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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Pressure falloff testing to characterize CO 2 plume and dry-out zone during CO 2 injection in saline aquifers

The fluid mobilities and the lateral extent of the CO 2 plume and the dry-out zone are among important unknowns required to effectively manage CO 2 storage in saline aquifers projects. These unknowns can be derived from pressure falloff testing which has been widely used to characterize the subsurface. This paper provides a three-region analytical solution for falloff pressure during the infinite-acting period of CO 2 injection in saline aquifers. The solution is derived using Laplace transformation technique of the governing diffusivity equations and the accompanying initial and boundary conditions. The exact solution is validated against numerical results obtained from a black-oil reservoir simulator for two cases. The solution is approximated into simple forms in real-time domain considering different time periods. Then, a graphical semi-log analysis technique based on the approximate solutions is provided as an interpretation tool of falloff tests where the fluid mobilities and the extents of dry-out zone and CO 2 plume can be inferred. The interpretation technique is applied to synthetic falloff test data to investigate its potential and reliability.

02 PETROLEUM↗
DOE OSTI2021

Pulsating Poiseuille flow of a cement slurry

In this paper, we investigate the pulsatile flow of a cement slurry in a pipe. The constitutive relation for the viscous stress tensor is based on the power-law model, where the shear viscosity not only depends on the shear rate but also the volume fraction of the cement particles. To solve for the volume fraction field, a convection-diffusion equation is used. The dimensionless form of the governing equations and the boundary conditions are solved numerically at different pulsatile cycles. Further, a parametric study is conducted to study the effect of different dimensionless numbers at different cycles on the velocity and the volume fraction profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Self-induced spin-orbit torques in metallic ferromagnets

Herein we present a phenomenological theory of spin-orbit torques in a metallic ferromagnet with spin-relaxing boundaries. The model is rooted in the coupled diffusion of charge and spin in the bulk of the ferromagnet, where we account for the anomalous Hall effects as well as the anisotropic magnetoresistance in the corresponding constitutive relations for both charge and spin sectors. The diffusion equations are supplemented with suitable boundary conditions reflecting the spin-sink capacity of the environment. In inversion-asymmetric heterostructures, the uncompensated spin accumulation exerts a dissipative torque on the order parameter, giving rise to a current-dependent linewidth in the ferromagnetic resonance with a characteristic angular dependence. We compare our model to recent spin-torque ferromagnetic resonance measurements, illustrating how rich self-induced spin-torque phenomenology can arise even in simple magnetic structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2021

Investigation of neutron noise in a micro-scale, natural circulation molten salt fission battery system

The Molten Salt Reactor (MSR) concept is a rapidly evolving Generation IV design that has recently attracted favorable attention due to the potential for reducing waste generation, realizing passive safety features, and seizing on the opportunity for cost effective economics. A specific novel micro-MSR, natural circulation, battery design concept involves placing all primary components within a single reactor vessel containment without the need of primary forced pumping. This small modular, integral design presents potential cost savings while producing safety, reliable, and transportable carbon-free power for decades. This investigation evaluates the neutron noise induced by density and flow fluctuations in a natural circulation MSR battery concept being developed at the University of Idaho. This study finds numerical solutions to the one-dimensional, one-group, coupled diffusion equations for a bare, homogeneous core, and evaluates neutron flux and delayed neutron precursor concentration noise due to core flow and fuel salt density perturbation sources. Noise analysis shows both point kinetic and space-dependent behavior is present despite the small size, low-flow, closely-coupled nature of the natural circulation MSR concept. Analysis shows that the low flow of the convective system results in noise behavior typically representative of a larger forced-circulation MSR system. Finally, noise techniques may also be useful in MSR core diagnostics as a non-invasive, low-cost radiological option to traditional monitoring schemes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Differentiating Thermal Conductances at Semiconductor Nanocrystal/Ligand and Ligand/Solvent Interfaces in Colloidal Suspensions

Infrared-pump, electronic-probe (IPEP) spectroscopy is used to measure heat flow into and out of CdSe nanocrystals suspended in an organic solvent, where the surface ligands are initially excited with an infrared pump pulse. Subsequently, the heat is transferred from the excited ligands to the nanocrystals and in parallel to the solvent. Parallel heat transfer in opposite directions uniquely enables us to differentiate the thermal conductances at the nanocrystal/ligand and ligand/solvent interfaces. Using a novel solution to the heat diffusion equation, we fit the IPEP data to find that the nanocrystal/ligand conductances range from 88 to 135 MW m –2 K –1 and are approximately 1 order of magnitude higher than the ligand/solvent conductances, which range from 7 to 26 MW m –2 K –1 . Transient nonequilibrium molecular dynamics (MD) simulations of nanocrystal suspensions agree with IPEP data and show that ligands bound to the nanocrystal by bidentate bonds have more than twice the per-ligand conductance as those bound by monodentate bonds.

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

Calculating RF current condensation with consistent ray-tracing and island heating

Finally, by exploiting the nonlinear amplification of the power deposition of RF waves, current condensation promises new pathways to the stabilization of magnetic islands. We present a numerical analysis of current condensation, coupling a geometrical optics treatment of wave propagation and damping to a thermal diffusion equation solver in the island. Taking into account the island geometry and relativistic damping, previous analytical theory can be made more precise and specific scenarios can be realistically predicted. With this more precise description, bifurcations and associated hysteresis effects could be obtained in an ITER-like scenario at realistic parameter values. Moreover, it is shown that dynamically varying the RF wave launching angles can lead to hysteresis and help to avoid the nonlinear shadowing effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Wide bandwidth frequency-domain thermoreflectance: Volumetric heat capacity, anisotropic thermal conductivity, and thickness measurements

A frequency-domain thermoreflectance (FDTR) system with a frequency range of 1 Hz to 75 MHz is presented. The wide bandwidth of pump modulation frequency enables enhanced simultaneous thermophysical property characterization of bulk and thin-film materials. The wide bandwidth FDTR system is demonstrated with simultaneous measurements of in-plane and cross-plane thermal conductivities of sapphire and muscovite mica, thickness and thermal conductivity of gold/titanium thin films, and isotropic thermal conductivity and volumetric heat capacity of lithium niobate and silicon. Thermophysical properties measured with FDTR are compared to literature values or independent measurements for verification. Finally, at low frequencies, a negligible error was achieved in the numerical integration of the heat diffusion equation with a Gauss–Legendre quadrature method solved with 400 integration points and an upper integration limit of 16weff, where weff is the effective pump–probe 1/e2 radius.

Ziade, Elbara (ORCID:0000000249020497)↗
DOE OSTI2022

Thermal model for time-domain thermoreflectance experiments in a laser-flash geometry

Time-domain thermoreflectance (TDTR) is a well-established pump–probe method for measuring thermal conductivity and interface conductance of multilayers. Interpreting signals in a TDTR experiment requires a thermal model. In standard front/front TDTR experiments, both pump and probe beams typically irradiate the surface of a multilayer. As a result, existing thermal models for interpreting thermoreflectance experiments assume that the pump and probe beams both interact with the surface layer. Here, we present a frequency-domain solution to the heat-diffusion equation of a multilayer in response to nonhomogeneous laser heating. This model allows analysis of experiments where the pump and probe beams irradiate opposite sides of a multilayer. We call such a geometry a front/back experiment to differentiate such experiments from standard TDTR experiments. As an example, we consider a 60nm amorphous Si film. We consider how signals differ in a front/front vs front/back geometry and compare thermal model predictions to experimental data.

Peng, Wanyue↗
DOE OSTI2025

Designing a validation experiment for radio frequency condensation

Abstract Theoretical studies have suggested that nonlinear effects can lead to ‘radio frequency (RF) condensation’, where an initially broad current profile can coalesce in islands when they reach sufficient width. In suitable conditions, RF condensation can ‘self-focus’ the driven current to the center of an island, improving stabilization efficiency and reducing control complexity. In unsuitable conditions, the effect can prematurely deplete the RF energy before it reaches the island center, impairing stabilization. It is predicted that the RF condensation effect can significantly impact reactor-scale tokamaks. This paper presents a set of simulations investigating the conditions under which RF condensation might be encountered in present-day tokamaks. For concreteness, the calculations use equilibrium reconstructions for two shots from DIII-D and AUG. The Current Condensation Amid Magnetic Islands (OCCAMI) simulation code has been used for this investigation. The code takes as its input a numerically specified axisymmetric EFIT equilibrium solution, and it perturbatively constructs a 3D field with an island embedded at the appropriate rational surface. In the OCCAMI code, the GENRAY code is used for ray tracing and for calculating the power deposition along a ray trajectory, and GENRAY is coupled self-consistently to a solution of the thermal diffusion equation in the island. The simulation results described in the paper illuminate the conditions required for experimental validation of the theory of RF condensation. The simulations also provide an explanation of why the effect was not noticed in experiments prior to the publication of theoretical papers on the subject.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Coupled Charge and Radiation Transport Processes in Thermophotovoltaic and Thermoradiative Cells

Accurate modeling of charge transport and both thermal and luminescent radiation is crucial to the understanding and design of radiative thermal energy converters. Charge-carrier dynamics in semiconductors are well-described by the Poisson-drift-diffusion equations, and thermal radiation in emitter-absorber structures can be computed using multilayer fluctuational electrodynamics. These two types of energy flows interact through radiation absorption or luminescence and charge-carrier generation or recombination. However, past research has typically only assumed limited interaction, with thermal radiation absorption as an input for charge-carrier models to predict device performance. To examine this assumption, we develop a fully coupled iterative model of charge and radiation transport in semiconductor devices, and we use our model to analyze near-field and far-field GaSb thermophotovoltaic and thermoradiative systems. By comparing our results with past methods that do not consider cross-influences between charge and radiation transport, we find that a fully coupled approach is necessary to accurately model photon recycling and near-field enhancement of external luminescence. As these effects can substantially alter device performance, our modeling approach can aid in the design of efficient thermophotovoltaic and thermoradiative systems.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Charge Transport and Space-Charge Formation in Cd 1 - x Zn x Te 1 - y Se y Radiation Detectors

The electron- and hole-transport properties in cadmium zinc telluride selenide (CZTS) crystals are studied using a laser-induced transient-current technique with pulsed and dc bias. The internal electric field profile and velocity of surface recombination are determined by Monte Carlo simulations of electron and hole transient currents combined with a numerical solution of the drift-diffusion equation coupled with Poisson’s equation. Electron and hole drift mobilities of μe = 830 cm 2 /Vs and μh = 40 cm 2 /Vs, respectively, are determined. We also develop a simple technique for evaluating surface recombination directly from measured current waveforms without the need for numerical simulation. The good quality of the prepared detector at pulsed bias, with electron- and hole-mobility-lifetime products of (μτ)e = 1.9 × 10 -3 cm 2 /V and (μτ)h = 1.4 × 10 -4 cm 2 /V, respectively, are observed. The formation of a positive space charge, originating from hole injection combined with a recombination level, is found. We observe a significant position dependence of the lifetime of electrons and holes in dc bias due to hole injection. The experiment is successfully fitted by a simple model dominated by a single deep recombination level with an energy of E t =E C -0.73eV; concentration of 7.3 × 10 11 cm -3 ; and electron- and hole-capture cross sections of 3.5 × 10 -14 cm 2 and 6.5 × 10 -14 cm 2 , respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Onset of anharmonicity and thermal conductivity in SnSe

The anharmonicity in SnSe is investigated through the analysis of moments of 119Sn nuclear resonant inelastic x-ray scattering and ab initio molecular dynamics calculations. Experimental evidences show that the anharmonic behavior started around 300 K, substantially lower than the usually suggested structural transition at 800 K. Both theory and experiments reveal substantial lifetime broadening and frequency renormalization of the optical phonons. Additionally, thermal conductivities calculated from the temporal energy moment using the Einstein diffusion equation are in good agreement with previous experiments. The abrupt increase of the thermal power near 800 K is driven by an electronic factor and not by the enhanced anharmonicity due to structural change.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Kinetic-energy dissipation and fluctuations in strongly damped heavy-ion collisions within the stochastic mean-field approach

Microscopic mean-field approaches have been successful in describing the most probable reaction outcomes in low-energy heavy-ion reactions. However, those approaches are known to severely underestimate dispersions of observables around the average values that has limited their applicability. Recently it has been shown that a quantal transport approach based on the stochastic mean-field (SMF) theory significantly improves the description, while its application has been limited so far to fragment mass and charge dispersions. In this work, we extend the quantal transport approach based on the SMF theory for relative kinetic energy dissipation and angular momentum transfer in low-energy heavy-ion reactions. Based on the SMF concept, analytical expressions are derived for the radial and tangential friction and associated diffusion coefficients. Those quantal transport coefficients are calculated microscopically in terms of single-particle orbitals within the time-dependent Hartree-Fock (TDHF) approach. As the first application of the proposed formalism, we consider the radial linear momentum dispersion, neglecting the coupling between radial and angular momenta. We analyze the total kinetic energy (TKE) distribution of binary reaction products in the 136 Xe + 208 Pb reaction at E c.m. =526 MeV and compare with experimental data. From time evolution of single-particle orbitals in TDHF, the radial diffusion coefficient is computed on a microscopic basis, while a phenomenological treatment is introduced for the radial friction coefficient. By solving the quantal diffusion equation for the radial linear momentum, the dispersion of the radial linear momentum is obtained, from which one can construct the TKE distribution. We find that the calculations provide a good description of the TKE distribution for strongly damped events with large energy losses, TKEL ≳ 150 MeV. However, the calculations underestimate the TKE distribution for smaller energy losses. Further studies are needed to improve the technical details of calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗