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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Generalizing the self-healing diffusion Monte Carlo approach to finite temperature: A path for the optimization of low-energy many-body bases

A statistical method is derived for the calculation of thermodynamic properties of many-body systems at low temperatures. This method is based on the self-healing diffusion Monte Carlo method for complex functions [F. A. Reboredo, J. Chem. Phys. 136, 204101 (2012)] and some ideas of the correlation function Monte Carlo approach [D. M. Ceperley and B. Bernu, J. Chem. Phys. 89, 6316 (1988)]. In order to allow the evolution in imaginary time to describe the density matrix, we remove the fixed-node restriction using complex antisymmetric guiding wave functions. In the process we obtain a parallel algorithm that optimizes a small subspace of the many-body Hilbert space to provide maximum overlap with the subspace spanned by the lowest-energy eigenstates of a many-body Hamiltonian. We show in a model system that the partition function is progressively maximized within this subspace. We show that the subspace spanned by the small basis systematically converges towards the subspace spanned by the lowest energy eigenstates. Possible applications of this method for calculating the thermodynamic properties of many-body systems near the ground state are discussed. The resulting basis can also be used to accelerate the calculation of the ground or excited states with quantum Monte Carlo.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Experimental study of the nonlinear diffusion of a magnetic field and skin explosion of cylindrical conductors

The paper presents the results of an experimental study of the skin explosion of cylindrical conductors of diameter 1–3 mm (copper, aluminum, titanium, steel 3, and stainless steel) at a peak magnetic field of 200–600 T. The experiments were carried out on the MIG pulsed power generator at a current of up to 2.5 MA and a current rise time of 100 ns. The surface explosion of a conductor was identified by the appearance of a flash of extreme ultraviolet radiation. A minimum magnetic induction has been determined below which no plasma is generated at the conductor surface. For copper, aluminum, steel 3, titanium, and stainless steel, the minimum magnetic induction has been estimated to be (to within 10%) 375, 270, 280, 220, and 245 T, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗
DOE OSTI2017

Shock Dynamics in Stellar Outbursts. I. Shock Formation

Wave-driven outflows and non-disruptive explosions have been implicated in pre-supernova outbursts, supernova impostors, luminous blue variable eruptions, and some narrow-line and superluminous supernovae. To model these events, we investigate the dynamics of stars set in motion by strong acoustic pulses and wave trains, focusing on nonlinear wave propagation, shock formation, and an early phase of the development of a weak shock. We identify the shock formation radius, showing that a heuristic estimate based on crossing characteristics matches an exact expansion around the wave front and verifying both with numerical experiments. Our general analytical condition for shock formation applies to one-dimensional motions within any static environment, including both eruptions and implosions. We also consider the early phase of shock energy dissipation. We find that waves of super-Eddington acoustic luminosity always create shocks, rather than damping by radiative diffusion. Therefore, shock formation is integral to super-Eddington outbursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On the Effects of Pickup Ion-driven Waves on the Diffusion Tensor of Low-energy Electrons in the Heliosphere

The effects of Alfvén cyclotron waves generated due to the formation in the outer heliosphere of pickup ions on the transport coefficients of low-energy electrons is investigated here. To this end, parallel mean free path (MFP) expressions are derived from quasilinear theory, employing the damping model of dynamical turbulence. These are then used as inputs for existing expressions for the perpendicular MFP and turbulence-reduced drift coefficient. Using outputs generated by a two-component turbulence transport model, the resulting diffusion coefficients are compared with those derived using a more typically assumed turbulence spectral form, which neglects the effects of pickup ion-generated waves. It is found that the inclusion of pickup ion effects greatly leads to considerable reductions in the parallel and perpendicular MFPs of 1–10 MeV electrons beyond ∼10 au, which are argued to have significant consequences for studies of the transport of these particles.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Variable boundary II heat conduction

Computer program for solving both transient and steady-state heat transfer problems is presented. Specific applications of computer program are described. Formulation for individual nodes of solid medium for heat balance is presented. Diffusion equation is solved for all nodes simultaneously at finite increments of time.

Gramer, J.↗
NASA NTRS1977

Model for premixed turbulent flames

Flame speed and mean degree of reaction profiles are determined for an ensemble mean one-dimensional premixed turbulent flame. Elementary closure models are employed to treat turbulence closure problems arising from the nonlinear convective and chemical reaction rate terms in the time-dependent fuel conservation equation; a constant eddy diffusivity is useed to model correlation terms derived from the convective term. An expression developed for fuel consumption is solved to obtain the turbulent flame speed and the variation of mean fuel mass fraction across the flame. Comparison of the numerical analysis with experimental data indicates the importance of the quenching of the combustion reaction in high-wavenumber eddies.

Gouldin, F. C.↗
NASA NTRS1981

High power densities from high-temperature material interactions

Thermionic energy conversion (TEC) and metallic-fluid heat pipes (MFHPs), offering unique advantages in terrestrial and space energy processing by virtue of operating on working-fluid vaporization/condensation cycles that accept great thermal power densities at high temperatures, share complex materials problems. Simplified equations are presented that verify and solve such problems, suggesting the possibility of cost-effective applications in the near term for TEC and MFHP devices. Among the problems discussed are: the limitation of alkali-metal corrosion, protection against hot external gases, external and internal vaporization, interfacial reactions and diffusion, expansion coefficient matching, and creep deformation.

Morris, J. F.↗
NASA NTRS1985

Upwind relaxation algorithms for the Navier Stokes equations

The development of upwind relaxation algorithms for obtaining efficient steady-state solutions to the compressible Navier-Stokes equations is described. The method is second-order accurate spatially and naturally disipative, using third-order flux splitting of the pressure and convective terms and second-order central differencing for shear and heat flux terms. A line Gauss-Seidel relaxation approach, shown to be unconditionally stable for model convection and diffusion equations, is used. The algorithm is demonstrated for several flows using the thin-layer form of the equations, including the problem of shock-induced separation over a flat plate.

Thomas, J. L.↗
NASA NTRS1988

Joints in fiber-reinforced aluminum alloys

Problems of joining involving structural components made of aluminum alloys with boron-fiber reinforcements are discussed, giving attention to a diffusion-welding process. The tension characteristics and weight factors in the case of various types of welded joints are considered. Diffusion-welding equipment used in the experimental investigation is described. The strength characteristics obtained in various cases of weld and component design are examined, taking into account static and dynamic stresses.

Gerber, K.↗
NASA NTRS1995

Momentum deposition on Wolf-Rayet winds: Nonisotropic diffusion with effective gray opacity

We derive the velocity and mass-loss rate of a steady state Wolf-Rayet (WR) wind, using a nonisotropic diffusion approximation applied to the transfer between strongly overlapping spectral lines. Following the approach of Friend & Castor (1983), the line list is assumed to approximate a statistically parameterized Poisson distribution in frequency, so that photon transport is controlled by an angle-dependent, effectively gray opacity. We show the nonisotropic diffusion approximation yields good agreement with more accurate numerical treatments of the radiative transfer, while providing analytic insight into wind driving by multiple scattering. We illustrate, in particular, that multiple radiative momentum deposition does not require that potons be repeatedly reflected across substantial distances within the spherical envelope, but indeed is greatest when photons undergo a nearly local diffusion, e.g., through scattering by many lines closely spaced in frequency. Our results reiterate the view that the so-called 'momentum problem' of Wolf-Rayet winds is better characterized as an 'opacity problem' of simply identfying enough lines. One way of increasing the number of thick lines in Wolf-Rayet winds is to transfer opacity from saturated to unsaturated lines, yielding a steeper opacity distribution than that found in OB winds. We discuss the implications of this perspective for extending our approach to W-R wind models that incorporate a more fundamental treatment of the ionization and excitation processes that determine the line opacity. In particular, we argue that developing statistical descriptions of the lines to allow an improved effective opacity for the line ensemble would offer several advantages for deriving such more fundamental W-R wind models.

Gayley, Kenneth G.↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS1999

A Gas-Kinetic Method for Hyperbolic-Elliptic Equations and Its Application in Two-Phase Fluid Flow

A gas-kinetic method for the hyperbolic-elliptic equations is presented in this paper. In the mixed type system, the co-existence and the phase transition between liquid and gas are described by the van der Waals-type equation of state (EOS). Due to the unstable mechanism for a fluid in the elliptic region, interface between the liquid and gas can be kept sharp through the condensation and evaporation process to remove the "averaged" numerical fluid away from the elliptic region, and the interface thickness depends on the numerical diffusion and stiffness of the phase change. A few examples are presented in this paper for both phase transition and multifluid interface problems.

Xu, Kun↗
NASA NTRS1995

Design of a three-layer antireflection coating for high efficiency indium phosphide solar cells using a chemical oxide as first layer

It is well known that the behavior of III-V compound based solar cells is largely controlled by their surface, since the majority of light generated carriers (63% for GaAs and 79% for InP) are created within 0.2 mu m of the surface of the illuminated cell. Consequently, the always observed high surface recombination velocity (SRV) on these cells is a serious limiting factor for their high efficiency performance, especially for those with p-n junction made by either thermal diffusion or ion implantation. A good surface passivation layer, ideally a grown oxide as opposed to a deposited one, will cause a significant reduction in the SRV without adding interface problems, thus improving the performance of III-V compound based solar cells. Another significant benefit to the overall performance of the solar cells can be achieved by a substantial reduction of their large surface optical reflection by the use of a well designed antireflection (AR) coating. In this paper, we demonstrate the effectiveness of using a chemically grown thermally and chemically stable oxide, not only for surface passivation but also as an integral part of a 3-layer AR coating for thermally diffused p+n InP solar cells. A phosphorus-rich interfacial oxide, In(PO3)3, is grown at the surface of the p+ emitter using an etchant based on HNO3, o-H3PO4 and H2O2. This oxide has the unique properties of passivating the surface as well as serving as an efficient antireflective layer yielding a measured record high AMO open-circuit voltage of 890.3 mV on a thermally diffused InP(Cd,S) solar cell. Unlike conventional single layer AR coatings such as ZnS, Sb2O3, SiO or double layer AR coatings such as ZnS/MgF2 deposited by e-beam or resistive evaporation, this oxide preserves the stoichiometry of the InP surface. We show that it is possible to design a three-layer AR coating for a thermally diffused InP solar cell using the In(PO3)3 grown oxide as the first layer and Al2O3 and MgF2 as the second and third layers respectively, so as to yield an overall theoretical reflectance of less than 2%. Since chemical oxides are readily grown on III-V semiconductors materials, the technique of using the grown oxide layer to both passivate the surface as well as serve as the first of a multilayer AR coating should work well for all III-V compound-based solar cells.

Moulot, Jacques↗
NASA NTRS1995

Design of a Three-Layer Antireflection Coating for High Efficiency Indium Phosphide Solar Cells Using a Chemical Oxide as First Layer

It is well known that the behavior of III-V compound based solar cells is largely controlled by their surface, since the majority of light generated carriers (63% for GaAs and 79% for InP) are created within 0.2 microns of the illuminated surface of the cell. Consequently, the always observed high surface recombination velocity (SRV) on these cells is a serious limiting factor for their high efficiency performance, especially for those with the p-n junction made by either thermal diffusion or ion implantation. A good surface passivation layer, ideally, a grown oxide as opposed to a deposited one, will cause a significant reduction in the SRV without adding interface problems, thus improving the performance of III-V compound based solar cells. Another significant benefit to the overall performance of the solar cells can be achieved by a substantial reduction of their large surface optical reflection by the use of a well designed antireflection (AR) coating. In this paper, we demonstrate the effectiveness of using a chemically grown, thermally and chemically stable oxide, not only for surface passivation but also as an integral part of a 3- layer AR coating for thermally diffused p(+)n InP solar cells. A phosphorus-rich interfacial oxide, In(PO3)3, is grown at the surface of the p(+) emitter using an etchant based on HNO3, o-H3PO4 and H2O2. This oxide has the unique properties of passivating the surface as well as serving as a fairly efficient antireflective layer yielding a measured record high AM0, 25 C, open-circuit voltage of 890.3 mV on a thermally diffused InP(Cd,S) solar cell. Unlike conventional single layer AR coatings such as ZnS, Sb2O3, SiO or double layer AR coatings such as ZnS/MgF2 deposited by e-beam or resistive evaporation, this oxide preserves the stoichiometry of the InP surface. We show that it is possible to design a three-layer AR coating for a thermally diffused InP solar cell using the In(PO3)3 grown oxide as the first layer and Al2O3, MgF2 or ZnS, MgF2 as the second and third layers respectively, so as to yield an overall theoretical reflectance of less than 2%. Since chemical oxides are readily grown on III-V semiconductor materials, the technique of using the grown oxide layer to both passivate the surface as well as serve as the first of a multilayer AR coating, should work well for essentially all III-V compound-based solar cells.

Moulot, Jacques↗
DOE OSTI2024

GraMeR: Gra ph Me ta R einforcement learning for multi-objective influence maximization

Influence maximization (IM) is a combinatorial problem of identifying a subset of seed nodes in a network (graph), which when activated, provide a maximal spread of influence in the network for a given diffusion model and a budget for seed set size. IM has numerous applications such as viral marketing, epidemic control, sensor placement and other network-related tasks. However, its practical uses are limited due to the computational complexity of current algorithms. Recently, deep reinforcement learning has been leveraged to solve IM in order to ease the computational burden. However, there are serious limitations in current approaches, including narrow IM formulation that only consider influence via spread and ignore self-activation, low scalability to large graphs, and lack of generalizability across graph families leading to a large running time for every test network. In this work, we address these limitations through a unique approach that involves: (1) Formulating a generic IM problem as a Markov decision process that handles both intrinsic and influence activations; (2)incorporating generalizability via meta-learning across graph families. There are previous works that combine deep reinforcement learning with graph neural network, but this work solves a more realistic IM problem and incorporates generalizability across graphs via meta reinforcement learning. Extensive experiments are carried out in various standard networks to validate performance of the proposed Graph Meta Reinforcement learning (GraMeR) framework. Finally, the results indicate that GraMeR is multiple orders faster and generic than conventional approaches when applied on small to medium scale graphs.

97 MATHEMATICS AND COMPUTING↗