Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

First lithium experiments in HIDRA and evidence of helium retention during quasi-steady-state stellarator plasma operations

Recent experiments in Hybrid Illinois Device for Research and Applications (HIDRA) have had operational discharges between t discharge = 60 and 1000 s using electron cyclotron resonant heating (ECRH) of the plasma. This means that quasi-steady-state plasma discharges reach conditions to study long-pulse plasma material interactions (PMIs). The newly commissioned HIDRA-Material Analysis Test-stand PMI diagnostic is used to place a drop of lithium onto a heated tungsten surface, transfer the sample in-vacuo and expose it in a helium plasma. Helium is of interest as there is an open question to whether lithium will be able to remove helium ash in real fusion devices. The introduction of the W-Li sample in HIDRA resulted in evaporation of lithium into the helium plasma during a 600 s pulse and caused a reduction of over 90% in neutral pressure during the discharge. It was also observed that the plasma density and temperature increased by over 2.5 times. Using spectroscopy and a helium collisional radiative model, the peak temperature and density of the helium plasma can be monitored during the discharge. During lithium evaporation, as significant lithium ionization occurs, there is a 85% drop in the HIDRA vessel neutral pressure, despite a constant flow rate of He gas. This reduction in neutral pressure is supported by spectroscopy data with corresponding reductions in He I line intensities (587 nm, 667 nm, 706 nm, and 728 nm), as well as those of other impurities. At one point in the discharge a lithium plasma is created, as indicated by an increase in Li + emission and a complete reduction in He + emission, but the electron density jumps from ne = 3 × 10 18 m –3 to over ne = 8 × 10 18 m –3 while the core temperature stays relatively constant between Te = 16 eV and 20 eV. Once lithium has completely evaporated from the sample and the majority of the ionized lithium has diffused from the plasma to the vessel walls, pressure and spectroscopy data paired with He collisional radiative model calculations shows a re-establishment of a helium plasma in a low recycling regime. In this regime, the density drops down to ne = 2 × 10 18 m –3 and the electron temperature increases from T e = 20 eV to over T e = 50 eV indicating an increase in helium heating efficiency. This is also indicated by the He+ emission re-establishing and having a higher intensity. Here, we show the results from the first lithium campaign in HIDRA. In the presence of lithium, and in particular when lithium ions are present, the helium disappears from the plasma via an as of yet unknown complex relationship that needs to be further studied. The most likely explanation is that the lithium ions are distributed around the vessel and able to trap helium to the surface turning HIDRA into a large gettering surface. These results have potential implications on future plasma facing component design using liquid lithium for impurity and recycling control using limiters and divertors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2014

Explanation of the local galactic cosmic ray energy spectra measured by Voyager 1. I. Protons

Almost exactly 100 yr after the original discovery of cosmic rays, the V1 spacecraft has observed, for the first time, the local interstellar medium energy spectra of cosmic ray H, He, C/O nuclei at nonrelativistic kinetic energies, after leaving the heliosphere modulation region on 2012 August 25. We explain these observations by modeling the propagation of these particles in the local Galactic environment with an updated steady-state spatial diffusion model including all particle momentum losses with the local interstellar gas (Coulomb/ionization, pion production, adiabatic deceleration, and fragmentation interactions). Excellent agreement with the V1 cosmic ray H observations is obtained if the solar system resides within a spatially homogeneous layer of distributed cosmic ray sources injecting the same momentum power law ∝p {sup –s} with s = 2.24 ± 0.12. The best fit to the V1 H observations also provides an estimate of the characteristic break kinetic energy T{sub C} = 116 ± 27 MeV, representing the transition from ionization/Coulomb energy losses at low energies to pion production and adiabatic deceleration losses in a Galactic wind at high energies. As the determined value is substantially smaller than 217 MeV in the absence of adiabatic deceleration, our results prove the existence of a Galactic wind in the local Galactic environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Geometries of edge and mixed dislocations in bcc Fe from first principles calculations

We use density functional theory (DFT) to compute the core structures of $a_0[100](010)$ edge, $a_0[100](011)$ edge, $a_0/2[\bar{1}\bar{1}1](1\bar{1}0)$ edge, and $a_0/2[111](1\bar{1}0)$ $71^{\circ}$ mixed dislocations in body-centered cubic (bcc) Fe. The calculations are performed using flexible boundary conditions (FBC), which effectively allow the dislocations to relax as isolated defects by coupling the DFT core to an infinite harmonic lattice through the lattice Green function (LGF). We use the LGFs of the dislocated geometries in contrast to most previous FBC-based dislocation calculations that use the LGF of the bulk crystal. The dislocation LGFs account for changes in the topology of the crystal in the core as well as local strain throughout the crystal lattice. A simple bulk-like approximation for the force constants in a dislocated geometry leads to dislocation LGFs that optimize the core structures of the $a_0[100](010)$ edge, $a_0[100](011)$ edge, and $a_0/2[111](1\bar{1}0)$ $71^{\circ}$ mixed dislocations. This approximation fails for the $a_0/2[\bar{1}\bar{1}1](1\bar{1}0)$ dislocation however, so in this case we derive the LGF from more accurate force constants computed using a Gaussian approximation potential. The standard deviations of the dislocation Nye tensor distributions quantify the widths of the dislocation cores. The relaxed cores are compact, and the local magnetic moments on the Fe atoms closely follow the volumetric strain distributions in the cores. We also compute the core structures of these dislocations using eight different classical interatomic potentials, and quantify symmetry differences between the cores using the Fourier coefficients of their Nye tensor distributions. Most of the core structures computed using the classical potentials agree well with the DFT results. Furthermore the DFT core geometries provide benchmarking for classical potential studies of work-hardening, as well as substitutional and interstitial sites for computing solute-dislocation interactions that serve as inputs for mesoscale models of solute strengthening and solute diffusion near dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Origin of Rapid Delithiation In Secondary Particles Of LiNi 0.8 Co 0.15 Al 0.05 O 2 and LiNi y Mn z Co 1− y − z O 2 Cathodes

Abstract Most research on the electrochemical dynamics in materials for high‐energy Li‐ion batteries has focused on the global behavior of the electrode. This approach is susceptible to misleading analyses resulting from idiosyncratic kinetic conditions, such as surface impurities inducing an apparent two‐phase transformation within LiNi 0.8 Co 0.15 Al 0.05 O 2 . Here, nano‐focused X‐ray probes are used to measure delithiation operando at the scale of secondary particle agglomerates in layered cathode materials during charge. After an initial latent phase, individual secondary particles undergo rapid, stochastic, and largely uniform delithiation, which is in contrast with the gradual increase in cell potential. This behavior reproduces across several layered oxides. Operando X‐ray microdiffraction (‐XRD) leverages the relationship between Li content and lattice parameter to further reveal that rate acceleration occurs between Li‐site fraction ( x Li ) ≈0.9 and ≈0.5 for LiNi 0.8 Co 0.15 Al 0.05 O 2 . Physics‐based modeling shows that, to reproduce the experimental results, the exchange current density ( i 0 ) must depend on x Li , and that i 0 should increase rapidly over three orders of magnitude at the transition point. The specifics and implications of this jump in i 0 are crucial to understanding the charge‐storage reaction of Li‐ion battery cathodes.

25 ENERGY STORAGE↗
DOE OSTI2017

A reduced order model to analytically infer atmospheric CO 2 concentration from stomatal and climate data

To address questions related to the acceleration or deceleration of the global hydrological cycle or links between the carbon and water cycles over land, reliable data for past climatic conditions based on proxies are required. In particular, the reconstruction of palaeoatmospheric CO 2 content (C a ) is needed to assist the separation of natural from anthropogenic C a variability and to explore phase relations between C a and air temperature T a time series. Both T a and C a are needed to fingerprint anthropogenic signatures in vapor pressure deficit, a major driver used to explain acceleration or deceleration phases in the global hydrological cycle. Current approaches to C a reconstruction rely on a robust inverse correlation between measured stomatal density in leaves (ν) of many plant taxa and C a . There are two methods that exploit this correlation: The first uses calibration curves obtained from extant species assumed to represent the fossil taxa, thereby restricting the suitable taxa to those existing today. The second is a hybrid eco-hydrological/physiological approach that determines C a with the aid of systems of equations based on quasi-instantaneous leaf-gas exchange theories and fossil stomatal data collected along with other measured leaf anatomical traits and parameters. In this contribution, a reduced order model (ROM) is proposed that derives C a from a single equation incorporating the aforementioned stomatal data, basic climate (e.g. temperature), estimated biochemical parameters of assimilation and isotope data. The usage of the ROM is then illustrated by applying it to isotopic and anatomical measurements from three extant species. The ROM derivation is based on a balance between the biochemical demand and atmospheric supply of CO 2 that leads to an explicit expression linking stomatal conductance to internal CO 2 concentration (C i ) and C a . The resulting expression of stomatal conductance from the carbon economy of the leaf is then equated to another expression derived from water vapor gas diffusion that includes anatomical traits. When combined with isotopic measurements for long-term C i /C a , C a can be analytically determined and is interpreted as the time-averaged C a that existed over the life-span of the leaf. Key advantages of the proposed ROM are: 1) the usage of isotopic data provides constraints on the reconstructed atmospheric CO 2 concentration from ν, 2) the analytical form of this approach permits direct links between parameter uncertainties and reconstructed C a , and 3) the time-scale mismatch between the application of instantaneous leaf-gas exchange expressions constrained with longer-term isotopic data is reconciled through averaging rules and sensitivity analysis. The latter point was rarely considered in prior reconstruction studies that combined models of leaf-gas exchange and isotopic data to reconstruct C a from ν. The proposed ROM is not without its limitations given the need to a priori assume a parameter related to the control on photosynthetic rate. Our work here further explores immanent constraints for the aforementioned photosynthetic parameter.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

PROPERTIES OF NEARBY STARBURST GALAXIES BASED ON THEIR DIFFUSE GAMMA-RAY EMISSION

The physical relationship between the far-infrared and radio fluxes of star-forming galaxies has yet to be definitively determined. The favored interpretation, the 'calorimeter model', requires that supernova generated cosmic-ray (CR) electrons cool rapidly via synchrotron radiation. However, this cooling should steepen their radio spectra beyond what is observed, and so enhanced ionization losses at low energies from high gas densities are also required. Further, evaluating the minimum energy magnetic field strength with the traditional scaling of the synchrotron flux may underestimate the true value in massive starbursts if their magnetic energy density is comparable to the hydrostatic pressure of their disks. Gamma-ray spectra of starburst galaxies, combined with radio data, provide a less ambiguous estimate of these physical properties in starburst nuclei. While the radio flux is most sensitive to the magnetic field, the GeV gamma-ray spectrum normalization depends primarily on gas density. To this end, spectra above 100 MeV were constructed for two nearby starburst galaxies, NGC 253 and M82, using Fermi data. Their nuclear radio and far-infrared spectra from the literature are compared to new models of the steady-state CR distributions expected from starburst galaxies. Models with high magnetic fields, favoring galaxy calorimetry, are overall better fits to the observations. These solutions also imply relatively high densities and CR ionization rates, consistent with molecular cloud studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Gravitational encounters and the evolution of galactic nuclei. IV. Captures mediated by gravitational-wave energy loss

Direct numerical integrations of the two-dimensional Fokker–Planck equation are carried out for compact objects orbiting a supermassive black hole at the center of a galaxy. As in Papers I–III, the diffusion coefficients incorporate the effects of the lowest-order post-Newtonian corrections to the equations of motion. In addition, terms describing the loss of orbital energy and angular momentum due to the 5/2-order post-Newtonian terms are included. In the steady state, captures are found to occur in two regimes that are clearly differentiated in terms of energy, or semimajor axis; these two regimes are naturally characterized as “plunges” (low binding energy) and “EMRIs,” or extreme-mass-ratio inspirals (high binding energy). The capture rate, and the distribution of orbital elements of the captured objects, are presented for two steady-state models based on the Milky Way: one with a relatively high density of remnants and one with a lower density. In both models, but particularly in the second, the steady-state f-bar (E) and the distribution of orbital elements of the captured objects are substantially different than if the Bahcall–Wolf energy distribution were assumed. The ability of classical relaxation to soften the blocking effects of the Schwarzschild barrier is quantified. These results, together with those of Papers I–III, suggest that a Fokker–Planck description can adequately represent the dynamics of collisional loss cones in the relativistic regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Understanding Process–Structure Relationships during Lamination of Halide Perovskite Interfaces

Fabrication of halide perovskite (HP) solar cells typically involves the sequential deposition of multiple layers to create a device stack, which is limited by the thermal and chemical incompatibility of top contact layers with the underlying HP semiconductor. One emerging strategy to overcome these restrictions on material selection and processing conditions is lamination, where two half-stacks are independently processed and then diffusion bonded to complete the device. Lamination reduces the processing constraints on the top side of the solar cell to allow new device designs, expanded use of deposition methods, and self-encapsulation of devices. While laminated perovskite solar cells with high efficiencies and novel interlayer combinations have been demonstrated, there is a limited understanding of how the lamination process parameters affect the diffusion-bond quality and material properties of the resulting HP layer. In this study, we systematically vary temperature, pressure, and time during lamination and quantify the resulting impacts on bonded area, grain domain size, and photoluminescence. A design of experiments is performed, and statistical analysis of the experimental results is used to quantitatively evaluate the resulting process–structure–property relationships. The lamination temperature is found to be the key parameter controlling these properties. Furthermore, a temperature of 150 °C enables successful bonding over 95% of the substrate area and also results in increases in apparent grain domain size and photoluminescence intensity. Based on these insights, the lamination temperature of functional perovskite solar cell devices is varied, demonstrating the importance of the resulting bond quality on device performance metrics.

14 SOLAR ENERGY↗
DOE OSTI2024

Review of the second charged-particle transport coefficient code comparison workshop

We report the results of the second charged-particle transport coefficient code comparison workshop, which was held in Livermore, California on 24–27 July 2023. This workshop gathered theoretical, computational, and experimental scientists to assess the state of computational and experimental techniques for understanding charged-particle transport coefficients relevant to high-energy-density plasma science. Data for electronic and ionic transport coefficients, namely, the direct current electrical conductivity, electron thermal conductivity, ion shear viscosity, and ion thermal conductivity were computed and compared for multiple plasma conditions. Additional comparisons were carried out for electron–ion properties such as the electron–ion equilibration time and alpha particle stopping power. Overall, 39 participants submitted calculated results from 18 independent approaches, spanning methods from parameterized semi-empirical models to time-dependent density functional theory. In the cases studied here, we find significant differences—several orders of magnitude—between approaches, particularly at lower temperatures, and smaller differences—roughly a factor of five—among first-principles models. We investigate the origins of these differences through comparisons of underlying predictions of ionic and electronic structure. The results of this workshop help to identify plasma conditions where computationally inexpensive approaches are accurate, where computationally expensive models are required, and where experimental measurements will have high impact.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Exploring NaCl-PuCl 3 molten salts with machine learning interatomic potentials and graph theory

Actinide molten salts are the basis of the liquid fuels used in molten salt reactors. Due to the inherent difficulties associated with high temperature and hazardous conditions, experimental investigations of fundamental properties of these materials are usually challenging. In this work, we describe the structure and transport of NaCl-PuCl 3 mixtures using computational techniques. Three compositions were considered (16, 25, and 36 mol% PuCl 3 ) over a temperature range (730 – 1257K) using ab initio molecular dynamics, which provided the necessary data sets for training machine learned interatomic potentials. Further, molecular dynamics simulations based on these potentials were then used to determine structure and transport properties. A substantial change was noted in the structure factor when increasing the PuCl 3 content from 25 to 36 mol%. This change is linked to the aggregation of larger Pu 3+ clusters. In addition, the similarity of the atomic environments of metal cations in molten salt systems to their solid states counterparts was investigated using an unsupervised learning technique. Finally, graph theory was employed to explore the structure and size of actinide networks. Consistent with the structure factor, a dense Pu 3+ intermolecular structure is observed within the 36 mol% PuCl 3 mixture. The structure of cation-cation inter-junctions is also discussed. In all cases, the diffusion of Pu 3+ is significantly lower than that of Na + and Cl - .

36 MATERIALS SCIENCE↗
DOE OSTI2016

Study on air ingress phenomena during a depressurization accident of a VHTR - 18632

A depressurization accident is the one of the design-basis accidents of a Very-High-Temperature Reactor (VHTR). When a depressurization accident occurs, air is expected to enter into the reactor pressure vessel from the breach and oxidize in-core graphite structures. Therefore, it is important to understand the mixing processes of different kinds of gases in the stable and unstable stratified fluid layers. In particular, it is also important to examine the influence of localized natural convection and molecular diffusion on the mixing process from a safety viewpoint. Therefore, in order to predict or analyze the air ingress phenomena during a depressurization accident, it is important to develop a method for the prevention of air ingress during an accident. We carried out experiments and numerical analysis using three-dimensional (3D) CFD code to obtain the mixing process of two-component gases and the flow characteristics of localized natural convection. This study also investigated a control method for the natural circulation of air through the injection of helium gas. The numerical model consists of a storage tank and a reverse U-shaped vertical rectangular passage. They are separated by a horizontal partition plate. One sidewall of the high-temperature side vertical passage is heated and the other sidewall is cooled. The low-temperature vertical passage is cooled by ambient air. The experimental procedures and numerical analysis process are as follows. Firstly, the storage tank is filled with heavy gas and the reverse U-shaped vertical passage is filled with a light gas. In the vertical passage of the high-temperature side, localized natural convection is generated by the temperature difference between the vertical walls. These flow characteristics were obtained from the experiments and steady-state analysis. The experiment and analysis began after the partition plate was opened. The results obtained from the experiments were quantitatively simulated using 3D numerical analysis. The two component gases were mixed via molecular diffusion and natural convection. After some time elapsed, natural circulation occurred through the reverse U-shaped vertical passage. These flow characteristics are the same as those of phenomena generated in the passage between a permanent reflector and a pressure vessel wall of the VHTR. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Alcove formation in dissolving cliffs driven by density inversion instability

We demonstrate conditions that give rise to cave-like features commonly found in dissolving cliffsides with a minimal two-phase physical model. Alcoves that are wider at the top and tapered at the bottom, with sharp-edged ceilings and sloping floors, are shown to develop on vertical solid surfaces dissolving in aqueous solvents. As evident from descending plumes, sufficiently large indentations evolve into alcoves as a result of the faster dissolution of the ceiling due to a solutal Rayleigh–Bénard density inversion instability. In contrast, defects of size below the boundary layer thickness set by the critical Rayleigh number smooth out, leading to stable planar interfaces. Furthermore, the ceiling recession rate and the alcove opening area evolution are shown to be given to first-order by the critical Rayleigh number. By tracking passive tracers in the fluid phase, we show that the alcoves are shaped by the detachment of the boundary layer flow and the appearance of a pinned vortex at the leading edge of the indentations. The attached boundary layer past the developing alcove is then found to lead to rounding of the other sides and the gradual sloping of the floor.

58 GEOSCIENCES↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

A trigonal nodal SP3 method with mesh refinement capabilities - Development and verification

The neutronics model of the nodal reactor dynamics code DYN3D developed for 3D analyses of steady states and transients in Light-Water Reactors has been extended by a simplified P{sub 3} (SP{sub 3}) neutron transport option - to overcome the limitations of the diffusion approach at regions with significant anisotropy effects. To provide a method being applicable to reactors with hexagonal fuel assemblies and to furthermore allow flexible mesh refinement, the nodal SP{sub 3} method has been developed on the basis of a flux expansion in triangular-z geometry. In this paper, the derivation of the trigonal SP{sub 3} method is presented in a condensed form and a verification of the methodology on quasi-pin level is performed by means of two single-assembly test examples. The corresponding pin-wise few-group cross sections were obtained by the deterministic lattice code HELIOS. The power distributions were calculated using both the trigonal DYN3D diffusion and SP{sub 3} solver and compared to the HELIOS reference solutions. Close to regions with non-negligible flux gradients, e.g., caused by the presence of a strong absorbing material, the power distribution calculated by DYN3D-SP{sub 3} shows a significant improvement in comparison to the diffusion method. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Accurate determination of electronic transport properties of silicon wafers by nonlinear photocarrier radiometry with multiple pump beam sizes

In this paper, photocarrier radiometry (PCR) technique with multiple pump beam sizes is employed to determine simultaneously the electronic transport parameters (the carrier lifetime, the carrier diffusion coefficient, and the front surface recombination velocity) of silicon wafers. By employing the multiple pump beam sizes, the influence of instrumental frequency response on the multi-parameter estimation is totally eliminated. A nonlinear PCR model is developed to interpret the PCR signal. Theoretical simulations are performed to investigate the uncertainties of the estimated parameter values by investigating the dependence of a mean square variance on the corresponding transport parameters and compared to that obtained by the conventional frequency-scan method, in which only the frequency dependences of the PCR amplitude and phase are recorded at single pump beam size. Simulation results show that the proposed multiple-pump-beam-size method can improve significantly the accuracy of the determination of the electronic transport parameters. Comparative experiments with a p-type silicon wafer with resistivity 0.1–0.2 Ω·cm are performed, and the electronic transport properties are determined simultaneously. The estimated uncertainties of the carrier lifetime, diffusion coefficient, and front surface recombination velocity are approximately ±10.7%, ±8.6%, and ±35.4% by the proposed multiple-pump-beam-size method, which is much improved than ±15.9%, ±29.1%, and >±50% by the conventional frequency-scan method. The transport parameters determined by the proposed multiple-pump-beam-size PCR method are in good agreement with that obtained by a steady-state PCR imaging technique.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

New Empirical Formulation for the Sublimational Breakup of Graupel and Dendritic Snow

Abstract Ice fragments are generated by sublimation of ice particles in subsaturated conditions in natural clouds. Conceivably, such sublimational breakup would be expected to cause ice multiplication in natural clouds. Any fragment that survives will grow to become ice precipitation that may sublimate and fragment further. As a first step toward assessing this overlooked process, a formulation is proposed for the number of ice fragments from sublimation of ice particles for an atmospheric model. This is done by amalgamating laboratory observations from previously published studies. The concept of a “sublimated mass activity spectrum” for the breakup is applied to the dataset. The number of ice fragments is determined by the relative humidity over ice and the initial size of the parent ice particles. The new formulation applies to dendritic crystals and heavily rimed particles only. Finally, a thought experiment is performed for an idealized scenario of subsaturation with in-cloud descent. Scaling analysis yields an estimate of an ice enhancement ratio of about 5 (10) within a weak deep convective downdraft of about 2 m s −1 , for an initial monodisperse population of dendritic snow (graupel) particles of 3 L −1 and 2 mm. During descent, there is a dynamic equilibrium between continual emission of fragments and their depletion by sublimation. A simplified bin microphysics parcel model exhibits this dynamical quasi equilibrium, consistent with the thought experiment. The fragments have average lifetimes of around 90 and 70 s for dendrites and graupel, respectively. Sublimational breakup is predicted to cause significant secondary ice production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Type I planet migration in a magnetized disk. I. Effect of large-scale vertical and azimuthal field components

We study the effects of a large-scale, ordered magnetic field in protoplanetary disks on Type I planet migration using a combination of numerical simulations in 2D and 3D and a linear perturbation analysis. Steady-state models of such disks require the inclusion of magnetic diffusivity. To make progress using ideal MHD, we focus on simplified field configurations, involving purely vertical (B{sub z}) and azimuthal (B{sub φ}) field components and a combination of the two. For each of the models we calculate the locations of the relevant resonances and of the turning points, which delineate the propagation regions of the MHD waves that transport angular momentum from the planet to the disk. We use both numerical and semianalytic methods to evaluate the cumulative back torque acting on the planet, and explore the effect of spatial gradients in the disk’s physical variables on the results. We conclude that, under realistic (3D) circumstances, a large-scale magnetic field can slow down the inward migration that characterizes the underlying unmagnetized disk—by up to a factor of ∼2 when the magnetic pressure approaches the thermal pressure—but it cannot reverse it. A previous inference that a pure-B{sub ϕ} field whose amplitude decreases fast enough with radius leads to outward migration applies only in 2D. In fact, we find that, in 3D, a pure-B{sub ϕ} disk undergoes a rapid transition to turbulence on account of a magnetorotational instability that is triggered by the planet-induced appearance of a weak B{sub z} component.

79 ASTRONOMY AND ASTROPHYSICS↗