Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLARIFY”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Trisodium Phosphate Phases and Solubility in Alkaline Solutions Relevant to Radioactive Waste Processing

The U.S. Department of Energy’s Hanford Site faces significant challenges in managing millions of gallons of legacy radioactive waste, where phosphate precipitation can obstruct pipelines during retrieval and processing. To resolve long-standing inconsistencies in reported solubility and clarify factors governing trisodium phosphate hydrates, we examined the Na3PO4:NaOH:H2O system. Powder and single-crystal X-ray diffraction revealed that commercial precursors undergo transformations that produce multiple hydrates, including three previously unreported phases comprising an ordered polymorph of Na3(PO4)·12H2O·1/6NaOH, Na3PO4·5H2O, and Na3PO4·9H2O. Computational modeling indicated that the ordered dodecahydrate is more stable than its disordered counterpart, suggesting kinetic persistence of structural disorder. Solubility measurements were conducted with solutions prepared either from anhydrous Na3PO4 or from Na3(PO4)·12H2O·1/6NaOH and revealed that release of interstitial NaOH from the hydrate precursor elevated solution alkalinity, thereby significantly reducing phosphate solubility relative to solutions prepared with anhydrous Na3PO4 across 20–44 °C. These results begin to reconcile inconsistencies in prior solubility data and clarify how phase composition dictates phosphate precipitation under alkaline conditions.

Graham, Trenton R. (ORCID:0000000189078004)↗
DOE OSTI2022

Exotic Carbonate Mineralization Recovered from a Deep Basalt Carbon Storage Demonstration

Mitigating climate change require transformational advances for carbon dioxide removal, including geologic carbon sequestration in reactive subsurface environments. The Wallula Basalt Carbon Storage Pilot Project demonstrated that CO 2 injected into >800 m deep Columbia River Basalt Group flow top reservoirs mineralizes on month-year timescales. Herein, we present new optical petrography, micro-computed X-ray tomography, and electron microscopy results sidewall cores collected two years after CO 2 injection. As no other anthropogenic carbonates from geologic carbon storage field studies have been recovered, this world-unique sample suite provides unparalleled insight for subsurface carbon mineralization products and paragenesis. Chemically-zoned nodules with Ca/Mn-rich cores and Fe-dominant outer rims are prominent examples of the neoformed carbonate assemblages with ankerite-siderite compositions and exotic divalent cation correlations. Further, paragenetic insights for the timing of aragonite, silica, and fibrous zeolites are clarified based on mineral texture and spatial relationships, along with time-resolved downhole fluid sampling. Collectively, these results clarify the mineralogy, chemistry, and paragenesis of carbon mineralization, providing insight into the ultimate fate and transport of CO 2 in reactive mafic-ultramafic reservoirs

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Exploring Ligand-to-Metal Charge-Transfer States in the Photo-Ferrioxalate System Using Excited-State Specific Optimization

The photo-ferrioxalate system (PFS), [Fe(III)(C 2 O 4 )] 3– , more than an exact chemical actinometer, has been extensively applied in wastewater and environment treatment. Despite many experimental efforts to improve clarity, important aspects of the mechanism of ferrioxalate photolysis are still under debate. In this paper, we employ the recently developed WΓ-CASSCF to investigate the ligand-to-metal charge-transfer states that are key to ferrioxalate photolysis. This investigation provides a qualitative picture of these states and key potential energy surface features related to the photolysis. Furthermore, our theoretical results are consistent with the prompt charge-transfer picture seen in recent experiments and clarify some features that are not visible in experiments. Two ligand-to-metal charge-transfer states contribute to the photolysis of ferrioxalate, and the avoided crossing barrier between them is low compared with the initial photoexcitation energy. Our data also clarify that one Fe–O bond cleaves first, followed by the C–C bond and the other Fe–O bond.

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

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Realism of Lagrangian Large Eddy Simulations Driven by Reanalysis Meteorology: Tracking a Pocket of Open Cells Under a Biomass Burning Aerosol Layer

An approach to drive Lagrangian large eddy simulation (LES) of boundary layer clouds with reanalysis data is presented and evaluated using satellite (Spinning Enhanced Visible and Infrared Imager, SEVIRI) and aircraft (Cloud-Aerosol-Radiation Interactions and Forcing, CLARIFY) measurements. The simulations follow trajectories of the boundary layer flow. They track the formation and evolution of a pocket of open cells (POC) underneath a biomass burning aerosol layer in the free troposphere. The simulations reproduce the evolution of observed stratocumulus cloud morphology, cloud optical depth, and cloud drop effective radius, and capture the timing of the cloud state transition from closed to open cells seen in the satellite imagery on the three considered trajectories. They reproduce a biomass burning aerosol layer identified by the in-situ aircraft measurements above the inversion of the POC. Entrainment of aerosol from the biomass burning layer into the POC is limited to the extent of having no impact on cloud- or boundary layer properties, in agreement with the CLARIFY observations. The two-moment bin microphysics scheme used in the simulations reproduces the in-situ cloud microphysical properties reasonably well. A two-moment bulk microphysics scheme reproduces the satellite observations in the non-precipitating closed-cell state, but overestimates liquid water path and cloud optical depth in the precipitating open-cell state due to insufficient surface precipitation. A boundary layer cold and dry bias occurring in LES can be counteracted by reducing the grid aspect ratio and by tightening the large scale wind speed nudging towards the surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Oxygen redox activities governing high-voltage charging reversibility of Ni-rich layered cathodes

A schematic illustration indicative of the improved structural reversibility through stabilizing oxidized oxygen by Al doping. The chemical reactions and phase transitions at high voltages are generally considered to determine the electrochemical properties of high-voltage layered cathodes such as Ni-rich rhombohedral oxides. Even if significantly higher SOCs (states-of-charge) are utilized above the capability of transition metal redox (primarily Ni and Co), the effect of oxygen redox on Ni-rich rhombohedral oxides still looks mysterious thereby necessitating research that can clarify the relationship between redox reactions and phase transitions. Here, we performed a comprehensive and comparative study of the cationic and anionic redox reactions, as well as the structural evolution of a series of commercial Ni-rich layered oxides with and without Al doping. We combined the results from X-ray spectroscopy, operando electrochemical mass spectrometry, and neutron diffraction with electrochemical properties and thereby revealed the different oxygen redox activities associated with structural and electrochemical degradations. We reveal that Al doping suppresses the irreversible oxygen release but enhances the lattice oxygen oxidization. With this modulated oxygen redox activity, the Ni-rich layered oxides’ notorious H2–H3 structural phase transition becomes highly reversible. Our findings disentangle the different oxygen redox activities during high-voltage cycling and clarify the role of dopants in the Ni-rich layered oxides in terms of structural and electrochemical stability, shedding light on the future direction of optimizing layered cathode materials for safer high energy-density secondary batteries.

Lee, Gi-Hyeok↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Examination of probability distribution of mixture fraction in LES/FDF modelling of a turbulent partially premixed jet flame

An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.

42 ENGINEERING↗
DOE OSTI2021

Laser scattering by submicron droplets formed during the electrical explosion of thin metal wires

Abstract This paper presents the results of studying dispersed media formation during the electrical explosion of thin metal wires in vacuum by using low-current generators (~1–10 kA). Particular attention is paid to the analysis of the composition and structure of the corresponding explosion products as well as to the problem of their visualization using simultaneous laser interferometry and shadow imaging at two wavelengths (1.064 µ m and 0.532 µ m). Our findings indicate the important role in the visualization of the explosion products that belongs to multiple scattering by submicron droplets of dense condensed matter, which are mixed with metal vapor. The hypothesis on the existence of submicron droplets in the products of exploding metal wires correlates with the results obtained by soft x-ray radiography combined with a laser probing technique. Taking into account the multiple scattering by submicron droplets, it is possible to significantly clarify the parameters of the explosion products visualized via laser probing techniques as well as to gain a deeper insight into the physics behind the electrical wire explosion.

Physics↗
DOE OSTI2021

Prediction of crystal structures and motifs in the Fe–Mg–O system at Earth’s core pressures

Abstract Fe, Mg, and O are among the most abundant elements in terrestrial planets. While the behavior of the Fe–O, Mg–O, and Fe–Mg binary systems under pressure have been investigated, there are still very few studies of the Fe–Mg–O ternary system at relevant Earth’s core and super-Earth’s mantle pressures. Here, we use the adaptive genetic algorithm (AGA) to study ternary Fe x Mg y O z phases in a wide range of stoichiometries at 200 GPa and 350 GPa. We discovered three dynamically stable phases with stoichiometries FeMg 2 O 4 , Fe 2 MgO 4, and FeMg 3 O 4 with lower enthalpy than any known combination of Fe–Mg–O high-pressure compounds at 350 GPa. With the discovery of these phases, we construct the Fe–Mg–O ternary convex hull. We further clarify the composition- and pressure-dependence of structural motifs with the analysis of the AGA-found stable and metastable structures. Analysis of binary and ternary stable phases suggest that O, Mg, or both could stabilize a BCC iron alloy at inner core pressures.

Wang, Renhai↗
DOE OSTI2021

Static response and Love numbers of Schwarzschild black holes

We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild-(A)dS background in arbitrary dimensions. Here, we then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers — in particular that they vanish for all types of perturbation in four spacetime dimensions — but also leads to new results. For instance, we find that neutral Schwarzschild black holes polarize in the presence of an electromagnetic background in any number of spacetime dimensions except four. Moreover, we calculate for the first time black hole Love numbers for vector-type gravitational perturbations in higher dimensions and find that they generically do not vanish. Along the way, we shed some light on an apparent discrepancy between previous results in the literature, and clarify some aspects of the matching between perturbative calculations of static response on a Schwarzschild background and the point-particle effective theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Toward holistic understanding of the ITER-like resonant magnetic perturbation (RMP) ELM control on KSTAR

Here, KSTAR has clarified a set of unresolved 3D physics issues utilizing the ITER-like in-vessel, three-row, resonant magnetic perturbation (RMP) configurations. Since RMP-driven, edge-localized-modes (ELMs)-crash control elevates the divertor heat flux peak through its impact on edge plasma parameters and transport, a series of intentionally misaligned RMP configurations (IMCs) have been explored to investigate the relationship between RMP ELM control and divertor heat fluxes, while searching for an ideal IMC that could be favorable in both aspects. First of all, the contrasting influence of kink vs anti-kink phasing on the ELM-crash suppression has been articulated, demonstrating the synergistic benefit of 'kink' phasing on ELM-crash-suppression. On the other hand, the three-row IMC in the anti-kink phasing becomes more insensitive to the ELM-crashes at the sub-marginal level of RMP, consistent with theory. Meanwhile, the divertor 'wetted' area of ELM-crash-suppression gets narrower than that of ELM-crash-mitigation, suggesting that ELM-crash-mitigation remains advantageous over ELM-crash-suppression in terms of time-averaged divertor thermal loading. In comparison, based on a set of two-row IMCs, no evidence of divertor heat flux broadening was found during ELM-crash-suppression, supporting a hypothesis that the dispersal of the divertor heat flux in three-row IMCs cannot be driven by helically structured two-row RMPs alone. Among ITER-like three-rows, lower two-row RMPs have been found to be much more effective in suppressing the ELM-crashes than upper two-row RMPs. Although it is quite preliminary, the up/down asymmetric dependence of RMP coupling may be generically attributed to lower-single-null plasmas. Such a holistic understanding of RMP-driven, ELM-crash-control in KSTAR is expected not only to elucidate various subtle points in the vicinity of ELM-crash-suppression, but also to clarify the relevant divertor thermal loading issues for ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Ecosystem as Super-Organ/ism, Revisited: Scaling Hydraulics to Forests under Climate Change

Synopsis Classic debates in community ecology focused on the complexities of considering an ecosystem as a super-organ or organism. New consideration of such perspectives could clarify mechanisms underlying the dynamics of forest carbon dioxide (CO2) uptake and water vapor loss, important for predicting and managing the future of Earth’s ecosystems and climate system. Here, we provide a rubric for considering ecosystem traits as aggregated, systemic, or emergent, i.e., representing the ecosystem as an aggregate of its individuals or as a metaphorical or literal super-organ or organism. We review recent approaches to scaling-up plant water relations (hydraulics) concepts developed for organs and organisms to enable and interpret measurements at ecosystem-level. We focus on three community-scale versions of water relations traits that have potential to provide mechanistic insight into climate change responses of forest CO2 and H2O gas exchange and productivity: leaf water potential (Ψcanopy), pressure volume curves (eco-PV), and hydraulic conductance (Keco). These analyses can reveal additional ecosystem-scale parameters analogous to those typically quantified for leaves or plants (e.g., wilting point and hydraulic vulnerability) that may act as thresholds in forest responses to drought, including growth cessation, mortality, and flammability. We unite these concepts in a novel framework to predict Ψcanopy and its approaching of critical thresholds during drought, using measurements of Keco and eco-PV curves. We thus delineate how the extension of water relations concepts from organ- and organism-scales can reveal the hydraulic constraints on the interaction of vegetation and climate and provide new mechanistic understanding and prediction of forest water use and productivity.

Zoology↗
DOE OSTI2020

Effective viscosities in a hydrodynamically expanding boost-invariant QCD plasma

Background: The near-equilibrium properties of a QCD plasma can be encoded into transport coefficients such as bulk and shear viscosity. In QCD, the ratio of these transport coefficients to entropy density, ζ/$\textit{s}$ and η/$\textit{s}$, depends nontrivially on the plasma's temperature. This is unlike in conformal systems where they take constant values such as η/$\textit{s}$ = 1/(4π). Purpose: In this work, we show that in a (0 + 1)D boost-invariant fluid with no transverse expansion, a temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ can be described by an equivalent effective viscosity $\langleζ/s\rangle_{\text{eff}}$ or $\langleη/s\rangle_{\text{eff}}$. This effective viscosity combines the actual temperature-dependent ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ with the temperature profile of the fluid. Further, we extend the concept of effective viscosity in systems with transverse expansion and discuss how effective viscosities can be used to identify families of ζ/$\textit{s(T)}$ and η /$\textit{s(T)}$ that lead to similar hydrodynamic evolution. Methods: The Navier-Stokes relativistic hydrodynamic equations are used to provide a first definition of effective viscosity, in (0 + 1)D and (1 + 1)D. In the (0 + 1)D case, the analysis is extended to Israel-Stewart-type second-order hydrodynamics to clarify the effect of higher-order hydrodynamics corrections on the effective viscosity. Results: In a boost-invariant fluid with no transverse expansion [(0 + 1)D], the effective viscosity is expressed as a simple integral of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ over temperature, with a weight determined by the speed of sound of the fluid. The result is general for any equation of state with a moderate temperature dependence of the speed of sound, including the QCD equation of state. This definition of effective viscosity can be used to identify infinite families of ζ/$\textit{s(T)}$ or η/$\textit{s(T)}$ that produce essentially indistinguishable temperature profiles. In a boost-invariant cylindrical system [(1 + 1)D], a similar definition of effective viscosity is obtained in terms of characteristic trajectories in time and transverse direction. This leads to an infinite number of constraints on an infinite functional space for ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$. Realistic examples are presented by using a finite number of constraints on a finite functional space. Conclusions: The definition of effective viscosity in a (0 + 1)D system clarifies how infinite families of ζ/$\textit{s(T)}$ and η/$\textit{s(T)}$ can result in nearly identical hydrodynamic temperature profiles. By extending the study to a boost-invariant cylindrical [(1 + 1)D] fluid, we identify an approximate but more general definition of effective viscosity that highlights the potential and limits of the concept of effective viscosity in fluids with limited symmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Global patterns of forest autotrophic carbon fluxes

Carbon (C) fixation, allocation, and metabolism by trees set the basis for energy and material flows in forest ecosystems and define their interactions with Earth’s changing climate. However, we lack a cohesive synthesis on how forest carbon fluxes vary globally with respect to climate and one another. Here, we draw upon 1,319 records from the Global Forest Carbon Database (ForC), representing all major forest types and the nine most significant autotrophic carbon fluxes, to comprehensively explore how C cycling in mature, undisturbed forests varies with latitude and climate on a global scale. We show that, across all flux variables analyzed, C cycling decreases continuously with absolute latitude – a finding that confirms multiple previous studies but contradicts the idea that net primary productivity of temperate forests rivals that of tropical forests. C flux variables generally displayed similar trends across latitude and multiple climate variables, with no differences in allocation detected at this global scale. Temperature variables in general, and mean annual temperature and temperature seasonality in particular, were the best univariate predictors of C flux, explaining 19 - 71% of variation in the C fluxes analyzed. The effects of temperature were modified by moisture availability, with C flux reduced under hot and dry conditions and sometimes under very high precipitation. C fluxes increased with growing season length, but this was never the best univariate predictor. Within the growing season, the influence of climate on C cycling was small but significant for a number of flux variables. These findings clarify how forest C flux varies with latitude and climate on a global scale. In a period of accelerating climatic change, this improved understanding of the fundamental climatic controls on forest C cycling sets a foundation for understanding patterns of change.

Banbury Morgan, Rebecca↗