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 163 records · Page 9

Formation of a conical distribution and intense ion heating in the presence of hydrogen cyclotron waves

In the considered investigation, it is assumed that the field aligned currents are responsible for producing electrostatic harmonic cyclotron waves (EHC). Using a one-dimensional simulation model in which the electron velocity distribution is maintained by a constant injection of the initial distribution, it is shown that, in contrast to earlier initial value simulations, EHC waves grow to a large amplitude, resulting in the formation of an anisotropic ion velocity distribution. Both the heating rate and the anisotropy are in reasonable agreement with the quasi-linear theory, taking into account the cyclotron resonance. The results show that the saturation is due to the combined effects of wave induced diffusion in an electron velocity space and the heating of ions perpendicularly. Both these effects reduce the growth rate.

Okuda, H.↗

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING↗

The anisotropic radiative transfer problem in optically thick, strongly magnetized plasma - A comparison of results

The accuracy of four different methods for solving the coherent radiative transfer problem in a strongly magnetized plasma is compared. Specific attention is given to the case of a semi-infinite homogeneous plasma with a magnetic field perpendicular to the surface. The four approaches evaluated were: diffusion approximation schemes; the generalized H-function method of Silant'ev (1982) the Feautrier method of Nagel (1981) and the integral equation approach of Meszares and Bonazzola (1981). It is found that the coupled diffusion approximation schemes of Nagel (1981) and Kaminker et al. (1982) provided satisfactory results and were more efficient in many cases than direct numerical methods.

Pavlov, G. G.↗

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. 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.

Methane-d↗

Mirror instability and origin of morningside auroral structure

Auroral optical imagery shows marked differences between auroral features of the evening and morning sectors: the separation between diffuse and discrete auroras in the evening sector is not distinct in the morning sector, which is dominated by auroral patches and multiple banded structures aligned along some direction. Plasma distribution function signatures also show marked differences: downward electron beams and inverted-V signatures prefer the evening sector, while the electron spectra on the morning sector are similar to the diffuse aurora. A theory of morningside auroras consistent with these features was constructed. The theory is based on modulation of the growth rates of electron cyclotron waves by the mirror instability, which is in turn driven by inward-convected ions that have become anisotropic. This modulation produces alternating bands of enhanced and reduced electron precipitation which approximate the observed multiple auroral bands and patches of the morning sector.

Chiu, Y. T.↗

Spectral ordering techniques for incomplete LU preconditoners for CG methods

The effectiveness of an incomplete LU (ILU) factorization as a preconditioner for the conjugate gradient method can be highly dependent on the ordering of the matrix rows during its creation. Detailed justification for two heuristics commonly used in matrix ordering for anisotropic problems is given. The bandwidth reduction and weak connection following heuristics are implemented through an ordering method based on eigenvector computations. This spectral ordering is shown to be a good representation of the heuristics. Analysis and test cases in two and three dimensional diffusion problems demonstrate when ordering is important, and when an ILU decomposition will be ordering insensitive. The applicability of the heuristics is thus evaluated and placed on a more rigorous footing.

Clift, Simon S.↗

Mirror instability and the origin of morningside auroral structure

Auroral optical imagery shows marked differences between auroral features of the evening and morning sectors: The separation between diffuse and discrete auroras in the evening sector is not distinct in the morning sector, which is dominated by auroral patches and multiple banded structures aligned along some direction. Plasma distribution function signatures also show marked differences: downward electron beams and inverted-V signatures prefer the evening sector, while the electron spectra on the morning sector are similar to the diffuse aurora. We have constructed a theory of morningside auroras consistent with these features. The theory is based on modulation of the growth rates of electron cyclotron waves by the mirror instability, which is in turn driven by inward-convected ions that have become anisotropic. This modulation produces alternating bands of enhanced and reduced electron precipitation which approximate the observed multiple auroral bands and patches of the morning sector.

Chiu, Y. T.↗

Turbulent solutions of the equations of fluid motion

Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. Initial three-dimensional cosine velocity fluctuations and periodic boundary conditions are used in most of the work considered. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations, the initially nonrandom flow develops into an apparently random turbulence. Thus randomness or turbulence can arise as a consequence of the structure of the Navier-Stokes equations. The cases considered include turbulence which is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A mean shear is present in some cases. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components, and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.

Deissler, R. G.↗

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↗

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.↗

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.↗

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↗

Tracking Dendritic Growth in Hydrogen-Based Hematite Reduction via Computer Vision

The reduction of hematite to metallic iron using hydrogen (H2) as a reducing agent presents a promising pathway for decarbonizing steel production. In this study, we employ a combination of in situ confocal scanning laser microscopy (CSLM) and advanced computer vision techniques to quantitatively analyze dendritic growth of ferrite during H2-based reduction of iron oxide at high temperatures. A workflow integrating Watershed Image Segmentation (WIS) and Lucas-Kanade Optical Flow (LKOF) is developed to extract both global and local kinetic information from time-resolved micrograph sequences. H2 reduction experiments conducted at 1400 degrees C and 1500 degrees C demonstrate a clear correlation between temperature and reduction rate, as evidenced by accuracy of fitted Johnson-Mehl-Avrami-Kolmogorov (JMAK) parameters. Optical flow analysis further elucidates the anisotropic and branched nature of dendritic growth, providing spatially resolved velocity fields that correlate well with global transformation kinetics. The proposed methodology demonstrates strong agreement with experimental measurements and literature values, offering a robust framework for automated image-based analysis to study kinetics through microstructural evolution in the reduction of iron ore, and likely other reaction-diffusion phenomena.

08 HYDROGEN↗

Atomic to molecular hydrogen transition in interstellar clouds

The photodestruction of molecular hydrogen is reexamined in the context of recent observational data. Improvements in the theory of radiation transfer are made by using a more accurate self-shielding integral, treating anisotropic scattering by dust grains, and including photodissociation from excited rotational states. The total column density of hydrogen required to obtain 10% conversion decreases with a parameter, epsilon, the ratio of the formation rate to the destruction rate, as epsilon to the -m power, with m in the range 1.4-1.6. Satisfactory agreement with the data for thicker diffuse clouds is obtained for epsilon of about 0.00006.

Federman, S. R.↗

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗