Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary condition”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Coming Closer to High Frequency Gravitational Wave Detection With MAGO

In the last years, low frequency gravitational waves (GWs) have been consistently measured by the LIGOVirgo collaboration, but little to no attention has been paid to higher frequencies GWs in the range of 10 kHz to 100 MHz, at which confirmation for current theories or even new physics could be hidden. The MAGO 2.0 project aims at filling this gap in the parameter space using superconducting radio-frequency (SRF) cavities. Exploiting the excellent Q-factors of these resonators, we plan to detect tiny harmonic deformations induced by GWs which change the boundary conditions of the oscillating electromagnetic field. We present the results of the first cold tests ran at DESY and FNAL using the cavity prototype built 20 years ago at the end of the MAGO collaboration, characterizing the RF spectrum, Q-factor and surface resistance. Additionally we introduce the mechanical vibration spectrum characterization and the RF response of the cavity with the injection of a “fake GW” signal using piezoelectric actuators.

Marconato, Giovanni [Hamburg U.]↗
DOE OSTI2025

Surrogate-driven Variance-based Sensitivity Analysis of Thermal Storage Tanks in Integrated Energy Systems

Sensitivity analysis and uncertainty quantification are essential steps for enhancing the accuracy of computational models by identifying and mitigating uncertainties. This study focuses on these steps for the Thermal Energy Delivery System at Idaho National Laboratory, specifically targeting the thermocline tank. Using a Modelica/Dymola simulation model, the study perturbed various design parameters and boundary conditions, including shape factor, porosity, outlet temperature, inlet mass flow rate, and system pressure, to predict and quantify uncertainty in the tank’s ax- ial temperature. A dataset of over 1,000 simulations was generated, and surrogate models were developed using the pyMAISE (Michigan Artificial Intelligence Standard Environment) library, which is an Automatic Machine Learning library for nuclear engineering applications. The optimal model, a feedforward neural network with two hidden layers, achieved an R2 score above 0.99 and a mean absolute error below 1 Kelvin. Sensitivity analyses using Sobol indices and Fourier amplitude sensitivity testing methods on this surrogate model revealed that the inlet mass flow rate at initial timestamps and porosity significantly impacts predicted temperatures across all sensors and time steps.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Physics and Components syntax to enable a systems-based approach to multiphysics

Simulations in MOOSE have traditionally used kernel and boundary condition classes to describe the equations. Downstream applications leveraged a system called Actions to define a pre-packaged discretization of the equations they solve. Unfortunately, the Action base class was very limited, and most applications implemented the same concepts in their Actions. This led to an increased maintenance burden and a reduction in coupling opportunities, save for the use of MultiApps which renders each input mostly independent. With the introduction of multi-system capabilities in MOOSE, there is growing interest in defining entire simulations of complex multiphysics systems in a single input file. By introducing a new Physics system, with its dedicated syntax and a new base class providing wide-ranging capabilities, we are now able to define multiple equations in a single input file in a compact and user-friendly way. With new interactions between Physics and the Component system, these equations can be defined on each component of a complex system. In this talk, we will present the capabilities of these new systems, their interactions, and how to define complex systems multiphysics simulations with Physics and Components.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

MOOSE Web Server Interface: A Message-based External Interface for Multiphysics Simulations

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a C++ toolkit designed to streamline the development of finite element and finite volume applications. It offers an interface for input-based coupling of these applications to create multiscale, multiphysics models. We introduce a new capability that enables external applications to integrate with MOOSE-based applications in situ via a web server using HTTP requests. An example of this integration is provided, where a MOOSE thermal-fluids solve has a boundary condition that is driven by an external Python application. Additionally, the coupling of the Python-based OpenMC depletion solver with the Cardinal application is demonstrated. A multiphysics model of a pressurized water reactor, incorporating neutronics, heat conduction, thermal-fluids, and depletion, is presented to showcase this new Cardinal capability that is enabled by the MOOSE web server capability.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

MOOSE Web Server Interface: A Message-based External Interface for Multiphysics Simulations

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a C++ toolkit designed to streamline the development of finite element and finite volume applications. It offers an interface for input-based coupling of these applications to create multiscale, multiphysics models. We introduce a new capability that enables external applications to integrate with MOOSE-based applications in situ via a web server using HTTP requests. An example of this integration is provided, where a MOOSE thermal-fluids solve has a boundary condition that is driven by an external Python application. Additionally, the coupling of the Python-based OpenMC depletion solver with the Cardinal application is demonstrated. A multiphysics model of a pressurized water reactor, incorporating neutronics, heat conduction, thermal-fluids, and depletion, is presented to showcase this new Cardinal capability that is enabled by the MOOSE web server capability.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Concept of Dynamic Heat Insulation for Rotating Detonation Engines

This work introduces a new class of materials concept to dynamically reduce instantaneous heat fluxes in Rotating Detonation Engine (RDE) combustor chamber walls. The high-frequency and high-amplitude surface heat fluxes observed in RDEs arise from large instantaneous temperature differences between the detonation shockwave and chamber wall surface. These temperature gradients drive substantial energy loss and reduce the chamber gas pressure, ultimately limiting the cycle’s thermodynamic efficiency. This work introduces a concept for dynamically insulating the combustion chamber surfaces using surface layers or coatings with low thermal time scale. With such coatings, the surface temperature may follow the fluctuations of the cyclic detonation wave temperature, thus reducing the instantaneous heat flux therefore cycle-mean heat flux. To analyze these cyclic thermal phenomena, a one-dimensional analytical conduction solver was utilized with the capability to handle multilayered structures. Parametric modeling was performed using transient heat flux boundary conditions representative of a hydrogen–air RDE across a broad range of coating thermal properties and engine conditions. The coating effectiveness scaled with the product of thermal time constant and detonation wave frequency and the results were non-dimensionalized to guide future materials development. This strategy may offer benefits in increasing material survivability, reducing cooling requirements, and enhancing pressure gain.

heat transfer↗
DOE OSTI2025

A Domain-Decomposed A-ϕ Formulation Based on Lagrange Multipliers for Low-Frequency Problems

A domain-decomposed A-ϕ formulation based on Lagrange multipliers is proposed to simulate low-frequency elec- tromagnetic problems. This method partitions the computational domain into smaller subdomains, allowing each subdomain to be independently formulated using Lagrange multipliers as Dirichlet boundary conditions, while ensuring continuity of the fields across the interfaces. A mixed finite element method, utilizing both vector and scalar basis functions, is employed to discretize the formulation, resulting in a global system to be solved. The proposed method is validated using TEAM Problem 7 at 50 Hz, demonstrating its effectiveness in handling complex geometries and addressing the low-frequency breakdown issues commonly encountered in traditional finite element methods.

Hossain, Amzad↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2026

Thermal Performance of Spandrel Assemblies in Glazed Wall Systems: Laboratory Test Design – Challenges and Test Results

Accurate thermal performance calculation procedures for opaque spandrel areas in curtain wall and window wall systems are essential for rating systems when comparing spandrel systems. However, there is a lack of consensus in thermal modeling needed for accurately characterizing heat transfer through spandrel assemblies due to the complex arrangement of materials and structural components. Several studies indicate that conventional 2D thermal simulations may overestimate R-values by 30% compared to physical testing and 3D simulations. Detailed simulations and well-curated laboratory test data are necessary to build confidence in simulation models, which will later be used to develop correlations to improve widely used conventional 2D thermal simulations. This study aims to experimentally test heat transfer through various spandrel assemblies to validate 3D simulation models. Also, the challenges of conducting a thorough testing design along with the solutions would be documented. The team developed a design for testing spandrel assemblies, making appropriate modifications to the existing heat, air, and moisture (HAM) chamber to accommodate the testing needs. Two moveable baffles were designed and fabricated to guide airflow direction parallel to the test article surface. The data acquisition capabilities in the chamber were upgraded to add more than two hundred sensors to the climate and indoor side of the chamber. The goal is to provide a quality dataset for validating complex 3D modeling simulations, which will be used to develop improved thermal simulation techniques that more accurately represent the thermal behavior of spandrel assemblies and their integration within the building envelope. This paper will summarize the results for the boundary conditions of the testing and the temperature variation across different locations of the spandrel assemblies.

Kunwar, Niraj [ORNL] (ORCID:0000000263457652)↗
DOE OSTI2020

Large-eddy simulations of the Northeastern US coastal marine boundary layer

In this study, large eddy simulations (LES) of offshore boundary layers near the Nantucket coast are performed using Nalu-Wind. The marine boundary layer conditions are chosen to match the predominantly unstable and neutral conditions measured by the Cape Wind platform. The appropriate domain, resolution, and boundary condition settings required for the LES are established through this work. Differences between stable and unstable cases are found in the wind speed profiles, averaged statistics, and wind spectra, and explained in terms of stratification effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Boundary layer moisture variability at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic observatory during marine conditions

Abstract. Boundary layer moisture variability at the Eastern North Atlantic (ENA) site during marine conditions is examined at monthly and daily timescales using 5 years of ground-based observations and output from the European Center for Medium range Weather Forecast (ECMWF) reanalysis model. The annual cycle of the mixed-layer total water budgets is presented to estimate the relative contribution of large-scale advection, local moisture tendency, entrainment, and precipitation to balance the moistening due to surface latent heat flux on monthly timescales. When marine conditions prevail, advection of colder and dry air from the north acts as an important moisture sink (∼ 50 % of the overall budget) during fall and winter driving the seasonality of the budget. Entrainment and precipitation contribute to the drying of the boundary layer (∼ 25 % and ∼ 15 % respectively), and the local change in moisture contributes to a small residual part. On a daily temporal scale, moist and dry mesoscale columns of vapor (∼ 10 km) are analyzed during 10 selected days of precipitating stratocumulus clouds. Adjacent moist and dry columns present distinct mesoscale features that are strongly correlated with clouds and precipitation. Dry columns adjacent to moist columns have more frequent and stronger downdrafts immediately below the cloud base. Moist columns have more frequent updrafts, stronger cloud-top cooling, and higher liquid water path and precipitation compared to the dry columns. This study highlights the complex interaction between large-scale and local processes controlling the boundary layer moisture and the importance of spatial distribution of vapor to support convection and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗
DOE OSTI2021

Study of Stratus-Lowering Marine-Fog Events Observed During C-FOG

Two stratus lowering marine fog events observed on 28 September and 04 October 2018 during the Coastal Fog (C-Fog) field campaign that took place offshore of Eastern Canada during 1 September to 6 October 2018 are described. In-situ, profiling and remote sensing observations were made at selected land sites in eastern Newfoundland (NL) and Nova Scotia (NS) as well as aboard the research vessel (R/V) Hugh R. Sharp that cruised in adjoining coastal waters. Synoptic scale analysis showed that both fog episodes were an outcome of the interaction between synoptic-scale surface level low-pressure systems and a contiguous high-pressure system. At the same time, back trajectories revealed that the bulk of the fog layer is formed due to differential advection. Here, the diameter of the fog droplets at the surface gradually decreased from the centre of the fog layer to its leading/trailing edges. The bimodal fog droplet diameter distribution with peaks at 5-10 µm and 20-25 µm provided clues on droplet collision and coalescence processes. The observed difference between microphysical variables and droplet distribution between the two fog events and within same fog layer may have been governed by atmospheric boundary layer conditions (e.g., humidity conditions and turbulence) prevailed in the fog layer. Overall, it is concluded that the life cycle of observed stratus-lowering coastal-fog episodes is dependent on synoptic conditions as well as atmospheric boundary layer characteristics such as stability, cloud top cooling and entrainment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Methodology to determine the coupling of continental clouds with surface and boundary layer height under cloudy conditions from lidar and meteorological data

Abstract. The states of coupling between clouds and surface or boundary layer have been investigated much more extensively for marine stratocumulus clouds than for continental low clouds, partly due to more complex thermodynamic structures over land. A manifestation is a lack of robust remote sensing methods to identify coupled and decoupled clouds over land. Following the idea for determining cloud coupling over the ocean, we have generalized the concept of coupling and decoupling to low clouds over land, based on potential temperature profiles. Furthermore, by using ample measurements from lidar and a suite of surface meteorological instruments at the U.S. Department of Energy's Atmospheric Radiation Measurement Program's Southern Great Plains site from 1998 to 2019, we have developed a method to simultaneously retrieve the planetary boundary layer (PBL) height (PBLH) and coupled states under cloudy conditions during the daytime. The new lidar-based method relies on the PBLH, the lifted condensation level, and the cloud base to diagnose the cloud coupling. The coupled states derived from this method are highly consistent with those derived from radiosondes. Retrieving the PBLH under cloudy conditions, which has been a persistent problem in lidar remote sensing, is resolved in this study. Our method can lead to high-quality retrievals of the PBLH under cloudy conditions and the determination of cloud coupling states. With the new method, we find that coupled clouds are sensitive to changes in the PBL with a strong diurnal cycle, whereas decoupled clouds and the PBL are weakly related. Since coupled and decoupled clouds have distinct features, our new method offers an advanced tool to separately investigate them in climate systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Passive Safety System Reliability Analysis: Methodology, Application, and Results

Many advanced reactor designs rely on passive systems to provide assurance that certain safety functions, such as reactivity control and reactor heat removal, are maintained during transient scenarios. The performance of these systems is sensitive to boundary and environmental conditions and can lead to functional failure, which is the inability of the system to perform as intended even without any physical failures. Additionally, conditions can cause the passive systems to operate in an intermediate or degraded state. These operational modes of passive systems and their overall reliability can be difficult to capture in traditional probabilistic frameworks where discrete operating modes are assumed or time-dependent boundary conditions are not accounted for. Novel strategies and methods are needed to be able to capture the reliability of these passive systems and their impact on facility safety. The current abstract provides a summary of such a methodology and examples of its use in a case study, the GE Hitachi (GEH) PRISM probabilistic safety assessment (PSA), and the U.S. Department of Energy (USDOE) Versatile Test Reactor (VTR) PSA.

Bucknor, M.↗
DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗