Search NASA⌕ Search

SEARCH · Search NASA

Results for “submodeling”

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 235 records · Page 13

System Modeling of Metabolic Heat Regenerated Temperature Swing Adsorption (MTSA) Subassembly for Prototype Design

This paper describes modeling methods for the three core components of a Metabolic heat regenerated Temperature Swing Adsorption (MTSA) subassembly: the sorbent bed, a sublimation (cooling) heat exchanger (SHX), and a condensing icing (warming) heat exchanger (CIHX). The primary function of the MTSA, removing carbon dioxide from a ventilation loop, is performed via the sorbent bed. The CIHX is used to heat the sorbent bed for desorption and to remove moisture from the ventilation loop while the SHX is alternately employed to cool the sorbent bed via sublimation of a spray of water at low pressure to prepare the reconditioned bed for the next cycle. This paper describes a system level model of the MTSA as developed in Thermal Desktop and SINDA/FLUINT including assumptions on geometry and physical phenomena, modeling methodology and relevant pa ra mete rizatio ns. Several areas of particular modeling interest are discussed. In the sorbent bed, capture of the translating CO2 saturation front and associated local energy and mass balance in both adsorbing and desorbing modes is covered. The CIHX poses particular challenges for modeling in SINDA/FLUINT as accounting for solids states in fluid submodels are not a native capability. Methods for capturing phase change and latent heat of ice as well as the transport properties across a layer of low density accreted frost are developed. This extended modeling capacity is applicable to temperatures greater than 258 K. To extend applicability to the minimum device temperature of 235 K, a method for a mapped transformation of temperatures from below the limit temperatures to some value above is given along with descriptions for associated material property transformations and the resulting impacts to total heat and mass transfer. Similar considerations are shown for the SHX along with assumptions for flow mechanics and resulting model methods for sublimation in a flow.

Bower, Chad↗

Short-Term Forecasting of Radiation Belt and Ring Current

A computer program implements a mathematical model of the radiation-belt and ring-current plasmas resulting from interactions between the solar wind and the Earth s magnetic field, for the purpose of predicting fluxes of energetic electrons (10 keV to 5 MeV) and protons (10 keV to 1 MeV), which are hazardous to humans and spacecraft. Given solar-wind and interplanetary-magnetic-field data as inputs, the program solves the convection-diffusion equations of plasma distribution functions in the range of 2 to 10 Earth radii. Phenomena represented in the model include particle drifts resulting from the gradient and curvature of the magnetic field; electric fields associated with the rotation of the Earth, convection, and temporal variation of the magnetic field; and losses along particle-drift paths. The model can readily accommodate new magnetic- and electric-field submodels and new information regarding physical processes that drive the radiation-belt and ring-current plasmas. Despite the complexity of the model, the program can be run in real time on ordinary computers. At present, the program can calculate present electron and proton fluxes; after further development, it should be able to predict the fluxes 24 hours in advance

Fok, Mei-Ching↗

Combustion Stability Analyses for J-2X Gas Generator Development

The National Aeronautics and Space Administration (NASA) is developing a liquid oxygen/liquid hydrogen rocket engine for upper stage and trans-lunar applications of the Ares vehicles for the Constellation program. This engine, designated the J-2X, is a higher pressure, higher thrust variant of the Apollo-era J-2 engine. Development was contracted to Pratt & Whitney Rocketdyne in 2006. Over the past several years, development of the gas generator for the J-2X engine has progressed through a variety of workhorse injector, chamber, and feed system configurations. Several of these configurations have resulted in injection-coupled combustion instability of the gas generator assembly at the first longitudinal mode of the combustion chamber. In this paper, the longitudinal mode combustion instabilities observed on the workhorse test stand are discussed in detail. Aspects of this combustion instability have been modeled at the NASA Marshall Space Flight Center with several codes, including the Rocket Combustor Interaction Design and Analysis (ROCCID) code and a new lumped-parameter MatLab model. To accurately predict the instability characteristics of all the chamber and injector geometries and test conditions, several features of the submodels in the ROCCID suite of calculations required modification. Finite-element analyses were conducted of several complicated combustion chamber geometries to determine how to model and anchor the chamber response in ROCCID. A large suite of sensitivity calculations were conducted to determine how to model and anchor the injector response in ROCCID. These modifications and their ramification for future stability analyses of this type are discussed in detail. The lumped-parameter MatLab model of the gas generator assembly was created as an alternative calculation to the ROCCID methodology. This paper also describes this model and the stability calculations.

Hulka, J. R.↗

Modeling Of Metabolic Heat Regenerated Temperature Swing Adsorption (MTSA) Subassembly For Prototype Design

This paper describes modeling methods for the three core components of a Metabolic heat regenerated Temperature Swing Adsorption (MTSA) subassembly: a sorbent bed, a sublimation (cooling) heat exchanger (SHX), and a condensing icing (warming) heat exchanger (CIHX). The primary function of the MTSA, removing carbon dioxide from a space suit Portable Life Support System (PLSS) ventilation loop, is performed via the sorbent bed. The CIHX is used to heat the sorbent bed for desorption and to remove moisture from the ventilation loop while the SHX is alternately employed to cool the sorbent bed via sublimation of a spray of water at low pressure to prepare the reconditioned bed for the next cycle. This paper describes subsystem heat a mass transfer modeling methodologies relevant to the description of the MTSA subassembly in Thermal Desktop and SINDA/FLUINT. Several areas of particular modeling interest are discussed. In the sorbent bed, capture of the translating carbon dioxide (CO2) front and associated local energy and mass balance in both adsorbing and desorbing modes is covered. The CIHX poses particular challenges for modeling in SINDA/FLUINT as accounting for solids states in fluid submodels are not a native capability. Methods for capturing phase change and latent heat of ice as well as the transport properties across a layer of low density accreted frost are developed. This extended modeling capacity is applicable to temperatures greater than 258 K. To extend applicability to the minimum device temperature of 235 K, a method for a mapped transformation of temperatures from below the limit temperatures to some value above is given along with descriptions for associated material property transformations and the resulting impacts to total heat and mass transfer. Similar considerations are given for the SHX along with functional relationships for areal sublimation rates as limited by flow mechanics in t1he outlet duct.

Bower, Chad E.↗

The Spatial Standard Observer

The spatial standard observer is a computational model that provides a measure of the visibility of a target in a uniform background image or of the visual discriminability of two images. Standard observers have long been used in science and industry to quantify the discriminability of colors. Color standard observers address the spectral characteristics of visual stimuli, while the spatial standard observer (SSO), as its name indicates, addresses spatial characteristics. The SSO is based on a model of human vision. The SSO was developed in a process that included evaluation of a number of earlier mathematical models that address optical, physiological, and psychophysical aspects of spatial characteristics of human visual perception. Elements of the prior models are incorporated into the SSO, which is formulated as a compromise between accuracy and simplicity. The SSO operates on a digitized monochrome still image or on a pair of such images. The SSO consists of three submodels that operate sequentially on the input image(s): 1. A contrast model, which converts an input monochrome image to a luminance contrast image, wherein luminance values are expressed as excursions from, and normalized to, a mean; 2. A contrast-sensitivity-filter model that includes an oblique-effect filter (which accounts for the decline in contrast sensitivity at oblique viewing angles); and 3. A spatial summation model, in which responses are spatially pooled by raising each pixel to the power beta, adding the results, and raising the sum to the 1/b power. In this model, b=2.9 was found to be a suitable value. The net effect of the SSO is to compute a numerical measure of the perceptual strength of the single image, or of the visible difference (denoted the perceptual distance) between two images. The unit of a measure used in the SSO is the just noticeable difference (JND), which is a standard measure of perceptual discriminability. A target that is just visible has a measure of 1 JND.

Watson, Andrew B.↗

Partial Model of Insulator/Insulator Contact Charging

Two papers present a two-phase equilibrium model that partly explains insulator/ insulator contact charging. In this model, a vapor of ions within a gas is in equilibrium with a submonolayer of ions of the same species that have been adsorbed on the surface of an insulator. The surface is modeled as having localized states, each with a certain energy of adsorption for an ion. In an earlier version of the model described in the first paper, the ions do not interact with each other. Using the grand canonical ensemble, the chemical potentials of both vapor and absorbed phases are derived and equated to determine the vapor pressure. If a charge is assigned to the vapor particles (in particular, if single ionization is assumed), then the surface charge density associated with adsorbed ions can be calculated as a function of pressure. In a later version of the model presented in the second paper, the submodel of the vapor phase is extended to include electrostatic interactions between vapor ions and adsorbed ones as well as the screening effect, at a given distance from the surface, of ions closer to the surface. Theoretical values of this model closely match preliminary experimental data on the discharge of insulators as a function of pressure.

Hogue, Michael↗

Analyses of Longitudinal Mode Combustion Instability in J-2X Gas Generator Development

The National Aeronautics and Space Administration (NASA) and Pratt & Whitney Rocketdyne are developing a liquid oxygen/liquid hydrogen rocket engine for future upper stage and trans-lunar applications. This engine, designated the J-2X, is a higher pressure, higher thrust variant of the Apollo-era J-2 engine. The contract for development was let to Pratt & Whitney Rocketdyne in 2006. Over the past several years, development of the gas generator for the J-2X engine has progressed through a variety of workhorse injector, chamber, and feed system configurations on the component test stand at the NASA Marshall Space Flight Center (MSFC). Several of the initial configurations resulted in combustion instability of the workhorse gas generator assembly at a frequency near the first longitudinal mode of the combustion chamber. In this paper, several aspects of these combustion instabilities are discussed, including injector, combustion chamber, feed system, and nozzle influences. To ensure elimination of the instabilities at the engine level, and to understand the stability margin, the gas generator system has been modeled at the NASA MSFC with two techniques, the Rocket Combustor Interaction Design and Analysis (ROCCID) code and a lumped-parameter MATLAB(TradeMark) model created as an alternative calculation to the ROCCID methodology. To correctly predict the instability characteristics of all the chamber and injector geometries and test conditions as a whole, several inputs to the submodels in ROCCID and the MATLAB(TradeMark) model were modified. Extensive sensitivity calculations were conducted to determine how to model and anchor a lumped-parameter injector response, and finite-element and acoustic analyses were conducted on several complicated combustion chamber geometries to determine how to model and anchor the chamber response. These modifications and their ramification for future stability analyses of this type are discussed.

Hulka, J. R.↗

Open-Source Software for Modeling of Nanoelectronic Devices

The Nanoelectronic Modeling 3-D (NEMO 3-D) computer program has been upgraded to open-source status through elimination of license-restricted components. The present version functions equivalently to the version reported in "Software for Numerical Modeling of Nanoelectronic Devices" (NPO-30520), NASA Tech Briefs, Vol. 27, No. 11 (November 2003), page 37. To recapitulate: NEMO 3-D performs numerical modeling of the electronic transport and structural properties of a semiconductor device that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantum-mechanical behavior of the device resolved to the atomistic level of granularity. NEMO 3-D solves the applicable quantum matrix equation on a Beowulf-class cluster computer by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. A prior upgrade of NEMO 3-D incorporated a capability for a strain treatment, parameterized for bulk material properties of GaAs and InAs, for two tight-binding submodels. NEMO 3-D has been demonstrated in atomistic analyses of effects of disorder in alloys and, in particular, in bulk In(x)Ga(1-x)As and in In(0.6)Ga(0.4)As quantum dots.

Oyafuso, Fabiano↗

Numerical Modeling of Nanoelectronic Devices

Nanoelectronic Modeling 3-D (NEMO 3-D) is a computer program for numerical modeling of the electronic structure properties of a semiconductor device that is embodied in a crystal containing as many as 16 million atoms in an arbitrary configuration and that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantummechanical behavior of the device resolved to the atomistic level of granularity. The system of electrons in the device is represented by a sparse Hamiltonian matrix that contains hundreds of millions of terms. NEMO 3-D solves the matrix equation on a Beowulf-class cluster computer, by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. In a recent update of NEMO 3-D, a new strain treatment, parameterized for bulk material properties of GaAs and InAs, was developed for two tight-binding submodels. The utility of the NEMO 3-D was demonstrated in an atomistic analysis of the effects of disorder in alloys and, in particular, in bulk In(x)Ga(l-x)As and in In0.6Ga0.4As quantum dots.

Klimeck, Gerhard↗

An Integrated Framework for Model-Based Distributed Diagnosis and Prognosis

Diagnosis and prognosis are necessary tasks for system reconfiguration and fault-adaptive control in complex systems. Diagnosis consists of detection, isolation and identification of faults, while prognosis consists of prediction of the remaining useful life of systems. This paper presents a novel integrated framework for model-based distributed diagnosis and prognosis, where system decomposition is used to enable the diagnosis and prognosis tasks to be performed in a distributed way. We show how different submodels can be automatically constructed to solve the local diagnosis and prognosis problems. We illustrate our approach using a simulated four-wheeled rover for different fault scenarios. Our experiments show that our approach correctly performs distributed fault diagnosis and prognosis in an efficient and robust manner.

Bregon, Anibal↗

Next Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models possible.

Arnold, William↗

Next-Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models possible

Arnold, William R., Sr.↗

Next Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models easier.

Arnold, William R., Sr.↗

Using SpaceClaimTD Direct for Modeling Components with Complex Geometries for the Thermal Desktop-Based Advanced Stirling Radioisotope Generator Model

SpaceClaim/TD Direct has been used extensively in the development of the Advanced Stirling Radioisotope Generator (ASRG) thermal model. This paper outlines the workflow for that aspect of the task and includes proposed best practices and lessons learned. The ASRG thermal model was developed to predict component temperatures and power output and to provide insight into the prime contractor's thermal modeling efforts. The insulation blocks, heat collectors, and cold side adapter flanges (CSAFs) were modeled with this approach. The model was constructed using mostly TD finite difference (FD) surfaces/solids. However, some complex geometry could not be reproduced with TD primitives while maintaining the desired degree of geometric fidelity. Using SpaceClaim permitted the import of original CAD files and enabled the defeaturing/repair of those geometries. TD Direct (a SpaceClaim add-on from CRTech) adds features that allowed the "mark-up" of that geometry. These so-called "mark-ups" control how finite element (FE) meshes are to be generated through the "tagging" of features (e.g. edges, solids, surfaces). These tags represent parameters that include: submodels, material properties, material orienters, optical properties, and radiation analysis groups. TD aliases were used for most tags to allow analysis to be performed with a variety of parameter values. "Domain-tags" were also attached to individual and groups of surfaces and solids to allow them to be used later within TD to populate objects like, for example, heaters and contactors. These tools allow the user to make changes to the geometry in SpaceClaim and then easily synchronize the mesh in TD without having to redefine the objects each time as one would if using TDMesher. The use of SpaceClaim/TD Direct helps simplify the process for importing existing geometries and in the creation of high fidelity FE meshes to represent complex parts. It also saves time and effort in the subsequent analysis.

Power Systems↗

Using SpaceClaim/TD Direct for Modeling Components with Complex Geometries for the Thermal Desktop-Based Advanced Stirling Radioisotope Generator Model

SpaceClaim/TD Direct has been used extensively in the development of the Advanced Stirling Radioisotope Generator (ASRG) thermal model. This paper outlines the workflow for that aspect of the task and includes proposed best practices and lessons learned. The ASRG thermal model was developed to predict component temperatures and power output and to provide insight into the prime contractors thermal modeling efforts. The insulation blocks, heat collectors, and cold side adapter flanges (CSAFs) were modeled with this approach. The model was constructed using mostly TD finite difference (FD) surfaces solids. However, some complex geometry could not be reproduced with TD primitives while maintaining the desired degree of geometric fidelity. Using SpaceClaim permitted the import of original CAD files and enabled the defeaturing repair of those geometries. TD Direct (a SpaceClaim add-on from CRTech) adds features that allowed the mark-up of that geometry. These so-called mark-ups control how finite element (FE) meshes were generated and allowed the tagging of features (e.g. edges, solids, surfaces). These tags represent parameters that include: submodels, material properties, material orienters, optical properties, and radiation analysis groups. TD aliases were used for most tags to allow analysis to be performed with a variety of parameter values. Domain-tags were also attached to individual and groups of surfaces and solids to allow them to be used later within TD to populate objects like, for example, heaters and contactors. These tools allow the user to make changes to the geometry in SpaceClaim and then easily synchronize the mesh in TD without having to redefine these objects each time as one would if using TD Mesher.The use of SpaceClaim/TD Direct has helped simplify the process for importing existing geometries and in the creation of high fidelity FE meshes to represent complex parts. It has also saved time and effort in the subsequent analysis.

Power Systems↗

Thermal Modeling Method Improvements for SAGE III on ISS

The Stratospheric Aerosol and Gas Experiment III (SAGE III) instrument is the fifth in a series of instruments developed for monitoring aerosols and gaseous constituents in the stratosphere and troposphere. SAGE III will be delivered to the International Space Station (ISS) via the SpaceX Dragon vehicle. A detailed thermal model of the SAGE III payload, which consists of multiple subsystems, has been developed in Thermal Desktop (TD). Many innovative analysis methods have been used in developing this model; these will be described in the paper. This paper builds on a paper presented at TFAWS 2013, which described some of the initial developments of efficient methods for SAGE III. The current paper describes additional improvements that have been made since that time. To expedite the correlation of the model to thermal vacuum (TVAC) testing, the chambers and GSE for both TVAC chambers at Langley used to test the payload were incorporated within the thermal model. This allowed the runs of TVAC predictions and correlations to be run within the flight model, thus eliminating the need for separate models for TVAC. In one TVAC test, radiant lamps were used which necessitated shooting rays from the lamps, and running in both solar and IR wavebands. A new Dragon model was incorporated which entailed a change in orientation; that change was made using an assembly, so that any potential additional new Dragon orbits could be added in the future without modification of the model. The Earth orbit parameters such as albedo and Earth infrared flux were incorporated as time-varying values that change over the course of the orbit; despite being required in one of the ISS documents, this had not been done before by any previous payload. All parameters such as initial temperature, heater voltage, and location of the payload are defined based on the case definition. For one component, testing was performed in both air and vacuum; incorporating the air convection in a submodel that was only built for the in-air cases allowed correlation of all testing to be done in a single model. These modeling improvements and more will be described and illustrated in the paper.

Liles, Kaitlin↗

Atmospheric Aerosol and Thermal Structure in the Boundary Layer Over the Los Angeles Basin

A field study using a mobile lidar was recently conducted in the L. A. Basin, California, to (1) examine the relationship between the vertical aerosol and the thermal structure, and (2) map the vertical aerosol structure in the atmospheric boundary layer over the basin. These data are needed for use in the development of a mixing-depth submodel required for photochemical air Quality simulation models. Toward these ends, a series of lidar aerosol measurements in conjunction with balloon and aircraft temperature soundings were taken at a site in El Monte, and in a mobile mode along a 90-mile freeway loop between El Monte, Santa Monica, and Long Beach. The lidar data are presented in the form of time-height and distance-height cross sections. The results indicate that, although aerosol concentrations are frequently present above the base of the marine inversion, these are generally in stratified layers in contrast to the more uniform nature of the lower convective layer, permitting the mixing depth to be distinguished on this basis. The lidar-derived mixing depths are well correlated (within 100 m) with daytime temperature inversions. Other significant features shown by the lidar data include large Basin-wide mixing-depth variations, waves with amplitudes of 200-300 m and wavelengths of 1000-1500 m on the lower aerosol layer, and apparent aerosol "chimneys" with overrunning in the vicinity of convergence zones.

Johnson, Warren B.↗

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗