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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Designing High-Efficiency Thin Silicon Solar Cells Using Parabolic-Pore Photonic Crystals

We demonstrate the efficacy of wave-interference-based light trapping and carrier transport in parabolic-pore photonic-crystal, thin-crystalline silicon (c -Si) solar cells to achieve above 29% power conversion efficiencies. Using a rigorous solution of Maxwell’s equations through a standard finite-difference time domain scheme, we optimize the design of the vertical-parabolic-pore photonic crystal (PhC) on a 10 μm -thick (c -Si) solar cell to obtain a maximum achievable photocurrent density (MAPD) of 40.6 mA/cm 2 beyond the ray-optical, Lambertian light-trapping limit. For a slanted-parabolic-pore PhC that breaks x-y symmetry, improved light trapping occurs due to better coupling into parallel-to-interface refraction modes. We achieve the optimum MAPD of 41.6 mA/cm 2 for a tilt angle of 10° with respect to the vertical axis of the pores. This MAPD is further improved to 41.72 mA/cm 2 by introducing a 75-nm SiO 2 antireflective coating on top of the solar cell. We use this MAPD and the associated charge-carrier generation profile as input for a numerical solution of Poisson’s equation coupled with semiconductor drift-diffusion equations using a Shockley-Read-Hall and Auger recombination model. Using experimentally achieved surface recombination velocities of 10 cm/s, we identify semiconductor doping profiles that yield power conversion efficiencies over 29%. Practical considerations of additional upper-contact losses suggest efficiencies close to 28%. This improvement beyond the current world record is largely due to an open-circuit voltage approaching 0.8 V enabled by reduced bulk recombination in our thin silicon architecture while maintaining a high short-circuit current through wave-interference-based light trapping.

14 SOLAR ENERGY↗
DOE OSTI2022

Mesoscopic finite-size effects of unconventional electron transport in PdCoO 2

A wide range of unconventional transport phenomena has recently been observed in single-crystal delafossite metals. Here, we present a theoretical framework to elucidate electron transport using a combination of first-principles calculations and numerical modeling of the anisotropic Boltzmann transport equation. Using PdCoO 2 as a model system, we study different microscopic electron and phonon scattering mechanisms and establish the mean free path hierarchy of quasiparticles at different temperatures. We treat the anisotropic Fermi surface explicitly to numerically obtain experimentally-accessible transport observables, which bridge between the “diffusive,” “ballistic,” and “hydrodynamic” transport regime limits. Further, we illustrate that the distinction between the “quasiballistic” and “quasihydrodynamic” regimes is challenging and often needs to be quantitative in nature. From first-principles calculations, we populate the resulting transport regime plots and demonstrate how the Fermi surface orientation adds complexity to the observed transport signatures in micrometer-scale devices. Our work provides key insights into microscopic interaction mechanisms on open hexagonal Fermi surfaces and establishes their connection to the macroscopic electron transport in finite-size channels.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Interaction between two conducting spheres in weakly ionized collisional plasma

Charging of two conducting spheres in a weakly ionized collisional plasma flow is considered. The spheres are arranged along the flow, and the plasma is assumed to consist of two ion species with the charges equal in magnitude but opposite in sign. The problem is analyzed with allowance for the external electric field, charging of the spheres due to the sedimentation of plasma ions on them, the fields of the sphere charges, the space charge field, and the processes of recombination and molecular diffusion. The interaction between the spheres and plasma is studied by numerically solving a time-dependent problem in a bispherical coordinate system by the finite difference method. The steady-state values of the sphere charges and the distributions of the space charge and ion densities in the ambient plasma are found as functions of the plasma parameters and the distance between the spheres. The electrostatic forces acting on the spheres are determined, and the effects of the external field, the space charge fields, and the fields of the sphere charges are comparatively analyzed. It is shown that, for the considered plasma parameters, the main electrostatic effect in the interaction between two spheres is their mutual approach in the external field due to the difference in their charges (one sphere catches up with the other). Due to the friction force with the neutral gas, this mutual approach is much slower than all other processes in the system. For widely spaced spheres, the results coincide with the solution obtained previously for a solitary sphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

An Efficient Parallel-in-Time Method for Optimization with Parabolic PDEs

To solve optimization problems with parabolic PDE constraints, often methods working on the reduced objective functional are used. They are computationally expensive due to the necessity of solving both the state equation and a backward-in-time adjoint equation to evaluate the reduced gradient in each iteration of the optimization method. In this study, we investigate the use of the parallel-in-time method PFASST in the setting of PDE-constrained optimization. In order to develop an efficient fully time-parallel algorithm, we discuss different options for applying PFASST to adjoint gradient computation, including the possibility of doing PFASST iterations on both the state and the adjoint equations simultaneously. Here, we also explore the additional gains in efficiency from reusing information from previous optimization iterations when solving each equation. Numerical results for both a linear and a nonlinear reaction-diffusion optimal control problem demonstrate the parallel speedup and efficiency of different approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Direct Simulations of Coupled Transport and Reaction on Nano-Scale X-Ray Computed Tomography Images of Platinum Group Metal-Free Catalyst Cathodes

The nano/micro-scale geometry of polymer electrolyte fuel cell (PEFC) catalyst layers critically affects cell performance. The small length scales and complex structure of these composite layers make it challenging to analyze cell performance and physics at the particle scale by experiment. We present a computational method to simulate transport and chemical reaction phenomena at the pore/particle-scale and apply it to a PEFC cathode with platinum group metal free (PGM-free) catalyst. Here, we numerically solve the governing equations for the physics with heterogeneous oxygen diffusion coefficient and proton conductivity evaluated using the actual electrode structure and ionomer distribution obtained using nano-scale resolution X-ray computed tomography (nano-CT). Using this approach, the oxygen concentration and electrolyte potential distributions imposed by the oxygen reduction reaction are solved and the impact of the catalyst layer structure on performance is evaluated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemical Characterization of Dissolved Oxides in Molten FLiNaK

Molten fluoride salts have long been the subject of investigation given their application in a variety of industrial processes. Oxides are a common exogeneous impurity within these systems, and can lead to many operational and regulatory concerns. In recent studies, the electro-oxidation of dissolved metallic oxides has been attributed to a single step oxygen evolution reaction despite voltametric evidence suggesting otherwise. Here, we first use square wave voltammetry to confirm that the oxidation of dissolved oxides in fluorides proceeds via a two-electron transfer reaction and then demonstrate that the presence of a peroxide-mediated redox reaction better explains the electrochemical data. Numerical simulations were then used to support the determination of improved diffusion coefficients to enable electroanalytical measurements of the oxide concentration across at a variety of temperatures. A linear Arrhenius relationship was observed when the experimental data was corrected for ohmic resistance effects. Concentration measurements using the corrected square wave voltammetry data demonstrated high accuracy across a large range of concentrations, while the uncorrected data showed plateauing and high errors. In total, this work serves to rectify previously misconstrued electrochemical data of oxides in molten salts and demonstrate how accurate, in situ concentration measurements can be achieved in real-world systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Determining Exact RANS Operators with the Macroscopic Forcing Method (Final Report)

This report contains a compilation of key results and findings from the ACT project “Determining Exact RANS Operators with the Macroscopic Forcing Method.” The Macroscopic Forcing Method (MFM), a numerical tool for determining closure operators, is used to measure eddy diffusivity moments in Rayleigh-Taylor (RT) instability. It is first applied to low-Atwood 2D RT instability; that work is then extended to 3D RT at different finite Atwood numbers. It is found that nonlocality is important for modeling the mean scalar transport closure operator in RT mixing. Additionally, work is done to improve the statistical convergence of MFM for chaotic problems like RT mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

UO 2 Grain Growth: Developing Phase Field Models for Pore Dragging, Solute Dragging and Anisotropic Grain Boundary Energies

A detailed phase field model for the effect of pore drag on grain growth kinetics was implemented in MARMOT. The model takes into consideration both the curvature-driven grain boundary motion and pore migration by surface diffusion. As such, the model accounts for the interaction between pore and grain boundary kinetics, which tends to retard the grain growth process. Our 2D and 3D simulations demonstrate that the model capture all possible pore-grain boundary interactions proposed in theoretical models. For high enough surface mobility, the pores move along with the migrating boundary as a quasi-rigid-body, albeit hindering its migration rate compared to the pore-free case. For less mobile pores, the migrating boundary can separate from the pores. For the pore-controlled grain growth kinetics, the model predicts a strong dependence of the growth rate on the number of pores, pore size, and surface diffusivity in agreement with theroretical models. An evolution equation for the grain size that includes these parameters was derived and showed to agree well with numerical solution. It shows a smooth transition from boundary-controlled kinetics to pore-controlled kinetics as the surface diffusivity decreases or the number of pores or their size increases. This equation can be utilized in BISON to give accurate estimate for the grain size evolution. This will be accomplished in the near future. The effect of solute drag and anisotropy of grain boundary on grain growth will be investigated in future studies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An on-demand, drop-on-drop method for studying enzyme catalysis by serial crystallography

Serial femtosecond crystallography has opened up many new opportunities in structural biology. In recent years, several approaches employing light-inducible systems have emerged to enable time-resolved experiments that reveal protein dynamics at high atomic and temporal resolutions. However, very few enzymes are light-dependent, whereas macromolecules requiring ligand diffusion into an active site are ubiquitous. In this work we present a drop-on-drop sample delivery system that enables the study of enzyme-catalyzed reactions in microcrystal slurries. The system delivers ligand solutions in bursts of multiple picoliter-sized drops on top of a larger crystal-containing drop inducing turbulent mixing and transports the mixture to the X-ray interaction region with temporal resolution. We demonstrate mixing using fluorescent dyes, numerical simulations and time-resolved serial femtosecond crystallography, which show rapid ligand diffusion through microdroplets. The drop-on-drop method has the potential to be widely applicable to serial crystallography studies, particularly of enzyme reactions with small molecule substrates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

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

An Investigation into Solution Verification for CFD-DEM

This report presents the study of the convergence behavior of the computational fluid dynamicsdiscrete element method (CFD-DEM) method, specifically National Energy Technology Laboratory’s (NETL) open source MFiX code (MFiX-DEM) with a diffusion based particle-tocontinuum filtering scheme. In particular, this study focused on determining if the numerical method had a solution in the high-resolution limit where the grid size is smaller than the particle size. To address this uncertainty, fixed particle beds of two primary configurations were studied: i) fictitious beds where the particles are seeded with a random particle generator, and ii) instantaneous snapshots from a transient simulation of an experimentally relevant problem. Both problems considered a uniform inlet boundary and a pressure outflow. The CFD grid was refined from a few particle diameters down to 1/6 th of a particle diameter. The pressure drop between two vertical elevations, averaged across the bed cross-section was considered as the system response quantity of interest. A least-squares regression method was used to extrapolate the grid-dependent results to an approximate “grid-free” solution in the limit of infinite resolution. The results show that the diffusion based scheme does yield a converging solution. However, the convergence is more complicated than encountered in simpler, single-phase flow problems showing strong oscillations and, at times, oscillations superimposed on top of globally non-monotonic behavior. The challenging convergence behavior highlights the importance of using at least four grid resolutions in solution verification problems so that (over-determined) regression-based extrapolation methods may be applied to approximate the grid-free solution. The grid-free solution is very important in solution verification and VVUQ exercise in general as the difference between it and the reference solution largely determines the numerical uncertainty. By testing different randomized particle configurations of the same general problem (for the fictitious case) or different instances of freezing a transient simulation, the numerical uncertainties appeared to be on the same order of magnitude as ensemble or time averaging uncertainties. By testing different drag laws, almost all cases studied show that model form uncertainty in this one, very important closure relation was larger than the numerical uncertainty, at least with a reasonable CFD grid, roughly five particle diameters. In this study, the diffusion width (filtering length scale) was mostly set at a constant of six particle diameters. A few exploratory tests were performed to show that similar convergence behavior was observed for diffusion widths greater than approximately two particle diameters. However, this subject was not investigated in great detail because determining an appropriate filter size is really a validation question which must be determined by comparison to experimental or highly accurate numerical data. Future studies are being considered targeting solution verification of transient simulations as well as validation of the filter size with direct numerical simulation data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Accretion of Gas Giants Constrained by the Tidal Barrier

After protoplanets have acquired sufficient mass to open partial gaps in their natal protostellar disks, residual gas continues to diffuse onto horseshoe streamlines under the effect of viscous dissipation, and to meander in and out of the planets' Hill sphere. Within the Hill sphere, the horseshoe streamlines intercept gas flow in circumplanetary disks. The host star's tidal perturbation induces a barrier across the converging streamlines' interface. Viscous transfer of angular momentum across this tidal barrier determines the rate of mass diffusion from the horseshoe streamlines onto the circumplanetary disks, and eventually the accretion rate onto the protoplanets. Here, we carry out a series of numerical simulations to test the influence of this tidal barrier on superthermal planets. In weakly viscous disks, the protoplanets' accretion rate steeply decreases with their masses above the thermal limit. As their growth timescale exceeds the gas depletion timescale, their masses reach asymptotic values comparable to that of Jupiter. In relatively thick and strongly viscous disks, protoplanets' asymptotic masses exceed several times that of Jupiter. Two-dimensional numerical simulations show that such massive protoplanets strongly excite the eccentricity of nearby horseshoe streamlines, destabilize orderly flow, substantially enhance the diffusion rate across the tidal barrier, and elevate their growth rate until their natal disk is severely depleted. In contrast, eccentric streamlines remain stable in three-dimensional simulations. Based on the upper falloff in the observed mass distribution of known exoplanets, we suggest that their natal disks had relatively low viscosity (α ~ 10 –3 ) and modest thickness (H/R ~ 0.03–0.05).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Interdiffusion and phase formation in the Fe–TiO{sub 2} thin-film system

The interaction of magnetron-sputtered metal iron with titanium-oxide films upon isothermal vacuum annealing is studied by X-ray phase analysis, secondary-ion mass spectrometry, atomic-force microscopy, and mathematical simulation. A mechanism for the formation of complex oxides at grain boundaries is suggested. The mechanism is based on the reaction diffusion of metal iron into titanium oxide. A quantitative model of reaction interdiffusion in two-layer polycrystalline metal–oxide film systems with limited component solubility is developed. From numerical analysis of the experimental distributions of the metal concentrations in the Fe–TiO{sub 2} film system, the individual diffusion coefficients are determined. It is found that, under the conditions of vacuum annealing at 1073 K, the diffusion coefficients of iron and titanium are 8.0 × 10{sup –13} and 3.0 × 10{sup –15} cm{sup 2} s{sup –1}, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A conservative nonlocal convection–diffusion model and asymptotically compatible finite difference discretization

Here in this paper, we first propose a nonlocal convection–diffusion model, in which the convection term is constructed in a special upwind manner so that mass conservation and maximum principle are maintained in any space dimension. The well-posedness of the proposed nonlocal model and its convergence to the classical local convection–diffusion model are established. A quadrature-based finite difference discretization is then developed to numerically solve the nonlocal problem and it is shown to be consistent and unconditionally stable. We further demonstrate that the numerical scheme is asymptotically compatible, that is, the approximate solutions converge to the exact solution of the corresponding local problem when δ → 0 and h → 0. Numerical experiments are also performed to complement the theoretical analysis.

42 ENGINEERING↗
DOE OSTI2018

A Legendre spectral quadrature tau method for the multi-term time-fractional diffusion equations

In this study, a Legendre spectral tau method is revisited to handle the multi-term time-fractional diffusion equations (MTT-FDEs). An error estimate and rigorous convergence analysis are carried out using some critical theorems. The applicability and accuracy of the solution method are demonstrated by a numerical example to verify the theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Equilibration in fermionic systems

Highlights: •Local equilibration in finite fermion systems is based on a nonlinear diffusion equation. •The fermionic diffusion equation is solved analytically in closed form. •The exact solutions agree with the numerical results. •Particle-number and energy conservation are shown to be fulfilled. •The model is applied to hadrons at low energies, and relativistic quarks. -- Abstract: The time evolution of a finite fermion system towards local statistical equilibrium is investigated using analytical solutions of a nonlinear partial differential equation that had been derived earlier from the Boltzmann collision term. The solutions of this fermionic diffusion equation are rederived in closed form, evaluated exactly for simplified initial conditions, and applied to hadron systems at low energies in the MeV-range, as well as to quark systems at relativistic energies in the TeV-range where antiparticle production is abundant. Conservation laws for particle number including created antiparticles, and for the energy are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Transport studies in high-performance field reversed configuration plasmas

A significant improvement of field reversed configuration (FRC) lifetime and plasma confinement times in the C-2 plasma, called High Performance FRC regime, has been observed with neutral beam injection (NBI), improved edge stability, and better wall conditioning [Binderbauer et al., Phys. Plasmas 22, 056110 (2015)]. A Quasi-1D (Q1D) fluid transport code has been developed and employed to carry out transport analysis of such C-2 plasma conditions. The Q1D code is coupled to a Monte-Carlo code to incorporate the effect of fast ions, due to NBI, on the background FRC plasma. Numerically, the Q1D transport behavior with enhanced transport coefficients (but with otherwise classical parametric dependencies) such as 5 times classical resistive diffusion, classical thermal ion conductivity, 20 times classical electron thermal conductivity, and classical fast ion behavior fit with the experimentally measured time evolution of the excluded flux radius, line-integrated density, and electron/ion temperature. The numerical study shows near sustainment of poloidal flux for nearly 1 ms in the presence of NBI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Processes of diffusion in a limited ensemble of charged particles in a static magnetic field

The dynamics of a limited ensemble of charged particles in a spatially uniform static magnetic field has been studied analytically and numerically. Calculations have been performed for particles with various masses and charges in a wide range of the parameters of the systems under study. It has been shown that the transverse diffusion coefficient in an ensemble of particles in the magnetic field for weakly nonideal systems is described by the Townsend formula. Analytical estimates have been proposed for the analysis of the effective diffusion radius of the ensemble of particles owing to their thermal motion. The proposed relations have been tested by the numerical simulation of the problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗