Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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

Experimental study on kinetic oxidation of graphite IG-110 by steam

Graphite is proposed for use in High-temperature Gas-cooled Reactors (HTGRs) as the fuel matrix, neutron moderator/reflector, and core structural material. One important property of nuclear grade graphite is their resistance to oxidation in high-temperature environment. Extensive investigation has been performed in the literature for graphite oxidation by air. However, available experimental data are still limited for graphite oxidation by steam under conditions comparable to a postulated steam ingress accident in HTGRs. In this study, the oxidation rate of graphite IG-110 by steam was measured at temperatures from 850 to 1100 °C with the steam partial pressure varying from 0.5 to 20.0 kPa and the hydrogen partial pressure varying from 0 to 2.0 kPa. Further analysis confirms the oxidation process in this present study is dominated by the chemical kinetics, which lends credit to the data for being used to develop numerical models. It was observed that the increase of the kinetic oxidation rate with the steam partial pressure tends to become less apparent if the steam partial pressure keeps increasing. In addition, it was found that the partitioning of hydrogen inhibits the graphite-steam reaction process even with the steam partial pressure up to 20.0 kPa. However, this inhibiting effect starts to become saturated when the hydrogen partial pressure exceeds 1.0 kPa. The oxidation rates were fitted to the conventional Langmuir-Hinshelwood (LH) and Boltzmann-enhanced Langmuir-Hinshelwood (BLH) models by a multivariable optimization algorithm. The BLH model exhibits a better accuracy than the LH model within the specified experimental conditions. The predicted oxidation rate using the BLH model shows a mean relative difference of about 24% with the maximum difference of about 55% when compared with our experimental data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

In the last decades, the Arctic climate has changed dramatically, with the loss of multiyear sea ice one of the clearest consequences. These changes have occurred on relatively rapid timescales, and both accurate short-term Arctic prediction (e.g., 10 days to three months) and climate projection of future Arctic scenarios present ongoing challenges. Here we describe a representation of the Arctic ocean and sea ice in a ultrahigh resolution simulation in which the horizontal grid mesh reduces from 8 km at the equator to 2 km at the poles (UH8to2) for the years 2017-2020. We find the simulation reproduces observed distributions of seasonal sea-ice thickness and concentration realistically, although concentration is biased low in the spring and summer and low biases in thickness are found in the central and eastern basins in the fall. Volume, fresh water, and heat transports through key passages are realistic, lying within observationally determined ranges. Climatological comparisons reveal that the UH8to2 Atlantic Water is shallower, warmer, and saltier than the World Ocean Atlas 2018 climatology for 2005-2017 in the eastern basin. Our analysis suggests that these biases, combined with a lack of stratification in the upper 100 m of the simulated ocean, contribute to the winter biases in modeled sea ice thickness. This relationship between biases in the sea ice and ocean points to a potential positive feedback within the model, illuminating challenges for long term model predictive power in a changing Arctic climate.

58 GEOSCIENCES↗
DOE OSTI2022

Deep-learning-based workflow for boundary and small target segmentation in digital rock images using UNet++ and IK-EBM

We report three-dimensional (3D) X-ray micro-computed tomography (μCT) has been widely used in petroleum engineering because it can provide detailed pore structural information for a reservoir rock, which can be imported into a pore-scale numerical model to simulate the transport and distribution of multiple fluids in the pore space. The partial volume blurring (PVB) problem is a major challenge in segmenting raw μCT images of rock samples, which impacts boundaries and small targets near the resolution limit. We developed a deep-learning (DL)-based workflow for accurate and fast partial volume segmentation. The DL model's performance depends primarily on the training data quality and model architecture. This study employed the entropy-based-masking indicator kriging (IK-EBM) to segment 3D Berea sandstone images as training datasets. The comparison between IK-EBM and manual segmentation using a 3D synthetic sphere pack, which had a known ground truth, showed that IK-EBM had higher accuracy on partial volume segmentation. We then trained and tested the UNet++ model, a state-of-the-art supervised encoder-decoder model, for binary (i.e., void and solid) and four-class segmentation. We compared the UNet++ with the commonly used U-Net and wide U-Net models and showed that the UNet++ had the best performance in terms of pixel-wise and physics-based evaluation metrics. Specifically, boundary-scaled accuracy demonstrated that the UNet++ architecture outperformed the regular U-Net architecture in the segmentation of pixels near boundaries and small targets, which were subjected to the PVB effect. Feature map visualization illustrated that the UNet++ bridged the semantic gaps between the feature maps extracted at different depths of the network, thereby enabling faster convergence and more accurate extraction of fine-scale features. The developed workflow significantly enhances the performance of supervised encoder-decoder models in partial volume segmentation, which has extensive applications in fundamental studies of subsurface energy, water, and environmental systems.

02 PETROLEUM↗
DOE OSTI2020

A crisis for the verification and validation of turbulence simulations

Three algorithms have been proposed for solution of the Rayleigh–Taylor turbulent mixing problem. They are based upon three different physical principles governing the Euler equations for fluid flow. The principles serve to select the physically relevant solution from among many nonunique solutions. The admissibility principle is in dispute. The three different algorithms, expressing the three physical admissibility principles can be formulated in terms of the three energy dissipation rates or the entropy production rates, as selected by the size of the sub grid scale coefficients. These have maximal values or less than maximal values. The resulting solutions are markedly different. We find strong validation evidence that supports the maximum rate principle, based on a review of prior results and on new results presented here. We review experimental data used for validation and sufficient to discriminate among the three. We present a new analysis of this data. We show that the hypothesized long wave length perturbations in the initial conditions are not significant, so that validation can be based on this data in a straight forward manner. One of the algorithms is labeled direct numerical simulation, but is not, and as a consequence, the two algorithms with less than maximal SGS coefficients are variants of one another. Finally recommendations for the numerical modeling of the deflagration to detonation transition in type Ia supernova are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Passive heat removal in horizontally oriented micro-HTGRs

There are novel, horizontally oriented high-temperature gas-cooled micro-reactor (HTGR) designs under consideration because of their portability, simplicity, and reliable operation. These features of the micro-HTGR design can meet the fission battery attributes, such as economic, standardize, installed, unattended and reliable. Although there have been HTGR related safety studies, micro-HTGRs when horizontally oriented are expected to exhibit significantly different thermal physics. During the unavailability of helium circulation, the internal reactor core is designed to cool by block-to-block conduction and radiation, and the reactor vessel surface is cooled by the ambient air. This scenario is anticipated during the transport of the micro-HTGR in a shipping container. The conduction and radiation between the prismatic micro-HTGR blocks in the core can be influenced by variances in the thermal contacts. This work investigated the conduction within a simulated horizontal HTGR core. An experimental setup was used to validate a numerical model. The experimental setup consisted of a hexagonal assembly with scaled prismatic blocks placed within a high-temperature vacuum environment. The gaps between the blocks were well controlled and monitored. The experimental setup was designed in such a way that the temperature variation in the axial direction was minimal, such that the experiment could be observed as a 2D (r, ) heat transfer problem. The pressurized tube was equipped with Infra-Red (IR) transparent windows, and an IR camera was used to capture the spatio-temporal temperature evolution in the hexagonal block assembly during the cooldown experiments. The experimental scenario was computationally modeled with a finite element analysis (FEA) program. Once validated, the computational model was used to investigate the impact of gap conductance on overall decay heat removal. Using a conservative estimate for gap conductance value (100 W/m 2 -K) between the prismatic blocks, there is a negligible increase in temperature observed during decay heat generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Post-detonation fireball thermometry via femtosecond-picosecond coherent anti-Stokes Raman Scattering (CARS)

Accurate knowledge of post-detonation fireball temperatures is important for understanding device performance and for validation of numerical models. Such measurements are difficult to make even under controlled laboratory conditions. Here, temperature measurements were performed in the fireball of a commercial detonator (RP-80, Teledyne RISI). The explosion and fragments were contained in a plastic enclosure with glass windows for optical access. A hybrid femtosecond-picosecond (fs-ps) rotational coherent anti-Stokes Raman scattering (CARS) instrument was used to perform gas-phase thermometry along a one-dimensional measurement volume in a single laser shot. The 13-mm-thick windows on the explosive-containment housing introduced significant nonlinear chirp on the fs lasers pulses, which reduced the Raman excitation bandwidth and did not allow for efficient excitation of high-J Raman transitions populated at flame temperatures. To overcome this, distinct pump and Stokes pulses were used in conjunction with spectral focusing, achieved by varying the relative timing between the pump and Stokes pulses to preferentially excite Raman transitions relevant to flame thermometry. Light scattering from particulate matter and solid fragments was a significant challenge and was mitigated using a new polarization scheme to isolate the CARS signal. Fireball temperatures were measured 35–40 mm above the detonator, 12–25 mm radially outward from the detonator centerline, and at 18 and 28 µs after initiation. At these locations and times, significant mixing between the detonation products and ambient air had occurred thus increasing the nitrogen-based CARS thermometry signal. Initial measurements show a distribution of fireball temperatures in the range 300–2000 K with higher temperatures occurring 28 µs after detonation.

42 ENGINEERING↗
DOE OSTI2026

Complex exposure-burial history and Pleistocene sediment recycling in the dead sea rift with implications for the age of the Acheulean site of ‘Ubeidiya

Here, we present early Pleistocene burial ages of the ‘Ubeidiya Formation sediments in the Jordan Valley, a segment of the Dead Sea Rift Valley. A minimum age of ∼1.1 Ma is constrained both by paleomagnetic analysis and U-Pb dating of Melanopsis shells. Simple cosmogenic burial ages (i.e., one very long exposure period followed by one period of burial) calculated from the ratios of 26 Al to 10 Be and 10 Be concentrations indicate ages of ∼3 Ma, contradicting the geological and paleomagnetic constraints as well as a reasonable age of the ‘Ubeidiya archeological site, as it contains human remains. A more sophisticated way of treating the results, by combining numerical modeling of cosmogenic nuclide build up during repeated burial-exposure cycles, paleomagnetic analysis (indicating reverse polarity) and a minimum burial age of 1.1 Ma set by U-Pb dating of Melanopsis shells, suggests two most probable time slots (1.19-1.77 and 1.93- 2.14 Ma) for the absolute age of the ‘Ubeidiya Formation. However, the skewed distribution of cosmogenic isotope burial ages, with a median age of 2.05 Ma, indicates a much higher probability of the older ages, most likely >1.9 Ma. The exposure-burial history that emerges from the model implies recycling of sediments previously deposited and buried in the rift valley, between 4.5 and 3.2 Ma, and subsequently exhumed, eroded, transported and then redeposited along the ‘Ubeidiya paleo lake shoreline.

Acheulean↗
DOE OSTI2020

Performance and wake flow characterization of a 1:8.7-scale reference USDOE MHKF1 hydrokinetic turbine to establish a verification and validation test database

As hydrokinetic turbine technologies continue to advance towards commercialization, public datasets on the performance characteristics for these devices and their flow field effects are invaluable to advance our understanding of these technologies and to validate analytical and numerical models. The Applied Research Laboratory at The Pennsylvania State University (ARL Penn State) collaborated with Sandia National Laboratories and the University of California at Davis to design, fabricate (at a 1:8.7 scale), and experimentally test a novel hydrokinetic turbine rotor design to provide an open platform and dataset for further study and development. The water tunnel test of this three-bladed, horizontal-axis rotor recorded power production, blade loading, the near-wake flow, cavitation effects, and noise generation. These state-of-the-art measurements demonstrate much of the complex physics associated with the flow through an unducted, horizontal-axis turbine, and they elucidate the performance characteristics and flow field effects at an unprecedented fidelity, accuracy and resolution. Measurements of powering coefficients (power, torque and thrust) as a function of tip-speed-ratio were performed. The dataset also includes unsteady measurements of driveshaft loading, blade strain, tower pressures, and radiated noise. Detailed flow mapping using laser Doppler velocimetry, and planar and stereo particle image velocimetry includes measurements of mean velocity and Reynolds stresses. Although the wake measurements are limited to less than half a diameter, they reveal the complex flow patterns in the near-wake structure of the rotor. The full database, available at the United States Department of Energy’s marine and hydrokinetic data repository, includes tunnel and model Computer Aided Design geometry files and inflow data sufficient for a “Model-the-Test” computational Verification and Validation study.

13 HYDRO ENERGY↗
DOE OSTI2023

A comparison of the power potential for surface- and seabed-deployed tidal turbines in the San Juan Archipelago, Salish Sea, WA

The San Juan archipelago lies along the axis of tidal movement between Straits of Juan de Fuca and Strait of Georgia in the Salish Sea. The amplitude of the tidal exchange produces significant tidal currents between the islands, as well as in Rosario and Haro Straits. These are of interest as a future source of electrical power generation, given the islands’ dependence on electricity supply by a subsea cable from the mainland. Here, we evaluate the tidal current energy potential in this region through a re-analysis of measurements collected by the National Ocean Service (NOS) and a high-resolution numerical model. Given the considerable variations in water depth and vertical velocity profiles across candidate tidal energy sites, we consider the trade-offs between tidal turbines deployed from a floating platform and those anchored to the seabed. Measurement re-analysis indicates several locations that could support tidal current power generation by MW-scale turbines with an acceptable balance between turbine size, rated power, and capacity factor. Even for relatively large (30 m) turbine diameters, surface-deployed turbines would be expected to produce up to 30 % more electricity than the same turbine deployed near the seabed due to vertical shear, with this difference increasing for smaller diameter turbines. A benchmark comparison at the measurement locations with the highest power generation potential shows relatively good model fidelity, though, even with relatively high resolution, time-average power density disagreements of ± 50 % persist throughout the water column. Nonetheless, the simulation identifies several locations with power generation potential more than twice as high as locations in the measurement re-analysis. These sites were either not surveyed by NOS or excluded due to data quality issues. Altogether, these results demonstrate the complimentary of models and measurements for tidal energy site assessment and suggest that tidal currents could be an important source of electricity generation in the San Juan archipelago.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Spatially varying seasonal modulation to tidal stream energy potential due to mixed tidal regimes in the Aleutian Islands, AK

We provide an assessment of the tidal stream energy resource of the Aleutian Islands, Alaska via a validated barotropic tidal numerical model of the region. Eight island passes are identified as energy “hotspots”. The annual mean kinetic energy fluxes, KEF , calculated at each pass vary from 1000 to 11,000 MW, while the annual available energy, AAE , varies from 5 to 42 MWh m −2 . Notable seasonal modulation to monthly power density averages and ranges are noted at some passes and not others. Seasonal adjustment is linked to the semi-annual solar declination cycle which enhances (dampens) diurnal (D 1 ) tidal amplitudes in summer/winter (spring/fall) as well as the time-varying phase lag between D 1 and semidiurnal (D 2 ) fortnightly tidal cycles. Annual variability in monthly mean power density scales with the tidal current form factor, F u , with the largest seasonal change occurring for F u > 1 (D 1 dominated tide). The spread in power density over a month is on average smaller for passes with mixed tides (F u = 1 ) than those with D 1 or D 2 dominance, as changes to fortnightly phase lag become influential to net power density ranges when tides are mixed. This study outlines overlooked, but relevant, long-term modulation to tidal streams in regions with mixed tides.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Underground sources of drinking water chemistry changes in response to potential CO 2 leakage

The purpose of this study was to quantify changes to underground sources of drinking water (USDW) quality in response to potential CO 2 leakage from geologic CO 2 sequestration (GCS) reservoirs. Here, the developed a framework of combined laboratory experiments and reactive transport simulations and used this framework to evaluate the Ogallala aquifer overlying the Farnsworth Unit (FWU), an active GCS site, as a case study. Using chemical reaction parameters obtained from laboratory experiments and numerical simulations, site-specific mechanisms of CO 2 -water-sediment interactions at the USDW aquifer were interpreted. Long-term risks of potential CO 2 leakage were then evaluated with field-scale numerical models using the regional hydrogeological characteristics and reaction parameters obtained from our experiments and simulations. Results suggest that carbonate mineral impurity and cation exchange are key mechanisms for interactions between CO 2 and the aquifer sediment. Additionally, for a large leakage rate of 0.1 % injection from one leaky well, the leakage plume might impact an area of 300 m in diameter and significantly affect the local water quality by changing pH and cation concentrations (e.g., Zn, Ba and Sr). After leakage ceases, the zone of impacted fluids would not migrate significantly in subsequent decades due to a low regional groundwater flowrate (for this case study). The relatively small area of impact might not be detected in a monitoring well given the broader spacing in a typical field scenario. Effective early leakage detection may require additional tools, e.g., borehole CO 2 movement, four-dimensional seismicity, CO 2 soil flux, samples from deeper aquifers, etc., to ensure effective leakage detection and long-term safety of GCS projects.

58 GEOSCIENCES↗
DOE OSTI2023

Identifying snowfall elevation patterns by assimilating satellite-based snow depth retrievals

Precipitation in mountain regions is highly variable and poorly measured, posing important challenges to water resource management. Traditional methods to estimate precipitation include in-situ gauges, Doppler weather radars, satellite radars and radiometers, numerical modeling and reanalysis products. Each of these methods is unable to adequately capture complex orographic precipitation. Here, we propose a novel approach to characterize orographic snowfall over mountain regions. We use a particle batch smoother to leverage satellite information from Sentinel-1 derived snow depth retrievals and to correct various gridded precipitation products. This novel approach is tested using a simple snow model for an alpine basin located in Trentino Alto Adige, Italy. Here, we quantify the precipitation biases across the basin and found that the assimilation method (i) corrects for snowfall biases and uncertainties, (ii) leads to cumulative snowfall elevation patterns that are consistent across precipitation products, and (iii) results in overall improved basin-wide snow variables (snow depth and snow cover area) and basin streamflow estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Numerical analyses of high temperature dense, granular flows coupled to high temperature flow property measurements for solar thermal energy storage

High temperature particle flow properties necessary to predict granular flow behavior for solar thermal energy storage applications were measured and calculated for Carbobead CP 30/60 up to 800 °C. Here, the measured properties included elastic and shear moduli, particle-particle coefficients of static sliding and rolling friction, and particle-particle coefficients of restitution. Poisson’s ratio was calculated with elastic and shear moduli. The flow properties were used as inputs for a numerical model using the discrete element method to examine granular flows along an inclined plane at high temperature. The flow behavior was strongly influenced by the coefficients of static friction, which impacted the particle residence time, shear effects from the side walls, and particle flow mass flux. An 8.7%, 15.6%, and 8.5% increase and 37.9% decrease in steady state mass flow rate was observed for 200 °C, 400 °C, 600 °C, and 800 °C, respectively, when compared to room temperature simulations. A 52%, 59%, and 33% decrease in the time to reach steady state was observed for 200 °C, 400 °C, and 600 °C, respectively, while a 53% increase in time was observed for 800 °C. A significant delay in the flow development at 800 °C was observed due to significantly higher frictional forces.

14 SOLAR ENERGY↗
DOE OSTI2022

Significant improvement of conversion efficiency by passivation of low-angle grain boundaries in flexible low-cost single-crystal-like GaAs thin-film solar cells directly deposited on metal tape

We develop an easily implementable yet effective process to achieve high power conversion efficiency in newly-developed flexible single-crystal-like GaAs thin-film photovoltaic solar cells on cheap metal tape. We study various techniques on the fabricated devices to passivate the defect states associated with dangling bonds at low-angle grain boundaries (LA-GBs) existing in the single-crystal-like GaAs material. We report that sulfur by trioctylphosphine sulfide (TOP:S) solution treatment can partially but effectively passivate the defect states at the LA-GBs to increase all the solar cell performance parameters of open-circuit voltage VOC (33.2%), short-circuit current JSC (19.5%), and fill factor FF (25.8%), resulting in the improvement of the power conversion efficiency more than double. The material and device characterizations and numerical modeling are performed to propose a mechanism, also suggesting that the power conversion efficiency higher than 20%, which is an important target milestone in emerging solar cell technologies, can be achieved. This work presents the first effective passivation of the bulk defects in GaAs at the device level and can also offer a solution for other photonic and electronic devices using non-single-crystalline GaAs films for flexible electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Substrate matters: Coupled phonon modes of a spherical particle on a substrate probed with EELS

Using vibrational electron energy loss spectroscopy (vib-EELS) combined with numerical modeling, we investigate the physical mechanisms governing the phonon coupling between a spherical particle sustaining multipolar surface phonon modes and an underlying thin film. Depending upon their dielectric composition, a variety of hybrid phonon modes arise in the EEL spectrum due to the interaction between polarization charges in the particle and film. Mirror charge effects and phonon mode hybridization are the active mechanisms acting on dielectric and metallic-type films, respectively. Furthermore, processes beyond dipole-dipole interactions are required to describe the sphere-film coupling.

Localized surface phonon response↗
DOE OSTI2021

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effect of ionization waves on dust chain formation in a DC discharge

An interesting aspect of complex plasma is its ability to self-organize into a variety of structural configurations and undergo transitions between these states. A striking phenomenon is the isotropic-to-string transition observed in electrorheological complex plasma under the influence of a symmetric ion wake field. Such transitions have been investigated using the Plasma Kristall-4 (PK-4) microgravity laboratory on the International Space Station. Recent experiments and numerical simulations have shown that, under PK-4-relevant discharge conditions, the seemingly homogeneous direct current discharge column is highly inhomogeneous, with large axial electric field oscillations associated with ionization waves occurring on microsecond time scales. here, a multi-scale numerical model of the dust–plasma interactions is employed to investigate the role of the electric field in the charge of individual dust grains, the ion wake field and the order of string-like structures. Results are compared with those for dust strings formed in similar conditions in the PK-4 experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗