Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

NASA GPM GV Science Implementation

Pre-launch algorithm development & post-launch product evaluation: The GPM GV paradigm moves beyond traditional direct validation/comparison activities by incorporating improved algorithm physics & model applications (end-to-end validation) in the validation process. Three approaches: 1) National Network (surface): Operational networks to identify and resolve first order discrepancies (e.g., bias) between satellite and ground-based precipitation estimates. 2) Physical Process (vertical column): Cloud system and microphysical studies geared toward testing and refinement of physically-based retrieval algorithms. 3) Integrated (4-dimensional): Integration of satellite precipitation products into coupled prediction models to evaluate strengths/limitations of satellite precipitation producers.

Petersen, W. A.↗
NASA NTRS2003

The Nonlinear Coupling of Electromagnetic Ion Cyclotron and Lower Hybrid Waves in the Ring Current Region: The Magnetic Storm May 1-7 1998

The excitation of lower hybrid waves (LHWs) is a widely discussed mechanism of interaction between plasma species in space, and is one of the unresolved questions of magnetospheric multi-ion plasmas. In this paper we present the morphology, dynamics, and level of LHW activity generated by electromagnetic ion cyclotron (EMIC) waves during the May 2-7, 1998 storm period on the global scale. The LHWs were calculated based on our newly developed self-consistent model that couples the system of two kinetic equations: one equation describes the ring current (RC) ion dynamic, and another equation describes the evolution of EMIC waves. It is found that the LHWs are excited by helium ions due to their mass dependent drift in the electric field of EMIC waves. The level of LHW activity is calculated assuming that the induced scattering process is the main saturation mechanism for these waves. The calculated LHWs electric fields are consistent with the observational data.

Khazanov, G. V.↗
NASA NTRS1990

The role of high-speed plasma flows in plasmaspheric refilling

A model of time-dependent one-stream interhemispheric plasma flow is used to investigate plasmaspheric refilling. In the model the coupled time-dependent hydrodynamic equations (continuity, momentum and energy) of a two-ion, H(+) and O(+), quasi-neutral, currentless plasma are solved for a closed geomagnetic field line. For the present set of calculations an L = 2 field line was used. A steady state solution was found and used as the initial condition in subsequent simulations of the effect of density depletions. Density depletions were modeled by reducing the densities by an arbitrary factor above 2500 km altitude, while keeping the velocities and temperatures unchanged. Shock structures develop which move up the field line, meeting and reflecting at the equator. The reflected shocks are absorbed by the dense neutral atmosphere when they reach the model flux tube boundaries. After the upwelling streams collide, refilling occurs from the equator downward for 2 to 3 hours; thereafter there is steady refilling from the ionospheres upward.

Guiter, S. M.↗
NASA NTRSDate not supplied

Challenges in Observing, Modeling, and Forecasting the June 2023 Smoke Event Over the Northeast United States

In the first week of May 2023, boreal Canada began an early start to the biomass burning season when fires erupted across Alberta and western Saskatchewan. Smoke entered the troposphere and spread across the United States as biomass burning emissions in Canada quadrupled the prior maximum since the onset of the MODIS satellite record. Adding to the already existing smoke, wildfires ignited in Quebec on June 2, 2023, producing a heavy plume of aerosol over the northeastern part of the continent. Meteorological factors, including an area of low pressure situated, nearly stationary, over Maine, transported the smoke into a densely populated corridor containing the cities of Washington, DC, Baltimore, Philadelphia, and New York City. Air quality alerts became widespread as near-surface levels of fine particulate matter (PM2.5) exceeded harmful levels due the low altitude of the smoke plume. Forecasting the smoke transport and impacts of such events requires a complex combination of observed initial conditions for meteorology and aerosols, knowledge of the emissions from wildfires, and a state-of-the-art Earth System model that couples these components together. Using a case-study perspective with the Goddard Earth Observing System (GEOS), challenges associated with near real time forecasting of the smoke plume are presented. Analyzed and forecasted aerosol optical depth will be evaluated using available satellite and AERONET observations, while surface aerosol will be assessed relative to observations of PM2.5. Lacking real time estimates of biomass burning emissions for inputs to our forecast model, particular attention is given to the use of day-old biomass burning emissions throughout the forecast period covering the first ten days in June 2023, and the resulting reduction in forecast skill will be quantified.

Allison Collow↗
NASA NTRSDate not supplied

Challenges in Observing, Modeling, and Forecasting the June 2023 Smoke Event over the Northeast United States

In the first week of May 2023, boreal Canada began an early start to the biomass burning season when fires erupted across Alberta and western Saskatchewan. Smoke entered the troposphere and spread across the United States as biomass burning emissions in Canada quadrupled the prior maximum since the onset of the MODIS satellite record. Adding to the already existing smoke, wildfires ignited in Quebec on June 2, 2023, producing a heavy plume of aerosol over the northeastern part of the continent. Meteorological factors, including an area of low pressure situated, nearly stationary, over Maine, transported the smoke into a densely populated corridor containing the cities of Washington, DC, Baltimore, Philadelphia, and New York City. Air quality alerts became widespread as near-surface levels of fine particulate matter (PM2.5) exceeded harmful levels due the low altitude of the smoke plume. Forecasting the smoke transport and impacts of such events requires a complex combination of observed initial conditions for meteorology and aerosols, knowledge of the emissions from wildfires, and a state-of-the-art Earth System model that couples these components together. Using a case-study perspective with the Goddard Earth Observing System (GEOS ), challenges associated with near real time forecasting of the smoke plume are presented. Analyzed and forecasted aerosol optical depth will be evaluated using available satellite and AERONET observations, while surface aerosol will be assessed relative to observations of PM2.5. Lacking real time estimates of biomass burning emissions for inputs to our forecast model, particular attention is given to the use of day-old biomass burning emissions throughout the forecast period covering the first ten days in June 2023, and the resulting reduction in forecast skill will be quantified.

Allison Collow↗
NASA NTRS1978

Effects of nonmethane hydrocarbons in the atmosphere

An investigation was conducted to determine whether nonmethane hydrocarbons (NMHC) are abundant enough to have a significant impact upon the ambient photochemistry. The vertical distribution of C2H6, C2H2, C3H8, C4H10, and C5H12 in the altitude range from 0 to 40 km was calculated in this connection. A one-dimensional steady state model with coupled photochemistry and vertical transport was employed in the investigation. The calculations imply that measurable quantities of relatively unreactive NMHC, especially C2H6 and C2H2, may be present in the upper troposphere and stratosphere. The results indicate, however, that NMHC are not likely to have a large impact on the background photochemistry of the troposphere, although local effects near source regions are probable. The findings support the current practice of many modelers who neglect NMHC in their calculations.

Chameides, W. L.↗
NASA NTRS2001

The Role of Mesoscale Variability on Plankton Dynamics in the North Atlantic

The intensive field observational phase of JGOFS in the North Atlantic Ocean has shown the importance of oceanic mesoscale variability on biogeochemical cycles and on the strength of the ocean biological pump. Mesoscale physical dynamics govern the major time/space scales of bulk biological variability (biomass, production, and export). Mesoscale eddies seem to have a strong impact on the ecosystem structure and functioning, but observational evidence is rather limited. For the signature of the mesoscale features to exist in the ecosystem, the comparison of temporal scales of formation and evolution of mesoscale features and reaction of the ecosystem is a key factor. Biological patterns are driven by active changes in biological source and sink terms rather than simply by passive turbulent mixing. A first modelling assessment of the regional balances between horizontal and vertical eddy-induced nutrient supplies in the euphotic zone shows that the horizontal transport predominates over the vertical route in the subtropical gyre, whereas the reverse holds true for the other biogeochemical provinces of the North Atlantic. Presently, despite some difference in numbers, the net impact of modelled eddies yields an enhancement of the biological productivity in most provinces of the North Atlantic Ocean. Key issues remaining include variation on the mesoscale of subsurface particle and dissolved organic matter remineralization, improved knowledge of the ecological response to patterns of variability, synopticity in mesoscale surveys along with refining measures of biogeochemical time/space variability. Eventual success of assimilation of in situ and satellite data, still in its infancy in coupled physical/biogeochemical models, will be crucial to achieve JGOFS synthesis in answering which data are most informative, standing stocks or rates, and which ones are relevant. Depending on which end of the spectrum quantification of the effect of mesoscale features on production and community structure is required, complementary strategies are offered. Either one may choose to increase resolution of models up to the very fine mesoscale features scale (a few kms) for the high end, or to include a parametric representation of eddies for the low end.

Garcon, Veronique C.↗
NASA NTRS2002

The Role of Mesoscale Variability on Plankton Dynamics in the North Atlantic

The intensive field observational phase of JGOFS in the North Atlantic Ocean has shown the importance of oceanic mesoscale variability on biogeochemical cycles and on the strength of the ocean biological pump. Mesoscale physical dynamics govern the major time/space scales of bulk biological variability (biomass, production and export). Mesoscale eddies seem to have a strong impact on the ecosystem structure and functioning, but observational evidence is rather limited. For the signature of the mesoscale features to exist in the ecosystem, the comparison of temporal scales of formation and evolution of mesoscale features and reaction of the ecosystem is a key factor. Biological patterns are driven by active changes in biological source and sink terms rather than simply by passive turbulent mixing. A first modelling assessment of the regional balances between horizontal and vertical eddy-induced nutrient supplies in the euphotic zone shows that the horizontal transport predominates over the vertical route in the subtropical gyre, whereas the reverse holds true for the other biogeochemical provinces of the North Atlantic. Presently. despite some difference in numbers, the net impact of modelled eddies yields an enhancement of the biological productivity in most provinces of the North Atlantic Ocean. Key issues remaining include variation on the mesoscale of subsurface particle and dissolved organic matter remineralization, improved knowledge of the ecological response to patterns of variability, synopticity in mesoscale surveys along with refining measures of biogeochemical time/space variability. Eventual success of assimilation of in situ and satellite data, still in its infancy in coupled physical/biogeochemical models, will be crucial to achieve JGOFS synthesis in answering which data are most informative, standing stocks or rates, and which ones are relevant. Depending on which end of the spectrum quantification of the effect of mesoscale features on production and community structure is required, complementary strategies are offered. Either one may choose to increase resolution of models up to the very fine mesoscale features scale (a few kms) for the high end, or to include a parametric representation of eddies for the low end.

Garcon, Veronique, C.↗
NASA NTRS2012

Validation of Atmosphere/Ionosphere Signals Associated with Major Earthquakes by Multi-Instrument Space-Borne and Ground Observations

The latest catastrophic earthquake in Japan (March 2011) has renewed interest in the important question of the existence of pre-earthquake anomalous signals related to strong earthquakes. Recent studies have shown that there were precursory atmospheric/ionospheric signals observed in space associated with major earthquakes. The critical question, still widely debated in the scientific community, is whether such ionospheric/atmospheric signals systematically precede large earthquakes. To address this problem we have started to investigate anomalous ionospheric / atmospheric signals occurring prior to large earthquakes. We are studying the Earth's atmospheric electromagnetic environment by developing a multisensor model for monitoring the signals related to active tectonic faulting and earthquake processes. The integrated satellite and terrestrial framework (ISTF) is our method for validation and is based on a joint analysis of several physical and environmental parameters (thermal infrared radiation, electron concentration in the ionosphere, lineament analysis, radon/ion activities, air temperature and seismicity) that were found to be associated with earthquakes. A physical link between these parameters and earthquake processes has been provided by the recent version of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model. Our experimental measurements have supported the new theoretical estimates of LAIC hypothesis for an increase in the surface latent heat flux, integrated variability of outgoing long wave radiation (OLR) and anomalous variations of the total electron content (TEC) registered over the epicenters. Some of the major earthquakes are accompanied by an intensification of gas migration to the surface, thermodynamic and hydrodynamic processes of transformation of latent heat into thermal energy and with vertical transport of charged aerosols in the lower atmosphere. These processes lead to the generation of external electric currents in specific regions of the atmosphere and the modifications, by dc electric fields, in the ionosphere-atmosphere electric circuit. We retrospectively analyzed temporal and spatial variations of four different physical parameters (gas/radon counting rate, lineaments change, long-wave radiation transitions and ionospheric electron density/plasma variations) characterizing the state of the lithosphere/atmosphere coupling several days before the onset of the earthquakes. Validation processes consist in two phases: A. Case studies for seven recent major earthquakes: Japan (M9.0, 2011), China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and, Chile (M8.8, 2010) and B. A continuous retrospective analysis was preformed over two different regions with high seismicity- Taiwan and Japan for 2003-2009. Satellite, ground surface, and troposphere data were obtained from Terra/ASTER, Aqua/AIRS, POES and ionospheric variations from DEMETER and COSMIC-I data. Radon and GPS/TEC were obtaining from monitoring sites in Taiwan, Japan and Italy and from global ionosphere maps (GIM) respectively. Our analysis of ground and satellite data during the occurrence of 7 global earthquakes has shown the presence of anomalies in the atmosphere. Our results for Tohoku M9.0 earthquake show that on March 7th, 2011 (4 days before the main shock and 1 day before the M7.2 foreshock of March 8, 2011) a rapid increase of emitted infrared radiation was observed by the satellite data and an anomaly was developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. Similar approach for analyzing atmospheric and ionospheric parameters has been applied for China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and Chile (M8.8, 2010) eahquakes. Results have revealed the presence of related variations of these parameters implying their connection with the earthquake process. The second phase (B) of this validation included 102 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. False alarm ratio for false positives is less then 10% and has been calculated for the same month of the earthquake occurrence for the entire period of analysis (2003-2009). The commonalities for detecting atmospheric/ionospheric anomalies are: i.) Regularly appearance over regions of maximum stress (i.e., along plate boundaries); ii.) Anomaly existence over land and sea; and iii) association with M>5.9 earthquakes not deeper than 100km. Due to their long duration over the same region these anomalies are not consistent with a meteorological origin. Our initial results from the ISTF validation of multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area, one to seven (average) days prior to the largest earthquakes, and suggest that it could be explained by a coupling process between the observed physical parameters and the pre-earthquake preparation processes.

Ouzounov, Dimitar↗
NASA NTRS2017

Comparing Proxy and Model Estimates of Hydroclimate Variability and Change over the Common Era

Water availability is fundamental to societies and ecosystems, but our understanding of variations in hydroclimate (including extreme events, flooding, and decadal periods of drought) is limited because of a paucity of modern instrumental observations that are distributed unevenly across the globe and only span parts of the 20th and 21st centuries. Such data coverage is insufficient for characterizing hydroclimate and its associated dynamics because of its multidecadal to centennial variability and highly regionalized spatial signature. High-resolution (seasonal to decadal) hydroclimatic proxies that span all or parts of the Common Era (CE) and paleoclimate simulations from climate models are therefore important tools for augmenting our understanding of hydroclimate variability. In particular, the comparison of the two sources of information is critical for addressing the uncertainties and limitations of both while enriching each of their interpretations. We review the principal proxy data available for hydroclimatic reconstructions over the CE and highlight the contemporary understanding of how these proxies are interpreted as hydroclimate indicators. We also review the available last-millennium simulations from fully coupled climate models and discuss several outstanding challenges associated with simulating hydroclimate variability and change over the CE. A specific review of simulated hydroclimatic changes forced by volcanic events is provided, as is a discussion of expected improvements in estimated radiative forcings, models, and their implementation in the future. Our review of hydroclimatic proxies and last-millennium model simulations is used as the basis for articulating a variety of considerations and best practices for how to perform proxy-model comparisons of CE hydroclimate. This discussion provides a framework for how best to evaluate hydroclimate variability and its associated dynamics using these comparisons and how they can better inform interpretations of both proxy data and model simulations.We subsequently explore means of using proxy-model comparisons to better constrain and characterize future hydroclimate risks. This is explored specifically in the context of several examples that demonstrate how proxy-model comparisons can be used to quantitatively constrain future hydroclimatic risks as estimated from climate model projections.

Common Era↗
NASA NTRS2000

NASA/GSFC Research Activities for the Global Ocean Carbon Cycle: A Prospectus for the 21st Century

There are increasing concerns that anthropogenic inputs of carbon dioxide into the Earth system have the potential for climate change. In response to these concerns, the GSFC Laboratory for Hydrospheric Processes has formed the Ocean Carbon Science Team (OCST) to contribute to greater understanding of the global ocean carbon cycle. The overall goals of the OCST are to: 1) detect changes in biological components of the ocean carbon cycle through remote sensing of biooptical properties, 2) refine understanding of ocean carbon uptake and sequestration through application of basic research results, new satellite algorithms, and improved model parameterizations, 3) develop and implement new sensors providing critical missing environmental information related to the oceanic carbon cycle and the flux of CO2 across the air-sea interface. The specific objectives of the OCST are to: 1) establish a 20-year time series of ocean color, 2) develop new remote sensing technologies, 3) validate ocean remote sensing observations, 4) conduct ocean carbon cycle scientific investigations directly related to remote sensing data, emphasizing physiological, empirical and coupled physical/biological models, satellite algorithm development and improvement, and analysis of satellite data sets. These research and mission objectives are intended to improve our understanding of global ocean carbon cycling and contribute to national goals by maximizing the use of remote sensing data.

Gregg, W. W.↗
NASA NTRS1993

Finite element calculations of transonic flutter in cascades

Aeroelastic stability and response of turbomachinery cascades operating in the transonic regime are explored using a new computational approach. The two-dimensional unsteady Euler equations are used to model the transonic cascade flows, and the cascade blades are modeled as thin plates. To minimize the discretization incompatibility, both the fluid and the solid domains are discretized using finite element procedures. The fluid and the structural models are coupled at the element level by introducing a channel theory, which is analogous to a strip theory in the classical sense. The system of space-discretized equations is marched forward in time using an explicit four-stage Runge-Kutta scheme. Results of free vibration analyses of typical blades are presented and compared to published finite element and Rayleigh-Ritz results. Aeroelastic stability and response characteristics of several transonic cascades are investigated using the new quasi-3D model. These results demonstrate the present 3D finite element blade model is capable of capturing aeroelastic instabilities involving higher modes, including (plate-type) camber bending modes.

Hsiao, Chingteng↗
NASA NTRS1993

Optical system storage design with diffractive optical elements

Optical data storage systems are gaining widespread acceptance due to their high areal density and the ability to remove the high capacity hard disk from the system. In magneto-optical read-write systems, a small rotation of the polarization state in the return signal from the MO media is the signal which must be sensed. A typical arrangement used for detecting these signals and correcting for errors in tracking and focusing on the disk is illustrated. The components required to achieve these functions are listed. The assembly and alignment of this complex system has a direct impact on cost, and also affects the size, weight, and corresponding data access rates. As a result, integrating these optical components and improving packaging techniques is an active area of research and development. Most designs of binary optic elements have been concerned with optimizing grating efficiency. However, rigorous coupled wave models for vector field diffraction from grating surfaces can be extended to determine the phase and polarization state of the diffracted field, and the design of polarization components. A typical grating geometry and the phase and polarization angles associated with the incident and diffracted fields are shown. In our current stage of work, we are examining system configurations which cascade several polarization functions on a single substrate. In this design, the beam returning from the MO disk illuminates a cascaded grating element which first couples light into the substrate, then introduces a quarter wave retardation, then a polarization rotation, and finally separates s- and p-polarized fields through a polarization beam splitter. The input coupler and polarization beam splitter are formed in volume gratings, and the two intermediate elements are zero-order elements.

Kostuk, Raymond K.↗
NASA NTRS1975

Propagation predictions and studies using a ray tracing program combined with a theoretical ionospheric model

Radio wave propagation predictions are described in which modern comprehensive theoretical ionospheric models are coupled with ray-tracing programs. In the computer code described, a network of electron density and collision frequency parameters along a band about the great circle path is calculated by specifying the transmitter and receiver geographic coordinates, time, the day number, and the 2800-MHz solar flux. The ray paths are calculated on specifying the frequency, mode, range of elevation angles, and range of azimuth angles from the great circle direction. The current program uses a combination of the Penn State MKI E and F region models and the Mitra-Rowe D and E region model. Application of the technique to the prediction of satellite to ground propagation and calculation of oblique incidence propagation paths and absorption are described. The implications of the study to the development of the next generation of ionospheric models are discussed.

Lee, M. K.↗
NASA NTRS1991

A dynamical explanation for the asymmetry in zonally averaged column abundances of ozone between northern and southern springs

The observed zonally averaged column ozone shows a maximum at 90 deg N during the northern winter and spring and at 60 deg S throughout the southern winter and spring. This asymmetry is explained in the context of a zonally averaged model with coupled radiation, dynamics, and chemistry, together with consistently parameterized planetary wave driving and wave transport. It is shown that in the presence of weak wave driving, the penetration of the tropospheric circulation into the lower stratosphere and the characteristics of ozone chemistry are such that they produce a column ozone maximum at subpolar latitudes. The effect of increased wave driving is to intensify the residual circulation and extend it farther poleward, resulting in an ozone maximum at the pole. The role of the mesospheric drag is to further enhance these column ozone maxima. Model calculations show that the positions of the observed column ozone maxima are consistent with intensities of wave driving in the two hemispheres derived from data.

Hou, Arthur Y.↗
NASA NTRS1996

Integrated Control Using the SOFFT Control Structure

The need for integrated/constrained control systems has become clearer as advanced aircraft introduced new coupled subsystems such as new propulsion subsystems with thrust vectoring and new aerodynamic designs. In this study, we develop an integrated control design methodology which accomodates constraints among subsystem variables while using the Stochastic Optimal Feedforward/Feedback Control Technique (SOFFT) thus maintaining all the advantages of the SOFFT approach. The Integrated SOFFT Control methodology uses a centralized feedforward control and a constrained feedback control law. The control thus takes advantage of the known coupling among the subsystems while maintaining the identity of subsystems for validation purposes and the simplicity of the feedback law to understand the system response in complicated nonlinear scenarios. The Variable-Gain Output Feedback Control methodology (including constant gain output feedback) is extended to accommodate equality constraints. A gain computation algorithm is developed. The designer can set the cross-gains between two variables or subsystems to zero or another value and optimize the remaining gains subject to the constraint. An integrated control law is designed for a modified F-15 SMTD aircraft model with coupled airframe and propulsion subsystems using the Integrated SOFFT Control methodology to produce a set of desired flying qualities.

Halyo, Nesim↗
NASA NTRS2000

Towards a Comprehensive Dynamic-chemistry Assimilation for Eos-Chem: Plans and Status in NASA's Data Assimilation Office

In order to support the EOS-Chem project, a comprehensive assimilation package for the coupled chemical-dynamical system is being developed by the Data Assimilation Office at NASA GSFC. This involves development of a coupled chemistry/meteorology model and of data assimilation techniques for trace species and meteorology. The model is being developed using the flux-form semi-Lagrangian dynamical core of Lin and Rood, the physical parameterizations from the NCAR Community Climate Model, and atmospheric chemistry modules from the Atmospheric Chemistry and Dynamics branch at NASA GSFC. To date the following results have been obtained: (i) multi-annual simulations with the dynamics-radiation model show the credibility of the package for atmospheric simulations; (ii) initial simulations including a limited number of middle atmospheric trace gases reveal the realistic nature of transport mechanisms, although there is still a need for some improvements. Samples of these results will be shown. A meteorological assimilation system is currently being constructed using the model; this will form the basis for the proposed meteorological/chemical assimilation package. The latter part of the presentation will focus on areas targeted for development in the near and far terms, with the objective of Providing a comprehensive assimilation package for the EOS-Chem science experiment. The first stage will target ozone assimilation. The plans also encompass a reanalysis (ReSTS) for the 1991-1995 period, which includes the Mt. Pinatubo eruption and the time when a large number of UARS observations were available. One of the most challenging aspects of future developments will be to couple theoretical advances in tracer assimilation with the practical considerations of a real environment and eventually a near-real-time assimilation system.

Pawson, Steven↗
NASA NTRS2020

INVESTIGATING THE FEASIBILITY AND STABILITY FOR MODELING ACOUSTIC WAVE SCATTERING USING A TIME-DOMAIN BOUNDARY INTEGRAL EQUATION WITH IMPEDANCE BOUNDARY CONDITION

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. This work considers using either an impedance or an admittance (inverse of impedance) boundary condition to allow for acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces. Three acoustic liner models are discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and the Broadband Impedance Model. In both the Helmholtz and Three-Parameter models, liner impedance is specified at a given frequency, whereas the Broadband model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for acoustic liners are derived for each model and coupled with a time-domain boundary integral equation. The scattering solution is obtained iteratively using a boundary element method with constant spatial and third-order temporal basis functions. Time-domain boundary integral equations are unfortunately prone to numerical instabilities due to resonant frequencies resulting from non-trivial solutions in the interior domain. When reformulated with the Burton-Miller method, the instabilities are eliminated. Using a Burton-Miller reformulation, the stability of the boundary element method assuming a liner boundary condition is assessed using eigenvalue analysis. The stability of each liner model is discussed, and it is shown that the Three-Parameter and Broadband models are sufficient for modeling an acoustic liner on the surface of scattering bodies. The Helmholtz model demonstrates strict limitations for stability, whereas the Three-Parameter and Broadband models are stable for most cases. Also included in this work is an assessment of the spatial accuracy of the time-domain boundary element method with respect to the surface element basis functions, as well as a performance study of the numerical algorithm.

computational mathematics↗