Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Instabilities in uranium plasma and the gas-core nuclear rocket engine

The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P approximates 0.00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.

Tidman, D. A.↗
NASA NTRS1982

A numerical solution of the supersonic flow over a rearward facing step with transverse non-reacting hydrogen injection

The present investigation represents an application of computational fluid dynamics to a problem associated with the flow in the combustor region of a supersonic combustion ramjet engine (scramjet). The governing equations are considered, taking into account the Navier-Stokes equations, a molecular viscosity calculation, the molecular thermal conductivity, molecular diffusion, and a turbulence model. The employed numerical solution is patterned after the explicit, time-dependent, unsplit, predictor-corrector, finite-difference method given by MacCormack (1969). The calculation is concerned with the supersonic flow over a rearward-facing step with transverse H2 injection at conditions germane to the combustor region of a scramjet engine. The H2 jet acts as an effective body which essentially shields the primary flow from the rearward-facing step, thus substantially changing the wave pattern in the primary flow.

Berman, H. A.↗
NASA NTRS1994

The growth of zeolites A, X and mordenite in space

Zeolites are a class of crystalline aluminosilicate materials that form the backbone of the chemical process industry worldwide. They are used primarily as adsorbents and catalysts and support to a significant extent the positive balance of trade realized by the chemical industry in the United States (around $19 billion in 1991). The magnitude of their efforts can be appreciated when one realizes that since their introduction as 'cracking catalysts' in the early 1960's, they have saved the equivalent of 60 percent of the total oil production from Alaska's North Slope. Thus the performance of zeolite catalysts can have a profound effect on the U.S. economy. It is estimated that a 1 percent increase in yield of the gasoline fraction per barrel of oil would represent a savings of 22 million barrels of crude oil per year, representing a reduction of $400 million in the United States' balance of payments. Thus any activity that results in improvement in zeolite catalyst performance is of significant scientific and industrial interest. In addition, due to their 'stability,' uniformity, and, within limits, their 'engineerable' structures, zeolites are being tested as potential adsorbents to purify gases and liquids at the parts-per-billion levels needed in today's electronic, biomedical, and biotechnology industries and for the environment. Other exotic applications, such as host materials for quantum-confined semiconductor atomic arrays, are also being investigated. Because of the importance of this class of material, extensive efforts have been made to characterize their structures and to understand their nucleation and growth mechanisms, so as to be able to custom-make zeolites for a desired application. To date, both the nucleation mechanics and chemistry (such as what are the 'key' nutrients) are, as yet, still unknown for many, if not all, systems. The problem is compounded because there is usually a 'gel' phase present that is assumed to control the degree of supersaturation, and this gel undergoes a continuous 'polymerization' type reaction during nucleation and growth. Generally, for structure characterization and diffusion studies, which are useful in evaluating zeolites for improving yield in petroleum refining as well as for many of the proposed new applications (e.g., catalytic membranes, molecular electronics, chemical sensors) large zeolites (greater than 100 to 1000 times normal size) with minimum lattice defects are desired. Presently, the lack of understanding of zeolite nucleation and growth precludes the custom design of zeolites for these or other uses. It was hypothesized that the microgravity levels achieved in an orbiting spacecraft could help to isolate the possible effects of natural convection (which affects defect formation) and minimize sedimentation, which occurs since zeolites are twice as dense as the solution from which they are formed. This was expected to promote larger crystals by allowing growing crystals a longer residence time in a high-concentration nutrient field. Thus it was hypothesized that the microgravity environment of Earth orbit would allow the growth of large, more defect-free zeolite crystals in high yield.

Sacco, Albert, Jr.↗
NASA NTRS2004

Parameter Estimation in Atmospheric Data Sets

In this study the structure tensor technique is used to estimate dynamical parameters in atmospheric data sets. The structure tensor is a common tool for estimating motion in image sequences. This technique can be extended to estimate other dynamical parameters such as diffusion constants or exponential decay rates. A general mathematical framework was developed for the direct estimation of the physical parameters that govern the underlying processes from image sequences. This estimation technique can be adapted to the specific physical problem under investigation, so it can be used in a variety of applications in trace gas, aerosol, and cloud remote sensing. As a test scenario this technique will be applied to modeled dust data. In this case vertically integrated dust concentrations were used to derive wind information. Those results can be compared to the wind vector fields which served as input to the model. Based on this analysis, a method to compute atmospheric data parameter fields will be presented. .

Wenig, Mark↗
NASA NTRS2004

Estimation of Dynamical Parameters in Atmospheric Data Sets

In this study a new technique is used to derive dynamical parameters out of atmospheric data sets. This technique, called the structure tensor technique, can be used to estimate dynamical parameters such as motion, source strengths, diffusion constants or exponential decay rates. A general mathematical framework was developed for the direct estimation of the physical parameters that govern the underlying processes from image sequences. This estimation technique can be adapted to the specific physical problem under investigation, so it can be used in a variety of applications in trace gas, aerosol, and cloud remote sensing. The fundamental algorithm will be extended to the analysis of multi- channel (e.g. multi trace gas) image sequences and to provide solutions to the extended aperture problem. In this study sensitivity studies have been performed to determine the usability of this technique for data sets with different resolution in time and space and different dimensions.

Wenig, Mark O.↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

A finite-volume HLLC-based scheme for compressible interfacial flows with surface tension

Shock waves are often used in experiments to create a shear flow across liquid droplets to study secondary atomization. Similar behavior occurs inside of supersonic combustors (scramjets) under startup conditions, but it is challenging to study these conditions experimentally. In order to investigate this phenomenon further, a numerical approach is developed to simulate compressible multiphase flows under the effects of surface tension forces. The flow field is solved via the compressible multicomponent Euler equations (i.e., the five equation model) discretized with the finite volume method on a uniform Cartesian grid. The solver utilizes a total variation diminishing (TVD) third-order Runge–Kutta method for time-marching and second order TVD spatial reconstruction. Surface tension is incorporated using the Continuum Surface Force (CSF) model. Fluxes are upwinded with a modified Harten–Lax–van Leer Contact (HLLC) approximate Riemann solver. An interface compression scheme is employed to counter numerical diffusion of the interface. The present work includes modifications to both the HLLC solver and the interface compression scheme to account for capillary force terms and the associated pressure jump across the gas–liquid interface. A simple method for numerically computing the interface curvature is developed and an acoustic scaling of the surface tension coefficient is proposed for the non-dimensionalization of the model. The model captures the surface tension induced pressure jump exactly if the exact curvature is known and is further verified with an oscillating elliptical droplet and Mach 1.47 and 3 shock-droplet interaction problems. The general characteristics of secondary atomization at a range of Weber numbers are also captured in a series of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2011

Magnetic Reconnection: Theoretical and Observational Perspectives: Preface

Magnetic reconnection is a fundamental plasma-physical process by which energy stored in a magnetic field is converted, often explosively, into heat and the kinetic energy of the charged particles that constitute the plasma. It occurs in a variety of astrophysical settings, ranging from the solar corona to pulsar magnetospheres and winds, as well as in laboratory fusion experiments, where it is responsible for sawtooth crashes. First proposed by R.G. Giovanelli in the late I 940s as the mechanism responsible for solar flares, magnetic reconnection was invoked at the beginning of the space age to explain not just solar flares but also the transfer of energy, mass, and momentum from the solar wind to Earth's magnetosphere and the subsequent storage and release of the transferred energy in the magnetotai\. During the half century or so that has followed the seminal theoretical works by J.W. Dungey, P.A. Sweet, E.N. Parker, and H.E. Petschek, in-situ measurements by Earth-orbiting satellites and remote-sensing observations of the solar corona have provided a growing body of evidence for the occurrence of reconnection at the Sun, in the solar wind, and in the near-Earth space environment. The last thirty years have also seen the development of laboratory reconnection experiments at a number of institutions. In parallel with the efforts of experimentalists in both space and laboratory plasma physics, theorists have investigated, analytically and with the help of increasingly powerful MHD, hybrid, and kinetic numerical simulations, the structure of the diffusion region, the factors controlling the rate, onset, and cessation of reconnection, and the detailed physics that enables the demagnetization of the ions and electrons and the topological reconfiguration of the magnetic field. Moreover, the scope of theoretical reconnection studies has been extended well beyond solar system and laboratory plasmas to include more exotic astrophysical plasma systems whose strong (10(exp 14)-10(exp 15) G) magnetic fields require that models of reconnection in these systems incorporate quantum electrodynamical, special relativistic, and radiative effects. The papers collected in this topical issue of Space Science Reviews cover different aspects of recent theoretical and observational work on magnetic reconnection in solar and space physics, astrophysics, and laboratory plasma physics. They derive from presentations given at a workshop on magnetic reconnection held in the Yosemite National Park, February 8-12,2010. The intent of the workshop was to stimulate, through a combination of tutorial talks, shorter focused talks, and extensive informal discussions, an interdisciplinary dialogue among members of the different research communities working on the problem of magnetic reconnection. One of the motivating considerations for holding the workshop was its relevance to NASA's Magnetospheric Multiscale (MMS) mission, scheduled for launch in 2014. The four identically instrumented MMS spacecraft are designed to study reconnect ion in Earth's magnetosphere and, specifically, to probe the electron diffusion region in order to determine the microphysical processes that enable the change in the topology of the magnetic field. Building on the achievements of the multi spacecraft Cluster and THEMIS missions, MMS will use the magnetosphere as an astrophysical plasma laboratory in which to test, through in-situ measurement of the plasma, energetic particles, and electric and magnetic fields, various models and theories that have emerged during the past twenty years, a period of extraordinarily productive theoretical and observational work.

Lewis, W. S.↗
DOE OSTI2015

TiO 2 Microboxes with Controlled Internal Porosity for High‐Performance Lithium Storage

Abstract Titanium dioxide (TiO 2 ) is considered a promising anode material for high‐power lithium ion batteries (LIBs) because of its low cost, high thermal/chemical stability, and good safety performance without solid electrolyte interface formation. However, the poor electronic conductivity and low lithium ion diffusivity of TiO 2 result in poor cyclability and lithium ion depletion at high current rates, which hinder them from practical applications. Herein we demonstrate that hierarchically structured TiO 2 microboxes with controlled internal porosity can address the aforementioned problems for high‐power, long‐life LIB anodes. A self‐templating method for the synthesis of mesoporous microboxes was developed through Na 2 EDTA‐assisted ion exchange of CaTiO 3 microcubes. The resulting TiO 2 nanorods were organized into microboxes that resemble the microcube precursors. This nanostructured TiO 2 material has superior lithium storage properties with a capacity of 187 mAh g −1 after 300 cycles at 1 C and good rate capabilities up to 20 C.

Gao, Xuehui↗
DOE OSTI2016

Actuator selection and placement for localized feedback flow control

The selection and placement of actuators and sensors to control compressible viscous flows is addressed by developing a novel methodology based upon the eigensystem structural sensitivity of the linearized evolution operator appropriate for linear feedback control. Forward and adjoint global modes are used to construct a space of possible perturbations to the linearized operator, which yields a small optimization problem for selecting the parameters that best achieve the control objective, including where they should be placed. The method is demonstrated by informing actuation to suppress amplification of the instabilities in boundary layer separation in a high-subsonic diffuser. Complete stabilization is observed in the separated shear layer for short downstream distances at modest Reynolds number. Higher Reynolds numbers and longer distances are expected to be more challenging to stabilize; here the control informed by the procedure still substantively suppresses amplification of instabilities. It is also demonstrated that more complex actuator–sensor selections may not yield superior controllers.

Mechanics↗
DOE OSTI2020

Fundamental gap of fluorographene by many-body GW and fixed-node diffusion Monte Carlo methods

Fluorographene (FG) is a promising graphene-derived material with a large bandgap. Currently existing predictions of its fundamental gap (Δ f ) and optical gap (Δ opt ) significantly vary when compared with experiment. We provide here an ultimate benchmark of Δ f for FG by many-body GW and fixed-node diffusion Monte Carlo (FNDMC) methods. Both approaches independently arrive at Δf ≈ 7.1 ± 0.1 eV. In addition, the Bethe–Salpeter equation enabled us to determine the first exciton binding energy, Eb = 1.92 eV. We also point to the possible misinterpretation problem of the results obtained for gaps of solids by FNDMC with single-reference trial wave functions of Bloch orbitals. We argue why instead of Δ opt , in the thermodynamic limit,

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

A modified many-body dissipative particle dynamics model for mesoscopic fluid simulation: methodology, calibration, and application for hydrocarbon and water

The many-body dissipative particle dynamics (mDPD) is a prominent mesoscopic multiphase model for fluid transport in mesoconfinement. However, it has been a long-standing challenge for mDPD (and other multiphase-enabled DPD models) to accurately predict real-fluid static and dynamic properties simultaneously. We have developed a modified mDPD model to overcome the issue and a rigorous calibration approach that uses reference data, including experimental and/or molecular dynamics (MD) simulations to parameterise the modified mDPD for real fluids. We choose heptane as a representative example of hydrocarbon in source rocks to demonstrate the model's capability to accurately predict the equation of state (EOS), free surface tension, diffusivity, and viscosity. Our timing test shows that the modified mDPD is 400–500 times faster than its MD counterpart for simulating bulk heptane in equivalent volumes. To further demonstrate the robustness of the model, we revisited the benchmark problem of mesoscopic modelling of water, in which all the previous works on DPD reported only a limited portion of the water properties. Here, we show that the modified mDPD can provide accurate modelling of water static and dynamic properties and an EOS that matches the experimental data to a large range of confinement pressure.

42 ENGINEERING↗
DOE OSTI2022

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Properties of Medium Nuclei from First Principles

During this project we have extended the Auxiliary Field Diffusion Monte Carlo (AFDMC) to calculate several properties of light/medium nuclei. In particular, we have included and tested several models of three-body forces and the associated approximations required, and the comparison with other exact calculations for simpler problems demonstrated the accuracy of the AFDMC method. In this project we have also developed new subroutines to solve for the ground-state of open shell nuclei. This case requires very complex wave function with respect to the case of closed-shell nuclei. We have used the AFDMC code to calculate many properties of nuclei up to the Oxygen, including energies, radii, spatial and momentum distributions and others.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Collaborative Research: Multi-scale Dynamics of Kinetic Turbulence in Weakly Collisional, High-Beta Plasmas

The transport of energy and momentum and the heating of plasma particles by waves and turbulence are key ingredients in many problems at the frontiers of heliospheric and astrophysics research. This includes the heating and acceleration of the solar wind; the observational appearance of black-hole accretion flows on event-horizon scales; and the properties of the hot, diffuse plasmas that fill dark-matter halos. This work addressed the physics of plasma waves and turbulence in weakly collisional, weakly magnetized (high-beta) plasmas. In this regime, deviations from local thermodynamic equilibrium (i.e., pressure anisotropies) and the instabilities they excite can dramatically change the properties of waves and turbulence from those predicted by fluid (i.e., magnetohydrodynamic) and linear kinetic theories. Through a combination of analytic theory and first-principles kinetic simulations, we have shown that: (1) ion-acoustic waves of sufficient amplitude drive pressure anisotropies large enough to trigger kinetic Larmor-scale instabilities, which in turn interrupt the otherwise strong collisionless damping such waves suffer; (2) pressure anisotropies generated by fluctuating magnetic fields in strong Alfvénic turbulence produce a viscous scale comparable to the driving scale and a cascade with significant non-local interactions; (3) strong Alfvénic turbulence in an expanding high-beta plasma self-organizes to maintain critical balance, even as kinetic Larmor-scale instabilities regulate the plasma viscosity and reduce the elasticity of magnetic-field lines; (4) imbalanced Alfvénic turbulence leads to a "helicity barrier" that ultimately triggers the cyclotron heating of ions, in a manner consistent with spacecraft measurements taken in the solar wind; and (5) relativistic synchrotron-radiating plasmas develop kinetically unstable pressure anisotropy, which changes the nature of their cooling in a way that can couple the thermodynamics of ions and electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Thermal preparation of foods in space-vehicle environments

Convection is the primary heat transfer mechanism for most foods heated in an earth-based environment. In contrast, in the low-gravity environment of space flight, the primary heat transfer mechanism is conduction (or radiation in the absence of a conducting medium). Conduction heating is significantly slower and less efficient than convection heating. This fact poses a problem for food heating during space flight. A numerical model has been developed to evaluate this problem. This model simulates the food-heating process for Skylab. The model includes the effect of a thermally controlled on/off heat flux. Parametric studies using this model establish how the required heating time is affected by: the thermal diffusivity of the nutrient materials, the power level of the heater, the initial food temperatures, and the food container dimensions.

Bannerot, R. B.↗
NASA NTRS1979

Modeling of premixing-prevaporizing fuel-air mixing passages

The development of a computer program for the analytical prediction of the distribution of liquid and vapor fuel in the premixing-prevaporizing passage by the direct injection method is described. The technical approach adopted for this program is to separate the problem into three parts each with its own computer code. These three parts are: calculation of the two-dimensional or axisymmetric air flow; calculation of the three-dimensional fuel droplet evaporation; and calculation of the fuel vapor diffusion. This method of approach is justified because premixing passages operate at lean equivalence ratios. Hence, a weak interaction assumption can be made wherein the airflow can affect the fuel droplet behavior but the fuel droplet behavior does not affect the airflow.

Anderson, O. L.↗
NASA NTRS1983

Kinetic effects on trace element partitioning

A potential difficulty regarding the application of trace element abundance data to problems in igneous petrogenesis is related to the possibility that partitioning of trace elements among minerals and melt may not be an equilibrium process. It may rather be controlled by kinetic processes, specifically the interplay between crystal growth and diffusion of the trace elements in crystals and melt. The present investigation has the objective to quantify kinetic effects on trace element partitioning. Crystals of olivine, pyroxene, and plagioclase were grown rapidly from synthetic melts resembling naturally-occurring mafic melts and doped with La, Sm, and Lu to quantify the magnitude of the effect of crystal growth rate on concentration of excluded trace elements. Olivine melt partitioning failed to show any appreciable kinetic effects, while diopside/melt distribution coefficients increased by 70 percent when the crystals grew very rapidly.

Lindstrom, D. J.↗