Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model types”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale Control of Generic Second Order Traffic Models by Driver-Assist Vehicles

We study the derivation of generic high order macroscopic traffic models from a follow-the-leader particle description via a kinetic approach. First, we recover a third order traffic model as the hydrodynamic limit of an Enskog-type kinetic equation. Next, we introduce in the vehicle interactions a binary control modeling the automatic feedback provided by driver-assist vehicles and we upscale such a new particle description by means of another Enskog-based hydrodynamic limit. The resulting macroscopic model is now a generic second order model (GSOM), which contains in turn a control term inherited from the microscopic interactions. We show that such a control may be chosen so as to optimize global traffic trends, such as the vehicle flux or the road congestion, constrained by the GSOM dynamics. By means of numerical simulations, we investigate the effect of this control hierarchy in some specific case studies, which exemplify the multiscale path from the vehiclewise implementation of a driver-assist control to its optimal hydrodynamic design.

GSOM↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2021

Thermal conductivity of α -U with point defects

Here, we develop a theoretical model for thermal conductivity of α-U that combines density functional theory calculations and the coupled electron–phonon Boltzmann transport equation. The model incorporates both electron and phonon contributions to thermal conductivity and achieves good agreement with experimental data over a wide temperature range. The dominant scattering mechanism governing thermal transport in α-U at different temperatures is examined. By including phonon–defect and electron–defect scatterings in the model, we study the effect of point defects including U-vacancy, U-interstitial, and Zr-substitution on the thermal conductivity of α-U. The degradation of anisotropic thermal conductivity due to point defects as a function of defect concentration, defect type, and temperature is reported. This model provides insights into the impact of defects on both phonon and electron thermal transport. It will promote the fundamental understanding of thermal transport in α-U and provide a ground for investigation of coupled electron–phonon transport in metallic materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Monitoring Operational States of a Nuclear Reactor Using Seismoacoustic Signatures and Machine Learning

Monitoring nuclear reactors is an important safety and security task with growing requirements. We explore the possibility of using seismic and acoustic data for inferring the power level of an operating reactor. Continuous data recorded at a single seismoacoustic station that is located about 50 m away from a research reactor was visualized and analyzed. The data show a clear correlation between seismoacoustic features and reactor main operational states. We designed a workflow that includes two machine learning (ML) models to classify the reactor operational states (OFF, transition, and ON) and estimate reactor power levels (10%, 30%, 50%, 70%, and 90%). We applied and compared five ML algorithms for the reactor OFF-transition-ON and four approaches for the power level classification. We also compared the performance of ML models trained with seismic-only, acoustic-only, and both types of data. Five-fold cross validations were implemented to assure a thorough evaluation of the model performances. Additionally, the results show the extreme boosting gradient algorithm worked best for the first model, whereas random forests performed best for the second model. Combining seismic and acoustic data leads to better performance than using a single type of data. Seismic data contributed more than acoustic data for both models. We reached an accuracy of 0.98 for reactor OFF and ON. The accuracies for the transition state and power levels are less optimal with a minimum accuracy of 0.66. However, our results suggest seismic and acoustic data contain useful information about the transition state as well as power levels. Seismic and acoustic data could be integrated with other observations to improve monitoring performance.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

AI Applications to Physics Experiments at Jefferson Lab

We survey how AI/ML is being deployed across Jefferson Lab's experimental and accelerator programs. In EPSCI, Hydra applies computer vision to automate real-time data-quality monitoring across all four experimental halls, replacing manual inspection of hundreds to thousands of histograms per shift. AIEC (AI Experiment Controls) uses ML to stabilize drift chamber gains and is now part of standard CEBAF production running, while AI Optimized Polarization (AIOP) targets autonomous control of polarized targets and photon beam angular alignment. In CASA, cavity fault classification models identify faulted cavities and trip types from waveform data with ~85% and ~78% agreement to labeled data, respectively, and are deployed in production; a separate effort applies LLMs and hybrid search to make the CEBAF operations logbook AI-ready. QCD-focused work includes transformer- and GAN-based generative models for particle-level event simulation, with distributed GAN training scaling studies on Polaris. Additional efforts span ML-on-FPGA for the EIC and a new Data Science Department coordinating anomaly detection, uncertainty quantification, and HPC-scalable ML lab-wide. Collectively, these projects illustrate AI's growing role in improving efficiency across JLab's nuclear physics mission.

Mei, Xinxin [Thomas Jefferson National Accelerator↗
DOE OSTI2026

Effects of slope steepness and cross-slope wind speed on fire spreading behavior for various vegetation

Wildland fire behavior is significantly influenced by environmental factors such as slope steepness, wind speed, and fuel type. Understanding these interactions is critical for improving predictive models and fire management. This study explores how slope steepness and cross-slope wind speed influence fire spread dynamics in various fuel bed types. Simulations are conducted using a physics-based wildland fire model, HIGRAD/FIRETEC, across six slope angles (0–50 %), four cross-slope wind speeds (4–10 m s –1 ), and three fuel bed types (grass, shrubland, and forest). Representative cases are additionally compared with FARSITE fireline evolution. Fire behavior is categorized into distinct propagation types based on spread characteristics and analyzed. The fire propagation angle, which indicates deviation from the wind direction, generally increases with steeper slopes and decreases with stronger cross-slope winds. Secondary upslope propagation is observed in shrubland under moderate slopes, while secondary downwind propagation occurs in all fuel beds at higher wind speeds. These findings highlight fire spread characteristics that differ from predictions by traditional models like Rothermel’s. By capturing complex propagation patterns and dynamics, this study demonstrates the value of a physics-based, atmosphere-fire coupled model for accurate wildland fire prediction. Incorporating secondary propagations and the influence of fuel bed complexities into predictive models can improve the accuracy of fire spread forecasts, enabling more effective fire management and risk mitigation efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

First Principle based Thermal Conductivity Modeling of Metal Alloys with Applications on Uranium Alloys

Thermal conductivity is an important materials property related to heat transport, which is essential to many applications, ranging from thermoelectrics to nuclear reactor materials. High- quality thermal conductivity data is critical to these materials and their associated technologies, including the nuclear fuel materials like traditional oxide fuels and metallic uranium (U) fuels, e.g. U-Zr and U-Mo alloys. And thermal conductivity modeling is widely used to interpolate or extrapolate experimental data, to obtain high-quality thermal conductivity data over wide temperature and composition range, and to understand the impacts of different factors including defects and microstructures. However, a general thermal conductivity model for metal alloys which can work regrading different phase components is still missing. Also, for U alloys, a thermal conductivity model working with different phases is needed. Therefore, in this work, we developed thermal conductivity models for metal alloys, based on ab-initio calculations, semi-classical physics rules, and limited experimental data. The goals of these models are to help obtain high- quality thermal conductivity data within a little experimental input as possible, have models that are extendable to varying microstructures, including different types of irradiation effects, and provide mechanistic understanding of heat transfer in the modeled alloys. In this work, our model solves several challenges in the development. The DFT-BTE approach is applied to decrease the reliance on experimental data and make the model extendable to composition changes and defects. A practical and efficient way to combine the DFT inputs with physics rules is pointed out in this work too. A staged approach, which starts from simple cases of elemental metals, then extends to solid solutions, different compound phases, and multi-phase mixtures, is presented to work with the metal alloys in different phase components. Our models are demonstrated respectively on aU for elemental metal model, on U-Zr and U-Mo alloys in aU temperature range for multi-phase mixture model, on high-temperature U-Zr, U-Nb, and U-Mo alloys in the body-centered cubic phase for concentrated solid solution model, and on irradiated U-Mo alloys for irradiated metal alloy model. All models show great agreement with experimental data. In these demonstrations, our model shows its advantages compared to previous empirical fitting model. Our model requires fewer experimental data, due to the inputs from DFT. The quantitative insights into the different physical factors are provided in our model, as it incorporates the electron and phonon scattering mechanisms. Our model also can be extended to incorporate the contributions of point defects, grain boundaries, and noble gas bubbles, to integrate their effects on heat transfer. This model can serve as both a foundation for understanding the more complex thermal conductivity of realistic U alloy fuels and a useful tool to guide further modeling of thermal conductivity to aid materials and device design.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Capabilities of multivariate Bayesian inference toward seismic hazard assessment

Multivariate Bayesian analysis can bring significant benefits to seismic hazard analysis: Its multivariate feature enables computing scalar and vector hazard without making any approximations; Correlations between intensity measures are implicitly modeled, permitting direct simulation of ground motion selection tools such as the conditional mean spectrum and the generalized conditioning intensity measure; and Its updating feature enables a seamless integration of new ground motion data into the hazard results. Here, we first develop a multivariate Bayesian ground motion model through the NGA-West2 database. The model functional form considers fault-type, magnitude, and distance dependencies, and also the linear and the rock intensity dependent site response. We use a hybrid Markov Chain Monte Carlo sampling to perform Bayesian inference consisting of Gibbs step and a multilevel Metropolis-Hastings step. We then perform several checks on the model and note that its performance is satisfactory. Finally, we illustrate the merits of this multivariate Bayesian analysis, which include: ground motion model updating with ground motion data recorded in the last four years not part of the NGA-West2 database; computation of scalar and vector seismic hazard using the un-updated and updated ground motion models for Los Angeles, CA; and simulation of the conditional mean spectrum under scalar and vector IM conditioning while accounting for different sources of aleatoric and epistemic uncertainties.

58 GEOSCIENCES↗
DOE OSTI2021

Search for heavy resonances decaying into a pair of Z bosons in the $\ell ^+\ell ^-\ell '^+\ell '^-$ and $\ell ^+\ell ^-\nu {{\bar{\nu }}}$ final states using 139 $\mathrm {fb}^{-1}$ of proton–proton collisions at $\sqrt{s} = 13$ TeV with the ATLAS detector

A search for heavy resonances decaying into a pair of Z bosons leading to $\ell ^+\ell ^-\ell '^+\ell '^-$ and $\ell ^+\ell ^-\nu {{\bar{\nu }}}$ final states, where ℓ stands for either an electron or a muon, is presented. The search uses proton–proton collision data at a centre-of-mass energy of 13 TeV collected from 2015 to 2018 that corresponds to the integrated luminosity of 139 fb -1 recorded by the ATLAS detector during Run 2 of the Large Hadron Collider. Different mass ranges spanning 200 GeV to 2000 GeV for the hypothetical resonances are considered, depending on the final state and model. In the absence of a significant observed excess, the results are interpreted as upper limits on the production cross section of a spin-0 or spin-2 resonance. The upper limits for the spin-0 resonance are translated to exclusion contours in the context of Type-I and Type-II two-Higgs-doublet models, and the limits for the spin-2 resonance are used to constrain the Randall–Sundrum model with an extra dimension giving rise to spin-2 graviton excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Systematic approach to define a prototype building model: Convenience Store

The US Department of Energy supports the development of commercial building energy codes and standards. To support commercial building energy research activities, continuous efforts have been made to convert 30 prototype building models, which cover 80% of commercial building floor space in the United States, into OpenStudio Standards prototype buildings. In this paper, we provide a systematic approach to define the building and system characteristics of a convenience store, which has the largest energy use intensity among the commercial building types not covered by current prototype models. This study will be used as a basis for developing the prototype convenience store building model.

Yoon, Yeobeom↗
DOE OSTI2020

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

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

Inferring Plant Acclimation and Improving Model Generalizability With Differentiable Physics‐Informed Machine Learning of Photosynthesis

Net photosynthesis (A N ) is a key component of the global carbon cycle influencing climate feedback over decadal scales. Although plant acclimation to environmental changes can modify A N , traditional vegetation models in Earth system models (ESMs) often rely on plant functional type (PFT)-specific parameterizations or simplified acclimation assumptions limiting generalizability across time, space, and PFTs. In this study, we developed a differentiable photosynthesis model to learn the environmental dependencies of V c,max25 (maximum carboxylation rate at 25°C, representing photosynthetic capacity), as this genre of hybrid physics-informed machine learning can seamlessly train neural networks and process-based equations together. Compared to PFT-specific parameterization of V c,max25 , learning the environment dependencies of key photosynthetic parameters improved model spatiotemporal generalizability. Applying environmental acclimation to V c,max25 led to substantial variations in global mean A N indicating the need to address acclimation in ESMs. The model effectively captured multivariate observations (V c,max25 , A N , and stomatal conductance (g s )) simultaneously with multivariate constraints, improving generalization across space and PFTs. It also learned sensible acclimation relationships of V c,max25 to different environmental conditions. The model explained more than 54%, 57%, and 62% of the variance of A N , g s , and V c,max25 , respectively, presenting a first global-scale spatial test benchmark of A N and g s . These results highlight the potential for differentiable modeling to enhance process-based modules in ESMs and effectively leverage information from large, multivariate data sets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

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

Asymptotic consistency of the WSINDy algorithm in the limit of continuum data

In this work we study the asymptotic consistency of the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy) in the identification of differential equations from noisy samples of solutions. We prove that the WSINDy estimator is unconditionally asymptotically consistent for a wide class of models that includes the Navier–Stokes, Kuramoto–Sivashinsky and Sine–Gordon equations. We thus provide a mathematically rigorous explanation for the observed robustness to noise of weak-form equation learning. Conversely, we also show that, in general, the WSINDy estimator is only conditionally asymptotically consistent, yielding discovery of spurious terms with probability one if the noise level exceeds a critical threshold σ c . We provide explicit bounds on σ c in the case of Gaussian white noise and we explicitly characterize the spurious terms that arise in the case of trigonometric and/or polynomial libraries. Furthermore, we show that, if the data is suitably denoised (a simple moving average filter is sufficient), then asymptotic consistency is recovered for models with locally-Lipschitz, polynomial-growth nonlinearities. Our results reveal important aspects of weak-form equation learning, which may be used to improve future algorithms. We demonstrate our findings numerically using the Lorenz system, the cubic oscillator, a viscous Burgers-growth model and a Kuramoto–Sivashinsky-type high-order PDE.

asymptotic consistency↗
DOE OSTI2021

The Initial–Final Mass Relation for Hydrogen-deficient White Dwarfs

The initial–final mass relation represents the total mass lost by a star during the entirety of its evolution from the zero age main sequence to the white-dwarf cooling track. The semiempirical initial–final mass relation (IFMR) is largely based on observations of DA white dwarfs, the most common spectral type of white dwarf and the simplest atmosphere to model. We present a first derivation of the semiempirical IFMR for hydrogen-deficient (non-DA) white dwarfs in open star clusters. We identify a possible discrepancy between the DA and non-DA IFMRs, with non-DA white dwarfs ≈0.07 M {sub ⊙} less massive at a given initial mass. Such a discrepancy is unexpected based on theoretical models of non-DA formation and observations of field white dwarf mass distributions. If real, the discrepancy is likely due to enhanced mass loss during the final thermal pulse and renewed post-AGB evolution of the star. However, we are dubious that the mass discrepancy is physical and instead is due to the small sample size, to systematic issues in model atmospheres of non-DAs, and to the uncertain evolutionary history of Procyon B (spectral type DQZ). A significantly larger sample size is needed to test these assertions. In addition, we also present Monte Carlo models of the correlated errors for DA and non-DA white dwarfs in the initial–final mass plane. We find the uncertainties in initial–final mass determinations for individual white dwarfs can be significantly asymmetric, but the recovered functional form of the IFMR is grossly unaffected by the correlated errors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Evaluation of Mesoscale Convective Systems in Climate Simulations: Methodological Development and Results from MPAS-CAM over the U.S.

In this study, a process-oriented approach is developed to evaluate warm-season mesoscale convective system (MCS) precipitation and their favorable large-scale meteorological patterns (FLSMPs) over the U.S. This approach features a novel observation-driven MCS-tracking algorithm using infrared brightness temperature and precipitation feature at 12, 25 and 50 km resolution and metrics to evaluate the model large-scale environment favorable for MCS initiation. The tracking algorithm successfully reproduces the observed MCS statistics from a reference 4-km radar MCS database. To demonstrate the utility of the new methodologies in evaluating MCS in climate simulations with mesoscale resolution, the process-oriented approach is applied to two climate simulations produced by the Variable-Resolution Model for Prediction Across Scales coupled to the Community Atmosphere Model physics, with refined horizontal grid spacing at 50 km and 25 km over North America. With the tracking algorithm applied to simulations and observations at equivalent resolutions, the simulated number of MCS and associated precipitation amount, frequency and intensity are found to be consistently underestimated in the Central U.S., particularly from May to August. The simulated MCS precipitation shows little diurnal variation and lasts too long, while MCS precipitation area is too large and intensity is too weak. The model is able to simulate four types of observed FLSMP associated with frontal systems and low-level jets (LLJ) in spring, but the frequencies are underestimated because of low-level dry bias and weaker LLJ. Precipitation simulated under different FLSMPs peak during daytime, in contrast to the observed nocturnal peak. Implications of these findings for future model development and diagnostics are discussed.

54 ENVIRONMENTAL SCIENCES↗