SEARCH · Search NASA

Results for “Model initialization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Enhancement of Feedback Efficiency by Active Galactic Nucleus Outflows via the Magnetic Tension Force in the Inhomogeneous Interstellar Medium

By performing three-dimensional magnetohydrodynamics simulations of subrelativistic jets and disk winds propagating into the magnetized inhomogeneous interstellar medium (ISM), we investigate the magnetic effects on the active galactic nucleus feedback. Our simulations reveal that the magnetic tension force promotes the acceleration of the dense gas clouds, since the magnetic field lines, which are initially straight, bend around the gas clouds. In the jet models, the velocity dispersion of the clouds increases with an increase in the initial magnetic fields. The increment of the kinetic energy of the clouds is proportional to the initial magnetic fields, implying that the magnetic tension force increases the energy conversion efficiency from the jet to the gas clouds. Through simulations of the mildly collimated disk wind and the funnel-shaped disk wind, we confirm that such an enhancement of the energy conversion efficiency via the magnetic fields appears even if the energy is injected via the disk winds. The enhancement of the acceleration of the dense part of the magnetized ISM via the magnetic tension force will occur wherever the magnetized inhomogeneous matter is blown away.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Mobility Data and Models Informing Smart Cities

Using emerging data platforms, new mobility technologies, and travel demand models (TDMs), researchers, industry, and communities seek to improve the quality of transportation while maximizing the energy efficiency, equity, and safety of transportation services. As transportation may soon reach over 30% of U.S. energy consumption and with urban areas representing an increasing proportion of the U.S. population (>80% since 2010), a critical need exists to engage in urban data science-informed approaches to enhancing mobility. The objective of this study is to explore and document how aspiring Smart Cities are using data and models to inform mobility and energy initiatives within Smart City programs and in so doing identify gaps in knowledge and processes guiding Smart City mobility investment strategies, programs, projects, and pilots. A primary focus of the Smart Cities studied was the creation of an integrated data sharing environment approach. Most of these systems are being developed in parallel with multiple new data analysis tools, while regional metropolitan planning organizations continue to slowly evolve TDMs to take into account impacts of long-term strategies for emerging mobility technologies and services. Smart City initiatives in the United States have keen interests in leveraging knowledge and research on the mobility benefits and risks of automated, connected, efficient/electric, and shared on-demand mobility services; and understanding the related energy, environmental, economic, and societal impacts of these shifts. The results serve to identify key gaps in data, knowledge, and methods required to advance energy efficient urban mobility innovation, and to enable research and analysis collaboration between Smart Cities and the U.S. Department of Energy's efforts enabling new Systems and Modeling for Accelerated Research in Transportation (SMART) Mobility.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Fast Reactor Physics Model Verification Studies using ARC and PyARC Workflows

PyARC was recently developed at Argonne National Laboratory to automate many of the tasks required in the ARC (Argonne Reactor Computation) fast reactor simulation workflow, from input file generation, code execution, data transfer between ARC codes, and output postprocessing. PyARC will likely be the path forward to train new users of the ARC codes with the goal of wide adoption by the national laboratories, academia, and industry. In particular, for the ANL-JAEA collaboration under the Civil Nuclear Working Group (CNWG) project agreement NE-01, PyARC will be used to model the Joyo and EBR-II reactors for comparisons with measured data and calculated results from JAEA (Task 3: Fast Reactor Fuel and Core). As an additional avenue for verification and validation, this report investigates the use of PyARC towards a variety of existing ARC-based reactor models, in order to understand its efficacy in replicating the behavior of base ARC codes and better understand any limitations within modeling realistic fast reactor problems. To this end, PyARC was used to model the Joyo MKI, RBEC Benchmark-M, PRISM Mod-B, and EBR-II Run 138B cores, and its results were compared to those from existing ARC-based models. It was found that for hexagonal-based geometries PyARC was able to replicate the behavior of ARC codes to within 10 pcm for small reactor cores, and ~150pcm difference in eigenvalue for larger cores. These discrepancies are attributed primarily to differences in local mesh refinement options between ARC and PyARC, which currently cannot be resolved with PyARC’s latest version (1.6.0). In some of these cases, PyARC was used to model steady-state problems with initial core compositions originating from a prior REBUS depletion calculation. While PyARC was not designed to support such steady-state calculations, workarounds were applied to replicate the behavior of ARC-based calculations as closely as possible. Thus, these results demonstrate the wide extent to which they can be applied to fast reactor problems while still providing immense benefit to the user in terms of automating and standardizing common routines within the fast reactor analysis workflow. This study concluded that PyARC will be suitable for modeling the steady-state conditions of the EBR-II and Joyo fast reactors as part of the CNWG project agreement.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

DOE OSTI2023

Uncertainty Aware Deep Learning for Fault Prediction Using Multivariate Time Series Signals

The superconducting radio-frequency cavities are a crucial component of the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab. When a cavity faults, beam delivery to experimental end users is disrupted. Prediction of cavity faults prior to onset is essential to reduce operation and maintenance costs. In this work, a parallel long short-term memory (LSTM)-convolution neural network (CNN)-based deep learning (DL) model is proposed to predict impending faults using pre-fault signals. Further, we introduce an uncertainty quantification approach using Monte Carlo dropout with the LSTM-CNN model to ascertain confidence in the prediction. The model was tested using multivariate time series signals from stable cavity operations and before faults. Initial results show that on the test dataset, the model can identify impending faults before their onset with an average 10-fold cross validation accuracy of 97.39% and a standard deviation of 0.12% using a 100-ms time window. It is also observed that the model performs better as the prediction time moves closer to the fault onset. For additional context, we compare the performance of the model with three machine-learning-based (ML) fault prediction models. Our proposed parallel LSTM-CNN-based DL method shows better performance than the ML-based methods.

Rahman, Md Monibor↗

DOE OSTI2021

Interpreting the Diurnal Cycle of Clouds and Precipitation in the ARM GoAmazon Observations: Shallow to Deep Convection Transition

The Green Ocean Amazon (GoAmazon) 2014/5 field campaign data are used to study the diurnal cycle of clouds and precipitation. Through a careful classification of days with shallow cumulus, congestus and deep convection, we investigate the major differences among locally generated convection regimes and the most important environmental factors governing the shallow-to-deep convection transition. On shallow cumulus days, a greater sensible heat flux drives deeper boundary layer growth, which entrains drier free-tropospheric air and lowers the relative humidity, thus leading to a significantly higher cloud base than those on days with deeper convection. Congestus and deep convection regimes exhibit distinct cloud top height distributions with noticeable differences in the vertical wind shear in the mid-troposphere, suggesting an important role of wind shear in limiting the vertical extent of convection. On deep convection days, with preexisting nocturnal convection or cold-pools from external disturbances, the timing of peak surface precipitation (12:00–13:00 LST) tends to be in-phase with the diurnal variation in surface fluxes. However, it takes longer for local deep convection to develop without these disturbances. A plume model with thermodynamic and dynamical constraints is developed to explore the relative importance of various convection-controlling factors. Initial cloud-base vertical velocity and buoyancy are important in helping parcels ascend to the level of free convection (LFC). After parcels reach the LFC, entrainment of environmental air and lower free troposphere humidity become crucial in determining cloud top. Entrainment rate differentiates among convection regimes, which may be tied to the cloud size distribution at cloud base.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Development of Explainable, Knowledge-Guided AI Models to Enhance the E3SM Land Model Development and Uncertainty Quantification

Focal Area(s): (2)Predictive modeling using AI techniques and AI-derived model components; use of AI and other tools to design a prediction system comprising of a hierarchy of models. (3) Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge- guided AI. Science Challenge: The Energy Exascale Earth System Model (E3SM) is a fully coupled, state-of-the-science Earth system model that uses code optimized for DOE's advanced computers to address the most critical scientific questions facing our nation and society (Golaz et al., 2019). The E3SM Land model (ELM) is designed to understand how the changes in terrestrial land surfaces will interact with other Earth system components and has been used to understand hydrologic cycles, biogeophysics, and ecosystem dynamics. In spite of great successes, the ELM has several known issues that restrain rapid improvements. For example, the ELM uses equilibrium models to simulate dynamic land-climate interactions and it requires long model spin-up time to identify suitable initial conditions for transient simulations. The ELM lacks built-in uncertainty mechanisms that can improve the robustness of model predictions. The ELM is a holistic, deterministic model system with a rigid design, and in many situations, it is hard to modify the ELM system to incorporate new theory/hypothesis and new data across scales to address emerging science problems (such as predicting the impacts of water cycle extremes). In addition, The ELM is technically optimized for traditional CPU-centric computers and it cannot fully utilize the current and incoming leadership computers for model simulations and uncertainty quantification (UQ). The success of artificial intelligence (AI) has inspired scientists to use AI models to discover intrinsic features from simulation data (Chattopadhyay et al., 2020) and observational data (Reichstein et al., 2019) to gain further process understanding of Earth science problems. However, autonomous AI model training through deep learning usually requires a huge amount of annotated data. To overcome the limitations from the data and computing resources, knowledge-guided AI models are necessary where human-knowledge is ingested in model construction (Banino et al., 2018) and training process (Silver et al., 2016) for efficient learning. Herein, we present a new way that leverages the process understanding from the ELM to guide AI model development for the ELM enhancement and UQ. We hope this study can inspire further Earth and environmental system model developments and transformations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Simulated hydrological dynamics and coupled iron redox cycling impact methane production in an Arctic soil" in the Journal of Geophysical Research-Biogeosciences. We simulated biogeochemical cycling in arctic soils using the PFLOTRAN geochemical model combined with measurements from previous NGEE Arctic incubations of polygonal permafrost soils in northern Alaska (Zheng et al., 2018). Simulated iron cycling, carbon dioxide production, and methane production were compared with incubation measurements and the parameterized model was then used to simulate coupled iron and carbon cycling over repeated oxic-anoxic cycles at different levels of carbon substrate availability and pH. The most recent data version (2.0) in the archive incorporates changes to the model and simulations as suggested by reviewers during the manuscript review process. These changes include an updated parameterization of the model; a new set of simulations omitting the iron cycle for direct evaluation of how iron cycle processes affect modeled outcomes; and a set of simulations testing different scenarios of carbon substrate availability in addition to scenarios of initial soil pH. This archive contains simulation code, model output, and analysis code for PFLOTRAN simulations. All scripts are python except the batch script for submitting multiprocessor jobs. Note that the model also requires compiled versions of the Alquimia interface and the NGEE Arctic fork of the PFLOTRAN geochemical simulator (see the README_INSTALL document for basic instructions). The Output directory contains eight data files in netCDF format generated by the model. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 10-year research effort (2012-2022) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Fast Decomposition of Three-Component Spectra of Fluorescence Quenching by White and Grey Methods of Data Modeling

‘White’ and ‘grey’ methods of data modeling have been employed to resolve the heterogeneous fluorescence from a fluorophore mixture of 9-cyanoanthracene (CNA), 10-chloro-9-cyanoanthracene (ClCNA) and 9,10-dicyanoanthracene (DCNA) into component individual fluorescence spectra. The three-component spectra of fluorescence quenching in methanol were recorded for increasing amounts of lithium bromide used as a quencher. The associated intensity decay profiles of differentially quenched fluorescence of single components were modeled on the basis of a linear Stern-Volmer plot. These profiles are necessary to initiate the fitting procedure in both ‘white’ and ‘grey’ modeling of the original data matrices. ‘White’ methods of data modeling, called also ‘hard’ methods, are based on chemical/physical laws expressed in terms of some well-known or generally accepted mathematical equations. The parameters of these models are not known and they are estimated by least squares curve fitting. ‘Grey’ approaches to data modeling, also known as hard-soft modeling techniques, make use of both hard-model and soft-model parts. In practice, the difference between ‘white’ and ‘grey’ methods lies in the way in which the ‘crude’ fluorescence intensity decays of the mixture components are estimated. In the former case they are given in a functional form while in the latter as digitized curves which, in general, can only be obtained by using dedicated techniques of factor analysis. In the paper, the initial values of the Stern-Volmer constants of pure components were evaluated by both ‘point-by-point’ and ‘matrix’ versions of the method making use of the concept of wavelength dependent intensity fractions as well as by the rank annihilation factor analysis applied to the data matrices of the difference fluorescence spectra constructed in two ways: from the spectra recorded for a few excitation lines at the same concentration of a fluorescence quencher or classically from a series of the spectra measured for one selected excitation line but for increasing concentration of the quencher. The results of multiple curve resolution obtained by all types of the applied methods have been scrutinized and compared. In addition, the effect of inadequacy of sample preparation and increasing instrumental noise on the shape of the resolved spectral profiles has been studied on several datasets mimicking the measured data matrices. .

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2018

Four Perspectives on Human Bias in Visual Analytics

Mixed-initiative visual analytic systems steer analytical models and adapt views by making inferences from users' behavioral patterns with the system. Because such systems rely on incorporating implicit and explicit user feedback, they are particularly susceptible to the injection and propagation of human biases. To ultimately guard against the potentially negative effects of systems biased by human users, we must first qualify what we mean by the term bias. Thus, in this paper we describe four different perspectives on human bias that are particularly relevant to visual analytics. We discuss the interplay of human and computer system biases, particularly their roles in mixed-initiative systems. Given that the term bias is used to describe several different concepts, our goal is to encourage mindful consideration of which perspective on bias researchers take in their work to facilitate a common language in research and development efforts going forward.

visual analytics, mixed-initiative systems, cognit↗

DOE OSTI2017

Modeling and Validation of Irradiation Damage in Ni-based Alloys for Long-Term LWR Applications

Grizzly is a multiphysics simulation code under development at Idaho National Laboratory to simulate aging mechanisms and their effects on the integrity of critical LWR components, including Ni-Cr alloy components. Predictions of the LWR component integrity rely on the ductile-to-brittle transition temperature (DBTT) curve to determine the toughness of the alloy after a period of exposure to in-service reactor conditions. The macroscopic DBTT curve is governed by the microstructure evolution of the alloy, thus, a capability to predict the microstructure evolution in response to thermal ageing is being implemented into Grizzly. Within this report section, the implementation of the Avrami equation to estimate the phase fraction of the aged long range ordered precipitates is discussed. The initial evaluation of current crystal plasticity models within Grizzly to model these Ni-Cr alloys is also presented.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Skeletal Kinetics Reduction for Astrophysical Reaction Networks

A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to Type Ia supernovae (SNe Ia). The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of 56 Ni in SNe Ia is produced for different temperatures, densities, and proton-to-neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of α-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and 44 Ti and 56 Ni abundances by each model. The consequences of ye ≠ 0.5 in the initial composition are also explored where ye is the electron fraction. The simulated results show that 56 Ni production decreases by decreasing ye as expected, and that the 43 Sc is a key isotope in proton and neutron channels toward 56 Ni production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the 56Ni abundance in SNe Ia for ye ≲ 0.5 initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Identifying arbitrary parameter zonation using multiple level set functions

In this paper, we extended the analytical level set method [1, 2] for identifying a piece-wisely heterogeneous (zonation) binary system to the case with an arbitrary number of materials with unknown material properties. In the developed level set approach, starting from an initial guess, the material interfaces are propagated through iterations such that the residuals between the simulated and observed state variables (hydraulic head) is minimized. We derived an expression for the propagation velocity of the interface between any two materials, which is related to the permeability contrast between the materials on two sides of the interface, the sensitivity of the head to permeability, and the head residual. We also formulated an expression for updating the permeability of all materials, which is consistent with the steepest descent of the objective function. The developed approach has been demonstrated through many examples, ranging from totally synthetic cases to a case where the flow conditions are representative of a groundwater contaminant site at the Los Alamos National Laboratory. These examples indicate that the level set method can successfully identify zonation structures, even if the number of materials in the model domain is not exactly known in advance. Although the evolution of the material zonation depends on the initial guess field, inverse modeling runs starting with different initial guesses fields may converge to the similar final zonation structure. These examples also suggest that identifying interfaces of spatially distributed heterogeneities is more important than estimating their permeability values.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Characterizing initial- and final-state effects of relativistic nuclear collisions

Here, the multiphase transport model (AMPT) is used to study the final-state effects on the symmetric correlations (SC), asymmetric correlations (ASC), normalized symmetric correlations (NSC), and normalized asymmetric correlations (NASC) in Au+Au collisions at 200 GeV. The correlators' sensitivity to nonflow effects associated with long- and short-range nonflow correlations are also discussed using the HIJING model. The results indicate that SC, ASC, NSC, and NASC can give accompanying constraints for initial- and final-state effects. In addition, conducting further detailed experimental measurements spanning a broad range of collision systems and beam energies will serve as an additional constraint for the theoretical models' calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

A near-stationary subspace for ridge approximation

Response surfaces are common surrogates for expensive computer simulations in engineering analysis. However, the cost of fitting an accurate response surface increases exponentially as the number of model inputs increases, which leaves response surface construction intractable for high-dimensional, nonlinear models. Here we describe ridge approximation for fitting response surfaces in several variables. A ridge function is constant along several directions in its domain, so fitting occurs on the coordinates of a low-dimensional subspace of the input space. We review essential theory for ridge approximation – e.g., the best mean-squared approximation and an optimal low-dimensional subspace – and we prove that the gradient-based active subspace is near-stationary for the least-squares problem that defines an optimal subspace. Motivated by the theory, we propose a computational heuristic that uses an estimated active subspace as an initial guess for a ridge approximation fitting problem. We show a simple example where the heuristic fails, which reveals a type of function for which the proposed approach is inappropriate. We then propose a simple alternating heuristic for fitting a ridge function, and we demonstrate the effectiveness of the active subspace initial guess applied to an airfoil model of drag as a function of its 18 shape parameters.

42 ENGINEERING↗

DOE OSTI2024

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗

DOE OSTI2017

Initial 1.5D BISON Simulation Verification and Validation CASL FY17 Letter Report

This milestone report describes the current efforts to improve the usability of the recently developed 1.5 dimensional capability in BISON, the initial verification of this BISON capability, and initial validation of the 1.5 dimensional models against existing 2 dimensional axisymmetric full length light water reactor (LWR) fuel rods. The 1.5 dimensional (1.5D) capability in BISON was developed to enable the simulation of LWR fuel rods in a less computationally demanding simulation formulation while retaining key physics of interest. BISON was originally developed to solve 2D axisymmetric and 3D models, and this original development focus allows BISON to examine a wider range of fuel types and geometrically complexity than can be solved with lower dimensional frameworks. As the set of problems which BISON is used to solve expands, the software design requirements continue to evolve. Coupling physics in a higher dimensional geometry is not always required to resolve the physics of interest, such as the modeling of radial heat transfer and radial displacement. BISON, built on the dimension-agnostic MOOSE framework, is able to solve problems in 3D, 2D, and 1D with only minor modifications to the material definition code so that the appropriate geometrical representation can be applied in different problems of interest. The motivation for the development and validation of the BISON 1.5D capability is the CASL program's need for a simplified and computationally efficient version of BISON for use in the large multi-physics VERA simulations which include hundreds of individual fuel rods. The original development of the BISON 1.5D simulation capability is described in the CASL FY17 Letter Report from January 2017, and this report documents the continuing work identified in the previous report. The improvements to the usability of the 1.5D capability through input file simplifications, the results of the verification simulations, the results of the initial validation against the 2D assessment cases in BISON, and the future development needs are included within this report.

1.5D Capability↗

DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗