Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional Distribution”

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,189 records · Page 66

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Review of Dew Point Evaporative Cooling Technology for Air Conditioning Applications

Indirect evaporative cooling has the potential to significantly improve the natural environment. It follows from a significant reduction in electricity consumption in the hot period, and hence lower operating costs for cooling systems. This paper presents the current state of knowledge and research directions on dew point indirect evaporative cooling. It was found that researchers focus on the development of dew point indirect evaporative coolers (DPIEC) by improving its design, geometry, water distribution, and new porous materials implementation. To evaluate the performance of new types of DPIEC, different methods are used by the scientists. Finally, optimized devices are studied in terms of their performance in different systems, like hybrid and desiccant systems, considering different climate conditions. Potential directions of development of evaporative technologies were indicated, such as increasing the coefficient of performance of solid desiccant evaporative cooling systems, developing novel geometry, and efficient water distribution, including development of porous materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Validation of two-layer model for underexpanded hydrogen jets

Previous studies have shown that the two-layer model more accurately predicts hydrogen dispersion than the conventional notional nozzle models without significantly increasing the computational expense. However, the model was only validated for predicting the concentration distribution and has not been adequately validated for predicting the velocity distributions. In the present study, particle imaging velocimetry (PIV) was used to measure the velocity field of an underexpanded hydrogen jet released at 10 bar from a 1.5 mm diameter orifice. The two-layer model was the used to calculate the inlet conditions for a two-dimensional axisymmetric CFD model to simulate the hydrogen jet downstream of the Mach disk. The predicted velocity spreading and centerline decay rates agreed well with the PIV measurements. The predicted concentration distribution was consistent with data from previous planar Rayleigh scattering measurements used to verify the concentration distribution predictions in an earlier study. The jet spreading was also simulated using several widely used notional nozzle models combined with the integral plume model for comparison. These results show that the velocity and concentration distributions are both better predicted by the two-layer model than the notional nozzle models to complement previous studies verifying only the predicted concentration profiles. Thus, this study shows that the two-layer model can accurately predict the jet velocity distributions as well as the concentration distributions as verified earlier. Though more validation studies are needed to improve confidence in the model and increase the range of validity, the present work indicates that the two-layer model is a promising tool for fast, accurate predictions of the flow fields of underexpanded hydrogen jets.

08 HYDROGEN↗

PacWave Grid Integration Study: Transient and Dynamic Conditions (Final Report)

This report describes the results of PSCAD simulations that were performed in 2020 to assess the impacts of PacWave generation on CLPUD's 12.47-kV distribution and 69-kV subtransmission systems. PacWave South (PacWave) is a wave energy test facility planned by Oregon State University. PacWave shore facilities will be located south of Seal Rock on the Oregon coast. PacWave is expected to be operational in 2022 and will connect up to 10 MW of generation to CLPUD's Seal Rock distribution feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Aerosol responses to precipitation along North American air trajectories arriving at Bermuda

North American pollution outflow is ubiquitous over the western North Atlantic Ocean, especially in winter, making this location a suitable natural laboratory for investigating the impact of precipitation on aerosol particles along air mass trajectories. We take advantage of observational data collected at Bermuda to seasonally assess the sensitivity of aerosol mass concentrations and volume size distributions to accumulated precipitation along trajectories (APT). The mass concentration of particulate matter with aerodynamic diameter less than 2.5 µm normalized by the enhancement of carbon monoxide above background (PM 2.5 /ΔCO) at Bermuda was used to estimate the degree of aerosol loss during transport to Bermuda. Results for December–February (DJF) show that most trajectories come from North America and have the highest APTs, resulting in a significant reduction (by 53 %) in PM 2.5 /ΔCO under high-APT conditions (> 13.5 mm) relative to low-APT conditions (< 0.9 mm). Moreover, PM 2.5 /ΔCO was most sensitive to increases in APT up to 5 mm (–0.044 µg m –3 ppbv –1 mm –1 ) and less sensitive to increases in APT over 5 mm. While anthropogenic PM 2.5 constituents (e.g., black carbon, sulfate, organic carbon) decrease with high APT, sea salt, in contrast, was comparable between high- and low-APT conditions owing to enhanced local wind and sea salt emissions in high-APT conditions. The greater sensitivity of the fine-mode volume concentrations (versus coarse mode) to wet scavenging is evident from AErosol RObotic NETwork (AERONET) volume size distribution data. A combination of GEOS-Chem model simulations of the 210 Pb submicron aerosol tracer and its gaseous precursor 222 Rn reveals that (i) surface aerosol particles at Bermuda are most impacted by wet scavenging in winter and spring (due to large-scale precipitation) with a maximum in March, whereas convective scavenging plays a substantial role in summer; and (ii) North American 222 Rn tracer emissions contribute most to surface 210 Pb concentrations at Bermuda in winter (~75 %–80 %), indicating that air masses arriving at Bermuda experience large-scale precipitation scavenging while traveling from North America. A case study flight from the ACTIVATE field campaign on 22 February 2020 reveals a significant reduction in aerosol number and volume concentrations during air mass transport off the US East Coast associated with increased cloud fraction and precipitation. These results highlight the sensitivity of remote marine boundary layer aerosol characteristics to precipitation along trajectories, especially when the air mass source is continental outflow from polluted regions like the US East Coast.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable Predictive Control and Optimization for Grid Integration of Large-Scale Distributed Energy Resources: Preprint

Integration of a large number of distributed energy resources (DERs) into the power grid needs a scalable power balancing method. We formulate the power balancing problem as a look-ahead optimization problem to be solved sequentially by a power distribution system aggregator based on a model predictive control (MPC) framework. Solving large-scale look-ahead control problem requires proper configuration of the control steps. In this paper, to solve large-scale control problems, we propose a variable time granularity where control time steps nearby the current control step have finer resolutions. The aggregator objective includes maximization of power production revenue and minimization of power purchasing expense, renewable power curtailment, and mileage costs for energy storage and electric vehicle (EV) charging stations while satisfying system capacity and operational constraints. The control problem is formulated as a mixed-integer linear program (MILP) and solved using the XpressMP solver. We perform simulations considering a copper plate representation of a large distribution network consisting of 2507 devices (controllable DERs) including curtailable photovoltaics (PVs), energy storage batteries, EV charging stations, and buildings with heating, ventilation, and air conditioning units (HVACs). We show the effectiveness of the proposed approach in managing DERs interactively for maximum energy trading profit and local supply-demand power balancing. Finally, we demonstrate that the proposed method outperformed other benchmark controllers regarding computation time without compromising operational performance.

DER↗

Cellulose nanofibrils and nanocrystals in confined flow: Single-particle dynamics to collective alignment revealed through scanning small-angle x-ray scattering and numerical simulations

Nanostructured materials made through flow-assisted assembly of proteinaceous or polymeric nanosized fibrillar building blocks are promising contenders for a family of high-performance biocompatible materials in a wide variety of applications. Optimization of these processes relies on improving our knowledge of the physical mechanisms from nano- to macroscale and especially understanding the alignment of elongated nanoparticles in flows. Here, we study the full projected orientation distributions of cellulose nanocrystals (CNCs) and nanofibrils (CNFs) in confined flow using scanning microbeam SAXS. For CNCs, we further compare with a simulated system of dilute Brownian ellipsoids, which agrees well at dilute concentrations. However, increasing CNC concentration to a semidilute regime results in locally arranged domains called tactoids, which aid in aligning the CNC at low shear rates, but limit alignment at higher rates. Similarly, shear alignment of CNF at semidilute conditions is also limited owing to probable bundle or flock formation of the highly entangled nanofibrils. Finally, this work provides a quantitative comparison of full projected orientation distributions of elongated nanoparticles in confined flow and provides an important stepping stone towards predicting and controlling processes to create nanostructured materials on an industrial scale.

36 MATERIALS SCIENCE↗

Particle Size Distribution Dynamics Can Help Constrain the Phase State of Secondary Organic Aerosol

Particle phase state is a property of atmospheric aerosols that has important implications for the formation, evolution, and gas/particle partitioning of secondary organic aerosol (SOA). In this work, we use a size-resolved chemistry and microphysics model (SOM-TOMAS), updated to include an explicit treatment of particle phase state, to constrain the bulk diffusion coefficient (Db) of SOA produced from ??- pinene ozonolysis. By leveraging data from laboratory experiments performed in the absence of a seed and under dry conditions, we find that the Db for SOA can be constrained (1-7 ×10-15 cm2 s-1 in these experiments) by simultaneously reproducing the time-varying SOA mass concentrations and the evolution of the particle size distribution. Another version of our model that used the predicted SOA composition to calculate the glass transition temperature, viscosity, and, ultimately, Db (~10-15 cm2 s-1) of the SOA was able to reproduce the mass and size distribution measurements when we included oligomer formation (oligomers accounted for about a fifth of the SOA mass). Our work highlights the potential of a size resolved SOA model to constrain the particle phase state of SOA by utilizing historical measurements of the evolution of the particle size distribution.

He, Yicong↗

Dominant Wave Energy Systems and Conditional Wave Resource Characterization for Coastal Waters of the United States

Opportunities and constraints for wave energy conversion technologies and projects are evaluated by identifying and characterizing the dominant wave energy systems for United States (US) coastal waters using marginal and joint distributions of the wave energy in terms of the peak period, wave direction, and month. These distributions are computed using partitioned wave parameters generated from a 30 year WaveWatch III model hindcast, and regionally averaged to identify the dominant wave systems contributing to the total annual available energy ( A A E ) for eleven distinct US wave energy climate regions. These dominant wave systems are linked to the wind systems driving their generation and propagation. In addition, conditional resource parameters characterizing peak period spread, directional spread, and seasonal variability, which consider dependencies of the peak period, direction, and month, are introduced to augment characterization methods recommended by international standards. These conditional resource parameters reveal information that supports project planning, conceptual design, and operation and maintenance. The present study shows that wave energy resources for the United States are dominated by long-period North Pacific swells (Alaska, West Coast, Hawaii), short-period trade winds and nor’easter swells (East Coast, Puerto Rico), and wind seas (Gulf of Mexico). Seasonality, peak period spread, and directional spread of these dominant wave systems are characterized to assess regional opportunities and constraints for wave energy conversion technologies targeting the dominant wave systems.

16 TIDAL AND WAVE POWER↗

Powertrain distribution based on total cost of ownership for privately owned vehicles and TNC in the Chicago metropolitan area

From a Total Cost of Ownership (TCO) standpoint, the lowest-cost powertrain option depends on driving conditions and number of miles driven over the lifetime of the vehicle. Studies looking at the value proposition of powertrain options traditionally assess the energy consumption based on standardized drive cycles. TCO is then calculated using average yearly Vehicle Miles Travelled (VMT). This approach which looks at TCO at the aggregate level, does not take into consideration the variations in driving patterns and the spread of VMT that exist in the vehicle population. In this study, we first determine the driving pattern of each vehicle in the Chicago metropolitan area. We then consider different powertrain options for each vehicle, calculate its associated TCO and provide a powertrain distribution based on lowest TCO. Vehicles under consideration are Privately Owned Vehicles (POV) and vehicles used by Transportation Network Companies (TNC).

Freyermuth, Vincent↗

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗

Experimental and Theoretical Evaluation of Feed Flow Collar Design for Shell Fed Hollow Fiber Membrane Modules

An experimental and theoretical study of module collar design is presented here. Hollow fiber membranes are prepared by dip coating a poly(vinylidene) (PVDF) support with a polydimethylsiloxane (PDMS) gutter layer and a Pebax 2533 selective layer. Fiber bundles with a well-defined fiber packing are prepared using a 3D printed module. A parallel fiber bundle consisting of 4-9 uniformly spaced fibers is created with printed tabs that align the fibers and create a tubesheet. The tabs are sealed within a printed case that possesses a series of external ports for gas introduction and removal. Uniquely, both port location and the use of a collar to assist fluid distribution in the shell can be varied for the same fiber bundle. Experimental measurements are compared to computational fluid dynamics (CFD) simulations. The experimental module design allows high-fidelity representation of the fiber bundle and module case in the simulations. Comparisons between experiment and simulation are in good agreement over a broad range of experimental conditions. The detrimental effect of having ports located too close, leading to stagnation regions, is captured as well as the beneficial effects of using a collar for shell-side fluid distribution around the fiber bundle. Such results help validate the use of CFD to develop high-performance module designs.

Tran, Thien↗

Coupled wave-current modeling for hydrodynamic load analysis of macroalgae cultivation farms

Large-scale cultivation of macroalgae is one of the most promising biofuel sources that could reduce our consumption of fossil fuels and would be particularly economically competitive when grown for the co-production of additional goods such as food and textiles. Current production of biofuel from algal biomass is limited by labor costs and lack of data on potential impacts of cultivation and harvesting on marine and coastal environments as well as a comprehensive characterization of the impact of environmental conditions on macroalgal feedstock and their bioproduct and biofuel yield and quality. Awarded by the United States (U.S.) Department of Energy Advanced Research Projects Agency-Energy (ARPA-E) Macroalgae Research Inspiring Novel Energy Resources (MARINER) program, which seeks to enable the U.S. as a global leader in the production of marine biomass, the Marine Biological Laboratory (MBL) has been developing a test system for tropical seaweed cultivation in the Gulf of Mexico and the Caribbean. An integral step in designing macroalgae cultivation farms is to adequately characterize the hydrodynamic climate at potential growth sites. In regions like the Gulf of Mexico and the Caribbean, which are highly exposed to tropical cyclones, it is critical to not only consider the day-to-day hydrodynamic conditions but to also assess the risk of extreme sea states to ensure the survivability of future macroalgae farms. Under these considerations, the Pacific Northwest National Laboratory (PNNL) is providing modeling support for MBL to inform the design and siting of the farm systems that have been proposed for their ARPA-E MARINER project. This report describes the development of a high-resolution coupled storm surge and wave model to simulate the hydrodynamics and wave climate at proposed macroalgae cultivation sites selected by MBL in Florida and Puerto Rico. Model results, including model validation, water level, current distributions, and sea states, are discussed for both selected sites in Florida and Puerto Rico coast. These simulations provide an accurate insight into the hydrodynamic conditions that a macroalgae farm is likely to experience during its operational lifetime, including current information to support fine-scale hydrodynamic load modeling, risk analysis and system design.

09 BIOMASS FUELS↗

Full time-dependent SOLPS-ITER simulation of the SPARC tokamak: actuator design for particle and divertor condition control *

Abstract This paper presents the application of full time-dependent SOLPS-ITER simulations for actuator design in the SPARC tokamak. This study employs both the EIRENE module, a neutral solver, and the B2.5 plasma module in a time-dependent mode. This is in contrast to most SOLPS simulations, which focus on steady-state solutions, where the neutral distribution is evolved without any time limit or for a time step of 1 ⋅ 10 − 3 second, which is several orders of magnitude larger than the fluid plasma time step. The time-dependent EIRENE was tested with a fixed B2.5 background and compared with a simple conductance based model in a simplified pump chamber geometry. This comparison aimed to verify the reliability of the neutral relaxation timescale derived from the time-dependent EIRENE. Subsequently, a full time-dependent simulation was performed in a realistic geometry, with the Monte-Carlo neutral time step synchronized with the plasma fluid time step. The numerical setup of the code, including relative time steps and the size of the census data used to store Monte-Carlo particle information is considered. The full-time dependent simulations are then applied to inform the design of the SPARC louver structure, which affects divertor plasma parameters by regulating the neutral conductance from the divertor to the pump. The response of the plasma and neutral parameters was captured on a timescale that enables the design of the actuator to consider time-dependent control capability. It was found that changing the louver opacity has an equivalent effect as varying the gas throughput via puff actuation. Therefore, equivalent divertor plasma conditions can be obtained from both actuators, while the neutral pressure distribution in the pump and divertor differs for each actuator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Habitat specialization and edge effects of soil microbial communities in a fragmented landscape

Abstract Soil microorganisms play outsized roles in nutrient cycling, plant health, and climate regulation. Despite their importance, we have a limited understanding of how soil microbes are affected by habitat fragmentation, including their responses to conditions at fragment edges, or “edge effects.” To understand the responses of soil communities to edge effects, we analyzed the distributions of soil bacteria, archaea, and fungi in an experimentally fragmented system of open patches embedded within a forest matrix. In addition, we identified taxa that consistently differed among patch, edge, or matrix habitats (“specialists”) and taxa that showed no habitat preference (“nonspecialists”). We hypothesized that microbial community turnover would be most pronounced at the edge between habitats. We also hypothesized that specialist fungi would be more likely to be mycorrhizal than nonspecialist fungi because mycorrhizae should be affected more by different plant hosts among habitats, whereas specialist prokaryotes would have smaller genomes (indicating reduced metabolic versatility) and be less likely to be able to sporulate than nonspecialist prokaryotes. Across all replicate sites, the matrix and patch soils harbored distinct microbial communities. However, sites where the contrasts in vegetation and pH between the patch and matrix were most pronounced exhibited larger differences between patch and matrix communities and tended to have edge communities that differed from those in the patch and forest. There were similar numbers of patch and matrix specialists, but very few edge specialist taxa. Acidobacteria and ectomycorrhizae were more likely to be forest specialists, while Chloroflexi, Ascomycota, and Glomeromycota (i.e., arbuscular mycorrhizae) were more likely to be patch specialists. Contrary to our hypotheses, nonspecialist bacteria were not more likely than specialist bacteria to have larger genomes or to be spore‐formers. We found partial support for our mycorrhizal hypothesis: arbuscular mycorrhizae, but not ectomycorrhizae, were more likely to be specialists. Overall, our results indicate that soil microbial communities are sensitive to edges, but not all taxa are equally affected, with arbuscular mycorrhizae in particular showing a strong response to habitat edges. In the context of increasing habitat fragmentation worldwide, our results can help inform efforts to maintain the structure and functioning of the soil microbiome.

Winfrey, Claire C. [Department of Ecology and Evol↗

Temperature response and brine availability to heated boreholes in bedded salt

Abstract There is a growing need for disposal of high‐level nuclear waste. To reduce uncertainty associated with brine availability to repository excavations in salt formations, a collaboration between Sandia, Los Alamos, and Lawrence Berkeley National Laboratories is performing a series of borehole‐scale coupled process tests. Here, we report on the first round of the Brine Availability Test in Salt (BATS) project, a “shakedown” experiment called Phase 1s. Experimental testing included placing a resistive heater, a 260‐W radiative heater, and a 750‐W radiative heater within previously drilled horizontal boreholes at the Waste Isolation Pilot Plant (WIPP) while monitoring temperature and water inflow. The experiments successfully achieved the targeted temperature of 120 °C when using the 750‐W radiative heater. Simulations using FEHM (Finite Element Heat and Mass transfer code) and TOUGH‐FLAC (Transport Of Unsaturated Groundwater and Heat–Fast Lagrangian Analysis of Continua) were able to accurately predict the coupled thermo–hydro–mechanical–chemical response of salt, matching the observed temperature and brine production. Due to the extremely low permeability of salt, these systems take many years to reach steady state when perturbed by mining activities. Long‐term numerical simulations are used to develop the initial pressure and saturation conditions. The inclusion of a damaged rock zone with higher permeability around the borehole also affects the saturation and pressure distributions and plays an important role in dissipating the potential for thermal pressurization. Knowledge gained from this round of experimentation and simulation will be used to conduct the next BATS project experiment in new boreholes at WIPP.

58 GEOSCIENCES↗

Agent-based simulation and child protection systems: Rationale, implementation, and verification

Simulation models are an important tool used in health care and other disciplines to support operational research and decision-making. In the child protection literature, simulation models are an under-utilized source of research evidence. Here, in this paper, we describe the rationale for and the development of an agent-based simulation of a child protection system in the US. Using the investigation, prevention service, and placement histories of 600,000 children served in an urban child welfare system, we walk the reader through the development of a prototype known as OSPEDALE. The governing equations built into OSPEDALE probabilistically simulate the onset of investigations. Then, drawing from empirical survival distributions, the governing equations trace the probability of subsequent interactions with the system (recurrence of maltreatment, service referrals, and placement) conditional on the characteristics of children, their assessed risk level, and prior child protection system involvement. As an initial test of OSPEDALE's utility, we compare empirical admission counts with counts generated from OSPEDALE. Though the verification step is admittedly simple, the comparison shows that OSPEDALE replicates the empirical count of new admissions closely enough to justify further investment in OSPEDALE. Management of public child protection systems is increasingly research evidence-dependent. The emphasis on research evidence as a decision-support tool has elevated evidence acquired through randomized clinical trials. Though important, the evidence from clinical trials represents only one type of research evidence. Properly specified, simulation models are another source of evidence with real-world relevance.

60 APPLIED LIFE SCIENCES↗

Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations

Here, we present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/R, where R is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.

Monte Carlo↗