SEARCH · Engineering Papers

Results for “Multiscale Modeling and 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.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Improvement of the $\mathrm{BISON U_3Si_2}$ modeling capabilities based on multiscale developments to modeling fission gas behavior

Uranium silicide (U 3 Si 2 ) is a concept explored as a potential alternative to UO 2 fuel used in light water reactors (LWRs) since it may improve accident tolerance and economics due to its higher thermal conductivity and increased uranium density. U 3 Si 2 has been previously used in research reactors in the form of dispersion fuel, but operated at lower temperatures than commercial LWRs. The research reactor data illustrated that significant gaseous swelling occurs as the fuel burnup increases. Therefore, it is imperative to understand the fission gas behavior of U 3 Si 2 under higher temperature LWR operating conditions. In this work, molecular dynamics and phase-field modeling techniques are used to reduce the uncertainty in select modeling assumptions made in developing the fission gas behavior model for U 3 Si 2 in the BISON fuel performance code. These lower length scale informed models are then utilized in the validation of BISON U 3 Si 2 modeling capabilities to simulate the ATF-1 experiments irradiated in the Advanced Test Reactor (ATR). Sensitivity analysis (SA) and uncertainty quantification (UQ) are included as part of the validation process to identify where further experiments and lower length scale modeling would be beneficial. Here, the multiscale modeling approach utilized in this work can be applied to new fuel concepts being explored for both LWRs and advanced reactors (e.g., uranium nitride, uranium carbide).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRSDate not supplied

Thermoviscoelastic modelling of high strain thin-ply composites by means of multiscale plate and beam model

High strain thin-ply (HS-TPC) technology is being increasingly adopted for high-performanceaerospace applications. Albeit many of these structures such as deployable booms can bemodeled as one-dimensional beam problems, there is a lack of themoviscoelastic beam mod-els to efficiently and accurately simulate HS-TPC. This work will use mechanics of struc-ture genome (MSG) to construct linear thermoviscoelastic beam models that can homogenizethree-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integralformulation for thermorheologically simple materials subject to finite temperature changeswith the restriction that the strain is small. A lenticular boom is used as a numerical exampleto verify the MSG-based linear thermoviscoelastic beam model against MSG-based shell/platemodel, which has already been validated against experimental data provided by NASA, anddirect numerical simulations performed in a finite element commercial package.

Orzuri Rique Garaizar↗

DOE OSTI2024

Coupled SAM/Griffin Model of a Reference Fluoride-Salt-Cooled High-Temperature Reactor for Multi-Physics Simulations

A multi-physics coupled simulation model of a reference pebble bed fluoride-salt-cooled high-temperature reactor (PB-FHR) has been developed with SAM and Griffin computer codes for transient safety analysis. The reference problem of a prototypical reactor design serves as the foundation for the U.S. NRC (Nuclear Regulatory Commission) to verify the adequacy of computer codes and evaluation models for specific reactor types. In this work, the previously developed SAM model for PB-FHR primary loop has been updated for the coupled simulation. The updated SAM PB-FHR model includes a 2-D axial symmetric core region and external core components in 0-D/1-D. In addition to the primary loop, a detailed model of the RCCS (reactor cavity cooling system) is added. The 2-D and 1-D domains are tightly coupled using the single-solve approach developed in SAM. In the pebble bed region, the SAM multiscale explicit pebble model is applied to calculate the pebble and TRISO fuel kernel temperatures. The Griffin model used in this work is based on a model developed at Idaho National Laboratory in collaboration with the U.S. NRC. The Griffin neutronics model and SAM thermal hydraulics model is coupled with the Comprehensive Reactor Analysis Bundle (CRAB or alternately BlueCRAB) application. Both steady-state and transient scenarios are simulated to demonstrate the model's suitability for multi-physics simulations of PB-FHR transients.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

The impact of capillary heterogeneity on CO 2 flow and trapping across scales

Capillary heterogeneity has been identified over the last decade as a key control on subsurface CO 2 flow behavior during geological CO 2 sequestration. These heterogeneities can be formed in all sedimentary rocks, ranging from slight variations in the sand grain sizes to extensive sequences of interbedded sands, shales, and limestones. Capillary heterogeneity has been largely, although not entirely, overlooked in subsurface flow modeling because it is assumed to only directly influence fluid redistribution over scales of centimeters to meters. However, even small-scale fluid movements can result in dramatic impacts on the mobility and trapping of the CO 2 over kilometers. Therefore, neglecting capillary heterogeneity at multiple scales could potentially lead to errors in modeling and predicting field-scale plume migration. In this review paper, we aim to provide a consistent overview to (1) establish that capillary heterogeneity can have a major impact on CO 2 plume migration, (2) establish the respective length scales at which capillary heterogeneity matters, and (3) provide guidance for numerical modeling. This review covers pertinent literature and extracts key observations from the core to the field scales. Experimental studies have shown that millimeter-decimeter scale capillary heterogeneity can cause the so-called capillary heterogeneity trapping in addition to pore-scale residual trapping. Even at such a small scale, capillary heterogeneity can already lead to complex upscaled constitutive relationships, such as flow-rate dependent and anisotropic relative permeability, which affects field-scale CO 2 migration even when field-scale heterogeneities are present. Under gravity-dominated flow regimes, centimeter-meter scale capillary heterogeneity can entrap a significant amount of CO 2 at field scale, not just after imbibition but also during drainage. In certain cases, the presence of capillary heterogeneity can even completely stop the vertical movement of the CO 2 plume, hence greatly reducing leakage risks. At meter-kilometer scale, the influence of capillary heterogeneity is more pronounced and can hinder or redirect CO 2 migration in both lateral and vertical directions. The impact of capillary heterogeneity across multiple spatial scales poses a great challenge in modeling CO 2 migration at field scale, because it is practically impossible to build a field-scale earth model with grid blocks at millimeter scale. We recommend a hierarchical modeling approach to address this challenge. At field scale, earth models are built to capture geological features and heterogeneities in high but still practical grid resolutions. For each facies or rock type of the field-scale model, high- resolution meter-scale “conceptual” models are built with millimeter-scale grid blocks to capture representative fine-scale bedding geometries and heterogeneities in various environments of deposition, bridging the gap from subcore scale to the size of a field-scale simulation grid block. Upscaling is then used to preserve the smaller-scale flow dynamics of various rock types in field-scale simulations. Here, future work is needed to (1) refine, improve, and validate the hierarchical modeling approach; (2) build libraries of fine-scale bedding models for facies in various environments of deposition; (3) quantify multiscale capillary heterogeneity effects under subsurface uncertainties; (4) gain learning from different storage formations; and (5) establish best practices that balance accuracy and computational speed.

Capillary heterogeneity↗

DOE OSTI2022

BattPhase—A Convergent, Non-Oscillatory, Efficient Algorithm and Code for Predicting Shape Changes in Lithium Metal Batteries Using Phase-Field Models: Part I. Secondary Current Distribution

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. This article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

25 ENERGY STORAGE↗

NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

We present the results of a combined experimental and phase-field modeling study of pattern formation during freeze casting. Those studies use unidirectional freezing of simple binary liquid mixtures of water and sugars (sucrose and trehalose) in a temperature gradient, which suffice to produce hierarchical templated structures similar to those observed in more complex multi-component freeze-cast systems including lamellae, undulated ridges, and more exotic “jellyfish-like” substructures. Multiscale 3D phase-field simulations reproduce remarkably well those structures quantitatively and identify key properties of the ice-water interface that control their formation. They further reveal that lamellae form as a result of a novel symmetry-breaking secondary instability of partially faceted cellular structures and pinpoint additional secondary instability mechanisms giving rise to smaller-scale substructures.

Kaihua Ji↗

NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

Kaihua Ji↗

DOE OSTI2020

Multiscale Approach to Fast ModSim for Laser Processing of Metals for Future Nuclear Deterrence Environments

Predicting performance of parts produced using laser-metal processing remains an out- standing challenge. While many computational models exist, they are generally too computationally expensive to simulate the build of an engineering-scale part. This work develops a reduced order thermal model of a laser-metal system using analytical Green's function solutions to the linear heat equation, representing a step towards achieving a full part performance prediction in an "overnight" time frame. The developed model is able to calculate a thermal history for an example problem 72 times faster than a traditional FEM method. The model parameters are calibrated using a non-linear solution and microstructures and residual stresses calculated and compared to a non-linear case. The calibrated model shows promising agreement with a non-linear solution.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Physics-Informed Learning Machines for Multiscale and Multiphysics Problems (PHILMS) (Technical Report)

The research work at University of California Santa Barbara (UCSB) resulted in several new developments in the areas of scientific machine learning, numerical analysis, and practical methods for data-driven modeling, prediction, reductions, and simulation. Many of the projects were carried out in collaboration with members of the national laboratories at Sandia National Laboratories (SNL), Pacific Northwestern National Laboratories (PNNL), and other institutions. Results included developing new scientific machine learning methods, related theory and mathematical frameworks for analysis and training, data-driven numerical solvers, and related tools and software for scientific computation. During the support period, over 16+ papers were submitted for publication, and 4 open-source software packages were developed and released (available at http://atzberger.org/). In addition, 7+ students and 2 post-docs were mentored in collaboration with the laboratory staff for future careers in academia, government labs, and industry.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Multiscale Modeling of Radiation Damage in UO 2 under Accelerated Burnup Conditions

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Development of a MOOSE thermal model of the MPC-32 canister and HI-STORM overpack

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. Safe management of spent nuclear fuel (SNF) is a key aspect of the back end of the nuclear fuel cycle, and SNF dry storage systems are becoming a popular, effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution for maintaining SNF (~60 years) prior to final disposal. This project aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack, using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. This study will investigate and demonstrate the thermal-hydraulics capabilities of the MOOSE framework, including natural circulation, heat transfer, porous flows, etc. The ultimate goal of the project is to verify whether MOOSE tools (including Pronghorn) can be used to study the thermal performance of the SNF dry cask storage system. This study provides reliable and inclusive solving strategy for dry cask problems. The detailed information about the solving scheme and the governing equations related to the physics of the system is provided in the report. The results for thermal-hydraulic analysis of the HI-STORM system is produced with using open source modules of the MOOSE framework. This results highlights the flexibility and modularity of the MOOSE which makes it a unique candidate for the frameworks and code packages. Therefore, integration of the MOOSE to UNF ST&DARDS will improve the thermal-hydraulic capability of the system while providing distinctive features to users.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Implementation of a Quasi-3D Multiscale Modeling Framework to ACME

Representation of unresolved cloud processes has been one of the most challenging problems in numerical modeling of the atmosphere. During the last decades, considerable efforts have been spent to better represent the collective effects of those processes in terms of resolved scales (see Arakawa 2004 for a review). Still, the progress on this problem has been unacceptably slow as pointed out by Randall et al. (2003). The large part of the difficulty is inherent in the artificial separation of atmospheric processes into the explicitly simulated resolved scales and the highly parameterized subgrid scales. Arakawa (2004) and Arakawa et al. (2011) suggested two plausible approaches to eliminate the problem of artificial separation ofatmospheric processes: “parameterize everything” and “resolve everything”.

58 GEOSCIENCES↗

NASA NTRS2010

Multi-Scale Modeling of Global of Magnetospheric Dynamics

To understand the role of magnetic reconnection in global evolution of magnetosphere and to place spacecraft observations into global context it is essential to perform global simulations with physically motivated model of dissipation that is capable to reproduce reconnection rates predicted by kinetic models. In our efforts to bridge the gap between small scale kinetic modeling and global simulations we introduced an approach that allows to quantify the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions near reconnection sites. We utilized the high resolution global MHD code BATSRUS and incorporate primary mechanism controlling the dissipation in the vicinity of reconnection sites in terms of kinetic corrections to induction and energy equations. One of the key elements of the multiscale modeling of magnetic reconnection is identification of reconnection sites and boundaries of surrounding diffusion regions where non-MHD corrections are required. Reconnection site search in the equatorial plane implemented in our previous studies is extended to cusp and magnetopause reconnection, as well as for magnetotail reconnection in realistic asymmetric configurations. The role of feedback between the non-ideal effects in diffusion regions and global magnetosphere structure and dynamics will be discussed.

Kuznetsova, M. M.↗

DOE OSTI2022

Load‐Balancing Intense Physics Calculations to Embed Regionalized High‐Resolution Cloud Resolving Models in the E3SM and CESM Climate Models

Abstract We design a new strategy to load‐balance high‐intensity sub‐grid atmospheric physics calculations restricted to a small fraction of a global climate simulation's domain. We show why the current parallel load balancing infrastructure of Community Earth System Model (CESM) and Energy Exascale Earth Model (E3SM) cannot efficiently handle this scenario at large core counts. As an example, we study an unusual configuration of the E3SM Multiscale Modeling Framework (MMF) that embeds a binary mixture of two separate cloud‐resolving model grid structures that is attractive for low cloud feedback studies. Less than a third of the planet uses high‐resolution (MMF‐HR; sub‐km horizontal grid spacing) relative to standard low‐resolution (MMF‐LR) cloud superparameterization elsewhere. To enable MMF runs with Multi‐Domain cloud resolving models (CRMs), our load balancing theory predicts the most efficient computational scale as a function of the high‐intensity work's relative overhead and its fractional coverage. The scheme successfully maximizes model throughput and minimizes model cost relative to precursor infrastructure, effectively by devoting the vast majority of the processor pool to operate on the few high‐intensity (and rate‐limiting) high‐resolution (HR) grid columns. Two examples prove the concept, showing that minor artifacts can be introduced near the HR/low‐resolution CRM grid transition boundary on idealized aquaplanets, but are minimal in operationally relevant real‐geography settings. As intended, within the high (low) resolution area, our Multi‐Domain CRM simulations exhibit cloud fraction and shortwave reflection convergent to standard baseline tests that use globally homogenous MMF‐LR and MMF‐HR. We suggest this approach can open up a range of creative multi‐resolution climate experiments without requiring unduly large allocations of computational resources.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020: Part 3

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: Task2c: Explicit modeling of pebble transient temperature response. In this simulation, the 400 MWth Pebble-Bed Modular Reactor (PBMR) design, PBMR- 400, experiences a 20-second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling.; Extended scope part 1: An assessment of the computational efficiency of the Discontinuous Finite Element Method (DFEM) heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the Finite Element Method (FEM) solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems.; Extended scope part 2: The GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables.; Extended scope part 3: Demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat-up transient and used in the GapHeatTransferInterface model. A full-domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench-marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗

NASA NTRS2001

Building Blocks for Reliable Complex Nonlinear Numerical Simulations

This chapter describes some of the building blocks to ensure a higher level of confidence in the predictability and reliability (PAR) of numerical simulation of multiscale complex nonlinear problems. The focus is on relating PAR of numerical simulations with complex nonlinear phenomena of numerics. To isolate sources of numerical uncertainties, the possible discrepancy between the chosen partial differential equation (PDE) model and the real physics and/or experimental data is set aside. The discussion is restricted to how well numerical schemes can mimic the solution behavior of the underlying PDE model for finite time steps and grid spacings. The situation is complicated by the fact that the available theory for the understanding of nonlinear behavior of numerics is not at a stage to fully analyze the nonlinear Euler and Navier-Stokes equations. The discussion is based on the knowledge gained for nonlinear model problems with known analytical solutions to identify and explain the possible sources and remedies of numerical uncertainties in practical computations. Examples relevant to turbulent flow computations are included.

Yee, H. C.↗