Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion Couple”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1997

Wave Driven Exothermic Heating in the Mesopause Region

A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.

Hickey, Michael P.↗
NASA NTRS1975

Diurnal variations in the thermosphere. I - Theoretical formulation

A nonlinear perturbation theory is formulated for the solution of the multicomponent equations of energy, mass, and momentum conservation in the atmosphere. The theory is three-dimensional and includes the effects of heat conduction and advection, viscosity, ion drag, and diffusion. The theory is described as a superposition of mathematical modes obtained by expanding the physical quantities into vector and spherical harmonics. The coupling between the various modes, both linear and nonlinear, is included. The theory provides a basis for the treatment of the thermosphere and its interaction with the lower atmosphere, where 'mode coupling' is most important. As an example, a comparison is presented between one-dimensional and three-dimensional calculations of the fundamental mode of the diurnal component in the thermosphere. Coupling between the lowest modes is considered to describe the physical conditions of the lower thermosphere where inertia and Coriolis forces become dominant over the ion-drag and viscous forces. In this region, the latitude structures of the temperature, wind field, and diffusively controlled oxygen are shown to change significantly.

Harris, I.↗
DOE OSTI2025

Detection of atomic oxygen and its electronic coherence decays using time-resolved ultrafast coherent Raman scattering

We report the detection of atomic oxygen and quantitative measurements of its electronic Raman coherence decays in flames and low-temperature plasmas using time-resolved hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (CARS). Atomic oxygen was detected using the Raman transitions between the spin–orbit coupled triplet ground states. Atomic oxygen was generated in an H 2 /O 2 /Ar diffusion flame and an O 2 /Ar pulsed plasma discharge. Single exponential decays were observed for the O( 3 P 2 )–O( 3 P 1 ) Raman transition at 158.3 cm −1 and the O( 3 P 2 )–O( 3 P 0 ) Raman transition at 227 cm −1 . From the decay measurements, the atomic O Raman linewidths were obtained from 25 to 150 Torr in non-equilibrium plasma and at 760 Torr in a flame. Enhanced signal-to-noise ratios (SNRs) of atomic oxygen and atomic to molecular oxygen signal contrasts were obtained by taking advantage of electronic triplet coherence beating. Enhancement of up to seven times in the atomic O SNR was observed. Furthermore, we also found that the dephasing rates of O 2 (v = 0–3, N = 37) were similar, which provides evidence for the assumption that vibrational excitation does not influence the dephasing of diatomic molecular rotational CARS transitions.

Atomic and molecular spectroscopy↗
NASA NTRS1974

Solutions of the Navier-Stokes equations for vortex breakdown

Steady solutions of the Navier-Stokes equations, in terms of velocity and pressure, for breakdown in an unconfined viscous vortex are obtained numerically using the artificial compressibility technique of Chorin combined with an ADI finite-difference scheme. Axisymmetry is assumed and boundary conditions are carefully applied at the boundaries of a large finite region in an axial plane while resolution near the axis is maintained by a coordinate transformation. The solutions, which are obtained for Reynolds numbers up to 200 based on the free-stream axial velocity and a characteristic core radius, show that breakdown results from the diffusion and convection of vorticity away from the vortex core which, because of the strong coupling between the circumferential and axial velocity fields in strongly swirling flows, can lead to stagnation and reversal of the axial flow near the axis.

Grabowski, W. J.↗
NASA NTRS1974

Study of glass preforms for glass fiber optics applications (study of space processing of ceramic materials)

The feasibility, and technical and economic desirability was studied of space processing of glass preforms for optical fiber transmission applications. The results indicate that space processing can produce glass preforms of equal quality at lower cost than earth bound production, and can produce diameter modulation in the glass preform which promotes mode coupling and lowers the dispersion. The glass composition can be modified through the evaporative and diffusion processes, and graded refractive index profiles can be produced. A brief summary of the state of the art in optical fiber transmission is included.

Wang, F. F. Y.↗
NASA NTRS1986

Observations of ionospheric magnetospheric coupling - DE and Chatanika coincidences

Observations from several experiments on board the Dynamics Explorer 1 and 2 (DE 1 and 2) spacecraft and ground-based radar measurements from the Chatanika radar are combined in order to examine the details of ionospheric/magnetospheric coupling in the local evening sector. DE 1 and DE 2 were in coplanar polar orbits that provided measurements almost simultaneously in time and magnetically coincident with the Chatanika radar from L = 3 to L = 17. The coupling processes are inferred from the density, temperature, composition, and angular distributions of the low-energy plasma observed from the E region of the ionosphere to magnetospheric altitudes of 2.5 earth radii. Plasma characteristics of the plasmasphere, main trough, auroral zone, and polar cap can be studied in this data set. The observations imply that as L increases, the dominant coupling mechanism between the ionosphere and magnetosphere in the measured energy range changes from equilibrium diffusion to perpendicular acceleration and finally to parallel acceleration.

Green, J. L.↗
NASA NTRS1988

Cosmic-ray-modified stellar winds. II - A perturbation approach

A perturbation method is developed to describe the modification of a stellar wind flow with termination shock by the Galactic cosmic rays. The perturbation parameter epsilon used in the analysis is the ratio of the galactic cosmic-ray pressure to the thermal gas pressure at a large distance from the star. The positive galactic cosmic-ray pressure gradient serves to brake the outflowing stellar wind gas, particularly just upstream of the termination shock of the wind. A one-fluid polytropic model is used to describe the thermal stellar wind gas, with the coupling between the cosmic rays and the thermal gas being determined by an average hydrodynamical diffusion coefficient kappa describing the scattering of cosmic rays by magnetohydrodynamic waves traveling in the background flow. The analysis takes into account cosmic-ray modifications of the critical point of the wind, the thermal gas entropy constants on both sides of the shock, and the fluid velocity profile.

Ko, C. M.↗
NASA NTRS1991

Dynamics and control of hypersonic vehicles - The integration challenge for the 1990's

The dynamic characteristics of hypersonic vehicles are reviewed, with special attention to the interactions between the airframe, engine, and structural dynamics. Based on a preliminary investigation of the dynamics of a generic vehicle configuration similar to the X-30 with Scramjet propulsion, an assessment of these interactions is presented. The control effectors include aerodynamic pitch-control surfaces, as well as engine fuel flow and diffuser area ratio. The study configuration is statically instable in pitch, and exhibits strong airframe/engine/elastic coupling in the attitude dynamics and engine responses. This strong coupling will require a highly integrated airframe-engine control system, and the performance of the attitude control system will be contingent upon the ability to adequately deal with the structural aeroelastic response and engine dynamics.

Schmidt, David K.↗
NASA NTRS1992

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗
NASA NTRS1995

Single-drop reactive extraction/extractive reaction with forced convective diffusion and interphase mass transfer

An algorithm has been developed for time-dependent forced convective diffusion-reaction having convection by a recirculating flow field within the drop that is hydrodynamically coupled at the interface with a convective external flow field that at infinity becomes a uniform free-streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet, or reactions can take place in both phases. The algorithm has been implemented, and for comparison results are shown here for the case of no reaction in either phase and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Kleinman, Leonid S.↗
NASA NTRS2005

Removal of organic contaminants by RO and NF membranes

Rejection characteristics of organic and inorganic compounds were examined for six reverse osmosis (RO) membranes and two nanofiltration (NF) membranes that are commercially available. A batch stirred-cell was employed to determine the membrane flux and the solute rejection for solutions at various concentrations and different pH conditions. The results show that for ionic solutes the degree of separation is influenced mainly by electrostatic exclusion, while for organic solutes the removal depends mainly upon the solute radius and molecular structure. In order to provide a better understanding of rejection mechanisms for the RO and NF membranes, the ratio of solute radius (r(i,s)) to effective membrane pore radius (r(p)) was employed to compare rejections. An empirical relation for the dependence of the rejection of organic compounds on the ratio r(i,s)/r(p) is presented. The rejection for organic compounds is over 75% when r(i,s)/r(p) is greater than 0.8. In addition, the rejection of organic compounds is examined using the extended Nernst-Planck equation coupled with a steric hindrance model. The transport of organic solutes is controlled mainly by diffusion for the compounds that have a high r(i,s)/r(p) ratio, while convection is dominant for compounds that have a small r(i,s)/r(p) ratio. c2005 Elsevier B.V. All rights reserved.

NASA Discipline Life Sciences Technologies↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
DOE OSTI2024

Morphology selection in dealloying: A phase field study of the coupling among kinetic mechanisms

A phase field model is used to investigate how the corrosion rate and morphology of a dealloying metal is controlled by the kinetic rates of bulk solid diffusion, interfacial diffusivity, and interface dissolution. A specific focus is investigating the manner by which variations in these kinetic parameters synergistically influence the onset of a morphological instability at the solid-liquid interface, which is known to precede rapid dealloying corrosion towards highly topologically complex structures such as bicontinuous porous solids. This work demonstrates that this instability can be suppressed by either an increase of interdiffusion in the solid alloy, or by a decreasing rate of dissolution at the solid-liquid interface. When the instability is suppressed, the result is a relatively planar, passivated dealloying front. Importantly, however, coupled changes in kinetic parameters can also promote the morphological instability and sustain fine-scale ligament and pore formation. Finally, the relevance of the solid-liquid interfacial diffusion for determining the rate of dealloying and ligament formation is highlighted. By elucidating the mechanistic influence of bulk and interfacial kinetic factors on the dealloyed morphology, this work helps to rationalize the contrasting morphologies observed in various dealloying systems.

Diffusion↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2025

A compartmentalized model of multiphase chemical kinetics

There are significant challenges in predicting multiphase chemical kinetics due to the complex coupling of reaction and mass transport across a phase boundary (i.e., interface). Here, we describe a framework for predicting multiphase kinetics that embeds the elementary kinetic steps of reaction, solvation, and diffusion into a coarse grain spatial description of two phases. The model is constructed to bridge the short-timescale interfacial dynamics observed in molecular simulations with the longer timescales observed in kinetic experiments. A simple set of governing differential equations is derived, which, when solved numerically or analytically, yield accurate predictions of multiphase kinetics in microdroplets. Although the equations are formulated for gas-liquid reactions, the underlying conceptual framework is general and can be applied to transformations in other two-phase systems (solid-liquid, liquid-liquid, etc.).

Chemical kinetics and dynamics↗
NASA NTRS1994

Erosional dynamics, flexural isostasy, and long-lived escarpments: A numerical modeling study

Erosional escarpments common features of high-elevation rifted continets. Fission track data suffest that these escarpments form by base level lowering and/or marginal uplift during rifting, followed by lateral retreat of an erosion front across tens to hundreds of kioometers. Previous modeling studies have shown that this characteristic pattern of denudation can have a profound impact upon marginal isostatic uplift and the evolution of offshore sedimentary basins. Yet at present there is only a rudimentary understanding of the geomorphic mechanisms capable of driving such prolonged escarpment retreat. In this study we present a nonlinear, two-dimensional landscape evolution model tha tis used to asses the necessary and sufficient conditions for long-term retreat of a rift-generated escarpment. The model represents topography as a grid of cells, with drainage networkds evolving as water flows across the grid in the direction of steepest descent. The model accounts for sediment production by weathering, fluvial sediment transport, bedrock channel erosion, and hillslope sediment transport by diffusive mechanisms and by mass failure. Numerical experiments presented explore the effects of different combinations of erosion processes and of dynamic coupling between denudation and flexural isostatic uplift. Model results suggest that the necessary and sufficient conditions for long-term escarpment retreat are (1) incising bedrock channels in which the erosion rate increases with increasing drainage area, so that the channels steepen and propagate headward; (2) a low rate of sediment production relative to sediment transport efficiency, which promotes relief-generating processes over diffusive ones; (3) high continental elevation, which allows greater freedom for fluvial dissection; and (4) any process, including flexural isostatic uplift, that helps to maintain a drainage divide near an escarpment crest. Flexural isostatic uplift also facilitates escarpment, thereby increasing channel gradients and accelerating erosion which in turn generates additional isostatic uplift. Of all the above conditions, high continental elevation is common ot most rift margin escarpments and may ultimately be the most important factor.

Tucker, Gregory E.↗