Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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 865 records · Page 48

The Inner Workings of Magnetic Reconnection: Diffusion Region in the Balance

The question of whether the micro scale controls the macroscale or vice-versa remains one of the most challenging problems in plasmas. A particular topic of interest within this context is collisionless magnetic reconnection, where both points of views are espoused by different groups of researchers. This presentation will focus on this topic. We will begin by analyzing the properties of electron diffusion region dynamics both for guide field and anti-parallel reconnection, and how they can be scaled to different inflow conditions. As a next step, we will study typical temporal variations of the microscopic dynamics with the objective of understanding the potential for secular changes to the macroscopic system. The research will be based on a combination of analytical theory and numerical modeling.

Hesse, Michael↗

CMB spectral distortions as solutions to the Boltzmann equations

We propose to re-interpret the cosmic microwave background spectral distortions as solutions to the Boltzmann equation. This approach makes it possible to solve the second order Boltzmann equation explicitly, with the spectral y distortion and the momentum independent second order temperature perturbation, while generation of μ distortion cannot be explained even at second order in this framework. We also extend our method to higher order Boltzmann equations systematically and find new type spectral distortions, assuming that the collision term is linear in the photon distribution functions, namely, in the Thomson scattering limit. As an example, we concretely construct solutions to the cubic order Boltzmann equation and show that the equations are closed with additional three parameters composed of a cubic order temperature perturbation and two cubic order spectral distortions. The linear Sunyaev-Zel'dovich effect whose momentum dependence is different from the usual y distortion is also discussed in the presence of the next leading order Kompaneets terms, and we show that higher order spectral distortions are also generated as a result of the diffusion process in a framework of higher order Boltzmann equations. The method may be applicable to a wider class of problems and has potential to give a general prescription to non-equilibrium physics.

79 ASTRONOMY AND ASTROPHYSICS↗

Heating of the solar chromosphere and corona. I - Generalized inhomogeneous wave equation for magnetoacoustic motions

The generalized inhomogeneous wave equation that governs magnetoacoustic, vortical, and thermal motions in compressible fluids and that is applicable to the problem of heating of the solar chromosphere and corona is obtained. The effects of kinematic and bulk viscosity, heat conduction, Joule dissipation, and magnetic diffusivity are included. Under the usual assumptions, the generalized wave equation reduces to the well-known equations of Lighthill, Kulsrud, Phillips, and others. The major problems encountered in applying Lighthill's (1952) mechanism to sound generation in turbulent media are reviewed for both the subsonic and supersonic cases.

Anand, S. P. S.↗

New Results from Fermi-LAT and Their Implications for the Nature of Dark Matter and the Origin of Cosmic Rays

The measured spectrum is compatible with a power law within our current systematic errors. The spectral index (-3.04) is harder than expected from previous experiments and simple theoretical considerations. "Pre-Fermi" diffusive model requires a harder electron injection spectrum (by 0.12) to fit the Fermi data, but inconsistent with positron excess reported by Pamela if it extends to higher energy. Additional component of electron flux from local source(s) may solve the problem; its origin, astrophysical or exotic, is still unclear. Valuable contribution to the calculation of IC component of diffuse gamma radiation.

Moiseev, Alexander↗

The cosmic-ray shock structure problem for relativistic shocks

The time asymptotic behaviour of a relativistic (parallel) shock wave significantly modified by the diffusive acceleration of cosmic-rays is investigated by means of relativistic hydrodynamical equations for both the cosmic-rays and thermal gas. The form of the shock structure equation and the dispersion relation for both long and short wavelength waves in the system are obtained. The dependence of the shock acceleration efficiency on the upstream fluid spped, long wavelength Mach number and the ratio N = P sub co/cP sub co+P sub go)(Psub co and P sub go are the upstream cosmic-ray and thermal gas pressures respectively) are studied.

Webb, G. M.↗

Plasma transport in an Eulerian AMR code

A plasma transport model has been implemented in an Eulerian AMR radiation-hydrodynamics code, xRage, which includes plasma viscosity in the momentum tensor, viscous dissipation in the energy equations, and binary species mixing with consistent species mass and energy fluxes driven by concentration gradients, ion and electron baro-diffusion terms and temperature gradient forces. The physics basis, computational issues, numeric options, and results from several test problems are discussed. The transport coefficients are found to be relatively insensitive to the kinetic correction factors when the concentrations are expressed with the molar fractions and the ion mass differences are large. The contributions to flow dynamics from plasma viscosity and mass diffusion were found to increase significantly as scale lengths decrease in an inertial confinement fusion relevant Kelvin-Helmholtz instability mix layer. The mixing scale lengths in the test case are on the order of 100 μm and smaller for viscous effects to appear and 10 μm or less for significant ion species diffusion, evident over durations on the order of nanoseconds. The temperature gradient driven mass flux is seen to deplete a high Z tracer ion at the ion shock front. The plasma transport model provides the generation of the atomic mix per unit of interfacial area between two species with no free parameters. The evolution of the total atomic mix then depends also on an accurate resolution or estimate of the interfacial area between the species mixing by plasma transport. High resolution simulations or a more Lagrangian-like treatment of species interfaces may be required to distinguish plasma transport and numerical diffusion in an Eulerian computation of complex and dynamically evolving mix regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An extended laser flash technique for thermal diffusivity measurement of high-temperature materials

Knowledge of thermal diffusivity data for high-temperature materials (solids and liquids) is very important in analyzing a number of processes, among them solidification, crystal growth, and welding. However, reliable thermal diffusivity versus temperature data, particularly those for high-temperature liquids, are still far from complete. The main measurement difficulties are due to the presence of convection and the requirement for a container. Fortunately, the availability of levitation techniques has made it possible to solve the containment problem. Based on the feasibility of the levitation technology, a new laser flash technique which is applicable to both levitated liquid and solid samples is being developed. At this point, the analysis for solid samples is near completion and highlights of the technique are presented here. The levitated solid sample which is assumed to be a sphere is subjected to a very short burst of high power radiant energy. The temperature of the irradiated surface area is elevated and a transient heat transfer process takes place within the sample. This containerless process is a two-dimensional unsteady heat conduction problem. Due to the nonlinearity of the radiative plus convective boundary condition, an analytic solution cannot be obtained. Two options are available at this point. Firstly, the radiation boundary condition can be linearized, which then accommodates a closed-form analytic solution. Comparison of the analytic curves for the temperature rise at different points to the experimentally-measured values will then provide the thermal diffusivity values. Secondly, one may set up an inverse conduction problem whereby experimentally obtained surface temperature history is used as the boundary conditions. The thermal diffusivity can then be elevated by minimizing the difference between the real heat flux boundary condition (radiation plus convection) and the measurements. Status of an experimental study directed at measuring the thermal diffusivity of high-temperature solid samples of pure Nickel and Inconel 718 superalloys are presented. Preliminary measurements showing surface temperature histories are discussed.

Shen, F.↗

Greenhouse Gas Source Attribution: Measurements Modeling and Uncertainty Quantification

In this project we have developed atmospheric measurement capabilities and a suite of atmospheric modeling and analysis tools that are well suited for verifying emissions of green- house gases (GHGs) on an urban-through-regional scale. We have for the first time applied the Community Multiscale Air Quality (CMAQ) model to simulate atmospheric CO 2 . This will allow for the examination of regional-scale transport and distribution of CO 2 along with air pollutants traditionally studied using CMAQ at relatively high spatial and temporal resolution with the goal of leveraging emissions verification efforts for both air quality and climate. We have developed a bias-enhanced Bayesian inference approach that can remedy the well-known problem of transport model errors in atmospheric CO 2 inversions. We have tested the approach using data and model outputs from the TransCom3 global CO 2 inversion comparison project. We have also performed two prototyping studies on inversion approaches in the generalized convection-diffusion context. One of these studies employed Polynomial Chaos Expansion to accelerate the evaluation of a regional transport model and enable efficient Markov Chain Monte Carlo sampling of the posterior for Bayesian inference. The other approach uses de- terministic inversion of a convection-diffusion-reaction system in the presence of uncertainty. These approaches should, in principle, be applicable to realistic atmospheric problems with moderate adaptation. We outline a regional greenhouse gas source inference system that integrates (1) two ap- proaches of atmospheric dispersion simulation and (2) a class of Bayesian inference and un- certainty quantification algorithms. We use two different and complementary approaches to simulate atmospheric dispersion. Specifically, we use a Eulerian chemical transport model CMAQ and a Lagrangian Particle Dispersion Model - FLEXPART-WRF. These two models share the same WRF assimilated meteorology fields, making it possible to perform a hybrid simulation, in which the Eulerian model (CMAQ) can be used to compute the initial condi- tion needed by the Lagrangian model, while the source-receptor relationships for a large state vector can be efficiently computed using the Lagrangian model in its backward mode. In ad- dition, CMAQ has a complete treatment of atmospheric chemistry of a suite of traditional air pollutants, many of which could help attribute GHGs from different sources. The inference of emissions sources using atmospheric observations is cast as a Bayesian model calibration problem, which is solved using a variety of Bayesian techniques, such as the bias-enhanced Bayesian inference algorithm, which accounts for the intrinsic model deficiency, Polynomial Chaos Expansion to accelerate model evaluation and Markov Chain Monte Carlo sampling, and Karhunen-Lo %60 eve (KL) Expansion to reduce the dimensionality of the state space. We have established an atmospheric measurement site in Livermore, CA and are collect- ing continuous measurements of CO 2 , CH 4 and other species that are typically co-emitted with these GHGs. Measurements of co-emitted species can assist in attributing the GHGs to different emissions sectors. Automatic calibrations using traceable standards are performed routinely for the gas-phase measurements. We are also collecting standard meteorological data at the Livermore site as well as planetary boundary height measurements using a ceilometer. The location of the measurement site is well suited to sample air transported between the San Francisco Bay area and the California Central Valley.

54 ENVIRONMENTAL SCIENCES↗

Finite-element solution for a combined radiative-conductive analysis with mixed diffuse-specular surface characteristics

The finite-element method has been applied to solve a combined radiative-conductive heat transfer problem for a large space telescope similar to those used in orbiting satellites. The derivation of the underlying matrices and associated solution algorithm for a 2-D triangular element is presented in detail. The resulting expressions for this triangular element typify such an analysis, which yields constitutive matrices when the heat equation is cast in the matrix form. The relevant matrices include those pertaining to thermal conductance, internal heat generation, radiative exchanges, and all possible external thermal loadings. Emphasis is placed on the treatment of non-linear radiative interchange between surfaces in an enclosure having mixed diffuse-specular surface characteristics. Essential differences in governing equations describing these distinctive surface characteristics are identified. Concluding remarks are drawn from an example simulating a Cassegrainian space telescope.

Lee, H.-P.↗

Decelerating Flows in TeV Blazars: A Resolution to the BL Lac -- FR I Unification Problem

TeV emission from BL Lacertae (BL) objects i s commonly modeled as Synchrotron-Self Compton (SSC) radiation from relativistically moving homogeneous plasma blobs. In the context of these models, the blob Lorentz factors needed to reproduce the corrected for absorption by the diffuse IR background (DIRB) TeV emission are large ($\delta \gtrsim 50$). The main reason for this is that stronger beaming eases the problem of the lack of $\sim$ IR-UV synchrotron seed photons needed to produce the de-absorbed $\sim $ few TeV peak of the spectral energy distribution (SED). However, such high Doppler factors are in strong disagreement with the unified scheme, according to which BLs are FR I radio galaxies with their jets closely aligned to the line of sight. Here, motivated by the detection of sub-luminal velocities in the sub-pc scale jets of the best studied TeV blazars, MKN 421 and MKN 501, we examine the possibility that the relativistic flows in the TeV BLs decelerate. In this case, the problem of the missing seed photons is solved because of Upstream Compton (UC) scattering, a process in which the upstream energetic electrons from the fast base of the flow 'see' the synchrotron seed photons produced in the slow part of the flow relativistically beamed. Modest Lorentz factors ($\Gamma \sim 15$), decelerating down to values compatible with the recent radio interferometric observations, reproduce the $\sim $ few TeV peak energy of these sources. Furthermore, such decelerating flows are shown to be in agreement with the BL - FR I unification, naturally reproducing the observed BL/FR I broadband luminosity ratios.

Kazanas, Demosthenes↗

Suppressing 1 H Spin Diffusion in Fast MAS Proton Detected Heteronuclear Correlation Solid-State NMR Experiments

Fast magic angle spinning (MAS) and indirect detection by high gyromagnetic ratio (γ) nuclei such as proton or fluorine are increasingly utilized to obtain 2D heteronuclear correlation (HETCOR) solid-state NMR spectra of spin-1/2 nuclei by using cross polarization (CP) for coherence transfer. However, one major drawback of CP HETCOR pulse sequences is that 1 H spin diffusion during the back X→ 1 H CP transfer step results in relayed correlations. This problem is particularly pronounced for the indirect detection of very low-γ nuclei such as 89 Y, 103 Rh, 109 Ag and 183 W where long contact times on the order of 10 to 30 ms are necessary for optimal CP transfer. In this work, we propose two methods that eliminate relayed correlations and allow more reliable distance information to be obtained from 2D HETCOR NMR spectra. The first method uses Lee-Goldburg (LG) CP during the X→ 1 H back-transfer step to suppress 1 H spin diffusion. We determine LG conditions compatible with fast MAS frequencies (ν rot ) of 40- 95 kHz and show that 1 H spin diffusion can be efficiently suppressed at low effective radiofrequency (RF) fields (ν 1,eff << 0.5ν rot ) and also at high effective RF fields (ν 1,eff >> 2ν rot ). We describe modified Hartmann-Hahn LG-CP match conditions compatible with fast MAS and suitable for indirect detection of moderate-γ nuclei such as 13 C, and low-γ nuclei such as 89 Y. The second method uses D-RINEPT (dipolar refocused insensitive nuclei enhanced by polarization transfer) during the X→ 1 H back-transfer step of the HETCOR pulse sequence. The effectiveness of these methods for acquiring HETCOR spectra with reduced relayed signal intensities is demonstrated with 1 H{ 13 C} HETCOR NMR experiments on L-histidine•HCl•H 2 O and 1 H{ 89 Y} HETCOR NMR experiments on an organometallic yttrium complex.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NOR-USA Scientific Traverse of East Antarctica: Science and Logistics on a Three-Month Expedition Across Antarctica's Farthest Frontier

Dr. Albert's current research is centered on transfer processes in porous media, including air-snow exchange in the Polar Regions and in soils in temperate areas. Her research includes field measurements, laboratory experiments, and theoretical modeling. Mary conducts field and laboratory measurements of the physical properties of natural terrain surfaces, including permeability, microstructure, and thermal conductivity. Mary uses the measurements to examine the processes of diffusion and advection of heat, mass, and chemical transport through snow and other porous media. She has developed numerical models for investigation of a variety of problems, from interstitial transport to freezing of flowing liquids. These models include a two-dimensional finite element code for air flow with heat, water vapor, and chemical transport in porous media, several multidimensional codes for diffusive transfer, as well as a computational fluid dynamics code for analysis of turbulent water flow in moving-boundary phase change problems.

Albert, Mary R.↗

Diffusion, atomic transport, and ordering in Al-Zr alloys: FCC and liquid phases

Additive manufacturing of materials with controlled microstructure demands knowledge of atomic scale properties near the solid-liquid transition state. Many of these properties are not affordable by experimental techniques and computer modeling is the possible solution to the problem. In this paper, we present the results of an extended atomistic study of intrinsic atomic transport due to vacancy diffusion in FCC and L12 solid phases and diffusion in the liquid phase of Al-Zr alloys. A deceleration of the overall self-diffusion was observed when Zr was added to Al. The effect was stronger in the solid and weaker in the liquid. Additionally, the effect was strongly temperature dependent in the solid phases, but not in the liquid. Atomic transport was chemically biased: transport of Zr atoms was significantly slower than that of Al atoms, and this bias effect was stronger in the solid phases. The overall diffusion and chemical ordering processes in the liquid state were five to six orders in magnitude faster than in the solid. Chemical short-range order parameters in the liquid saturated at values close to those in the ordered L12 structure of Al3Zr. Chemical and structural ordering in the solid phases was negligible over the modeled microsecond time scale. Here, the results are discussed in view of optimizing additive manufacturing parameters for the controlled formation of metastable L1 2 precipitates.

36 MATERIALS SCIENCE↗

Density Large Deviations for Multidimensional Stochastic Hyperbolic Conservation Laws

We investigate the density large deviation function for a multidimensional conservation law in the vanishing viscosity limit, when the probability concentrates on weak solutions of a hyperbolic conservation law. When the mobility and diffusivity matrices are proportional, i.e. an Einstein-like relation is satisfied, the problem has been solved in Bellettini and Mariani (Bull Greek Math Soc 57:31–45, 2010). When this proportionality does not hold, we compute explicitly the large deviation function for a step-like density profile, and we show that the associated optimal current has a non trivial structure. We also derive a lower bound for the large deviation function, valid for a more general weak solution, and leave the general large deviation function upper bound as a conjecture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Homology modeling and docking of AahII-Nanobody complexes reveal the epitope binding site on AahII scorpion toxin

Highlights: • AahII specific nanobodies share common CDRs sequence patterns. • Residues forming sequences' pattern interact specifically with AahII. • CDR1 Tyr27, Tyr29 and Ser30 polar residues are involved in AahII neutralization. • Interacting AahII Arg62/NbAahII10 CDR3 Asp112 charged residues are crucial. • AahII-NbAahII10 trapping involves a distinct toxin 3D orientation. Scorpion envenoming and its treatment is a public health problem in many parts of the world due to highly toxic venom polypeptides diffusing rapidly within the body of severely envenomed victims. Recently, 38 AahII-specific Nanobody sequences (Nbs) were retrieved from which the performance of NbAahII10 nanobody candidate, to neutralize the most poisonous venom compound namely AahII acting on sodium channels, was established. Herein, structural computational approach is conducted to elucidate the Nb-AahII interactions that support the biological characteristics, using Nb multiple sequence alignment (MSA) followed by modeling and molecular docking investigations (RosettaAntibody, ZDOCK software tools). Sequence and structural analysis showed two dissimilar residues of NbAahII10 CDR1 (Tyr27 and Tyr29) and an inserted polar residue Ser30 that appear to play an important role. Indeed, CDR3 region of NbAahII10 is characterized by a specific Met104 and two negatively charged residues Asp115 and Asp117. Complex dockings reveal that NbAahII17 and NbAahII38 share one common binding site on the surface of the AahII toxin divergent from the NbAahII10 one's. At least, a couple of NbAahII10 – AahII residue interactions (Gln38 – Asn44 and Arg62, His64, respectively) are mainly involved in the toxic AahII binding site. Altogether, this study gives valuable insights in the design and development of next generation of antivenom.

60 APPLIED LIFE SCIENCES↗

Metastability of stratified magnetohydrostatic equilibria and their relaxation

Motivated by explosive releases of energy in fusion, space and astrophysical plasmas, we consider the nonlinear stability of stratified magnetohydrodynamic equilibria against two-dimensional interchanges of straight magnetic-flux tubes. We demonstrate that, even within this restricted class of dynamics, the linear stability of an equilibrium does not guarantee its nonlinear stability: equilibria can be metastable. We show that the minimum-energy state accessible to a metastable equilibrium under non-diffusive two-dimensional dynamics can be found by solving a combinatorial optimisation problem. These minimum-energy states are, to good approximation, the final states reached by our simulations of destabilised metastable equilibria for which turbulent mixing is suppressed by viscosity. To predict the result of fully turbulent relaxation, we construct a statistical mechanical theory based on the maximisation of Boltzmann's mixing entropy. This theory is analogous to the Lynden-Bell statistical mechanics of collisionless stellar systems and plasma, and to the Robert–Sommeria–Miller theory of two-dimensional vortex turbulence. Our theory reproduces well the results of our numerical simulations for sufficiently large perturbations to the metastable equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the crust of neutron stars with neural-network quantum states

An accurate description of low-density nuclear matter is crucial for explaining the physics of neutron star crusts. In the density range between approximately 0.01 fm −3 and 0.1 fm −3 , matter transitions from neutron-rich nuclei to various higher-density pasta shapes, before ultimately reaching a uniform liquid. In this work, we introduce a variational Monte Carlo method based on a neural Pfaffian-Jastrow quantum state, which allows us to model the transition from the liquid phase to neutron-rich nuclei microscopically. At low densities, nuclear clusters dynamically emerge from the microscopic interactions among protons and neutrons, which we model based on pionless effective field theory. Our variational Monte Carlo approach represents a significant improvement over the state-of-the-art auxiliary-field diffusion Monte Carlo method, which is severely hindered by the fermion-sign problem in this low-density regime and cannot capture the onset of clusters. In addition to computing the energy per particle of symmetric nuclear matter and pure neutron matter, we analyze an intermediate isospin-asymmetry configuration to elucidate the formation of nuclear clusters. We also provide evidence that the presence of such nuclear clusters influences the amount of protons in the crust compared to protons in beta-equilibrated, neutrino-transparent matter.

Nuclear astrophysics↗