SEARCH · Engineering Papers

Results for “diffusion length”

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,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Effect of wave localization on plasma instabilities

The Anderson model of wave localization in random media is involved to study the effect of solar wind density turbulence on plasma processes associated with the solar type III radio burst. ISEE-3 satellite data indicate that a possible model for the type III process is the parametric decay of Langmuir waves excited by solar flare electron streams into daughter electromagnetic and ion acoustic waves. The threshold for this instability, however, is much higher than observed Langmuir wave levels because of rapid wave convection of the transverse electromagnetic daughter wave in the case where the solar wind is assumed homogeneous. Langmuir and transverse waves near critical density satisfy the Ioffe-Reigel criteria for wave localization in the solar wind with observed density fluctuations -1 percent. Numerical simulations of wave propagation in random media confirm the localization length predictions of Escande and Souillard for stationary density fluctations. For mobile density fluctuations localized wave packets spread at the propagation velocity of the density fluctuations rather than the group velocity of the waves. Computer simulations using a linearized hybrid code show that an electron beam will excite localized Langmuir waves in a plasma with density turbulence. An action principle approach is used to develop a theory of non-linear wave processes when waves are localized. A theory of resonant particles diffusion by localized waves is developed to explain the saturation of the beam-plasma instability. It is argued that localization of electromagnetic waves will allow the instability threshold to be exceeded for the parametric decay discussed above.

Levedahl, William Kirk↗

NASA NTRS2003

Thermal Characteristics and Structure of Fully-Modulated, Turbulent Diffusion Flames in Microgravity

Turbulent jet diffusion flames are studied in microgravity and normal gravity under fully-modulated conditions for a range of injection times and a 50% duty cycle. Diluted ethylene was injected through a 2-mm nozzle at a Reynolds number of 5,000 into an open duct, with a slow oxidizer co-flow. Microgravity tests are conducted in NASA's 2.2 Second Drop Tower. Flames with short injection times and high duty cycle exhibit a marked increase in the ensemble-averaged flame length due to the removal of buoyancy. The cycle-averaged centerline temperature profile reveals higher temperatures in the microgravity flames, especially at the flame tip where the difference is about 200 K. In addition, the cycle-averaged measurements of flame radiation were about 30% to 60% greater in microgravity than in normal gravity.

Hermanson, J. C.↗

NASA NTRS1974

Experimental determination of turbulence in a GH2-GOX rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a gaseous rocket combustion chamber have been determined from the experimental measurements of the tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and Spalding's numerical method for solving the conservation equations of fluid mechanics was used to calculate these quantities. Taylor's theory was extended to consider the inhomogeneity of the turbulence field in the axial direction of the combustion chamber. An exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the maximum value of the intensity of turbulence is about 15% and the Lagrangian correlation coefficient drops to about 0.12 in one inch of the chamber length.

Tou, P.↗

NASA NTRS1998

Critical Slowing Down in Time-to-Extinction: An Example of Critical Phenomena in Ecology

We study a model for two competing species that explicitly accounts for effects due to discreteness, stochasticity and spatial extension of populations. The two species are equally preferred by the environment and do better when surrounded by others of the same species. We observe that the final outcome depends on the initial densities (uniformly distributed in space) of the two species. The observed phase transition is a continuous one and key macroscopic quantities like the correlation length of clusters and the time-to-extinction diverge at a critical point. Away from the critical point, the dynamics can be described by a mean-field approximation. Close to the critical point, however, there is a crossover to power-law behavior because of the gross mismatch between the largest and smallest scales in the system. We have developed a theory based on surface effects, which is in good agreement with the observed behavior. The course-grained reaction-diffusion system obtained from the mean-field dynamics agrees well with the particle system.

Gandhi, Amar↗

NASA NTRS2000

Smoke-Point Properties of Non-Buoyant Round Laminar Jet Diffusion Flames

The laminar smoke-point properties of non-buoyant round laminar jet diffusion flames were studied emphasizing results from long-duration (100-230 s) experiments at microgravity carried out in orbit aboard the space shuttle Columbia. Experimental conditions included ethylene- and propane-fueled flames burning in still air at an ambient temperature of 300 K, pressures of 35-130 kPa, jet exit diameters of 1.6 and 2.7 mm, jet exit velocities of 170-690 mm/s, jet exit Reynolds numbers of 46-172, characteristic flame residence times of 40-302 ms, and luminous flame lengths of 15-63 mm. Contrary to the normal-gravity laminar smoke point, in microgravity, the onset of laminar smoke-point conditions involved two flame configurations: closed-tip flames with soot emissions along the flame axis and open-tip flames with soot emissions from an annular ring about the flame axis. Open-tip flames were observed at large characteristic flame residence times with the onset of soot emissions associated with radiative quenching near the flame tip: nevertheless, unified correlations of laminar smoke-point properties were obtained that included both flame configurations. Flame lengths at laminar smoke-point conditions were well correlated in terms of a corrected fuel flow rate suggested by a simplified analysis of flame shape. The present steady and non-buoyant flames emitted soot more readily than non-buoyant flames in earlier tests using ground-based microgravity facilities and than buoyant flames at normal gravity, as a result of reduced effects of unsteadiness, flame disturbances, and buoyant motion. For example, present measurements of laminar smoke-point flame lengths at comparable conditions were up to 2.3 times shorter than ground-based microgravity measurements and up to 6.4 times shorter than buoyant flame measurements. Finally, present laminar smoke-point flame lengths were roughly inversely proportional to pressure to a degree that is a somewhat smaller than observed during earlier tests both at microgravity (using ground-based facilities) and at normal gravity.

Urban, D. L.↗

NASA NTRS2000

Smoke-Point Properties of Nonbuoyant Round Laminar Jet Diffusion Flames

The laminar smoke-point properties of non-buoyant round laminar jet diffusion flames were studied emphasizing results from long-duration (100-230 s) experiments at microgravity carried out in orbit aboard the space shuttle Columbia. Experimental conditions included ethylene- and propane-fueled flames burning in still air at an ambient temperature of 300 K, pressures of 35-130 kPa, jet exit diameters of 1.6 and 2.7 mm, jet exit velocities of 170-690 mm/s, jet exit Reynolds numbers of 46-172, characteristic flame residence times of 40-302 ms, and luminous flame lengths of 15-63 mm. Contrary to the normal-gravity laminar smoke point, in microgravity the onset of laminar smoke-point conditions involved two flame configurations: closed-tip flames with soot emissions along the flame axis and open-tip flames with soot emissions from an annular ring about the flame axis. Open-tip flames were observed at large characteristic flame residence times with the onset of soot emissions associated with radiative quenching near the flame tip: nevertheless, unified correlations of laminar smoke-point properties were obtained that included both flame configurations. Flame lengths at laminar smoke-point conditions were well correlated in terms of a corrected fuel flow rate suggested by a simplified analysis of flame shape. The present steady and nonbuoyant flames emitted soot more readily than non-buoyant flames in earlier tests using ground-based microgravity facilities and than buoyant flames at normal gravity, as a result of reduced effects of unsteadiness, flame disturbances, and buoyant motion. For example, present measurements of laminar smokepoint flame lengths at comparable conditions were up to 2.3 times shorter than ground-based microgravity measurements and up to 6.4 times shorter than buoyant flame measurements. Finally, present laminar smoke-point flame lengths were roughly inversely proportional to pressure to a degree that is a somewhat smaller than observed during earlier tests both at microgravity (using ground-based facilities) and at normal gravity,

Urban, D. L.↗

DOE OSTI2023

New constraints on warm dark matter from the Lyman- α forest power spectrum

The forest of Lyman-α absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable—the 1D flux power spectrum (FPS)—should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard Λ⁢CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the graphics processing unit-accelerated code cholla to study the evolution of the Lyman-α forest under a wide range of physically motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts 4.0≲z≲5.2 [E. Boera et al., Revealing reionization with the thermal history of the intergalactic medium: New constraints from the Ly⁢α flux power spectrum, Astrophys. J. 872, 101 (2019)] yields a lower-limit m WDM >3.1 keV (95% C.L.) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of m WDM =4.5 keV. In conclusion, we find that the suppression of the FPS from free-streaming saturates at k≳0.1 s km -1 because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1973

Thermotransport in liquid aluminum-copper alloys

A thermotransport study was made on a series of liquid aluminum-copper alloys which contained from trace amounts to 33 weight percent copper. The samples in the form of narrow capillaries were held in known temperature gradient of thermotransport apparatus until the stationary state was reached. The samples were analyzed for the concentration of copper along the length. Copper was observed to migrate to the colder regions in all the samples. The heat of transport, Q*, was determined for each composition from a plot of concentration of copper versus reciprocal absolute temperature. The value of Q* is the highest at trace amounts of copper (4850 cal/gm-atom), but decreases with increasing concentration of copper and levels off to 2550 cal/gm-atom at about 25 weight percent copper. The results are explained on the basis of electron-solute interaction and a gas model of diffusion.

Bhat, B. N.↗

DOE OSTI2025

Ohmic contacts to nitrogen-incorporated n-type ultrananocrystalline diamond film grown on intrinsic single crystal diamond substrate

Nitrogen-incorporated ultrananocrystalline diamond (n-UNCD) films offer tremendous potential for diverse electronic applications. However, the absence of a reliable Ohmic contact at room temperature limits their practical integration and broader applicability. Here, in this study, we investigate the room temperature specific contact resistivity (ρ c ) of Ti/Pt/Au metal stack deposited on n-UNCD films grown on an intrinsic single crystal diamond substrate using a microwave plasma chemical vapor deposition system. We employ a circular transfer length model (c-TLM) and find the room temperature ρ c to be ∼4.67 × 10 −5 Ω cm 2 , which is among the lowest reported value for n-UNCD films. High temperature vacuum annealing conducted at 700 and 800 °C results in an initial improvement, followed by a minor degradation in ρ c values, respectively. The electrical contacts remain highly Ohmic for all measurements. Furthermore, cross-sectional transmission electron microscopy analysis suggests formation of conductive titanium carbide layer with no significant inter metallic diffusion. Overall, the electrical contacts demonstrate robust thermal stability, both of which are critical for attaining high-performance nanocrystalline diamond-based electronic devices.

Low resistance contacts↗

DOE OSTI2024

Deuterium Permeation Through a Self-Supported Palladium-Silver Membrane in Helium Gas Mixtures

Palladium membranes and membrane reactors can separate and purify tritium from impurities in the plasma exhaust processing system for the fusion energy fuel cycle. Membranes can also act as a continuous separation method to remove tritium from helium streams in the breeder blanket tritium extraction system, such as from the purge gas of solid breeders. To investigate the potential of these membranes for these applications, we performed a deuterium permeation campaign with a self-supported palladium-silver (Pd-25Ag wt%) tube of 0.15 m length, 76 µ m wall thickness, and 3.0 mm inner diameter. A gas mixture of 3.95% D 2 and a balance of He was delivered to the inside of the tube and D permeated radially outwards through the membrane into a vacuum chamber. Permeation experiments were conducted at 300 ° C, 350 ° C, 400° C, and 450°C, from 100 to 1000 sccm total flow rate, and with total pressures of 90, 150, 190, and 250 kPa. Further, deuterium permeation was consistently lower than predicted from diffusion-limited permeation models, thus we developed a transport model that included gas-phase mass transfer and surface reactions to model experimental results. The dissociation constant was optimized to fit the developed model to experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Deep learning assisted jet tomography for the study of Mach cones in QGP

Abstract Mach cones are expected to form in the expanding quark-gluon plasma (QGP) when energetic quarks and gluons traverse the hot medium at a velocity faster than the speed of sound in high-energy heavy-ion collisions. The shape of the Mach cone and the associated diffusion wake are sensitive to the initial jet production location and the propagation direction of the parton shower relative to the radial flow because of the distortion caused by the collective expansion of the QGP and the large density gradient. The shape of jet-induced Mach cones and their distortions in heavy-ion collisions provide a unique and direct probe of the dynamical evolution and the equation of state of QGP. However, it is difficult to identify the Mach cone and the diffusion wake in current experimental measurements of final hadron distributions because they are averaged over all possible initial jet production locations and parton-shower propagation directions. To overcome this difficulty, we develop a deep learning assisted jet tomography which uses the full information of the final hadrons from jets to localize the initial jet production positions. This method can help to constrain the initial regions of jet production in heavy-ion collisions and enable a differential study of Mach-cones with different path lengths and orientations relative to the radial flow of the QGP in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2004

Sensitivity Studies for In-Situ Automated Tape Placement of Thermoplastic Composites

This modeling effort seeks to improve the interlaminate bond strength of thermoplastic carbon composites produced by the in-situ automated tape placement (ATP) process. An existing high productivity model is extended to lower values of the Peclet number that correspond to the present operating conditions of the Langley ATP robot. (The Peclet number is the dimensionless ratio of inertial to diffusive heat transfer.) In sensitivity studies, all of the process and material parameters are individually varied. The model yields the corresponding variations in the effective bonding time (EBT) referred to the glass transition temperature. According to reptation theory, the interlaminate bond strength after wetting occurs is proportional to the one-fourth power of EBT. The model also computes the corresponding variations in the thermal input power (TIP) and the mass and volumetric process rates. Process studies show that a 10 percent increase in the consolidation length results in a 20 percent increase in EBT and a 5 percent increase in TIP. A surprising result is that a 10 K decrease in the tooling temperature results in a 25 percent increase in EBT and an 8 percent increase in TIP. Material studies show that a 10 K decrease in glass transition temperature results in an 8 percent increase in EBT and a 8 percent decrease in TIP. A 20 K increase in polymer degradation temperature results in a 23 percent increase in EBT with no change in TIP.

Costen, Robert C.↗

NASA NTRS2001

The Cool Flames Experiment

A space-based experiment is currently under development to study diffusion-controlled, gas-phase, low temperature oxidation reactions, cool flames and auto-ignition in an unstirred, static reactor. At Earth's gravity (1g), natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles via the Arrhenius temperature dependence of the reaction rates. Natural convection is important in all terrestrial cool flame and auto-ignition studies, except for select low pressure, highly dilute (small temperature excess) studies in small vessels (i.e., small Rayleigh number). On Earth, natural convection occurs when the Rayleigh number (Ra) exceeds a critical value of approximately 600. Typical values of the Ra, associated with cool flames and auto-ignitions, range from 104-105 (or larger), a regime where both natural convection and conduction heat transport are important. When natural convection occurs, it alters the temperature, hydrodynamic, and species concentration fields, thus generating a multi-dimensional field that is extremely difficult, if not impossible, to be modeled analytically. This point has been emphasized recently by Kagan and co-workers who have shown that explosion limits can shift depending on the characteristic length scale associated with the natural convection. Moreover, natural convection in unstirred reactors is never "sufficiently strong to generate a spatially uniform temperature distribution throughout the reacting gas." Thus, an unstirred, nonisothermal reaction on Earth does not reduce to that generated in a mechanically, well-stirred system. Interestingly, however, thermal ignition theories and thermokinetic models neglect natural convection and assume a heat transfer correlation of the form: q=h(S/V)(T(bar) - Tw) where q is the heat loss per unit volume, h is the heat transfer coefficient, S/V is the surface to volume ratio, and (T(bar) - Tw ) is the spatially averaged temperature excess. This Newtonian form has been validated in spatially-uniform, well-stirred reactors, provided the effective heat transfer coefficient associated with the unsteady process is properly evaluated. Unfortunately, it is not a valid assumption for spatially-nonuniform temperature distributions induced by natural convection in unstirred reactors. "This is why the analysis of such a system is so difficult." Historically, the complexities associated with natural convection were perhaps recognized as early as 1938 when thermal ignition theory was first developed. In the 1955 text "Diffusion and Heat Exchange in Chemical Kinetics", Frank-Kamenetskii recognized that "the purely conductive theory can be applied at sufficiently low pressure and small dimensions of the vessel when the influence of natural convection can be disregarded." This was reiterated by Tyler in 1966 and further emphasized by Barnard and Harwood in 1974. Specifically, they state: "It is generally assumed that heat losses are purely conductive. While this may be valid for certain low pressure slow combustion regimes, it is unlikely to be true for the cool flame and ignition regimes." While this statement is true for terrestrial experiments, the purely conductive heat transport assumption is valid at microgravity (mu-g). Specifically, buoyant complexities are suppressed at mu-g and the reaction-diffusion structure associated with low temperature oxidation reactions, cool flames and auto-ignitions can be studied. Without natural convection, the system is simpler, does not require determination of the effective heat transfer coefficient, and is a testbed for analytic and numerical models that assume pure diffusive transport. In addition, mu-g experiments will provide baseline data that will improve our understanding of the effects of natural convection on Earth.

Pearlman, Howard↗

NASA NTRS2019

Surface Roughness-Induced Spectral Degradation of Multi-Spaceborne Solar Diffusers Due to Space Radiation Exposure

Solar diffusers (SDs) have often been used as the onboard calibrators for the radiometric calibration of reflective solar band imaging sensors. After being spaceborne, the reflectance of SDs is observed to degrade with spectral dependence due to exposure to solar UV and energetic particle radiation. Long-term spectral reflectance data of SDs onboard multiple LEO imaging sensors, such as the Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua and the Visible Infrared Imaging Radiometer Suite (VIIRS) on SNPP, are analyzed. The reflectance of SDs on these three instruments degrades faster for the shorter wavelength (0.4-0.6 μm) bands than the longer wavelength bands. The Surface Roughness-induced Rayleigh Scattering (SRRS) model is applied to simulate the SD degradation on these instruments, and the growth of the surface roughness parameter of the SDs is derived. It is determined that the change of surface roughness scale length is ~tens of nanometers. To show the consistency of roughness growth rates among the SDs on Terra/Aqua MODIS and SNPP VIIRS instruments, the functional dependences of the growth rates are characterized according to the SD exposure time and the stage of surface roughness. It is also found that the flattening or reverse in the growth trend of the surface roughness for these three SDs occurred around the same interval between October 2013 and October 2015. The confirmation of the applicability of SRRS model with the long-term spectral reflectance data from three independent spaceborne SDs facilitates a better understanding of the origin and physical processes of the SD degradation.

Visible Infrared Imaging Radiometer Suite (VIIRS)↗

NASA NTRSDate not supplied

Thermodynamic and Diffusion Model Estimates on Metamorphic Temperatures and Timescales for Basaltic Eucrite GRA 98098

Introduction: HED meteorites are thought to rep-resent igneous rocks from Vesta’s basaltic crust and preserve evidence of early crustal metamorphism. Determining the temperatures and timescales of thermal metamorphism is important for reconstructing crustal evolution in the early solar system. Here, we study basaltic eucrite Graves Nunataks (GRA) 98098, which has been identified as a highly metamorphosed eucrite [1]. We present new estimates on metamorphic temperatures determined via thermodynamic modeling as well as the initial results from diffusion models constraining timescales of thermal metamorphism. Sample Description: GRA 98098 is an unbrecciated eucrite with a granoblastic plagioclase and pyroxene mineralogy. Millimeter to cm-long lathes of tridymite cross-cut and poikilitically enclose plagioclase and pyroxene [1,this work]. Pyroxene grains have exsolved into Ca-rich (~Wo38En29Fs33) and Ca-poor (~Wo4.5En36Fs59.5) lamellae. Both unzoned and zoned plagioclase grains are observed. Unzoned plagioclase grains are found solely with tridymite laths. These grains have ~An92 compositions. The cores of the zoned plagioclase grains have the same composition and thin, relatively sodic rims (~An67), (Fig. 1). The bulk sample is unusually enriched in highly in-compatible elements and has one of the most fractionated REE patterns reported [1]. Maximum metamorphic temperatures of 985±78°C have been estimated using two-pyroxene thermometry [2]. Methods: Thermodynamic modeling. Thermodynamic models were constructed using the software Perple_X, which employs a Gibbs free energy minimization in order to determine the most stable phase assemblage for a given bulk rock composition [3]. Bulk composition was calculated using mineral com-positions acquired via EMPA (this study) and the observed abundancies present in the thin section. Two bulk compositions were estimated; 1) includes all phases present in the thin section, (assumes that all phases are present during metamorphism), 2) excludes tridymite from the bulk calculation (assumes that tridymite was not present during metamorphism). In order to determine whether metamorphic equilibria was achieved and estimate temperatures of metamorphism, we compared measured pyroxene compositions with thermodynamically predicted compositions [4, Fig. 2]. Diffusion Modeling. Several time-temperature de-pendent diffusion profiles were calculated in order to determine the best match for XAn chemical profiles observed at the edges of the zoned plagioclase (Fig. 3). We assumed that the start condition was a stepwise gradient at the plagioclase/pyroxene interface. We also assumed an average diffusion coefficient (D) and a constant temperature using the equation in [5]. D was determined for two temperatures (T = 1060ºC; near eucrite solidus [6] and T = 985ºC; metamorphism reported in [4]) and then XAn was calculated as a function of distance from plagioclase core to rim using an error function solution to Fick’s second law. Results: Thermodynamic model results are summarized in Fig. 2. For a bulk composition that includes all phases in the thin section, pyroxene endmember compositions plot in the following temperature ranges: Fs ~660-860ºC, En~1000ºC & 1150ºC, and Wo~760-900ºC (Fig. 2a). For a bulk composition that excludes tridymite from the peak metamorphic assemblage (i.e., the bulk composition minus the contribution from tridymite), a temperature range could not be determined for the Fs component of pyroxene. For Wo, T~760-900ºC and En, T~1000ºC & 1150ºC (Fig. 2b). Fig. 3 summarizes the diffusion model results. For T = 1060°C & 985°C, the most appropriate time interval was estimated based on which diffusion curve most matched (solid lines, Fig. 3) the EMPA data. For T = 1060°C, the best looking match was t = 500 ka. For T = 985°C, the best match was t = 7 Ma. Discussion and future work: Temperature estimates from thermodynamic models are not conclusive because the temperature ranges determined for pyroxene endmember stability do not overlap (colored fields in Fig. 2), thus implying that there is disequilibrium between pyroxene crystals and the bulk composition considered [4]. Thus, additional exploration is needed to define a metamorphically equilibrated do-main that accurately records peak temperature. The utility of defining metamorphically equilibrated do-mains to improve the accuracy and level of detail elucidated regarding the petrogenetic history of metamorphose samples has been demonstrated previously [4,7]. We suggest that in the case of Fig. 2a, the thin section composition is not representative of the length scales over which metamorphic equilibrium was achieve and in the case of Fig. 2b, the assumption that tridymite was not present during metamorphism was incorrect. However, results from thermodynamic models can provide insight into the relative timing of mineral and compositional textures. For example from texture alone, it is unclear whether tridymite was igneous in origin and represents the last bits of melt in a crystallizing magma chamber, or if it formed during (and possibly initiated) open system thermal metamorphism. The latter could be consistent with a partial melt hypothesis [8,9] while the former implies that simple fractional crystallization can yield the textures present in GRA 98098. The lack of coincidence be-tween pyroxene endmember compositions in Fig. 2b suggest that the bulk composition minus tridymite was not the assemblage in equilibrium with the pyroxene, suggesting that tridymite was present during metamorphism and formed during igneous crystallization. We conclude that the development of the Na-rich plagioclase rims likely occurred during or immediately after peak thermal metamorphism, because eucrites of similar metamorphic grade and texture have unzoned plagioclase (~An92) [2,4,8], and Na zoning is only observed in the plagioclase not included in the tridymite. This suggests that the zoning formed after tridymite formation, and therefore after igneous crystallization. Thus, the timescales calculated via diffusion modeling possibly represent the time interval over which thermal metamorphism occurred. Cooling rates approximated for the Vestan crust predict that the crust cooled below 300°C around 35-40 Ma after formation[10]. This is consistent with our modeling results that predict formation of the Na rich plagioclase rims occurring at higher temperatures over a period of 0.5 to 7 Ma years. Future work. Additional thermodynamic modeling work will focus on selecting an equilibrated bulk rock domain in which to elucidate metamorphic conditions. Diffusion models currently provide a minimum time-scale, since diffusion slows down as the system cools. Future work with will focus on integrating cooling into the diffusion models and constraining the depth at which thermal metamorphism occurs because it could be used to determine whether the range of time-scales calculated for thermal metamorphism are consistent with the geologic environment.

J S Gorce↗

DOE OSTI2020

Connecting Hindered Transport in Porous Media across Length Scales: From Single-Pore to Macroscopic

Hindered mass transport is widely observed in various porous media, including filtration membranes and porous rocks. However, there is no universal model capable of predicting transport in porous media due to the heterogeneity of porous structures and the complexity of the underlying microscopic mechanisms, including electrostatic and hydrodynamic effects. Here, we used a highly ordered porous medium (i.e. a silica inverse opal film) as a model system to directly explore the effects of geometric parameters (i.e. pore size, pore throat size and tracer particle size) and microscopic interaction parameters (e.g. controlled by ionic strength) on nanoparticle transport in porous environments using three-dimensional single-particle tracking. We found a linear scaling relation between the macroscopic diffusion coefficient and microscopic diffusion behaviors and a combination of parameters associated with microscopic features and phenomena, including both geometric effects and particle-level interactions. The proportionality coefficient relating micro and macro behaviors was complex and related to the connectivity of the matrix and pore-size variation, which could lead to tortuous diffusion pathways, hindering macroscopic transport. Furthermore, these findings provide important insights into nanoparticles and macromolecules transport in various porous media.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2000

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗

DOE OSTI2021

Faster Form of Electron Magnetic Reconnection with a Finite Length X-Line

Observations in Earth’s turbulent magnetosheath downstream of a quasiparallel bow shock reveal a prevalence of electron-scale current sheets favorable for electron-only reconnection where ions are not coupled to the reconnecting magnetic fields. In small-scale turbulence, magnetic structures associated with intense current sheets are limited in all dimensions. And since the coupling of ions are constrained by a minimum length scale, the dynamics of electron reconnection is likely to be 3D. Here, both 2D and 3D kinetic particle-in-cell simulations are used to investigate electron-only reconnection, focusing on the reconnection rate and associated electron flows. A new form of 3D electron-only reconnection spontaneously develops where the magnetic X-line is localized in the out-of-plane (z) direction. The consequence is an enhancement of the reconnection rate compared with two dimensions, which results from differential mass flux out of the diffusion region along z, enabling a faster inflow velocity and thus a larger reconnection rate. This outflow along z is due to the magnetic tension force in z just as the conventional exhaust tension force, allowing particles to leave the diffusion region efficiently along z unlike the 2D configuration.

79 ASTRONOMY AND ASTROPHYSICS↗