Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model order reduction”

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 307 records · Page 17

Direct statistical simulation of the Lorenz96 system in model reduction approaches

Direct statistical simulation (DSS) of nonlinear dynamical systems bypasses the traditional route of accumulating statistics by lengthy direct numerical simulations by solving the equations that govern the statistics themselves. DSS suffers, however, from the curse of dimensionality as the statistics (such as correlations) generally have higher dimensions than the underlying dynamical variables. Here we investigate two approaches to reduce the dimensionality of DSS, illustrating each method with numerical experiments with the Lorenz96 dynamical system. The forms of DSS chosen here involve approximate closures at second and third order in the equal-time cumulants. We demonstrate significant reduction in computational effort that can be achieved without sacrificing the accuracy of DSS. The methods developed here can be applied to turbulent fluid and magnetohydrodynamical systems. Published by the American Physical Society 2025

Li, Kuan↗

Green synthesis of sulfur‐doped g‐C 3 N 4 nanosheets for enhanced removal of oxytetracycline under visible‐light irradiation and reduction of its N ‐nitrosodimethylamine formation potential

Abstract BACKGROUND In this study, sulfur‐doped g‐C 3 N 4 (SCN) nanosheets were synthesized using thiourea as a precursor, and then used to catalyze the degradation of oxytetracycline (OTC) under visible‐light irradiation. OTC was chosen as target contaminant because of it being an important precursor of N ‐nitrosodimethylamine (NDMA). RESULTS The removal of OTC was influenced by its initial concentration, initial pH, SCN dosage and water matrices. OTC (≤10 mg L −1 ) was almost completely removed in 40 min under the conditions of an initial pH of 7.0 and a SCN dose of 1.0 g L −1 , and its degradation was well described by the pseudo‐first‐order kinetics model. Water matrices would also influence the degradation of OTC in the system. Electron paramagnetic resonance analysis and trapping experiments confirmed that the function of the radicals for OTC removal followed the order of • O 2 − > • OH > h + . Eight transformation products were determined during the degradation of OTC, it being removed through four pathways: methyl oxidization, demethylation, decarbonylation and secondary alcohol oxidation. SCN was stable in the system, as the removal efficiency of OTC was not significantly reduced even after three recycles. In addition, 55.5% of NDMA formation potential of OTC was reduced after SCN catalytic oxidation. CONCLUSIONS The results imply that the visible‐light‐driven SCN catalytic oxidation is an effective technology not only for the removal of OTC but also for the reduction of NDMA formation potential. © 2021 Society of Chemical Industry

Dou, Yicheng↗

Pt Particle Size Affects Both the Charge Separation and Water Reduction Efficiencies of CdS–Pt Nanorod Photocatalysts for Light Driven H 2 Generation

Decreasing the metal catalyst size into nanoclusters or even single atom is an emerging direction of developing more efficient and cost-effective photocatalytic systems. Because the catalyst particle size affects both the catalyst activity and light driven charge separation efficiency, their effects on the overall photocatalytic efficiency are still poorly understood. Herein, using a well-defined semiconductor–metal heterostructure with Pt nanoparticle catalysts selectively grown on the apexes of CdS nanorods (NRs), we study the effect of the Pt catalyst size on light driven H 2 generation quantum efficiency (QE H 2 ). With the increase of the Pt catalyst size from 0.7 ± 0.3 to 3.0 ± 0.8 nm, the QEH 2 of CdS–Pt increases from 0.5 ± 0.2% to 38.3 ± 5.1%, by nearly 2 orders of magnitude. Transient absorption spectroscopy measurement reveals that the electron transfer rate from the CdS NR to the Pt tip increases with the Pt diameter following a scaling law of d 5.6 , giving rise to the increase of electron transfer efficiency at larger Pt sizes. Additionally, the observed trend can be understood by a simplified kinetic model that assumes the overall efficiency is the product of the quantum efficiencies of charge separation (including hole transfer, electron transfer, and hole scavenging) and water reduction steps, and for CdS–Pt NRs, the quantum efficiencies of electron transfer and water reduction steps increase with the Pt sizes. Our findings suggest the importance of improving the quantum efficiencies of both charge separation and catalysis in designing efficient semiconductor–metal hybrid photocatalysts, especially in the regime of small metal particle sizes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Entropic formation of a thermodynamically stable colloidal quasicrystal with negligible phason strain

Quasicrystals have been discovered in a variety of materials ranging from metals to polymers. Yet, why and how they form is incompletely understood. In situ transmission electron microscopy of alloy quasicrystal formation in metals suggests an error-and-repair mechanism, whereby quasiperiodic crystals grow imperfectly with phason strain present, and only perfect themselves later into a high-quality quasicrystal with negligible phason strain. The growth mechanism has not been investigated for other types of quasicrystals, such as dendrimeric, polymeric, or colloidal quasicrystals. Soft-matter quasicrystals typically result from entropic, rather than energetic, interactions, and are not usually grown (either in laboratories or in silico) into large-volume quasicrystals. Consequently, it is unknown whether soft-matter quasicrystals form with the high degree of structural quality found in metal alloy quasicrystals. Here, we investigate the entropically driven growth of colloidal dodecagonal quasicrystals (DQCs) via computer simulation of systems of hard tetrahedra, which are simple models for anisotropic colloidal particles that form a quasicrystal. Using a pattern recognition algorithm applied to particle trajectories during DQC growth, we analyze phason strain to follow the evolution of quasiperiodic order. As in alloys, we observe high structural quality; DQCs with low phason strain crystallize directly from the melt and only require minimal further reduction of phason strain. Finally, we also observe transformation from a denser approximant to the DQC via continuous phason strain relaxation. Our results demonstrate that soft-matter quasicrystals dominated by entropy can be thermodynamically stable and grown with high structural quality––just like their alloy quasicrystal counterparts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Modeling of Cr 2 O 3 -doped UO 2 Creep and Fracture (Rev.1)

The pursuit of UO 2 fuel with various dopants (e.g., Cr 2 O 3 and Al 2 O 3 ) for improved accident tolerance due to larger grain sizes and suspected reductions in fission gas release and creep requires advanced modeling and simulation tools. Such tools enable rapid multiscale development and enhanced understanding of material behavior in regimes for which experimental measurements may be lacking. This work builds on previous multiscale modeling efforts targeting the fission gas behavior of Cr 2 O 3 -doped UO 2 , in order to begin exploring the effects of dopants on mechanical properties, with a particular focus on creep and fracture. The results indicate that dislocation climb is the most dominate creep mechanism at almost all temperatures, and sensitivity analyses provide evidence in further support of this claim. The limited experimental data generated by the U.S. Department of Energy’s Advanced Fuels Campaign (AFC) on the tensile strength of Cr 2 O 3 -doped UO 2 indicate that fracture may be less severe in doped specimens. The results on creep and fracture are preliminary, given the large degree of uncertainty associated with the models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulation and Experimental Validation of an Integrated Heat Pump – Thermal Energy Storage Using a Room-Temperature Phase Change Material

As the dependence on electrical heat pumps (HPs) and intermittent renewables increases, grid strains are expected to grow. This necessitates an energy storage system to reduce the mismatch between energy supply and demand. Thus, a proposed dual-mode commercially available 14.1 kW HP was integrated with a single 22°C phase change material (PCM) thermal storage system (TES) to load-shift both cooling and heating loads. The HP-TES system was manufactured and experimentally tested using a novel test matrix based on AHRI 210/240 psychrometric conditions. Furthermore, transient dual-mode system-level HP-TES models were developed in Modelica and validated using the experimental test conditions. Base HP cooling and heating experimental tests at ambient temperatures of 35°C and −8.3°C show that the modified HP-TES maintained the rated system capacity and performance. The HP-TES discharge provided approximately 30% and 50% reductions in cooling and heating demand, respectively. The transient HP-TES models predicted system capacity and total power input for discharge and recharge operating modes within ±4% mean percentage error, and recharge power input within ±2%, with maximum errors occurring at the equipment startup. During system operation, the sources of model deviations are first-order polynomial fits of the PCM digital scanning calorimetry (DSC) data and unaccounted supercooling in the PCM during solidification. Nonetheless, the model predictions agree with the experimental tests, demonstrating the availability of robust, accurate, and validated transient models that can be used for further validation and the development of system controls.

25 ENERGY STORAGE↗

The Installation of Direct Water-Cooling Systems to Reduce Cooling Energy Requirements for High-Performance Computing Centers

A large cluster of High-Performance Computing (HPC) equipment at the Lawrence Livermore National Laboratory in California was retrofitted with an Asetek cooling system. The Asetek system is a hybrid scheme with water-cooled cold plates on high-heat-producing components in the information technology (IT) equipment, and with the remainder of the heat being removed by conventional air-cooling systems. In order to determine energy savings of the Asetek system, data were gathered and analyzed two ways: using top-down statistical models, and bottom-up engineering models. The cluster, “Cabernet”, rejected its heat into a facilities cooling water loop which in turn rejected the heat into the same chilled water system serving the computer-room air handlers (CRAHs) that provided the air-based cooling for the room. Because the “before” and “after” cases both reject their heat into the chilled water system, the only savings is due to reduction in CRAH fan power. The top-down analysis showed a 4% overall energy savings for the data center (power usage effectiveness (PUE) —the ratio of total data center energy to IT energy— dropped from 1.60 to 1.53, lower is better); the bottom-up analysis showed a 3% overall energy savings (PUE from 1.70 to 1.66) and an 11% savings for the Cabernet system by itself (partial PUE of 1.51). Greater savings, on the order of 15-20%, would be possible if the chilled water system was not used for rejecting the heat from the Asetek system. About 37% of the heat from the Cab system was rejected to the cooling water, lower than at other installations.

Earni, Shankar↗

Quantitative Account of the Bonding Properties of a Rubredoxin Model Complex [Fe(SCH3)4]q, q = -2, -1, +2, +3

Iron-sulfur clusters play important roles in biology as parts of electron transfer chains and catalytic cofactors. Here, we report a detailed computational analysis of a structural model of the simplest natural iron-sulfur cluster of rubredoxin and its cationic counterparts. Specifically, we report results for the ground and low-lying electronically excited states of the complex [Fe(SCH3)4]2-/1-/2+/3+, using Multi-Reference (CASSCF, MRCISD), and Coupled Cluster [CCSD(T)] methodology in order to provide accurate adiabatic reduction energies, dissociation energies and insights into the bonding analysis. The nature of the Fe-S chemical bond and the magnitude of the ionization potentials in the anionic and cationic [Fe(SCH3)4] complexes offer a physical rationale for the relative stabilization, structure and speciation of these complexes. Anionic and cationic complexes present different types of chemical bonds: prevalently ionic in [Fe(SCH3)4]2-/1- complexes and covalent in [Fe(SCH3)4]2+/3+ complexes. The ionic bonds result in an energy gain for the transition [Fe(SCH3)4]2-®[Fe(SCH3)4]- (i.e., FeII®FeIII) of 1.5 eV, while the covalent bonds result in an energy loss for the transition [Fe(SCH3)4]2+®[Fe(SCH3)4]3+ of 16.6 eV, almost half of the IP of Fe2+. The ionic vs covalent bond character influences the Fe-S bond strength and length, i.e., ionic Fe-S bonds are longer than covalent ones by about 0.2 Å (for FeII) and 0.04 Å (for FeII). Finally, the average Fe-S heterolytic bond strength is 6.7 eV (FeII) and 14.6 (FeIII) eV at the RCCSD(T) level of theory.

Tzeli, Demeter↗

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗

Examining transport and integrated modeling predictive capabilities for negative-triangularity scenarios

This paper investigates the predictive capabilities of TGYRO and TGLF models in assessing the performance of negative triangularity (NT) plasmas compared to positive triangularity (PT) plasmas in fusion devices. TGYRO predicts kinetic profiles, while TGLF analyzes turbulent transport. The study reveals that TGYRO reasonably predicts NT profiles similar to PT, although it overpredicts the high-power scenarios where there is increased experimental MHD activity. TGLF analysis finds reduced linear growth rates in NT and altered flux spectra relative to PT. Additionally, the TGLF SAT0 saturation model is observed to predict high-k transport and a reduction of particle transport with the electron temperature gradient. These findings are further corroborated by core-pedestal modeling using the Stability Transport Equilibrium Pedestal workflow, showing stronger confinement improvements in NT, particularly at higher power densities for the SAT0 saturation model. Furthermore, the study underscores the importance of accurately capturing turbulence saturation mechanisms for NT in order to project its performance accurately in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

System Modeling of the HTTR and Economic Dispatch Model of the Secondary System

High Temperature Gas-cooled Reactors (HTGRs) can be used for the generation of electricity and their process heat can be used to improve the efficiency of chemical processes such as hydrogen production. The JAEA-operated High-Temperature engineering Test Reactor (HTTR-GT/H2) is exploring using the reactor for electricity and hydrogen production. A RELAP5-3D model of the HTTR-GT/H2 secondary system has been developed using design information. The various components and heat exchangers in the secondary system were modeled and results were compared to the design conditions. The results for the sole-power generation mode were shown to fit the design conditions very well. The largest temperature difference was on the order of 7 K, and the largest pressure difference was on the order of 0.05 MPa. The results for the hydrogen cogeneration mode did not match the design conditions nearly as well. The largest temperature difference was about 39 K and the largest pressure difference was about 0.27 MPa at the compressor outlet. The larger differences for the hydrogen cogeneration mode are attributed to the various complex components and the flow being split in the secondary loop. A transient reduction in heat removal capability of the secondary system was investigated. Reactor temperatures are anticipated to rise as a result. The core reactivity response due to this increase in temperature is investigated and is expected to add negative reactivity to the reactor. An economic dispatch model was developed for a nuclear-driven iodine-sulfur cycle system to determine hydrogen sale prices that would make such a system profitable. The study focuses on the development of the economic model and the role that input data plays on final calculated values. It was found that the input electricity prices, whether using historical data or a host of synthetic time histories, produce significantly different breakeven hydrogen sale prices. As such, great care should be used in these economic dispatch analyses to select reasonable input assumptions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Scheduler Modeling of Distributed Energy Resources for Providing Ancillary Services

Distribution energy resources (DERs) have been integral components of modern power systems, and their capability to provide grid services has been widely studied. To promote the deployment of these resources in providing grid services in real-world utility operations, this paper proposes a day-ahead scheduler model for a distribution system connected DER plant. A certain amount of generation capacity of this DER plant is reserved for frequency services, and some ancillary services for the distribution system-including peak load reduction, voltage regulation, and power factor control-are integrated into the model. The model is tested on a real-world distribution system. From the simulation results, the energy and reserve schedule of the solar photovoltaic (PV) unit and battery energy storage system (BESS) can be determined, and voltage and power factor are well maintained. Additionally, in order to demonstrate the specific characteristics of the co-located and hybrid operation modes for the PV and BESS, these two modes are analyzed both theoretically and through real-time simulation. Simulation results show that most of PV's variability is transferred to the net power in the co-located mode, whereas it is transferred to the BESS in the hybrid mode. This proposed scheduler model and the comparison of co-located and hybrid modes can provide practical guidance for the applications of DER plant in the real-world utility.

14 SOLAR ENERGY↗

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shape restricted additive hazards models: Monotone, unimodal, and U‐shape hazard functions

We consider estimation of the semiparametric additive hazards model with an unspecified baseline hazard function where the effect of a continuous covariate has a specific shape but otherwise unspecified. Such estimation is particularly useful for a unimodal hazard function, where the hazard is monotone increasing and monotone decreasing with an unknown mode. A popular approach of the proportional hazards model is limited in such setting due to the complicated structure of the partial likelihood. Our model defines a quadratic loss function, and its simple structure allows a global Hessian matrix that does not involve parameters. Thus, once the global Hessian matrix is computed, a standard quadratic programming method can be applicable by profiling all possible locations of the mode. However, the quadratic programming method may be inefficient to handle a large global Hessian matrix in the profiling algorithm due to a large dimensionality, where the dimension of the global Hessian matrix and number of hypothetical modes are the same order as the sample size. We propose the quadratic pool adjacent violators algorithm to reduce computational costs. The proposed algorithm is extended to the model with a time‐dependent covariate with monotone or U‐shape hazard function. In simulation studies, our proposed method improves computational speed compared to the quadratic programming method, with bias and mean square error reductions. We analyze data from a recent cardiovascular study.

Mathematical & Computational Biology↗

Computational Modeling of Multi-Pass Rolling Parameters Effect on Resulting Fuel Foil Shape

A focus of the U.S. Department of Energy is to improve production yield and reduce the cost of Low Enriched Uranium (LEU)-molybdenum alloy (U-10Mo) monolithic fuel plates that will be replacing High Enriched Uranium (HEU) oxide dispersion fuels used currently in the United States High Performance Research Reactors (USHPRR). One area of improvement is lowering the high transverse waviness and longitudinal waviness currently present within rolled foils prior to cladding to produce fuel plates. Traditional rolling manufacturing techniques for other metal foils use winders to pull and straighten the foil as it is rolled back and forth to the final thickness. This approach cannot be used to roll thin U-10Mo foils (0.008-0.025” thick) because only small castings can be rolled due to nuclear criticality safety concerns. As a result, the fuel foils are too short (1 m in length) to use traditional winders. Therefore, it is crucial to identify other rolling parameters (i.e., roller friction, axial tension load, roller diameter, and roll pass reduction percent) that might reduce transverse waviness and longitudinal waviness in the rolled fuel foil and develop a high-yield, low-cost multi-pass rolling manufacturing process. This report documents a systematic finite element modeling study to investigate the effects of numerous rolling parameters to reduce resulting transverse waviness and longitudinal waviness in the fuel foil during multi-pass rolling of U-10Mo foils. The rolling of a U-10Mo plate with initial dimensions of 1”x1”x 0.048” is modeled using Abaqus CAE. This rolling is modeled to undergo eight 20% reduction roll passes to a final fuel foil thickness of 0.01”. The elastic-plastic constitutive model of the U-10Mo alloy was input to the fuel foil rolling model. The rollers were modeled as rigid bodies. A comparison of rolling friction coefficients of 0.3 and 0.7 over a wide range of applied axial tension loads were investigated in order to evaluate the effect of using a lubricant during rolling. The effect of roller diameter on the resulting transverse waviness and longitudinal waviness of the fuel foil over a wide range of axial tension loads were also investigated by modeling rollers 7/8” and 3.75” in diameter. The results of this systematic finite element method study will aid manufacturers in producing low transverse waviness and reduced longitudinal waviness in U-10Mo fuel foils.

U-10Mo, FEA, Rolling, Residual Stress, Fuel Foil↗

Optimized Tandem Catalyst Patterning for CO 2 Reduction Flow Reactors

Tandem catalysis involves two or more catalysts arranged in proximity within a single reaction vessel, with the aim of synergistically aligning the catalysts’ reaction pathways to maximize overall system performance. This study presents a proof of concept showing the integration of continuum transport modeling with design optimization in a simplified two-dimensional flow reactor setup for electrochemical CO 2 reduction. Ag catalysts provide the CO 2 ⟶ CO reaction capability, and Cu catalysts provide the CO ⟶ high-value products reaction capability. Given a set of input parameters, the optimization algorithm uses adjoint methods to modify the Ag/Cu surface patterning in order to maximize the current density toward high-value products, such as ethylene. The optimized designs yield significant performance enhancement especially at more negative applied voltages (i.e., stronger surface reactions) and for larger numbers of patterning sections. For an applied voltage of −1.7 V vs. SHE, the 12-section optimized design increases the current density toward ethylene by up to 65% compared to the unoptimized 2-section design. For the optimized cases, observed differences in the production and consumption of CO (the key intermediate species) and minimized zones of low CO reactant surface concentration on Cu sections explain the improved reactor performance.

CO2 reduction↗

Grounding our Understanding of the Impacts of Boreal Forest Expansion on Shallow Cumulus Clouds with a Simple Modeling Framework

Abstract The expansion of the boreal forest poleward is a potentially important driver of feedbacks between the land surface and Arctic climate. A growing body of work has highlighted the importance of differences in evaporative resistance between different possible future Arctic land covers, which in turn alters humidity and cloudiness in the boundary layer, for these feedbacks. While thus far this problem has been studied primarily with complex Earth system models, we turn to a locally focused, idealized model capable of diagnosing and testing the sensitivity of first-order processes connecting vegetation, the atmospheric boundary layer, and low clouds in this critical region. This allows us to benchmark the mechanisms and results at the center of predictions from larger-scale simulations. A surface dominated by broadleaf trees, characterized by higher albedo and lower surface evaporative resistance, drives cooling and moistening of the boundary layer relative to a surface of needleleaf trees, characterized by lower albedo and higher surface evaporative resistance. Differences in evaporative resistance between these hypothetical Arctic vegetation covers are of equal importance to changes in albedo for the initial response of the boundary layer to boreal expansion, even with our idealized approach. However, compensation between the elevation of the lifting condensation level (LCL) and more rapid growth of the mixed layer over higher evaporative resistance surfaces can minimize changes in the favorability of shallow clouds over different land cover types under some conditions. We then perform two tests on the sensitivity of this compensating effect, to changes in water availability, represented first by a reduction in boundary layer humidity and then by both a reduction in humidity and soil moisture available to our vegetation surface. Finally, given the importance of this potential LCL–mixed-layer height compensation in our idealized modeling results, we look to determine its relevance in observational data from a field campaign in boreal Finland. These observations do confirm that such a coupling plays an important role in cumulus-topped boundary layers over a needleleaf forest surface. While our results confirm some underlying mechanisms at the center of prior work with Earth system models, they also provide motivation for future work to constrain the impact of boreal forest expansion. This will include both large eddy simulations to examine the impact of processes and feedbacks not resolved by a mixed-layer model, as well as a more systematic evaluation and comparison of relevant observations at the site in Finland and sites from prior boreal field campaigns. Significance Statement Clouds and vegetation are both important components of the climate system that interact across a range of scales. These interactions are central to understanding how changes at the land surface feedback on climate. For example, if a forest expands or recedes, diagnosing how that will impact clouds will determine whether you predict warming or cooling temperatures from that shift in the forest area. These predictions are often made with complex Earth system models, but we look to a more idealized representation of the land–atmosphere system to diagnose how shallow clouds should respond to changes in surface properties with different scenarios of boreal forest expansion at a more foundational level. This both grounds our understanding of previous analysis and provides helpful direction for future studies of this relevant and impactful land cover change.

Meteorology & Atmospheric Sciences↗

Slow manifold reduction for plasma science

The classical Chapman–Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand breakdown of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham–Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation n e = Z i n i .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗