Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled simulation”

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

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

At least 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Large eddy simulation of wind farm performance in horizontally and vertically staggered layouts

This numerical investigation employs Large Eddy Simulation (LES) coupled with Actuator Disk Model (ADM) to evaluate wind farm layout optimization strategies. The study presents a systematic analysis of aligned, horizontal staggering, vertical staggering, and mixed (combination of horizontal and vertical) staggering configurations, aiming to establish optimal design parameters for enhanced power production. The investigation examines key performance metrics including mean velocity distributions, turbulence intensity characteristics, and power generation efficiency. Results demonstrate better performance of both horizontal and vertical staggering patterns compared to conventional aligned configurations, with horizontal staggering exhibiting notably higher power output than vertical arrangements. Our findings also suggest that mixed configurations, incorporating both horizontal and vertical staggering, can offer optimal performance characteristics. As a result, this research advances the understanding of wake interactions in complex wind farm layouts and provides design guidelines for maximizing wind farm power generation efficiency through strategic turbine positioning.

17 WIND ENERGY↗
DOE OSTI2015

On the impact of temperatures up to 200°C in clay repositories with bentonite engineer barrier systems: A study with coupled thermal, hydrological, chemical, and mechanical modeling

One of the most important design variables for a geological nuclear waste repository is the temperature limit up to which the engineered barrier system (EBS) and the natural geologic environment can be exposed. Up to now, almost all design concepts that involve bentonite-backfilled emplacement tunnels have chosen a maximum allowable temperature of about 100°C. Such a choice is largely based on the consideration that in clay-based materials illitization and the associated mechanical changes in the bentonite (and perhaps the clay host rock) could affect the barrier attributes of the EBS. However, existing experimental and modeling studies on the occurrence of illitization and related performance impacts are not conclusive, in part because the relevant couplings between the thermal, hydrological, chemical, and mechanical (THMC) processes have not been fully represented in the models. This work presents a fully coupled THMC simulation of a nuclear waste repository in a clay formation with a bentonite-backfilled EBS for 1000 years. Two scenarios were simulated for comparison: a case in which the temperature in the bentonite near the waste canister can reach about 200°C and a case in which the temperature in the bentonite near the waste canister peaks at about 100°C. The model simulations demonstrate some degree of illitization in both the bentonite buffer and the surrounding clay formation. Other chemical alterations include the dissolution of K-feldspar and calcite, and precipitation of quartz, chlorite, and kaolinite. In general, illitization in the bentonite and the clay formation is enhanced at higher temperature. However, the quantity of illitization is affected by many chemical factors and therefore varies a great deal. The most important chemical factors are the concentration of K in the pore water as well as the abundance and dissolution rate of K-feldspar; less important ones are the concentration of sodium and the quartz precipitation rate. In our modeling scenarios, the calculated decrease in smectite volume fraction in bentonite ranges from 1 to 8% of the initial volume fraction of smectite in the 100°C scenario and 1–27% in the 200°C scenario. Chemical changes in the 200°C scenario could also lead to a reduction in swelling stress up to 15–18% whereas those in the 100°C scenario result in about 14–15% reduction in swelling stress for the base case scenario. Model results also show that the 200°C scenario results in a much higher total stress than the 100°C scenario, mostly due to thermal pressurization. While cautions should be taken regarding the model results due to some limitations in the models, the modeling work is illustrative in light of the relative importance of different processes occurring in EBS bentonite and clay formation at higher than 100°C conditions, and could be of greater use when site specific data are available.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗
DOE OSTI2016

A hybrid Boundary Element Unstructured Transmission-line (BEUT) method for accurate 2D electromagnetic simulation

Time domain electromagnetic simulation tools have the ability to model transient, wide-band applications, and non-linear problems. The Boundary Element Method (BEM) and the Transmission Line Modeling (TLM) method are both well established numerical techniques for simulating time-varying electromagnetic fields. The former surface based method can accurately describe outwardly radiating fields from piecewise uniform objects and efficiently deals with large domains filled with homogeneous media. The latter volume based method can describe inhomogeneous and non-linear media and has been proven to be unconditionally stable. Furthermore, the Unstructured TLM (UTLM) enables modelling of geometrically complex objects by using triangular meshes which removes staircasing and unnecessary extensions of the simulation domain. The hybridization of BEM and UTLM which is described in this paper is named the Boundary Element Unstructured Transmission-line (BEUT) method. It incorporates the advantages of both methods. The theory and derivation of the 2D BEUT method is described in this paper, along with any relevant implementation details. The method is corroborated by studying its correctness and efficiency compared to the traditional UTLM method when applied to complex problems such as the transmission through a system of Luneburg lenses and the modelling of antenna radomes for use in wireless communications. - Graphical abstract:.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A hybridizable discontinuous Galerkin method for modeling fluid–structure interaction

This study presents a novel application of the hybridizable discontinuous Galerkin (HDG) finite element method to the multi-physics simulation of coupled fluid–structure interaction (FSI) problems. Recent applications of the HDG method have primarily been for single-physics problems including both solids and fluids, which are necessary building blocks for FSI modeling. Utilizing these established models, HDG formulations for linear elastostatics, a nonlinear elastodynamic model, and arbitrary Lagrangian–Eulerian Navier–Stokes are derived. The elasticity formulations are written in a Lagrangian reference frame, with the nonlinear formulation restricted to hyperelastic materials. With these individual solid and fluid formulations, the remaining challenge in FSI modeling is coupling together their disparate mathematics on the fluid–solid interface. This coupling is presented, along with the resultant HDG FSI formulation. Verification of the component models, through the method of manufactured solutions, is performed and each model is shown to converge at the expected rate. The individual components, along with the complete FSI model, are then compared to the benchmark problems proposed by Turek and Hron [1]. The solutions from the HDG formulation presented in this work trend towards the benchmark as the spatial polynomial order and the temporal order of integration are increased.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Structure-induced microalloying effect in multicomponent alloys

In this study, the microalloying effect on glass-forming ability (GFA) has been investigated from the structural aspect, by performing synchrotron radiation x-ray diffraction and absorption measurements coupled with simulations in the NiNbZr ternary system. We propose a new parameter which counts the fraction of the fivefold symmetries in all clusters and find it is strongly associated with the GFA. In particular, this structural parameter has the highest value in a composition where the best GFA is achieved. The present work provides an in-depth understanding of microalloying-induced high GFAs in multicomponent alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Electron Accelerator Shielding Design of KIPT Neutron Source Facility

The Argonne National Laboratory of the United States and the Kharkov Institute of Physics and Technology of the Ukraine have been collaborating on the design, development and construction of a neutron source facility at Kharkov Institute of Physics and Technology utilizing an electron-accelerator-driven subcritical assembly. The electron beam power is 100 kW using 100-MeV electrons. The facility was designed to perform basic and applied nuclear research, produce medical isotopes, and train nuclear specialists. The biological shield of the accelerator building was designed to reduce the biological dose to less than 5.0e-03 mSv/h during operation. The main source of the biological dose for the accelerator building is the photons and neutrons generated from different interactions of leaked electrons from the electron gun and the accelerator sections with the surrounding components and materials. The Monte Carlo N-particle extended code (MCNPX) was used for the shielding calculations because of its capability to perform electron-, photon-, and neutron-coupled transport simulations. The photon dose was tallied using the MCNPX calculation, starting with the leaked electrons. However, it is difficult to accurately tally the neutron dose directly from the leaked electrons. The neutron yield per electron from the interactions with the surrounding components is very small, similar to 0.01 neutron for 100-MeV electron and even smaller for lower-energy electrons. This causes difficulties for the Monte Carlo analyses and consumes tremendous computation resources for tallying the neutron dose outside the shield boundary with an acceptable accuracy. To avoid these difficulties, the SOURCE and TALLYX user subroutines of MCNPX were utilized for this study. The generated neutrons were banked, together with all related parameters, for a subsequent MCNPX calculation to obtain the neutron dose. The weight windows variance reduction technique was also utilized for both neutron and photon dose calculations. Two shielding materials, heavy concrete and ordinary concrete, were considered for the shield design. The main goal is to maintain the total dose outside the shield boundary less than 5.0e-03 mSv/h during operation. The shield configuration and parameters of the accelerator building were determined and are presented in this paper. Copyright (C) 2016, Published by Elsevier Korea LLC on behalf of Korean Nuclear Society.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Analysis of soot formation from aviation fuels in laminar counterflow flames

Combustion emissions from aviation contribute to the formation of condensation trail (contrail) that can lead to the formation of anthropogenic cirrus clouds. Ice particles that form contrails are observed to have a linear correlation with soot particle number density. Synthetic aviation fuels (SAFs) offer a promising route to mitigate the production of soot particles while also increasing energy security. Although studies have focused on combustion and spray behavior, the detailed investigation of soot formation processes for different jet fuels and their impact on models for computational fluid dynamics (CFD) applications is not well understood. Moreover, experimental measurements of soot for canonical flames using Synthetic aviation fuels (SAF) for model validation remain scarce. To address this, we use employed the Lawrence Livermore National Laboratory (LLNL) detailed soot model based on the discrete sectional method. Additionally, we develop two reduced chemical mechanisms for Jet-A and Alcohol-to-Jet (C1) that are suitable for turbulent flame simulations and couple them with the Hybrid Method of Moments (HMOM). The detailed and reduced model frameworks are validated against experimental measurements of soot volume fraction (ƒ ν ) from a counterflow burner experiment previously reported in the literature. Given the good agreement between modeling results and experimental measurements for the (1) spatial distribution of ƒ ν and (2) the non-linear variation of peak ƒ ν with strain rate, we further investigate the modeled sub-processes (nucleation, condensation, surface growth, and oxidation) using the LLNL model to analyze the assumptions in the reduced model framework. Furthermore, the results indicate a significant contribution from resonant radicals to the surface growth of soot particles, which are not accounted for in the current implementation of HMOM and could help reconcile soot predictions by the reduced model with observations.

Counterflow↗
DOE OSTI2021

Concentration-Dependent Solvation Structure and Dynamics of Aqueous Sulfuric Acid Using Multinuclear NMR and DFT

Sulfuric acid is a ubiquitous compound for industrial processes and aqueous sulfate solutions also play a critical role in electrolytes for many prominent battery chemistries. While the thermodynamic literature for it is quite well-developed, comprehensive studies of the solvation structure, and particularly molecular-scale dynamical and transport properties, are less available. This study applies a multinuclear NMR approach to the elucidation of the solvation structure and dynamics over a wide temperature (-10 – 50?C) and concentration (0 – 18 M) range, combining 17O shift, linewidth, and T1 relaxation measurements, 33S shift and linewidth measurements, and 1H PFG-NMR measurements of the proton self-diffusivity. In conjunction, these results indicate a crossover between two regimes of solvation structure and dynamics, occurring above the concentration associated with the deep eutectic point (~ 4.5 M), with the high-concentration regime dominated by strong water-sulfate correlation. This description was borne out in detail by the activation energy trends with increasing concentration derived from the relaxation of both the H2O/H3O+ and H2SO4/HSO4-/SO42- 17O resonances and the 1H self-diffusivity. However, the 17O chemical shift difference between the H2O/H3O+ and H2SO4/HSO4-/SO42- resonances across the entire temperature range is nevertheless strikingly linear. A computational approach coupling MD simulations and DFT NMR shift calculations to reproduce this trend is presented, which will be the subject of further development. This combination of multinuclear, dynamical NMR and computational methods, and the results furnished by this study, will provide a platform for future studies on battery electrolytes where aqueous sulfate chemistry plays a central role in the solution structure.

Bazak, Jonathan D.↗
DOE OSTI2020

Selecting Adsorbents to Separate Diverse Near-Azeotropic Chemicals

Industrial separations of near-azeotropic chemicals, species with very similar boiling points, are energy- and capital-intensive. Adsorption-based processes can energy-efficiently separate near-azeotropic mixtures provided suitable adsorbent materials can be found. Among the full diversity of industry-relevant molecules, millions of these mixtures exist, meaning that discovery of mixture-specific adsorbents by direct experiment is infeasible. We show that vast numbers of adsorbents and adsorbing molecules can be explored in a powerful way by coupling atomistic simulations with machine learning. This concept is demonstrated by describing the adsorption of ~54000 industry-relevant chemicals in an experimentally derived set of thousands of metal–organic framework materials. Here, our results identify thousands of near-azeotropic mixtures that can be efficiently separated using adsorption and open possibilities for creating adsorption processes for complex mixtures with many components.

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

Role of Well Operations and Multiphase Geomechanics in Controlling Fault Stability During CO 2 Storage and Enhanced Oil Recovery

CO 2 injection in active oil fields is a technology proposed for industrializing carbon storage operations. However, monitoring CO 2 migration, oil recovery, and mechanical deformation within caprock and on faults is a challenging problem in large-scale storage-enhanced oil recovery operations due to intersecting faults, multiphase flow effects, and wells with complex production-injection schedule. We develop and demonstrate a methodology based on our coupled modeling framework to monitor the movement of CO 2 , hydrocarbons and water, and the associated evolution in mechanical stability of faults during CO 2 storage-enhanced oil recovery in the Farnsworth Unit oil field in Texas, United States. The methodology honors geological, geophysical, and production-injection data acquired in the field over six decades. Differential depletion from hydraulically isolated fault compartments followed by water and CO 2 injection-induced overpressure causes volumetric contraction and expansion of the reservoir, and changes in the total and effective stresses in the overburden-reservoir-underburden complex. CO 2 migrates upward to accumulate near top of the geologic structure, and water migrates downward to pressurize the faults. Three-dimensional changes in the pressure and stress fields in the system lead to changes in the shear and effective normal tractions on three major faults compartmentalizing the field. Evolution in fault tractions is used to compute the evolution in the Coulomb failure function of the faults to quantify induced-slip tendency under production and injection. We explain the spatial heterogeneity and time variability of Coulomb failure function in terms of well location heterogeneity and well rate variability.

58 GEOSCIENCES↗
DOE OSTI2020

The Sensitivity of Hyporheic Exchange to Fractal Properties of Riverbeds

Hyporheic exchange in riverbeds is driven by current-bed topography interactions. Because riverbeds exhibit topographic roughness across scales, from individual grains to bedforms and bars, they can exhibit fractal patterns. This study analyzed the influence of fractal properties of riverbed topography on hyporheic exchange. A set of synthetic fractal riverbeds with different scaling statistics was used as inputs to sequentially coupled numerical simulations of turbulent channel flow and hyporheic flow. In the analysis, the maximum power spectrum (dune size) and the fractal dimension (topographic complexity) were considered as independent variables and we then investigated how interfacial fluxes and hyporheic travel times are functionally related to these variables. As the maximum power spectrum increases (i.e., dune height to flow depth ratio), the average interfacial flux increases logarithmically whereas it increases exponentially with an increase in fractal dimension. Hyporheic exchange is more sensitive to additional roughness (larger fractal dimensions) than to bedform size (larger maximum power). Our results imply that fractal properties of riverbeds are crucial to predicting hyporheic exchange. The predictive relationships we propose could be integrated with reduced complexity, large-scale models. Further, they can also be used to design artificial topographies that target hyporheic ecosystem services.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hyporheic Exchange in Sand Dunes Under a Freely Deforming River Water Surface

Bedform-driven hyporheic exchange is conditioned by the head gradients at the sediment-water interface. Local exchange phenomena between the surface and the subsurface are often driven by the dynamic forces related to the velocity distributions around sediment waves. In open channels, the static forces are represented by the water depth, yet most computational fluid dynamics models use a rigid-lid approximation. Here, we investigated whether and when the deformation of the river’s free-surface influences bedform-driven hyporheic exchange. This was done through simulations of coupled open channel and hyporheic flows with the air-water interface modeled either as a free-surface or a rigid-lid across increasing subcritical Froude numbers from ~0.05 to ~0.95. The normalized hyporheic flux was higher when considering free-surface deformation across most of the range of Froude numbers considered. When the Froude number was larger than 0.6 and smaller than 0.85, both hydrostatic and nonhydrostatic-driven fluxes increased significantly compared to the rigid-lid approach and the total flux was about five times that predicted with the rigid-lid. These results indicate that bedform-driven hyporheic fluxes are underestimated in most studies that use a rigid-lid assumption in subcritical flows.

58 GEOSCIENCES↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Use of Shallow Ice Core Measurements to Evaluate and Constrain 1980–1990 Global Reanalyses of Ice Sheet Precipitation Rates

Abstract Sea‐level rise (SLR) projections by Earth System Models (ESMs) depend on ice sheet surface mass balances. Accurate, global atmosphere reanalyses would be ideal for providing equilibrated ice sheet model initial conditions in fully coupled ESM simulations. Here we present the first evaluation of 1980–1990 global reanalysis precipitation over Greenland and Antarctica that uses independent observations of net accumulation rates derived from shallow ice cores. Precipitation distributions from both the European Centre for Medium‐Range Weather Forecast's Reanalysis (ERA5) and the Modern‐Era Retrospective Analysis for Research and Applications (MERRA‐2) are highly correlated with contemporaneous co‐located net accumulation rates from Greenland ( r 2 > 0.95) and West Antarctica ( r 2 > 0.7). Three other commonly used reanalyses (WFDE5, CRUNCEP, and GSWP3) exhibit significantly weaker correlations on one or both ice sheets. Our findings imply that ESMs should use ERA5 or MERRA‐2 in data‐forced simulations to validate ice sheet model dynamics and precondition firn for SLR projections.

Schneider, Adam↗
DOE OSTI2023

Sensitivity of Arctic Surface Temperature to Including a Comprehensive Ocean Interior Reflectance to the Ocean Surface Albedo Within the Fully Coupled CESM2

Abstract Almost all current climate models simplify the ocean surface albedo (OSA) by assuming the reflected solar energy without the ocean interior contribution. In this study, an improved ocean surface albedo scheme is incorporated into the Community Earth System Model version 2 (CESM2) to assess the sensitivity of Arctic surface temperature to including ocean interior reflectance to the OSA. Fully coupled CESM2 simulations with and without ocean interior reflectance are subsequently performed, we focus on the analysis of Arctic surface temperature responses. Incorporating ocean interior reflectance increases absorbed solar radiation and warms the ocean, enhancing seasonal heat storage and release across the Arctic Ocean, and increasing sea ice reduction and positive climate feedbacks that elevates Arctic surface temperature. Seasonal variations in air‐surface temperature differences induce changes in turbulent heat flux patterns, concurrently modifying dynamic advection and moisture processes that affect boundary layer humidity and low clouds, especially in winter. Based on partitioning physical processes in the thermodynamic energy equation, surface air warming is induced primarily through positive heating anomalies of vertical advection, latent heat release, and longwave radiative forcing. Through an examination of the surface energy budget, skin temperature warming is driven predominantly by increased downward longwave radiation, positive surface albedo feedback in summer, and increased conductive heat transport from the ocean particularly in winter. Significant effects of ocean interior reflectance on the Arctic Ocean, including sea surface warming and sea ice reduction, justify the importance of ocean interior reflectance in climate models for better understanding of ongoing Arctic climate changes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2017

Structural responses of metallic glasses under neutron irradiation

Seeking nuclear materials that possess a high resistance to particle irradiation damage is a long-standing issue. Permanent defects, induced by irradiation, are primary structural changes, the accumulation of which will lead to structural damage and performance degradation in crystalline materials served in nuclear plants. In this work, structural responses of neutron irradiation in metallic glasses (MGs) have been investigated by making a series of experimental measurements, coupled with simulations in ZrCu amorphous alloys. It is found that, compared with crystalline alloys, MGs have some specific structural responses to neutron irradiation. Although neutron irradiation can induce transient vacancy-like defects in MGs, they are fully annihilated after structural relaxation by rearrangement of free volumes. In addition, the rearrangement of free volumes depends strongly on constituent elements. In particular, the change in free volumes occurs around the Zr atoms, rather than the Cu centers. This implies that there is a feasible strategy for identifying glassy materials with high structural stability against neutron irradiation by tailoring the microstructures, the systems, or the compositions in alloys. Furthermore, this work will shed light on the development of materials with high irradiation resistance.

36 MATERIALS SCIENCE↗