Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Experimental and theoretical analysis of the rate of solvent equilibration in the hanging drop method of protein crystal growth

The principles of the hanging-drop method of crystal growth are discussed, and the rate of water evaporation in a water droplet (containing protein, buffer, and a precipitating agent) suspended above a well containing a double concentration of precipitating agent is investigated theoretically. It is shown that, on earth, the rate of evaporation may be determined from diffusion theory and the colligative properties of solutions. The parameters affecting the rate of evaporation include the temperature, the vapor pressure of water, the ionization constant of the salt, the volume of the drop, the contact angle between the droplet and the coverslip, the number of moles of salt in the droplet, the number of moles of water and salt in the well, the molar volumes of water and salt, the distance from the droplet to the well, and the coefficient of diffusion of water vapor through air. To test the theoretical equations, hanging-drop experiments were conducted using various reagent concentrations in 25-microliter droplets and measuring the evaporation times at 4 C and 25 C. The results showed good agreement with the theory.

Fowlis, William W.↗
NASA NTRS2018

Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Comparisons with experimental data showed a sensitivity of the current prediction to the modeling of the sidewalls. This was found to be common among various computational studies in the literature where periodic boundary conditions were used in the spanwise direction, as was the case in the present work. Thus, although the experiment was quasi-two-dimensional, the present simulation was determined to be two-dimensional. Quantities present in the exact equation of the Reynolds-stress transport, i.e., production, molecular diffusion, turbulent transport, pressure diffusion, pressure strain, dissipation, and turbulent mass flux were calculated. Reynolds-stress budgets were compared with past large-eddy simulation and direct numerical simulation datasets in the undisturbed portion of the turbulent boundary layer to validate the current approach. The budgets in SBLI showed the growth in the production term for the primary normal stress and energy transfer mechanism was led by the pressure strain term in the secondary normal stresses. The pressure diffusion term, commonly assumed as negligible by turbulence model developers, was shown to be small but non-zero in the normal stress budgets, however it played a key role in the primary shear stress budget.

shock wave/boundary-layer interactio↗
NASA NTRSDate not supplied

Edge Based Viscous Method for Node-Centered Formulations

This paper presents a novel, efficient, conservative, edge-based method for evaluation of mean flow viscous fluxes and turbulence-model diffusion terms of the Reynolds-averaged Navier-Stokes equations on tetrahedral grids. The new method is implemented in a practical, node-centered, finite-volume computational fluid dynamics solver. The baseline finite-volume scheme that is equivalent to a second-order accurate finite-element Galerkin approximation of viscous stresses is reformulated. The order of operations to compute the cell-based Green-Gauss gradients is changed to combine the operations by edge, which leads to an equivalent formulation on tetrahedral grids, improves efficiency, and preserves the compact discretization stencil based on the nearest neighbors. The computational results presented in this paper verify the implementation of this edge-based method by comparing its accuracy and iterative convergence with those of the well verified and validated baseline formulation. Efficiency gains for residual and Jacobian evaluations result in significant reduction of time to solution. This novel edge-based formulation on tetrahedra can be seamlessly combined with the baseline formulation on cells of other types for computing solutions on mixed-element grids.

Edge Based↗
NASA NTRS1993

Numerical vorticity capturing for vortex-solid body interaction problems

A new numerical method based on adding a term to the Euler/Navier-Stokes equations has been demonstrated it can effectively treat vortex-dominated flows using low-order numerical schemes and coarse grids. The numerical diffusion introduced by convection schemes can be eliminated by modifying the velocity field as a result of solving the modified flow governing equations. The modification of the velocity field, which conserves the total vorticity, essentially convects the vorticity toward its local extreme to offset the numerical diffusion. This method preserves the vortex structure even when vortices travel on coarse grid region. The method is implemented with two distinctive codes: a Navier-Stokes flow solver based on vorticity-velocity formulation and an Euler code based on velocity-pressure formulation. Problems presented in this paper include airfoil dynamic stall, vortex airfoil interaction, and vortex-fuselage interactions.

Wang, Clin M.↗
NASA NTRS1983

An EBIC equation for solar cells

When an electron beam of a scanning electron microscope (SEM) impinges on an N-P junction, the generation of electron-hole pairs by impact ionization causes a characteristic short circuit current I(sc) to flow. The I(sc), i.e., EBIC (electron beam induced current) depends strongly on the configuration used to investigate the cell's response. In this paper the case where the plane of the junction is perpendicular to the surface is considered. An EBIC equation amenable to numerical computations is derived as a function of cell thickness, source depth, surface recombination velocity, diffusion length, and distance of the junction to the beam-cell interaction point for a cell with an ohmic contact at its back surface. It is shown that the EBIC equation presented here is more general and easier to use than those previously reported. The effects of source depth, ohmic contact, and diffusion length on the normalized EBIC characteristic are discussed.

Luke, K. L.↗
NASA NTRS1983

Day and night models of the Venus thermosphere

A model atmosphere of Venus for altitudes between 100 and 178 km is presented for the dayside and nightside. Densities of CO2, CO, O, N2, He, and O2 on the dayside, for 0800 and 1600 hours local time, are obtained by simultaneous solution of continuity equations. These equations couple ionospheric and neutral chemistry and the transport processes of molecular and eddy diffusion. Photodissociation and photoionization J coefficients are presented to facilitate the incorporation of chemistry into circulation models of the Venus atmosphere. Midnight densities of CO2 CO, O, N2, He, and N are derived from integration of the continuity equations, subject to specified fluxes. The nightside densities and fluxes are consistent with the observed airglow of NO and O2(1 Delta). The homopause of Venus is located near 133 km on both the dayside and nightside.

Massie, S. T.↗
NASA NTRS2022

Prediction of Turbulent Diffusing Flows Using FUN3D

The FUN3D code was used to perform Reynolds-averaged Navier-Stokes (RANS) simulations to compute subsonic flow in an S-duct diffuser and transonic flow in a two-dimensional diffuser using the Speziale-Sarkar-Gatski/Launder-Rodi-Reece (SSG/LRR) Reynolds stress model (RSM). For comparison purposes, additional simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST) turbulence models. Each model was run with and without the quadratic constitutive relation (QCR) for computing the turbulent stresses. It was shown that the simulations that utilized the RSM had better overall predictions of the diffusive flow fields compared to the simulations that utilized the one and two equation turbulence models.

Inlet Distortion↗
DOE OSTI2024

High-accuracy method for modeling nucleation and growth of particles

State-of-the-art numerical models describing the kinetics of aerosol particle nucleation and growth from a cooling vapor primarily use a nodal method, in which particles that are smaller than the critical size are omitted from consideration because they are thermodynamically unfavorable. This omission is based on the assumption that most newly formed particles are above the critical size, so that subcritical-size particles are not important to take into account. Due to the nature of the nodal method, it suffers from numerical diffusion, which can cause an artificial broadening of the cluster size distribution leading to a significant overestimation of the number of large-size particles. To address these issues, we propose a more accurate numerical method that explicitly models particles of all sizes, and uses a special numerical scheme that substantially reduces the numerical diffusion and provides high solution accuracy and numerical stability. We extensively compare this novel method to the commonly used nodal solver of the general dynamic equation (GDE) for particle growth and demonstrate that it offers GDE solutions with higher accuracy with low numerical diffusion. Incorporating small subcritical clusters into the solution is crucial for: 1) more precise determination of the entire particle size distribution function and 2) wider applicability of the model to experimental studies with non-monotonic temperature variations leading to particle evaporation. The computational code implementing this numerical method in Python is available upon request.

42 ENGINEERING↗
DOE OSTI2026

High-Fidelity, Low-Dissipation/Symmetry-Preserving Numerical Scheme for Solving the Euler Equations with Unstructured, Metric-Based Mesh Adaptation

This work presents an overview of a high-fidelity compressible Euler solver that utilizes the continuous Galerkin (CG) method with added artificial numerical diffusion for stabilization to solve a variety of unsteady and steady benchmark inviscid flow problems. This work shows that discretizing the Euler equations with this CG approach and first order basis functions produces a cost-effective stencil as well as simple well-posed boundary conditions. We show through convergence testing with manufactured solutions that the reduced stencil of CG, combined with the low amount of artificial diffusion required when using the stabilization method outlined in this work, leads to stable and highly accurate results for a variety of unsteady and steady applications. When combined with the adaptive mesh refinement approach used for many of the cases in this work, our results show that the flow solver achieves even more accurate results. A variety of inviscid flow cases are presented in this work, including transient 2D cases with complex shock structures and several steady 3D airfoils sections with a constant span.

Doetsch, Kevin [ORNL] (ORCID:0000000267051705)↗
NASA NTRS2002

Dynamics of Diffusion Flames in von Karman Swirling Flows Studied

Von Karman swirling flow is generated by the viscous pumping action of a solid disk spinning in a quiescent fluid media. When this spinning disk is ignited in an oxidizing environment, a flat diffusion flame is established adjacent to the disk, embedded in the boundary layer (see the preceding illustration). For this geometry, the conservation equations reduce to a system of ordinary differential equations, enabling researchers to carry out detailed theoretical models to study the effects of varying strain on the dynamics of diffusion flames. Experimentally, the spinning disk burner provides an ideal configuration to precisely control the strain rates over a wide range. Our original motivation at the NASA Glenn Research Center to study these flames arose from a need to understand the flammability characteristics of solid fuels in microgravity where slow, subbuoyant flows can exist, producing very small strain rates. In a recent work (ref. 1), we showed that the flammability boundaries are wider and the minimum oxygen index (below which flames cannot be sustained) is lower for the von Karman flow configuration in comparison to a stagnation-point flow. Adding a small forced convection to the swirling flow pushes the flame into regions of higher strain and, thereby, decreases the range of flammable strain rates. Experiments using downward facing, polymethylmethacrylate (PMMA) disks spinning in air revealed that, close to the extinction boundaries, the flat diffusion flame breaks up into rotating spiral flames (refs. 2 and 3). Remarkably, the dynamics of these spiral flame edges exhibit a number of similarities to spirals observed in biological systems, such as the electric pulses in cardiac muscles and the aggregation of slime-mold amoeba. The tail of the spiral rotates rigidly while the tip executes a compound, meandering motion sometimes observed in Belousov-Zhabotinskii reactions.

Nayagam, Vedha↗
DOE OSTI2021

Efficient quadrature rules for finite element discretizations of nonlocal equations

In this paper we design efficient quadrature rules for finite element discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first showthat the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient finite element assembly.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Reflectance and unsteady temperature distribution of diffuse reflectors

The applicability of a previously derived solution of the Kubelka-Munk differential equations for the reflectance of an absorbing and scattering material for the case of a zero absorption coefficient to the description of the reflectance of a weakly absorbing, highly scattering material is examined. A solution for the transient heating of a one-dimensional volume reflector is obtained, and the error incurred in the use of the zero absorption coefficient solution is assessed for early times appropriate to the entry heat protection system application. It is found that the small difference in reflectance causes a significant difference in internal temperature.

Weston, K. C.↗
NASA NTRS1990

Improved Analysis Of Heat-Pulse Data

Equations derived to improve analysis of experimental data in flash method for measurement of thermal diffusivity and heat capacity. In flash method, pulse of radiant energy from flashlamp or other source deposited on front face of small, thin specimen, and temperature of rear face monitored as function of time. New equations account for both losses and exponentially decaying pulses.

Vining, C. B.↗
DOE OSTI2025

Resolving turbulent magnetohydrodynamics: a hybrid operator-diffusion framework

We present a hybrid machine learning framework that combines physics-informed neural operators (PINOs) with score-based generative diffusion models to simulate the full spatio-temporal evolution of two-dimensional, incompressible, resistive magnetohydrodynamic turbulence across a broad range of Reynolds numbers (Re). The framework leverages the equation-constrained generalization capabilities of PINOs to predict coherent, low-frequency dynamics, while a conditional diffusion model stochastically corrects high-frequency residuals, enabling accurate modeling of fully developed turbulence. Trained on a comprehensive ensemble of high-fidelity simulations with Re ϵ {100, 250, 500, 750, 1000, 3000, 10000}, the approach achieves state-of-the-art accuracy in regimes previously inaccessible to deterministic surrogates. At Re = 1000 and 3000, the model faithfully reconstructs the full spectral energy distributions of both velocity and magnetic fields late into the simulation, capturing non-Gaussian statistics, intermittent structures, and cross-field correlations with high fidelity. At extreme turbulence levels (Re = 10 000), it remains the first surrogate capable of recovering the high-wavenumber evolution of the magnetic field, preserving large-scale morphology and enabling statistically meaningful predictions.

Diffusion-Integrated Neural Operators↗
NASA NTRSDate not supplied

Prediction of Turbulent Diffusing Flows Using FUN3D

The FUN3D code was used to perform Reynolds-averaged Navier-Stokes (RANS) simulations to compute subsonic flow in an S-duct diffuser and transonic flow in a two-dimensional diffuser using the Speziale-Sarkar-Gatski/Launder-Rodi-Reece (SSG/LRR) Reynolds stress model (RSM). For comparison purposes, additional simulations were run with the one-equation Spalart-Allmaras (SA) and the two-equation Menter Shear-Stress Transport (SST) turbulence models. Each model was run with and without the quadratic constitutive relation (QCR) for computing the turbulent stresses. It was shown that the simulations that utilized the RSM had better overall predictions of the diffusive flow fields compared to the simulations that utilized the one and two equation turbulence models.

Inlet Distortion↗
NASA NTRS1975

Steady-state solutions of a diffusive energy-balance climate model and their stability

A diffusive energy-balance climate model, governed by a nonlinear parabolic partial differential equation, was studied. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. Models similar to the main one are considered, and the number of their steady states was determined. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The stability under small perturbations of the main model's climates was investigated. A stability criterion is derived, and its application shows that the present climate and the deep freeze are stable, whereas the model's glacial is unstable. The dependence was examined of the number of steady states and of their stability on the average solar radiation.

Ghil, M.↗