Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and 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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Structure of an Atmospheric River over Australia and the Southern Ocean. Part I: Tropical and Midlatitude Water Vapor Fluxes

An atmospheric river (AR) impacting Australia-Tasmania and the Southern Ocean during the Austral summer on 28-29 January 2018 during the SOCRATES campaign is analyzed using both a modeling and observational approach. Gulfstream-V dropsonde measurements and GPM radar analyses were used in conjunction with Weather Research and Forecasting model simulations with water vapor tracers to investigate the relative contributions of tropical and mid-latitude moisture sources to the AR. Moisture associated with a monsoonal tropical depression became entrained into a mid-latitude frontal system that extended to 60°S over the Southern Ocean reaching the associated low pressure system 850 km off the east coast of Antarctica – effectively connecting the tropics and the polar region. Tropical moisture contributed to about 50% of the precipitable water within the AR. The tropical contribution to precipitation decreased with latitude, from > 70% over the Australian continent, to ~50% off the Australian coast, to less than 10% poleward of 60°S. Precipitation in the AR conformed to a seeder-feeder model, the seeder precipitation forming primarily from tropical sourced moisture above the altitude of the 0°C isotherm, and the feeder precipitation sourced from mid-latitude sourced moisture below the altitude of the 0°C isotherm. The integrated vapor transport (IVT) through the core of the AR reached above 500 kg m-1s-1 during 1200 UTC 28 January to 0600 UTC 29 January 29, 1.29 times the average amount of water carried by the world’s largest river, the Amazon. The high IVT strength might be attributed to the higher water vapor content associated with the warmer temperatures in Austral summer.

Rauber, Robert M.↗
DOE OSTI2024

Towards accelerating particle-resolved direct numerical simulation with neural operators

In this paper, we present our ongoing work aimed at accelerating a particle-resolved direct numerical simulation model designed to study aerosol–cloud–turbulence interactions. The dynamical model consists of two main components—a set of fluid dynamics equations for air velocity, temperature, and humidity, coupled with a set of equations for particle (i.e., cloud droplet) tracing. Rather than attempting to replace the original numerical solution method in its entirety with a machine learning (ML) method, we consider developing a hybrid approach. We exploit the potential of neural operator learning to yield fast and accurate surrogate models and, in this study, develop such surrogates for the velocity and vorticity fields. We discuss results from numerical experiments designed to assess the performance of ML architectures under consideration as well as their suitability for capturing the behavior of relevant dynamical systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Biomass char gasification kinetic rates compared to data, including ash effects

This investigation provides both experimental and theoretical analyses of biomass char gasification kinetics over 1150–1350 °C in the environment of CO 2 , H 2 O and combinations of these reactants. The experimental data include continuous measurements of particle mass, external and internal temperatures, shape, and size with periodic but not continuous measurements of char porosity, pore size distribution, and pore surface area. The model describes heterogeneous char gasification for poplar wood, corn stover and switchgrass in different shapes and industrially relevant sizes (12.5 mm and smaller). Nonlinear least-squares regression produces optimized power-law model parameters for global kinetics that describe biomass char gasification with respect to both CO 2 and H 2 O separately and in combination. A single set of kinetic parameters describes reaction rates for all three fuels. Model simulations agree with measured data at all stages of char conversion even though the observed mass loss rates change significantly with conversion. In conclusion, the change in rate with conversion depends strongly on ash content, and this investigation provides a simple and theoretically based reaction rate expression that includes the effect of inert ash occupying an increasing portion of the reactive surface with increased conversion.

09 BIOMASS FUELS↗
DOE OSTI2022

Produced Water and Waste Heat-aided Blowdown Water Treatment: Using Chemical and Energy Synergisms for Value Creation

The project objective was to develop a cooling blowdown water (BDW) treatment process utilizing produced water (PW) and low-grade heat to maximize water reuse and saleable by-product generation while reducing chemical and energy footprints of the treatment. The proposed treatment process consists of mixing, softening, organics and suspended solids removal, reverse osmosis (RO), thermal desalination, and brine electrolysis. BDW samples collected from a local coal-fired power plant and PW samples from two shale gas production wells were used in this study. Each treatment unit was first designed and tested to quantify its treatment efficiency, and its chemical and energy requirements. In addition, a process model was developed and model simulations were conducted based on the experimental results and literature data to optimize the treatment process. A techno-economic analysis was conducted to quantify chemical and energy savings as well as production of 10-lb brine as a saleable product. With the field-collected BDW and PW samples, mixing experiments determined a volumetric mixing ratio 10:1 (BDW:PW) resulted in the best performance of multivalent ions removal and largest chemical savings for softening. Softening of the BDW/PW mixtures using alkaline chemicals (Na 2 CO 3 and NaOH) achieved 95%-100% removal of scaling-forming cations (Ca, Mg, Fe, Ba, Sr) and 60% of silicon, and 10% of total organic carbon (TOC). The mixing and softening treatments yielded an effluent with total dissolved solids (TDS) concentration of 23 g/L. Activated carbon (AC) filtration removed TOC to a low level (< 3 mg/L) and further removed remaining scale-forming divalent metals and silica from the softened water. The AC filtration resulted in a slight reduction of TDS from 23 g/L to 20 g/L, leaving behind only mostly monovalent ions (i.e., sodium and chloride) in the filtered water. These pretreatments yielded a feed water that met the criteria of the downstream reverse osmosis (RO) to prevent membrane fouling. A cross-flow RO system was used to further concentrate the TDS of the AC effluent. Various factors including TDS, pH, and applied pressure were examined and optimal conditions were determined for the co-treatment process. An integrated process consisting of mixing, softening, AC filtration and RO was used to treat a continuous flow (0.25 – 1.2 L/min, or 0.07 – 0.32 gpm) and successfully generated RO permeate as product water (TDS < 0.5 g/L) for reuse in cooling operation, and a concentrate (TDS ~ 45 g/L) to be further treated in a thermal desalination unit. These flow rates meet the FOA’s criterion of 0.01 – 1 gpm. Overall, the co-treatment of BDW/PW allowed shorter ramp-up time compared to treatment of BDW alone. It resulted in 40% and 55% savings of Na 2 CO 3(s) and NaOH, respectively, compared to treating the BDW and PW individually for the same level of softening. The co-treatment also resulted in a 29% energy saving compared to treatment of BDW only for the level of TDS concentration. A thermal desalination system was designed using CFD simulations and manufactured in the WVU Innovation Hub for further treatment of the RO concentrate to generate 10-lb brine. The system has a design flow rate of 2 gpm and has been successfully tested. A bench-scale brine electrolysis system was developed for on-site generation of chlorine/hypochlorite (Cl 2 /OCl - ) and caustic soda (NaOH) as useful chemicals for the co-treatment process. Using salt solutions (0.5 M and 1 M), the system achieved faradaic efficiencies of 93%-97% and 70%-77% for caustic soda and chlorine/hypochlorite generation, respectively. An economic analysis showed that the electricity costs for on-site generation of these chemicals were significantly lower than the chemical prices offered by suppliers. An industrial-scale process model consisting of mixing, softening, AC filtration, RO, thermal desalination, and brine electrolysis was developed using the Aspen Plus V9 in conjunction with Aspen Custom Modeler V9. The model serves as a solvable Aspen Plus model and as basis to form the costing infrastructure. In addition, techno-economic analysis considering capital, operating, and transportation costs was conducted. An optimization solution showed that produced water for mixing is still advantageous in low quantities. The optimum solution approaches a leveled cost of water (LCW) of 2 $/m 3 which becomes cost competitive with nominal water treatment prices.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

ION PRODUCTION AND MITIGATION IN DC HIGH-VOLTAGE PHOTO-GUNS

One of the biggest obstacles to operating a GaAs polarized electron source with a long charge lifetime is the mitigation of ion back-bombardment. Several techniques exist to either clear ions from the accelerator or to mitigate ion damage of the photocathode. Predicting the effectiveness of these techniques requires sophisticated simulation models of electron impact ionization within the photo-gun. In this work, the effectiveness of applying a positive anode bias voltage to mitigate ion damage and increase the charge lifetime of the GaAs photocath?ode was studied over three run periods at the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab). The charge lifetime with the biased anode configuration was 1.80±0.51 times longer than the lifetime of the usual grounded anode configuration. Sim?ulations of ionization within the CEBAF photo-gun and adjacent beamline were made using General Particle Tracer (GPT) and a new C++ custom element to predict and explain the substantial improvement in charge lifetime. The experimental results and the development of the ionization custom element, with its use in simulations of ion back-bombardment with the biased anode, are described in detail.

Yoskowitz, Joshua↗
DOE OSTI2022

An Evolutionary Algorithm to Personalize Stool-Based Colorectal Cancer Screening

Fecal immunochemical testing (FIT) is an established method for colorectal cancer (CRC) screening. Measured FIT-concentrations are associated with both present and future risk of CRC, and may be used for personalized screening. However, evaluation of personalized screening is computationally challenging. In this study, a broadly applicable algorithm is presented to efficiently optimize personalized screening policies that prescribe screening intervals and FIT-cutoffs, based on age and FIT-history. We present a mathematical framework for personalized screening policies and a bi-objective evolutionary algorithm that identifies policies with minimal costs and maximal health benefits. The algorithm is combined with an established microsimulation model (MISCAN-Colon), to accurately estimate the costs and benefits of generated policies, without restrictive Markov assumptions. The performance of the algorithm is demonstrated in three experiments. In Experiment 1, a relatively small benchmark problem, the optimal policies were known. The algorithm approached the maximum feasible benefits with a relative difference of 0.007%. Experiment 2 optimized both intervals and cutoffs, Experiment 3 optimized cutoffs only. Optimal policies in both experiments are unknown. Compared to policies recently evaluated for the USPSTF, personalized screening increased health benefits up to 14 and 4.3%, for Experiments 2 and 3, respectively, without adding costs. Generated policies have several features concordant with current screening recommendations. The method presented in this paper is flexible and capable of optimizing personalized screening policies evaluated with computationally-intensive but established simulation models. It can be used to inform screening policies for CRC or other diseases. For CRC, more debate is needed on what features a policy needs to exhibit to make it suitable for implementation in practice.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

ComStock/COVID-19: U.S. Commercial Building Stock Analysis of COVID-19 Mitigation Strategies

ComStock is a DOE platform, with NREL as the lead developer, that models the U.S. commercial building stock using OpenStudio/EnergyPlus to perform physics-based building energy modeling simulations. The ASHRAE Epidemic Taskforce (ASHRAE-ETF) has released several COVID-19 mitigation strategies for commercial buildings. ComStock was used to analyze the impact of four strategies on the existing commercial building stock. These strategies are upgrading to MERV-13 filters, disabling demand control ventilation (DCV), increasing minimum outdoor air to 100%, and HVAC flushing mode operation. Results show the impact of each measure on energy consumption, electric consumption, gas consumption and peak electric.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Faster recovery of North Atlantic tropical cyclone-induced cold wakes in recent decades

Intense winds associated with tropical cyclones (TCs) generate surface ocean cooling in their wakes, which can persist for several weeks in their aftermath. While multi-decadal observations of the sea surface have shown a substantial warming of the ocean, long-term changes in cold wake recovery time remain largely unknown. Here we find a trend toward faster recovery of TC cold wakes in the Atlantic main development region (MDR) since 2001. This is due primarily to a decrease in the strength of the North Atlantic trade winds, which reduces evaporative cooling of the ocean. The faster damping of TC cold wakes has led to a significant increase in the intensification of subsequent TCs that encounter lingering wakes from prior TCs, with a magnitude that is about 9% of that from long-term warming of the ocean. Finally, earth system model simulations indicate that the observed decrease in the cold wake recovery time will likely continue into the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year -1 (1Tg = 10 12 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH 4 year -1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of -34 to -46 Tg CH 4 year -1 , we estimated that the global net land methane emissions are 149–176 Tg CH 4 year -1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

58 GEOSCIENCES↗
DOE OSTI2022

Orographic enhancement of rainfall over the Congo Basin

Abstract The Congo rainforest located in central equatorial Africa is an important, yet understudied part of the globe surrounded by complex orographic features. A primitive understanding of precipitation processes such as mesoscale convective dynamics magnifies uncertainties in the future climate projections of the hydrological cycle over the Congo. Furthermore, the effects of orography, which is an important forcing for convection and precipitation, are poorly resolved by climate models, and ill‐conceptualized over the Congo. To address this knowledge gap, perturbed orographic forcing experiments are conducted using the Weather Research and Forecasting (WRF) mesoscale numerical model in a high‐resolution convection‐permitting model setup. The model simulated selected dates in November 2014. The thunderstorms and rainfall simulated in the control run for the case study analyzed in this article compared reasonably well to satellite‐derived brightness temperature and rainfall data. The results from this case study show that the dynamical impact of increasing the height of the East African Highlands is the blocking of the lower‐tropospheric tropical easterlies. This weakening of the lower‐tropospheric zonal winds increases the windshear over the Congo Basin resulting in slower propagating, more intense mesoscale convective systems with enhanced rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Persisting volcanic ash particles impact stratospheric SO 2 lifetime and aerosol optical properties

Volcanic ash is often neglected in climate simulations because ash particles are assumed to have a short atmospheric lifetime, and to not participate in sulfur chemistry. After the Mt. Kelut eruption in 2014, stratospheric ash-rich aerosols were observed for months. Here we show that the persistence of super-micron ash is consistent with a density near 0.5 g cm -3 , close to pumice. Ash-rich particles dominate the volcanic cloud optical properties for the first 60 days. We also find that the initial SO 2 lifetime is determined by SO 2 uptake on ash, rather than by reaction with OH as commonly assumed. About 43% more volcanic sulfur is removed from the stratosphere in 2 months with the SO 2 heterogeneous chemistry on ash particles than without. This research suggests the need for re-evaluation of factors controlling SO 2 lifetime in climate model simulations, and of the impact of volcanic ash on stratospheric chemistry and radiation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗
DOE OSTI2021

Techno‐economic feasibility analysis of engineered energycane‐based biorefinery co‐producing biodiesel and ethanol

Abstract High feedstock cost and low oil yields per unit of land from temperate oilseed crops limit the growth of commercial‐scale biodiesel production. Recently, highly productive crops, such as sugarcane and energycane, have been engineered to accumulate triacylglycerides (TAGs) that allow the production of far more industrial vegetable oil than previously possible. A proof‐of‐concept suggests that biodiesel production from engineered energycane will be possible. However, before making efforts for scale‐up, it is critical to understand the commercial feasibility and economic competitiveness of this process. This study performs techno‐economic analysis of a unique biorefinery processing energycane to co‐produce biodiesel and ethanol. Comprehensive process simulation models were developed for two scenarios: (i) biodiesel from TAGs and ethanol from fermentation of sugars in juice and (ii) biodiesel from TAGs and ethanol from fermentation of sugars in juice and hydrolysis of carbohydrates in bagasse. Based on the target levels, the analysis was performed for energycane containing 0%, 5%, and 7.7% TAGs (d.b.). The biodiesel from engineered energycane was found economically viable and competitive to soybean biodiesel. Although the capital investment is higher compared to the soybean biodiesel plant, the biodiesel production costs ($0.66–$0.9/L) were lower than soybean biodiesel ($0.91/L). Biorefinery‐scenario‐1 processing energycane containing 7.7% TAG produces biodiesel with profitability (IRR 7.84) slightly lower than soybean biodiesel (IRR 8.3), but yields five times of biodiesel per unit land and is self‐sustainable for energy requirements. The surplus electricity can displace fossil electricity and provide environmental benefits. Monte Carlo simulation indicated that biorefinery is profitable with a 29%–65% probability (NPV > 0) which is largely controlled by feedstock composition and biodiesel market price. It is important to note that energycane can be grown on the marginal rainfed lands in S.E. USA, where soybean would not be viable. Biodiesel from engineered energycane would therefore be complementary to soydiesel in the United States.

09 BIOMASS FUELS↗
DOE OSTI2024

Model associated with: "Thermodynamic control on the decomposition of organic matter across different electron acceptors"

This model data package is associated with the publication “Thermodynamic control on the decomposition of organic matter across different electron acceptors” submitted to Soil Biology and Biochemistry (Zheng et al., 2023; https://doi.org/10.1016/j.soilbio.2024.109364).In this research, a thermodynamic modeling framework is built to flexibly incorporate both organic matter (OM) molecules and electron acceptors for estimating potential free energy release from various redox reactions and to further predict reaction rates based on Microbial Transition State Theory. The model package includes scripts for thermodynamic modeling and postprocessing. Input Fourier-transform ion cyclotron resonance (FTICR) data are from a previous experimental study (Boye et al., 2018), and model outputs are free energy predictions and stoichiometric coefficients associated with all possible redox reactions.This data package is associated with the project GitHub repository found at MM_bioenergetic_modeling.This data package contains four folders (Input_FTICR, Model, Output, and Output_processing), a file-level metadata (FLMD) csv, and a data dictionary (dd) csv. Please see Zheng_bioenergetic_modeling_flmd.csv for a list of all files contained in this data package and descriptions for each. The Zheng_bioenergetic_modeling_dd.csv file describes the csv column headers. The “Model” folder contains scripts to run energy balance calculations for each electron acceptor. The “Output” folder contains csv files with stoichiometric information from model simulations. And the "Output_processing" folder contains scripts for reaction rate calculations and to generate plots.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Control System of Multi-Port Autonomous Reconfigurable Solar Power Plant (MARS) & HIL Platforms for Design

Multi-port autonomous reconfigurable solar power plant (MARS) provides an attractive alternative to connect photovoltaic (PV) and energy storage systems (ESSs) to high-voltage direct current (HVdc) links and high-voltage alternating current (ac) grids. In this paper, a unique hierarchical control system of MARS is proposed and evaluated. To evaluate the control system and associated algorithms in early-stage research of complex architectures like MARS, it is important to develop unique suitable hardware-in-the-loop (HIL) platforms. In this paper, the HIL platforms for MARS to evaluate the performance of the hierarchical control system and the control algorithms implemented are presented. Further, they help with the design process of control systems. The real-time simulation models and algorithms that are utilized for MARS in the HIL platforms are also discussed in the paper. The HIL experiments of the control system of MARS showcase the capability to provide continuity of operation under faults and frequency support to the power grid during loss of generation. They also showcase the stability of the proposed hierarchical control system of MARS.

14 SOLAR ENERGY↗
DOE OSTI2023

Response of Extreme Rainfall to Atmospheric Warming and Wetting: Implications for Hydrologic Designs Under a Changing Climate

Understanding the processes of rainfall extremes and their response to anthropogenic climate change is pivotal for improved adaptation of unprecedented flood hazards around the world. Here we take the record-breaking 20 July 2021 storm over central China as an example. We investigate the response of this particular storm to atmospheric warming (i.e., increase in air temperature) and wetting (i.e., increase in atmospheric moisture content) based on a series of convection-permitting model simulations. Our results show non-monotonic changes of the space-time rainfall variability to either increased temperature or atmospheric moisture content. The most extreme rain rate is produced when relative humidity is increased by 20~40 % or temperature is increased by less than 2°C. The non-monotonic rainfall response is more clearly revealed at fine spatial (100-1000 km 2 ) and temporal scales (less than 6 hours) rather than over the entire domain (~10 4 km 2 ) and aggregated over the storm duration (around two days). This is mainly attributable to the distinct feedbacks from atmospheric dynamics (i.e., moisture convergence and interaction with regional topography) rather than regulated by thermodynamic changes alone. Furthermore, atmospheric warming poses notable changes in the vertical structure of storm cells, contributing to reduced areal reduction factors at small spatial scales and short durations, while atmospheric wetting additionally modifies storm evolution properties. Our modeling analyses challenge the existing practices for hydrologic designs under a changing climate, highlighting particular vulnerability for cities or small basins to short-duration rainfall extremes and the resultant flash flood hazards.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The LAKE model input dataset for three Arctic lakes

This dataset contains meteorological data collected for three Arctic lakes and compiled to satisfy input requirements of the LAKE 2.0 model. The dataset was generated to act as a benchmarking dataset for future model-data inter-comparisons. The LAKE 2.0 model simulates temperatures within the water later and the sedimentary layer of a lake. The LAKE2.0. is an open-source code and available to download via this weblike http://tesla.parallel.ru/Viktor/LAKE/-/wikis/LAKE-model (last visit July 14, 2021). The meteorological data are required to simulate the surface energy balance at the surface of a lake. This dataset includes a compilation of the meteorological data pulled from multiple data streams, including National Oceanic and Atmospheric Administration (NOAA) climate data, Circumarctic Lakes Observation Network (CALON) data, and the United States Geological Survey (USGS) data. The data were compiled for three Arctic lakes: FoxDen (66.55877, -164.45670), Atqasuk (70.452497, -156.951984), and Toolik (68.63150, -149.60740). Each meteorological data is in comma-delimited format (file extension ‘.dat’) and includes eight columns: Temperature [K], Pressure [Pa], longwave downward radiation [W/m2], shortwave downward radiation [W/m2], “U” wind speed [m/s], ”V” wind speed [m/s], humidity [kg/kg], precipitation [m/s]. In addition to the meteorological data file, we included setup and driver files. The Toolik lake is the deepest out of three lakes and has inflowing and outflowing groundwater data. InflowOutflowREADME.txt has more information about inflow and outflow flies. The other two lakes are much shallower and modeled as a closed system (i.e. no water inflow or outflow).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of Thermal Scattering Law and Cross Sections of Hydrogen in Paraffin Wax

Paraffin wax is frequently used as a neutron moderator and shielding material. The main component of paraffin wax is straight-chain alkanes (n-alkanes). The deposition of paraffin wax is primarily attributed to the crystallization of n-alkanes. It is important to gain a deeper understanding of the mechanisms underlying the behavior of paraffin wax, which would impact the thermal scattering Law (TSL) and cross sections and affect the analysis of neutronic and critical systems. In this work, a classical molecular dynamics (CMD) simulation model was used in LAMMPS to create the TSL and cross sections at room temperature and pressure. To generate the required velocity autocorrelation functions (VACF), previously published data were used to validate the approach and models by comparing them with key model parameters. The phonon density of state (DOS) was calculated using Fourier transformation of the normalized VACF. This DOS was used as the primary input to estimate the TSL (S($a$, $β$)) and cross sections of hydrogen in paraffin wax. The TSL and cross sections of hydrogen were estimated using the Full Law Analysis Scattering System Hub (FLASSH) code. The cross section of hydrogen in paraffin wax is consistent with other hydrocarbon materials such as polyethylene with deviations due to structure in the lowest energy region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗