SEARCH · Search NASA

Results for “Coupled process modeling”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Wind Tunnel Testing of the NASA-DFRC Flutterometer using a Two DOF Wing Section

Flutter of an aeroelastic structure is potentially destructive aeroelastic instability. This phenomenon has motivated research within the aeroelastic community to develop methods that can accurately predict aeroelastic instabilities. The Flutterometer method used herein, and as developed by NASA DFRC, is based upon the mu method which has been coupled with wavelet filtering processes in estimating aeroelastic models from flight data. The approach leads to a methodology to predict the occurrence of flutter boundaries, and may prove to reliably predict flutter boundaries during flight tests. An analytical model is used as the first estimate of the aeroelastic structural dynamics, and uncertainty operators are introduced into the system to model variations between the theoretical system and the physical system. The modelling uncertainties are then updated from experimental data. Although the model used did not work well with this particular experiment, a sensitivity analysis was additionally performed and improvements suggested.

Strganac, Thomas W.↗

NASA NTRS2016

Stratospheric Aerosol--Observations, Processes, and Impact on Climate

Interest in stratospheric aerosol and its role in climate have increased over the last decade due to the observed increase in stratospheric aerosol since 2000 and the potential for changes in the sulfur cycle induced by climate change. This review provides an overview about the advances in stratospheric aerosol research since the last comprehensive assessment of stratospheric aerosol was published in 2006. A crucial development since 2006 is the substantial improvement in the agreement between in situ and space-based inferences of stratospheric aerosol properties during volcanically quiescent periods. Furthermore, new measurement systems and techniques, both in situ and space based, have been developed for measuring physical aerosol properties with greater accuracy and for characterizing aerosol composition. However, these changes induce challenges to constructing a long-term stratospheric aerosol climatology. Currently, changes in stratospheric aerosol levels less than 20% cannot be confidently quantified. The volcanic signals tend to mask any nonvolcanically driven change, making them difficult to understand. While the role of carbonyl sulfide as a substantial and relatively constant source of stratospheric sulfur has been confirmed by new observations and model simulations, large uncertainties remain with respect to the contribution from anthropogenic sulfur dioxide emissions. New evidence has been provided that stratospheric aerosol can also contain small amounts of nonsulfatematter such as black carbon and organics. Chemistry-climate models have substantially increased in quantity and sophistication. In many models the implementation of stratospheric aerosol processes is coupled to radiation and/or stratospheric chemistry modules to account for relevant feedback processes.

Kresmer, Stefanie↗

NASA NTRS2022

VULCAN-CFD Theory Manual: Ver. 7.2.0

VULCAN-CFD Theory Manual↗

NASA NTRS2022

VULCAN-CFD User Manual: Ver. 7.2.0

VULCAN-CFD offers a comprehensive set of capabilities to enable the simulation of continuum flowfields from subsonic to hypersonic conditions. The governing equations that are employed include allowances for both chemical and thermal nonequilibrium processes, coupled with a wide variety of turbulence models for both Reynolds-averaged and large eddy simulations. The software package can simulate two-dimensional, axisymmetric, or three-dimensional problems on structured multiblock meshes or fully unstructured meshes. A parabolic (i.e., space-marching) treatment can also be used for any subset of a structured mesh that can accommodate this solution strategy. The flow solver provides a significant level of geometric flexibility for structured grid simulations by allowing for arbitrary face-to-face C(0) continuous and non-C(0) continuous block interface connectivities. The unstructured grid paradigm allows for mixed element unstructured meshes that contain any combination of tetrahedral, prismatic, pyramidal, and hexahedral cell elements. The flow solver is also fully parallelized using MPI (Message Passing Interface) libraries in a data-parallel fashion, allowing for efficient simulations on modern High Performance Computing (HPC) systems. This document provides information related to the installation and execution of the VULCAN-CFD software package. A detailed description of the physical and numerical models available in the software are provided in the VULCAN-CFD Theory Manual.

VULCAN-CFD User Manual↗

NASA NTRS1994

Erosional dynamics, flexural isostasy, and long-lived escarpments: A numerical modeling study

Erosional escarpments common features of high-elevation rifted continets. Fission track data suffest that these escarpments form by base level lowering and/or marginal uplift during rifting, followed by lateral retreat of an erosion front across tens to hundreds of kioometers. Previous modeling studies have shown that this characteristic pattern of denudation can have a profound impact upon marginal isostatic uplift and the evolution of offshore sedimentary basins. Yet at present there is only a rudimentary understanding of the geomorphic mechanisms capable of driving such prolonged escarpment retreat. In this study we present a nonlinear, two-dimensional landscape evolution model tha tis used to asses the necessary and sufficient conditions for long-term retreat of a rift-generated escarpment. The model represents topography as a grid of cells, with drainage networkds evolving as water flows across the grid in the direction of steepest descent. The model accounts for sediment production by weathering, fluvial sediment transport, bedrock channel erosion, and hillslope sediment transport by diffusive mechanisms and by mass failure. Numerical experiments presented explore the effects of different combinations of erosion processes and of dynamic coupling between denudation and flexural isostatic uplift. Model results suggest that the necessary and sufficient conditions for long-term escarpment retreat are (1) incising bedrock channels in which the erosion rate increases with increasing drainage area, so that the channels steepen and propagate headward; (2) a low rate of sediment production relative to sediment transport efficiency, which promotes relief-generating processes over diffusive ones; (3) high continental elevation, which allows greater freedom for fluvial dissection; and (4) any process, including flexural isostatic uplift, that helps to maintain a drainage divide near an escarpment crest. Flexural isostatic uplift also facilitates escarpment, thereby increasing channel gradients and accelerating erosion which in turn generates additional isostatic uplift. Of all the above conditions, high continental elevation is common ot most rift margin escarpments and may ultimately be the most important factor.

Tucker, Gregory E.↗

NASA NTRS2023

Convective Impact on the Global Lower Stratospheric Water Vapor Budget

Water vapor in the stratosphere is primarily controlled by temperatures in the tropical upper troposphere and lower stratosphere. However, the direct impact of deep convection on the global lower stratospheric water vapor budget is still an actively debated issue. Two complementary modeling approaches are used to investigate the convective impact in boreal winter and summer. Convective influence is diagnosed by tracing trajectories through convective cloud top altitude fields derived from global rainfall and brightness temperature data. Backward trajectory model simulations coupled with a detailed treatment of cloud microphysical processes indicate that convection moistens the global lower stratosphere by approximately 0.3 ppmv in boreal winter and summer 2010. The diurnal peak in convection is responsible for about half of the total convective moistening during winter and nearly all of the convective moistening during summer. Deep convective clouds overshooting the tropopause have relatively minor effect on global lower stratospheric water vapor. A forward trajectory model coupled with a simplified cloud module is used to estimate the relative magnitude of the interannual variability of the convective impact. Combining the results from the two models, we find that the convective impact on the global lower stratospheric water vapor during 2006-2016 is approximately 0.3 ppmv with year-to-year variations of up to 0.1 ppmv. An important mechanism of convective hydration of the lower stratosphere is via the detrainment of saturated air and ice into the tropical uppermost troposphere and the subsequent upward transport of some of these moist air parcels across relatively warm and subsaturated tropopause.

CONVECTION↗

NASA NTRS1982

Modelling of auroral electrodynamical processes: Magnetosphere to mesosphere

Research conducted on auroral electrodynamic coupling between the magnetosphere and ionosphere-atmosphere in support of the development of a global scale kinetic plasma theory is reviewed. Topics covered include electric potential structure in the evening sector; morning and dayside auroras; auroral plasma formation; electrodynamic coupling with the thermosphere; and auroral electron interaction with the atmosphere.

Chiu, Y. T.↗

NASA NTRS1972

Magnetospheric electrons

Coupling of source, transport, and sink processes produces a fairly accurate model for the macroscopic structure and dynamics of magnetospheric electrons. Auroral electrons are controlled by convective transport from a plasma sheet source coupled with a precipitation loss due to whistler and electrostatic plasma turbulence. Outer and inner zone electrons are governed by radial diffusion transport from convection and acceleration sources external to the plasmapause and by parasitic precipitation losses arising from cyclotron and Landau interactions with whistler and ion cyclotron turbulence.

Coroniti, F. V.↗

NASA NTRS2019

The Connected Isotopic Water Cycle in the Community Earth System Model Version 1

Because of the pervasive role of water in the Earth system, the relative abundances of stable isotopologues of water are valuable for understanding atmospheric, oceanic, and biospheric processes, and for interpreting paleoclimate proxy reconstructions. Isotopologues are transported by both largescale and turbulent flows, and the ratio of heavy to light isotopologues changes due to fractionation that can accompany condensation and evaporation processes. Correctly predicting the isotopic distributions requires resolving the relationships between largescale ocean and atmospheric circulation and smallerscale hydrological processes, which can be accomplished within a coupled climate modeling framework. Here we present the water isotopeenabled version of the Community Earth System Model version 1 (iCESM1), which simulates global variations in water isotopic ratios in the atmosphere, land, ocean, and sea ice. In a transient Last Millennium simulation covering the 850-2005 period, iCESM1 correctly captures the latetwentiethcentury structure of δ(exp 18)O and δD over the global oceans, with more limited accuracy over land. The relationship between salinity and seawater δ(exp 18)O is also well represented over the observational period, including interbasin variations. We illustrate the utility of coupled, isotopeenabled simulations using both Last Millennium simulations and freshwater hosing experiments with iCESM1. Closing the isotopic mass balance between all components of the coupled model provides new confidence in the underlying depiction of the water cycle in CESM, while also highlighting areas where the underlying hydrologic balance can be improved. The iCESM1 is poised to be a vital community resource for ongoing model development with both modern and paleoclimate applications.

Brady, E.↗

NASA NTRSDate not supplied

Improved Treatment of Snow Over Ice Sheets in the NASA GISS Climate Model: Towards Ice Sheet–Climate Coupling

Representing the interactions between ice sheets and climate is essential for more accurate prediction of climate change and sea level rise. Ice sheets interact with the overlying atmosphere via the accumulation of snow and its compaction into firn, then ice, as well as the melting of surface snow and ice and the creation of runoff water. Getting an adequate representation of heat transfer, compaction, and melting processes is essential for an accurate representation of snow on land ice in global climate models. We are implementing an improved snow model on top of land ice as part of an effort to couple the NASA GISS climate model with the PISM ice sheet model. The new snow model includes additional layers and processes that are not currently incorporated (e.g., liquid water retention, percolation and refreezing, and snow densification), and mass and energy transfer methods that are consistent with both static ice sheets (with implicit iceberg fluxes) and interactive ice sheets (with explicit dynamics). We are tuning the densification scheme of this snow model with temperature and density data from common FirnCover and SumUp observations at locations in the accumulation zone of Greenland, and we compare the resulting density profiles to other SumUp density profiles in Greenland and Antarctica. We will assess the impact of this new snow model in climate model simulations with a static ice sheet compared with the previous (simpler) 2-layer snow model. Finally, we aim to use the non-coupled simulations as a baseline to assess the impact of dynamic coupling with an interactive ice sheet model.

snow↗

NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

NASA NTRS2013

The NASA-Goddard Multi-Scale Modeling Framework - Land Information System: Global Land/atmosphere Interaction with Resolved Convection

The present generation of general circulation models (GCM) use parameterized cumulus schemes and run at hydrostatic grid resolutions. To improve the representation of cloud-scale moist processes and landeatmosphere interactions, a global, Multi-scale Modeling Framework (MMF) coupled to the Land Information System (LIS) has been developed at NASA-Goddard Space Flight Center. The MMFeLIS has three components, a finite-volume (fv) GCM (Goddard Earth Observing System Ver. 4, GEOS-4), a 2D cloud-resolving model (Goddard Cumulus Ensemble, GCE), and the LIS, representing the large-scale atmospheric circulation, cloud processes, and land surface processes, respectively. The non-hydrostatic GCE model replaces the single-column cumulus parameterization of fvGCM. The model grid is composed of an array of fvGCM gridcells each with a series of embedded GCE models. A horizontal coupling strategy, GCE4fvGCM4Coupler4LIS, offered significant computational efficiency, with the scalability and I/O capabilities of LIS permitting landeatmosphere interactions at cloud-scale. Global simulations of 2007e2008 and comparisons to observations and reanalysis products were conducted. Using two different versions of the same land surface model but the same initial conditions, divergence in regional, synoptic-scale surface pressure patterns emerged within two weeks. The sensitivity of largescale circulations to land surface model physics revealed significant functional value to using a scalable, multi-model land surface modeling system in global weather and climate prediction.

cumulus schemes↗

NASA NTRS1977

A three-D circulation model study of the radiative-dynamic coupling within the stratosphere

The role of radiative transfer processes in the general circulation of the stratosphere was studied by comparing the thermal and dynamical structure of the stratosphere as simulated by two controlled numerical experiments performed with the aid of a spectral 3-dimensional quasi-geostrophic circulation model. In the first experiment, a detailed radiative transfer model is employed to treat the longwave radiative transfer in the stratosphere. In the second experiment, longwave radiative transfer is accounted for by employing the Newtonian cooling approximation. It is found that the exchange of longwave radiation between the troposphere and lower stratosphere has a net heating effect on the lower stratosphere. It is shown that this heating effect contributes partly to the maintenance of the warm high-latitude belt in the lower stratosphere during winter and spring seasons. The strong temperature dependence of the Newtonian cooling coefficient plays an important role in determining the zonal temperatures and has a significant influence on the transmissivity of stratosphere to propagating planetary scale waves.

Ramanathan, V.↗

NASA NTRS1998

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗

NASA NTRS2006

Phosphate Reactions as Mechanisms of High-Temperature Lubrication

One of the major problems preventing the operation of advanced gas turbine engines at higher temperatures is the inability of currently used liquid lubricants to survive at these higher temperatures under friction and wear conditions. Current state-of-the-art organic liquid lubricants rapidly degrade at temperatures above 300 C; hence some other form of lubrication is necessary. Vapor-phase lubrication is a promising new technology for high-temperature lubrication. This lubrication method employs a liquid phosphate ester that is vaporized and delivered to bearings or gears; the vapor reacts with the metal surfaces, generating a solid lubricious film that has proven very stable at high temperatures. In this study, solid lubricious films were grown on cast-iron foils in order to obtain reaction and diffusion rate data to help characterize the growth mechanism. A phenomenological mathematical model of the film deposition process was derived incorporating transport and kinetic parameters that were coupled to the experimental data. This phenomenological model can now be reliably used as a predictive and scale-up tool for future vapor-phase lubrication studies.

Nagarajan, Anitha↗

NASA NTRS2011

Physical Validation of GPM Retrieval Algorithms Over Land: An Overview of the Mid-Latitude Continental Convective Clouds Experiment (MC3E)

The joint NASA Global Precipitation Measurement (GPM) -- DOE Atmospheric Radiation Measurement (ARM) Midlatitude Continental Convective Clouds Experiment (MC3E) was conducted from April 22-June 6, 2011, centered on the DOE-ARM Southern Great Plains Central Facility site in northern Oklahoma. GPM field campaign objectives focused on the collection of airborne and ground-based measurements of warm-season continental precipitation processes to support refinement of GPM retrieval algorithm physics over land, and to improve the fidelity of coupled cloud resolving and land-surface satellite simulator models. DOE ARM objectives were synergistically focused on relating observations of cloud microphysics and the surrounding environment to feedbacks on convective system dynamics, an effort driven by the need to better represent those interactions in numerical modeling frameworks. More specific topics addressed by MC3E include ice processes and ice characteristics as coupled to precipitation at the surface and radiometer signals measured in space, the correlation properties of rainfall and drop size distributions and impacts on dual-frequency radar retrieval algorithms, the transition of cloud water to rain water (e.g., autoconversion processes) and the vertical distribution of cloud water in precipitating clouds, and vertical draft structure statistics in cumulus convection. The MC3E observational strategy relied on NASA ER-2 high-altitude airborne multi-frequency radar (HIWRAP Ka-Ku band) and radiometer (AMPR, CoSMIR; 10-183 GHz) sampling (a GPM "proxy") over an atmospheric column being simultaneously profiled in situ by the University of North Dakota Citation microphysics aircraft, an array of ground-based multi-frequency scanning polarimetric radars (DOE Ka-W, X and C-band; NASA D3R Ka-Ku and NPOL S-bands) and wind-profilers (S/UHF bands), supported by a dense network of over 20 disdrometers and rain gauges, all nested in the coverage of a six-station mesoscale rawinsonde network. As an exploratory effort to examine land-surface emissivity impacts on retrieval algorithms, and to demonstrate airborne soil moisture retrieval capabilities, the University of Tennessee Space Institute Piper aircraft carrying the MAPIR L-band radiometer was also flown during the latter half of the experiment in coordination with the ER-2. The observational strategy provided a means to sample the atmospheric column in a redundant framework that enables inter-calibration and constraint of measured and retrieved precipitation characteristics such as particle size distributions, or water contents- all within the umbrella of "proxy" satellite measurements (i.e., the ER-2). Complimenting the precipitation sampling framework, frequent and coincident launches of atmospheric soundings (e.g., 4-8/day) then provided a much larger mesoscale view of the thermodynamic and winds environment, a data set useful for initializing cloud models. The datasets collected represent a variety cloud and precipitation types including isolated cumulus clouds, severe thunderstorms, mesoscale convective systems, and widespread regions of light to moderate stratiform precipitation. We will present the MC3E experiment design, an overview of operations, and a summary of preliminary results.

Petersen, Walter A.↗

NASA NTRS1979

A coupled radiative-convective-photochemical model of the stratosphere

A one-dimensional atmospheric model is described. The neutral photochemical processes of the troposphere and stratosphere are presented by a reaction set of 87 kinetic processes and 28 photolytic reactions, with vertical transport being taken into account by an eddy diffusion parameter. The vertical temperature profiles are coupled with the chemical calculations by means of a radiative-convective code. The model has the capability of including effects of multiple scattering and diurnal variations in species concentrations. The numerical solution methodology is described in detail, and current input data are given along with results obtained in a sample calculation.

Boughner, R. E.↗

NASA NTRS1997

Coupling Processes Between Atmospheric Chemistry and Climate

This is the first semi-annual report for NAS5-97039 summarizing work performed for January 1997 through June 1997. Work in this project is related to NAS1-20666, also funded by NASA ACMAP. The work funded in this project also benefits from work at AER associated with the AER three-dimensional isentropic transport model funded by NASA AEAP and the AER two-dimensional climate-chemistry model (co-funded by Department of Energy). The overall objective of this project is to improve the understanding of coupling processes between atmospheric chemistry and climate. Model predictions of the future distributions of trace gases in the atmosphere constitute an important component of the input necessary for quantitative assessments of global change. We will concentrate on the changes in ozone and stratospheric sulfate aerosol, with emphasis on how ozone in the lower stratosphere would respond to natural or anthropogenic changes. The key modeling tools for this work are the AER two-dimensional chemistry-transport model, the AER two-dimensional stratospheric sulfate model, and the AER three-wave interactive model with full chemistry.

Ko, M. K. W.↗