Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed 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 559 records · Page 31

Shrubs Strongly Influence Snow Properties in Two Subarctic Watersheds

Understanding changes in snow distribution in permafrost ecosystems is fundamental to predicting their response to future climate change. The expansion of tall shrubs into tundra ecosystems can trap snow and insulate permafrost ecosystems during the winter, but the overall insulation effect is dependent upon many ecosystem properties. To study shrub–snow–ground interactions, small temperature sensors were deployed at two research sites on the Seward Peninsula of Alaska, USA, during the 2019–2020 winter. Snow temperatures were used to extrapolate multiple metrics, including freezing n-factors, the snow insulation effect, snow cover duration, and the length of the snowmelt period. Statistical and spatial analysis showed that shrub patches were a dominant control on all snow metrics. Within shrub patches, average ground temperatures were 2.1°C warmer, snow persisted 50 days longer, snow insulation was double, and a longer, later spring snowmelt period occurred compared to nonshrubby areas. Site-level differences contributed relatively little to variation in snow metrics, indicating that shrub presence is an overarching driver of snow–ground interactions at the locations examined. Shrub expansion, which is anticipated under climate change, will strongly impact future permafrost distribution and Arctic energy, water, and carbon cycles through snow–shrub–ground feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Transition Core Planning and Safety Analyses in Support of LEU Fuel Conversion of the University of Missouri Research Reactor (MURR)

The University of Missouri Research Reactor (MURR®) is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is working with the National Nuclear Security Administration (NNSA) Office of Material Management and Minimization (M3) Reactor Conversion Program to convert from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel. The M3 Reactor Conversion USHPRR Project objectives include the development of LEU fuel element designs that will ensure safe reactor operations and to maintain the existing experimental performance of each facility. The work is being conducted through many inter-related activities being completed by four Project Pillars: Fuel Qualification (FQ), Fuel Fabrication (FF), Reactor Conversion (RC), and Cross Cutting (CC). A new type of LEU fuel based on an alloy of uranium-10 wt% molybdenum (U-10Mo) is expected to allow the conversion of those USHPRR, like MURR, requiring higher density fuels. The very-high-density LEU U-10Mo monolithic fuel is currently undergoing irradiation testing and post-irradiation examination under a planned and documented fuel qualification effort. The FQ Pillar will document fuel property and fuel performance data and qualify the fuel for use in these reactors. The FF Pillar is fabricating fuel for ongoing and future irradiation tests, as well as conducting fabrication demonstrations to validate or update preliminary fabrication assumptions. The FF Pillar is also working to develop and install commercial manufacturing capacity with the U-10Mo monolithic fuel to produce prototypic fuel. Working with the RC Pillar at Argonne, MURR has progressed through a preliminary fuel element design using preliminary data for the proposed monolithic alloy of U-10Mo. Analyses were completed in previous work that found for typical equilibrium operations with the preliminary LEU fuel element design, in conjunction with a power uprate to 12 MW and appropriate changes to the MURR Limiting safety system settings (LSSS), MURR will have adequate margins to safety for steady-state operations and postulated transient accidents and will have experimental performance in key locations that meets or exceeds current operations with HEU fuel. The purpose of this work is to develop a sequence of transition cycles that will enable MURR to transition from operation with the reactor core loaded with fresh LEU fuel elements only to typical equilibrium operations with mixed-burnup cores following conversion while meeting operational requirements on safety and experimental performance. It is expected that the use of fresh LEU fuel at conversion and subsequent low burnup of the LEU fuel elements that will initially be available for use following conversion will result in critical control blade positions that will substantially change the axial power distribution in the core and the neutron flux available in key experimental locations relative to equilibrium LEU operations. Given the constraints of MURR safety margins, operational practices, and production and research, a novel method has been developed to identify a transition sequence that minimizes the time MURR operates atypically compared to the current prototypic cycles using HEU fuel. The proposed transition sequence moves quickly to the same sort of equilibrium cycles for the LEU fuel that have already been evaluated in documented preliminary safety analyses. Although shifting the neutron flux peak to the lower half of the core during initial cycles with LEU at 12 MW reduces the experiment performance in some key locations relative to current HEU operations at 10 MW, all LEU cores provide an average performance that meets or exceeds that of HEU. An LEU cycle is reached that meets or exceeds the level of experimental performance predicted for current HEU and equilibrium LEU operations in more than 450 key locations identified by a reactor specialist at MURR by the 23rd cycle following conversion and that afterwards will enable MURR to consistently meet its experimental performance requirements. The proposed transition sequence only requires the fabrication of 34 fresh LEU elements in the first year of operation and does not exceed the anticipated availability of fresh elements that can be produced by the fuel fabricator. By the third year after conversion, 22 fresh LEU elements will be required each year, which is the same as expected for equilibrium LEU operations and the same as current operations with HEU fuel. The proposed transition sequence thus combines a relatively short time period before equilibrium burnup is achieved, a temporary increase of fuel elements needed annually relative to typical operations that are within the production capabilities of the fuel fabricator, and demonstrates comparable experimental performance of the LEU cores relative to current HEU operations. Further measures may be taken to reduce any initial experimental performance penalty even further, where possible, by repositioning certain experiments to leverage the increased performance in the lower axial experimental positions in the initial cycles following conversion or leaving the experiments in the irradiation facilities longer in order to achieve the required neutron fluence. This analysis may require refinement depending on the experimental facilities in use at the time of conversion. Nonetheless, the results presented here, including the experimental performance, core burnup, and critical control blade positions throughout the transition cycles, show that the proposed transition cycle fuel management patterns are consistent with what is expected and desired for MURR operation with LEU U-10Mo fuel. Detailed core power distributions from the neutronics models were also used to evaluate safety margins during steady-state operations for the selected transition cycles and the equilibrium LEU core. It is shown that there are adequate safety margins for both steady-state operations and postulated accident scenarios. For the steady-state operations with the preliminary LEU fuel element design the analysis predicts at least 2.49 MW margin to the onset of flow instability at the LSSS power of 15 MW. Considering the LSSS power is 125% of full license power, the margin to OFI is sufficient. In addition, the critical heat flux ratio at LSSS power is well above the requirement of CHFR > 2.0 from NUREG-1537 for all considered cases. For postulated transient accidents, the minimum margin to the fuel temperature safety limit is at least 109 °C. In summary, the proposed sequence of core loadings for MURR operations following conversion to LEU fuel and a power uprate to 12 MW provides sufficient safety margins for both steady-state operations and postulated transient accidents during a proposed sequence of transition cycles to equilibrium operations. Analysis has shown that there are some local experimental performance penalties during the initial cycles. Although there are local shifts in the experimental performance, on average all LEU cores at 12 MW have equal or higher performance than HEU at 10 MW. Temporary adjustments are being planned that will produce suitable experimental performance during these cycles. The results indicate that for the equilibrium LEU core the experimental performance exceeds that of current HEU operations in all key locations while also demonstrating sufficient safety margins.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Distributed Load Shedding Application Architecture and Bi-Level Predictive Estimator Algorithm

Increasing penetrations of distributed renewables are decreasing the effectiveness of traditional decentralized under-frequency load shedding (UFLS) schemes. As more distribution circuits begin to back-feed the transmission system, operation of traditional UFLS may exacerbate frequency instability. This paper presents the conceptual framework for a data-rich environment to coordinate UFLS across multiple distribution providers based on the laminar coordination framework in order to ensure optimal adaptive setting of UFLS relays. Communication and control are enabled through a distributed implementation of the IEC 61968-1 Common Information Model message bus structure. In addition to the proposed architecture, a novel adaptive UFLS scheme informed by a bi-level state estimator to create optimal relay setpoints is introduced. Initial simulation results are presented for the IEEE 14-bus test system on scenarios leading to mis-operation of traditional UFLS.

Anderson, Alexander A.↗

Local Strain and Polarization Mapping in Ferrielectric Materials

CuInP 2 S 6 (CIPS) is a van der Waals material that has attracted attention because of its unusual properties. Recently, a combination of density functional theory (DFT) calculations and piezoresponse force microscopy (PFM) showed that CIPS is a uniaxial quadruple-well ferrielectric featuring two polar phases and a total of four polarization states that can be controlled by external strain. In this study, we combine DFT and PFM to investigate the stress-dependent piezoelectric properties of CIPS, which have so far remained unexplored. The two different polarization phases are predicted to differ in their mechanical properties and the stress sensitivity of their piezoelectric constants. This knowledge is applied to the interpretation of ferroelectric domain images, which enables investigation of local strain and stress distributions. The interplay of theory and experiment produces polarization maps and layer spacings which we compare to macroscopic X-ray measurements. We found that the sample contains only the low-polarization phase and that domains of one polarization orientation are strained, whereas domains of the opposite polarization direction are fully relaxed. The described nanoscale imaging methodology is applicable to any material for which the relationship between electromechanical and mechanical characteristics is known, providing insight on structural, mechanical, and electromechanical properties down to ~10 nm length scales.

36 MATERIALS SCIENCE↗

A novel construct for scaling groundwater–river interactions based on machine-guided hydromorphic classification

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. Predictive numerical models are needed to understand responses of river corridors to environmental change and to support sustainable watershed management. We posit that systematic hydromorphic classification provides a scaling construct that facilitates extrapolation of outputs from local-scale mechanistic models to reduced-order models applicable at reach and watershed scales. This in turn offers the potential to improve large-scale predictions of river corridor hydrobiogeochemical processes. Here we present a new machine-guided hydromorphic classification methodology that addresses the key requirements of this objective, and we demonstrate its application to a segment of the Columbia River in the northwestern United States. The resulting hydromorphic classes form spatially coherent and physically interpretable hydromorphic units that exhibit distinct behaviors in terms of distributions of subsurface residence times (a primary control on critical biogeochemical reactions). This approach forms the basis of ongoing research that is evaluating the formulation of reduced-order models and transferability of results to other river reaches and larger scales.

54 ENVIRONMENTAL SCIENCES↗

Toxin-Antitoxin Systems Reflect Community Interactions Through Horizontal Gene Transfer

Bacterial evolution through horizontal gene transfer (HGT) reflects their community interactions. In this way, HGT networks do well at mapping community interactions, but offer little toward controlling them—an important step in the translation of synthetic strains into natural contexts. Toxin–antitoxin (TA) systems serve as ubiquitous and diverse agents of selection; however, their utility is limited by their erratic distribution in hosts. Here we examine the heterogeneous distribution of TAs as a consequence of their mobility. By systematically mapping TA systems across a 10,000 plasmid network, we find HGT communities have unique and predictable TA signatures. We propose these TA signatures arise from plasmid competition and have further potential to signal the degree to which plasmids, hosts, and phage interact. To emphasize these relationships, we construct an HGT network based solely on TA similarity, framing specific selection markers in the broader context of bacterial communities. This work both clarifies the evolution of TA systems and unlocks a common framework for manipulating community interactions through TA compatibility.

59 BASIC BIOLOGICAL SCIENCES↗

Depth-resolved cathodoluminescence and surface photovoltage spectroscopies of gallium vacancies in β-Ga 2 O 3 with neutron irradiation and forming gas anneals

The gallium vacancy is one of the dominant native point defects in β-Ga 2 O 3 , one that, together with its complexes, can have a major effect on free carrier densities and transport in this wide bandgap semiconductor. We used a combination of depth-resolved cathodoluminescence spectroscopy and surface photovoltage spectroscopy to identify the optical and energy-level properties of these defects as well as how their defect densities and spatial distributions vary with neutron irradiation and temperature-dependent-forming gas anneals. These studies reveal optical signatures that align closely with theoretical energy-level predictions. Likewise, our optical techniques reveal variations in these defect densities that are consistent with hydrogen passivation of gallium vacancies as a function of temperature and depth from the free Ga 2 O 3 surface. These techniques can help guide the understanding and control of dominant native point defects in Ga 2 O 3 .

42 ENGINEERING↗

Editorial: Selective Controls on Microbial Energy Metabolisms: From the Microscale to the Macroscale

A mechanistic and predictive understanding of the genotypic and phenotypic controls on microbial element cycling remains a grand challenge in microbiology. While the composition and concentration of carbon sources, electron donors, and electron acceptors are known to influence microbial community structure by selecting for microbial sub-populations with distinct catabolic and respiratory pathways, selective inhibitors and trace nutrient availability modulate the activity of metabolic enzymes. This in turn influences the distribution of microbial sub-populations with distinct metabolic and respiratory traits, which are then structured by complex multi-dimensional environmental gradients that influence the composition, gene content and element cycling activity of microbiomes. While some parameters are known to select for different respiratory activities (e.g., lanthanides as essential nutrients for methanotrophs, molybdate as a specific inhibitor of sulfate reduction, carbon:nitrogen ratio and concentration as a control on the end-products of nitrate respiration), there are others to discover, and demonstrating how selective parameters operate and mediate element cycling across scales requires multi-disciplinary research in both the lab and field.

59 BASIC BIOLOGICAL SCIENCES↗

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

An unsteady lifting surface method for single rotation propellers

The mathematical formulation of a lifting surface method for evaluating the steady and unsteady loads induced on single rotation propellers by blade vibration and inflow distortion is described. The scheme is based on 3-D linearized compressible aerodynamics and presumes that all disturbances are simple harmonic in time. This approximation leads to a direct linear integral relation between the normal velocity on the blade (which is determined from the blade geometry and motion) and the distribution of pressure difference across the blade. This linear relation is discretized by breaking the blade up into subareas (panels) on which the pressure difference is treated as approximately constant, and constraining the normal velocity at one (control) point on each panel. The piece-wise constant loads can then be determined by Gaussian elimination. The resulting blade loads can be used in performance, stability and forced response predictions for the rotor. Mathematical and numerical aspects of the method are examined. A selection of results obtained from the method is presented. The appendices include various details of the derivation that were felt to be secondary to the main development in Section 1.

Williams, Marc H.↗

Machine Learning for Optimized Polarization at Jefferson Lab

Polarized cryo-targets and polarized photon beams are widely used in experiments at Jefferson Lab. Traditional methods for maintaining the optimal polarization involve manual adjustments throughout data taking by human shift takers. This may introduce some level of inconsistency simply due to the wide variety of experience and expertise of the shift takers themselves. Implementing machine learning-based control systems can improve the stability of the polarization without relying on human intervention. The cryo-target polarization is influenced by temperature, microwave energy, the distribution of paramagnetic radicals, as well as operational conditions including the radiation dose. Diamond radiators are used to generate linearly polarized photons from a primary electron beam. The energy spectrum of these photons can drift over time due to changes in the primary electron beam conditions and diamond degradation. As a first step towards automating the continuous optimization and control processes, uncertainty aware surrogate models have been developed to predict the polarization based on historical data. This talk will provide an overview of the use cases and models developed, highlighting the collaboration between data scientists and physicists at Jefferson Lab.

Jeske, Torri [Thomas Jefferson National Accelerato↗

Residual Stress Modeling and Advanced Diffraction Measurements of 347H Steel Weldments

The 347H stainless steel is a primary high-temperature material for many energy and power generation industries. Stress relief cracking (SRC) has been a particular concern in welding of this material. The residual stress induced by welding and its evolution during post-welding heat treatment (PWHT) and subsequent operating and service conditions is one of the primary factors contributing to SRC. The lifetime of welded structure components is also controlled by the precipitation kinetics that accompanies PWHT, stress relaxation process, and long-term aging and complex synergistic factors. Various theories have been proposed in the past to explain SRC. However, a widely accepted approach to predicting the entire damage evolution and the resulting performance reduction is still lacking. This study is to demonstrate a reliable solution for two critical issues that affect the predictions. First, the residual stress distribution obtained from both simulation and neutron diffraction is compared, which increases the accuracy of mechanical analysis simulation model and therefore builds a solid basis for the lifetime prediction model. Second, the generation, evolution, and annihilation of precipitates are monitored by the synchrotron diffraction experiment. Preliminary results demonstrate the critical importance of precipitation kinetics on the residual stress distribution/redistribution during heat treatments.

Yang, Yi↗

Active Structural Acoustic Control of Interior Noise on a Raytheon 1900D

An active structural acoustic control system has been demonstrated on a Raytheon Aircraft Company 1900D turboprop airliner. Both single frequency and multi-frequency control of the blade passage frequency and its harmonics was accomplished. The control algorithm was a variant of the popular filtered-x LMS implemented in the principal component domain. The control system consisted of 21 inertial actuators and 32 microphones. The actuators were mounted to the aircraft's ring frames. The microphones were distributed uniformly throughout the interior at head height, both seated and standing. Actuator locations were selected using a combinatorial search optimization algorithm. The control system achieved a 14 dB noise reduction of the blade passage frequency during single frequency tests. Multi-frequency control of the first 1st, 2nd and 3rd harmonics resulted in 10.2 dB, 3.3 dB and 1.6 dB noise reductions respectively. These results fall short of the predictions which were produced by the optimization algorithm (13.5 dB, 8.6 dB and 6.3 dB). The optimization was based on actuator transfer functions taken on the ground and it is postulated that cabin pressurization at flight altitude was a factor in this discrepancy.

Palumbo, Dan↗

A Technique for Rapidly Deploying a Concentration Gradient with Applications to Microgravity

The latter half of the last century has seen rapid advancements in semiconductor crystal growth powered by the demand for high performance electronics in myriad applications. The reduced gravity environment of space has also been used for crystal growth tests, especially in instances where terrestrial growth has largely been unsuccessful. While reduced gravity crystal growth affords some control of the gravity parameter, many crystals grown in space, to date, have structural flaws believed to result from convective motions during the growth phase. The character of these instabilities is not well understood but is associated with thermal and solutal density variations near the solidification interface in the presence of residual gravity and g-jitter. In order to study these instabilities in a separate, controlled space experiment, a concentration gradient would first have to be artificially established in a timely manner as an initial condition. This is generally difficult to accomplish in a microgravity environment because the momentum of the fluid injected into a test cell tends to swirl around and mix in the absence of a restoring force. The use of magnetic fields to control the motion and position of liquids has received growing interest in recent times. The possibility of using the force exerted by a non-uniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for space applications. This paper describes a technique for quickly establishing a linear or exponential fluid concentration gradient using a magnetic field in place of gravity to stabilize the deployment. Also discussed is a photometric technique for measuring the concentration profile using light attenuation. Results of the ground-based experiments indicate that the concentration distribution is within 3% of the predicted value. Although any range of concentations can be realized, photometric constraints are discussed which impose some limitations on measurements.

Leslie, Fred↗

Formation of the Potential Jump Over the Geomagnetically Quiet Sunlit Polar Cap Region

We consider the formation of a potential drop over the Earth's polar cap duringgeomagnetically quiet daytime. The observed potential drop is primarily defined by the hydrogen,photoelectron, and polar rain fluxes ratios and depends strongly on the energy distribution of thephotoelectron flux. Polar rain is an essential component of the model required for plasma quasineutrality.The potential distribution along the magnetic field line has two regions, with a small, gradual, potentialdrop of 3–4 V and a potential jump. The value of the potential jump depends on the hydrogen ion tophotoelectron flux ratio and is also controlled by polar rain electrons. With quasineutrality required at itsupper boundary, the jump only occurs in the presence of polar rain and its location depends on the polarrain flux. Model predictions compare well with FAST observations presented by Kitamura et al.(2012, https://doi.org/10.1029/2011JA017459).

Khazanov, George V.↗

Manufacturing control systems and logic for prognosis of defects in composite materials

Presented are manufacturing control systems for composite-material structures, methods for assembling/operating such systems, and transfer molding techniques for predicting and ameliorating void conditions in fiber-reinforced polymer panels. A method for forming a composite-material construction includes receiving a start signal indicating a fiber-based preform is inside a mold cavity, and transmitting a command signal to inject pressurized resin into the mold to induce resin flow within the mold cavity and impregnate the fiber-based preform. An electronic controller receives, from a distributed array of sensors attached to the mold, signals indicative of pressure and/or temperature at discrete locations on an interior face of the mold cavity. The controller determines a measurement deviation between a calibrated baseline value and the pressure and/or temperature values for each of the discrete locations. If any one of the respective measurement deviations exceeds a calibrated threshold, a void signal is generated to flag a detected void condition.

36 MATERIALS SCIENCE↗

Data from: "Lab and Field Warming Similarly Advance Germination Date and Limit Germination Rate for High and Low Elevation Provenances of Two Widespread Subalpine Conifers"

This archive contains data used to draw conclusions in “Lab and Field Warming Similarly Advance Germination Date and Limit Germination Rate for High and Low Elevation Provenances of Two Widespread Subalpine Conifers”, by Kueppers et al. 2017. Cone collection and field experiments were conducted on Niwot Ridge, Colorado, with the latter completed within the Alpine Treeline Warming Experiment (ATWE). Laboratory germination experiments were conducted at the University of California Merced. Geospatial files are also provided for additional geographical context.There are four data file formats in this package: comma-separated values (.csv), Microsoft Excel (.xlsx), keyhole markup language (.kml), and ESRI shapefile (.shp). The latter two are geospatial file formats. .csv files can be opened using any simple text editor program such as TextEdit (Mac) and Notepad (Windows), and .xlsx files can be opened using Microsoft Excel or Google Sheets. The .kml file contains coordinate data, and is compatible with Google Earth and Google Maps. ESRI shapefiles (.shp) are compatible with any geospatial software able to read the file type, such as QGIS and ESRI’s ArcGIS Suite. This data user’s guide is available in as a .pdf, and can be opened using software such as Adobe Acrobat Reader.-----------------------------------------------------------------------------------------------------------------------------------Accurately predicting upslope shifts in subalpine tree ranges with warming requires understanding how future forest populations will be affected by climate change, as these are the seed sources for new tree line and alpine populations. Early life history stages are particularly sensitive to climate and are also influenced by genetic variation among populations. We tested the climate sensitivity of germination and initial development for two widely distributed subalpine conifers, using controlled-environment growth chambers with one temperature regime from subalpine forest in the Colorado Rocky Mountains and one 5 °C warmer, and two soil moisture levels. We tracked germination rate and timing, rate of seedling development, and seedling morphology for two seed provenances separated by ~300 m elevation. Warming advanced germination timing and initial seedling development by a total of ~2 weeks, advances comparable to mean differences between provenances. Advances were similar for both provenances and species; however, warming reduced the overall germination rate, as did low soil moisture, only for Picea engelmannii. A three-year field warming and watering experiment planted with the same species and provenances yielded responses qualitatively consistent with the lab trials. Together these experiments indicate that in a warmer, drier climate, P. engelmannii germination, and thus regeneration, could decline, which could lead to declining subalpine forest populations, while Pinus flexilis forest populations could remain robust as a seed source for upslope range shifts.

54 ENVIRONMENTAL SCIENCES↗

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗