Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Little influence of Arctic amplification on mid-latitude climate

Observations and model simulations show enhanced warming in the Arctic under increasing greenhouse gases, a phenomenon known as the Arctic amplification (AA), that is likely caused by sea-ice loss. AA reduces meridional temperature gradients linked to circulation, thus mid-latitude weather and climate changes have been attributed to AA, often on the basis of regression analysis and atmospheric simulations. However, other modelling studies show only a weak link. This inconsistency may result from deficiencies in separating the effects of AA from those of natural variability or background warming. In this work, using coupled model simulations with and without AA, we show that cold-season precipitation, snowfall and circulation changes over northern mid-latitudes come mostly from background warming. AA and sea-ice loss increase precipitation and snowfall above ~60° N and reduce meridional temperature gradients above ~45° N in the lower–mid troposphere. However, minimal impact on the mean climate is seen below ~60° N, with weak reduction in zonal wind over 50°–70° N and 150–700 hPa, mainly over the North Atlantic and northern central Asia. These results suggest that the climatic impacts of AA are probably small outside the high latitudes, thus caution is needed in attributing mid-latitude changes to AA and sea-ice loss on the basis of statistical analyses that cannot distinguish the impact of AA from other correlated changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Transient Analysis of a Micro-reactor using the DireWolf Code Suite

Transient analyses of heat pipe micro-reactors are necessary to ensure that hypothetical accident scenarios do not comprise reactor safety. Due to its small size and reliance on heat-pipes for cooling, the micro-reactor design introduced in this paper is a tightly coupled system which requires multi-physics tools to accurately model transient events. Idaho National Laboratory’s DireWolf code suite based on the MOOSE framework is tailor-built to model heat-pipe reactors. This paper demonstrates DireWolf’s ability to simulate the coupled thermal-neutronics transient behavior of a heat-pipe micro-reactor. The transient events presented here include an inadvertent rotation of all control drums simultaneously and a sudden complete ejection of a single control drum. A detailed description of each event is provided along with simulation results, including time dependent power and temperature distributions, and discussion. This is a Westinghouse Electric Company (WEC) led publication.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Insights on Simulating Summer Warming of the Great Lakes: Understanding the Behavior of a Newly Developed Coupled Lake-Atmosphere Modeling System

The Laurentian Great Lakes are the world's largest freshwater system and regulate the climate of the Great Lakes region, which has been increasingly experiencing climatic, hydrological, and ecological changes. An accurate mechanistic representation of the Great Lakes thermal structure in Regional Climate Models (RCMs) is paramount to studying the climate of this region. Currently, RCMs have primarily represented the Great Lakes through coupled one-dimensional (1D) column lake models; this approach works well for small inland lakes but is unable to resolve the realistic hydrodynamics of the Great Lakes and leads to inaccurate representations of lake surface temperature (LST) that influence regional climate and weather patterns. This work overcomes this limitation by developing a fully two-way coupled modeling system using the Weather Research and Forecasting model and a three-dimensional (3D) hydrodynamic model. The coupled model system resolves the interactive physical processes between the atmosphere, lake, and surrounding watersheds; and validated against a range of observational data. The model is then used to investigate the potential impacts of lake-atmosphere coupling on the simulated summer LST of Lake Superior. By evaluating the difference between our two-way coupled modeling system and our observation-driven modeling system, we find that coupled-lake atmosphere dynamics can lead to a higher LST during June-September through higher net surface heat flux entering the lake in June and July and a lower net surface heat flux entering the lake in August and September. The unstratified water in June distributes the entering surface heat flux throughout the water column leading to a minor LST increase, while the stratified waters of July create a conducive thermal structure for the water surface to warm rapidly under the higher incoming surface heat flux. This research provides insight into the coupled modeling system behavior, which is critical for enhancing our predictive understanding of the Great Lakes climate system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Electrode Mesoscale as a Collection of Particles: Coupled Electrochemical and Mechanical Analysis of NMC Cathodes

Battery electrodes are composed of polydisperse particles and a porous, composite binder domain. These materials are arranged into a complex mesostructure whose morphology impacts both electrochemical performance and mechanical response. We present image-based, particle-resolved, mesoscale finite element model simulations of coupled electrochemical-mechanical performance on a representative NMC electrode domain. Beyond predicting macroscale quantities such as half-cell voltage and evolving electrical conductivity, studying behaviors on a per-particle and per-surface basis enables performance and material design insights previously unachievable. Voltage losses are primarily attributable to a complex interplay between interfacial charge transfer kinetics, lithium diffusion, and, locally, electrical conductivity. Mesoscale heterogeneities arise from particle polydispersity and lead to material underutilization at high current densities. Particle-particle contacts, however, reduce heterogeneities by enabling lithium diffusion between connected particle groups. While the porous composite binder domain (CBD) may have slower ionic transport and less available area for electrochemical reactions, its high electrical conductivity makes it the preferred reaction site late in electrode discharge. Mesoscale results are favorably compared to both experimental data and macrohomogeneous models. This work enables improvements in materials design by providing a tool for optimization of particle sizes, CBD morphology, and manufacturing conditions.

42 ENGINEERING↗
DOE OSTI2020

CyDER: A Cyber Physical Co-simulation Platform for Distributed Energy Resources in Smartgrids

The CyDER project aimed at developing an open-source, modular and scalable co-simulation platform for power grids with large shares of Distributed Energy Resources (DERs). The project partners are the Lawrence Berkeley National Lab (LBNL), Lawrence Livermore National Lab (LLNL), PG&E, SolarCity, and ChargePoint. The prime recipient is LBNL; SolarCity and ChargePoint were partners for the project’s first two years. Increased DER integration introduces a number of challenges in power grid operation including a more dynamic interaction between the transmission grid and distribution grids, and increased modeling complexity. Although specialized software exists to precisely model different components of the power system, it is far from trivial to integrate all various models and perform a holistic simulation. Instead of replicating all models in a common simulation program, a commonly accepted approach to tackle this model diversity is to couple third-party simulators and models through a co-simulation platform that coordinates information exchange among the various components. Following this line of research, this project’s objective was to develop a co-simulation platform based on a widely accepted industrial standard called Functional Mockup Interface (FMI). Within this process, the project developed models compliant with the FMI standard, called Functional Mockup Units (FMUs), and used them to perform various operational and planning power system analyses. Relying and building upon an industrial standard is the main differentiation of this project compared with previous or parallel efforts in the co-simulation area. Particular emphasis was put on delivering software utilities to facilitate setting up and running co-simulations by end-users. Furthermore, a strong aspect of this project is demonstrating that co-simulation techniques can be used to perform Hardware-in-the-Loop (HIL) simulations that couple software components (e.g., simulated models) with hardware components (e.g., real devices such PV systems and batteries). The long-term goal of CyDER project is to help establish FMI as a powerful standard for co-simulation and promote adoption by electric utilities and other interested stakeholders. The main accomplishments of the project include the development of several FMUs including distribution and transmission grid models, PV inverters with Volt/Var/Watt controllers, batteries, and predictive optimal controllers. Additionally, a unique software package was developed, called SimulatorToFMU, which is capable of exporting any Python-driven simulator or Python script as an FMU. This is an important contribution towards establishing FMI as one of the main co-simulation standards, because more and more third-party programs for sub-system modeling and simulation are delivered with Python APIs. The CyDER platform was used to perform PV hosting capacity analyses in real utility feeders with and without smart inverter controls, battery storage, and EV charging. Smart inverter controls include conventional Volt/Var/Watt controls for reactive power support and active power curtailment, but also predictive controls that optimize the charging and discharging profile of the battery connected on the DC side in order to minimize the customer’s economic benefit. Finally, an important result of this project is delivering an experimental setup that consists of residential-scale PV inverters with battery storage, a real-time grid simulator with an ideal voltage source as grid emulator, and micro Phasor Measurement Units (PMUs). All these components and additional software modules are coupled to one another using the FMI standard and can be co-simulated with the CyDER platform.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Coupled SAM/Griffin Model of a Reference Pebble Bed High-Temperature Gas Cooled Reactor for Multi-Physics Simulations

An effort has been dedicated to developing a reference model for multi-physics coupled simulation of the pebble bed high-temperature gas-cooled reactor (PB-HTGR) with SAM and Griffin computer codes for safety analysis and licensing purpose. 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, a SAM model of the HTR-PM reactor has been developed based on publicly available design information and the multi-dimensional Pronghorn model developed by Idaho National Laboratory. The SAM HTR-PM model includes a multi-dimensional core region and 0- D/1-D fluid components. Additionally, a simplified air RCCS loop is modeled for decay heat removal. The Griffin model is based on a recent work by Idaho National Laboratory. The coupling between the models is achieved through the MOOSE MultiApp system. Both steady-state and transient scenarios were simulated to demonstrate the coupled model’s capability for multi-physics simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Intrinsic bulk viscosity of the one-component plasma

The intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at Γ ≈ 1, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average interparticle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: Γ ≈ 10 –2 –10 2 . Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. As a result, the generalized shear and bulk viscosity coefficients are found to have a similar magnitude for ω ≳ 2⁢ω p at strongly coupled conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivity of Simulated Polar Climate to Improved Partitioning of Spectral Solar Irradiance Between Visible and Near‐Infrared Bands

The solar radiative processes that contribute to Earth's surface and atmospheric energy budgets are strongly dependent on wavelength. For example, snow and water vapor become more absorptive as incident radiation shifts from visible (VIS) to near-infrared (NIR) wavelengths. Some earth system models (ESMs) aggregate solar radiation into just two bands (VIS and NIR) partitioned at 0.7 μm for transmission between the atmospheric and surface components. In the widely used radiative transfer model RRTMG_SW this partition is near the center of the overlap spectral band spanning 0.625–0.778 μm, whose flux is often approximated as being evenly divided between the VIS and NIR surface bands. Using a hyperspectral radiative transfer model, we show that the fractional downwelling surface flux within the overlap band is usually distributed about 55.5%:44.5% VIS:NIR. This improved approximation shifts as much as 4.27 W m −2 from the NIR to the VIS band, leading to an instantaneous decrease in surface absorption of up to 0.91 W m −2 over snow-covered surfaces. Century-long fully coupled ESM simulations show surface absorption over snow decreases by over 1.6 W m −2 . The coupled response in sea ice regions amplifies the initial forcing ten-fold, and increases seasonal sea ice area by up to 4.9%. These results highlight the importance of accurately representing the spectral distribution of solar radiation in the cryosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assembling Multiphysics Nuclear Reactor Simulations Using the MOOSE Framework

The Multiphysics Object Oriented Simulation Environment (MOOSE) [1] is an open-source, parallel finite element framework which provides the foundation for many advanced modeling and simulation tools developed under the Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program [2] for the analysis of advanced reactors. The MOOSE framework provides the common foundational capability on which many NEAMS codes for reactor analysis are built. The MOOSE framework also includes several systems to assemble unique workflows and couplingamong MOOSE-based applications. In particular, the MultiApp and Transfer Systems are widely used to assemble different MOOSE-based or MOOSE-wrapped physics applications together to perform loosely or tightly coupled multiphysics simulations. The National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB) [3] hosts publicly available nuclear reactor multiphysics simulation examples which leverage MOOSE’s MultiApp System to meet the modeling needs of different reactor types. The flexibility and robustness of coupling provided by MOOSE permits rapid development of coupled physics models for a wide range of reactor types and events

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A process-based evaluation of biases in extratropical stratosphere–troposphere coupling in subseasonal forecast systems

Abstract. Two-way coupling between the stratosphere and troposphere is recognized as an important source of subseasonal-to-seasonal (S2S) predictability and can open windows of opportunity for improved forecasts. Model biases can, however, lead to a poor representation of such coupling processes; drifts in a model's circulation related to model biases, resolution, and parameterizations have the potential to feed back on the circulation and affect stratosphere–troposphere coupling. We introduce a set of diagnostics using readily available data that can be used to reveal these biases and then apply these diagnostics to 22 S2S forecast systems. In the Northern Hemisphere, nearly all S2S forecast systems underestimate the strength of the observed upward coupling from the troposphere to the stratosphere, downward coupling within the stratosphere, and the persistence of lower-stratospheric temperature anomalies. While downward coupling from the lower stratosphere to the near surface is well represented in the multi-model ensemble mean, there is substantial intermodel spread likely related to how well each model represents tropospheric stationary waves. In the Southern Hemisphere, the stratospheric vortex is oversensitive to upward-propagating wave flux in the forecast systems. Forecast systems generally overestimate the strength of downward coupling from the lower stratosphere to the troposphere, even as most underestimate the radiative persistence in the lower stratosphere. In both hemispheres, models with higher lids and a better representation of tropospheric quasi-stationary waves generally perform better at simulating these coupling processes.

Garfinkel, Chaim I. (ORCID:000000017258666X)↗
DOE OSTI2021

Assessment of Fast Reactor Hot Channel Factor Calculation Capability in Griffin and NekRS

The DOE-NE Advanced Modeling and Simulation (NEAMS) program develops, demonstrates, and deploys predictive computational modeling tools for nuclear reactor design and analysis. Under the Fast Reactor Applications activity within the Application Drivers Technical Area of NEAMS, the ability to perform high fidelity “hot channel factor” (HCF) coupled multiphysics simulations has been identified as a valuable simulation capability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Multiphysics analysis system for heat pipe cooled micro-reactors employing PRAGMA-OpenFOAM-ANLHTP

A multiphysics analysis system for neutronics/thermo-mechanical/heat-pipe analysis of heat pipe cooled micro-reactors was developed using the PRAGMA code as the neutronics engine. PRAGMA, which used to be a GPU-based continuous-energy MC code for power reactor applications, now has an extended geometry package to handle geometries with unstructured meshes generated by the ANSYS Design-Modeler and Meshing. The NVIDIA ray tracing engine OptiX was exploited for efficient neutron transport on unstructured geometry. On the multiphysics side, the open-source CFD tool OpenFOAM and one-dimensional heat pipe analysis code ANLHTP were adopted. The manager-worker system based on the MPI dynamic process management (DPM) model enables efficient coupling of codes employing different parallelization schemes. With all the features, the multiphysics analysis of the one-sixth symmetrical MegaPower 2D core was performed and it demonstrated the benefit of the tight integration of the three-way coupling system and one-to-one geometry coupling strategy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

pystorms: A simulation sandbox for the development and evaluation of stormwater control algorithms

Advances in cyber–physical technologies have enabled real-time sensing and adaptive control of stormwater infrastructure. These smart stormwater systems allow for inexpensive, minimally-invasive stormwater control interventions in lieu of new construction. Importantly, however promising the area of smart stormwater control, there still remain barriers – for experts and novices alike – to access shared tools and methods for investigating, developing, and contributing to it. In an effort to make smart stormwater control research more methodical and accessible, we present pystorms, an open-source Python-based simulation sandbox that facilitates the quantitative evaluation and comparison of control strategies. pystorms consists of a collection of real world-inspired smart stormwater control scenarios on which any number of control strategies can be applied and tested via an accompanying Python programming interface and coupled stormwater simulator. pystorms provides a framework for the rigorous and efficient evaluation of smart stormwater control methodologies across diverse watersheds with only a few lines of code.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessment of Lost Circulation Material Particle-Size Distribution on Fracture Sealing: A Numerical Study

Summary Lost circulation materials (LCMs) are essential to combat fluid loss while drilling and may put the whole operation at risk if a proper LCM design is not used. The focus of this research is understanding the function of LCMs in sealing fractures to reduce fluid loss. One important consideration in the success of fracture sealing is the particle-size distribution (PSD) of LCMs. Various studies have suggested different guidelines for obtaining the best size distribution of LCMs for effective fracture sealing based on limited laboratory experiments or field observations. Hence, there is a need for sophisticated numerical methods to improve the LCM design by providing some predictive capabilities. In this study, computational fluid dynamics (CFD) and discrete element methods (DEM) numerical simulations are coupled to investigate the influence of PSD of granular LCMs on fracture sealing. Dimensionless variables were introduced to compare cases with different PSDs. We validated the CFD-DEM model in reproducing specific laboratory observations of fracture-sealing experiments within the model boundary parameters. Our simulations suggested that a bimodally distributed blend would be the most effective design in comparison to other PSDs tested here.

Engineering↗
DOE OSTI2023

Natural Fracture Pattern Development Using 3D Coupled Numerical Model: Mechanics and Diagenesis

In the earth's crust, while tectonic processes can mechanically create natural fractures, cementation can simultaneously grow on the surfaces of the fractures. This crystallization can affect the growth of the fractures. To more accurately address the modeling of natural fracture development, both mechanical deformation and chemical diagenesis should be carried out in an integrated manner. Here, we achieve this with a three-dimensional (3D) simulator that couples mechanical and chemical processes for natural fracture pattern development in layer confined situations. Tectonic loading over geologic time is mimicked by displacement driven boundary conditions in a 3D displacement discontinuity method (DDM) scheme. Quartz diagenesis competes with mechanical fracture opening, resulting in three different cementation conditions - open, bridging or fully infilled – which define fracture element stiffness that is used as feedback within the 3D-DDM code during propagation. An important milestone of this work is that this coupled structural diagenetic model can reproduce power-law aperture population statistics which cannot be generated using mechanics alone. More accurately estimated mechanical fracture apertures and their cementation condition is a critical step toward better permeability estimation and fluid path analysis for subsurface reservoirs.

58 GEOSCIENCES↗
DOE OSTI2024

Anthropogenic effects on flood hazards in a hyper-arid watershed: The 2015 Atacama floods

An unprecedented precipitation event in the hyper-arid Atacama Desert of Northern Chile occurred in March 2015. Geomorphic alterations to the river channel and the coastal zone, coupled with the exceptional magnitude of the rainfall, caused catastrophic damage and loss of life. On the coast of the El Salado watershed, legacy mine tailings infilled the watershed-ocean connection, while the river channel was altered both by tailings and urbanization. The consequences of this event resulted from the coupling of anthropogenic geomorphic changes with an unusual climate event. Lack of field data, complex geomorphology and sediment loads influenced by human activity make analysing floods in these regions especially challenging. The objective of this work is to improve our understanding of the factors that control flood hazards by using numerical simulations to reconstruct the 2015 flood in El Salado. We carry out unsteady two-dimensional simulations fully coupled with the sediment concentration to identify the influence of tailing deposits, considering high-resolution data of the pre- and post-2015 flood topography. In conclusion, the results highlight the importance of specific event-based studies, using models that can help designing better strategies for climate change adaptation and risk mitigation, while providing information for risk reduction and channel restoration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The influence of permeability anisotropy in the upper ocean crust on advective heat transport by a ridge-flank hydrothermal system

Here, in this study, we highlight the importance of permeability anisotropy on the hydrogeological regime of a ridge-flank hydrothermal system. Our study site, North Pond, is a marine sediment pond on ~8 Ma seafloor in the North Atlantic, and represents a low-temperature, end-member ridge-flank hydrothermal system. Previous simulations of North Pond elucidated long-standing hypotheses concerning hydrothermal fluid and heat transport in the upper volcanic crust but failed to fully explain observed patterns of seafloor heat flux in this area. Here we use variography, a geostatistical method, to quantify relations between seafloor heat-flux measurements, and coupled numerical simulations of fluid and heat flow to simulate the hydrogeologic regime. Directional variography shows that heat-flux observations are correlated along-strike of the regional crustal fabric. Three-dimensional simulations that include permeability anisotropy are able to replicate seafloor heat-flux patterns across North Pond. The simulations that result in the best match to thermal data incorporate permeability anisotropy in the horizontal plane. We find that the feedback between permeability anisotropy and the asymmetric geometry of North Pond combine to promote advective removal of heat and mass within the crustal aquifer. These findings suggest that permeability anisotropy in the oceanic crust may influence ridge-flank hydrothermal circulation more broadly.

58 GEOSCIENCES↗