Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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 847 records · Page 47

Constraints on the Physical Properties of GW190814 through Simulations Based on DECam Follow-up Observations by the Dark Energy Survey

On 2019 August 14, the LIGO and Virgo Collaborations detected gravitational waves from a black hole and a 2.6 solar mass compact object, possibly the first neutron star–black hole merger. In search of an optical counterpart, the Dark Energy Survey (DES) obtained deep imaging of the entire 90% confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger. Objects with varying brightness were detected by the DES Pipeline, and we systematically reduced the candidate counterparts through catalog matching, light-curve properties, host-galaxy photometric redshifts, Southern Astrophysical Research spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to full light-curve simulations of supernovae and kilonovae as they would appear in the DECam observations. Because the source class of the merger was uncertain, we utilized an agnostic, three-component kilonova model based on tidally disrupted neutron star (NS) ejecta properties to quantify our detection efficiency of a counterpart if the merger included an NS. We find that, if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than 10 -8.56 , and has a velocity between 0.18c and 0.21c are disfavored at the 2σ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the 4σ level in 95% of future follow-up observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

On the characterization of biases arising from methods and approximations used for sensitivity analyses

The preliminary assessment of the performance and safety of advanced reactors, as well as the identification of research and development needs, rely on computer simulations due to the lack of operational experience. Those simulations should be based on verified and validated computational tools, calculation schemes, and nuclear data libraries, and should be accompanied of a proper characterization of the involved uncertainties, providing confidence in the computational-based predictions. Thus, sensitivity and uncertainty studies, together with integral experiments, play an essential role in that process. Sensitivities are then a critical element and guidelines about how to produce accurate enough sensitivities depending on the subsequent analyses to be performed are of interest. The present study addresses the sensitivities biases arising from the use of different nuclear data libraries, computational methods and the assumption of modelling simplifications. Not only sensitivities for multiplication factor but also for safety-relevant reactivity responses have been analyzed. In the study, propagated uncertainties are also included for the considered parameters to illustrate the impact of the sensitivity profiles obtained for each case on the final uncertainty. Then, relevant recommendations are given for sensitivity analyses of highly complex systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Sensitivity Studies of Premixed Flame Response to Transverse, High Frequency Disturbances

This study analyzes the response of premixed flames to high frequency transverse acoustic disturbances. A significant body of experiment, modeling and simulation efforts has been performed to understand the response of the spatially integrated heat release of premixed flames to low-frequency disturbances. These studies quantified this sensitivity through the Flame Transfer Function (FTF), as the global heat release is the key quantity influencing thermoacoustic stability for acoustically compact flames. Recent efforts have expanded these models to consider high-frequency acoustic disturbances, and the sensitivity of premixed flames through the Rayleigh integral (RI). The full Rayleigh Integral, as opposed to simply the spatially integrated heat release, must be considered in these cases where the flame is not acoustically compact. Prior modeling studies considered several effects - stability of axisymmetric flames centered in a cylindrical combustion chamber and the effect of flame offset from the combustor axis and the effect of mean flame asymmetries in the flame centered in the combustor. These studies showed that the flame's stability was sensitive to all these parameters, as well as to the transverse mode in question. In this paper, we quantify these different parametric sensitivities of premixed flame stability in the presence of high-frequency transverse acoustic modes - specifically sensitivity to flame length and nozzle location within the acoustic mode shape for 2 different transverse modes. The sensitivity of the RI to these control parameters provides an improved understanding of which flame parameters are critical to improve the overall thermoacoustic stability of the system.

Acharya, Vishal S.↗

Towards a Unified Setup to Simulate Mid‐Latitude and Tropical Mesoscale Convective Systems at Kilometer‐Scales

Abstract Mesoscale convective systems (MCSs) are the main source of precipitation in the tropics and parts of the mid‐latitudes and are responsible for high‐impact weather worldwide. Studies showed that deficiencies in simulating mid‐latitude MCSs in state‐of‐the‐art climate models can be alleviated by kilometer‐scale models. However, whether these models can also improve tropical MCSs and whether we can find model settings that perform well in both regions is understudied. We take advantage of high‐quality MCS observations collected over the Atmospheric Radiation Measurement (ARM) facilities in the US Southern Great Plains (SGP) and the Amazon basin near Manaus (MAO) to evaluate a perturbed physics ensemble of simulated MCSs with 4 km horizontal grid spacing. A new model evaluation method is developed that enables to distinguish biases stemming from spatiotemporal displacements of MCSs from biases in their reflectivity and cloud shield. Amazon MCSs are similarly well simulated across these evaluation metrics than SGP MCSs despite the challenges anticipated from weaker large‐scale forcing in the tropics. Generally, SGP MCSs are more sensitive to the choice of model microphysics, while Amazon cases are more sensitive to the planetary boundary layer (PBL) scheme. Although our tested model physics combinations had strengths and weaknesses, combinations that performed well for SGP simulations result in worse results in the Amazon basin and vice versa. However, we identified model settings that perform well at both locations, which include the Thompson and Morrison microphysics coupled with the Yonsei University (YSU) PBL scheme and the Thompson scheme coupled with the Mellor‐Yamada‐Janjic PBL scheme.

54 ENVIRONMENTAL SCIENCES↗

Advancing Multiscale Simulation of Plasma-Surface Interfaces

We report the development of an atomistic-informed, surface-state-dependent predictive model for particle exchange in a carbon-tungsten plasma-surface interface. The predictive model uses machine learning (ML) techniques to learn the energy and angular distributions for particle exchange and rate functions for surface state evolution from molecular dynamics simulations of cumulative bombardment of tungsten by energetic carbon ions. Each predictive component is sensitive to the energy and trajectory of incident plasma species and the surface state. The surface state is represented by a set of surface state descriptors, which were derived from the atomistic surface state for each independent carbon bombardment event. These descriptors are representative of the composition and degree of amorphization of the outermost angstrom of surface material and were chosen to optimize predictive performance for particle exchange at the interface. The distributions for particle exchange (reflection/sputtering) are demonstrated to vary with each surface state descriptor, motivating the development of surface-state-dependent particle exchange models for plasma simulations. The performance of various ML methods was compared, including polynomial quantile regression, artificial neural networks, k-nearest neighbors, and random forest algorithms, with polynomial regression performing the best for interpolation and extrapolation of learned relationships. In addition to the particle exchange model, a neutral network was developed and used to identify data sufficiency throughout surface descriptor space, which will enable real-time feedback during future data production to ensure data is produced where it is most needed, and we provide commentary on improvements to the data production workflow for future endeavors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data-Enabled Predictive Control for Building HVAC Systems

Model predictive control is widely used as a control technology for the computation of optimal control inputs of building heating, ventilating, and air conditioning (HVAC) systems. However, both the benefits and widespread adoption of model predictive control (MPC) are hindered by the effort of model creation, calibration, and accuracy of the predictions. In this paper, we apply the data-enabled predictive control (DeePC) algorithm for designing controls for building HVAC systems. The algorithm solely depends on input/output data from the system to predict future state trajectories without the need for system identification. The algorithm relies on the idea that a vector space of all input–output trajectories of a discrete-time linear time-invariant (LTI) system is spanned by time-shifts of a single measured trajectory, given the input signal is persistently exciting. Closed-loop simulations using EnergyPlus are performed to demonstrate the approach. The simulated building modeled in EnergyPlus is a modified commercial large office prototype building served by an air handling unit-variable air volume HVAC system. Temperature setpoints of zones are used as control variables to minimize the HVAC energy cost of the building considering a time-of-use electricity rate structure. Furthermore, sensitivity analysis is conducted to gain insights into the effect of parameter tuning on DeePC performance. Simulation results are used to illustrate the performance of the algorithm and compare the algorithm with model-based MPC and occupancy-based setpoint controller. Overall, DeePC achieves similar performance compared to MPC for lower engineering effort.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simulation-guided analysis of resonant soft X-ray scattering for determining the microstructure of triblock copolymers

Resonant soft X-ray scattering (RSoXS) probes structure with chemical sensitivity that is useful for determining the morphology of multiblock copolymers. However, the hyperspectral scattering data produced by this technique can be challenging to interpret. Here, in this study, we use computational scattering simulations to extract the microstructure of a model triblock copolymer from the energy-dependent scattering from RSoXS. An ABC triblock terpolymer formed from poly(4-methylcaprolactone) (P4MCL), poly(2,2,2-trifluoroethylacrylate) (PTFEA), and poly (dodecylacrylate) (PDDA), P4MCL-block-PTFEA-block-PDDA, was synthesized as the model triblock system. Through quantitative evaluation of simulated scattering data from a physics-informed set of candidate structure models against experimental RSoXS data, we find the best agreement with hexagonally packed core–shell cylinders. This result is also consistent with electron-density reconstruction from hard X-ray scattering data evaluated against electron-density maps generated with the same model set. These results demonstrate the utility of simulation-guided scattering analysis to study complex microstructures that are challenging to image by microscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Physics-informed machine learning for fault-leakage reduced-order modeling

Geologic carbon storage (GCS) is a promising technology for mitigating CO 2 emissions. The overall success of GCS depends on safe operations that are informed by risk assessment and have proper mitigation plans in place. Performing quantitative probabilistic risk assessment for a GCS site using traditional reservoir simulators can be challenging due to the high computational costs. To overcome this challenge, the US Department of Energy’s National Risk Assessment Partnership (NRAP) project has developed an integrated assessment modeling approach that utilizes computationally efficient reduced-order models (ROM) for simulating various parts of a GCS storage site to quantify uncertainty. Here, in this study, we develop a reduced-order model for fault leakage risk assessment. We use a deep learning approach to build the reduced-order model. We perform a sensitivity analysis and find that the deep learning model yields high accuracy with a much smaller computational cost than full-physics simulation. We also evaluate the performance of the model in scenarios where simulations are not possible to run, providing analysis not previously performed in fault-leakage ROM analyses. Based on a sensitivity analysis of the model, we suggest a simplified conceptual model for fault leakage and site monitoring.

58 GEOSCIENCES↗

Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing

Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.

58 GEOSCIENCES↗

LOCOMOTIVES - Comprehensive Impact and Cost Assessment Framework of Carbon Lowering Approaches for the US Rail Freight System

The goal of this project is to develop a tool to aid railroads and other stakeholders assess and approach the decarbonization of freight rail operations by identifying new, viable low-carbon energy storage and conversion systems for future locomotive systems and how they should be deployed on the existing US freight rail network. In the first quarter, the project focused on collecting data, establishing a simulation workflow, and engaging industry through the creation of the Industry Advisory Board (IAB). In the second quarter, the project focused on selecting fuel pathways and powertrain technologies, setting performance targets, conducting a techno-economic analyses, and developing the simulation framework that would serve as the backbone of the future toolhead. The third quarter involved developing an industry-oriented interactive dashboard powered by a five-step sequential framework, as well as holding industry advisory board meetings as per the initial technology-to-market plan. In the remaining project quarters, the NUFRIEND dashboard were fine-tuned with the help of IAB member feedback and in-depth scenario analyses were conducted to support the techno-economic analysis of energy sources. Additionally, dashboard documentation, project insights, and open-source code on GitHub were prepared and released. Throughout the project, the team completed testing and analysis of all model components, integrated all initial test scenarios, and conducted stakeholder engagement. Lower-carbon drop-in fuels can be deployed as admixtures and are considered uniform across the network at a desired penetration rate, while hydrogen and battery-electric technology deployment poses a more complex problem as they require significant investments to be made in the siting of refueling/charging facilities and the replacement of locomotive fleets. Thus, strategies for locating and sizing refueling/charging facilities on a railroad’s network to meet their energy demands were developed to inform deployment decisions. To address this challenge, the Northwestern University Freight Rail Infrastructure & Energy Network Decarbonization (NUFRIEND) framework presents a five-step sequential framework to select O-D paths, locate facilities, reroute flows, size facilities, and evaluate the deployment for alternative energy sources that require locomotive powertrains to be converted and new refueling infrastructure to be deployed. The NUFRIEND Framework is an industry-oriented tool for simulating the deployment of new energy technologies across the US freight rail network. The framework provides a comprehensive network-level optimization and scenario simulation tool for decarbonizing the freight rail sector, addressing the uncertainties surrounding technological developments by supporting sensitivity analyses for different operational and technological parameters through a transparent and flexible input module. It offers practical alternatives to diesel locomotives and can be applied for any railroad considering the specific network structure and freight demand, outputting evaluation metrics for the associated emissions and costs relative to diesel operations. A number of relevant simulation scenarios were run and analyzed for key insights on the value of different alternative technologies for freight rail decarbonization. The project developments and findings have been presented at numerous conferences and events.

08 HYDROGEN↗

Interpreting radial correlation Doppler reflectometry using gyrokinetic simulations

A linear response, local model for the DBS amplitude applied to gyrokinetic simulations shows that radial correlation Doppler reflectometry measurements (RCDR, Schirmer et al 2007 Plasma Phys. Control. Fusion 49 1019) are not sensitive to the average turbulence radial correlation length, but to a correlation length that depends on the binormal wavenumber $k_{\perp}$ selected by the Doppler backscattering (DBS) signal. Nonlinear gyrokinetic simulations show that the turbulence naturally exhibits a nonseparable power law spectrum in wavenumber space, leading to a power law dependence of the radial correlation length with binormal wavenumber $l_r \sim C k_{\perp}^{-\alpha} (\alpha \approx 1)$ which agrees with the inverse proportionality relationship between the measured $l_r$ and $k_{\perp} $ observed in experiments (Fernández-Marina et al 2014 Nucl. Fusion 54 072001). This new insight indicates that RCDR characterizes the eddy aspect ratio in the perpendicular plane to the magnetic field. It also motivates future use of a nonseparable turbulent spectrum to quantitatively interpret RCDR and potentially other turbulence diagnostics. The radial correlation length is only measurable when the radial resolution at the cutoff location $W_\text n$ satisfies $W_\text n \ll l_r$, while the measurement becomes dominated by $W_\text n$ for $W_\text n \gg l_r$. This suggests that $l_r$ is likely to be inaccessible for electron-scale DBS measurements (${k_{\perp}}{\rho_{\text s}} \gt 1$). The effect of $W_\text n$ on ion-scale radial correlation lengths could be nonnegligible.

radial correlation Doppler reflectometry↗

The Influence of Environment on Post-Detonation Chemistry and Debris Formation (Abbreviated Final Report: 20-SI-006)

Predicting, responding to, or interpreting the chemical record preserved in debris derived from nuclear events can be challenging due to chemical fractionation. Chemical fractionation is where different species of the evolving radionuclide inventory segregate and/or are lost from the system over the timescales of debris formation. Both historic data and recent research suggest that the interaction and character of the local environment may exert controls on chemical fractionation by influencing the cooling and evolution of the associated fireball as well as the composition of the vapor term and resultant speciation. Prior to this work, an integrated platform permitting dynamic and concurrent consideration of physical and chemical evolution of early time post-detonation event environments did not exist. Our work merged historic data and experimental approaches to support development of a computational framework able to simulate fundamental processes (e.g., entrainment of local environment, oxidation chemistry, and cooling time scales) that may perturb the radionuclide inventory captured in post-detonation debris. Work with historic debris confirmed that entrained environmental material affect debris composition, structure, and radionuclide incorporation. Complementary work utilizing a readily controllable and tunable benchtop setup (a plasma flow reactor) simulated the late cooling of a nuclear fireball (e.g., T < 6000 K) and bounded the sensitivity of actinide speciation and particle size distribution to variations in oxygen concentration and cooling rates. Concurrent laser ablation and laser heating experiments were used to investigate the chemistry and physics of processes occurring in vaporized and/or rapidly heated actinides and other elements in the presence of oxygen. A more computationally efficient microphysical model was developed for predicting and evolving size distributions of particles forming from mixed vapor terms and simulating particle formation processes under a variety of extreme conditions. Continued study of historic nuclear event film confirmed that shockwave data and physics codes agree to within the uncertainty of the data. Good agreement was achieved for thermal emission from an airburst, however the paucity of low-temperature molecular opacity data for mixtures of air, bomb debris, entrained dirt, and water vapor complicate agreement for more elaborate scenarios. A multiphysics code (ALE3D) was modified to bring the necessary physics and chemistry, including these new data and insights, onto a single platform. Code development included improved initialization of large physical systems, modernization of chemistry capabilities, and modifications to enable inclusion of particle transport.

07 ISOTOPE AND RADIATION SOURCES↗

Integrating Experiments and Simulations to Reveal Anisotropic Growth Mechanisms and Interfaces of a One-Dimensional Zeolite

Zeolites are nanoporous crystalline materials critical for diverse industrial applications, yet their growth mechanisms are poorly understood. Here, this study presents a novel integrated framework combining experimental synthesis, high-resolution imaging, coarse-grained molecular dynamics simulations, and computer vision to uncover the mechanisms of growth of SSZ-24, a 1D channel zeolite. We demonstrate how synthesis conditions, such as temperature and reactant concentration, govern crystal anisotropy and surface roughness with growth dynamics differing markedly by crystallographic orientation. Along the channels, growth involves minimal energy barriers and rapid nucleation, resulting in rough surfaces. In contrast, growth perpendicular to the channels requires cooperative molecular organization and is highly sensitive to thermodynamic and kinetic conditions, yielding smooth anisotropic surfaces under low driving forces. By simulating transmission electron microscopy (TEM) images, we bridge molecular-scale simulations with experimental observations, identifying distinct growth mechanisms along different crystal planes. This work offers molecular-level insights into zeolite crystallization, advancing the rational design of nanoporous materials. The integration of cross-disciplinary methodologies establishes a transformative framework for optimizing zeolite synthesis, with implications for broader classes of materials.

Bertolazzo, Andressa A. [Univ. of Utah, Salt Lake ↗

Efficiency and Output Power Sensitivity to Control Parameter Variations in an LCC-Series Wireless Power Transfer System

In wireless power transfer (WPT) systems, output power (Pout) can be controlled through input parameters such as, input voltage (Vin) , duty cycle (d) , and switching frequency (fsw) . The impact of these control parameters on the overall system efficiency has not been explored in previous studies. This paper focuses on the efficiency (η) and Pout sensitivity to these control parameter variations in an LCC-Series tuned WPT system. There are multiple combinations of Vin, d , and fsw that can achieve a specified target Pout . However, it is important to determine the optimal combination of these control variables that will yield the maximum η for a target Pout . A MATLAB/Simulink model is built to calculate and measure η and Pout . In order to achieve this result, several simulations are run where Vin,d , and fsw are swept across a wide range to cover all possible combinations. MATLAB's Sensitivity Analysis toolbox is used to evaluate how the parameters influence the model's output. Finally, an experimental prototype is implemented. The optimal values acquired using the sensitivity analysis is then plugged into the hardware prototype and compared with the simulation. The results obtained confirm the validity of the proposed method.

Harave, Sudarshan↗

Improving GCM Representation of Convective Cloud Microphysics by Using ARM Raman Lidar and Cloud Radar Observations

The overall objective of this Department of Energy (DOE) Atmospheric System Research (ASR) funded project is to improve the representation of convective cloud microphysics in global climate models (GCMs) and check it by comparing the model simulations with observations. We have derived the cloud ice water content by synthesizing ARM Raman lidar (RL) and cloud radar observations at the Atmospheric Radiation Measurement (ARM) sites. Noting that the simulated anvil clouds in terms of their macro- and micro-physical properties are sensitive to the parameterization of convective microphysical processes, observed ice water content in anvil clouds provides a useful constraint on these parameterizations. We have improved the convective microphysics parameterization scheme by (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. We have examined the impact of improved convective cloud microphysics parameterization on the simulated global climate from GCM simulations.

54 ENVIRONMENTAL SCIENCES↗

Data-Driven Discovery of Linear Molecular Probes with Optimal Selective Affinity for PFAS in Water

Approaches to tackle the wide and growing variety of highly persistent per- and polyfluoroalkyl substances (PFAS) are of pressing global need because of their detrimental human health effects, such as cancer, birth defects, and hormone imbalance. Sensitive, selective, and easy-to-use real-time sensors to monitor and detect PFAS and sorbents to extract them are critical to meeting government-mandated environmental concentrations. In this work, we combine all-atom molecular dynamics simulations, enhanced sampling, deep representational learning, and Bayesian optimization to perform high-throughput virtual screening for highly sensitive and selective molecular probes. Our molecular design space consists of 3850 linear hydrocarbon chains with varying degrees of halogenation with and without amine- and phosphine-based headgroups. By employing a data-driven search process, we efficiently explore the molecular design space to optimize the sensitivity to perfluorooctanesulfonic acid (PFOS) as a prototypical PFAS analyte and selectivity relative to a sodium dodecyl sulfate (SDS) interferent. We calculate 504 Gibbs free energies of probe-analyte and probe-interferent interactions and identify probes with PFOS association free energies of up to (-ΔG PFOS ) = 9.8 ± 0.2 kJ/mol and selectivities relative to SDS of (-ΔΔG PFOS–SDS ) = 3.1 ± 1.5 kJ/mol. A C 11 Br 23 P(CH 3 ) 2 probe containing 11 backbone brominated carbons and a tertiary phosphine headgroup possesses the most sensitive binding constant to PFOS within the defined search space of K b PFOS = 177.4 ± 12.7, and a semibrominated probe C 5 H 11 C 7 Br 14 N(CH 3 ) 2 containing 12 backbone carbons and a tertiary amine headgroup possesses the highest selectivity relative to SDS of K b PFOS /K b SDS = 4.6 ± 1.7. A retrospective analysis of our data to extract interpretable design rules reveals that the sensitivity of linear hydrogenated probes increases by approximately 1 kJ/mol per C–C bond. The addition or removal of halogen atoms and amine or phosphine headgroups produces nonmonotonic changes in both sensitivity and selectivity with changes to the sensitivity of up to 2.5 kJ/mol. Finally, this work places empirical limitations on the performance of a wide range of linear probes for PFOS detection and offers a generic strategy for high-throughput computational screening to promote selective and sensitive binding.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dataset for manuscript "Equipartition and the temperature of maximum density of TIP4P/2005 water"

We simulate TIP4P/2005 water in the temperature range of 257 K to 318 K with time-steps 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density-temperature behavior obtained using 0.25 or 0.50 fs are in excellent agreement with each other but differ from those obtained using time-steps that have been shown earlier to lead to a breakdown of equipartition. The temperature of maximum density (TMD) is 277.15 K with time-step 0.25 or 0.50 fs, but is shifted to progressively lower values for longer time-steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water-water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid-vapor phase envelope. We present a simple physically transparent reasoning to highlight the separation of the time-scales between translational and rotational motion. We also develop a metric, Chi, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating Chi is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time-steps (less than or equal to 0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct-ensemble sampling is fundamental to ensure reproducibility across codes and simulation alogrithms. The included dataset provides the raw data used in the preparation of the graphs noted in the manuscript.

36 MATERIALS SCIENCE↗

Validation of SAS4A/SASSYS-1 for predicting steady-state single-phase natural circulation

The validation of system analysis codes for nuclear reactor systems is required for the development and application of these computational tools. Designed as a comprehensive system analysis code for advanced nuclear reactors, SAS4A/SASSYS-1 requires validation of its physics model for capturing single-phase natural circulation behavior. To support the validation of SAS4A/SASSYS-1, high-precision experiments are performed capturing steady-state single-phase natural circulation on a scaled facility with comprehensive instrumentation. In this work, dedicated tests are performed quantifying the critical modeling parameters, and a single-phase natural circulation benchmark dataset is obtained with well-documented uncertainty and comprehensive facility description. The validation is then performed against the dataset examining the capability of SAS4A/SASSYS-1 in simulating steady-state single-phase natural circulation. The experimental facility is modeled in the candidate code. Solution verification is performed using Richardson-extrapolation-based estimators which quantify and restrict numerical errors from discretization. Input uncertainty provided by the benchmark dataset is forward propagated through the candidate code, quantifying the output uncertainty in a Monte Carlo approach. The composition of the output uncertainty is also quantified through a variance-based sensitivity analysis. With the uncertainty quantified for each individual condition, a detailed comparison between the simulation results and experimental data is performed covering the whole dataset. The results show consistent agreement for all primary parameters. The current validation activity provides a valuable benchmark dataset for the validation of system analysis codes in capturing single-phase natural circulation and demonstrates satisfactory prediction capability of SAS4A/SASSYS1 for steady-state single-phase natural circulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗