Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional generative models”

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 109 records · Page 6

Extracellular matrix protein DMP1 suppresses osteogenic differentiation of Mesenchymal Stem Cells

Highlights: • Extracellular matrix protein DMP1 suppresses MSCs osteogenic differentiation. • DMP1 works not only in osteoblasts and osteocytes, but also in MSCs. Mesenchymal Stem Cells (MSCs) are self-renewing and multipotent stem cells which was investigated for diverse clinical applications. However, complex mechanism of MSCs fate determination is still not fully disclosed. Extracellular matrix (ECM) proteins contribute to maintain MSCs stemness by providing extracellular microenvironment. Increasing evidences show that ECM proteins could also regulate the fate of MSCs directly. Dentin matrix protein 1 (DMP1) is an ECM protein enrich in bone tissue and terminal cells, which well-known in promoting osteoblasts and osteocytes maturation, and facilitate mineralization. Recently, our experiment indicated that DMP1 was also expressed in MSCs of long bone. In present study, it is found that DMP1 expressed in Prx1 positive MSCs. And, DMP1 is down-regulated in early osteoblasts and up-regulated again in mature osteoblasts. DMP1 conditional knockout mice model under Prx1cre was generated to explore whether DMP1 regulates MSCs osteogenic differentiation. Specific ablation of DMP1 in Prx1 positive MSCs increased bone mass in vivo and promoted osteoblasts activity in vitro. This study provides a new understanding of DMP1's function in regulation of osteogenesis: not only an enhancer of bone formation, but also a negative regulator of MSCs differentiation in bone.

60 APPLIED LIFE SCIENCES↗

Improved subseasonal-to-seasonal precipitation prediction of climate models with nudging approach for better initialization of Tibetan Plateau-Rocky Mountain Circumglobal wave train and land surface conditions

Abstract Reliable subseasonal-to-seasonal (S2S) precipitation prediction is highly desired due to the great socioeconomical implications, yet it remains one of the most challenging topics in the weather/climate prediction research area. As part of the Impact of Initialized Land Temperature and Snowpack on Sub-seasonal to Seasonal Prediction (LS4P) project of the Global Energy and Water Exchanges (GEWEX) program, twenty-one climate models follow the LS4P protocol to quantify the impact of the Tibetan Plateau (TP) land surface temperature/subsurface temperature (LST/SUBT) springtime anomalies on the global summertime precipitation. We find that nudging towards reanalysis winds is crucial for climate models to generate atmosphere and land surface initial conditions close to observations, which is necessary for meaningful S2S applications. Simulations with nudged initial conditions can better capture the summer precipitation responses to the imposed TP LST/SUBT spring anomalies at hotspot regions all over the world. Further analyses show that the enhanced S2S prediction skill is largely attributable to the substantially improved initialization of the Tibetan Plateau-Rocky Mountain Circumglobal (TRC) wave train pattern in the atmosphere. This study highlights the important role that initial condition plays in the S2S prediction and suggests that data assimilation technique (e.g., nudging) should be adopted to initialize climate models to improve their S2S prediction.

54 ENVIRONMENTAL SCIENCES↗

Apparent Monte Carlo Source Convergence Problem with BWR Fuel Depleted with Partial Control Blade Insertion

Oak Ridge National Laboratory (ORNL) and the US Nuclear Regulatory Commission (NRC) have initiated a multiyear project to investigate application of burnup credit (BUC) to boiling-water reactor (BWR) fuel in storage and transportation casks and spent fuel pools (SFPs). Ideally, a single combination of all operating parameters including moderator density and axial burnup profiles, control blade exposure, and others can be identified to allow a simple, bounding approach to BWR BUC. One aspect of this investigation involves examining the effect of control blade insertion during operation on discharged assembly reactivity. In the course of these calculations, an apparent source convergence problem was discovered for fuel that has experienced partial blade insertion for the entire lifetime of the assembly. The source convergence problem appears not to impact k{sub eff}, but it manifests in predictions of the axial fission density profile. This phenomenon is of interest in identifying areas of concern with BWR BUC since it has not been studied as thoroughly as BUC for pressurized-water (PWR) reactor fuel. There may also be interest in this configuration for further study of a realistic source Monte Carlo convergence problem as opposed to some earlier manufactured source convergence problems through careful selection of initial conditions. This summary presents the depletion conditions used to generate the problem, a description of the models, a background based on well-behaved calculations, the identification of the problem, and attempts made to resolve the problem. Only a brief description of the depletion conditions and models is possible within the limited scope of this summary. Finally, conclusions are presented. This summary identifies a source convergence problem in the process of investigating BUC for BWR spent fuel. The depletion scenario which created the depleted fuel and burnable absorber compositions included a control blade inserted in the bottom 14 nodes throughout the entire depletion. This insertion is slightly over half the length of the assembly, and it appears to create two regions of nearly equal reactivity. As expected, one region is near the top of the assembly, and the other region is near the middle of the assembly. This node includes significant {sup 239}Pu generation as a result of rodded operation, but it is also the top node of the full lattice zone. While the calculated k{sub eff} values appear to be statistically equivalent, the axial neutron production profile varies significantly among independent simulations of the same model. A range of standard approaches has been used to attempt to resolve the source convergence problem and determine the true axial neutron production profile. This profile is of interest in cask criticality safety analysis because it influences the interaction among assemblies loaded in the cask. The forthcoming SCALE 6.2 release will include source convergence testing based on Shannon entropy, which may assist in identifying these types of problems in the future. It should be noted that this source convergence problem results from the depletion conditions and the final distribution of fissile and absorber material within the discharged fuel assembly. The results, while entirely generated within the SCALE code system, are representative of any Monte Carlo program attempting to calculate k{sub eff} for such a system. The most important conclusion of this work is that the BUC criticality safety community must be aware of the potential for this type of source convergence problem. Analysts must approach such partially rodded scenarios with caution and critically examine the results. Fortunately, it is extremely unlikely that BWR assemblies would experience the bladed operation necessary to create this particular source convergence problem. These results may also have some bearing on the analysis of PWR BUC for new plant designs such as AP1000 that incorporate significant rodded operations. In a broader context, the source convergence problem identified here could have been magnified if it had occurred in a Monte Carlo depletion simulation. It may be possible to identify this type of problem in a depletion calculation by comparing the apparently converged fluxes between steps to ensure that no oscillations or other serious shifts are occurring. The depletion calculations performed in this work were performed with a two-dimensional deterministic code and are therefore not subject to this problem. These results may also raise questions about the validity of 2D depletion models with extruded geometry in the 3D k{sub eff} calculation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computational Characterization and Boundary Condition Models of a Jet Interaction Fluidic Oscillator

Fluidically oscillating jet actuators that generate sweeping jets when supplied with a pressurized fluid have been used to mitigate separation and reduce drag in a range of flow control applications. Their implementation in future flight vehicles would require fundamental understanding and accurate predictive techniques of the physics of their internal flow and jet formation. The present investigations focus on high-fidelity, time-accurate simulations to characterize the flow physics of the actuator in quiescent conditions. An important element of the present simulations is to demonstrate the ability of the computational fluid dynamics (CFD) solver to predict the jet characteristics and provide a basis for the development of improved boundary conditions (BC) without entirely resolving the geometrical features of the fluidic device. The CFD-predicted oscillation frequencies of the engendered jets were found to be in excellent agreement with experiments, even on two-dimensional meshes. The study revealed that three-dimensional simulations are required to capture some of the flow features of the sweeping jet such as the double peak in time-averaged velocity distributions downstream of the actuator’s orifice that were measured in experiments. Furthermore, several approaches for modeling the actuator were implemented and assessed in quiescent conditions. The evaluation of a boundary condition at the device throat, based on the phase-averaged flow variables, provides the basis for devising surface-based boundary conditions. The influence and necessity of including turbulent characteristics as part of the boundary conditions have also been identified.

42 ENGINEERING↗

Multimodal imaging and machine learning to enhance microscope images of shale

A machine learning based image processing workflow is presented to enhance shale source rock microscopic images obtained using diverse imaging platforms. Images were acquired from a 30 μm diameter cylindrical Vaca Muerta shale sample using both nondestructive Transmission X-Ray Microscopy (TXM, alternately referred to as nano computed tomography) and destructive Focused Ion Beam-Scanning Electron Microscopy (FIB-SEM). Output cross-sectional images from each modality were aligned using a combination of manual and automated registration techniques to create a registered image dataset. We then apply this dataset for two image processing tasks: prediction of image cross sections with SEM-like resolution from nondestructive TXM data and repair of charged region artifacts (localized accumulation of electrons) within SEM images. The image processing algorithms for both tasks use deep learning models, specifically image-to-image Convolutional Neural Networks (CNNs) and conditional Generative Adversarial Networks (cGANs). In the image enhancement tasks, we are able to achieve significant qualitative and quantitative improvement in TXM images. Here, the best model reaches an average Peak Signal to Noise Ratio (PSNR) of 15.8 dB. Conditioning on TXM data is also shown to reduce artifacts from SEM charging, achieving an average PSNR of 25.8 dB. Furthermore, our results suggest that properly trained and validated networks are capable of significant enhancement of images obtained using nondestructive techniques, thereby improving interpretation of two- and three-dimensional images while preserving samples for future use.

58 GEOSCIENCES↗

GeneratorSE: A Sizing Tool for Variable-Speed Wind Turbine Generators

This report documents a set of analytical models employed by the optimization algorithms within the GeneratorSE framework. The initial values and boundary conditions employed for the generation of the various designs and initial estimates for basic design dimensions, masses, and efficiency for the four different models of generators are presented and compared with empirical data collected from previous studies and some existing commercial turbines. These models include designs applicable for variable-speed, high-torque application featuring direct-drive synchronous generators and low-torque application featuring induction generators. In all of the four models presented, the main focus of optimization is electromagnetic design with the exception of permanent-magnet and wire-wound synchronous generators, wherein the structural design is also optimized. Thermal design is accommodated in GeneratorSE as a secondary attribute by limiting the winding current densities to acceptable limits. A preliminary validation of electromagnetic design was carried out by comparing the optimized magnetic loading against those predicted by numerical simulation in FEMM4.2, a finite-element software for analyzing electromagnetic and thermal physics problems for electrical machines. For direct-drive synchronous generators, the analytical models for the structural design are validated by static structural analysis in ANSYS.

17 WIND ENERGY↗

Using probabilistic solar power forecasts to inform flexible ramp product procurement for the California ISO

How can independent system operators (ISOs) take advantage of probabilistic solar forecasts to lower generation costs and improve reliability of power systems? We discuss one three-step approach for doing so, focusing on how such forecasts might help the California Independent System Operator (CAISO) prepare unexpected net load ramps, where net load equals gross demand minus wind and solar production. First, we enhance an existing solar forecasting system to provide well-calibrated hours-ahead probabilistic forecasts. We then relate the degree of uncertainty reflected in the forecasted prediction intervals (independent variables) to error distributions for net load ramp forecasts for the CAISO real-time market (dependent variable) using machine learning and quantile regression. Projected ramp forecast errors conditioned on solar uncertainty are translated into flexible ramp requirements that therefore reflect real-time meteorological and solar conditions, improving on typical ISO procedures. Detailed descriptions are provided on the quantile regression and kth-nearest neighbor categorization methods for accomplishing that translation. Finally, a multiple time-scale look-ahead market simulation model is applied to a 118-bus IEEE Reliability Test System, modified to represent the CAISO generation mix and demand distributions. The model runs quantify how solar-conditioned ramp requirements can, first, decrease operating costs by reducing requirements compared to often conservative unconditional methods and, second, decrease generation scarcity events and consequently improve reliability by increasing flexibility requirements at times when unconditional forecast-based requirements understate actual ramp uncertainty. Solar-conditioned ramp requirements are found to reduce generation operating costs by about 2% for the test system (which would be equivalent to over $\$100$ million per year for a CAISO-size system).

14 SOLAR ENERGY↗

Structural constraint integration in a generative model for the discovery of quantum materials

Billions of organic molecules have been computationally generated, yet functional inorganic materials remain scarce due to limited data and structural complexity. Here, in this work, we introduce Structural Constraint Integration in a GENerative model (SCIGEN), a framework that enforces geometric constraints, such as honeycomb and kagome lattices, within diffusion-based generative models to discover stable quantum materials candidates. SCIGEN enables conditional sampling from the original distribution, preserving output validity while guiding structural motifs. This approach generates ten million inorganic compounds with Archimedean and Lieb lattices, over 10% of which pass multistage stability screening. High-throughput density functional theory calculations on 26,000 candidates shows over 95% convergence and 53% structural stability. A graph neural network classifier detects magnetic ordering in 41% of relaxed structures. Furthermore, we synthesize and characterize two predicted materials, TiPd 0.22 Bi 0.88 and Ti 0.5 Pd 1.5 Sb, which display paramagnetic and diamagnetic behaviour, respectively. Our results indicate that SCIGEN provides a scalable path for generating quantum materials guided by lattice geometry.

36 MATERIALS SCIENCE↗

Mesoscale to Microscale Coupling for a Wind Ramp Case over Complex Terrain

Within the wind energy community, it has become increasingly important to simulate realistic weather events at high fidelity in order to improve wind turbine design, wind farm efficiency, and power forecasting. One of the major difficulties in simulating atmospheric events with numerical models at high fidelity is in accurately transitioning from the mesoscale to the microscale. The mesoscale flow, when input into the large-eddy simulation (LES) domain, will be free of resolved turbulence; needing sufficient time and space (known as fetch) for resolved-scale turbulence to develop. In order to reduce this fetch, several methods exist in which the temperature field, wind field, or both, are stochastically perturbed in order to accelerate the generation of turbulence. However, for cases involving complex terrain, it remains unclear if these perturbation processes are necessary, and if the interaction of the flow over the terrain, alone, is sufficient to generate adequate turbulence. In this study, a wind ramp event observed at the National Wind Technology Center (NWTC) outside of Boulder, CO is simulated first with a mesoscale model and then, using the simulated mesoscale solution as initial and boundary conditions, within a microscale LES solver. The mesoscale solution generated by the Weather Research and Forecasting (WRF) model accurately captures the wind ramp event with only a slight temporal shift and a small temperature bias when compared to the meteorological towers in the NWTC site. This information is then taken as the initial and boundary conditions within the Simulator for Wind Farm Applications (SOWFA) microscale LES model with varying terrain inputs. These inputs include the high-resolution (roughly 30 m grid spacing) terrain data from the Shuttle Radar Topography Mission (SRTM) (Figure 1a), as well as, three intermediate terrain resolution data sets produced by applying a uniform convolution window to the SRTM data with varying window sizes (Figure 1b-d). These are compared to the WRF terrain data interpolated to the high-resolution microscale grid (Figure 1e) in order to determine the impact of low-resolution (about 1 km grid spacing) terrain data on turbulence generation. The wind field is compared at several heights and distances along the dominant wind flow axis between each of the terrain cases. It is found that the development of turbulence across the domain is accelerated as terrain information increases in fidelity. Further, because this ramp event at the location of interest is largely influenced by the interaction of the short-wave trough with the terrain, as the topographic information decreases in fidelity the accuracy of the ramp event is decreased. To make generalizations and describe the benefits and shortcomings of this approach, several of the existing perturbation techniques are employed and compared to that of turbulence generation through topographic features.

17 WIND ENERGY↗

Continuous conditional generative adversarial networks for data-driven solutions of poroelasticity with heterogeneous material properties

Machine learning-based data-driven modeling can allow computationally efficient time-dependent solutions of PDEs, such as those that describe subsurface multiphysical problems. In this work, our previous approach (Kadeethum et al., 2021d) of conditional generative adversarial networks (cGAN) developed for the solution of steady-state problems involving highly heterogeneous material properties is extended to time-dependent problems by adopting the concept of continuous cGAN (CcGAN). The CcGAN that can condition continuous variables is developed to incorporate the time domain through either element-wise addition or conditional batch normalization. Moreover, this framework can handle training data that contain different timestamps and then predict timestamps that do not exist in the training data. As a numerical example, the transient response of the coupled poroelastic process is studied in two different permeability fields: Zinn & Harvey transformation and a bimodal transformation. The proposed CcGAN uses heterogeneous permeability fields as input parameters while pressure and displacement fields over time are model output. Our results show that the model provides sufficient accuracy with computational speed-up. This robust framework will enable us to perform real-time reservoir management and robust uncertainty quantification in poroelastic problems.

97 MATHEMATICS AND COMPUTING↗

Development of Application Function Blocks for Power-Hardware-in-the-Loop Testing of Grid-Connected Inverters: Preprint

Stability and accuracy of the PHIL platform has been a common problem and challenging issue for a long time, thus continuous efforts and research have been devoted to developing interface algorithms to improve upon this. Some of existing works exhibit promising solutions to compensate the errors in the PHIL platform. However, there is a need to harmonize various solutions into a standard framework that provides a general solution for various PHIL test beds and that can be reused with minimum modifications to avoid extra effort. A good way to do this is by using application function blocks (AFBs), which are model-based and can be developed as a standard library component, thus achieving satisfactory reusability and reconfigurability for various platforms. Since PHIL simulation has the flexibility of the software simulation, proper functions can be easily implemented in the digital real-time simulator (DRTS) to preprocess the reference signal before sending it to the power amplifier. Therefore, this paper aims to develop standard AFBs in DRTS for the PHIL interface to test grid-connected inverters. The PHIL interface algorithm is separated into five AFBs that will be presented in detail in this paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Machine Learning Based Prediction of Airflow Maldistribution in A-Type Heat Exchangers

Airflow maldistribution is one of the primary causes of performance degradation in air-to refrigerant heat exchangers (HX) and has been shown to decrease heat transfer by as much as 35%. As a result, many units are oversized to meet the target capacity, resulting in increased system cost and refrigerant charge. Several studies have explored how characteristics like package type and HX geometry impact the flow profile, but results are restricted to a limited range of parameters and cannot be extrapolated to new designs. In this work, a machine learning (ML) model is trained to predict the inlet flow profile of dry air entering A-type HXs across a broad range of geometries and conditions. Flow profiles are generated using a porous media CFD model and used to train an Artificial Neural Network (ANN) which exhibits maximum and average relative L2 norm errors of 0.48 and 0.05. Additionally, these predictions take less than a second to generate resulting in a speed up factor of 2.42E5 compared to CFD. Component-level simulations are conducted to determine the performance degradation resulting from the predicted airflow maldistribution profiles. The new ML model will enable rapid and accurate prediction of performance degradation resulting from airflow maldistribution in A-type HXs, allowing for more accurate and cost-effective HX design.

42 ENGINEERING↗

Evaluating the Impact of Managed EV Charging for Reliable Operation of Bulk Power Systems with High Non-Dispatchable Generation

The growth of electric vehicles (EVs) and variable-generation (VG) sources introduces new challenges for power-system operations. This study introduces a modeling framework and evaluates five EV charging strategies under projected 2040 grid conditions in the Evergy service territory with high non-dispatchable generation. Using realistic EV behavior and generation models, their impacts on system peak demand, ramp rate, and reserve capacity are evaluated. Results show that only the peak-avoidance strategy effectively reduces system peak demand, while decentralized strategies-particularly cost based dynamic charging-can exacerbate peaks due to synchronized user behavior. However, ramp-rate minimization strategy significantly reduce the stress on dispatchable generation achieving the lowest maximum absolute ramp rate (MARR) (56.81 MW) and lowest reserve requirement (2.39 GW). In contrast, unmanaged and TOU random strategies increase the stress on dispatchable generation sources with increased MARR and reserve requirements. These findings highlight the importance of coordinated, system-aware managed charging strategies to ensure reliable and affordable grid operation in the presence of EVs and VG sources.

14 - SOLAR ENERGY↗

OmicsMLMentor: A Web Application for Guided Machine Learning Analysis of Omics Data

Expression-based omics technologies (e.g. proteomics, metabolomics, transcriptomics, etc.) increasingly rely on supervised and unsupervised machine learning (ML) models to find key biomolecules distinguishing conditions, identify natural groupings in biological data, or generate predictions for outcomes of interest. Fitting ML models to omics data presents several challenges, including handling missing data, selecting a normalization method, choosing a valid model, and optimizing hyperparameters, all requiring statistical programming skills to address these challenges. Thus, the open-source web application SLOPE was designed to lower the barrier to ML modeling for omics data. SLOPE supports the fitting of 15 ML models (10 supervised and 5 unsupervised) tailored to omics datasets, such as proteomics, metabolomics, lipidomics, and transcriptomics. SLOPE offers several omics-specific features, including methods for handling missingness (imputation, conversion, removal), normalization tests, ranking of models based on the structure of a user’s data and user input, and optimal hyperparameter selections using cross-validation splits. By streamlining ML workflows for omics analysis, SLOPE address critical gaps in existing online web tools, facilitating a broader adoption of these models for omics research. Here, SLOPE is applied to data from a lignin exposure study to highlight the workflow for fitting both supervised and unsupervised models to data.

lipidomics↗

Evaluation of historical precipitation interannual variability in CMIP6 over the United States

Interannual precipitation variability profoundly influences society via its effects on agriculture, water resources, infrastructure, and disaster risks. In this study, we use daily in situ precipitation observations from the global historical climatology network-daily (GHCN-D) to assess the ability of 21 Coupled Model Intercomparison Project Phase 6 (CMIP6) models, including the 50-member fifth-generation Canadian Earth System Model single model initial-condition large ensemble (CanESM5_SMILE), to realistically simulate historical interannual precipitation variability trends within 17 regions of the contiguous United States (CONUS). We assess how accurately the CMIP6 simulations align with observational data across annual, summer, and winter periods, focusing on four key hydrometeorological metrics, including interannual precipitation variability, relative interannual precipitation variability (coefficient of variation), annual mean precipitation, and annual wet day frequency. Our findings reveal that CMIP6 ensemble members generally reproduce the spatial patterns of observed trends in annual mean precipitation. In most regions, models agree well with the signs of observed changes in annual mean precipitation, though discrepancies in trend magnitude are evident. Further, observed trends in winter mean precipitation broadly exhibit a spatial pattern similar to that of the observed annual mean. However, analysis of the CanESM5_SMILE shows that trends in precipitation variability may primarily be the result of model-simulated internal variability, suggesting caution in interpreting multi-model single-realization ensemble results. Challenges in accurately simulating interannual precipitation variability underscore the need for ongoing model refinement and validation to enhance climate projections, especially in regions vulnerable to extreme precipitation events.

54 ENVIRONMENTAL SCIENCES↗

Optimization control technology for building energy conservation

A simulation processor generates and stores a simulation model based on conditions associated with a physical structure, such as a building. A neural network processor implements a neural network, having an input layer coupled to receive sensor data from the structure and having an output layer coupled to supply control signals to the at least one electrically operable environmental control device. The neural network is trained using the simulation model. A particle swarm optimization processor programmed to receive the simulation results and perform particle swarm optimization, ascertains optimal parameters for controlling the at least one electrically operable environmental control device and supplies these optimal parameters to the neural network processor. The neural network processor uses the optimal parameters supplied by the particle swarm optimization processor to further train the neural network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analysis of contact conditions and microstructure evolution in shear assisted processing and extrusion using smoothed particle hydrodynamics method

Shear assisted processing and extrusion (ShAPE) is a solid-phase processing technique that adds an additional shear force as compared with a conventional extrusion approach. Recently, ShAPE has demonstrated the capability of extruding high-performance aluminum alloy 7075 (AA7075) tubes at speeds up to 12.2 m/min without surface tearing. However, the relationship among the ShAPE processing parameters, thermomechanical conditions, contact conditions, heat generation, and microstructure evolution remains primarily empirical because an insightful understanding of the associated physics is still lacking. To help elucidate these relationships, this work proposes a thermomechanical meshfree model for the first time for ShAPE processing of AA7075 using the smoothed particle hydrodynamics (SPH) method. The meshfree model is first validated thoroughly by experimental data in terms of material flow, die face temperature, and extrusion force with various processing parameters. The validated model is then used to analyze the steady-state contact conditions and heat generation rates during ShAPE processing. Distributions of the average grain size of AA7075 being extruded are calculated using the SPH model output. The meshfree model results reveal that extrusions conducted at lower temperatures and higher strain rates yield more refined grains and possibly higher material strength, which is also consistent with the experimental observations.

36 MATERIALS SCIENCE↗

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↗