Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed model predictive control”

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 667 records · Page 37

Virtual Design Studio (VDS) - Development of an Integrated Computer Simulation Environment for Performance Based Design of Very-Low Energy and High IEQ Buildings

Executive Summary The objective of this study was to develop a “Virtual Design Studio (VDS)”: a software platform for integrated, coordinated and optimized design of green building systems with low energy consumption, high indoor environmental quality (IEQ), and high level of sustainability. This VDS is intended to assist collaborating architects, engineers and project management team members throughout from the early phases to the detailed building design stages. It can be used to plan design tasks and workflow, and evaluate the potential impacts of various green building strategies on the building performance by using the state of the art simulation tools as well as industrial/professional standards and guidelines for green building system design. Engaged in the development of VDS was a multi-disciplinary research team that included architects, engineers, and software developers. Based on the review and analysis of how existing professional practices in building systems design operate, particularly those used in the U.S., Germany and UK, a generic process for performance-based building design, construction and operation was proposed. It distinguishes the whole process into five distinct stages: Assess, Define, Design, Apply, and Monitoring (ADDAM). The current VDS is focused on the first three stages. The VDS considers building design as a multi-dimensional process, involving multiple design teams, design factors, and design stages. The intersection among these three dimensions defines a specific design task in terms of “who”, “what” and “when”. It also considers building design as a multi-objective process that aims to enhance the five aspects of performance for green building systems: site sustainability, materials and resource efficiency, water utilization efficiency, energy efficiency and impacts to the atmospheric environment, and IEQ. The current VDS development has been limited to energy efficiency and IEQ performance, with particular focus on evaluating thermal performance, air quality and lighting environmental quality because of their strong interaction with the energy performance of buildings. The VDS software framework contains four major functions: 1) Design coordination: It enables users to define tasks using the Input-Process-Output flow approach, which specifies the anticipated activities (i.e., the process), required input and output information, and anticipated interactions with other tasks. It also allows task scheduling to define the work flow, and sharing of the design data and information via the internet. 2) Modeling and simulation: It enables users to perform building simulations to predict the energy consumption and IEQ conditions at any of the design stages by using EnergyPlus and a combined heat, air, moisture and pollutant simulation (CHAMPS) model. A method for co-simulation was developed to allow the use of both models at the same time step for the combined energy and indoor air quality analysis. 3) Results visualization: It enables users to display a 3-D geometric design of the building by reading BIM (building information model) file generated by design software such as SketchUp, and the predicted results of heat, air, moisture, pollutant and light distributions in the building. 4) Performance evaluation: It enables the users to compare the performance of a proposed building design against a reference building that is defined for the same type of buildings under the same climate condition, and predicts the percent of improvements over the minimum requirements specified in ASHRAE Standard 55-2010, 62.1-2010 and 90.1-2010. An approach was developed to estimate the potential impact of a design factor on the whole building performance, and hence can assist the user to identify areas that have most pay back for investment. The VDS software was developed by using C++ with the conventional Model, View and Control (MVC) software architecture. The software has been verified by using a simple 3-zone case building. The application of the VDS concepts and framework for building design and performance analysis has been illustrated by using a medium-sized, five story office building that received LEED Platinum Certification from USGBC.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Global particle buildup simulations with gas puff scan: application to WEST discharge

This paper deals with the distribution of sources, transport, and exhaust of particles in a tokamak. Knowledge and understanding of all the physical phenomena involved in the global particle buildup are necessary to study and predict density regimes and subsequently to develop optimized scenarios for tokamak operation in order to control heat and particle exhaust. Neutral particles and their interactions with plasma are central in this perspective. This paper discusses the impact of varying the intensity of particle fueling in 2D transport simulations of a WEST discharge. Simulations are performed with an updated version of SOLEDGE-HDG that allows a more realistic transport of neutrals using a self-consistent diffusive model based on charge exchange and ionization processes. New code capabilities allow the entire WEST poloidal cross section to be simulated in a realistic configuration for both geometry and the range of control parameters. A gas puff scan illustrates the main features of the sheath-limited, high-recycling, and detached regimes, such as the buildup of the temperature gradient and the pressure drop in the scrape-off layer (SOL), the target temperature falling to 1 eV, and the ionization source moving away from the targets, as well as the particle flux rollover. A crude estimate of wall erosion is also provided, showing the respective role of each plasma wall component in each of these regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Residual Stress Distribution in a Hydroformed Advanced High Strength Steel Component: Neutron Diffraction Measurements and Finite Element Simulations

Today’s automotive industry is witnessing increasing applications of advanced high strength steels (AHSS) combined with innovative manufacturing techniques to satisfy fuel economy requirements of stringent environmental regulations. The integration of AHSS in novel automotive structure design has introduced huge advantages in mass reduction while maintaining their structural performances, yet several concerns have been raised for this relatively new family of steels. One of those concerns is their potentially high springback after forming, which can lead to geometrical deviation of the final product from its designed geometry and cause difficulties during assembly. From the perspective of accurate prediction, control and compensation of springback, further understanding on the effect of residual stress in AHSS parts is urged. In this work, the residual stress distribution in a 980GEN3 steel part after hydroforming is investigated via experimental and numerical approaches. A non-destructive neutron diffraction technique is adopted to reveal the residual stress profiles across a 94-degree tube bending radius section of the hydroformed part. The finite element analysis (FEA) is also conducted to simulate the hydroforming and the subsequent unloading processes to predict the residual stress distribution within the same region. The correlation between experimental and simulation results is presented and the effectiveness of the FEA model for residual stress prediction is discussed. Findings in this study set the basis for springback analysis with improved accuracy and reliability.

Huang, Lu↗

Radial evolution of solar wind thermal electron distributions due to expansion and collisions

ISEE 3 electron observations near 1 AU show that the solar wind thermal electron temperature anisotropy, T(parallel)/T(perpendicular), is typically 1.0 to 1.5, with densest distributions most nearly isotropic, but is sometimes much higher when density is low. For a small observational subset characterized by high density and low bulk speed, T(perpendicular) can exceed T(parallel). Based on these and other observations, a simple model is proposed for radial evolution of thermal electrons in a structureless solar wind under the influence of Coulomb collisions and geometric expansion in a spiral interplanetary magnetic field. The model, which satisfactorily explains the 1 AU observations, shows that the evolution of electron temperature and T(parallel)/T(perpendicular) is controlled by plasma density, flow speed, and initial temperature, with density the most important factor. The model makes predictions for radial distances other than 1 AU and for high solar latitudes.

Phillips, J. L.↗

Machine learning based simultaneous control of air handling unit discharge air and condenser water temperatures set-point for minimized cooling energy in an office building

In this study, an artificial intelligence based real-time prediction and control model to optimize condenser water temperature and discharge air temperature (DAT) set-points in water-cooled air handling unit (AHU) system has been developed. EnergyPlus-MATLAB co-simulation has been conducted to analyze the developed model's effectiveness. Here, to develop artificial neural networks (ANN) model, embedded neural network objects in MATLAB was utilized. The developed model could decide an optimal temperature set-points based on outdoor air wet-bulb temperature to reflect the Korean climate context. As a result, the developed ANN prediction model showed the predictive performance of Cv(RMSE) of approximately 21%. Compared to the conventional fixed temperature algorithm, which fixes AHU DAT at 14°C and condenser water temperature at 32°C, the ANN based optimized control showed a 22% total cooling energy reduction. These results show that significant energy savings can be achieved by simultaneously controlling condenser water temperature and AHU DAT set-points considering Korean climatic characteristics using AI technologies such as ANN models.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Incorporating Parallel Computing into the Goddard Earth Observing System Data Assimilation System (GEOS DAS)

Atmospheric data assimilation is a method of combining actual observations with model forecasts to produce a more accurate description of the earth system than the observations or forecast alone can provide. The output of data assimilation, sometimes called the analysis, are regular, gridded datasets of observed and unobserved variables. Analysis plays a key role in numerical weather prediction and is becoming increasingly important for climate research. These applications, and the need for timely validation of scientific enhancements to the data assimilation system pose computational demands that are best met by distributed parallel software. The mission of the NASA Data Assimilation Office (DAO) is to provide datasets for climate research and to support NASA satellite and aircraft missions. The system used to create these datasets is the Goddard Earth Observing System Data Assimilation System (GEOS DAS). The core components of the the GEOS DAS are: the GEOS General Circulation Model (GCM), the Physical-space Statistical Analysis System (PSAS), the Observer, the on-line Quality Control (QC) system, the Coupler (which feeds analysis increments back to the GCM), and an I/O package for processing the large amounts of data the system produces (which will be described in another presentation in this session). The discussion will center on the following issues: the computational complexity for the whole GEOS DAS, assessment of the performance of the individual elements of GEOS DAS, and parallelization strategy for some of the components of the system.

Larson, Jay W.↗

Integration for Airborne Dust Prediction Systems and Vegetation Phenology to Track Pollen for Asthma Alerts in Public Health Decision Support Systems

The objective of the program is to assess the feasibility of combining a dust transport model with MODIS derived phenology to study pollen transport for integration with a public health decision support system. The use of pollen information has specifically be identified as a critical need by the New Mexico State Health department for inclusion in the Environmental Public Health Tracking (EPHT) program. Material and methods: Pollen can be transported great distances. Local observations of plan phenology may be consistent with the timing and source of pollen collected by pollen sampling instruments. The Dust REgional Atmospheric Model (DREAM) is an integrated modeling system designed to accurately describe the dust cycle in the atmosphere. The dust modules of the entire system incorporate the state of the art parameterization of all the major phases of the atmospheric dust life such as production, diffusion, advection, and removal. These modules also include effects of the particles size distribution on aerosol dispersion. The model was modified to use pollen sources instead of dust. Pollen release was estimated based on satellite-derived phenology of key plan species and vegetation communities. The MODIS surface reflectance product (MOD09) provided information on the start of the plant growing season, growth stage, and pollen release. The resulting deterministic model is useful for predicting and simulating pollen emission and downwind concentration to study details of phenology and meteorology and their dependencies. The proposed linkage in this project provided critical information on the location timing and modeled transport of pollen directly to the EPHT> This information is useful to support the centers for disease control and prevention (CDC)'s National EPHT and the state of New Mexico environmental public health decision support for asthma and allergies alerts.

Luvall, J. C.↗

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

residual stress↗

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

processing↗

On the Low Risk of Subsynchronous Resonance (SSR) in Type III Wind Turbines Operating in Grid-Forming Control

We have shown in a previous work that the risk of sub-synchronous resonance (SSR) between a wind power plant with Type III wind turbines and series-compensated transmission lines is low when the wind turbines in the plant are operated in grid-forming mode, instead of the standard grid-following mode. The fundamental mechanism behind the improved damping characteristics is explained in this paper by modeling the positive and negative sequence impedances of Type III wind turbines for GFM operation mode. It is discovered that the GFM control naturally acts against the negative resistance of Type III wind turbines that results from an interaction between the proportional gain of the rotor-side converter current controller and negative slip. The developed sequence impedance models and improved damping behavior are verified using PSCAD simulations of a 2.5-MW Type III GFM wind turbine. The modeling predictions are also supported by experimentally measuring the sequence impedance response of a 2.5-MW Type III wind turbine during operation in GFL and GFM modes.

grid-following turbine↗

LOW-TEMPERATURE SPECTROSCOPY OF THE {sup 12}C{sub 2}H{sub 2} (υ{sub 1} + υ{sub 3}) BAND IN A HELIUM BUFFER GAS

Buffer gas cooling with a {sup 4}He gas is used to perform laser-absorption spectroscopy of the {sup 12}C{sub 2}H{sub 2} (υ{sub 1} + υ{sub 3}) band at cryogenic temperatures. Doppler thermometry is first carried out to extract translational temperatures from the recorded spectra. Then, rotational temperatures down to 20 K are retrieved by fitting the Boltzmann distribution to the relative intensities of several ro-vibrational lines. The potential of our setup to tune the thermal equilibrium between translational and rotational degrees of freedom is also demonstrated. This can be used to reproduce in a controlled way the regime of non-local thermal equilibrium typically encountered in the interstellar medium. The underlying helium-acetylene collisional physics, relevant for modeling planetary atmospheres, is also addressed. In particular, the diffusion time of {sup 12}C{sub 2}H{sub 2} in the buffer cell is measured against the {sup 4}He flux at two separate translational temperatures; the observed behavior is then compared with that predicted by a Monte Carlo simulation, thus providing an estimate for the respective total elastic cross sections: σ{sub el}(100 K) = (4 ± 1) × 10{sup –20} m{sup 2} and σ{sub el}(25 K) = (7 ± 2) × 10{sup –20} m{sup 2}.

79 ASTRONOMY AND ASTROPHYSICS↗

CFD Simulations of a Regenerative Process for Carbon Dioxide Capture in Advanced Gasification Based Power Systems

This project describes the work carried out to prepare a highly reactive and mechanically strong MgO based sorbents and to develop a Population Balance Equations (PBE) approach to describe the evolution of the particle porosity distribution that is linked with Computational Fluid Dynamics (CFD) to perform simulations of the CO2 capture and sorbent regeneration. A large number of MgO-based regenerable sorbents were prepared using low cost and abundant dolomite as the base material. Among various preparation parameters investigated the potassium/magnesium (K/Mg) ratio was identified as the key variable affecting the reactivity and CO2 capacity of the sorbent. The optimum K/Mg ratio is about 0.15. The sorbent formulation HD52-P2 was identified as the “best” sorbent formulation and a large batch (one kg) of the sorbent was prepared for the detailed study. The results of parametric study indicate the optimum carbonation and regeneration temperatures are 360° and 500°C, respectively. The results also indicate that steam has a beneficial effect on the rate of carbonation and regeneration of the sorbent and that the reactivity and capacity of the sorbent decreases in the cycling process (sorbent deactivation). The results indicate that to achieve a high CO2 removal efficiency, the bed of sorbent should be operated at a temperature range of 370-410°C which also favors production of hydrogen through the WGS reaction. To describe the carbonation reaction kinetics of the MgO, the Variable Diffusivity shrinking core Model (VDM) was developed in this project, which was shown to accurately fit the experimental data. An important advantage of this model is that the changes in the sorbent conversion with time can be expressed in an explicit manner, which will significantly reduce the CFD computation time. A Computational Fluid Dynamic/Population Balance Equations (CFD/PBE) model was developed that accounts for the particle (sorbent) porosity distribution and a new version of the method of moments, called Finite size domain Complete set of trial functions Method Of Moments (FCMOM) was used to solve the population balance equations. The PBE model was implemented in a commercial CFD code, Ansys Fluent 13.0. The code was used to test the model in some simple cases and the results were verified against available analytical solution in the literature. Furthermore, the code was used to simulate CO2 capture in a packed-bed and the results were in excellent agreement with the experimental data obtained in the packed bed. The National Energy Laboratory (NETL) Carbon Capture Unit (C2U) design was used in simulate of the hydrodynamics of the cold flow gas/solid system (Clark et al.58). The results indicate that the pressure drop predicted by the model is in good agreement with the experimental data. Furthermore, the model was shown to be able to predict chugging behavior, which was observed during the experiment. The model was used as a base-case for simulations of reactive flow at elevated pressure and temperatures. The results indicate that by controlling the solid circulation rate, up to 70% CO2 removal can be achieved and that the solid hold up in the riser is one of the main factors controlling the extent of CO2 removal. The CFD/PBE simulation model indicates that by using a simulated syngas with a composition of 20% CO2, 20% H2O, 30% CO, and 30% H2, the composition (wet basis) in the reactor outlet corresponded to about 60% CO2 capture with and exit gas containing 65% H2. A preliminary base-case-design was developed for a regenerative MgO-based pre-combustion carbon capture process for a 500 MW IGCC power plant. To minimize the external energy requirement, an extensive heat integration network was developed in Aspen/HYSYS® to produce the steam required in the regenerator and heat integration. In this process, liquid CO2 produced at 50 atm can easily be pumped and sequestered or stored. The preliminary economic analyses indicate that the estimated cost of carbon v capture is in the range of $31-$44/ton, suggesting that a regenerative MgO-Based process can be a viable option for pre-combustion carbon dioxide capture in advanced gasification based power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Probing light quark Yukawa couplings through angularity distributions in Higgs boson decay

We propose to utilize angularity distributions in Higgs boson decay to probe light quark Yukawa couplings at e + e - colliders. Angularities $τ$ a are a class of 2-jet event shapes with variable and tunable sensitivity to the distribution of radiation in hadronic jets in the final state. Using soft-collinear effective theory (SCET), we present a prediction of angularity distributions from Higgs decaying to quark and gluon states at e + e - colliders to NNLL + $\mathcal{O}$(α s ) accuracy. Due to the different color structures in quark and gluon jets, the angularity distributions from H → $q\overline{q}$ and H → gg show different behaviors and can be used to constrain the light quark Yukawa couplings. We show that the upper limit of light quark Yukawa couplings could be probed to the level of ~ 15% of the bottom quark Yukawa coupling in the Standard Model in a conservative analysis window far away from nonperturbative effects and other uncertainties; the limit can be pushed to ≲ 7 – 9% with better control of the nonperturbative effects especially on gluon angularity distributions and/or with multiple angularities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A New Efficient Quantitative Multi-component Phase Field: Lattice Boltzmann Model for Simulating Ti6Al4V Solidified Dendrite Under Forced Flow

Ti6Al4V is a widely used, multi-component alloy in additive manufacturing, during which the fluid flow in the molten pool significantly affects the solidified dendrites. To predict and further control the microstructure, modeling and simulating the microstructure evolution play a critical role. In this study, a newly developed, efficient, quantitative multi-component phase-field (PF) model is coupled with a lattice Boltzmann (LB) model to simulate Ti6Al4V solidified dendrite evolution under fluid flow. The accuracy and convergence behavior of the model is validated by the Gibbs–Thomson relation at the dendrite tip. Single and multiple two-dimensional (2D) equiaxed dendrite evolution cases under forced flow were simulated. Results show that the dendrite pattern is influenced remarkably by the fluid flow. Underlying mechanisms of the asymmetrical evolution are revealed by discussing the interaction among the flow, composition distribution and dendrite morphology, quantitatively. The dendrite kinetics are also derived, which ascertains the relationship between tip velocity and undercooling and inlet velocity and is the foundation for larger-scale simulation. We believe that the coupled quantitative multi-component PF–LB framework employed in this study helps in investigating the solidified dendrite morphology evolution in a deep and quantitate manner.

36 MATERIALS SCIENCE↗

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spatial Trends and Drivers of Bedload and Suspended Sediment Fluxes in Global Rivers

Abstract Bedload is notoriously challenging to measure and model; its dynamics, therefore, remains largely unknown in most fluvial systems worldwide. We present results from a global scale bedload flux model as part of the WBMsed modeling framework that well predict the distribution of water discharge, suspended sediment and bedload. The sensitivity of bedload predictions to river slope, particle size, discharge, river width, and suspended sediment were analyzed, showing the model to be most responsive to spatial dynamics in river discharge and slope. The relationship between bedload and total sediment flux is analyzed globally, and for representative longitudinal river profiles (Amazon, Mississippi, and Lena Rivers). The results show that while the proportion of bedload decreases from headwaters to the coasts, there is considerable variability between basins and along river corridors. The topographic and hydrological longitudinal profiles of rivers are shown to be the key drivers of bedload trends, with fluctuations in slope controlling its more local dynamics. Estimates of water and sediment fluxes to global oceans from 2,067 largest river outlets (draining 67% of the global continental area) are provided. Estimated water discharge at 30,579 km 3 /y corresponds well to past estimates; however, sediment flux is higher. Total global particulate load of 17.8 Gt/y is delivered to global oceans, 14.8 Gt/y as washload, 1.1 Gt/y as bedload, and 2.6 Gt/y as suspended bed material. The largest 25 rivers are predicted to transport more than half of the total sediment flux to global oceans.

Cohen, Sagy↗

Spray formation processes of impinging jet injectors

A study examining impinging liquid jets has been underway to determine physical mechanisms responsible for combustion instabilities in liquid bi-propellant rocket engines. Primary atomization has been identified as an important process. Measurements of atomization length, wave structure, and drop size and velocity distribution were made under various ambient conditions. Test parameters included geometric effects and flow effects. It was observed that pre-impingement jet conditions, specifically whether they were laminar or turbulent, had the major effect on primary atomization. Comparison of the measurements with results from a two dimensional linear aerodynamic stability model of a thinning, viscous sheet were made. Measured turbulent impinging jet characteristics were contrary to model predictions; the structure of waves generated near the point of jet impingement were dependent primarily on jet diameter and independent of jet velocity. It has been postulated that these impact waves are related to pressure and momentum fluctuations near the impingement region and control the eventual disintegration of the liquid sheet into ligaments. Examination of the temporal characteristics of primary atomization (ligament shedding frequency) strongly suggests that the periodic nature of primary atomization is a key process in combustion instability.

Anderson, W. E.↗