SEARCH · Engineering Papers

Results for “flux analysis”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

DOE OSTI2021

OVERVIEW OF RECENT PU-238 PRODUCTION ACTIVITIES AT IDAHO NATIONAL LABORATORY

The Plutonium Fuel Services (PFS) program at Idaho National Laboratory (INL) is active in the qualification of irradiation targets containing Np-237 for irradiation in the Advanced Test Reactor (ATR) to produce Pu-238 for future NASA missions. INL qualified and loaded 7 targets in ATR’s South Flux Trap (SFT) for cycle 169A, which occurred in Spring 2021. This program was reinitiated after two baseline production targets in three positions validated significant production of Pu-238 [ref.1]. The validation model was followed by the PFS-1 experimental test in the ATRC (Critical) facility [ref.2]. This paper outlines the progress and status of the PFS program. The qualification effort, safety analysis, hardware status, and future activities for qualification of an updated target design for use in the ATR will be discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Assessment of MiniFuel Subcapsule Design Recommendations on Previous Experiments

MiniFuel describes the class of separate effects nuclear fuels irradiation experiments that have been conducted in the High Flux Isotope Reactor (HFIR) since 2018. These experiments comprise a stack of six fuel-bearing subcapsules contained in a stainless steel target housing that is in contact with HFIR coolant on its exterior. All MiniFuel targets have a near-standardized architecture, and the primary design variables that change between experiments are the radial gap size between the subcapsule and housing and the target fill gas composition. Finite element heat transfer models are used to determine the optimum gas composition and gap sizes, and recent studies were performed to identify model parameters that contribute the most uncertainty to fuel specimen temperature predictions. That work, which is referenced herein, also recommended a set of design modifications to the subcapsule internal architecture and assembly process. These modifications are intended to reduce fuel temperature uncertainty in future experiments. In this report, the subcapsule design modifications were retroactively applied to a previously conducted MiniFuel experiment to determine how these changes affect the established safety and performance envelope of the experimental capability. These effects were determined in two steps. First, the modifications were applied to the subcapsule design without any other changes to determine their isolated effect on the predicted fuel specimen temperatures. This portion of the analysis showed that fuel temperatures were modestly reduced because the implemented changes improved heat transfer efficacy. Next, traditional MiniFuel design activities (i.e., sizing the gas gaps and determining the fill gas composition) were reperformed, and they confirmed that the original desired fuel temperatures could be achieved while remaining within established safety limits. Therefore, this report demonstrates improved performance resulting from the subcapsule modifications, which mitigate uncertainty while meeting the objectives of past experiments. An additional benefit of the design changes is reduced sensitivity of the fuel temperature to the evolving flux spectrum in HFIR, leading to more stable temperatures and enhanced utility of MiniFuel as a separate effects irradiation platform.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Linear instabilities in the hot-ion regime in a high-field spherical tokamak

Abstract The present operating high-field compact spherical tokamak ST40 is an important step toward an ST-based fusion reactor (Gryaznevich and Asunta 2017 Fusion Eng. Des. 123 177; McNamara 2023 Nucl. Fusion 63 054002). Temperature and density profiles and their uncertainties in recent hot ion ST40 plasmas with central ion temperatures exceeding 8.6 keV, have been inferred using an integrated analysis of several diagnostics including line-of-sight integration and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. A linear gyrokinetic stability analysis has been carried out to identify the most unstable micro-instabilities in these hot ion plasmas. In one of these plasmas, it is found that linear growth rates of both ion- and electron-scale ( k θ ρ s ⩾ 0.2 where k θ is the poloidal wavenumber and ρ s is the ion gyro-radius with sound speed) modes decrease from the edge toward the core of the plasma (i.e. from ρ = 0.8 to 0.3 where ρ is the square root of the normalized magnetic flux), while the change of linear growth rates at k θ ρ s < 0.2 is non-monotonic. In particular, at ρ = 0.3 no unstable mode was found at k θ ρ s > 5 , and about an order of magnitude reduction in the maximum linear growth rate in the ion-scale wavenumber range of 0.2 ⩽ k θ ρ s < 1 is seen at ρ = 0.3 compared with the ρ = 0.8 location. It is found through parametric scans that while the ion temperature gradient (ITG) mode/trapped electron mode (TEM) and kinetic ballooning mode (KBM)/kinetic shear Alfvén (KSA) mode are more important in the plasma core, ubiquitous mode (UM) is found to be the important ion-scale instability at ρ ⩾ 0.6 . At ρ = 0.8, two instabilities are found at the electron scale, with one being the typical electrostatic electron temperature gradient (ETG) mode at ρ e scale and beyond and the other being an unidentified low-real-frequency instability (at the intermediate scale between electron and ion gyroradii) which can propagate in the ion direction and is driven by ETG. The existence of UM at the ion scale and the low-real-frequency instability and ETG mode at the electron scale at ρ = 0.8 is supported by an extensive linear gyrokinetic stability analysis of another ST40 hot ion plasma in a wider parametric range, taking into account the uncertainties in the experimental profiles.

ubiquitous mode↗

DOE OSTI2020

Impact of nuclear reactor radiation on the performance of AlN/sapphire surface acoustic wave devices

The performance of an AlN/sapphire surface acoustic wave (SAW) delay line device was characterized in real time under irradiation inside a nuclear reactor. Both its resonant frequency and transmission efficiency were observed to respond to a change in reactor power. The response follows an exponentially saturating behavior after a step power increase, followed by an exponentially decaying recovery after reactor shutdown. A sensitivity analysis based on the governing electro-mechanical equations shows that the frequency shift can be attributed to the softening of sapphire’s elastic constants under neutron radiation. Additionally, a kinetic rate equation is adopted to interpret device response and describe its microstructural evolution. These results suggest that the AlN/sapphire SAW device remains functional under irradiation, is sensitive to neutron and gamma ray fluxes, and offers an opportunity for remote sensing and in-situ measurement of material properties when exposed to nuclear reactor environment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Multi-harmonic Rutherford island theory

Rutherford island theory, which governs the nonlinear evolution of tearing modes in tokamak plasmas, is generalized to take into account situations in which the conventional one-harmonic approximation is not valid. The analysis incorporates non-inductive currents driven by radio frequency (RF) electromagnetic waves injected into the plasma. A multi-harmonic tearing mode dispersion relation is derived that takes the form of a nonlinear inhomogeneous matrix eigenvalue problem. The dispersion relation is solved in the so-called two-harmonic approximation, in which only the principal Fourier harmonic of the perturbed magnetic flux and its first overtone are included in the calculation. In the absence of RF current drive, the nonlinear behavior of a tearing mode predicted in the two-harmonic approximation does not differ substantially from that predicted in the one-harmonic approximation. On the other hand, RF current drive that is sufficiently localized in the vicinity of the O-points of the mode's magnetic island chain is capable of triggering bifurcations of the O-points (which is impossible in the one-harmonic approximation). However, the current drive is incapable of triggering bifurcations of the island X-points. This finding is significant because Bardóczi and Evans recently observed bifurcations of magnetic island chain O-points in the presence of RF current drive in the DIII-D tokamak but did not observe bifurcations of the X-points. Finally, the changes in the topology of the magnetic island flux-surfaces induced by RF current drive are found to facilitate the stabilization of the tearing mode.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Transient computational and experimental thermal analysis of graphite foam monoblock

A densified graphite foam is being explored for its applicability as plasma-facing material in fusion devices. Three different graphite foam monoblocks are constructed and tested in the Garching Large Divertor Sample Testing Facility. The monoblock samples consist of graphite foam press-fit to a single tube, graphite foam cubes brazed to a single, and graphite foam cubes press-fit to a single tube. The tube is composed of CuCrZr with a steel twisted tape. The monoblocks are exposed to the heat fluxes of 5, 6, and 8 MW/m 2 for 30 s to determine the maximum surface and body temperatures measured with thermocouple for each monoblock design. The press-fit monoblocks are exposed to 8 MW/m 2 for 15 s for 100 cycles to determine the effect of thermal cycling on the contact between the graphite foam and the tube. STAR-CCM+ is used to predict how much the contact between the foam and tubes varies as a result of thermal cycling. In addition, the 6 MW/m 2 loading is modeled in STAR-CCM+ to compare the transient cooldown curves of the computational results to the recorded temperatures at the surface and two different thermocouple locations. This comparison is used to validate the temperature-dependent thermal conductivity and specific heat capacity used in the models of the graphite foam. The computational modeling of the experiment has been used to hypothesize ways to use and improve upon the graphite foam as a suitable material for fusion applications.

Gehrig, Monica↗

DOE OSTI2021

An Analysis of the Impact of Particle Growth on Transport Losses in Polymer-Electrolyte Fuel Cells

Voltage cycling causes catalyst nanoparticles in polymer-electrolyte fuel cells to grow. The concomitant loss of interfacial area results in larger kinetic, transport, and possibly ohmic overpotentials. This paper uses recently published experimental data and mathematical models to investigate the evolutions of transport and ohmic resistances to platinum nanoparticles located on the surface and inside the micropores of carbon black supports. Resistance to oxygen transport rises as surface area declines primarily because the flux to each remaining larger platinum particle increases. The path lengths governing oxygen diffusion to surface and buried platinum sites also increase as nanoparticles grow. Platinum nanoparticles on the surface become relatively less favorable as voltage cycling proceeds because they grow faster than platinum in micropores. Because voltage cycling causes total interfacial area to decline, and the fraction located inside micropores to increase, nanoscale ohmic losses increase as a catalyst layer decays. The practical importance of this effect is difficult to discern because proton conductivity in carbon micropores is not well characterized.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2020

Molecular Optoelectronics (Final Technical Report)

The project developed new theoretical tools for analysis of nanodevices subjected to external optical and bias driving. Traditional theoretical techniques for molecular optoelectronics are either based on kinetic schemes or utilize quasiparticle description. The former fails to account for quantum mechanics of the system-bath interface and may fail qualitatively when used in open interacting systems. The latter is inconvenient in (usual for molecular optical spectroscopy) many body states description of the system. We developed a new technique - the Hubbard nonequilibrium Green's functions method - which allows to avoid shortcomings of standard theoretical approaches. The Hubbard NEGF was applied to analysis of experimental measurements of photoinduced current and bias-induced electroluminescence in STM molecular junctions. We also demonstrated practical usefulness of the Hubbard NEGF within newly proposed universal theory of current-induced forces for nonadiabatic molecular dynamics which generalizes celebrated Head-Gordon and Tully expression for electronic friction. Finally, following latest trends in laser techniques we studied local responses (local fluxes and local noise spectroscopy) of biased junctions as a source of information not accessible by total responses (total fluxes and noises) of the system.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2023

Global evaluation of terrestrial biogeochemistry in the Energy Exascale Earth System Model (E3SM) and the role of the phosphorus cycle in the historical terrestrial carbon balance

Abstract. The importance of carbon (C)–nutrient interactions to the prediction of future C uptake has long been recognized. The Energy Exascale Earth System Model (E3SM) land model (ELM) version 1 is one of the few land surface models that include both N and P cycling and limitation (ELMv1-CNP). Here we provide a global-scale evaluation of ELMv1-CNP using the International Land Model Benchmarking (ILAMB) system. We show that ELMv1-CNP produces realistic estimates of present-day carbon pools and fluxes. Compared to simulations with optimal P availability, simulations with ELMv1-CNP produce better performance, particularly for simulated biomass, leaf area index (LAI), and global net C balance. We also show ELMv1-CNP-simulated N and P cycling is in good agreement with data-driven estimates. We compared the ELMv1-CNP-simulated response to CO2 enrichment with meta-analysis of observations from similar manipulation experiments. We show that ELMv1-CNP is able to capture the field-observed responses for photosynthesis, growth, and LAI. We investigated the role of P limitation in the historical balance and show that global C sources and sinks are significantly affected by P limitation, as the historical CO2 fertilization effect was reduced by 20 % and C emission due to land use and land cover change was 11 % lower when P limitation was considered. Our simulations suggest that the introduction of P cycle dynamics and C–N–P coupling will likely have substantial consequences for projections of future C uptake.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Impact of Temperature and Optical Error on the Combined Optical and Thermal Efficiency of Solar Tower Systems for Industrial Process Heat: Preprint

Concentrating solar thermal (CST) power towers can provide high flux concentrations at commercial scale. As a result, CST towers exhibit potential for high-temperature solar industrial process heat (SIPH) applications. However, at higher operating temperatures, thermal radiation losses can be significant. This study explores the trade-off between thermal and optical losses for SIPH applications using a collection of three case studies at operating temperatures that range from 900-1,550 degrees C. We assume blackbody radiation to represent the thermal losses at the receiver and we use ray tracing to estimate the optical losses. The results show the impact of process temperature on the maximum attainable system efficiency, as well as the higher flux concentration requirements as the temperature increases.

concentrating solar power↗

DOE OSTI2025

Measurement of $\nu_\mu$ CC Interactions With Two-Proton Final State in MINERvA

This dissertation presents a measurement of charged–current (CC) muon–neutrino interactions with exactly two protons and no pions in the final state (CC~$2p\,0\pi$), using data collected by the MINERvA detector in the NuMI medium–energy beam at Fermilab. Such two–proton topologies are a sensitive probe of nuclear dynamics in the few–GeV regime, including multi–nucleon correlations (npnh, notably $2p2h$) and intranuclear final–state interactions (FSI) such as pion absorption and nucleon rescattering. A precise experimental characterization of these processes is essential both for neutrino–interaction theory and for reducing systematic uncertainties in oscillation experiments that rely on accurate modeling of neutrino–nucleus interactions. Events are selected by requiring a $\nu_\mu$ CC interaction with a reconstructed $\mu^-$ and two proton tracks originating from a common vertex in MINERvA’s finely segmented scintillator tracker, with no reconstructed mesons. Muon charge and momentum are constrained by matching to the MINOS Near Detector, while proton identification exploits energy–loss profiles and stopping–proton features. Backgrounds from pion–producing channels that enter the signal region through FSI or reconstruction effects are constrained with data–driven sidebands (Michel–electron and isolated–cluster “blob” samples) and tuned via a simultaneous fit across signal and sideband regions. To correct detector resolution and acceptance effects, the analysis employs iterative Bayesian unfolding with extensive validation: statistical pseudo–experiments, and robustness checks against generator systematic “universes” and additional strong shape warps. Single–differential cross sections are reported for three observables tailored to the two–proton final state: the opening–angle cosine $\cos\!\left(\theta_{pp}\right)$, the leading–proton momentum, and the subleading–proton momentum. Systematic uncertainties include contributions from neutrino flux, interaction modeling (e.g., npnh and resonance parameters, pion FSI), and detector response (calibration, reconstruction efficiencies). The resulting distributions provide targeted constraints on the interplay of multi–nucleon dynamics and FSI that shape CC~$2p\,0\pi$ final states on hydrocarbon. Comparisons to modern GENIE–based simulations highlight kinematic regions where model components require refinement. These measurements thus inform generator tuning and improve the reliability of neutrino–energy reconstruction strategies for current and future long–baseline oscillation programs.

Syrotenko, Vladyslav S. [Tufts U.]↗

DOE OSTI2023

The ecosystem wilting point defines drought response and recovery of a Quercus‐Carya forest

Abstract Soil and atmospheric droughts increasingly threaten plant survival and productivity around the world. Yet, conceptual gaps constrain our ability to predict ecosystem‐scale drought impacts under climate change. Here, we introduce the ecosystem wilting point (Ψ EWP ), a property that integrates the drought response of an ecosystem's plant community across the soil–plant–atmosphere continuum. Specifically, Ψ EWP defines a threshold below which the capacity of the root system to extract soil water and the ability of the leaves to maintain stomatal function are strongly diminished. We combined ecosystem flux and leaf water potential measurements to derive the Ψ EWP of a Quercus‐Carya forest from an “ecosystem pressure–volume (PV) curve,” which is analogous to the tissue‐level technique. When community predawn leaf water potential (Ψ pd ) was above Ψ EWP (=−2.0 MPa), the forest was highly responsive to environmental dynamics. When Ψ pd fell below Ψ EWP , the forest became insensitive to environmental variation and was a net source of carbon dioxide for nearly 2 months. Thus, Ψ EWP is a threshold defining marked shifts in ecosystem functional state. Though there was rainfall‐induced recovery of ecosystem gas exchange following soaking rains, a legacy of structural and physiological damage inhibited canopy photosynthetic capacity. Although over 16 growing seasons, only 10% of Ψ pd observations fell below Ψ EWP , the forest is commonly only 2–4 weeks of intense drought away from reaching Ψ EWP , and thus highly reliant on frequent rainfall to replenish the soil water supply. We propose, based on a bottom‐up analysis of root density profiles and soil moisture characteristic curves, that soil water acquisition capacity is the major determinant of Ψ EWP , and species in an ecosystem require compatible leaf‐level traits such as turgor loss point so that leaf wilting is coordinated with the inability to extract further water from the soil.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A Case for Electron-Astrophysics

Abstract The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the research theme of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. In addition, plasma electrons often play an important role for the spatial transfer of thermal energy due to the high heat flux associated with their velocity distribution. The regulation of this electron heat flux is likewise not understood. By focussing on these and other fundamental electron processes, the research theme of electron-astrophysics links outstanding science questions of great importance to the fields of space physics, astrophysics, and laboratory plasma physics. In this White Paper, submitted to ESA in response to the Voyage 2050 call, we review a selection of these outstanding questions, discuss their importance, and present a roadmap for answering them through novel space-mission concepts.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Understanding Precipitation Bias Sensitivities in E3SM-Multi-Scale Modeling Framework From a Dilution Framework

We investigate a set of Energy Exascale Earth System Model Multi-scale modeling framework (MMF) (E3SM-MMF) simulations that vary the dimensionality and momentum transport configurations of the embedded cloud-resolving models (CRMs), including unusually ambitious 3D configurations. Issues endemic to all MMF simulations include too much Intertropical Convergence Zone rainfall and too little over the Amazon. Systematic MMF improvements include more on-equatorial rainfall across the Warm Pool. Interesting sensitivities to the CRM domain are found in the regional time-mean precipitation pattern over the tropics. The 2D E3SM-MMF produces an unrealistically rainy region over the northwestern tropical Pacific; this is reduced in computationally ambitious 3D configurations that use 1,024 embedded CRM grid columns per host cell. Trajectory analysis indicates that these regional improvements are associated with desirably fewer tropical cyclones and less extreme precipitation rates. To understand why and how the representation of precipitation improved in 3D, we propose a framework that dilution is stronger in 3D. This viewpoint is supported by multiple indirect lines of evidence, including a delayed moisture-precipitation pickup, smaller precipitation efficiency, and amplified convective mass flux profiles and more high clouds. We also demonstrate that the effects of varying embedded CRM dimensionality and momentum transport on precipitation can be identified during the first few simulated days, providing an opportunity for rapid model tuning without high computational cost. Meanwhile the results imply that other less computationally intensive ways to enhance dilution within MMF CRMs may also be strategic tuning targets.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A complete spectroscopic catalogue of local galaxies in the northern spring sky: Gas properties and nuclear activity in different environments

With the aim of providing the complete demography of galaxies in the local Universe, including their nuclear properties, we present SPRING, a complete census of local galaxies limited to the spring quarter of the northern sky (10h < RA < 16h; 0° < Dec < 65°). The SPRING catalogue is a flux- and volume-limited sample (r < 17.7 mag, cz < 10 000 km s -1 ) of 30 597 galaxies, including the Virgo, Coma, and A1367 clusters. Images and spectra were individually examined to clear the sample from unwanted entries. To inspect possible secular and environmental dependencies of the various nuclear excitation properties (star-forming versus active nuclei), we performed a multi-dimensional analysis by dividing the total sample according to: (i) their position in the (NUV - i) versus M star diagram, (ii) the local galaxy density, (iii) the stellar mass, (iv) the halo mass of the group to which galaxies belong, and (v) the neutral hydrogen content. We present a new calibration of the optical diameter-based H I-deficiency parameter H I def , employing a reference sample of isolated galaxies extracted from SPRING. At intermediate distances between Virgo and Coma, we identify a ring-like structure of galaxies constituted by three large filaments, each with a length of approximately 20h -1 Mpc, mostly composed of blue-cloud galaxies with stellar masses M star ≲ 10 10 M ⊙ . The fraction of H I-deficient galaxies within the filament (~30%) suggests that filaments are a transitioning environment between lower- and higher-overdensity environments in terms of H I content, as we find a clear progression from field galaxies to the filament and cluster galaxies for an increasing H I def parameter. We classify the nuclear spectra according to the four-line Baldwin-Phillips-Terlevich (BPT) and the two-line EWHα versus [NII]/Hα (WHAN) diagnostic diagrams, and investigate the variation in the fraction of active nuclei hosts with stellar mass, as well as their colours and environments. We observe that the fraction of low-ionisation nuclear emitting regions (LINERs) is a steep function of stellar mass, for example, it is consistent with zero up to M star ≲ 10 9.5 M ⊙ and becomes ~40% for M star ≳ 101 0.5 M ⊙ , whereas, for M star ≲ 10 9 - 9.5 M ⊙ , almost the entire spectroscopic sample is constituted of galaxies with star-forming nuclei. We investigate whether the nuclear-excitation fractions depend predominantly on the stellar mass or, conversely, on the galaxy environment. In general, we observe that the mass dependence of the fraction of Seyfert nuclei is not very sensitive to the galaxy environment, whereas the fraction of star-forming nuclei is a steeper function of stellar mass in lower-density environments and in blue-cloud galaxies. We find that the fraction of LINERs depends on galaxy colour and, for M star ≳ 10 9.5 - 10 M ⊙ , increases in galaxies belonging to the green valley.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Etched Silicon Planar CRL Optics for the High-Energy X-ray Diffraction Beamlines 11-ID-B and 11-ID-C at the APS

We report the focusing of high-energy X-rays is most commonly achieved by commercially available compound refractive lenses (CRL), which are added together in larger structures. There is a need for high-energy X-ray beamlines to enable focusing in only one dimension down to a few micrometers to study layered materials and thin films. High-energy X-rays above 50 keV require often 100 lenses or more, therefore it is more efficient to etch the complete arrays of lenses into a monolithic substrate like a silicon wafer. Custom lens arrays were fabricated for the high-energy x-ray beamlines 11-ID-B, and 11-ID-C at the Argonne National Laboratory. The fabrication required tuning parameters of the Bosch process to optimize the etching of deep vertical sidewalls at a reasonable speed without overheating. The beamline 11-ID-B is equipped with lens arrays of 85-110 lenses for 59 keV and 190-250 lenses for 87 keV photon energy, focusing the beam vertically below 5.0 µm at the sample positions between 1600 mm to 2100 mm. 11-ID-C has received arrays of 435 lenses that focus 106 keV X-rays to a vertical height of 2.0 µm at the sample position. In both cases, the flux density gain in the focal spot is larger than a factor of 20, relative to using slits to cut the beam size. The lens setup is permanently installed and can be moved in and out of the beam during the experiment without further alignment. The lenses have enabled operando diffraction studies on interfaces in batteries and the development of pair distribution function (PDF) analysis using grazing incidence geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Quasi-Steady Simulation of Glaze Ice Accretion and Heat Transfer in the Supercooled Large Droplet Regime in Atmospheric Turbulent Air Flow

Abstract To avoid the use of a computationally intensive full unsteady simulation while providing an accurate solution, a quasi-steady simulation has been performed to study glaze and mixed ice accretion in the supercooled large droplet (SLD) regime in turbulent flow of atmospheric air. We have attempted to find the minimum time-step necessary to adequately simulate the icing process on an airfoil surface. Based on node displacement, a mesh morphing scheme has been adopted in the computations to account for the moving boundaries that are caused by the continuous ice buildup. We have modeled ice accretion on an airfoil surface for a time period of 232 s using several time steps. At each time-step we solved the steady-state conservation equations for the air and droplet phases and then used the results as initial conditions for the next time position. The magnitude of the time-step ranged from a least accurate two-shot simulation (where the time-step is 116 s) to a most accurate 2320-shot simulation (where the time-step is 0.1 s). In-between these two extreme time steps, we have performed a three-shot simulation (where the time-step is 77.33 s), a four-shot simulation (where the time-step is 58 s), a six-shot simulation (where the time-step is 38.67 s), a 46-shot simulation (where the time-step is 5 s), and a 232-shot simulation (where the time-step is 1 s). We have done so to find out the degree of accuracy (or inaccuracy) of the multishot simulation approach and to find out the appropriate time-step needed for a successful and valid quasi-steady simulation. A valid quasi-steady simulation needs to use a time-step that is small enough to reproduce the full time-dependent solution within a very small error band. We have found that both the 1 s and the 0.1 s time steps produce virtually identical results. This is the primary litmus test that proves the validity of the quasi-steady-state assumption. The results adopted in this paper are thus all based on the more conservative 0.1 s time-step. In the process of performing the simulations, remeshing was required in order to maintain the grid density in zones of high curvature to be able to capture the full physics in those zones. After successfully modeling glaze ice accretion over the airfoil surface using the 0.1 s quasi-steady simulation approach, the effects of supercooled large droplets (SLDs) impacting the surface have been examined and presented in terms of the variation of the local collection efficiency, the water film thickness, and the heat transfer rate. Examination of the variation of the angle that the ice horn makes with the airfoil chord line demonstrated a 20% improvement in angle prediction when the time-step is reduced from 116 s to 0.1 s. The analysis also reveals a 12% or 8.5% increase in the maximum collection efficiency, βmax, depending on whether the value of the liquid water content (LWC) has been doubled from 0.5 g/m3 to 1 g/m3, or the value of the freestream velocity has been doubled from 75 m/s to 150 m/s, respectively. Because of the need to monitor the local collection efficiency and convective heat fluxes at each shot (regardless of the number of shots employed), the approach adopted here was found to be effective in successively and successfully reproducing the curvature of the glaze ice horn.

Engineering↗

DOE OSTI2020

Deep abiotic weathering of pyrite

Pyrite is a ubiquitous iron sulfide mineral that is oxidized by trace oxygen. The mineral has been largely absent from global sediments since the rise in oxygen concentration in Earth’s early atmosphere. We analyzed weathering in shale, the most common rock exposed at Earth’s surface, with chemical and microscopic analysis. By looking across scales from 10 −9 to 10 2 meters, we determined the factors that control pyrite oxidation. Under the atmosphere today, pyrite oxidation is rate-limited by diffusion of oxygen to the grain surface and regulated by large-scale erosion and clast-scale fracturing. We determined that neither iron- nor sulfur-oxidizing microorganisms control global pyrite weathering fluxes despite their ability to catalyze the reaction. This multiscale picture emphasizes that fracturing and erosion are as important as atmospheric oxygen in limiting pyrite reactivity over Earth’s history.

58 GEOSCIENCES↗