Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion model”

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,297 records · Page 72

Application of a Reynolds stress turbulence model to a supersonic radiating hydrogen-air diffusion flame

A second-order differential Reynolds Stress turbulence model has been applied to the Favre-averaged Navier-Stokes equations for the study of supersonic flows with finite-rate chemistry and radiation. An assumed Beta Probability Density Function is applied to account for the chemical source terms and the radiative flux terms in the conservation equations. A seven-species, seven-reaction finite rate chemistry mechanism is used to simulate the combustion process. The tangent slab approximation is used in radiative flux formulation. A pseudo-gray gas model is used to represent the absorption-emission characteristics of the participating species. The turbulence/radiation interaction is achieved via a new formulation. The resulting formulation is validated by comparison with experimental data on reacting supersonic axisymmetric jets. Results obtained for specific conditions indicate that the effect of chemical reaction on the turbulence is significant. Also, the radiative heat transfer is enhanced by the turbulence.

Chandrasekhar, R.↗

Development and assessment of models for turbulent Rayleigh-Taylor mixing using the macroscopic forcing method

Reynolds-Averaged Navier Stokes (RANS) simulations are a popular method for designing ICF experiments, and accurate mixing models are crucial for these simulations to give good predictions. To this end, the present work seeks to demonstrate the Macroscopic Forcing Method (MFM) as a tool for both improving existing RANS models as well as assessing RANS model forms. First, MFM analysis from Lavacot et al. (Phys. Rev. Fluids, 2025) is used to develop the k–L–F model, an extension of the k–L model of Dimonte and Tipton (Phys. Fluids, 2006) that incorporates nonlocality through addition of a turbulent species flux transport equation. MFM is then applied to the k–L–F model along with the k–L and BHR–4 models to assess their forms and compare the model-implied eddy diffusivity moments to those measured from high-fidelity simulations. Furthermore, the analysis reveals that models incorporating nonlocality (k–L–F and BHR–4) match the high-fidelity simulation data better than purely local models (k–L), both in terms of mean fields and eddy diffusivity moments. However, all of the considered RANS models struggle to match temporal moments at high Atwood numbers, highlighting the importance of temporal nonlocality in these regimes and the need for additional improvement even among models incorporating nonlocality.

general physics↗

Diffusive shock acceleration - Comparison of a unified shock model to bow shock observations

A comparison is made between recent AMPTE/IRM observations of diffuse ions detected upstream of the earth's bow shock when the interplanetary magnetic field was nearly parallel to the solar wind direction and a known collisionless quasi-parallel shock model. These observations, which provide the proton spectrum for all of velocity space, give a direct measure of the shock acceleration efficiency and show how thermal solar winds are injected into the diffusive shock acceleration mechanism. The model accurately describes the proton spectrum; moreover, it predicts the shock structure, the complete particle spectrum, the relative velocity distributions of different ion species in the down stream region, and the enhancement of heavy ions over protons relative to the unshocked solar wind.

Ellison, Donald C.↗

Theoretical estimates of flammability bounds for thin condensed fuel diffusion flames in microgravity using detailed models of chemistry and radiation

Recently, U-shape flammability maps have been constructed showing minimal oxygen vs. flame strain for opposed flame spread in micro-gravity by Olson and Ferkul. The U-shape defines the limiting flammability bounds from radiative extinction and flame blow-off. Here, the minimum of the U corresponds to the minimum possible oxidizer concentration where burning can occur, and is an important quantity of interest for fire safety. While high strain extinction bounds have been well analyzed, low strain radiative extinction has not. To estimate low strain extinction, in this study an analytical theory is developed based on thin flame theory coupled with a heat and mass transfer model for solid fuels. A reaction progress variable based on the Damköhler number is adapted in the theory to account for incomplete combustion at high strain rates and enable the capturing of the full flammability map. The analytical model is compared to a one dimensional numerical model w/ detailed chemical kinetics and coupled radiation heat transfer in planar and spherical geometries. The flammability maps are then qualitatively compared to experimental extinguishment data compiled by Olson and Ferkul for cylindrical rods of PMMA showing similar trends. The results show the newly developed analytics capture the radiative extinction bound compared to the numerical model and qualitatively agrees with microgravity data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

A computational method for subsonic compressible flow in diffusers

Integral methods are used to model subsonic flow in planar diffusers operating with or without stall. Turbulent boundary layers are computed with a lag-entrainment method employing new correlations for shape factor and skin friction, which are based on a Coles-Van Driest velocity profile. The inviscid flow is modeled as one dimensional and is easily coupled with the boundary layer method to yield equations that can be solved by integrating in the downstream direction. Analysis of the equations describing flow in a diffuser shows that the boundary layer separation singularity may not be eliminated by viscous-inviscid interaction; a modified interaction scheme is used to guarantee a non-singular equation set. The nature of this singularity offers an explanation for the apparent universality of asymmetric stall in planar diffusers. Results for several flows demonstrate the accuracy and range of applicability of the method.

Childs, R. E.↗

The structure of hydrophobic gas diffusion electrodes.

The 'flooded agglomerate' model of the Teflon-bonded gas diffusion electrode is discussed. A mathematical treatment of the 'flooded agglomerate' model is given; it can be used to predict the performance of the electrode as a function of measurable physical parameters.

Giner, J.↗

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗

Interactions Enhance Ramp Reversal Memory in Locally Phase Separated Materials

The ramp-reversal memory (RRM) effect in metal–insulator transition metal oxides (TMOs), a non-volatile resistance change induced by repeated temperature cycling, has attracted considerable interest in neuromorphic computing and non-volatile memory devices. Our previous defect motion model successfully explained RRM in vanadium dioxide (VO 2 ), capturing observed critical temperature shifts and memory accumulation throughout the sample. However, this approach lacked interactions between metallic and insulating domains. Here, we extend our model by combining a correlated Random Field Ising Model with defect diffusion-segregation, enabling accurate hysteresis modeling while predicting the relationship between RRM and domain interactions. Our simulations demonstrate that the maximum RRM occurs when the turnaround temperature approaches the inflection point. This peak in RRM vs. turnaround temperature is consistent with prior transport measurements, as well as our own optical measurements reported here. Significantly, we find that increasing nearest-neighbor interactions enhances the maximum memory effect, thus providing a clear mechanism for optimizing RRM performance. Since our model employs minimal assumptions, we predict that RRM should be a widespread phenomenon in materials exhibiting patterned phase coexistence of electronic domains. This work not only advances fundamental understanding of memory behavior in TMOs but also establishes a much-needed theoretical framework for optimizing device applications.

36 MATERIALS SCIENCE↗

Analytical Study of Gravity Effects on Laminar Diffusion Flames

A mathematical model is presented for the description of axisymmetric laminar-jet diffusion flames. The analysis includes the effects of inertia, viscosity, diffusion, gravity and combustion. These mechanisms are coupled in a boundary layer type formulation and solutions are obtained by an explicit finite difference technique. A dimensional analysis shows that the maximum flame width radius, velocity and thermodynamic state characterize the flame structure. Comparisons with experimental data showed excellent agreement for normal gravity flames and fair agreement for steady state low Reynolds number zero gravity flames. Kinetics effects and radiation are shown to be the primary mechanisms responsible for this discrepancy. Additional factors are discussed including elipticity and transient effects.

Edelman, R. B.↗

Microstructural Modeling of the α + β Phase in Ti-6Al-4V: A Diffusion-Based Approach

Complex heat treatment operations and advanced manufacturing processes such as laser or electron-beam welding will see the metallic workpiece experience a considerable range of temperatures and heating/cooling rates. These intrinsic conditions will have a significant bearing upon the microstructure of the material, and in turn upon the thermo-mechanical properties. In this work, a diffusion-based approach to model the growth and shrinkage of precipitates in the alpha + beta field of the common titanium alloy Ti-6Al-4V is established. Further, the numerical model is extended using a JMA-type approach to explore the dependency of the beta-transus temperature on extremely high heating rates, whereby dissolution alone is too slow to accurately describe the alpha to beta-phase transformation. Experimental heat treatments at rates of 5, 50, and 500 °C/s were performed, and metallographic analysis of the samples was used to validate the numerical modeling framework predictions for lamellar shrinkage, while data from the literature has been used to evaluate the numerical modeling framework for the growth of equiaxed microstructures. The agreement between measurements and numerical predictions was found to be good.

36 MATERIALS SCIENCE↗

Numerical simulation of flow through inlets/diffusers with terminal shocks

Two-dimensional Euler and Navier-Stokes solutions of the flow through three inlet/diffuser configurations with terminal shock systems are reported. Calculations without bleed indicate that the terminal shock location is very sensitive to the outflow back pressure. For cases where there are no available experimental results, it becomes difficult to estimate the back pressure that will result in a terminal shock. Estimates based on quasi-one-dimensional analysis are not found adequate for complex two-dimensional flows. It is found that since the flow downstream of the terminal shock is subsonic, and what happens at the outflow boundary affects the flow inside the inlet, enough of the subsonic diffuser must be modeled to accurately predict the terminal shock region. The diffuser portion should be fairly long with the outflow boundary occurring in a region of more or less uniform flow to be able to prescribe a uniform back pressure. The third configuration studied was investigated with and without incorporating bleed in the code. It is found that the use of bleed stabilizes the shock location and allows solutions which without bleed result in unstarting of the inlet. Comparisons are made with available experimental data.

Talcott, N. A., Jr.↗

The effect of breaking gravity waves on the dynamics and chemistry of the mesosphere and lower thermosphere (invited review)

The influence of breaking gravity waves on the dynamics and chemical composition of the 60 to 110 km region is investigated with a two dimensional model that includes a parameterization of gravity wave momentum deposition and diffusion. The dynamical model is described by Garcia and Solomon (1983) and Solomon and Garcia (1983) and includes a complete chemical scheme for the mesosphere and lower thermosphere. The parameterization of Lindzen (1981) is used to calculate the momentum deposited and the turbulent diffusion produced by the gravity waves. It is found that wave momentum deposition drives a very vigorous mean meridional circulation, produces a very cold summer mesopause and reverse the zonal wind jets above about 85 km. The seasonal variation of the turbulent diffusion coefficient is consistent with the behavior of mesospheric turbulences inferred from MST radar echoes. The large degree of consistency between model results and various types of dynamical and chemical data supports very strongly the hypothesis that breaking gravity waves play a major role in determining the zonally-averaged dynamical and chemical structure of the 60 to 110 km region of the atmosphere.

Garcia, R. R.↗

Baffle aperture design study of hollow cathode equipped ion thrusters

A simple theoretical model which can be used as an aid in the design of the baffle aperture region of a hollow cathode equipped ion thruster was developed. An analysis of the ion and electron currents in both the main and cathode discharge chambers is presented. From this analysis a model of current flow through the aperture, which is required as an input to the design model, was developed. This model was verified experimentally. The dominant force driving electrons through the aperture was the force due to the electrical potential gradient. The diffusion process was modeled according to the Bolm diffusion theory. A number of simplifications were made to limit the amount of detailed plasma information required as input to the model to facilitate the use of the model in thruster design. This simplified model gave remarkably consistant results with experimental results obtained with a given thruster geometry over substantial changes in operating conditions. The model was uncertain to about a factor of two for different thruster cathode region geometries. The design usefulness was limited by this factor of two uncertainty and by the accuracy to which the plasma parameters required as inputs to the model were specified.

Brophy, J. R., Jr.↗

A model for the trap-assisted tunneling mechanism in diffused n-p and implanted n(+)-p HgCdTe photodiodes

A theoretical model for the trap-assisted tunneling process in diffused n-on-p and implanted n(+)-on-p HgCdTe photodiodes is presented. The model describes the traps and the trap characteristics: concentration, energy level, and capture cross sections. We have observed that the above two types of diodes differ in the voltage dependence of the trap-assisted tunneling current and dynamic resistance. Our model takes this difference into account and offers an explanation of the phenomenon. The good fit between measured and calculated DC characteristics of the photodiodes (for medium and high reverse bias and for temperatures from 65 to 140 K) supports the validity of the model.

Rosenfeld, David↗

Measuring Fault Zone and Host Rock Hydraulic Properties Using Tidal Responses

Pore pressure diffusion along faults influences induced seismicity and rupture mechanics. In situ hydraulic diffusivity measurements along faults are rare and generally lower than inferred from seismicity migration. Here we use the tidal response of deep geothermal boreholes to measure fault diffusivity and permeability. Initial interpretations of the observation with a homogeneous confined aquifer model result in diffusivities of 10 -3 –10 -1 m 2 /s. However, this model mixes signals from both the conduit and the host rock. We develop a model for tidal response with a fault passing through the aquifer based on the fault-guided fracture network and solve for hydraulic properties in both the fault and the host rock. The resulting fault permeability is 2 × 10 -14 –7 × 10 -14 m 2 (90% confidence interval[CI]) and fault diffusivity is 0.08–0.33 m 2 /s (90% CI), which is 2 orders of magnitude higher than the host rock diffusivity in some wells, thus highlighting the role of faults as fluid conduits.

58 GEOSCIENCES↗

Modeling Studies of the MODIS Solar Diffuser Attenuation Screen and Comparison with On-Orbit Measurements

The MODIS instrument relies on solar calibration to achieve the required radiometric accuracy. This solar calibration occurs as the TERRA spacecraft comes up over the North Pole. The earth underneath the spacecraft is still dark for approximately one minute and the sun is just rising over the earth's north polar regions. During this time the sun moves through about 4 degrees, the scan mirror rotates about 19 times and about 50 frames (exposures) are made of the white solar diffuser. For some of MODIS'S bands the brightness of the screen is reduced, to prevent detector saturation, by means of a pinhole screen, which produces approximately 600 pinhole images of the sun within the field of view of any one detector. Previous attempts at creating a detailed radiometric model of this calibration scenario produced intensity variations on the focal planes with insufficient detail to be useful. The current computational approach produces results, which take into account the motion of the sun and the scan mirror and produces variations, which strongly resemble the observed focal plane intensity variations. The computational approach and results and a comparison with actual observational data are presented.

Waluschka, Eugene↗