Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anisotropic 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 289 records · Page 16

Electromagnetic wave propagation based upon spectral-element methodology in dispersive and attenuating media

SUMMARY We build on mathematical equivalences between Maxwell’s wave equations for an electromagnetic medium and elastic seismic wave equations. This allows us to readily model Maxwell’s wave propagation in the spectral-element codes SPECFEM2D and SPECFEM3D, written for acoustic, viscoelastic and poroelastic seismic wave propagation, providing the ability to handle complex geometries, inherent to finite-element methods and retaining the strength of exponential convergence and accuracy due to the use of high-degree polynomials to interpolate field functions on the elements, characteristic to spectral-element methods (SEMs). Attenuation and dispersion processes related to the frequency dependence of dielectric permittivity and conductivity are also included using a Zener model, similar to shear attenuation in viscoelastic media or viscous diffusion in poroelastic media, and a Kelvin–Voigt model, respectively. Ability to account for anisotropic media is also discussed. Here, we limit ourselves to certain dielectric permittivity tensor geometries, in order to conserve a diagonal mass matrix after discretization of the system of equations. Doing so, simulation of Maxwell’s wave equations in the radar frequency range based on SEM can be solved using explicit time integration schemes well suited for parallel computation. We validate our formulation with analytical solutions. In 2-D, our implementation allows for the modelling of both a transverse magnetic (TM) mode, suitable for surface based reflection ground penetration radar type of applications, and a transverse electric (TE) mode more suitable for crosshole and vertical radar profiling setups. Two 2-D examples are designed to demonstrated the use of the TM and TE modes. A 3-D example is also presented, which allows for the full TEM solution, different antenna orientations, and out-of-plane variations in material properties.

58 GEOSCIENCES↗

Mechanical and Thermophysical Properties of 3D-Printed SiC-FY20

In the Transformational Challenge Reactor (TCR), the fuel blocks consist of an additively-manufactured silicon carbide (SiC) matrix and uranium nitride tristructural isotropic (UN TRISO) fuel particles, which are stacked to form fuel columns. The SiC matrix is manufactured using binderjet 3D printing followed by loading the TRISO fuel particles and the chemical vapor infiltration (CVI) process. Because the fuel matrix is a primary component of the TCR core and its response to mechanical and thermal loads during operation is one of the most influential factors on the integrity of TCR core, testing and evaluation have focused on producing mechanical and thermophysical properties data for the binderjet/CVI SiC. Mechanical and thermophysical properties were measured from various types of specimens printed for two or three orientations, which included equibiaxial flexural failure strength, elastic constants, thermal diffusivity and conductivity, density, and the coefficient of thermal expansion. Flexural failure strength datasets showed similar Weibull distributions regardless of sample variants including different orientations. The mean failure strengths of the 3D-printed SiC variants were in the range of 286–306 MPa, which are 22–27% lower than that of the CVD SiC. Thermophysical test results showed that specific heat and thermal expansion are not sensitive to the build directions of SiC samples, while thermal diffusivity is highly dependent on the build direction and can be correlated to the anisotropic character of the 3D-printed SiC. This report also includes discussions on the uniaxial tensile properties of the as-printed SiC before CVI and on ongoing efforts for irradiation effects studies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Cosmic-Ray Feedback on Bistable Interstellar Medium Turbulence

Abstract While cosmic rays ( E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Anisotropic thermal conductivity of thin polycrystalline oxide samples

This paper reports about the development of a modified laser-flash technique and relation to measure the in-plane thermal diffusivity of thin polycrystalline oxide samples. Thermal conductivity is then calculated with the product of diffusivity, specific heat and density. Design and operating features for evaluating in-plane thermal conductivities are described. The technique is advantageous as thin samples are not glued together to measure in-plane thermal conductivities like earlier methods reported in literature. The approach was employed to study anisotropic thermal conductivity in alumina sheet, textured kaolin ceramics and montmorillonite. Since it is rare to find in-plane thermal conductivity values for such anisotropic thin samples in literature, this technique offers a useful variant to existing techniques.

36 MATERIALS SCIENCE↗

Accelerating the solution of the S N equations with highly anisotropic scattering using the Fokker-Planck approximation

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. In this paper, we explore the idea of using the Fokker-Planck approximation as a preconditioner to accelerate the convergence of solution of forward-peaked transport problems. Finally, we observe that such preconditioning can successfully reduce the iteration count by up to four orders of magnitude, and solver runtime by up to three orders of magnitude when compared to diffusion synthetic acceleration.

42 ENGINEERING↗

Interfacial mixing in high-energy-density matter with a multiphysics kinetic model

We have extended here a recently developed multispecies, multitemperature Bhatnagar-Gross-Krook model [Haack et al., J. Stat. Phys. 168, 822 (2017)], to include multiphysics capabilities that enable modeling of a wider range of physical conditions. In terms of geometry, we have extended from the spatially homogeneous setting to one spatial dimension. In terms of the physics, we have included an atomic ionization model, accurate collision physics across coupling regimes, self-consistent electric fields, and degeneracy in the electronic screening. We apply the model to a warm dense matter scenario in which the ablator-fuel interface of an inertial confinement fusion target is heated, but for larger length and time scales and for much higher temperatures than can be simulated using molecular dynamics. Relative to molecular dynamics, the kinetic model greatly extends the temperature regime and the spatiotemporal scales over which we are able to model. In our numerical results we observe hydrogen from the ablator material jetting into the fuel during the early stages of the implosion and compare the relative size of various diffusion components (Fickean diffusion, electrodiffusion, and barodiffusion) that drive this process. We also examine kinetic effects, such as anisotropic distributions and velocity separation, in order to determine when this problem can be described with a hydrodynamic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microscale Effects from Global Hot Plasma Imagery

We have used a three-dimensional model of recovery phase storm hot plasmas to explore the signatures of pitch angle distributions (PADS) in global fast atom imagery of the magnetosphere. The model computes mass, energy, and position-dependent PADs based on drift effects, charge exchange losses, and Coulomb drag. The hot plasma PAD strongly influences both the storm current system carried by the hot plasma and its time evolution. In turn, the PAD is strongly influenced by plasma waves through pitch angle diffusion, a microscale effect. We report the first simulated neutral atom images that account for anisotropic PADs within the hot plasma. They exhibit spatial distribution features that correspond directly to the PADs along the lines of sight. We investigate the use of image brightness distributions along tangent-shell field lines to infer equatorial PADS. In tangent-shell regions with minimal spatial gradients, reasonably accurate PADs are inferred from simulated images. They demonstrate the importance of modeling PADs for image inversion and show that comparisons of models with real storm plasma images will reveal the global effects of these microscale processes.

Moore, T. E.↗

Meso-Scale Fuel Performance Modeling using the MARMOT Code

The in-pile behavior of light-water-reactor (LWR) fuels is governed by the physics from various time and spatial scales. Accordingly, multi-scale modeling and simulation are needed to accurately predict the fuel performance at the engineering scale. To improve the traditional fuel performance modeling, the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program targets at developing multi-scale, multiphysics modeling tools for in-pile fuel behavior prediction under various operating or accident scenarios. In traditional fuel performance modeling, the history of fuel is represented by a single parameter, fuel burn-up. Burn-up measures the energy that has been extracted from fuels. However, it is not an adequate measure of fuel history as it ignores the flux effect, which is critical for the irradiation-induced microstructural evolution in fuels. Moreover, the traditional fuel performance models are usually fitted to experimental data. Their capabilities of being extrapolated to new operating conditions and being extended to new types of fuels are strongly limited by the empirical or semi-empirical nature. To improve the burn-up based models, a microstructure based fuel approach, the MOOSE-BISON-MARMOT (MBM) tool-kit, is being developed under the fuel product line (FPL) of the NEAMS program for predictive fuel performance modeling. In the MBM approach, the evolution of fuel microstructure and changes in fuel properties are described by the MARMOT code. The transient fuel properties obtained from MARMOT are then fed into the BISON code for engineering scale fuel performance modeling. Both BISON and MARMOT are developed within the MOOSE framework, which was developed at Idaho National Laboratory. In this summary, the development and the capabilities of the MARMOT code will be introduced. In the past few years, several materials models have been developed including: (1) a grain growth model that describes the grain size evolution as a function of time, temperature and fuel porosity, considering the effects of impurity and porosity dragging and anisotropic grain boundary energy, (2) a fission gas release model that describes the fission gas diffusion in grain interior, fission gas bubble formation in grain interior and at grain boundaries, and fission gas release to the fuel-cladding gap, (3) a fracture model that describes the fracture initiation and propagation in grain interior and along grain boundaries, and (4) a thermal conductivity model that describes the change in thermal conductivity due to grain growth, gas bubble formation, fracture and radiation damage. Several other models are also being developed in MARMOT for hydride formation in Zr-based cladding, radiation damage and grain size subdivision. In most cases, these materials models are coupled with each other. For instance, to accurate predict the transient thermal conductivity, information on grain size, fission gas density, porosity and fracture is required. Further development is on going to improve all these models. These models are developed primarily for UO{sub 2} fuels used in LWRs. They can be extended to other fuels with possibly new development. All models in MARMOT will be rigorously verified and validated before being applied for fuel performance modeling. A version of MARMOT was recently released in September 2015 with an assessment report. (authors)

36 MATERIALS SCIENCE↗

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

Convective and morphological instability in vapor crystal growth

Theoretical and experimental work on the fluid dynamics of physical vapor transport is reported. It is shown that diffusion in viscous interaction with container walls leads to concentration gradients normal to the main transport direction. Consequently any convection threshold is removed. This coupling with convective instabilities makes the theoretical interpretation of earthbound (anisotropic) morphological stability studies on vapor-solid interfaces intractable. Hence low-gravity experiments are suggested that will allow for the establishment of morphological stability criteria under diffusion controlled conditions.

Rosenberger, F.↗

Nonlinear Landau damping in the ionosphere

A model which explains the nonresonant waves which produce the diffuse resonance observed near 3/2 f(H) by the Alouette and Isis topside sounders, where f(H) is the ambient electron cyclotron frequency, is presented. These waves are the result of plasma wave instabilities driven by anisotropic electron velocity distributions initiated by the high-power short-duration sounder pulse. Calculations of the nonlinear wave-particle coupling coefficients show that the diffuse resonance wave can be maintained by nonlinear Landau damping of the sounder-stimulated 2f(H) wave which is observed with a time duration longer than that of the diffuse resonance wave. The time duration of the diffuse resonance is determined by the transit time of the instability-generated and nonlinearly maintained diffuse resonance wave from the remote short-lived hot region back to the antenna. The model is consistent with the Alouette/Isis observations and it demonstrates the existence of nonlinear wave-particle interactions in the ionosphere.

Kiwamoto, Y.↗

A diffusion-limited reaction model for self-propagating Al/Pt multilayers with quench limits

A diffusion-limited reaction model was calibrated for Al/Pt multilayers ignited on oxidized silicon, sapphire, and tungsten substrates, as well as for some Al/Pt multilayers ignited as free-standing foils. The model was implemented in a finite element analysis code and used to match experimental burn front velocity data collected from several years of testing at Sandia National Laboratories. Moreover, both the simulations and experiments reveal well-defined quench limits in the total Al + Pt layer (i.e., bilayer) thickness. At these limits, the heat generated from atomic diffusion is insufficient to support a self-propagating wave front on top of the substrates. Quench limits for reactive multilayers are seldom reported and are found to depend on the thermal properties of the individual layers. Here, the diffusion-limited reaction model is generalized to allow for temperature- and composition-dependent material properties, phase change, and anisotropic thermal conductivity. Utilizing this increase in model fidelity, excellent overall agreement is shown between the simulations and experimental results with a single calibrated parameter set. However, the burn front velocities of Al/Pt multilayers ignited on tungsten substrates are over-predicted. Finally, possible sources of error are discussed and a higher activation energy (from 41.9 kJ/mol.at. to 47.5 kJ/mol.at.) is shown to bring the simulations into agreement with the velocity data observed on tungsten substrates. Finally, this higher activation energy suggests an inhibited diffusion mechanism present at lower heating rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A NUMERICAL TREATMENT OF ANISOTROPIC RADIATION FIELDS COUPLED WITH RELATIVISTIC RESISTIVE MAGNETOFLUIDS

We develop a numerical scheme for solving fully special relativistic, resistive radiation magnetohydrodynamics. Our code guarantees conservation of total mass, momentum, and energy. The radiation energy density and the radiation flux are consistently updated using the M-1 closure method, which can resolve an anisotropic radiation field, in contrast to the Eddington approximation, as well as the flux-limited diffusion approximation. For the resistive part, we adopt a simple form of Ohm's law. The advection terms are explicitly solved with an approximate Riemann solver, mainly the Harten-Lax-van Leer scheme; the HLLC and HLLD schemes are also solved for some tests. The source terms, which describe the gas-radiation interaction and the magnetic energy dissipation, are implicitly integrated, relaxing the Courant-Friedrichs-Lewy condition even in an optically thick regime or a large magnetic Reynolds number regime. Although we need to invert 4 Multiplication-Sign 4 matrices (for the gas-radiation interaction) and 3 Multiplication-Sign 3 matrices (for the magnetic energy dissipation) at each grid point for implicit integration, they are obtained analytically without preventing massive parallel computing. We show that our code gives reasonable outcomes in numerical tests for ideal magnetohydrodynamics, propagating radiation, and radiation hydrodynamics. We also applied our resistive code to the relativistic Petschek-type magnetic reconnection, revealing the reduction of the reconnection rate via radiation drag.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Glancing Angle Deposition in Gas Sensing: Bridging Morphological Innovations and Sensor Performances

Glancing Angle Deposition (GLAD) has emerged as a versatile and powerful nanofabrication technique for developing next-generation gas sensors by enabling precise control over nanostructure geometry, porosity, and material composition. Through dynamic substrate tilting and rotation, GLAD facilitates the fabrication of highly porous, anisotropic nanostructures, such as aligned, tilted, zigzag, helical, and multilayered nanorods, with tunable surface area and diffusion pathways optimized for gas detection. This review provides a comprehensive synthesis of recent advances in GLAD-based gas sensor design, focusing on how structural engineering and material integration converge to enhance sensor performance. Key materials strategies include the construction of heterojunctions and core–shell architectures, controlled doping, and nanoparticle decoration using noble metals or metal oxides to amplify charge transfer, catalytic activity, and redox responsiveness. GLAD-fabricated nanostructures have been effectively deployed across multiple gas sensing modalities, including resistive, capacitive, piezoelectric, and optical platforms, where their high aspect ratios, tailored porosity, and defect-rich surfaces facilitate enhanced gas adsorption kinetics and efficient signal transduction. These devices exhibit high sensitivity and selectivity toward a range of analytes, including NO2, CO, H2S, and volatile organic compounds (VOCs), with detection limits often reaching the parts-per-billion level. Emerging innovations, such as photo-assisted sensing and integration with artificial intelligence for data analysis and pattern recognition, further extend the capabilities of GLAD-based systems for multifunctional, real-time, and adaptive sensing. Finally, current challenges and future research directions are discussed, emphasizing the promise of GLAD as a scalable platform for next-generation gas sensing technologies.

Chemistry↗

A finite difference method of solving anisotropic scattering problems

A new method of solving radiative transfer problems is described including a comparison of its speed with that of the doubling method, and a discussion of its accuracy and suitability for computations involving variable optical properties. The method uses a discretization in angle to produce a coupled set of first-order differential equations which are integrated between discrete depth points to produce a set of recursion relations for symmetric and anti-symmetric angular sums of the radiation field at alternate depth points. The formulation given here includes depth-dependent anisotropic scattering, absorption, and internal sources, and allows arbitrary combinations of specular and non-Lambertian diffuse reflection at either or both boundaries. Numerical tests of the method show that it can return accurate emergent intensities even for large optical depths. The method is also shown to conserve flux to machine accuracy in conservative atmospheres

Barkstrom, B. R.↗

High-temperature properties of ceramic fibers and insulations for thermal protection of atmospheric entry and hypersonic cruise vehicles

Multilayer insulations (MIs) which will operate in the 500 to 1000 C temperature range are being considered for possible applications on aerospace vehicles subject to convective and radiative heating during atmospheric entry. The insulations described consist of ceramic fibers, insulations, and metal foils quilted together with ceramic thread. As these types of insulations have highly anisotropic properties, the total heat transfer characteristics must be determined. Data are presented on the thermal diffusivity and thermal conductivity of four types of MIs and are compared to the baseline Advanced Flexible Reusable Surface Insulation currently used on the Space Shuttle Orbiter. In addition, the high temperature properties of the fibers used in these MIs are discussed. The fibers investigated included silica and three types of aluminoborosilicate (ABS). Static tension tests were performed at temperatures up to 1200 C and the ultimate strain, tensile strength, and tensile modulus of single fibers were determined.

Kourtides, Demetrius A.↗