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 451 records · Page 25

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is a suite of tools for the analysis of porous materials and generation of material microstructures. From microstructural data, often obtained through X-ray microtomography, PuMA can determine a number of effective material properties and perform material response simulations. Version 2.2 includes capabilities for computing volume fractions, porosity, specific surface area, effective thermal and electrical conductivities, and continuum and rarefied diffusive tortuosity. PuMA can also simulate competitive diffusion/reaction processes at the micro-scale, such as surface oxidation. In this poster, recent advancements to the PuMA software are detailed, including the full refactoring of PuMA into v3.0, a new module to compute heat conduction in anisotropic materials, a particle method for simulating molecular beam experiments, a new finite-volume Laplace solver, complex fibrous material generation, woven material generation, and a coupling of PuMA with the DAKOTA software for advanced statistics.

PuMA↗

The First Widespread Solar Energetic Particle Event Observed By Solar Orbiter on 2020 November 29

Context. On 2020 November 29, the first widespread solar energetic particle (SEP) event of solar cycle 25 was observed at four widely separated locations in the inner (. 1AU) heliosphere. Relativistic electrons as well as protons with energies > 50MeV were observed by Solar Orbiter (SolO), Parker Solar Probe (PSP), the Solar Terrestrial Relations Observatory (STEREO)-A and multiple near-Earth spacecraft. The SEP event was associated with an M4.4 class X-ray flare and accompanied by a coronal mass ejection (CME) and an extreme ultraviolet (EUV) wave as well as a type II radio burst and multiple type III radio bursts. Aims. We present multi-spacecraft particle observations and place them in context with source observations from remote sensing instruments and discuss how such observations may further our understanding of particle acceleration and transport in this widespread event. Methods. Velocity dispersion analysis (VDA) and time shift analysis (TSA) were used to infer the particle release times at the Sun. Solar wind plasma and magnetic field measurements were examined to identify structures that influence the properties of the energetic particles such as their intensity. Pitch angle distributions and first-order anisotropies were analyzed in order to characterize the particle propagation in the interplanetary medium. Results. We find that during the 2020 November 29 SEP event, particles spread over more than 230° in longitude close to 1AU. The particle onset delays observed at the different spacecraft are larger as the flare–footpoint angle increases and are consistent with those from previous STEREO observations. Comparing the timing when the EUV wave intersects the estimated magnetic footpoints of each spacecraft with particle release times from TSA and VDA, we conclude that a simple scenario where the particle release is only determined by the EUV wave propagation is unlikely for this event. Observations of anisotropic particle distributions at SolO, Wind, and STEREO-A do not rule out that particles are injected over a wide longitudinal range close to the Sun. However, the low values of the first-order anisotropy observed by near-Earth spacecraft suggest that diffusive propagation processes are likely involved.

A Kollhoff↗

The j and k dependencies of the line coupling and line mixing effects: Theoretical studies of the relaxation matrices of N 2 -broadened CH 3 D

Line coupling and line mixing effects in parallel and perpendicular bands of CH 3 D perturbed by N 2 have been studied. The work focuses on exhibiting the j and k dependencies of these two processes. The calculations were based on a previously reported anisotropic intermolecular potential including both the long-range multipole, induction, and dispersion forces and a short-range atom-atom model. Throughout this work, it is shown that components with L 1 = 3 of the atom-atom model are dominant contributions to the diffusion operator. As a consequence, in comparison with other molecular systems such as the CH 3 Cl-N 2 and CH 3 I-N 2 , theoretically predicted line coupling and line mixing effects exhibit completely new j and k dependencies. In general, the theoretically calculated halfwidths and intra-doublets’ off-diagonal elements of the relaxation matrix are in reasonable agreement with measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure and lattice dynamics of 1⁢𝑇−VSe 2 : Origin of the three-dimensional charge density wave

To characterize in detail the charge density wave (CDW) transition of 1⁢𝑇−VSe 2 , its electronic structure and lattice dynamics are comprehensively studied by means of x-ray diffraction, muon spectroscopy, angle resolved photoemission (ARPES), diffuse and inelastic x-ray scattering, and state-of-the-art first-principles density functional theory calculations. Resonant elastic x-ray scattering does not show any resonant enhancement at either V or Se, indicating that the CDW peak at the 𝐾 edges describes a purely structural modulation of the electronic ordering. ARPES experiments identify (i) a pseudogap at 𝑇 > 𝑇⁢ CDW , which leads to a depletion of the density of states in the ML-M'L' plane at 𝑇 < 𝑇⁢ CDW , and (ii) anomalies in the electronic dispersion reflecting a sizable impact of phonons on it. A diffuse scattering precursor, characteristic of soft phonons, is observed at room temperature (RT) and leads to the full collapse of the low-energy phonon (𝜔 1 ) with propagation vector (0.25 0 −0.3) r.l.u. Here, we show that the frequency and linewidth of this mode are anisotropic in momentum space, reflecting the momentum dependence of the electron-phonon interaction (EPI), hence demonstrating that the origin of the CDW is, to a much larger extent, due to the momentum dependent EPI with a small contribution from nesting. The pressure dependence of the 𝜔 1 soft mode remains nearly constant up to 13 GPa at RT, with only a modest softening before the transition to the high-pressure monoclinic C2/m phase. The wide set of experimental data is well captured by our state-of-the art first-principles anharmonic calculations with the inclusion of van der Waals corrections in the exchange-correlation functional. The comprehensive description of the electronic and dynamical properties of VSe 2 reported here adds important pieces of information to the understanding of the electronic modulations in the family of transition-metal dichalcogenides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO 2 in fractured rocks

Analytical solutions with infinite exponential series are available to calculate the rate of diffusive transfer between low-permeability blocks and high-permeability zones in the subsurface. Truncation of these series is often employed by neglecting the early-time regime. Here in this paper, we present unified-form approximate solutions in which the early-time and the late-time solutions are continuous at a switchover time. The early-time solutions are based on three-term polynomial functions in terms of square root of dimensionless time, with the first coefficient dependent only on the dimensionless area-to-volume ratio. The last two coefficients are either determined analytically for isotropic blocks (e.g., spheres and slabs) or obtained by fitting the exact solutions, and they solely depend on the aspect ratios for rectangular columns and parallelepipeds. For the late-time solutions, only the leading exponential term is needed for isotropic blocks, while a few additional exponential terms are needed for highly anisotropic rectangular blocks. The optimal switchover time is between 0.157 and 0.229, with highest relative approximation error less than 0.2%. The solutions are used to demonstrate the storage of dissolved CO 2 in fractured reservoirs with low-permeability matrix blocks of single and multiple shapes and sizes. These approximate solutions are building blocks for development of analytical and numerical tools for hydraulic, solute, and thermal diffusion processes in low-permeability matrix blocks.

54 ENVIRONMENTAL SCIENCES↗

Kinetics modeling of precipitation with characteristic shape during post-implantation annealing

In this study, we investigated the precipitation with characteristic shape in the microstructure during post-implantation annealing via a theoretical modeling approach. The processes of precipitates formation and evolution during phase separation were based on a nucleation and growth mechanism of atomic diffusion. Different stages of the precipitation, including the nucleation, growth and coalescence, were distinctly revealed in the numerical simulations. In addition, the influences of ion dose, temperature and crystallographic symmetry on the processes of faceted precipitation were also demonstrated. To comprehend the kinetic mechanism, the simulation results were further analyzed quantitatively by the Kolmogorov-Johnson-Mehl-Avrami (KJMA) equation. The Avrami exponents obtained from the regression curves varied from 1.47 to 0.52 for different conditions. With the increase of ion dose and temperature, the nucleation and growth of precipitations were expedited in accordance with the shortened incubation time and the raised coefficient of growth rate. A miscellaneous shape of precipitates in various crystallographic symmetry systems could be simulated through this anisotropic model. From the analyses of the kinetics, more fundamental information about the nucleation and growth mechanism of faceted precipitation during post-implantation annealing was acquired for future application.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

ACCELERATION OF LOW-ENERGY IONS AT PARALLEL SHOCKS WITH A FOCUSED TRANSPORT MODEL

We present a test particle simulation on the injection and acceleration of low-energy suprathermal particles by parallel shocks with a focused transport model. The focused transport equation contains all necessary physics of shock acceleration, but avoids the limitation of diffusive shock acceleration (DSA) that requires a small pitch angle anisotropy. This simulation verifies that the particles with speeds of a fraction of to a few times the shock speed can indeed be directly injected and accelerated into the DSA regime by parallel shocks. At higher energies starting from a few times the shock speed, the energy spectrum of accelerated particles is a power law with the same spectral index as the solution of standard DSA theory, although the particles are highly anisotropic in the upstream region. The intensity, however, is different from that predicted by DSA theory, indicating a different level of injection efficiency. It is found that the shock strength, the injection speed, and the intensity of an electric cross-shock potential (CSP) jump can affect the injection efficiency of the low-energy particles. A stronger shock has a higher injection efficiency. In addition, if the speed of injected particles is above a few times the shock speed, the produced power-law spectrum is consistent with the prediction of standard DSA theory in both its intensity and spectrum index with an injection efficiency of 1. CSP can increase the injection efficiency through direct particle reflection back upstream, but it has little effect on the energetic particle acceleration once the speed of injected particles is beyond a few times the shock speed. This test particle simulation proves that the focused transport theory is an extension of DSA theory with the capability of predicting the efficiency of particle injection.

79 ASTRONOMY AND ASTROPHYSICS↗

Persistent low-energy phonon broadening near the charge-order q vector in the bilayer cuprate Bi 2 Sr 2 CaCu 2 O 8 + δ

Here, we report a persistent low-energy phonon broadening around q B ~ 0.28 reciprocal lattice unit (r.l.u.) along the Cu-O bond direction in the high- T c cuprate Bi 2 Sr 2 CaCu 2 O 8+δ (Bi-2212). We show that such broadening exists both inside and outside the conventional charge-density wave (CDW) phase, via temperature-dependent measurements in both underdoped and heavily overdoped samples. Combining inelastic hard x-ray scattering, diffuse scattering, angle-resolved photoemission spectroscopy, and resonant soft x-ray scattering at the Cu L 3 -edge, we exclude the presence of a CDW in the heavily overdoped Bi-2212 similar to that observed in the underdoped systems. Finally, we discuss the origin of such anisotropic low-energy phonon broadening and its potential precursory role to the CDW phase in the underdoped region.

36 MATERIALS SCIENCE↗

Freeze Tape Casting of Functionally Graded Porous Ceramics

Freeze tape casting is a means of making preforms of ceramic sheets that, upon subsequent completion of fabrication processing, can have anisotropic and/or functionally graded properties that notably include aligned and graded porosity. Freeze tape casting was developed to enable optimization of the microstructures of porous ceramic components for use as solid oxide electrodes in fuel cells: Through alignment and grading of pores, one can tailor surface areas and diffusion channels for flows of gas and liquid species involved in fuel-cell reactions. Freeze tape casting offers similar benefits for fabrication of optimally porous ceramics for use as catalysts, gas sensors, and filters.

Sofie, Stephen W.↗

Impact of electrode porosity architecture on electrochemical performances of 1 mm-thick LiFePO4 binder-free Li-ion electrodes fabricated by Spark Plasma Sintering

Thick electrodes with high active material loadings have been intensively studied over the last couple of decades in pursuit of achieving high energy density systems. To optimize and enhance the electrochemical performance of such electrodes, one has to control the pore morphology by, for example, varying the pore size and shape, and the level of porosity. In the present work, the fabrication of thick binder-free LiFePO4 (LFP) electrodes with two different pore sizes (12 and 20 mu m) and porosities (21 vol% and 44 vol%) using Spark Plasma Sintering (SPS) and templating approach is reported. The well-controlled porous architecture inside the thick electrodes is realized by fine-tuning experimental parameters. The impact of porosity architecture on electrochemical performance is quantified and correlated with the 3D tortuosity values determined from both micro-computed tomography and electrochemical impedance-based experimental methods. Based on the micro-computed tomography data analysis, estimated tortuosity values along X, Y, and Z axes reveal an anisotropic effect perpendicularly to the SPS compression axis (Z-direction). This is particularly profoundly observed in the samples with larger pores (20 mu m). The correlation between morphological properties and the rate capability performance is established indicating that the capacity loss happens mainly due to the Li-ion diffusion limitations.

Cathode material↗

Particle-in-cell modeling of electron beam generated plasma

Plasmas generated using energetic electron beams are well known for their low electron temperature (T e ) and plasma potential, which makes them attractive for atomic-precision plasma processing applications such as atomic layer etch and deposition. A 2-dimensional particle-in-cell model for an electron beam-generated plasma in argon confined by a constant applied magnetic field is described here in this article. Plasma production primarily occurs in the path of the beam electrons in the center of the chamber. The resulting plasma spreads out in the chamber through non-ambipolar diffusion with a short-circuiting effect allowing unequal electron and ion fluxes to different regions of the bounding conductive chamber walls. The cross-field transport of the electrons (and thus the steady-state characteristics of the plasma) are strongly impacted by the magnetic field. T e is anisotropic in the electron beam region, but low and isotropic away from the plasma production zone. The plasma density increases and the plasma becomes more confined near the region of production when the magnetic field strengthens. The magnetic field reduces both electron physical and energy transport perpendicular to the magnetic field. T e is uniform along the magnetic field lines and slowly decreases perpendicular to it. Electrons are less energetic in the sheath regions where the sheath electric field repels and confines the low-energy electrons from the bulk plasma. Even though electron and ion densities are similar in the bulk plasma due to quasi-neutrality, electron and ion fluxes on the grounded chamber walls are unequal at most locations. Electron confinement by the magnetic field weakens with increasing pressure, and the plasma spread out farther from the electron beam region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

In situ estimation of ice crystal properties at the South Pole using LED calibration data from the IceCube Neutrino Observatory

The IceCube Neutrino Observatory instruments about 1 km 3 of deep, glacial ice at the geographic South Pole. It uses 5160 photomultipliers to detect Cherenkov light emitted by charged relativistic particles. An unexpected light propagation effect observed by the experiment is an anisotropic attenuation, which is aligned with the local flow direction of the ice. We examine birefringent light propagation through the polycrystalline ice microstructure as a possible explanation for this effect. The predictions of a first-principles model developed for this purpose, in particular curved light trajectories resulting from asymmetric diffusion, provide a qualitatively good match to the main features of the data. This in turn allows us to deduce ice crystal properties. Since the wavelength of the detected light is short compared to the crystal size, these crystal properties include not only the crystal orientation fabric, but also the average crystal size and shape, as a function of depth. By adding small empirical corrections to this first-principles model, a quantitatively accurate description of the optical properties of the IceCube glacial ice is obtained. In this paper, we present the experimental signature of ice optical anisotropy observed in IceCube light-emitting diode (LED) calibration data, the theory and parameterization of the birefringence effect, the fitting procedures of these parameterizations to experimental data, and the inferred crystal properties.

54 ENVIRONMENTAL SCIENCES↗

Modeling Lithium Diffusion in Battery Cathodes Considering Chemo-Mechanically Induced Damage

This talk will present a 3D, continuum-level damage model for simulating Lithium diffusion within generated Li_{x}Ni_{0.5}Mn_{0.3}Co_{0.2} (NMC 532) secondary cathode particles. The primary motivation of the particle-level model is to inform cathode-particle design and determine charging profiles that reduce cathode fracture. The model considers NMC 532 secondary particles containing an agglomeration of anisotropic, randomly oriented grains. The model predicts that secondary-particle fracture is primarily due to non-ideal grain interactions with slight dependence on high-rate charge demands. The model predicts that small secondary-particles with large grains develop significantly less damage than larger secondary particles with small grains. Finally, the model predicts most of the chemo-mechanical damage accumulates in the first high-rate cycles. This chemo-mechanical “damage saturation” effect indicates that initial secondary-particle fracture occurs within the first few cycles, while long-term cathode degradation is not solely chemo-mechanically induced.

47 OTHER INSTRUMENTATION↗

Mechanochemical topological defects in an active nematic

We propose a reaction-diffusion system that converts topological information of an active nematic into chemical signals. We show that a curvature-activated reaction dipole is sufficient for creating a system that dynamically senses topology by producing a concentration field possessing local extrema coinciding with ±$\frac{1}{2}$ defects. The enabling term is analogous to polarization charge density seen in dielectric materials. We demonstrate the ability of this system to identify defects in both passive and active nematics. Our results illustrate that a relatively simple feedback scheme, expressed as a system of partial differential equations, is capable of producing chemical signals in response to inherently nonlocal structures in anisotropic media. Here, we posit that such coarse-grained systems can help generate testable hypotheses for regulated processes in biological systems, such as morphogenesis, and motivate the creation of bio-inspired materials that utilize dynamic coupling between nematic structure and biochemistry.

42 ENGINEERING↗

Tidal evolution of globular clusters. II - The effects of Galactic tidal field and diffusion

The tidal evolution of globular clusters subject to various degrees of the internal diffusion process is investigated. In cases of negligible diffusion, clusters are found to be tidally truncated to the theoretical tidal radius at perigalacticon. There is no apparent orbital phase dependence of the tidal radius for clusters with eccentric orbits. In clusters with moderately short two-body relaxation time scales, diffusion processes significantly modify the structure of the outer regions in such a way that the limiting radius may be comparable to the tidal radius at apogalacticon. The Galactical tidal torque induces isotropy in the velocity dispersion of the outer regions of the cluster. For relaxed clusters, the velocity dispersion may be isotropic in the core, anisotropic in the envelope and isotropic near the limiting radius. Disk shocking is also very efficient for isotropizing the orbits of stars in the outer cluster regions. Stars with direct orbits are less stable, so that prolonged tidal interaction can lead to apparent retrograde rotation in the outer regions of the cluster.

Oh, K. S.↗

Coherent-backscatter effect - A vector formulation accounting for polarization and absorption effects and small or large scatterers

Previous theoretical work on the coherent-backscatter effect in the context of speckle time autocorrelation has gone beyond the diffusion approximation and the assumption of isotropic (point) scatterers. This paper extends the theory to include the effects of polarization and absorption, and to give the angular line shape. The results are expressions for angular variations valid for small and large scatterers and linear and circular polarizations, in lossless or lossy media. Calculations show that multiple anisotropic scattering results in the preservation of incident polarization. Application to a problem in radar astronomy is considered. It is shown that the unusual radar measurements (high reflectivity and polarization ratios) of Jupiter's icy Galilean satellites can be explained by coherent backscatter from anisotropic (forward) scatterers.

Peters, Kenneth J.↗

Alfvénic turbulence beyond the ambipolar diffusion scale

We investigate the nature of the Alfvénic turbulence cascade in two-fluid magnetohydrodynamic (MHD) simulations in order to determine if turbulence is damped once the ion and neutral species become decoupled at a critical scale called the ambipolar diffusion scale (L{sub AD}). Using mode decomposition to separate the three classical MHD modes, we study the second-order structure functions of the Alfvén mode velocity field of both neutrals and ions in the reference frame of the local magnetic field. On scales greater than L{sub AD} we confirm that two-fluid turbulence strongly resembles single-fluid MHD turbulence. Our simulations show that the behavior of two-fluid turbulence becomes more complex on scales less than L{sub AD}. We find that Alfvénic turbulence can exist past L{sub AD} when the turbulence is globally super-Alfvénic, with the ions and neutrals forming separate cascades once decoupling has taken place. When turbulence is globally sub-Alfvénic and hence strongly anisotropic, with a large separation between the parallel and perpendicular decoupling scales, turbulence is damped at L{sub AD}. We also find that the power spectrum of the kinetic energy in the damped regime is consistent with a k{sup −4} scaling (in agreement with the predictions of Lazarian et al.).

79 ASTRONOMY AND ASTROPHYSICS↗

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗