Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Dry Reforming of Ethane and Butane with CO 2 over PtNi/CeO 2 Bimetallic Catalysts

Dry reforming is a potential process to convert CO 2 and light alkanes into syngas (H 2 and CO), which can be subsequently transformed to chemicals and fuels. Here in this work, PtNi bimetallic catalysts have been investigated for dry reforming of ethane and butane using both model surfaces and supported powder catalysts. The PtNi bimetallic catalyst shows an improvement in both activity and stability as compared to the corresponding monometallic catalysts. The formation of PtNi alloy and the partial reduction of Ce 4+ to Ce 3+ under reaction conditions are demonstrated by in-situ Ambient Pressure X-ray Photoemission Spectroscopy (AP-XPS), X-ray Diffraction (XRD) and X-ray Absorption Fine Structure (XAFS) measurements. A Pt-rich bimetallic surface is revealed by Diffuse Reflectance Infrared Fourier Transform Spectroscopy (DRIFTS) following CO adsorption. Combined in-situ experimental results and Density Functional Theory (DFT) calculations suggest that the Pt-rich PtNi bimetallic surface structure would weaken the binding of surface oxygenates/carbon species and reduce the activation energy for C-C bond scission, leading to an enhanced dry reforming activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

FunDiff: diffusion models over function spaces for physics-informed generative modeling

Recent advances in generative modeling-particularly diffusion models and flow matching-have been widely used for synthesizing discrete data such as images and videos. However, adapting these models to physical applications remains challenging, as the quantities of interest are continuous functions governed by complex physical laws. To address this, we introduce FunDiff, an efficient and robust framework for generative modeling in function spaces. FunDiff combines a latent diffusion process with a function autoencoder architecture to handle input functions with varying discretizations, generates continuous functions that can be evaluated at arbitrary locations, and seamlessly incorporate physical priors. These priors are enforced through architectural constraints or physics-informed loss functions, ensuring that generated samples satisfy fundamental physical laws. We theoretically establish minimax optimality guarantees for density estimation in function spaces, demonstrating that diffusion-based estimators achieve optimal convergence rates under suitable regularity conditions. We further demonstrate the practical effectiveness of FunDiff across diverse applications in fluid dynamics and solid mechanics. Empirical results indicate that our method can generate physically consistent samples with high fidelity to the target distribution, and exhibit robustness to noisy and low-resolution data.

Wang, Sifan [Yale University, New Haven, CT (Unite↗
DOE OSTI2017

The physics of long- and intermediate-wavelength asymmetries of the hot spot: Compression hydrodynamics and energetics

To achieve ignition with inertial confinement fusion (ICF), it is important to under- stand the effect of asymmetries on the hydrodynamics and energetics of the compres- sion. This paper describes a theoretical model for the compression of distorted hot spots, and quantitative estimates using hydrodynamic simulations. The asymmetries are categorized into low (Ι < 6) and intermediate (Ι < A < 40) modes by comparison of the wavelength with the thermal-diffusion scale length. Long-wavelength modes introduce substantial nonradial motion, whereas intermediate-wavelength modes in- volve more cooling by thermal ablation. We discover that for distorted hot spots, the measured neutron-averaged properties can be very different from the real hydro- dynamic conditions. This is because mass ablation driven my thermal conduction introduces flows in the Rayleigh–Taylor bubbles, this results in pressure variation, in addition to temperature variation between the bubbles and the neutron-producing region (~1 keV for intermediate modes). The differences are less pronounced for long-wavelength asymmetries since the bubbles are relatively hot and sustain fusion reactions. The yield degradation$-$ with respect to the symmetric$-$ results primarily from a reduction in the hot-spot pressure for low modes and from a reduction in burn volume for intermediate modes. It is shown that the degradation in internal energy of the hot-spot is equivalent for both categories, and is equal to the total residual energy in the shell including the bubbles. This quantity is correlated with the shell residual kinetic energy for low-modes, and includes the kinetic energy in the bubbles for mid-modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Mechanical behavior at high temperatures of highly oxygen- or hydrogen-enriched α and prior-β phases of Zirconium alloys - 2016-0063

During a hypothetical loss-of-coolant accident, zirconium alloy fuel claddings can be loaded by internal pressure and exposed to steam at high temperatures (HTs) up to 1,200 deg. C and then cooled and quenched in water. A significant fraction of the oxygen that reacts with the cladding during HT oxidation diffuses beneath the oxide through the metallic substrate. This induces a progressive transformation of the metallic β{sub Zr}-phase layer into an intermediate layer of α{sub Zr}(O) phase containing up to 7 weight % (wt.%) oxygen. Furthermore, in some specific conditions, the cladding may rapidly absorb a significant amount of hydrogen during steam exposition at HTs. As a β{sub Zr} stabilizer, hydrogen would mainly diffuse and concentrate up to several thousands of parts per million by weight (wppm) into the inner β{sub Zr}-phase layer. This study intends to provide new and more comprehensive data on the HT mechanical behavior of the α{sub Zr}(O) and (prior-)β{sub Zr} phases containing high contents of oxygen and hydrogen, respectively. Model samples produced from M5{sup R} and Zircaloy-4 cladding tubes homogeneously charged in oxygen (less than or equal to 6 wt.%) and hydrogen (less than or equal to 3,000 wppm), respectively, were prepared. Their mechanical behavior was determined under vacuum between 800 and 1,100 deg. C for the oxygen-enriched α{sub Zr} phase and in air between 700 and 20 deg. C after cooling from the β{sub Zr} temperature domain for the hydrogen-enriched (prior-)β{sub Zr} phase. The α{sub Zr} phase was substantially strengthened and embrittled by oxygen. The contribution of the α{sub Zr}(O) layer to the HT creep behavior of an oxidized fuel cladding tube subjected to internal pressure is evaluated by finite element analysis. Mechanical strength and ductility of the model (prior-)β{sub Zr} phase appear to be affected by hydrogen contents of 2000-3000 wppm in ways that depend on temperature. (authors)

36 MATERIALS SCIENCE↗
NASA NTRS2014

A Computational and Experimental Study of Coflow Laminar Methane/Air Diffusion Flames: Effects of Fuel Dilution, Inlet Velocity, and Gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4-air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (mu g) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90 cm/s. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton's method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

flames↗
NASA NTRS2015

Vibroacoustic Tailoring of a Rod-Stiffened Composite Fuselage Panel with Multidisciplinary Considerations

An efficient multi-objective design tailoring procedure seeking to improve the vibroacoustic performance of a fuselage panel while maintaining or reducing weight is presented. The structure considered is the pultruded rod stitched efficient unitized structure, a highly integrated composite structure concept designed for a noncylindrical, next-generation flight vehicle fuselage. Modifications to a baseline design are evaluated within a six-parameter design space including spacing, flange width, and web height for both frame and stringer substructure components. The change in sound power radiation attributed to a design change is predicted using finite-element models sized and meshed for analyses in the 500 Hz, 1 kHz, and 2 kHz octave bands. Three design studies are carried out in parallel while considering a diffuse acoustic field excitation and two types of turbulent boundary-layer excitation. Kriging surrogate models are used to reduce the computational costs of resolving the vibroacoustic and weight objective Pareto fronts. The resulting Pareto optimal designs are then evaluated under a static pressurization ultimate load to assess structural strength and stability. Results suggest that choosing alternative configurations within the considered design space can reduce weight and improve vibroacoustic performance without compromising strength and stability of the structure under the static load condition considered, but the tradeoffs are significantly influenced by the spatial characteristics of the assumed excitation field.

Allen, Albert R.↗
NASA NTRSDate not supplied

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗
NASA NTRSDate not supplied

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗
NASA NTRSDate not supplied

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗
NASA NTRSDate not supplied

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗
DOE OSTI2019

Physical and Thermodynamic Properties of Secondary Organic Materials for Modeling (Final Technical Report)

Secondary organic material constitutes more than 50% of the submicron component of atmospheric aerosol particles. This organic material is produced by oxidation of biogenic and anthropogenic volatile organic compounds. The physicochemical properties of these aerosol particles and their variability with atmospheric conditions influence many climate-relevant endpoints, such as rates of chemical reactions, production of secondary organic material, number-diameter distributions, direct radiative forcing, and cloud condensation nuclei. The results of the project research highlight a nonlabile-to-labile transition for material representative of anthropogenic secondary organic particulate matter at a threshold humidity range. This behavior differs markedly from materials representing biogenic aerosol sources, which lack such a transition. More specifically, the mass lability of atmospheric organic particulate matter, meaning its tendency to evaporate, is an important property in the mechanisms governing the climate effects. This property is affected by both the volatility and the molecular diffusivity of the constituent organic species. Via a series of laboratory experiments, the project directly determined the mass labilities of films of secondary organic material representative of anthropogenic and biogenic organic particulate matter (PM). The evaporation rates for anthropogenic organic PM increased above a threshold relative humidity (RH) between 20% and 30%, indicating rapid partitioning above a transition RH. Below the threshold, the characteristic time for equilibration was estimated at up to 1 week for a typically sized particle. In contrast, for films representing biogenic PM, no RH threshold was observed, suggesting equilibrium partitioning is rapidly obtained for all RHs. The effective molecular diffusivity for the biogenic case was at least 1000 times greater than that of the anthropogenic case. The results allow for an improved representation of anthropogenic and biogenic secondary organic aerosol formation in climate models. The results of the project also provide new measurements on how organic aerosol particles take up water and serve as cloud condensation nuclei, thus explaining the observed complex thermodynamic behaviors using a simple thermodynamic model. The project results reconciled the discrepancy between sub-saturation hygroscopicity and super-saturation cloud condensation nuclei activity for organic aerosol particles. The key to the explanation is the phase separation between hydrophobic and hydrophilic organics and the induced surface tension lowering. The results are essential for the assessment of the aerosol indirect effect on climate. In more detail, hygroscopic growth and cloud condensation nuclei activation are key processes for accurately modeling the climate impacts of organic particulate matter. Nevertheless, the microphysical mechanisms of these processes remain unresolved. The project results focus on complex thermodynamic behaviors, including humidity-dependent hygroscopicity, diameter-dependent cloud condensation nuclei activity, and liquid-liquid phase separation in the laboratory for biogenically derived secondary organic material representative of similar atmospheric organic particulate matter. These behaviors can be explained by the non-ideal mixing of water with hydrophobic and hydrophilic organic components. The non-ideality-driven liquid-liquid phase separation further enhances water uptake and induces lowered surface tension at high relative humidity, which results in a lower barrier to cloud condensation nuclei activation. By comparison, secondary organic material representing anthropogenic sources does not exhibit complex thermodynamic behavior. The combined results highlight the importance of detailed thermodynamic representations of the hygroscopicity and cloud condensation nuclei activity in models of the Earth’s climate system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗
DOE OSTI2014

Design of anti-ring back reflectors for thin-film solar cells based on three-dimensional optical and electrical modeling

The optical and electrical properties of a photonic-plasmonic nanostructure on the back contact of thin-film solar cells were investigated numerically through the three-dimensional (3D) finite-difference time-domain method and the 3D Poisson and drift-diffusion solver. The focusing effect and the Fabry-Perot resonances are identified as the main mechanisms for the enhancement of the optical generation rate as well as the short circuit current density. However, the surface topography of certain nanopattern structures is found to reduce the internal electrostatic field of the device, thus limiting charge collection. The optimized conditions for both optics and electronics have been analyzed in this paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Heterogeneous Two-Phase Pillars in Epitaxial NiFe 2 O 4 -LaFeO 3 Nanocomposites

Self-assembled epitaxial oxide nanocomposites have been explored for a wide range of applications, including multiferroic and magnetoelectric properties, plasmonics, and catalysis. These so-called “vertically aligned nanocomposites” form spontaneously during the deposition process when segregation into two phases is energetically favorable as compared to a solid solution. However, there has been surprisingly little work understanding the driving forces that govern the synthesis of these materials, which can include point defect energetics, surface diffusion, and interfacial energies. To explore these factors, La-Ni-Fe-O films have been synthesized by molecular beam epitaxy and it is shown that these phase segregate into spinel-perovskite nanocomposites. Using complementary scanning transmission electron microscopy and atom-probe tomography, the elemental composition of each phase is examined and found that Ni ions are exclusively found in the spinel phase. From correlative analysis, a model for the relative favorability of the Ni2+ and Ni3+ valences under the growth conditions is developed. It is shown that multidimensional characterization techniques provide previously unobserved insight into the growth process and complex driving forces for phase segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Colloidal State Machines as Smart Tracers for Chemical Reactor Analysis

A widely utilized tool in reactor analysis is passive tracers that report the residence time distribution, allowing estimation of the conversion and other properties of the system. Recently, advances in microrobotics have introduced powered and functional entities with sizes comparable to some traditional tracers. This has motivated the concept of Smart Tracers that could record the local chemical concentrations, temperature, or other conditions as they progress through reactors. Herein, the design constraints and advantages of Smart Tracers by simulating their operation in a laminar flow reactor model conducting chemical reactions of various orders are analyzed. It is noted that far fewer particles are necessary to completely map even the most complex concentration gradients compared with their conventional counterparts. Design criteria explored herein include sampling frequency, memory storage capacity, and ensemble number necessary to achieve the required accuracy to inform a reactor model. Cases of severe particle diffusion and sensor noise appear to bind the functional upper limit of such probes and require consideration for future design. The results of the study provide a starting framework for applying the new technology of microrobotics to the broad and impactful set of problems classified as chemical reactor analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Understanding coarsening of a post-corrosion microstructure in a molten salt by combining phase-field modeling and in situ tomography

Alloys corroding in molten salt have been observed to form bicontinuous, nanoporous microstructures via dealloying, which subsequently undergo coarsening due to facile transport in high-temperature conditions. In this work, we describe a methodology to elucidate the underlying transport mechanisms during coarsening of a bicontinuous microstructure via quantitative comparisons between phase-field simulations and four-dimensional in situ experiments, in this case X-ray nanotomography of the coarsening of a dealloyed 80 wt% Ni-20 wt% Cr microwire in molten KCl-MgCl 2 at 800°C. We conduct phase-field simulations initialized from experimental data to model coarsening via three different transport mechanisms: surface diffusion, solid bulk diffusion, and liquid bulk diffusion. These simulations reproduce key features of the experiment, such as the densification of the outer layer of the dealloyed wire and the reduction in radius over time. We quantitatively compare different microstructural characteristics between the simulations and experiment and extract temporal scaling factors that optimally match the time scales of the simulations to that of the experiment. This allows us to evaluate morphological similarity between the simulations and experiment and relate the experimental coarsening kinetics to fundamental material properties. We find that surface diffusion is most likely to be the dominant coarsening mechanism, and its kinetics imply a surface diffusivity of D S = 8.9 x 10 -20 m 3 /s, which is within the range of reported values for Ni-vacuum interfaces at 800°C. However, the difference between the experiment and the surface diffusion simulation increases substantially at late times, suggesting that other mechanisms, such as the dissolution of residual Cr, may be at play.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Reproducing curvature effects due to differential diffusion in tabulated chemistry for premixed flames

Tabulated chemistry reduces the cost of modeling reacting flows by transporting fewer scalars. It also has the benefit of introducing fewer terms that must be closed for Reynolds-averaged Navier-Stokes and large eddy simulation modeling. In this study of premixed lean hydrogen–air flames, an extension to an existing tabulated chemistry model was developed by including differential diffusion and thermal diffusion of multiple species, using a mixture-averaged diffusion model. The transport equations for the tabulation variables (a progress variable and a mixture fraction-like variable) have been formulated, and curvature effects are intrinsic to the transport equations derived from first principles. The mathematical derivation is general, such that other fuel/air mixtures could be examined. The new tabulated chemistry model was then validated using a range of configurations, including flames with curvature, extinction, and hot spots. The model validation indicated excellent agreement with detailed chemistry simulations for lean premixed hydrogen–air mixtures. Conditional means of the progress variable source term were shown to agree well with detailed chemistry simulations of two-dimensional freely propagating flames and three-dimensional turbulent flames. Finally, progress variable source terms three times higher than one-dimensional flat flame results were identified and validated with detailed chemistry, and super-adiabatic hot spots were also predicted in the turbulent flame configuration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Kinetic and Mechanistic Evaluation of Inorganic Arsenic Species Adsorption onto Humic Acid Grafted Magnetite Nanoparticles

In this study, humic acid coated magnetic iron oxide nanoparticles (HA-MNP) were synthesized, characterized, and studied for the removal of toxic inorganic arsenic species from aqueous media. The adsorption of As(III) and As(V) followed pseudo-second order kinetics and the observed data were accurately modeled employing Freundlich adsorption isotherm. Application of the Weber and Morris intraparticle diffusion model to the observed kinetic data suggest that the adsorption occurs in three distinct stages, fast, intermediate, and slow steps. We propose the initial stage is governed by surface association, followed by intraparticle diffusion of arsenic through the HA matrix and finally, chemical reaction or bonding between the arsenic species and HA functionality. The HAMNP nano-adsorbent absorbs > 95% of the inorganic arsenic species with an absorption capacity of 12.2-12.6 mg/g from aqueous media and is effective under a variety of conditions. Post arsenic adsorption characterization of the nanoparticles suggests that As(III) binds with the carboxylate group of HA through a proposed ester type linkage while electrophilic reactions can occur between the nucleophilic functional groups present in HA and the electrophilic arsenic atom in As(V). Finally, the results obtained demonstrated that HA-MNP are robust and have promise for effective As (III) and As(V) remediation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗