Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional diffusion”

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 757 records · Page 42

Molecular aspect ratio effect on axial thermal transport in solution-spun carbon nanotube fibers

Neat, densely packed, and highly aligned carbon nanotube fibers (CNTFs) have appealing room-temperature axial thermal conductivity (⁠k⁠) and thermal diffusivity (⁠α⁠) for applications in lightweight heat spreading, flexible thermal connections, and thermoelectric active cooling. Although CNTFs are regularly produced from different input carbon nanotubes (CNTs), prior work has not quantified how the CNT molecular aspect ratio r (i.e., molecular length-to-diameter ratio) influences k and in well-aligned, packed CNTFs. Here, we perform self-heated steady-state and three-omega thermal measurements at room temperature on CNTF suspended in vacuum. Our results show that k increases from 150 to 380 W / mK for viscosity-averaged molecular aspect ratios increasing from r = 960 to 5600 and nanotube diameters of ~2 nm, which we attribute to the effects of thermal resistances between CNT bundles. CNTFs made with varying volume fraction Φ of constituent high-r and low-r CNT have properties that fall within or below the typical macroscopic rule-of-mixtures bounds. The thermal diffusivity α scales with k, leading to a sample-averaged volumetric heat capacity of 1.5 ± 0.3 MJ / m 3 K⁠. This work's findings that fibers made from longer CNT have larger k and α at room temperature motivate further investigation into thermal transport in solution-spun CNTF.

42 ENGINEERING↗

Diffusion in Melts and Magmas

Diffusion results from random motion of particles and entities. Diffusion in melts and magmas is due to thermally excited random motion of atoms, ions, and clusters, and plays a critical role in magmatic and volcanic processes. In melts and magmas, diffusion is one of the two mechanisms of mass transfer; the other being bulk flow (referred to as convection or advection). When both are present, diffusion refers to the dispersive motion relative to the mean bulk flow in a given reference frame (Richter et al. 1998). Diffusion plays critical roles in controlling magma mixing (Watson 1982; Koyaguchi 1985, 1989; Lesher 1994; Huber et al. 2009; Guo and Zhang 2020), mineral growth and dissolution rates in magmas (e.g., Watson 1982; Harrison and Watson 1983; Zhang et al. 1989; Newcombe et al. 2014; Macris et al. 2018), bubble growth and dissolution rate in magmas (Sparks 1978; Proussevitch and Sahagian 1998; Liu and Zhang 2000; Zhang 2013), and elemental and isotope fractionation during mineral 2 growth and dissolution (Jambon 1980; Richter et al. 1999, 2003; Watson and Muller 2009; Chopra et al. 2012; Watkins et al. 2014, 2017; Holycross et al. 2016, 2018). As a result, diffusion also plays an essential role in explosive volcanic eruptions and magma crystallization. Furthermore, diffusion has important applications in geospeedometry (Lasaga 1983, 1998; Zhang 1994, 2008; Trail et al. 2016; Zhang and Xu 2016). Experimental investigation of diffusion in geologically relevant silicate melts began to flourish in the 1970’s when micro-analytical measurements of diffusion profiles became available. In addition to the vast number of papers published since then, numerous books and reviews are available for diffusion in silicate melts. Hofmann et al. (1974) edited a book titled “Geochemical Transport and Kinetics” published by Carnegie Institution of Washington. This was the first landmark book summarizing the field. Lasaga and Kirkpatrick (1981) edited a book “Kinetics of Geochemical Processes” as volume 8 of the Reviews in Mineralogy (later becoming Reviews in Mineralogy and Geochemistry) series. Zhang and Cherniak (2010) edited “Diffusion in Minerals and Melts” as volume 72 of Reviews of Mineralogy and Geochemistry series, in which one chapter focused on diffusion theory, five chapters on diffusion in silicate melts (Behrens 2010; Lesher 2010; Liang 2010; Zhang and Ni 2010; Zhang et al. 2010), and other chapters were on experimental, analytical, and computational methods, and diffusion in minerals. Several textbooks covered the principles and applications of diffusion theories (Kirkaldy and Young 1987; Shewmon 1989; Cussler 1997; Lasaga 1998; Zhang 2008; Vrentas and Vrentas 2016), and two classic books covered the mathematics of diffusion (Carslaw and Jaeger 1959; Crank 1975). In preparing for this review chapter, we thought carefully about what to cover for this vast field, and decided to briefly go through the fundamentals of diffusion (more complete review can be found in Chakraborty 1995; Zhang 2008, 2010) and solutions to often-encountered 3 diffusion problems, and then focus on post-2010 diffusion studies on silicate melts and magmas. Here, melts refer to (mostly natural) silicate liquid, and magmas refer to crystal-bearing and/or bubble-bearing melts in which the continuous phase is the melt. There is a large body of work on diffusion in glasses, especially in the materials science literature, which is not covered in this review.

Youxue Zhang↗

Turbulent nonpremixed cool flames: Experimental measurements, Direct Numerical Simulation, and manifold-based combustion modeling

Turbulence, low-temperature chemistry, and their interactions in the form of turbulent cool flames are critical to understanding and improving advanced engines. Design of such engines requires tractable simulations which in turn necessitate turbulent combustion models that can account for cool flames. While manifold-based turbulent combustion models are an attractive option for hot flames, their applicability to cool flames is not yet fully understood. This is partially due to the lack of turbulent cool flame experiments, which has made model validation difficult. This work addresses these points with a combined experimental and computational investigation of turbulent nonpremixed dimethyl ether cool flames. First, a Co-flow Axisymmetric Reactor-Assisted Turbulent (CARAT) burner is developed and tested. Turbulent cool flames are studied using the formaldehyde planar laser-induced fluorescence (PLIF), acetone PLIF, and planar Rayleigh scattering techniques. The acetone PLIF signals are converted into mixture fraction values, and measurements of time-averaged temperature are derived from the Rayleigh scattering signals by taking advantage of the similarities between cool flames and unburned mixtures. Second, a Direct Numerical Simulation (DNS) of the CARAT burner is conducted. The flame is shown to be sensitive to thermal boundary conditions such that the stabilization method is dependent on a 10 K change in inlet temperature. Comparisons of first- and second-order statistics of temperature, mixture fraction, and formaldehyde between the DNS and experiment show good agreement. The validity of manifold-based turbulent combustion models in turbulent cool flames is then explored by first studying the effective Lewis numbers in the flame using a local differential diffusion parameter analysis. Turbulent cool flames are shown to have effective unity Lewis number transport even at a much lower Reynolds number than typically required for hot flames. Good agreement is shown for temperature and formaldehyde between DNS conditional means and one-dimensional nonpremixed flame solutions. Additionally, the motion of the cool flame in mixture fraction space is shown to be well described by one-dimensional nonpremixed flame solutions. These points indicate that the fundamental physics of turbulent cool flames are analogous to turbulent hot flames, implying that the same reduced-order manifold-based modeling approaches utilized for turbulent hot flames can be utilized for turbulent cool flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Macromolecule Sorption and Diffusion in HEMA/MAA Hydrogels

Transient solute absorption and desorption concentration profiles were measured in a 70 wt % hydroxyethyl methacrylate (HEMA)/30 wt % methacrylic acid (MAA) anionic hydrogel using two-photon confocal microscopy. Dilute aqueous solutes included fluorescently labeled dextrans with molecular masses of 4, 10, and 20 kDa, and fluorescently labeled cationic avidin protein. Cross-linking densities with ethylene glycol dimethacrylate (EGDMA) varied from 0 to 1 wt % with polymer volume fractions increasing from 0.15 to 0.25. Average gel mesh sizes, determined from zero-frequency oscillatory shear storage moduli, ranged from about 3.6 to 8.4 nm over the cross-link ratios studied. All solutes exhibit Stokes-Einstein hydrodynamic radii obtained from measured free diffusion coefficients, D o , comparable to or larger than the average gel mesh size. In spite of considerable size exclusion, the studied solutes penetrate the gels indicating a range of mesh sizes available for transport. Transient uptake and release concentration profiles for FITC-dextrans follow simple diffusion theory with diffusion coefficients, D , essentially independent of loading or release characteristic of reversible absorption. Although strongly size-excluded, these solutes do not interact specifically with the polymer network. Diffusivities are accordingly predicted from a large-pore effective-medium (LPEM) model developed to account for solute size, hydrodynamic drag, and distribution of mesh sizes available for transport in the polymer network. For this class of solute, anrs of magnitude slower than that of dextran of similar hydrodynamic radius. Desorption of avidin is not cd using no adjustable parameters, diffusivities predicted from the new effective-medium model demonstrate good agreement with experiment. For the specific-interacting cationic protein, avidin, gel loading is 3 ordeomplete even after 2 weeks of extraction. On the basis of size alone, avidin is strongly size-excluded, yet it exhibits a partition coefficient of over 20. For the positively charged protein, we observed specific ion binding on the negatively charged carboxylate groups of MAA-decorated polymer strands in the larger mesh spaces. Finally, simple linear sorption kinetics gives an adsorption time constant of 5 min and a desorption time constant of about 20 days, suggesting nearly irreversible uptake of cationic avidin on the anionic gel matrix.

diffusion↗

What Can Plastic Flow Fields Tell Us About Heat Sources in Deformation Processing?

We characterize primary (shape-change) and secondary (friction) deformation, and associated temperature fields, in metal cutting and forming processes, using in situ imaging and simulation. The experimental configurations enable access to the deformation zones and die contact-interfaces, for measuring deformation, temperature, and frictional drag. Infra-red thermography reveals that the plastic strain-rate field is an excellent proxy for the deformation-induced heat sources. Both spatially confined and diffuse strain-rate fields occur, depending on the initial workpiece deformation state. When the strain rate is confined, as in pre-hardened material, the temperature modeling is much simplified, as the heat source is also now localized. However complex, microstructure-driven deformation modes, like sinuous flow in annealed metals, result in spatially diffuse strain-rate and body heat sources, more challenging to analyze. Furthermore, our unified measurements should be of value for accurately estimating the fraction of plastic dissipation that is converted into heat in large-strain deformation processes.

36 MATERIALS SCIENCE↗

Use of physics to improve solar forecast: Part III, impacts of different cloud types

Cloud-type impacts present a great challenge to solar forecasting due to diverse and complex cloud-radiation interactions. This third part of our paper sequence seeks to address this challenge by quantifying the forecast accuracies under eight cloud types: cumulus (Cu), stratified clouds (St), altocumulus (Ac), altostratus (As), cirrostratus/anvil (Cr), cirrus (Ci), congestus (Co), deep convective clouds (Dc) across four physics-informed persistence models reported in Part I. To generalize the cloud impacts, the eight cloud types are further grouped into three cloud categories based on their common features: weak convective clouds, stratiform clouds, and strong convective clouds. Here, the decade-long (2001 ~ 2014) collocated measurements of irradiances and cloud types at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Program South Great Plain (SGP) Central Facility site are used for model evaluation. Results reveal a clear performance hierarchy for global horizontal irradiance (GHI) and direct normal irradiance (DNI): best for weak convective clouds and cirrus, intermediate for stratiform clouds, and worst for strong convective clouds. Performance for diffuse horizontal irradiance (DHI) is less influenced by cloud types. Cloud albedo dominates all three irradiances for Dc, while both cloud albedo and cloud fraction are influential for other cloud types. A 12 %~33 % improvement in accuracy at 6-hour lead time compared to the benchmark smart model confirms the effectiveness of incorporating physics into the models for various cloud types; further improvements are expected by directly integrating cloud type information into forecasting models by modifying the physical formulation of cloud-radiation interaction, and/or using more advanced machine learning models.

14 SOLAR ENERGY↗

Ab initio prediction of fast non-equilibrium transport of nascent polarons in SrI 2 : a key to high-performance scintillation [First-principles study of hole polaron formation and migration in strontium iodide]

The excellent light yield proportionality of europium-doped strontium iodide (SrI 2 :Eu) has resulted in state-of-the-art γ-ray detectors with remarkably high-energy resolution, far exceeding that of most halide compounds. In this class of materials, the formation of self-trapped hole polarons is very common. However, polaron formation is usually expected to limit carrier mobilities and has been associated with poor scintillator light-yield proportionality and resolution. Here using a recently developed first-principles method, we perform an unprecedented study of polaron transport in SrI 2 , both for equilibrium polarons, as well as nascent polarons immediately following a self-trapping event. We propose a rationale for the unexpected high-energy resolution of SrI 2 . We identify nine stable hole polaron configurations, which consist of dimerised iodine pairs with polaron-binding energies of up to 0.5 eV. They are connected by a complex potential energy landscape that comprises 66 unique nearest-neighbour migration paths. Ab initio molecular dynamics simulations reveal that a large fraction of polarons is born into configurations that migrate practically barrier free at room temperature. Consequently, carriers created during γ-irradiation can quickly diffuse away reducing the chance for nonlinear recombination, the primary culprit for non-proportionality and resolution reduction. We conclude that the flat, albeit complex, landscape for polaron migration in SrI 2 is a key for understanding its outstanding performance. This insight provides important guidance not only for the future development of high-performance scintillators but also of other materials, for which large polaron mobilities are crucial such as batteries and solid-state ionic conductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vegetation clumping modulates global photosynthesis through adjusting canopy light environment

Abstract The spatial dispersion of photoelements within a vegetation canopy, quantified by the clumping index (CI), directly regulates the within‐canopy light environment and photosynthesis rate, but is not commonly implemented in terrestrial biosphere models to estimate the ecosystem carbon cycle. A few global CI products have been developed recently with remote sensing measurements, making it possible to examine the global impacts of CI. This study deployed CI in the radiative transfer scheme of the Community Land Model version 5 (CLM5) and used the revised CLM5 to quantitatively evaluate the extent to which CI can affect canopy absorbed radiation and gross primary production (GPP), and for the first time, considering the uncertainty and seasonal variation of CI with multiple remote sensing products. Compared to the results without considering the CI impact, the revised CLM5 estimated that sunlit canopy absorbed up to 9%–15% and 23%–34% less direct and diffuse radiation, respectively, while shaded canopy absorbed 3%–18% more diffuse radiation across different biome types. The CI impacts on canopy light conditions included changes in canopy light absorption, and sunlit–shaded leaf area fraction related to nitrogen distribution and thus the maximum rate of Rubisco carboxylase activity ( V cmax ), which together decreased photosynthesis in sunlit canopy by 5.9–7.2 PgC year −1 while enhanced photosynthesis by 6.9–8.2 PgC year −1 in shaded canopy. With higher light use efficiency of shaded leaves, shaded canopy increased photosynthesis compensated and exceeded the lost photosynthesis in sunlit canopy, resulting in 1.0 ± 0.12 PgC year −1 net increase in GPP. The uncertainty of GPP due to the different input CI datasets was much larger than that caused by CI seasonal variations, and was up to 50% of the magnitude of GPP interannual variations in the tropical regions. This study highlights the necessity of considering the impacts of CI and its uncertainty in terrestrial biosphere models.

59 BASIC BIOLOGICAL SCIENCES↗

Pathway to tune sputter coated Nb$_3$Sn with Zr

Superconducting radio frequency (SRF) cavities are essential elements of many modern-day particle accelerators such as, European XFEL (DESY), LCLS II-HE (SLAC) and PIP II (Fermilab). Nb$_3$Sn superconducting material promises significant potential to exceed the performance of conventional niobium based SRF cavities. We aim to advance the ongoing R$&$D efforts to achieve Nb$_3$Sn coated cavities with high material quality. In this work, we have developed to achieve a route for effectively reducing the intrinsic defects in magnetron sputter coated Nb$_3$Sn SRF cavities and eventually improve the overall RF performance. We have observed that Zr is mainly being incorporated as ZrO$_2$ precipitates of average dimensions ranging from 20-100nm. By increasing the concentration of Zr, we have three major effects in the physical properties. Firstly, we noted that the density of the surface and bulk voids as well as their average dimensions are dramatically reduced on increasing the Zr content in Nb$_3$Sn. These surface voids are believed to be detrimental to cavity applications due to increase the local dissipation of electric and magnetic fields. These results indicate that inclusion of Zr in Nb$_3$Sn sputtered coating will improve the material quality and might help to reduce the overall RF power dissipation. Secondly, we observed that increasing Zr concentration from 0 to 2.4% substantially improves both superconducting transition temperature (Tc) and upper critical magnetic field (Hc2). However, on further increasing the Zr concentration to 24%, we noticed a small decrement in the Hc2 value. These results infer that the increase in Zr concentration up to an optimal concentration will be beneficial to enhance the superconducting properties. Increase in Tc promises lesser RF losses and improvement in accelerating gradient, whereas the increase in Hc2 suggests more efficient flux pinning which could give rise to high quality factors in presence of multi-Tesla magnetic field. Additionally, we also observed that increasing Zr concentration possibly prevents the oxygen diffusion and leads to thinner formation of primary surface oxide (Nb$_2$O$_5$ and SnO). This study gives a route to tune Zr fraction in Nb3Sn to modify the physical properties that control the high field performance of Nb3Sn RF cavities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

KINETIC SIMULATIONS OF THE LOWEST-ORDER UNSTABLE MODE OF RELATIVISTIC MAGNETOSTATIC EQUILIBRIA

We present the results of particle-in-cell numerical pair plasma simulations of relativistic two-dimensional magnetostatic equilibria known as the “Arnold–Beltrami–Childress” fields. In particular, we focus on the lowest-order unstable configuration consisting of two minima and two maxima of the magnetic vector potential. Breaking of the initial symmetry leads to exponential growth of the electric energy and to the formation of two current layers, which is consistent with the picture of “X-point collapse” first described by Syrovatskii. Magnetic reconnection within the layers heats a fraction of particles to very high energies. After the saturation of the linear instability, the current layers are disrupted and the system evolves chaotically, diffusing the particle energies in a stochastic second-order Fermi process, leading to the formation of power-law energy distributions. The power-law slopes harden with the increasing mean magnetization, but they are significantly softer than those produced in simulations initiated from Harris-type layers. The maximum particle energy is proportional to the mean magnetization, which is attributed partly to the increase of the effective electric field and partly to the increase of the acceleration timescale. We describe in detail the evolving structure of the dynamical current layers and report on the conservation of magnetic helicity. These results can be applied to highly magnetized astrophysical environments, where ideal plasma instabilities trigger rapid magnetic dissipation with efficient particle acceleration and flares of high-energy radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

After the interaction: an efficiently star-forming molecular disk in NGC 5195

We present new molecular gas maps of NGC 5195 (alternatively known as M51b) from the Combined Array for Research in Millimeter Astronomy, including {sup 12}CO(1–0), {sup 13}CO(1–0), CN(1{sub 0,2}–0{sub 0,1}), CS(2–1), and 3 mm continuum. We also detected HCN(1–0) and HCO{sup +}(1–0) using the Onsala Space Observatory. NGC 5195 has a {sup 12}CO/{sup 13}CO ratio (R{sub 12/13} = 11.4 ± 0.5) consistent with normal star-forming galaxies. The CN(1–0) intensity is higher than is seen in an average star-forming galaxy, possibly enhanced in the diffuse gas in photo-dissociation regions. Stellar template fitting of the nuclear spectrum of NGC 5195 shows two stellar populations: an 80% mass fraction of old (≳10 Gyr) and a 20% mass fraction of intermediate-aged (≈1 Gyr) stellar populations. This provides a constraint on the timescale over which NGC 5195 experienced enhanced star formation during its interaction with M51a. The average molecular gas depletion timescale in NGC 5195 is 〈τ{sub dep}〉 = 3.08 Gyr, a factor of ≈2 larger than the depletion timescales in nearby star-forming galaxies, but consistent with the depletion seen in CO-detected early-type galaxies. While radio continuum emission at centimeter and millimeter wavelengths is present in the vicinity of the nucleus of NGC 5195, we find it is most likely associated with nuclear star formation rather than radio-loud AGN activity. Thus, despite having a substantial interaction with M51a ∼1/2 Gyr ago, the molecular gas in NGC 5195 has resettled and is currently forming stars at an efficiency consistent with settled early-type galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

ALMA reveals weak [N ii] emission in “typical” galaxies and intense starbursts at z = 5-6

We report interferometric measurements of [N ii] 205 μm fine-structure line emission from a representative sample of three galaxies at z = 5–6 using the Atacama Large (sub)Millimeter Array (ALMA). These galaxies were previously detected in [C ii] and far-infrared continuum emission and span almost two orders of magnitude in star formation rate (SFR). Our results show at least two different regimes of ionized interstellar medium properties for galaxies in the first billion years of cosmic time, separated by their L{sub [CII]}/L{sub [NII]} ratio. We find extremely low [N ii] emission compared to [C ii] (L{sub [CII]}/L{sub [NII]}=68{sub −28}{sup +200}) from a “typical” ∼L{sub UV}{sup ∗} star-forming galaxy, likely directly or indirectly (by its effect on the radiation field) related to low dust abundance and low metallicity. The infrared-luminous modestly star-forming Lyman-break galaxy (LBG) in our sample is characterized by an ionized-gas fraction (L{sub [CII]}/L{sub [NII]}≲20) typical of local star-forming galaxies and shows evidence for spatial variations in its ionized-gas fraction across an extended gas reservoir. The extreme SFR, warm and compact dusty starburst AzTEC-3 shows an ionized fraction higher than expected given its SFR surface density (L{sub [CII]}/L{sub [NII]}=22±8) suggesting that [N ii] dominantly traces a diffuse ionized medium rather than star-forming H ii regions in this type of galaxy. This highest redshift sample of [N ii] detections provides some of the first constraints on ionized and neutral gas modeling attempts and on the structure of the interstellar medium at z = 5–6 in “normal” galaxies and starbursts.

79 ASTRONOMY AND ASTROPHYSICS↗

A simple reaction-rate model for turbulent diffusion flames

A simple reaction rate model is proposed for turbulent diffusion flames in which the reaction rate is proportional to the turbulence mixing rate. The reaction rate is also dependent on the mean mass fraction and the mean square fluctuation of mass fraction of each reactant. Calculations are compared with experimental data and are generally successful in predicting the measured quantities.

Bangert, L. H.↗

Pulsar and diffuse contributions to the observed galactic gamma radiation

With the acquisition of satellite data on the energy spectrum of galactic gamma-radiation, it is clear that such radiation has a multicomponent nature. A calculation of the pulsar gamma ray emission spectrum is used together with a statistical analysis of recent data on 328 known pulsars to make a new determination of the pulsar contribution to galactic gamma ray emission. The contributions from diffuse interstellar cosmic ray induced production mechanisms to the total emission are then reexamined. It is concluded that pulsars may account for a significant fraction of galactic gamma ray emission.

Harding, A. K.↗

Analysis of Pt/SnO(sub x) during catalysis of CO oxidation

Temperature-programmed reduction using 6kPaH2 suggests that a sample consisting of 3 percent Pt supported directly on SnO2 is, under conditions of catalysis of CO oxidation used here, best represented as 3 percent Pt/SnO sub x, since the support is likely to partially reduced, probably in the vicinity of the metal/oxide interface. Catalytic measurements at 421 to 424 K show that this 3 percent Pt/SnO sub x is significantly more active per unit area of Pt than 6 percent Pt/SiO2 in catalyzing the oxidation of CO. In situ micro-FTIR reveals that while the latter has predominantly linearly bound CO on the surface under reaction conditions, the Pt/SnO sub x also has a species absorbing at 2168 cm(exp -1) which may be CO upon Pt in a positive oxidation state or weakly chemisorbed CO on zero-valent Pt. This may be directly involved in the low temperature oxidation of CO on the Pt/SnO sub x, since being weakly held the activation energy for its surface diffusion to the metal/oxide interface will be low; such mobile species could allow the high rates of surface transport and an increase in the fraction of the surface over which the CO oxidation occurs. FTIR also reveals carbonate-type species on the P/SnO sub c surface.

Sermon, Paul A.↗

Spinel-bearing, Al-rich chondrules in two chondrite finds from Roosevelt County, New Mexico - Indicators of nebular and parent body processes

Two rare spinel-bearing Al-rich chondrules are identified in chondrite finds from Roosevelt County, New Mexico-RC 071 (L4) and RC 072 (L5). These chondrules have unusual mineralogies dominated by highly and asymmetrically zoned Al-Cr-rich spinels. Two alternatives exist to explain the origin of this zoning-fractional crystallization or metamorphism. Fractional crystallization formed the zoning of the trivalent cations and caused a localized depletion in chromites around the large Al-Cr-rich spinels. Diffusive exchange and partitioning of Fe and Mg between olivine and spinel during parent-body metamorphism can explain the asymmetric zoning of these elements. The bulk compositions of the chondrules suggest affinities with the Na-Al-Cr-rich chondrules, as would be expected from the abundance of Al-Cr-rich spinels. The most important factors are the temperature to which the molten chondrule was heated and the cooling rate during crystallization. These two chondrules cooled rapidly from near the liquidus, as indicated by the zoning, occurrence and sizes of spinels, radiating chondrule textures and localized chromite depletions.

Mccoy, Timothy J.↗

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.↗

GMC Collisions as Triggers of Star Formation. IV. The Role of Ambipolar Diffusion

We investigate the role of ambipolar diffusion (AD) in collisions between magnetized giant molecular clouds (GMCs), which may be an important mechanism for triggering star cluster formation. Three-dimensional simulations of GMC collisions are performed using a version of the Enzo magnetohydrodynamics code that has been extended to include AD. The resistivities are calculated using the 31-species chemical model of Wu et al. (2015). We find that in the weak-field, 10 μG case, AD has only a modest effect on the dynamical evolution during the collision. However, for the stronger-field, 30 μG case involving near-critical clouds, AD results in the formation of dense cores in regions where collapse is otherwise inhibited. The overall efficiency of formation of cores with n{sub H}⩾10{sup 6} cm{sup −3} in these simulations is increases from about 0.2% to 2% once AD is included, comparable to observed values in star-forming GMCs. The gas around these cores typically has relatively slow infall at speeds that are a modest fraction of the free-fall speed.

79 ASTRONOMY AND ASTROPHYSICS↗