Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Analysis of X-Ray Microradiographs of Al-Au Interface Quench Profile using Modeling of Solidification Including Double-Diffusion and Convection in the Melt

Experimental data on Al-0.8Au horizontal solidification of a 1 mm thick specimen in a BN crucible shows the effect of growth rate on the solidification interface shape. For translation rates below 0.5 micron/s the interface maintains a plain and flat shape. When the translation rate is 3 to 5 micron/s or more, the interface appearance changes to two planar zones, with the zone closer to the bottom having higher inclination. The interface shapes were measured by first quenching in place during growth. X-ray microscopy shows the interface shape within the quenched sample by viewing through the side of the specimen. In order to provide theoretical explanation of the phenomena, numerical modeling was undertaken using finite element code FIDAP. Double diffusion convection in Al-0.8Au melt and crystal-melt interface curvature during directional solidification was analyzed numerically. Actual thermophysical properties of Al-0.8Au including the binary Al-Au phase diagram were used. Although convection in the sample is weak, for the slower translation rate convection and diffusion is sufficient for the redistribution of initial compositional stratification caused by gravity. When translation rate is raised, neither convection nor diffusion can provide proper mixing so that solidification temperatures differ significantly near the bottom within the bulk of the sample. As a result, the solid-liquid interface appears to have two planar zones with different inclination.

Bune, Andris V.↗
NASA NTRS1986

Thermophoretically augmented mass transfer rates to solid walls across laminar boundary layers

Predictions of mass transfer (heavy vapor and small particle deposition) rates to solid walls, including the effects of thermal (Soret) diffusion ('thermophoresis' for small particles), are made by numerically solving the two-dimensional self-similar forced convection laminar boundary-layer equations with variable properties, covering the particle size range from vapor molecules up to the size threshold for inertial (dynamical nonequilibrium) effects. The effect of thermophoresis is predicted to be particularly important for submicron particle deposition on highly cooled solid surfaces, with corresponding enhancement factors at atmospheric conditions being over a thousand-fold at T(w)/T(e) equal to about 0.6. As a consequence of this mass transfer mechanism, the particle size dependence of the mass transfer coefficient to a cooled wall will be much weaker than for the corresponding case of isothermal capture by Brownian-convective diffusion.

Gokoglu, S. A.↗
DOE OSTI2017

Hydrophilic domain structure in polymer exchange membranes: Simulations of NMR spin diffusion experiments to address ability for model discrimination

Here, we detail the development of a flexible simulation program (NMR_DIFFSIM) that solves the nuclear magnetic resonance (NMR) spin diffusion equation for arbitrary polymer architectures. The program was used to explore the proton ( 1 H) NMR spin diffusion behavior predicted for a range of geometrical models describing polymer exchange membranes. These results were also directly compared with the NMR spin diffusion behavior predicted for more complex domain structures obtained from molecular dynamics (MD) simulations. The numerical implementation and capabilities of NMR_DIFFSIM were demonstrated by evaluating the experimental NMR spin diffusion behavior for the hydrophilic domain structure in sulfonated Diels-Alder Poly(Phenylene) (SDAPP) polymer membranes. The impact of morphology variations as a function of sulfonation and hydration level on the resulting NMR spin diffusion behavior were determined. These simulations allowed us to critically address the ability of NMR spin diffusion to discriminate between different structural models, and to highlight the extremely high fidelity experimental data required to accomplish this. A direct comparison of experimental double-quantum-filtered 1 H NMR spin diffusion in SDAPP membranes to the spin diffusion behavior predicted for MD-proposed morphologies revealed excellent agreement, providing experimental support for the MD structures at low to moderate hydration levels.

36 MATERIALS SCIENCE↗
NASA NTRS1989

A high-order Lagrangian-decoupling method for the incompressible Navier-Stokes equations

A high-order Lagrangian-decoupling method is presented for the unsteady convection-diffusion and incompressible Navier-Stokes equations. The method is based upon: (1) Lagrangian variational forms that reduce the convection-diffusion equation to a symmetric initial value problem; (2) implicit high-order backward-differentiation finite-difference schemes for integration along characteristics; (3) finite element or spectral element spatial discretizations; and (4) mesh-invariance procedures and high-order explicit time-stepping schemes for deducing function values at convected space-time points. The method improves upon previous finite element characteristic methods through the systematic and efficient extension to high order accuracy, and the introduction of a simple structure-preserving characteristic-foot calculation procedure which is readily implemented on modern architectures. The new method is significantly more efficient than explicit-convection schemes for the Navier-Stokes equations due to the decoupling of the convection and Stokes operators and the attendant increase in temporal stability. Numerous numerical examples are given for the convection-diffusion and Navier-Stokes equations for the particular case of a spectral element spatial discretization.

Ho, Lee-Wing↗
DOE OSTI2018

Internal Resistive Barriers Related to Zinc Diffusion During the Growth of Inverted Metamorphic Multijunction Solar Cells

Majority carrier barriers at heterointerfaces are a common source of non-linear resistance that hinders concentrator solar cell performance. The source of a particular barrier is often unclear in a multijunction device with numerous heterointerfaces. In this work, we demonstrate Zn-dopant diffusion during inverted metamorphic multijunction (IMM) device growth to be one key cause of internal barrier formation. Using an inverted GaAs/GaAs tandem solar cell with a high temperature annealing layer grown in between each subcell, we simulate the annealing conditions of a multijunction growth in a simplified structure. Through analysis of the device by secondary ion mass spectrometry (SIMS) and electrochemical capacitance-voltage profiling, we show that annealing causes Zn to diffuse out of the top cell Ga 0.5 In 0.5 P back surface field (BSF) and accumulate in the GaAs base. Through equilibrium band modeling, we show that the resultant doping profile forms an energetic barrier to hole flow in the valence band, which correlates with fill factor losses in the current-voltage curves measured under concentration. When we, instead, employ a C-doped Al 0.2 Ga 0.8 As BSF layer in the top cell, we do not observe evidence of a heterojunction barrier. We attribute this difference to the reduced diffusivity of carbon, confirmed by SIMS, as well as more favorable valence band offsets between GaAs and Al 0.2 Ga 0.8 As. Finally, we compare 5-junction IMM cells with Al 0.2 Ga 0.8 As:C and Ga 0.5 In 0.5 P:Zn BSF layers in the GaAs third junction, respectively, and show a significantly improved concentrator device performance when Al 0.2 Ga 0.8 As:C is employed. Lastly, we demonstrate the importance of designing annealing tolerance into multijunction structures that are subjected to extended annealing during growth.

14 SOLAR ENERGY↗
DOE OSTI2018

Coiled Tube Gas Heaters For Nuclear Gas-Brayton Power Conversion

This project developed an alternative design for heat exchangers for application to heating supercritical carbon dioxide (S-CO 2 ) or air for power conversion. We have identified an annular coiled tube bundle configuration–where hot sodium enters tubes from multiple vertical inlet manifold pipes, flows in a spiral pattern radially inward and downward, and then exits into an equal number of vertical outlet manifold pipes–as a potentially attractive option. The S-CO 2 gas or air flows radially outward through the tube bundle. Coiled tube gas heaters (CTGHs) are expected to have excellent thermal shock, long-term thermal creep, in-service inspection, and reparability characteristics, compared to alternative options. CTGHs have significant commonality with modern nuclear steam generators. Extensive experience exists with the design, manufacture, operation, in-service inspection and maintenance of nuclear steam generators. The U.S. Nuclear Regulatory Commission also has extensive experience with regulatory guidance documented in NUREG 0800. CTGHs leverage this experience and manufacturing capability. The most important difference between steam generators and gas-Brayton cycles such as the S-CO 2 cycle is that the heat exchangers must operate with counter flow with high effectiveness to minimize the pinch-point temperature difference between the hot liquid coolant and the heated gas. S-CO 2 -cycle gas heaters also operate at sufficiently elevated temperatures that time dependent creep is important and allowable stresses are relatively low. Designing heat exchangers to operate in this regime requires configurations that minimize stresses and stress concentrations. The cylindrical tubes and cylindrical manifold pipes used in CTGHs are particularly effective geometries. The first major goal of this research project was to develop and experimentally validate a detailed, 3-D multi-phase (gas-solid-liquid) heat transport model for CTGHs, using methods similar to earlier UCB multi-scale models for PCHEs, which will enable optimization of CTGH designs with respect to the number of manifold pipes, S-CO 2 and sodium circulating power, and other design parameters. CTGHs are relatively novel because the tubes loaded in compression as well as the manifold pipes which form the tube sheets. The second major goal of the experimental and modeling effort was to optimize the assembly and heat treatment of the tube-to-tubesheet joints, using a novel tapered joint configuration, and to develop diffusion-bonded joints with high shear strength and creep resistance, as has been done in other applications for tapered plugs. The project would also conduct creep testing of these tubesheet joints in pressurized furnaces under prototypical conditions of stress and temperature, however this was not feasible due to numerous hinderances. Detailed microstructural investigations was conducted on the diffusion bond interface, while the bond strength was evaluated by tensile testing. The mechanical testing quantified the overall strength of the bond and measured the pull out force as well as long term behavior, which complemented the microstructural analysis and gave a complete understanding of the diffusion bond.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2019

Nanoparticle diffusion in sheared cellular blood flow

Using a multiscale blood flow solver, the complete diffusion tensor of nanoparticles (NPs) in sheared cellular blood flow is calculated over a wide range of shear rate and haematocrit. In the short-time regime, NPs exhibit anomalous dispersive behaviors under high shear and high haematocrit due to the transient elongation and alignment of the red blood cells (RBCs). In the long-time regime, the NP diffusion tensor features high anisotropy. Particularly, there exists a critical shear rate ( ${\sim}100~\text{s}^{-1}$ ) around which the shear-rate dependence of the diffusivity tensor changes from linear to nonlinear scale. Above the critical shear rate, the cross-stream diffusivity terms vary sublinearly with shear rate, while the longitudinal term varies superlinearly. The dependence on haematocrit is linear in general except at high shear rates, where a sublinear scale is found for the vorticity term and a quadratic scale for the longitudinal term. Through analysis of the suspension microstructure and numerical experiments, the nonlinear haemorheological dependence of the NP diffusion tensor is attributed to the streamwise elongation and cross-stream contraction of RBCs under high shear, quantified by a capillary number. The RBC size is shown to be the characteristic length scale affecting the RBC-enhanced shear-induced diffusion (RESID), while the NP submicrometre size exhibits negligible influence on the RESID. Based on the observed scaling behaviours, empirical correlations are proposed to bridge the NP diffusion tensor to specific shear rate and haematocrit. As a result, the characterized NP diffusion tensor provides a constitutive relation that can lead to more effective continuum models to tackle large-scale NP biotransport applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Spherical flame characteristics of heptane/iso-butanol mixture droplets

Experimental data and detailed numerical modelling are presented on the burning characteristics of a model gasoline/biofuel mixture consisting of n-heptane and iso-butanol. A droplet burning in an environment that minimises the influence of buoyant and forced convective flows in the standard atmosphere is used to promote one-dimensional gas transport to facilitate numerical modelling of the droplet burning process. The numerical model includes a detailed combustion kinetic mechanism, unsteady gas and liquid transport, multicomponent diffusion inside the droplet, variable properties, and non-luminous radiative heat transfer from the flame. The numerical simulation was validated by experimental measurements in the standard atmosphere which showed good agreement with the evolutions of droplet and flame diameters. The iso-butanol concentration had a strong effect on formation of particulates. Above ~20% (volume) iso-butanol, flame luminosity was significantly diminished anddecreased with increasing iso-butanol concentration, while CO2 emissions as a representative greenhouse gas were not strongly influenced by the iso-butanol loading. The soot shell was located near a 1350 K isotherm for concentrations up to 20% (volume) iso-butanol, suggesting this value as a possible soot inception temperature for the mixture droplet. The combustion rate decreased with increasing iso-butanol concentration which was attributed to iso-butanol's higher liquid density. No evidence of a low temperature burning regime, or of extinction, was found (in experiments and simulations) for the small droplet sizes investigated.

Dalili, Alireza↗
NASA NTRS2011

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.↗
DOE OSTI2018

A class of efficient difference method for time fractional reaction–diffusion equation

The fractional reaction–diffusion equation has profound physical background and theoretical connotation, and its numerical methods are of great scientific significance and practical value. In this paper, a class of explicit–implicit (E–I) difference method and implicit–explicit (I–E) difference method are proposed for one-dimensional time fractional reaction–diffusion equation. The methods are constructed by combining classical explicit scheme and classical implicit scheme. The existence, uniqueness and convergence of solutions are given for E–I and I–E schemes. Theoretical analysis and numerical experiments show that the E–I and I–E schemes are unconditionally stable, with second-order spatial accuracy and 2−α order temporal accuracy. Compared with the classical implicit difference method with the similar accuracy, the computation time of the presented methods is improved by nearly 41%. It is shown that the E–I difference method and I–E difference method are effective for solving the time fractional reaction–diffusion equation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Enhancement of diffusive transport in oscillatory flows

The theory of transport of passive scalars in oscillatory flows is reexamined. The differences between transport in standing and traveling waves are emphasized. Both Lagrangian and Eulerian diffusivities are calculated, and the conditions for their applicability are discussed. Numerical simulations are conducted to understand the expulsion of gradients from time-dependent eddies and the resulting transport. The results indicate that it is the Eulerian diffusivity that is of primary relevance for describing enhanced transport on spatial scales larger than that of the eddies.

Knobloch, E.↗
DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Towards Unifying the Planetary Boundary Layer and Shallow Convection in CAM5 with the Eddy-Diffusivity/Mass-Flux Approach

The modular structure of the boundary layer and convection parameterizations in atmospheric models have long been affecting the numerical representation of subgrid-scale motions and their mutual interactions. A promising alternative, the eddy-diffusivity/mass-flux approach (EDMF), has the potential for unifying the existing formulations into a consistent scheme and improving some of the long-standing issues. This study documents a step towards developing such a unified approach by implementing a stochastic multi-plume EDMF scheme into the Community Atmosphere Model (Version 5.0). Its performance in single-column mode is evaluated against the control parameterization and large-eddy simulation (LES) for two benchmark cases: marine and continental shallow convection. Overall, the results for the two parameterizations agree well with each other and with LES in terms of mean profiles of moist conserved variables and their vertical fluxes, as well as the updraft properties. However, systematic differences between the two schemes, especially for transient continental convection, are also documented. Using EDMF helps improve some of the parameterized features of shallow convection. In particular, for the highest tested vertical resolution, the EDMF cloud base and top heights and the vertical fluxes of energy and water are remarkably close to LES.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2014

First-Order Hyperbolic System Method for Time-Dependent Advection-Diffusion Problems

A time-dependent extension of the first-order hyperbolic system method for advection-diffusion problems is introduced. Diffusive/viscous terms are written and discretized as a hyperbolic system, which recovers the original equation in the steady state. The resulting scheme offers advantages over traditional schemes: a dramatic simplification in the discretization, high-order accuracy in the solution gradients, and orders-of-magnitude convergence acceleration. The hyperbolic advection-diffusion system is discretized by the second-order upwind residual-distribution scheme in a unified manner, and the system of implicit-residual-equations is solved by Newton's method over every physical time step. The numerical results are presented for linear and nonlinear advection-diffusion problems, demonstrating solutions and gradients produced to the same order of accuracy, with rapid convergence over each physical time step, typically less than five Newton iterations.

Mazaheri, Alireza↗
NASA NTRS1995

A case study of view-factor rectification procedures for diffuse-gray radiation enclosure computations

The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Taylor, Robert P.↗
DOE OSTI2025

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗
DOE OSTI2018

Convective carbon dioxide dissolution in a closed porous medium at low pressure

Motivated by the persistence of natural carbon dioxide ( $\text{CO}_{2}$ ) fields, we investigate the convective dissolution of $\text{CO}_{2}$ at low pressure (below 1 MPa) in a closed system, where the pressure in the gas declines as convection proceeds. This introduces a negative feedback that reduces the convective dissolution rate even before the brine becomes saturated. We analyse the case of an ideal gas with a solubility given by Henry’s law, in the limits of very low and very high Rayleigh numbers. The equilibrium state in this system is determined by the dimensionless dissolution capacity, $\unicode[STIX]{x1D6F1}$ , which gives the fraction of the gas that can be dissolved into the underlying brine. Analytic approximations of the pure diffusion problem with $\unicode[STIX]{x1D6F1}>0$ show that the diffusive base state is no longer self-similar and that diffusive mass transfer declines rapidly with time. Direct numerical simulations at high Rayleigh numbers show that no constant flux regime exists for $\unicode[STIX]{x1D6F1}>0$ ; nevertheless, the quantity $F/C_{s}^{2}$ remains constant, where $F$ is the dissolution flux and $C_{s}$ is the dissolved concentration at the top of the domain. Simple mathematical models are developed to predict the evolution of $C_{s}$ and $F$ for high-Rayleigh-number convection in a closed system. The negative feedback that limits convection in closed systems may explain the persistence of natural $\text{CO}_{2}$ accumulations over millennial time scales.

Mechanics↗