Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure modeling”

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

Catalytic Autoxidation for Depolymerization of Multilayer Plastic Films

Recycling multilayer plastic films is challenged by a diversity of polymers, prompting development of new recycling methods. For the depolymerization of mixed polymers like those in multilayer films, metal-catalyzed autoxidation offers a versatile chemical recycling method to deconstruct multiple polymers to useful oxygenates. Here, we demonstrate that catalytic autoxidation is effective for depolymerizing multilayer films across diverse chemistries. We investigated conditions for a model polyethylene substrate using a Co, Mn, and Br cocatalyst system, achieving full carbon closure with oxygenated small molecules contributing up to 48 mol% carbon. Subsequently, we characterized product distributions for several common polymers used in multilayer films using high-resolution mass spectrometry (HRMS) and developed analytical methods to quantify the resulting complex product streams. Optimized conditions for polyethylene were applied to 11 multilayer plastic films containing 10 different polymers, including films with nonpolymeric potential disrupters like aluminum foil and titanium dioxide, showing that catalytic autoxidation is effective across a broad range of polymer types and is resistant to disrupters and additives. The generation of CO 2 in these reactions overall suggests that both reaction engineering and modifications to the reaction conditions will be required to achieve higher yields of soluble oxygenated products.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nonadiabatic Dynamics of Photoexcited cis-Stilbene Using Ab Initio Multiple Spawning

The photochemistry of cis-stilbene proceeds through two pathways: cis-trans isomerization and ring closure to 4a,4b-dihydrophenanthrene (DHP). Despite serving for many decades as a model system for photoisomerization, the photodynamics of cis-stilbene is still not fully understood. In this work, we use ab initio multiple spawning on a SA-2-CASSCF(2,2) potential energy surface to simulate the nonadiabatic dynamics of isolated cis-stilbene. We find the cyclization (to DHP and cis-stilbene) and isomerization (to trans- and cis-stilbene) reaction coordinates to be orthogonal; branching between the two pathways is determined on the S1 excited state within 150 fs of photoexcitation. Furthermore ,trajectory basis functions (TBFs) undergoing cyclization decay rapidly to the ground state in 250 fs, while TBFs moving along the isomerization coordinate remain on the excited state longer, with the majority decaying between 300 and 500 fs. We observe three avoided crossing regions in the dynamics: two along the isomerization coordinate (displaying pyramidalization and migration of an ethylenic hydrogen or phenyl group), and one DHP-like conical intersection along the cyclization coordinate. The isomeric form of the vibrationally hot photoproducts (as determined by measurement 2 ps after photoexcitation) is determined within less than 50 fs of decay to the ground state mediated by passage through a conical intersection. Excess vibrational energy of ground state cis- and trans-stilbene is channelled into phenyl torsions (with mostly opposing directionality). Finally, our simulations are validated by direct comparison to experiment for the absorption spectrum, branching ratio of the three photoproducts (44:52:4 cis-stilbene:trans-stilbene:DHP), and excited state lifetime (520 ± 40 fs).

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

Feasibility Study for Disposal Control Rod Assemblies Using UNF-ST&DARDS As-Loaded Zion Dual Purpose Cask Models

This report documents an initial evaluation to support future use of disposal control rod assemblies (DCRAs) for post-closure criticality control in dual purpose canisters (DPCs). The work described herein is an extension of previous efforts performed by Walker using inputs generated by the Used Nuclear Fuel-Storage, Transportation & Disposal Analysis Resource and Data System (UNFST&DARDS) for the Zion site DPCs with as-loaded isotopic compositions. The results of this analysis demonstrate that there are multiple pathways to using DCRA for post-closure criticality control. Various configurations of DCRA material, diameter, number of rods per DCRA, and number and location of DCRAs within a DPC were shown to be effective in varying degrees for the set of DPCs analyzed by Walker. Because of the variations in DPC as-loaded isotopic compositions considered in an array of DCRA parametric sweeps, it can be concluded that a DPC-specific methodology is feasible (i.e., a one-size-fits-all approach may not be needed). Instead, the utility program created for this work can be expanded to develop capabilities to provide DPC-specific DCRA arrangements to limit cost and weight and to allow for operational considerations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Nuclide Inventory Benchmark for BWR Spent Nuclear Fuel: Challenges in Evaluation of Modeling Data Assumptions and Uncertainties

This work discusses challenges and approaches to uncertainty analyses associated with the development of a nuclide inventory benchmark for fuel irradiated in a boiling water reactor. The benchmark under consideration is being developed based on experimental data from the SFCOMPO international database. The focus herein is on how to address missing data in fuel design and operating conditions that are important for adequately simulating the time-dependent changes in fuel during irradiation in the reactor. The effects of modeling assumptions and uncertainties in modeling parameters on the calculated nuclide inventory were analyzed and quantified through computational models developed using capabilities in the SCALE code system. Particular attention was given to the impact of the power history and water coolant density on the calculated nuclide inventory, as well as to the effect of geometry modeling considerations not usually addressed in a nuclide inventory benchmark. These considerations include gap closure, channel bow, and channel corner radius, which do not usually apply to regular reactor operation but are relevant for assessing impacts of potential anomalous operating scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

On the Existence of the Plateau Emission in High-energy Gamma-Ray Burst Light Curves Observed by Fermi-LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) shows long-lasting high-energy emission in many gamma-ray bursts (GRBs), similar to X-ray afterglows observed by the Neil Gehrels Swift Observatory (Swift). Some LAT light curves (LCs) show a late-time flattening reminiscent of X-ray plateaus. We explore the presence of plateaus in LAT temporally extended emission analyzing GRBs from the second Fermi-LAT GRB Catalog from 2008 to 2016 May with known redshifts, and check whether they follow closure relations corresponding to four distinct astrophysical environments predicted by the external forward shock model. We find that three LCs can be fit by the same phenomenological model used to fit X-ray plateaus and show tentative evidence for the existence of plateaus in their high-energy extended emission. The most favorable scenario is a slow-cooling regime, whereas the preferred density profile for each GRBs varies from a constant-density interstellar medium to an r –2 wind environment. We also compare the end time of the plateaus in γ-rays and X-rays using a statistical comparison with 222 Swift GRBs with plateaus and known redshifts from 2005 January to 2019 August. Within this comparison, the case of GRB 090510 shows an indication of chromaticity at the end time of the plateau. Finally, we update the 3D fundamental plane relation among the rest-frame end time of the plateau, its correspondent luminosity, and the peak prompt luminosity for 222 GRBs observed by Swift. We find that these three LAT GRBs follow this relation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Aerosol-induced closure of marine cloud cells: enhanced effects in the presence of precipitation

Abstract. The Weather Research Forecasting (WRF) version 4.3 model is configured within a Lagrangian framework to quantify the impact of aerosols on evolving cloud fields. Kilometer-scale simulations utilizing meteorological boundary conditions are based on 10 case study days offering diverse meteorology during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA). Measurements from aircraft, the ground-based Atmosphere Radiation Measurement (ARM) site at Graciosa Island in the Azores, and A-Train and geostationary satellites are utilized for validation, demonstrating good agreement with the WRF-simulated cloud and aerosol properties. Higher aerosol concentration leads to suppressed drizzle and increased cloud water content in all case study days. These changes lead to larger radiative cooling rates at cloud top, enhanced vertical velocity variance, and increased vertical and horizontal wind speed near the base of the lower-tropospheric inversion. As a result, marine cloud cell area expands, narrowing the gap between shallow clouds and increasing cloud optical thickness, liquid water content, and the top-of-atmosphere outgoing shortwave flux. While similar aerosol effects are observed in lightly to non-raining clouds, they tend to be smaller by comparison. These simulations show a relationship between cloud cell area expansion and the radiative adjustments caused by liquid water path and cloud fraction changes. The adjustments are positive and scale as 74 % and 51 %, respectively, relative to the Twomey effect. While higher-resolution large-eddy simulations may provide improved representation of cloud-top mixing processes, these results emphasize the importance of addressing mesoscale cloud-state transitions in the quantification of aerosol radiative forcing that cannot be attained from traditional climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Full-Scale Retrieval and Transfer System Mockup for Pneumatic Conveyance of High-Level Calcine Waste - 19563

The Idaho National Laboratory reprocessed spent nuclear fuel from 1953 to 1994. The reprocessing used solvent extraction to separate uranium, resulting in a high-level radioactive liquid waste. A calcination process was used to convert the liquid waste to a high-level radioactive granular solid, referred to as calcine. The calcine is stored in stainless-steel bins contained within six concrete vaults called the Calcined Solids Storage Facilities (CSSFs). The Calcine Disposition Project is tasked with demonstrating the ability to retrieve the calcine from the bins of CSSF 1 and transferring it to the bins of CSSF 6. CSSF 1 was constructed to be a final calcine repository, with no forethought of retrieval. It represents the greatest technical hurdle to calcine extraction of all the CSSFs. The challenges inherent with accessing the calcine, safely extracting it, and securely transporting it are copious and complex. To overcome those challenges in a timely and objective manner, a full-scale mockup will be constructed. The mockup includes a full-size model of a concentric group of bins from CSSF 1 and a storage vessel representing CSSF 6. Mockup equipment will include an air compressor, filtration, blower/vacuum, and retrieval devices. In addition to defining process parameters, the mockup testing scope will demonstrate the amount of residual calcine that can be expected to be in the bins following extraction. This will be key to preparing for closure of CSSF 1 following removal of the calcine. Once testing is complete, the mockup will constitute a model of the calcine extraction and transport system and will be prepared to fulfill a second advantageous objective: the operators that will carry out the calcine transfer from CSSF 1 to CSSF 6 will be trained using the mockup. The control system will be identical, the major pieces of equipment will be representative, and the procedures will be analogous; hence, high-fidelity training can take place at the mockup. Once radiological operations commence at CSSF 1, the mock-up's commission will still not be completed. If any problems arise during calcine transport, the mockup will be ready to serve as a simulator to more effectively reproduce dilemmas and deduce solutions. The mock-up's mission spans the duration of the project. Its usefulness begins with helping to identify process criteria. Operations personnel will benefit greatly from its use as a training tool. And finally, the immediate availability of a simulator to reproduce problems encountered during radiological extraction and transfer will be invaluable. The knowledge gained during this process will also be of continued benefit as future work is done to remove various types of calcine from the remaining CSSFs. The methodology of safe and secure pneumatic transfer of radioactive material may be of benefit to others throughout the U.S. Department of Energy complex. (authors)

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

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗
DOE OSTI2021

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗
DOE OSTI2020

Safety Functions and Features, Events and Processes for the E-Area Performance Assessment

The DOE Technical Standard, “Disposal Authorization Statement and Tank Closure Documentation,” (DOE 2017) recommends the use of safety functions and features, events and processes (FEPs) to support development of conceptual models and identification of scenarios to be considered in a performance assessment (PA). The FEP process provides a means to describe how a PA considers and addresses the factors that could influence the performance of key barriers. Understanding the roles of barriers in terms of limiting migration helps to focus on how changes in the system could lead to a situation where those roles cannot be fulfilled and there is the potential for compromised performance. The FEPs screening and review process was used to identify FEPs that are relevant for the EArea Low-Level Waste Facility (LLWF) and specifically those FEPs that could have a detrimental impact on the effectiveness of a given safety function. For this PA, a default list of FEPs developed at the International Atomic Energy Agency (IAEA 2004) and an approach implemented for PAs at the Hanford and Idaho sites (Mehta et al. 2016, DOE-ID 2019) are used to identify processes and events that could influence the effectiveness of a given safety function for the E-Area LLWF (e.g., subsidence can impact the safety function of the cover system and lead to increased infiltration). The Hanford and Idaho PAs represent two of the most recent applications of this approach. The PA evaluates the potential impacts of changes in performance of different features of the system and demonstrates that the safety functions represent multiple and redundant barriers. Barrier analyses, assuming a safety function is not present, also test the robustness of the system in the event of the loss of one or more safety functions. Such evaluations also support a qualitative illustration of the concept of defense in depth. The safety concept for closure of the E-Area LLWF (generically referred to as “E-Area”) encompasses a variety of different features (i.e., administrative controls, natural site features, and engineered barriers) that reduce the potential impacts on human health and the environment from the residual waste that will remain after closure. These features can be represented as a collection of safety functions acting independently and as a system to provide for overall safety. In some applications, there have been attempts to assign numerical expectations to specific safety functions, but that is not the intent in this case. The concept of safety functions is used more qualitatively in two ways for this PA: 1. To illustrate the robustness of the E-Area design, operational practices and closure approach by documenting features that are and are not credited in different modeling cases. 2. To identify the roles of the different features and potential processes and events that could compromise the performance of safety features and need to be considered when developing the modeling approach. This report addresses both safety functions and FEPs for the E-Area PA.

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

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales

Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Relating Convection to GCM Grid-Scale Fields Using Cloud-Resolving Model Simulation of a Squall Line Observed during MC3E Field Experiment

In this study, a WRF (Weather Research and Forecasting) model is used as a cloud-resolving model to simulate a squall line observed on 20 May 2011 in the Southern Great Plains (SGP) of the United States. The model output is then used to examine the relationships between convective precipitation and coarse-grained variables averaged over a range of subdomain sizes equivalent to various global climate model horizontal resolutions. The objective is to determine to what extent convection within the subdomains can be related to these “large-scale” variables, thus that they can potentially serve as closure in convective parameterization. Results show that convective precipitation is well correlated with the vertical velocity at 500 hPa, column integrated moisture convergence and CAPE change due to large-scale advective forcing (dCAPE) for various subdomain sizes, but the correlation decreases with decreasing subdomain size. dCAPE leads convective precipitation for all subdomain sizes examined; however, the lead time decreases with decreasing subdomain size. Moisture convergence leads convective precipitation for subdomain sizes greater than 32 km but has no lead time for smaller subdomain sizes. Mid-tropospheric vertical velocity has no lead time or slightly lags convective precipitation. The lead/lag composite analysis with respect to maximum precipitation time indicates that peaks of large-scale variables increase with decreasing subdomain size. The peaks of 500 hPa vertical velocity and column integrated moisture convergence occur at the same time as maximum precipitation, but maximum dCAPE leads maximum precipitation by twelve minutes.

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

Surprising Charge-Radius Kink in the Sc Isotopes at N = 20

For this work, charge radii of neutron deficient 40 Sc and 41 Sc nuclei were determined using collinear laser spectroscopy. With the new data, the chain of Sc charge radii extends below the neutron magic number N=20 and shows a pronounced kink, generally taken as a signature of a shell closure, but one notably absent in the neighboring Ca, K, and Ar isotopic chains. Theoretical models that explain the trend at N=20 for the Ca isotopes cannot reproduce this puzzling behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Scale Separation Effects on Simulations of Plasma Turbulence

Abstract Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v A / c ) or the relative mass ratio of ions and electrons ( m i / m e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m i / m e ∼ 250 captures the asymptotic behavior of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Performance of Cloud 3D Solvers in Ice Cloud Shortwave Radiation Closure Over the Equatorial Western Pacific Ocean

Abstract For retrieving cloud optical properties from satellite images or computing these properties from climate model output, computationally efficient treatments of cloud horizontal inhomogeneity include the Monte Carlo Independent Column Approximation (McICA) and the Tripleclouds method. Computationally efficient treatment of cloud horizontal radiation exchanges includes the SPeedy Algorithm for Radiative TrAnsfer through CloUd Sides (SPARTACUS). As a test to derive properties from satellite images, we collocate Moderate Resolution Imaging Spectroradiometer (MODIS) cloud retrievals with near‐nadir Cloud and the Earth's Radiant Energy System (CERES) footprints in July 2008 over an equatorial western Pacific Ocean region to compare the performance of the McICA, Tripleclouds, and SPARTACUS solvers to the conventional plane‐parallel homogeneous (PPH) treatment. PPH overestimates cloud albedo, and the three solvers effectively reduce overestimation with root mean square error of shortwave upwelling irradiance decreasing between 15.72 and 18.53 W m −2 , or about 22%–25%. Although cloud top variability does not get fed into the simulations, all three solvers also reduce the effect of cloud top variability on cloud albedo. Entrapment (energy reflected downward from clouds) and horizontal radiation transfer have opposite effects on the SPARTACUS cloud albedo simulation. The net effect depends on the cloud vertical extent, the unawareness of which limits the performance of the SPARTACUS solver.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Plasma thermal transport with a generalized 8-moment distribution function

Moment equations that model plasma transport require an ansatz distribution function to close the system of equations. The resulting transport is sensitive to the specific closure used, and several options have been proposed in the literature. Two different 8-moment distribution functions can be generalized to form a single-parameter family of distribution functions. The transport coefficients resulting from this generalized distribution function can be expressed in terms of this free parameter. This provides the flexibility of matching the 8-moment model to some validating result at a given magnetization value, such as Braginskii’s transport, or the more recent results of Davies et al. [Physics of Plasma, 28, 012305 (2021)]. Here, this process can be thought of as solving for the 8-moment distribution function that matches the value of a transport coefficient given by a Chapman-Enskog expansion, while retaining the improved physical properties, such as finite propagation speeds and time dependence which belong to the hyperbolic moment models. Since the presented generalized distribution function only has a single free parameter, only a single transport coefficient can be matched at a time. However, this generalization process may be extended to provide multiple free parameters. The focus of this Brief Communication is on the dramatically improved thermal conductivity of the proposed model compared to the two base moment models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗