SEARCH · Engineering Papers

Results for “numerical relativity”

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.

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Interannual Fluctuations in Atmospheric Angular Momentum Simulated by the National Centers for Environmental Prediction Medium Range Forecast Model

An earlier study established the existence of globally coherent interannual fluctuations in atmospheric angular momentum (AAM), associated with the El Nino-Southern Oscillation (ENSO) cycle. In this paper, we pursue the origin and the structure of these fluctuations using an ensemble of experiments generated by the National Centers for Environmental Prediction, medium range forecast model version 9. In the control experiments, where the observed sea surface temperatures (SSTs) were used as the lower boundary conditions, the model captures the characteristic V-like structure in time-latitude plots of zonally averaged AAM, while experiments with climatological SSTs and those with either perpetual warm or cold ENSO conditions superimposed on the climatological SSTs failed to reproduce this structure. The numerical results indicate that these AAM structures are related to SST variations associated with transitions between different phases of the ENSO cycle and have both propagating and standing components. The largest zonal wind contribution from the levels studied (850, 500, and 200 hPa) is at 200 hPa, where the tropical convective outflow is the strongest. Composites of zonal wind and geopotential height show a clear relationship between the stages of the global AAM oscillation and the ENSO cycle. The strong similarity between the simulated and observed AAM series attests to the model's ability to realistically simulate the interannual response of the atmosphere to ENSO SST anomalies.

Mo, Kingtse C.↗

NASA NTRS2000

Exchange Forces in Dispersion Relations Investigated Using Circuit Relations

We propose a novel method to compute in an exact manner the left-hand cut discontinuity of the electron-Atom partial wave scattering amplitude in the complex energy plane within the static stage approximation. Zero energy dispersion relations for electron-Hydrogen scattering are computed numerically for illustration.

Vrinceanu, D.↗

NASA NTRS1996

Hyperbolic Methods for Surface and Field Grid Generation

This chapter describes the use of hyperbolic partial differential equation methods for structured surface grid generation and field grid generation. While the surface grid generation equations are inherently three dimensional, the field grid generation equations can be formulated in two or three dimensions. The governing equations are derived from orthogonality relations and cell area/volume constraints; and are solved numerically by marching from an initial curve or surface. The marching step size and marching distance can be prescribedly the user. Exact specifications of the side and outer boundaries are not possible with a one sweep marching scheme but limited control is achievable. Excellent orthogonality and grid clustering characteristics are provided by hyperbolic methods with one to two orders of magnitude savings in time over typical elliptic methods. Since hyperbolic grid generation methods do not require the exact specifications of the side and outer boundaries of a grid, these methods are particularly well suited for the overlapping grid approach for solving problems on complex configurations. Grid generation software based on hyperbolic methods and their applications on several complex configurations will be described.

Chan, William M.↗

NASA NTRS2000

Self-Scheduling Parallel Methods for Multiple Serial Codes with Application to WOPWOP

This paper presents a scheme for efficiently running a large number of serial jobs on parallel computers. Two examples are given of computer programs that run relatively quickly, but often they must be run numerous times to obtain all the results needed. It is very common in science and engineering to have codes that are not massive computing challenges in themselves, but due to the number of instances that must be run, they do become large-scale computing problems. The two examples given here represent common problems in aerospace engineering: aerodynamic panel methods and aeroacoustic integral methods. The first example simply solves many systems of linear equations. This is representative of an aerodynamic panel code where someone would like to solve for numerous angles of attack. The complete code for this first example is included in the appendix so that it can be readily used by others as a template. The second example is an aeroacoustics code (WOPWOP) that solves the Ffowcs Williams Hawkings equation to predict the far-field sound due to rotating blades. In this example, one quite often needs to compute the sound at numerous observer locations, hence parallelization is utilized to automate the noise computation for a large number of observers.

Long, Lyle N.↗

NASA NTRS2000

Numerical Study of Wake Vortex Behavior in Turbulent Domains with Ambient Stratification

A three-dimensional large eddy simulation model is used to investigate the sensitivity of ambient stratification with turbulence on the behavior of aircraft wake vortices. Modeled ambient turbulence levels range from very weak to moderate, and stratification levels range from strongly stable to unstable. The results of profound significance from this study are: 1) very little sensitivity between vortex linking time and the level of stratification, 2) the mean vortex separation remained nearly constant regardless of stratification and turbulence (at least prior to linking), 3) the wake vortices did not rise regardless of the level of stratification, and 4) for very strong stratification, the vortex stopped descending and quickly dissipated even before vortex linking could occur. These results are supported by experimental data and are contrary to conclusions from other numerical studies that assume laminar flow and/or relatively-low Reynolds numbers.

Switzer, George F.↗

NASA NTRS2004

EARLY IMPACT MELTING AND SPACE EXPOSURE HISTORY OF THE PAT91501 LCHONDRITE

Collisions probably occurred frequently in the early history of the asteroid belt. Their effects, which should be recorded in meteorites, must have included heating and melting along with shock alteration of mineral textures. Some non-chondritic meteorite types e.g., eucrites and IIE and IAB irons - do indeed give evidence of extensive impact heating more than 3.4 Gyr ago. The ordinary chondrites, in contrast, show little evidence of early impact heating. The Ar-Ar and Rb-Sr ages of ordinary chondrites that experienced intense shock are for the most part relatively young, many less than 1.5 Gyr. The numerous L-chondrites with Ar- Ar ages clustering near 0.5 Gy are a well-known example. One of them, the 105-kg Chico Lchondrite, shows the effects of unusually intense heating. It is approximately 60% impact melt and likely formed as a dyke beneath a large crater when the L-chondrite parent body underwent a very large impact approximately 0.5 Gyr ago. In rare instances, older shock dates are indicated for ordinary chondrites. Dixon et al show early impact resetting of Ar-Ar ages of a few LL-chondrites including MIL 99301 at 4.23 0.03 Gyr, but in none of these stones did shock lead to extensive melting. As of 2003, searches for chondritic melts attributable to early shock had turned up only the Shaw L-chondrite, which has an Ar-Ar age of approximately 4.42 Gyr. PAT91501 is an 8.55-kg L-chondrite containing vesicles and metal-troilite nodules. It is a unique, near-total impact melt, unshocked, depleted in siderophile and chalcophile elements, and contains only approximately 10% relic chondritic material. The authors conclude that PAT91501 crystallized rapidly and from a much more homogeneous melt than did Shaw. They suggest that PAT resembles Chico and likely formed as an impact melt vein within an impact crater. To define the history of PAT, we have determined its Ar-39-Ar-40 age and measured several radioactive and stable nuclides produced during its space exposure to cosmic rays.

Bogard, Donald D.↗

NASA NTRS2005

Multi-reference evaluation of uncertainty in earth orientation parameter measurements

Uncertainties in polar motion and length-of-day measurements are evaluated empirically using several data series from the space-geodetic techniques of the global positioning system (GPS), satellite laser ranging (SLR), and very long baseline interferometry (VLBI) during 1997-2002. In the evaluation procedure employed here, known as the three-corner hat (TCH) technique, the signal common to each series is eliminated by forming pair-wise differences between the series, thus requiring no assumed values for the 'truth' signal. From the variances of the differenced series, the uncertainty of each series can be recovered when reasonable assumptions are made about the correlations between the series. In order to form the pair-wise differences, the series data must be given at the same epoch. All measurement data sets studied here were sampled at noon (UTC); except for the VLBI series, whose data are interpolated to noon and whose UT1 values are also numerically differentiated to obtain LOD. The numerical error introduced to the VLBI values by the interpolation and differentiation is shown to be comparable in magnitude to the values determined by the TCH method for the uncertainties of the VLBI series. The TCH estimates for the VLBI series are corrupted by such numerical errors mostly as a result of the relatively large data intervals. Of the remaining data sets studied here, it is found that the IGS Final combined series has the smallest polar motion and length-of-day uncertainties.

Combination↗

NASA NTRS2008

Shear Thinning Near the Critical Point of Xenon

We measured shear thinning, a viscosity decrease ordinarily associated with complex liquids, near the critical point of xenon. The data span a wide range of reduced shear rate: 10(exp -3) < gamma-dot tau < 700, where gamma-dot tau is the shear rate scaled by the relaxation time tau of critical fluctuations. The measurements had a temperature resolution of 0.01 mK and were conducted in microgravity aboard the Space Shuttle Columbia to avoid the density stratification caused by Earth's gravity. The viscometer measured the drag on a delicate nickel screen as it oscillated in the xenon at amplitudes 3 mu,m < chi (sub 0) >430 mu, and frequencies 1 Hz < omega/2 pi < 5 Hz. To separate shear thinning from other nonlinearities, we computed the ratio of the viscous force on the screen at gamma-dot tau to the force at gamma-dot tau approximates 0: C(sub gamma) is identical with F(chi(sub 0), omega tau, gamma-dot tau )/F)(chi(sub 0, omega tau, 0). At low frequencies, (omega tau)(exp 2) < gamma-dot tau, C(sub gamma) depends only on gamma-dot tau, as predicted by dynamic critical scaling. At high frequencies, (omega tau)(exp 2) > gamma-dot tau, C(sub gamma) depends also on both x(sub 0) and omega. The data were compared with numerical calculations based on the Carreau-Yasuda relation for complex fluids: eta(gamma-dot)/eta(0)=[1+A(sub gamma)|gamma-dot tau|](exp - chi(sub eta)/3+chi(sub eta)), where chi(sub eta) =0.069 is the critical exponent for viscosity and mode-coupling theory predicts A(sub gamma) =0.121. For xenon we find A(sub gamma) =0.137 +/- 0.029, in agreement with the mode coupling value. Remarkably, the xenon data close to the critical temperature T(sub c) were independent of the cooling rate (both above and below T(sub c) and these data were symmetric about T(sub c) to within a temperature scale factor. The scale factors for the magnitude of the oscillator s response differed from those for the oscillator's phase; this suggests that the surface tension of the two-phase domains affected the drag on the screen below T(sub c).

Zimmerli, Gregory A.↗

NASA NTRS2009

Examining the Impacts of High-Resolution Land Surface Initialization on Model Predictions of Convection in the Southeastern U.S.

One of the most challenging weather forecast problems in the southeastern U.S. is daily summertime pulse convection. During the summer, atmospheric flow and forcing are generally weak in this region; thus, convection typically initiates in response to local forcing along sea/lake breezes, and other discontinuities often related to horizontal gradients in surface heating rates. Numerical simulations of pulse convection usually have low skill, even in local predictions at high resolution, due to the inherent chaotic nature of these precipitation systems. Forecast errors can arise from assumptions within physics parameterizations, model resolution limitations, as well as uncertainties in both the initial state of the atmosphere and land surface variables such as soil moisture and temperature. For this study, it is hypothesized that high-resolution, consistent representations of surface properties such as soil moisture and temperature, ground fluxes, and vegetation are necessary to better simulate the interactions between the land surface and atmosphere, and ultimately improve predictions of local circulations and summertime pulse convection. The NASA Short-term Prediction Research and Transition (SPORT) Center has been conducting studies to examine the impacts of high-resolution land surface initialization data generated by offline simulations of the NASA Land Informatiot~ System (LIS) on subsequent numerical forecasts using the Weather Research and Forecasting (WRF) model (Case et al. 2008, to appear in the Journal of Hydrometeorology). Case et al. presents improvements to simulated sea breezes and surface verification statistics over Florida by initializing WRF with land surface variables from an offline LIS spin-up run, conducted on the exact WRF domain and resolution. The current project extends the previous work over Florida, focusing on selected case studies of typical pulse convection over the southeastern U.S., with an emphasis on improving local short-term WRF simulations over the Mobile, AL and Miami, FL NWS county warning areas. Future efforts may involve examining the impacts of assimilating remotely-sensed soil moisture data, and/or introducing weekly greenness vegetation fraction composites (as opposed to monthly climatologies) into ol'fline NASA LIS runs. Based on positive impacts, the offline LIS runs could be transitioned into an operational mode, providing land surface initialization data to NWS forecast offices in real time.

Case, Jonathan L.↗

NASA NTRS2010

A Modeling and Verification Study of Summer Precipitation Systems Using NASA Surface Initialization Datasets

One of the most challenging weather forecast problems in the southeastern U.S. is daily summertime pulse-type convection. During the summer, atmospheric flow and forcing are generally weak in this region; thus, convection typically initiates in response to local forcing along sea/lake breezes, and other discontinuities often related to horizontal gradients in surface heating rates. Numerical simulations of pulse convection usually have low skill, even in local predictions at high resolution, due to the inherent chaotic nature of these precipitation systems. Forecast errors can arise from assumptions within parameterization schemes, model resolution limitations, and uncertainties in both the initial state of the atmosphere and land surface variables such as soil moisture and temperature. For this study, it is hypothesized that high-resolution, consistent representations of surface properties such as soil moisture, soil temperature, and sea surface temperature (SST) are necessary to better simulate the interactions between the surface and atmosphere, and ultimately improve predictions of summertime pulse convection. This paper describes a sensitivity experiment using the Weather Research and Forecasting (WRF) model. Interpolated land and ocean surface fields from a large-scale model are replaced with high-resolution datasets provided by unique NASA assets in an experimental simulation: the Land Information System (LIS) and Moderate Resolution Imaging Spectroradiometer (MODIS) SSTs. The LIS is run in an offline mode for several years at the same grid resolution as the WRF model to provide compatible land surface initial conditions in an equilibrium state. The MODIS SSTs provide detailed analyses of SSTs over the oceans and large lakes compared to current operational products. The WRF model runs initialized with the LIS+MODIS datasets result in a reduction in the overprediction of rainfall areas; however, the skill is almost equally as low in both experiments using traditional verification methodologies. Output from object-based verification within NCAR s Meteorological Evaluation Tools reveals that the WRF runs initialized with LIS+MODIS data consistently generated precipitation objects that better matched observed precipitation objects, especially at higher precipitation intensities. The LIS+MODIS runs produced on average a 4% increase in matched precipitation areas and a simultaneous 4% decrease in unmatched areas during three months of daily simulations.

Jonathan L. Case↗

NASA NTRS2008

Metrology Optical Power Budgeting in SIM Using Statistical Analysis Techniques

The Space Interferometry Mission (SIM) is a space-based stellar interferometry instrument, consisting of up to three interferometers, which will be capable of micro-arc second resolution. Alignment knowledge of the three interferometer baselines requires a three-dimensional, 14-leg truss with each leg being monitored by an external metrology gauge. In addition, each of the three interferometers requires an internal metrology gauge to monitor the optical path length differences between the two sides. Both external and internal metrology gauges are interferometry based, operating at a wavelength of 1319 nanometers. Each gauge has fiber inputs delivering measurement and local oscillator (LO) power, split into probe-LO and reference-LO beam pairs. These beams experience power loss due to a variety of mechanisms including, but not restricted to, design efficiency, material attenuation, element misalignment, diffraction, and coupling efficiency. Since the attenuation due to these sources may degrade over time, an accounting of the range of expected attenuation is needed so an optical power margin can be book kept. A method of statistical optical power analysis and budgeting, based on a technique developed for deep space RF telecommunications, is described in this paper and provides a numerical confidence level for having sufficient optical power relative to mission metrology performance requirements.

stellar interferometry↗

NASA NTRS2015

Construction of Polarimetric Radar-Based Reference Rain Maps for the Iowa Flood Studies Campaign

The Global Precipitation Measurement (GPM) Mission Iowa Flood Studies (IFloodS) campaign was conducted in central and northeastern Iowa during the months of April-June, 2013. Specific science objectives for IFloodS included quantification of uncertainties in satellite and ground-based estimates of precipitation, 4-D characterization of precipitation physical processes and associated parameters (e.g., size distributions, water contents, types, structure etc.), assessment of the impact of precipitation estimation uncertainty and physical processes on hydrologic predictive skill, and refinement of field observations and data analysis approaches as they pertain to future GPM integrated hydrologic validation and related field studies. In addition to field campaign archival of raw and processed satellite data (including precipitation products), key ground-based platforms such as the NASA NPOL S-band and D3R Ka/Ku-band dual-polarimetric radars, University of Iowa X-band dual-polarimetric radars, a large network of paired rain gauge platforms, and a large network of 2D Video and Parsivel disdrometers were deployed. In something of a canonical approach, the radar (NPOL in particular), gauge and disdrometer observational assets were deployed to create a consistent high-quality distributed (time and space sampling) radar-based ground "reference" rainfall dataset, with known uncertainties, that could be used for assessing the satellite-based precipitation products at a range of space/time scales. Subsequently, the impact of uncertainties in the satellite products could be evaluated relative to the ground-benchmark in coupled weather, land-surface and distributed hydrologic modeling frameworks as related to flood prediction. Relative to establishing the ground-based "benchmark", numerous avenues were pursued in the making and verification of IFloodS "reference" dual-polarimetric radar-based rain maps, and this study documents the process and results as they pertain specifically to efforts using the NPOL radar dataset. The initial portions of the "process" involved dual-polarimetric quality control procedures which employed standard phase and correlation-based approaches to removal of clutter and non-meteorological echo. Calculation of a scale-adaptive KDP was accomplished using the method of Wang and Chandrasekar (2009; J. Atmos. Oceanic Tech.). A dual-polarimetric blockage algorithm based on Lang et al. (2009; J. Atmos. Oceanic Tech.) was then implemented to correct radar reflectivity and differential reflectivity at low elevation angles. Next, hydrometeor identification algorithms were run to identify liquid and ice hydrometeors. After the quality control and data preparation steps were completed several different dual-polarimetric rain estimation algorithms were employed to estimate rainfall rates using rainfall scans collected approximately every two to three minutes throughout the campaign. These algorithms included a polarimetrically-tuned Z-R algorithm that adjusts for drop oscillations (via Bringi et al., 2004, J. Atmos. Oceanic Tech.), and several different hybrid polarimetric variable approaches, including one that made use of parameters tuned to IFloodS 2D Video Disdrometer measurements. Finally, a hybrid scan algorithm was designed to merge the rain rate estimates from multiple low level elevation angle scans (where blockages could not be appropriately corrected) in order to create individual low-level rain maps. Individual rain maps at each time step were subsequently accumulated over multiple time scales for comparison to gauge network data. The comparison results and overall error character depended strongly on rain event type, polarimetric estimator applied, and range from the radar. We will present the outcome of these comparisons and their impact on constructing composited "reference" rainfall maps at select time and space scales.

Radar↗

NASA NTRS2016

System and Method for Outlier Detection via Estimating Clusters

An efficient method and system for real-time or offline analysis of multivariate sensor data for use in anomaly detection, fault detection, and system health monitoring is provided. Models automatically derived from training data, typically nominal system data acquired from sensors in normally operating conditions or from detailed simulations, are used to identify unusual, out of family data samples (outliers) that indicate possible system failure or degradation. Outliers are determined through analyzing a degree of deviation of current system behavior from the models formed from the nominal system data. The deviation of current system behavior is presented as an easy to interpret numerical score along with a measure of the relative contribution of each system parameter to any off-nominal deviation. The techniques described herein may also be used to "clean" the training data.

Iverson, David J.↗

NASA NTRS2017

Long Duration Head Down Tilt Bed Rest and Spaceflight Effects on Neurocognitive Performance: Extent, Longevity and Neural Bases

We have recently completed a long duration head down tilt bed rest (HDBR) study in which we performed structural and functional magnetic resonance brain imaging to identify the relationships between changes in neurocognitive function and neural structural alterations in a spaceflight analog environment. We are also collecting the same measures in crewmembers prior to and following a six month International Space Station mission. We will present data demonstrating that bed rest resulted in functional mobility and balance deterioration with recovery post-HDBR. We observed numerous changes in brain structure, function, and connectivity relative to a control group which were associated with pre to post bed rest changes in sensorimotor function. For example, gray matter volume (GMv) increased in posterior parietal areas and decreased in frontal regions. GMv increases largely overlapped with fluid decreases and vice versa. Larger increases in precentral gyrus (M1)/ postcentral gyrus (S1+2) GMv and fluid decreases were associated with smaller balance decrements. Vestibular activation in the bilateral insular cortex increased with bed rest and subsequently recovered. Larger increases in vestibular activation in multiple brain regions were associated with greater decrements in balance and mobility. We found connectivity increases between left M1 with right S1+2 and the superior parietal lobule, and right vestibular cortex with the cerebellum. Decreases were observed between right Lobule VIII with right S1+2 and the supramarginal gyrus, right posterior parietal cortex (PPC) with occipital regions, and the right superior posterior fissure with right Crus I and II. Connectivity strength between left M1 and right S1+2/superior parietal lobule increased the most in individuals that exhibited the least balance impairments. In sum, we observed HDBR-related changes in measures of brain structure, function, and network connectivity, which correlated with indices of sensorimotor function. Recovery was observed post HDBR but remained incomplete at 12 days post-HDBR. Preliminary findings from our parallel ongoing flight study will be compared and contrasted with bed rest results during this presentation.

Seidler, R. D.↗

NASA NTRS2021

Quirky Circling Behavior in Mice Informs Research on Humans in Space

As interest in long duration effects of space habitation increases, understanding the behavior of model organisms living within the habitats engineered to fly them is vital for designing, validating, and interpreting future spaceflight studies. We previously conducted a detailed phenotypic analysis of mouse behavior during long duration (33-day) spaceflight in the NASA Rodent Habitat (Rodent Research-1; RR1 mission). Notably, we documented the emergence within the 8-10 days of launch of spontaneous ambulatory behavior in the form of circling or ‘race-tracking’ behavior in spaceflight but not in an identical ground control condition. Circling is unique to the NASA RH not having been reported to occur in other mouse habitats that provide limited opportunity for grasping, and/or are characterized by smaller habitable volumes viz., Italian Mice Drawer System (MDS) flown on ISS; Russian Block Obespecheniya Soderzhaniya (BOS) flown on Bion M-1, or JAXA Habitat Cage Unit (HCU) flown in the ISS Kibo module. Over time, circling became the primary dark cycle activity of FLT mice, occurring in individuals then as a coordinated group activity. Here we discuss possible interpretations of circling behavior including: (1) Stereotypic or abnormal repetitive behaviors (ARBs) that are unvarying, and apparently functionless behavior patterns documented numerous species in laboratory and zoo settings, possibility related to insufficient environmental enrichment and stress, (2) Rewarding effects of physical activity that are well-documented in terrestrial studies of rodents given the opportunity to run in wheels, and (3) Vestibular self-stimulation, i.e., the generation biologically-relevant amounts of vestibular sensory input to reduce the effects of microgravity. Affording mice the opportunity to grab and run in the RH resembles physical activities that the crew participate in routinely. Our approach is yielding an interesting analogue for better understanding human responses to spaceflight, and providing the opportunity to begin to address how physical movement influences responses to microgravity.

spaceflight↗

NASA NTRSDate not supplied

Utilization of Streamtubes to Analyze the Physical Interaction of a Dispersed Cloud with the CRM65 Hybrid Midspan Model

Recently there have been numerous efforts to identify the relevant icing physics related to the formation of complex three dimensional features, sometimes referred to as ’scallops’, on swept wings. However, much of the physics is still not well understood. This paper computationally investigates the interaction of the icing cloud with the 65 percent Common Research Model (CRM65). Both the interaction with an uniced model, and the interaction with a representative simulated three dimensional ice accretion are analyzed. Preliminary results suggest that the liquid water content increases near the aerodynamic body. For small droplets on the uniced geometry, the particle velocity vector becomes nearly parallel with the aerodynamic model resulting in a low value of collection efficiency. When a three dimensional feature with a length scale much smaller than the leading edge of the airfoil is introduced into the flow, the impingement of small particles can become significantly more perpendicular to the particle velocity vector in the region of these features. Since the liquid water content near the body has increased due to the interaction with the larger scale features of the aerodynamic model, i.e. the leading edge of the swept wing, extremely high collection efficiency is observed. These results suggest that three dimensional features are likely a significant physical driver that should be modeled in some capacity when simulating the impingement of a cloud on an aerodynamic model.

Icing↗

NASA NTRSDate not supplied

Utilization of Streamtubes to Analyze the Physical Interaction of a Dispersed Cloud with the CRM65 Hybrid Midspan Model

Icing↗

DOE OSTI2021

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗