Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface roughness”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Crucial role of thermal fluctuations and vertex corrections for the magnetic pseudogap

It is generally believed that in a two-dimensional metal, whose ground state is antiferromagnetically ordered with Q = (π,π), thermal (static) magnetic fluctuations give rise to precursor behavior above TN in which the spectral function of a hot fermion (the one for which k and k + Q are Fermi surface points) contains two peaks, separated by roughly the same energy as in the antiferromagnetically ordered state. The two peaks persist in some range of T > T N and eventually merge into a single peak at zero frequency. This behavior is obtained theoretically by departing from free fermions in a paramagnet and evaluating the dressed fermionic Green's function by summing up infinite series of diagrams with contributions from thermal magnetic fluctuations. We show, following that keeping vertex renormalization diagrams in these series is crucial as other terms only broaden the spectral function of a hot fermion but do not shift its maximum away from zero frequency. As the consequence, the magnetic pseudogap should be treated as an input for theories that neglect vertex corrections, such as, e.g., Eliashberg theory for magnetically mediated superconductivity. We also analyze the potential pseudogap behavior at T = 0. Here, we argue that it may exist, but only at a finite correlation length, and not as a precursor to antiferromagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

LatticeAnalytics: Strut-Level Visualization and Inspection of Additively Manufactured Lattice Structures

Additive manufacturing (AM) is revolutionizing the production of custom components with complex internal geometries, essential for high-performance applications in diverse fields such as medicine and defense. These AM parts optimize strength while minimizing weight by utilizing internal lattice structures consisting of large quantities of small interconnected struts. However, the complexity of these structures, combined with the challenges of using X-ray Computed Tomography (XCT) data, makes validation of part reliability difficult. This ultimately inhibits the development of novel parts for our collaborating material scientists. Here, we introduce LatticeAnalytics, a novel framework specifically designed for visual inspection of defects in these lattice structures. Our framework offers an end-to-end solution that includes the data management of XCT scans, enables remote access for geographically dispersed teams through a web-based dashboard, and incorporates novel visualizations. Our analysis is facilitated by a coarse alignment between the lattice’s nominal model, a spatial graph, and the XCT data. We employ a simple VR-based approach for fast and rough alignment, followed by an offline registration and identification of the struts. With the nodes and struts aligned and identified in the volume, our framework allows querying of subvolumes containing a single strut at multiple resolutions. This avoids computation over the entire lattice and also allow for easy parallelization of down-stream computations, such as strut-specific metrics. To depict a fast overview of the strut quality, we introduce two innovative visual encodings, crucial for our collaborators’ research in creating novel AM parts: the Contour View and the Roughness Map, which depict critical geometrical and surface features of individual struts in standardized two 2D views. We evaluated the integrated system through expert interviews. The feedback confirms the framework’s practicality and its effectiveness in enhancing current inspection workflows. It solves major bottlenecks for our collaborators, ultimately helping them create novel parts with advanced properties.

Miao, Haichao [Lawrence Livermore National Laborat↗
DOE OSTI2025

Utah FORGE: Well 16B(78)-32 Drill Core Fracture Analysis Images and Data

This dataset contains drilling core data from well 16B(78)-32, including PDF documents with flattened core images annotated by feature type and core interval, as well as spreadsheets detailing feature morphologies by depth, planar feature measurements, and planar feature orientations rotated to in situ conditions. Core was recovered from three intervals, one per stimulation stage, in the crystalline rocks affected by the stimulation of well 16A(78)-32. Seven core runs were conducted, yielding 135.8 feet of recovered core. Features in the core were categorized into planar fractures, semi-planar fractures, unbroken mineralized fractures, rough fractures, curviplanar fractures, concave-convex surfaces, and planar compositional features such as mylonite or dike-like structures. Planar features were measured while the core was positioned horizontally, with the core axis aligned to a downhole azimuth of 42 degrees. Planar core measurements from stimulations 2 and 3 that could be confidently correlated with FMI data were rotated to in situ orientations. This was done by rotating the planes along vertical and horizontal axes to match the azimuth and inclination data recorded in the directional survey of well 16B(78)-32, as well as applying an axial rotation to resemble the fracture orientations observed in the FMI log at corresponding depths. Coherent sets of planar fracture measurements were made by aligning the core within each 3-foot section of the dissected core barrel, and between adjacent 3-foot sections within a core run by matching rock fabrics, saw cuts and/or tool marks. Where coherent fracture measurements could not be made within a core run, data sets are denoted by a subscript (i.e. 2-Ta and 2-Tb both come from tangent core run number 2).

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Study of unripe and inferior banana flours pre-gelatinized by four different physical methods

This study aimed to prepare the pre-gelatinized banana flours and compare the effects of four physical treatment methods (autoclaving, microwave, ultrasound, and heat-moisture) on the digestive and structural characteristics of unripe and inferior banana flours. After the four physical treatments, the resistant starch (RS) content values of unripe and inferior banana flours were decreased from 96.85% (RS2) to 28.99–48.37% (RS2 + RS3), while C∞ and k values were increased from 5.90% and 0.039 min −1 to 56.22–74.58% and 0.040–0.059 min −1 , respectively. The gelatinization enthalpy (ΔHg) and I 1047/1022 ratio (short-range ordered crystalline structures) were decreased from 15.19 J/g and 1.0139 to 12.01–13.72 J/g, 0.9275–0.9811, respectively. The relative crystallinity decreased from 36.25% to 21.69–26.30%, and the XRD patterns of ultrasound (UT) and heat-moisture (HMT) treatment flours maintained the C-type, but those samples pre-gelatinized by autoclave (AT) and microwave (MT) treatment were changed to C + V-type, and heat-moisture (HMT) treatment was changed to A-type. The surface of pre-gelatinized samples was rough, and MT and HMT showed large amorphous holes. The above changes in structure further confirmed the results of digestibility. According to the experimental results, UT was more suitable for processing unripe and inferior banana flours as UT had a higher RS content and thermal gelatinization temperatures, a lower degree and rate of hydrolysis, and a more crystalline structure. The study can provide a theoretical basis for developing and utilizing unripe and inferior banana flours.

Zhang, Siwei↗
DOE OSTI2025

Controlling Oxidation of Nb in Oxygen Abundant Environments

Modern particle accelerators depend on Superconducting Radio Frequency (SRF) cavities made from high-purity niobium (Nb) to achieve optimal performance, including high quality factors and strong accelerating gradients. However, when exposed to air, niobium naturally forms a complex oxide layer that can introduce surface imperfections and carbon contamination. This work examines an alternative oxidation strategy under oxygen-rich conditions to better regulate the oxide formation process. The ultimate objective is to improve surface uniformity and cleanliness, thereby reducing defect density and enhancing performance. We used Confocal Microscopy, Scanning Electron Microscopy (SEM), and X-Ray Photoelectron Spectroscopy (XPS) to analyze surface changes. We used standard metrics, like Arithmetic Average Roughness (Ra) and Root Mean Square (Rq). Three oxidation methods were applied: short, extended, and a repeated HF and H₂O₂ oxidation process. Preliminary results show promising reductions in carbon contamination and surface defects.

Romero, Juan [Fermilab]↗
DOE OSTI2020

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗
DOE OSTI2021

Alteration of pool boiling heat transfer on metallic surfaces by in situ oxidation

The critical heat flux during pool boiling has been investigated for a range of applications including electrical power generation and thermal management. Reported experimental CHF values during pool boiling of water on flat metallic surfaces, however, show a large discrepancy across studies. Here, we address this discrepancy in CHF values by accounting for oxidation of metallic surfaces during boiling. We studied the effect of in situ oxidation on flat Cu and Ni surfaces by changing the duration that samples were held in saturated water before conducting boiling experiments. The morphology and chemical composition of surfaces after the boiling experiments were analyzed by atomic force microscopy and X-ray photoelectron spectroscopy, respectively. Cu surfaces showed gradually increasing CHF values as the duration in saturated water increased, which could be attributed to the increase in roughness due to the formation of Cu 2 O nanostructures. Conversely, Ni surfaces showed relatively stable CHF and morphology as a nearly flat layer of NiO formed, with one exception: formation of a highly wetting hydroxide, Ni(OH) 2 , on a Ni coupon held in saturated water for 24 h resulted in a uniquely high CHF value, signifying the importance of surface chemistry in addition to morphology. Finally, the fundamental mechanisms resulting in the wide spread of CHF values on metallic surfaces elucidated in this work will lead to more accurate estimation of CHF as well as a deeper mechanistic understanding of CHF values on engineered surfaces.

42 ENGINEERING↗
DOE OSTI2023

A quantitative method to determine the region not influenced by injected interstitial and surface effects during void swelling in ion-irradiated metals

We propose and demonstrate a microstructurally-based experimental method to quantitively determine the depth regions in self-ion-irradiated metals that are affected by the injected interstitial effect and various surface effects, focusing on the choice of safe analysis zones to minimize the impact of these phenomena. The goal is to define the depth ranges where extracted data can be confidently applied to ion-neutron correlations for reactor application. Since ion energies in the range of 1–5 MeV are most frequently employed by the radiation effects community, irradiations were conducted at four energies in this range, all proceeding at 475 °C. The experiment was conducted on relatively pure single crystal iron to focus only on physical phenomena, avoiding the influence of possibly confounding chemical or segregation processes. Care was also taken to minimize the influence of other physical factors such as crystalline orientation. It was shown that, at 475 °C, ion energies of ≤1 MeV were too shallow in penetration and did not yield a safe depth range, but irradiations at 2.5 MeV and above yielded useful safe zones with predicted swelling behavior becoming independent of ion energy. The surface-affected zone width was found to be roughly twice that of the void-denuded zone width and to be independent of accumulated displacement dose. The largest injected-interstitial effect arises from the injected interstitial depression of void nucleation and growth. The interstitial-affected region starts at about one half of the projected range and does not show any “spreading” of its influence in depth as the peak damage level increases from 50 to 100 dpa. This study provides some confidence that enhances the credibility of ion simulation when applied to prediction of void swelling in neutron environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Temporal contrast degradation from mid-spatial-frequency surface error on stretcher mirrors

Temporal contrast degradation due to mid-spatial-frequency error in chirped-pulse amplification stretcher optics is studied. Third-order cross-correlation measurements reveal a temporal peak that appears when using two different mirrors processed by magnetorheological finishing, despite an improvement in rms roughness compared to a third unprocessed mirror. Simulations based on measured power spectral density show how the actual impact on contrast is different from measurements using a typical bandwidth-limited third-order cross correlator. Furthermore, strategies are proposed to avoid this type of contrast degradation while exploiting computer numerically controlled polishing techniques for enhancement of surface figure and roughness.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Oxide derived Cu nanofibril assembly for enhanced nonenzymatic glucose sensing

In this report, we have prepared copper nanofibril assembly via thermal oxidation followed by electrochemical reduction processes, exhibiting superior glucose detection ability. The morphological analysis evidenced the formation of rough islands or nanofibril structure on the Cu surface depending on initial thermal oxidation temperature. The glucose detection performance was investigated by performing cyclic voltammetry and chronoamperometry with varying glucose concentration. A high sensitivity (4131.57 μA mM – 1 cm – 2 ), low detection limit (1.41 μM), wider linear range (0–3.9 mM) and long–term stability (30 days) have been recorded for electrode thermally oxidized at 400 °C followed by electrochemical reduction process (rCu_400). The sensitivity is almost three times higher in comparison to the planner Cu surface. This significantly enhanced glucose sensing ability rCu_400 has been attributed to the nanofibril morphology, origination of Cu (111) facet and formation of a stable oxide layer evidenced by scanning electron microscopy, X-ray diffraction and energy dispersive spectroscopy analysis. Despite higher sensitivity, rCu_400 electrode does not show any response to the chloride ion, dopamine, ascorbic acid and uric acid. Finally, these results indicate that the Cu nanofibril structure prepared via simple oxidation/reduction process can be an excellent candidate to be used as an electrode for glucose sensing application.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A biokinetic model for systemic sodium

This paper describes an updated biokinetic model for systemic sodium (Na), developed for use in a series of reports by the International Commission on Radiological Protection (ICRP) on occupational intake of radionuclides. In contrast to the ICRP's previous model for intake of radio-sodium by workers, the updated model depicts realistic directions of movement of Na in the body including recycling of activity between blood and tissues. The updated model structure facilitates extension of the baseline transfer coefficients for adults to different age groups and to special exposure scenarios such as transfer of radio-sodium from the mother to the foetus or the nursing infant. Dose coefficients for 22Na and 24Na based on the updated model generally do not differ greatly from those based on the ICRP's previous Na model when both models are connected to the ICRP's latest dosimetry system. The main exception is that the updated model yields roughly twofold higher dose coefficients for endosteal bone surface than does the previous model due to the dosimetrically cautious assumption in the updated model that exchangeable Na in bone resides on bone surface.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

3D Printed Parahydrophobic Surfaces as Multireaction Platforms

We report Parahydrophobic surfaces (PHSs) composed of arrays of cubic μ-pillars with a double scale of roughness and variable wettability were systematically obtained in one step and a widely accessible stereolithographic Formlabs 3D printer. The wettability control was achieved by combining the geometrical parameters (H = height and P = pitch) and the surface modification with fluoroalkyl silane compounds. Homogeneous distribution of F and Si atoms onto the pillars was observed by XPS and SEM-EDAX. A nano-roughness on the heads of the pillars was achieved without any post-treatment. The smallest P values lead to surfaces with static contact angles (CAs) >150° regardless of the H utilized. Interestingly, the relationship 0.6 ≤ H/P ≤ 2.6 obtained here was in good agreement with the H/P values reported for nano- and submicron pillars. Furthermore, experimental CAs, advancing and receding CAs, were consistent with the theoretical prediction from the Cassie–Baxter model. Structures covered with perfluorodecyltriethoxysilane with high H and short P lead to PHSs. Conversely, structures covered with perfluorodecyltrimethoxysilane exhibited a superhydrophobic behavior. Finally, several aqueous reactions, such as precipitation, coordination complex, and nanoparticle synthesis, were carried out by placing the reactive agents as microdroplets on the parahydrophobic pillars, demonstrating the potential application as chemical multi-reaction array platforms for a large variety of relevant fields in microdroplet manipulation, microfluidics systems, and health monitoring, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗
DOE OSTI2026

Scattering Processes from Quantum Simulation Algorithms for Scalar Field Theories

We provide practical simulation methods for scalar field theories on a quantum computer that yield improved asymptotics as well as concrete gate estimates for the simulation and physical qubit estimates using the surface code. We achieve these improvements through two optimizations. First, we consider a finite volume approach for estimating the elements of the S-matrix. This approach is appropriate in general for 1+1D and for certain low-energy elastic collisions in higher dimensions. Second, we implement our approach using a series of different fault-tolerant simulation algorithms for Hamiltonians formulated both in the field occupation basis and field amplitude basis. Our algorithms are based on either second-order Trotterization or qubitization. The cost of Trotterization in occupation basis scales as O ( λ N 7 | Ω | 3 / ( M 5 / 2 ϵ 3 / 2 ) ) where λ is the coupling strength, N is the occupation cutoff, | Ω | is the volume of the spatial lattice, M is the mass of the particles and ϵ is the uncertainty in the energy calculation used for the S -matrix determination. Qubitization in the field basis scales as O ( | Ω | 2 ( k 2 Λ + k M 2 ) / ϵ ) , where k is the cutoff in the field and Λ is a scaled coupling constant. We find in both cases that the bounds suggest physically meaningful simulations can be performed using on the order of 4 × 10 6 physical qubits and 10 12 T -gates which corresponds to roughly one day on a superconducting quantum computer with surface code and a cycle time of 100 ns. This places the simulation of scalar field theory within striking distance of the gate counts for the best available chemistry simulation results.

Hardy, Andrew [Toronto U.] (ORCID:0000000235817382↗
DOE OSTI2022

Experimental study of crystallization fouling by calcium carbonate: Effects of surface structure and material

We report mitigation of crystallization fouling on heat transfer surface is crucial for developing highly efficient desalination techniques, because the tenacious and porous crystal layer can drastically deteriorate system efficiency. Surface modification is a potential approach to decrease scale formation and promote crystal removal, but some of the mechanisms between crystallization and substrate properties are still not fully understood. In this study, we use dual-angle microscopic visualization and image processing to quantify three-dimensional crystallization behavior and analyze the impact of surface structure and material. The transition of crystal size and distance with roughness and wettability are identified from statistical analysis. The determination of crystal porosity and effective thermal conductivity is enabled by visualized geometry and measured deposition mass. This study shed lights on potential strategies of surface modification to mitigate scaling in desalination by enhancing scale removal and localizing crystal growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Influence of preheating on tungsten damage under pulsed hydrogen plasma exposure

Preheated tungsten plate behavior under multiple pulsed hydrogen plasma irradiations in the PW-7 plasma accelerator was investigated. The preheating temperature was 973 K, which is above the ductile to brittle transition temperature (DBTT ≈ 673 K). The plate was exposed to three consecutive series of 13 plasma pulses each. Preheating sufficiently suppressed surface cracking at the early irradiation stage: after 13 pulses, cracks were rare, localized, and did not form extensive fractures. During subsequent irradiation series, the area fraction of fine cracks increased by ∼ 2% due to accumulation of dislocations around the pre-existing defects. However, their width remained nearly constant, averaging 0.2–0.3 μm, with rare cracks up to 1 μm in size. Repeated plasma exposure resulted in pronounced surface roughening associated with multiple melting cycles, crater formation and blister rupture, with roughness reaching ∼ 204 nm and a recorded sample mass loss of ∼ 2.8 mg. These results indicate that controlling tungsten preheating well above DBTT effectively reduces damage. However, the surface condition during initial irradiation significantly influences its subsequent morphological evolution. Preheating to 773 K, selected regarding the minimum DBTT value, has previously been demonstrated to be an effective approach for limiting crack growth. A further increase in temperature of approximately ∼ 1000 K did not qualitatively improve results: crack still formed and widths did not decrease. At the same time, the absence of further crack coarsening indicates that the limiting effect of preheating has been reached. Considering the increased energy costs, further increasing the preheating temperature seems unjustified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics↗
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↗