Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale simulations”

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 181 records · Page 10

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.

Yee, H. C.↗

Physics-Reinforced Machine Learning Algorithms for Multiscale Closure Model Discovery

The central objective of this project was to address the challenge of modeling and simulating complex multiscale turbulence phenomena by leveraging physics-guided machine learning (PGML) and hybrid modeling approaches. By integrating physics-based methods with data-driven models, the research focused on achieving robust and scalable solutions for geophysical turbulence, enhancing numerical weather prediction and climate research tools. The project resulted in significant advancements in computational modeling paradigms, predictive tools for reduced-order modeling, and innovative algorithms for fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗

MrHyDE v.1.0

SAND2024-01324O MrHyDE, which stands for Multi-resolution Hybridized Differential Equations, is a general-purpose C++ package for the solution of coupled multiphysics and multiscale systems on massively parallel computing systems. MrHyDE is designed to enable moving beyond forward simulation for multiscale applications which includes optimization, control, uncertainty quantification, and stochastic inversion. The framework provides interfaces to several packages within the Trilinos framework and leverages automatic differentiation to enable adjoint capabilities for large-scale, gradient-based optimization. MrHyDE provides automated multiscale capabilities through a subgrid model interface and multiscale Dirichlet-to-Neumann maps. For extreme-scale applications, MrHyDE provides in situ data-compression algorithms to reduce memory requirements while maintaining performance. MrHyDE is a general-purpose, computational framework for the solution of multiscale and multiphysics applications. It uses a combination of structure-preserving, physics-compatible discretizations, fully implicit methods, multi-resolution schemes, or fully explicit methods. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

HI-STORM Overpack and MPC-32 Thermal-Hydraulic Model with MOOSE Framework

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. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and 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 to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work 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. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and 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 to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work 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. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

MuSIKAL: Multiphysics Simulations and Knowledge Discovery through AI/ML Technologies

Under the MuSiKAL project, we developed a framework for a coastal digital twin (DT) platform capable of integrating diverse data resources, configuring multiscale model simulations, performing SciML‐accelerated predictions, with applications primarily driven by storm surge and heavily rainfall events impacting the Gulf Coast of the U.S.

54 ENVIRONMENTAL SCIENCES↗

Multiscale Failure Analysis of a 3D Woven Unit Cell Containing Defects

Multiscale failure simulations have been performed for a Three-dimensional woven composite unit cell considering five, or more, length scales spanning the woven composite mesoscale to the sub-microscale voids. The multiscale recursive micromechanics approach, which enables recursive integration of general micromechanics theories over an arbitrary number of length scales, has been employed within the NASA Multiscale Analysis Tool. The multiscale model uses both the generalized method of cells and Mori-Tanaka micromechanics theories, and considers failure in the constituent materials using a simple damage model. Baseline results, containing distributed voids, are compared to uniaxial experimental data for an AS4 carbon fiber/ RTM6 epoxy matrix 3D orthogonal woven composite with good agreement in terms of global stiffness and global failure stress. The simulations demonstrate that the 3D woven composite exhibits damage tolerance through sustaining increasing axial load far beyond the first initiation of damage. The multiscale model is used to examine the nonlinear response of the material to other loading conditions. Case studies, motivated by X-ray computed tomography data, are presented on the effects of manufacturing induced voids and cracks.

3D woven↗

Multiscale Cloud System Modeling

The central theme of this paper is to describe how cloud system resolving models (CRMs) of grid spacing approximately 1 km have been applied to various important problems in atmospheric science across a wide range of spatial and temporal scales and how these applications relate to other modeling approaches. A long-standing problem concerns the representation of organized precipitating convective cloud systems in weather and climate models. Since CRMs resolve the mesoscale to large scales of motion (i.e., 10 km to global) they explicitly address the cloud system problem. By explicitly representing organized convection, CRMs bypass restrictive assumptions associated with convective parameterization such as the scale gap between cumulus and large-scale motion. Dynamical models provide insight into the physical mechanisms involved with scale interaction and convective organization. Multiscale CRMs simulate convective cloud systems in computational domains up to global and have been applied in place of contemporary convective parameterizations in global models. Multiscale CRMs pose a new challenge for model validation, which is met in an integrated approach involving CRMs, operational prediction systems, observational measurements, and dynamical models in a new international project: the Year of Tropical Convection, which has an emphasis on organized tropical convection and its global effects.

Tao, Wei-Kuo↗

Quantification of Numerical Uncertainty via Nonlinear Dynamical Approach

Motivations (Ensure a Higher Level of Confidence in the Predictability & Reliability of Numerical Simulation for Multiscale Complex Nonlinear Fluid Problems) - The last two decades have been an era when computation is ahead of analysis & when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems & non-traditional fields (Especially when computations offer the ONLY way of generating this type of data limited simulations). - At present some of the numerical uncertainties can be explained and minimized by traditional numerical analysis and standard CFD practices. However, such practices, usually based on linearized analysis, MIGHT NOT be sufficient for strongly nonlinear and/or stiff problems. - We need a good understanding of the nonlinear behavior of numerical schemes being used as an integral part of code verification, validation and certification.

HEC↗

Multiscale numerical investigation of ratchet growth damage effects in PBX 9502

This paper presents results of numerical experiments conducted on the high explosive PBX 9502 to investigate how recently observed grain-scale damage mechanisms of ratchet growth affect uniaxial compression measurements. Simulations are multiscale in the sense of directly resolving grains, pores, cracks, and grain-interfaces based upon scanning electron microscope (SEM) images of damaged and undamaged samples. The combined finite-discrete element method (FDEM) is utilized to resolve both grain-scale microfracture and elastoplastic deformation of solid grains. Pristine (undamaged) and damaged microstructures are compared in simulation of unconfined compression tests of the same material from the literature. Here, the simulation results show the observed microscale mechanisms of damage, specifically microfracture predominantly around and sometimes through grains and crack-associated pore growth, can well-explain the effective degradation of strength and stiffness observed in the laboratory measurements.

36 MATERIALS SCIENCE↗

Three-dimensional inhomogeneity of electron-temperature-gradient turbulence in the edge of tokamak plasmas

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Here our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport—in our case by ~40%—through multiscale interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: Presentation to DOE-NE [Slides]

The research goal of this project is to develop at Idaho National Laboratory (INL) a first-of-its-kind Multiphysics Object-Oriented Simulation Environment (MOOSE)-based, multiscale, multiphysics spark plasma sintering (SPS) modeling and simulation code application, termed “Freya.” Freya will simulate the thermo-mechanical-electrical aspects of the SPS fabrication process and will be paired with lower length scale sub-models, such as phase-field, to predict the resulting microstructure. SPS is an advanced manufacturing process that can be used to solve a variety of material manufacturing challenges; however, this process is an extremely challenging problem for modeling and simulation. The SPS process is inherently multiphysics and multi-scale, with the macroscale electro-thermo-mechanical behavior linked intricately to the microstructure evolution of the part being sintered. Accurate modeling and simulation tools, specifically geared towards the SPS process, are needed to predict the influence of the multiple variables involved in the manufacturing process. Modeling and simulation accuracy is achieved and demonstrated through comparison to multiple validation experiments. The validation efforts for Freya include both separate effects and complete multiphysics SPS process experiments. One of the key benefits this Laboratory Directed Research & Development (LDRD) project offers stems from the emphasis placed on experimental validation of the Freya models, both on the individual length scales and of the final coupled multiscale multiphysics simulations. Experimental validation of Freya’s multiscale coupling capability provides the technical credibility necessary for potential future industry and research partners to accept the simulation predictions.

36 MATERIALS SCIENCE↗

VOTCA: multiscale frameworks for quantum and classical simulations in soft matter

Many physical phenomena in liquids and soft matter are multiscale by nature and can involve processes with quantum and classical degrees of freedom occurring over a vast range of length- and timescales. Examples range from structure formation processes of complex polymers or even polymer blends (Svaneborg & Everaers, 2023) on the classical side to charge and energy transport and conversion processes (Lee et al., 2019) involving explicit electronic and, therefore, quantum information. The Versatile Object-oriented Toolkit for Coarse-graining Applications (VOTCA) provides multiscale frameworks built on a comprehensive set of methods for the development of classical coarse-grained potentials (VOTCA-CSG) as well as state-of-the art excited state electronic structure methods based on density-functional and many-body Green’s function theories, coupled in mixed quantum-classical models and used in kinetic network models (VOTCA-XTP).

97 MATHEMATICS AND COMPUTING↗