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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74

All sources NASA NTRS DOE OSTI

NASA NTRS1983

The physics of a single-event upset in integrated circuits: A review and critique of analytical models for charge collection

When an energetic particle (kinetic energy 0.5 MeV) originating from a radioactive decay or a cosmic ray transverse the active regions of semiconductor devices used in integrated circuit (IC) chips, it leaves along its track a high density electron hole plasma. The subsequent decay of this plasma by drift and diffusion leads to charge collection at the electrodes large enough in most cases to engender a false reading, hence the name single-event upset (SEU). The problem of SEU's is particularly severe within the harsh environment of Jupiter's radiation belts and constitutes therefore a matter of concern for the Galileo mission. The physics of an SEU event is analyzed in some detail. Owing to the predominance of nonlinear space charge effects and the fact that positive (holes) and negative (electrons) charges must be treated on an equal footing, analytical models for the ionized-charge collection and their corresponding currents as a function of time prove to be inadequate even in the simplest case of uniformly doped, abrupt p-n junctions in a one-dimensional geometry. The necessity for full-fledged computer simulation of the pertinent equations governing the electron-hole plasma therefore becomes imperative.

Vonroos, O.↗

NASA NTRS1989

Galerkin/least-squares procedures in computational fluid dynamics

The application of Galerkin/least-squares methodology to computational fluid dynamics is reviewed. Applications to advective-diffusive systems of equations and incompressible and compressible Euler and Navier-Stokes equations are described. Space-time finite-element methodology for time-dependent problems is emphasized and the status and role of so-called 'discontinuity-capturing operators' is reviewed.

Hughes, Thomas J. R.↗

NASA NTRS1991

Time-derivative preconditioning for viscous flows

A time-derivative preconditioning algorithm that is effective over a wide range of flow conditions from inviscid to very diffusive flows and from low speed to supersonic flows was developed. This algorithm uses a viscous set of primary dependent variables to introduce well-conditioned eigenvalues and to avoid having a nonphysical time reversal for viscous flow. The resulting algorithm also provides a mechanism for controlling the inviscid and viscous time step parameters to be of order one for very diffusive flows, thereby ensuring rapid convergence at very viscous flows as well as for inviscid flows. Convergence capabilities are demonstrated through computation of a wide variety of problems.

Choi, Yunho↗

NASA NTRS1991

Time-derivative preconditioning for viscous flows

Choi, Yunho↗

NASA NTRS2016

Preferential Heating of Oxygen 5+ Ions by Finite-Amplitude Oblique Alfven Waves

Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5+ ions by large-scale finite-amplitude Alfven waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5+ ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfven-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles theta less than or equal to 30 deg. The obliquely propagating Alfven pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.

Maneva, Yana G.↗

DOE OSTI2018

OPTIMIZATION OF FAST FISSION GAS RELEASE MODEL PARAMETERS USING MACHINE LEARNING ACCELERATED EVOLUTIONARY ALGORITHMS

Fission gas release models used in FRAPCON and FAST are based on equations that describe the physical transport phenomena that control release from the fuel pellet to the rod void volume of fission gas produced during irradiation. Although these models are physically based, many of the parameters that control the release such as the impact of burnup on diffusion and bubble formation and saturation are critical to the accurate prediction of fission gas release, but cannot be measured in any quantitative way. To overcome this problem, these parameters have been empirically derived to provide a best fit to the available data that includes rod puncture data and electron probe microanalysis (EMPA) and X-ray fluorescence (XRF) data of radial location of residual gas within the pellet. In general, a heuristic approach for deriving these parameters via practical implementation of expert analyses on just a few dependent variables has been shown to perform quite well in most situations. In this paper, we present an approach for applying machine learning to expedite model development. This approach is based on developing a deep artificial neural network which describes FRAPCON’s fission gas release models, and optimizes parameters by a differential evaluation algorithm. This approach allows us to quickly and accurately tune physical models based on expert judgment, and works as a human-in-the-loop approach, to assist the modeler in identifying and addressing regions of high uncertainty in a multi-parameter space. Results of the updated fission gas release model will be shown for all the FAST assessment data.

Johns, Jesse M.↗

NASA NTRS1998

Theoretical Evaluation of the Radiation Hazards from Cosmic Rays Within Space Vehicles

We may summarize our efforts as follows: a. Improvement of our calculations of the radial dose distribution from delta rays ejected in the passage of heavy ions through matter through the application of new data to a previous calculation by Kobetich and Katz (1968). Supplementing this calculation, we have found the radial distribution of electron energy spectra and the radial distribution of microdosimetric quantities (Cucinotta et al, 1996, 1997). b. Extension of the Katz model of cellular survival to bacteria, to lethal mutations in C. Elegans in vivo, to mutation induction in vitro, to thindown in radiobiology (observed experimentally at GSI, Darmstadt, and there called "Darmstadt hooks", predicted by Katz theory years before GSI was constructed). c. Coupling the Katz theory of RBE to the NASA theory of the diffusion of heavy ion beams in matter to yield predictions of the effects for monoenergetic heavy ion beams as well as range modulated beams used for cancer therapy. Here we have directed attention to the role of "ion-kill" (the effects produced by heavy ions passing through the nucleus of a cell), responsible for increased RBE, decreased OER, and reduced repair. We predict that the use of beams of heavy ions in cancer therapy will create late effect problems for fractionated therapy. We highlight also the damage by "ion-kill", from single heavy ions in the cosmic rays, to the central nervous system in space flight. d. The coupling of Katz theory and the NASA theory of heavy ion diffusion and penetration through matter, and knowledge of the space radiation environment, has been applied to design of shielding, to the cell damage in space flight.

Katz, Robert↗

DOE OSTI2022

Minimum stabilizing energy release for mixing processes

Diffusive operations, which mix the populations of different elements of phase space, can irreversibly transform a given initial state into any of a spectrum of different states from which no further energy can be extracted through diffusive operations. We call these ground states. Here, the lower bound of accessible ground-state energies represents the maximal possible release of energy. This lower bound, sometimes called the diffusively accessible free energy, is of interest in theories of instabilities and wave-particle interactions. On the other hand, the upper bound of accessible ground-state energies has escaped identification as a problem of interest. Yet, as demonstrated here, in the case of a continuous system, it is precisely this upper bound that corresponds to the paradigmatic “quasilinear plateau” ground state of the bump-on-tail distribution. Although for general discrete systems the complexity of calculating the upper bound grows rapidly with the number of states, using techniques adapted from treatments of the lower bound, the upper bound can in fact be computed directly for the three-state discrete system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1985

Pattern selection in a boundary-layer model of dendritic growth in the presence of impurities

Presently analyzed, in the context of a boundary-layer model, is the problem of pattern selection in dendritic growth in a situation where impurities are present in the undercooled liquid. It is found that the tip-velocity selection criterion that has been proposed recently for the geometrical model and the boundary-layer model of a pure substance can be extended, in a nontrivial way, to this more complex situation where two coupled diffusion fields (temperature and solute) determine the interface dynamics. This model predicts a sharp enhancement of tip velocity in good qualitative agreement with experiment. This agreement is consistent with the conjecture that a solvability condition can be used to determine the operating point of the dendrite in the full nonlocal problem.

Karma, A.↗

DOE OSTI2019

Countercations Control Local Specific Bonding Interactions and Nucleation Mechanisms in Concentrated Water-in-Salt Solutions

One of the continuing challenges presented in aqueous salt solutions is to understand ion association reactions driving dynamic demixing from solvation, complexation and solute clustering. The problems understanding this phenomenon are exacerbated in the highly-concentrated water-in-salt electrolytes, where the deficiency of water leads to a dramatic retardation of solvating water and formation of extended solvent-solute clustering networks. By probing microscopic dynamics of water and prenucleation clusters using quasi-elastic neutron scattering and proton nuclear magnetic resonance spectroscopy, we observed contrasting mechanistic specifics of ion-water mobilities in highly-concentrated Na+- vs. K+-based aluminate solutions (diffusion coefficients of 0.2 vs. 2.6?10-10m2s-1 at 293 K, respectively). The magnitude of the differences is far beyond the kosmotropic (Na+) and chaotropic (K+) effects typically understood at low ionic concentrations. The distinct crystallization mechanisms observed in analogous deuterated systems further imply that different prenucleation cluster dynamics and network packing morphologies can either frustrate or promote crystallization, as described by non-classical nucleation theory.

Wang, Hsiu-Wen↗

DOE OSTI2014

Limitation of parallel flow in double diffusive convection: Two- and three-dimensional transitions in a horizontal porous domain

Two-dimensional (2D) and three-dimensional (3D) numerical simulations of double diffusion natural convection in an elongated enclosure filled with a binary fluid saturating a porous medium are carried out in the present work. The Boussinesq approximation is made in the formulation of the problem, and Neumann boundary conditions for temperature and concentration are adopted, respectively, on vertical and horizontal walls of the cavity. The used numerical method is based on the control volume approach, with the third order quadratic upstream interpolation scheme in approximating the advection terms. A semi implicit method algorithm is used to handle the velocity-pressure coupling. To avoid the excessively high computer time inherent to the solution of 3D natural convection problems, full approximation storage with full multigrid method is used to solve the problem. A wide range of the controlling parameters (Rayleigh-Darcy number Ra, lateral aspect ratio Ay, Lewis number Le, and the buoyancy ration N) is investigated. We clearly show that increasing the depth of the cavity (i.e., the lateral aspect ratio) has an important effect on the flow patterns. The 2D perfect parallel flows obtained for small lateral aspect ratio are drastically destabilized by increasing the cavity lateral dimension. This yields a 3D fluid motion with a much more complex flow pattern and the usually considered 2D parallel flow model cannot be applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Generic expansion of the Jastrow correlation factor in polynomials satisfying symmetry and cusp conditions

Jastrow correlation factors play an important role in quantum Monte Carlo calculations. Together with an orbital based antisymmetric function, they allow the construction of highly accurate correlation wave functions. In this paper, a generic expansion of the Jastrow correlation function in terms of polynomials that satisfy both the electron exchange symmetry constraint and the cusp conditions is presented. In particular, an expansion of the three-body electron-electron-nucleus contribution in terms of cuspless homogeneous symmetric polynomials is proposed. The polynomials can be expressed in fairly arbitrary scaling function allowing a generic implementation of the Jastrow factor. It is demonstrated with a few examples that the new Jastrow factor achieves 85%–90% of the total correlation energy in a variational quantum Monte Carlo calculation and more than 90% of the diffusion Monte Carlo correlation energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Imaging properties and its improvements of scanning/imaging x-ray microscope

A scanning / imaging X-ray microscope (SIXM) system has been developed at SPring-8. The SIXM consists of a scanning X-ray microscope with a one-dimensional (1D) X-ray focusing device and an imaging (full-field) X-ray microscope with a 1D X-ray objective. The motivation of the SIXM system is to realize a quantitative and highly-sensitive multimodal 3D X-ray tomography by taking advantages of both the scanning X-ray microscope using multi-pixel detector and the imaging X-ray microscope. Data acquisition process of a 2D image is completely different between in the horizontal direction and in the vertical direction; a 1D signal is obtained with the linear-scanning while the other dimensional signal is obtained with the imaging optics. Such condition have caused a serious problem on the imaging properties that the imaging quality in the vertical direction has been much worse than that in the horizontal direction. In this paper, two approaches to solve this problem will be presented. One is introducing a Fourier transform method for phase retrieval from one phase derivative image, and the other to develop and employ a 1D diffuser to produce an asymmetrical coherent illumination.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2015

On origin and destruction of relativistic dust and its implication for ultrahigh energy cosmic rays

Dust grains may be accelerated to relativistic speeds by radiation pressure, diffusive shocks, and other acceleration mechanisms. Such relativistic grains have been suggested as primary particles of ultrahigh energy cosmic rays (UHECRs). In this paper, we first revisit the problem of acceleration by radiation pressure and calculate maximum grain velocities achieved. We find that grains can be accelerated to relativistic speeds with Lorentz factor γ<2 by powerful radiation sources, which is lower than earlier estimates showing that γ could reach ∼10. We then investigate different destruction mechanisms for relativistic grains traversing a variety of environments. In solar radiation, we find that the destruction by thermal sublimation and Coulomb explosions is important. We also quantify grain destruction due to electronic sputtering by ions and grain–grain collisions. Electronic sputtering is found to be rather inefficient, whereas the evaporation following grain–grain collisions is shown to be an important mechanism for which the a ∼ 0.01–1 μm grains would be destroyed after sweeping a gas column N{sub coll}∼5×10{sup 19}–5×10{sup 20} cm{sup −2}. Relativistic dust in the interstellar medium and intergalactic medium (IGM) would be disrupted by Coulomb explosions due to collisional charging after traversing a gas column N{sub coll}∼10{sup 17} cm{sup −2} unless grain material is very strong. We show that photoelectric emission by optical and ultraviolet background radiation is also significant for the destruction of relativistic dust in the IGM. The obtained results indicate that relativistic dust from galaxies would be destroyed before reaching the Earth's atmosphere and unlikely to account for UHECRs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Resonance broadened quasi-linear (RBQ) model for fast ion distribution relaxation due to Alfvénic eigenmodes

The burning plasma performance is limited by the confinement of the super-alfvénic fusion products such as alpha particles and the auxiliary heating ions capable of exciting the Alfvénic eigenmodes (AEs) (Gorelenkov et al 2014 Nucl. Fusion 54 125001). Here in this work the effect of AEs on fast ions is formulated within the quasi-linear (QL) theory generalized for this problem recently (Duarte 2017 PhD Thesis University of São Paulo, Brazil). The generalization involves the resonance line broadened interaction of energetic particles (EP) with AEs supplemented by the diffusion coefficients depending on EP position in the velocity space. A new resonance broadened QL code (or RBQ1D) based on this formulation allowing for EP diffusion in radial direction is built and presented in details. In RBQ1D applications we reduce the wave particle interaction (WPI) dynamics to 1D case when the particle kinetic energy is nearly constant. The diffusion equation for EP distribution evolution is then solved simultaneously for all particles along the angular momentum direction. We make initial applications of the RBQ1D to a DIII-D plasma with elevated q-profile where the beam ions show stiff transport properties (Collins et al (The DIII-D Team) 2016 Phys. Rev. Lett. 116 095001). In conclusion, AE driven fast ion profile relaxation is studied for validations of the QL approach in realistic conditions of beam ion driven instabilities in DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Multilevel Sequential Monte Carlo Samplers for Normalizing Constants

This article considers the sequential Monte Carlo (SMC) approximation of ratios of normalizing constants associated to posterior distributions which in principle rely on continuum models. Therefore, the Monte Carlo estimation error and the discrete approximation error must be balanced. A multilevel strategy is utilized to substantially reduce the cost to obtain a given error level in the approximation as compared to standard estimators. Two estimators are considered and relative variance bounds are given. The theoretical results are numerically illustrated for two Bayesian inverse problems arising from elliptic partial differential equations (PDEs). The examples involve the inversion of observations of the solution of (i) a 1-dimensional Poisson equation to infer the diffusion coefficient, and (ii) a 2-dimensional Poisson equation to infer the external forcing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Diffusion bonding of tungsten-vanadium-zirconium using vacuum hot pressing for the development of a low decay heat cladding solution for tungsten spallation targets

Tantalum has been used as cladding material for water-cooled solid tungsten targets at many leading spallation neutron production facilities thanks to its high neutron yield, manageable radiation damage behavior, and excellent corrosion/erosion resistance in radiation environments. However, from a safety hazard perspective, thermal neutron capture of tantalum in spallation environments causes a high specific decay heat in the target volume, which often becomes a limiting factor in increasing the beam power on the target. In this paper, we studied vacuum hot pressing (VHP) parameters to diffusion bond zirconium to tungsten to explore the feasibility of using zirconium alloys as an alternative cladding material to tantalum. Zirconium alloys have long been used as cladding material for early generation solid spallation targets, and nuclear fuel rods. In spallation environments zirconium has significantly lower decay heat with shorter decay time compared to tantalum. The hot isostatic pressing (HIP) of zirconium and tungsten is known to produce limited bonding quality due to the formation of the brittle ZrW2 intermetallic layer. To overcome this problem, placing a vanadium interlayer between tungsten and zirconium has been proposed by exploring parameter space in binary alloy phase diagrams. Under the VHP conditions, 860 ◦C at 70 MPa for 4 hours, Zr-V and V-W showed good diffusion bonding, which demonstrates the feasibility of a single step HIP process to make the zirconium alloy clad tungsten spallation volumes.

Lee, Yong Joong↗