Search NASA⌕ Search

SEARCH · Search NASA

Results for “POTENTIAL GRADIENT”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Applications of velocity potential function to acoustic duct propagation and radiation from inlets using finite element theory

A finite element velocity potential program was developed to study acoustic wave propagation in complex geometries. For irrotational flows, relatively low sound frequencies, and plane wave input, the finite element solutions showed significant effects of inlet curvature and flow gradients on the attenuation of a given acoustic liner in a realistic variable area turbofan inlet. The velocity potential approach can not be used to estimate the effects of rotational flow on acoustic propagation, since the potential acoustic disturbances propagate at the speed of the media in sheared flow. Approaches are discussed that are being considered for extending the finite element solution to include the far field, as well as the internal portion of the duct. A new matrix partitioning approach is presented that can be incorporated in previously developed programs to allow the finite element calculation to be marched into the far field. The partitioning approach provided a large reduction in computer storage and running times.

Baumeister, K. J.↗
DOE OSTI2025

Strategies to mitigate urban heat: Effects on overheating and cloud formation

This study evaluates the effectiveness of various urban heat island (UHI) overheating mitigation strategies and their associated impacts on urban cloud dynamics and thermal processes. This study shows cloud-resolving and urban-resolving modeling results estimating the impact of Houston-Galveston heat mitigation scenarios and other resilient strategies contemplated in the Climate Adaptation Plan and Resilient Houston reports. The simulated scenarios include high intensity green rooftops, rooftop photovoltaic solar panels, enhanced urban irrigation, white/cool roofs and roads, and street trees. We contrast the adaptation scenarios against a present baseline case, a no city scenario and a larger and denser city as projected by the Houston-Galveston Area Council for 2045. During the daytime, cooling strategies such as cool roofs, cool roads, and green roofs exhibit superior performance in mitigating urban overheating. At night, enhanced urban irrigation emerges as the most effective cooling intervention. Cooling strategies significantly reduce sensible heat flux partitioning during the day, a process that reduces the uplift of air, suppressing the formation of urban shallow cumulus clouds. The extent of urban cloud mixing ratio is reduced in proportion to the decrease in sensible heating. Under the BEP-Tree scenario, which includes wind effects and evapotranspiration driven by a stomatal conductance model, urban trees demonstrated negligible environmental cooling effects and minimal urban cloud modifications. In contrast, the scenarios with more urban cooling and higher latent heat fluxes led to suppressed urban clouds. The net cooling effect achieved by the heat mitigation strategies is influenced by a combination of indirect processes, including the reduction of downwelling longwave radiation flux, due to reduced cloud presence, while some warming is attributed to a modest increase in shortwave radiation that offsets the cooling benefit. Additionally, reduced heat dissipation, weakened thermal gradients, and diminished vertical mixing over urban areas further moderate the cooling potential. These findings highlight the pivotal role of clouds and moist atmospheric processes in shaping the UHI effect and offer insights for designing more effective urban cooling strategies.

albedo changes↗
DOE OSTI2025

A quantum eigenvalue solver based on tensor networks

Electronic ground states are of central importance in chemical simulations, but have remained beyond the reach of efficient classical algorithms except in cases of weak electron correlation or one-dimensional spatial geometry. We introduce a hybrid quantum-classical eigenvalue solver that constructs a wavefunction ansatz from a linear combination of matrix product states in rotated orbital bases, enabling the characterization of strongly correlated ground states with arbitrary spatial geometry. The energy is converged via a gradient-free generalized sweep algorithm based on quantum subspace diagonalization, with a potentially exponential speedup in the off-diagonal matrix element contractions upon translation into compact quantum circuits of linear depth in the number of qubits. Chemical accuracy is attained in numerical experiments for both a stretched water molecule and an octahedral arrangement of hydrogen atoms, achieving substantially better correlation energies compared to a unitary coupled-cluster benchmark, with orders of magnitude reductions in quantum resource estimates and a surprisingly high tolerance to shot noise. This proof-of-concept study suggests a promising new avenue for scaling up simulations of strongly correlated chemical systems on near-term quantum hardware.

chemistry↗
DOE OSTI2024

Decoupling of redox processes from soil saturation in Arctic tundra

Permafrost thaw in warming Arctic landscapes alters hydrology and saturation-driven biogeochemical processes. Models assume that aerobic respiration occurs in drained soils while saturated soils support methanogenesis; however, saturated soils maintain redox gradients that host a range of anaerobic metabolisms. We evaluated how redox potential and redox-active solutes vary with soil moisture in the active layer of permafrost-affected acidic and non-acidic tundra hillslopes. Oxidizing conditions persisted in highly permeable organic horizons of both unsaturated tussock tundra and saturated wet sedge meadows. Redox potential decreased with depth in all soils as increasing soil bulk density restricted groundwater flow and oxygen diffusion. High concentrations of dissolved iron, phosphate, and organic carbon coincided with redox boundaries below the soil surface in acidic tundra, indicating active iron redox cycling and potential release of adsorbed phosphate during iron (oxyhydr)oxide dissolution. In non-acidic tundra, weatherable minerals affected nutrient dynamics more than redox-driven iron cycling, especially in low-lying, saturated areas where thaw reached mineral soils. The role of thaw depth and the ability of saturated soils to maintain oxidizing conditions in organic surface layers highlight the importance of soil physical properties and hydrology in predicting biogeochemical processes and greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Electric fields and stationary drift flows in the island divertor SOL of Wendelstein 7-X

Abstract In the island divertor scrape-off layer (SOL) of the Wendelstein 7-X stellarator, 2D-resolved Langmuir probe measurements uncover a complex spatial structure of electron temperature, density, and floating potential. Depending on details of magnetic island size and position, significant poloidal gradients and non-monotonic radial profiles are observed. The resulting electric fields imply the presence of sheared stationary drift flows with poloidal velocities of a few km/s, with flow channels as narrow as one cm. Such flows are observed with a gas puff imaging diagnostic in qualitative agreement with the probe results and imply a significant transport channel in the SOL. The 2D distribution of plasma parameters is intricately associated to the unique magnetic structure of the island divertor SOL and highlights the necessity for more accurate magnetic field reconstructions.

Killer, Carsten (ORCID:0000000177473066)↗
DOE OSTI2024

Probing the state of hydrogen in δ − AlOOH at mantle conditions with machine learning potential

Hydrous and nominally anhydrous minerals are a fundamental class of solids of enormous significance to geophysics. They are the water carriers in the deep geological water cycle and impact structural, elastic, plastic, and thermodynamic properties and phase relations in Earth's forming aggregates (rocks). They play a critical role in the geochemical and geophysical processes that shape the planet. Their complexity has prevented predictive calculations of their properties, but progress in materials simulations ushered by machine-learning potentials is transforming this state of affairs. Here, we adopt a hybrid approach that combines deep learning potentials (DPs) with the strongly constrained and appropriately normed meta-generalized gradient approximation functional to simulate a prototypical hydrous system. We illustrate the success of this approach to simulate δ − AlOOH ( δ ) , a phase capable of transporting water down to near the core-mantle boundary of the Earth ( ∼ 2900 km depth and ∼ 135 GPa ) in subducting slabs. A high-throughput sampling of phase space using molecular dynamics simulations with DPs sheds light on the hydrogen-bond behavior and proton diffusion at geophysical conditions. These simulations provide a pathway for a deeper understanding of these crucial components that shape Earth's internal state. Published by the American Physical Society 2024

Luo, Chenxing (ORCID:0000000341166851)↗
DOE OSTI2025

A Near-Real-Time Model for Predicting Electricity Disruptions in Texas During Winter Storms

There has been an increase in extreme weather events, posing a threat to power grid systems, potentially influenced by factors such as population growth, changes in ecosystems, land cover, and land use in the service area, as well as the growth of certain vegetation types. This research seeks to develop a predictive model to mitigate potential damages caused by future winter storms. This research utilizes the Light Gradient Boosting Machine (LightGBM), incorporating the number of power outages experienced at the county level, geographic details, weather information, and lagged outage and lagged weather data. The developed models were broadly divided into two groups, with six models in each group - one group without optimization and another with optimization, totaling 12 trained models. For model optimization, Bayesian optimization was employed using Root Mean Squared Error (RMSE) as the objective function. In results, when comparing Group 2 (the optimized group) with Group 1 (the non-optimized group), it was found that optimization did not always lead to a reduction in RMSE and Mean Absolute Error (MAE). However, in terms of Mean Directional Accuracy (MDA), while all results in Group 1 were below the baseline accuracy of 0.33, all results in Group 2 exceeded 0.33, with some cases showing an increase of more than three times the baseline. The results indicated that, in the optimized model group, Population and Pressure were the most influential factors when using current weather data and geographical information. When using lagged data, lagged recorded outages and lagged Pressure emerged as the most significant factors. Among the 12 developed models, the L-1-2-O model showed the lowest RMSE and MAE, as well as the highest accuracy, with values of 390.62 households and 168.13 households, respectively. To normalize the RMSE and MAE values, each metric was divided by the average number of households among the counties in Texas. For the L-1-2-O model, the scaled RMSE was 0.88% and the scaled MAE was 0.38%. In terms of MDA, which indicates the accuracy of the prediction direction, the L-1-O model achieved the highest score of 0.41. Although this study focused on Texas, which suffered the greatest impact from the winter storms in 2021, with additional validation, the methodology used in this research could be applied to other regions.

Lee, Jangjae [Texas A & M Univ., College Station, ↗
DOE OSTI2022

Continental-scale niche differentiation of dominant topsoil archaea in drylands

Archaea represent a diverse group of microorganisms often associated with extreme environments. However, an integrated understanding of biogeographical patterns of the specialist Haloarchaea and the potential generalist ammonia-oxidizing archaea (AOA) across large-scale environmental gradients remains limited. Here, we hypothesize that niche differentiation determines their distinct distributions along environmental gradients. To test the hypothesis, we use a continental-scale research network including 173 dryland sites across northern China. Our results demonstrate that Haloarchaea and AOA dominate topsoil archaeal communities. As hypothesized, Haloarchaea and AOA show strong niche differentiation associated with two ecosystem types mainly found in China's drylands (i.e. deserts vs. grasslands), and they differ in the degree of habitat specialization. The relative abundance and richness of Haloarchaea are higher in deserts due to specialization to relatively high soil salinity and extreme climates, while those of AOA are greater in grassland soils. Our results further indicate a divergence in ecological processes underlying the segregated distributions of Haloarchaea and AOA. Haloarchaea are governed primarily by environmental-based processes while the more generalist AOA are assembled mostly via spatial-based processes. Our findings add to existing knowledge of large-scale biogeography of topsoil archaea, advancing our predictive understanding on changes in topsoil archaeal communities in a drier world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Leaf habit affects the distribution of drought sensitivity but not water transport efficiency in the tropics

Abstract Considering the global intensification of aridity in tropical biomes due to climate change, we need to understand what shapes the distribution of drought sensitivity in tropical plants. We conducted a pantropical data synthesis representing 1117 species to test whether xylem‐specific hydraulic conductivity (K S ), water potential at leaf turgor loss (Ψ TLP ) and water potential at 50% loss of K S (Ψ P50 ) varied along climate gradients. The Ψ TLP and Ψ P50 increased with climatic moisture only for evergreen species, but K S did not. Species with high Ψ TLP and Ψ P50 values were associated with both dry and wet environments. However, drought‐deciduous species showed high Ψ TLP and Ψ P50 values regardless of water availability, whereas evergreen species only in wet environments. All three traits showed a weak phylogenetic signal and a short half‐life. These results suggest strong environmental controls on trait variance, which in turn is modulated by leaf habit along climatic moisture gradients in the tropics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Leveraging Additive Manufacturing to Fabricate High Temperature Alloys With Co-Designed Mechanical Properties and Environmental Resistance

A paradigm shift in the traditional sequential design approaches is critically essential to create application-specific hierarchical and multifunctional materials with superior long-term performance for next-generation energy technologies involving extreme environments. In the current work, we aim to leverage the flexibility and geometric/compositional complexity offered by additive manufacturing to demonstrate this new approach by codesigning a compositionally graded Ni-based alloy for molten salts\sCO 2 heat exchangers to enable mitigation of environmental degradation of surfaces exposed to molten halide salts, while simultaneously suppressing the consequent deterioration in mechanical stability. Further, thermokinetic modeling describing the underlying physics of thermally- and environmentally induced spatiotemporal compositional and microstructural evolution will be employed to predict the parameter space of material deposition processes and precisely identify the required composition gradient. Preliminary corrosion and mechanical testing of the dual material demonstrated the potential of the material to replace existing solid solution strengthened materials for this application.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗
DOE OSTI2023

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗
NASA NTRS1980

The stability properties of cylindrical force-free fields - Effect of an external potential field

A large-scale potential field with an embedded smaller-scale force-free structure gradient x B equals alpha B is studied in cylindrical geometry. Cases in which alpha goes continuously from a constant value alpha 0 on the axis to zero at large r are considered. Such a choice of alpha (r) produces fields which are realistic (few field reversals) but not completely stable. The MHD-unstable wavenumber regime is found. Since the considered equilibrium field exhibits a certain amount of magnetic shear, resistive instabilities can arise. The growth rates of the tearing mode in the limited MHD-stable region of k space are calculated, showing time-scales much shorter than the resistive decay time.

Chiuderi, C.↗
NASA NTRS1986

Developmental testing of a programmable multizone furnace

A multizone furnace was evaluated for its potential utilization for process experimentation on board the Space Shuttle. A temperature gradient can be created through the use of a series of connected temperature zones and can be translated by the coordinated sequencing of zone temperatures. The Bridgman-Stockbarger thermal configuration for directional solidification was implemented so that neither the sample nor furnace was translated. The thermal behavior of the furnace was measured and characterized. Limitations due to both thermal and electronic (computer) factors are identified. The results indicate that the multizone design is limited to low temperature gradients because of the indirect furnace-to-sample thermal coupling needed to blend the discrete thermal zones. The multizone furnace design inherently consumes more power than a similar (two temperature) conventional Bridgman type directional solidification furnace because every zone must be capable of the high cooling rates needed to produce the maximum desired temperature drop. Typical achievable static temperature gradients for the furnace tested were between 6 and 75 C/in. The maximum gradient velocity was approximately 10 in./hr. Several aspects of the tested system could be improved, but the dependence of the multizone design on high heat loss will limit Space Shuttle applications in the form tested unless additional power is available. The multizone furnace offers great flexibility but requires a high level of operator understanding for full advantage to be obtained.

Ting, E. Y.↗
NASA NTRS1990

Thermocapillary convection in a model float-zone

A finite-element method has been used to study thermocapillary convection in a model of relevance to the float-zone, crystal-growth process. The geometry consists of a pair of horizontal isothermal surfaces with a Boussinesq liquid suspended between them. Because of the variation of surface tension with temperature, the temperature gradient along the free surface drives thermocapillary convection. The free surface is allowed to deform and its location is calculated along with the velocity and temperature fields. Cases in which thermocapillarity dominates and others in which it interacts with an unstable axial buoyancy gradient are treated. In order to obtain information of possible interest to potential microgravity applications, the influence of the Grashof number is investigated.

Neitzel, G. P.↗
NASA NTRS1992

Kinetic plasma processes occurring in the outer plasmasphere

One area of data analysis work that was begun under this contract is the fitting of the perpendicular velocity distributions of equatorially trapped ions with a Kappa function. This type of characterization of the trapped ions will be very useful for comparison with velocity distributions produced by the model. A second area of data analysis is to study data from consecutive passes when DE 1's apogee was near the magnetic equator and the spacecraft was often skimming along nearly the same L shell. In 1982 three such periods occurred in May, June, and July. For these consecutive events we have Kp histories, density measurements from a number of sources (Whistler data, DE SFR, ISEE SFR) and consecutive samples of ion pitch angle distributions along field lines. It is clear from this data how the pitch angle distributions evolve during a flux tube refilling event. Our modeling of the flow of plasma along closed field lines is following two basic tracks. The first is a study of the basic refilling process without the effect of wave-particle heating near the equator or the effect of large or abrupt field-aligned electric potential drops. This model includes the effects of Coulomb self-collisions and collisions with the O+ ions in the topside ionosphere. The second track is a study of the effects of wave produced pitch-angle scattering and perpendicular heating occurring near the magnetic equator, in connection with the development of large potential drops that result from electron heating and the development of density gradients.

Wilson, Gordon R.↗
NASA NTRS1994

Control theory based airfoil design for potential flow and a finite volume discretization

This paper describes the implementation of optimization techniques based on control theory for airfoil design. In previous studies it was shown that control theory could be used to devise an effective optimization procedure for two-dimensional profiles in which the shape is determined by a conformal transformation from a unit circle, and the control is the mapping function. The goal of our present work is to develop a method which does not depend on conformal mapping, so that it can be extended to treat three-dimensional problems. Therefore, we have developed a method which can address arbitrary geometric shapes through the use of a finite volume method to discretize the potential flow equation. Here the control law serves to provide computationally inexpensive gradient information to a standard numerical optimization method. Results are presented, where both target speed distributions and minimum drag are used as objective functions.

Reuther, J.↗
NASA NTRS1999

Detailed Studies on the Structure and Dynamics of Reacting Dusty Flows at Normal and Microgravity

Dusty reacting flows are of particular interest for a wide range of applications. Inert particles can alter the flammability and extinction limits of a combustible mixture. Reacting particles can release substantial amount of heat and can be used either for power generation or propulsion. Accumulation of combustible particles in air can result in explosions which, for example, can occur in grain elevators, during lumber milling and in mine galleries. Furthermore, inert particles are used as flow velocity markers in reacting flows, and their velocity is measured by non-intrusive laser diagnostic techniques. Despite their importance, dusty reacting flows have been less studied and understood compared to gas phase as well as sprays. The addition of solid particles in a flowing gas stream can lead to strong couplings between the two phases, which can be of dynamic, thermal, and chemical nature. The dynamic coupling between the two phases is caused by the inertia that causes the phases to move with different velocities. Furthermore, gravitational, thermophoretic, photophoretic, electrophoretic, diffusiophoretic, centrifugal, and magnetic forces can be exerted on the particles. In general, magnetic, electrophoretic, centrifugal, photophoretic, and diffusiophoretic can be neglected. On the other hand, thermophoretic forces, caused by steep temperature gradients, can be important. The gravitational forces are almost always present and can affect the dynamic response of large particles. Understanding and quantifying the chemical coupling between two phases is a challenging task. However, all reacting particles begin this process as inert particles, and they must be heated before they participate in the combustion process. Thus, one must first understand the interactions of inert particles in a combustion environment. The in-detail understanding of the dynamics and structure of dusty flows can be only advanced by considering simple flow geometries such as the opposed-jet, stagnation-type. In such configurations the imposed strain rate is well characterized, and the in-depth understanding of the details of the physico-chemical processes can be systematically obtained. A number of computational and experimental studies on spray and particle flows have been conducted in stagnation-type configurations. Numerically, the need for a hybrid Eulerian-Lagrangian approach has been identified by Continillo and Sirignano, and the use of such approach has allowed for the prediction of the phenomenon of droplet flow reversal. Gomez and Rosner have conducted a detailed study on the particle response in the opposed-jet configuration, and the particle thermophoretic diffusivities were determined experimentally. Sung, Law and co-workers have conducted numerical studies on the effect of strain rate and temperature gradients on the dynamics of inert particles, as a way of understanding potential errors in experimental LDV data that may arise from thermophoretic forces. This investigation is a combined experimental and numerical study on the details of reacting dusty flows. The specific tasks are: (1) Experimental determination of laminar flame speeds, and extinction strain rates of dusty flows at normal- and micro-gravity as functions of the particle type, particle initial diameter, particle initial number density, and gas phase chemical composition; (2) Detailed numerical simulation of the experiments. Results are compared with experiments and the adequacy of theoretical models is assessed; and (3) Provision of enhanced insight into the thermo-chemical coupling between the two phases.

Andac, M. Gurhan↗