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 1,063 records · Page 59

Self-induced temperature correction for inter-phase heat transfer in Euler-Lagrange point-particle simulation

In a two-way coupled Euler-Lagrange simulation, particles are approximated as point sources and their momentum and energy exchange with the surrounding flow are modeled through point-particle force and heat transfer models. As the particle size increases and approaches the Eulerian grid, the feedback force and heat transfer become large and strongly affect the local flow at the particle location. The fluid velocity and temperature computed in the EL simulation when interpolated to the particle location will be substantially different from the corresponding undisturbed values. In this work we will follow the approach pursued in Balachandar et al. (2018) [1] for self-induced velocity correction to develop an analogous correction procedure for self-induced temperature. We obtain analytical solutions for the self-induced thermal perturbation for a steady uniform cross flow past a Gaussian filtered heat source under both steady and unsteady conditions. The analytical results are extended to finite Peclet number using numerical simulations. The resulting quasi-steady and unsteady models are tested with a simple problem of exponential thermal evolution against the corresponding EL simulation results. We also provide a simple criterion for when the self-induced temperature correction is needed in an Euler-Lagrange simulation.

42 ENGINEERING↗

Coupled CFD and chemical-kinetics simulations of cellulosic-biomass enzymatic hydrolysis: Mathematical-model development and validation

Computational simulations of cellulosic-biomass enzymatic hydrolysis using a computational-fluid-dynamics (CFD) model coupled to a chemical-kinetics model are presented in this work. A time-step sub-cycling strategy was used to circumvent large disparities in time scales with respect to transport physics and chemical kinetics. The unknown reaction and transport parameters in the mathematical model were first fit to two canonical experiments: a well-mixed case that was limited by reaction kinetics and a settled-solids case that was limited by transport. Here, the model was then validated against a benchtop experiment for an intermediate-mixing scenario, for which simulation results were in reasonable agreement with experimentally measured cellulose conversions. Settling of cellulose particles and diffusion of enzymes determine overall conversion rates for lower mixing rates, while reaction rates dominate conversion at higher mixing rates for which greater homogenization of substrates is achieved.

09 BIOMASS FUELS↗

Modeling and Simulation of Deepwater Pipeline S-Lay With Coupled Dynamic Positioning

Dynamic position (DP) control and pipeline dynamics are the two main parts of the deepwater S-lay simulation model. In this study, a fully coupled analysis tool for deepwater S-lay deployment by dynamically positioned vessels is developed. The method integrates the major aspects related to numerical simulation, including coupled pipeline motion and roller contact forces. The roller–pipe interaction is incorporated in the S-lay pipeline model using a contact search method based on a lumped-mass (LM) formulation in global coordinates. A proportional-integration-differentiation (PID) controller and a Kalman filter are applied in the vessel motion equation to calculate the thrust allocation of the DP system in time domain. Numerical simulation results showed that the dynamic effects add a significant contribution to the tension, but have little influence on the maximum pipe stress and strain. The dynamic response of the coupled S-lay and DP pipeline deployment system increases the demand on the tensioner load carrying capability as well as the maximum DP thruster power.

Engineering↗

Root Dynamics Mitigate Warm and Dry Biases over the Central United States

The central United States frequently exhibits warm and dry biases in simulations of summertime conditions, a persistent feature that remains unresolved. While previous studies linked these biases to misrepresented surface energy exchanges, the role of belowground processes remains poorly understood. Here, we demonstrate that inadequate representation of root water uptake in land surface models contributes to this bias. Using both offline Noah-MP and coupled WRF-Noah-MP simulations with static and dynamic root water uptake schemes, we show that the inclusion of dynamic root processes reduces 2-m air temperature biases and enhances precipitation, primarily by increasing the rain rate of convective systems. Offline and coupled simulations further reveal that the cooling effects and precipitation increases are amplified through positive land-atmosphere feedback, active only in the coupled model. These findings highlight an important role of root in modulating land-atmosphere interactions and underscore the need to refine root-zone processes to improve regional atmospheric simulations.

Yang, Zhao (ORCID:0000000288027130)↗

Coupled Calculations in Helical Steam Generator: Validation on Legacy Data

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. Flow-induced vibration (FIV) is widespread problem in energy systems because they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occurs. Given the importance of reliable components in the nuclear industry, flow-induced vibration has long been a major concern in safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from flow-induced vibration and related failures. Advanced reactors, such as integral Pressurized Water Reactors (PWRs) considered for Small Modular Reactors (SMR), often rely on innovative component designs to meet cost and safety targets. One component that is the subject of advanced designs is the steam generator, some designs of which forego the usual shell-and-tube architecture in order to fit within the primary vessel. In addition to being more cost- and space-efficient, such steam generators need to be more reliable, since failure of the primary vessel represents a potential loss of coolant and a safety concern. A significant amount of data exists on flow-induced vibration in shell-and-tube heat exchangers, and heuristic methods are available to predict their occurrence based on a set of given assumptions. In contrast, advanced designs have far less data available. Advanced modeling and simulation based on coupled structural and fluid simulations have the potential to predict flow-induced vibration in a variety of designs, reducing the need for expensive experimental programs, especially at the design stage. Over the past five years, the Reactor Product Line has developed the integrated multi-physics code suite SHARP. The goal of developing such a tool is to perform multi-physics neutronics, thermal/fluid, and structural mechanics modeling of the components inside the full reactor core or portions of it with a user-specified fidelity. In particular SHARP contains high-fidelity single-physics codes Diablo for structural mechanics and Nek5000 for fluid mechanics calculations. Both codes are state-of-the-art, highly scalable tools that have been extensively validated. These tools form a strong basis on which to build a flow-induced vibration modeling capability. In this report we discuss one-way coupled calculations performed with Nek5000 and Diablo aimed at simulating available FIV experiments in helical steam generators in the turbulent buffeting regime. In this regime one-way coupling is judged sufficient because the pressure loads do not cause substantial displacements. It is also the most common source of vibration in helical steam generators at the low flows expected in integral PWRs. The legacy data is obtained from two datasets developed at Argonne and B&W.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupled 2D MHD and runaway electron fluid simulations of SPARC disruptions

Runaway electrons (REs) generated during disruption events in tokamaks can carry mega-Ampère level currents, potentially causing damage to plasma-facing components. Understanding RE evolution during disruption events is important for evaluating strategies to mitigate RE damage. Using two-dimensional toroidally symmetric magnetohydrodynamic (MHD) simulations in M3D-C1, which incorporates a fluid RE model evolved self-consistently with the bulk MHD fluid, we examine the seeding and avalanching of REs during disruptions in the SPARC tokamak – a compact, high-field, high-current device designed to achieve a fusion gain Q > 2 in deuterium–tritium plasmas. The M3D-C1 simulations of unmitigated disruptions demonstrate RE plateau formation and peaking of the final current density, which agree well with the results of lower-fidelity reduced RE fluid models. This work provides the first systematic comparison and benchmarking of different primary sources, including activated tritium beta decay and Compton scattering, in SPARC disruption simulations with self-consistent MHD and RE coupling.

Compton scattering↗

Future changes in regional precipitation simulated by a half-degree coupled climate model: Sensitivity to horizontal resolution

The global fully coupled half-degree Community Climate System Model Version 4 (CCSM4) was integrated for a suite of climate change ensemble simulations including five historical runs, five Representative Concentration Pathway 8.5 [RCP8.5) runs, and a long Pre-Industrial control run. This study focuses on precipitation at regional scales and its sensitivity to horizontal resolution. The half-degree historical CCSM4 simulations are compared to observations, where relevant, and to the standard 1° CCSM4. Both the halfdegree and 1° resolutions are coupled to a nominal 1° ocean. North American and South Asian/Indian monsoon regimes are highlighted because these regimes demonstrate improvements due to higher resolution, primarily because of better-resolved topography. Agriculturally sensitive areas are analyzed and include Southwest, Central, and Southeast U.S., Southern Europe, and Australia. Both mean and extreme precipitation is discussed for convective and large-scale precipitation processes. Convective precipitation tends to decrease with increasing resolution and large-scale precipitation tends to increase. Improvements for the half-degree agricultural regions can be found for mean and extreme precipitation in the Southeast U.S., Southern Europe, and Australian regions. Climate change responses differ between the model resolutions for the U.S. Southwest/Central regions and are seasonally dependent in the Southeast and Australian regions. Both resolutions project a clear drying signal across Southern Europe due to increased greenhouse warming. As a result, differences between resolutions tied to the representation of convective and large-scale precipitation play an important role in the character of the climate change and depend on regional influences.

58 GEOSCIENCES↗

Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes

Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.

Astronomy & Astrophysics↗

Demonstration of Fully Coupled Calculations in Helical Steam Generator: Toward Predictions of Fluid Elastic Instability

The NEAMS program aims to develop an integrated multi-physics simulation capability “pellet-to-plant” for the design and analysis of future generations of nuclear power plants. In particular, the Reactor Product Line code suite's multi-resolution hierarchy is being designed to ultimately span the full range of length and time scales present in relevant reactor design and safety analyses, as well as scale from desktop to petaflop computing platforms. Flow-induced vibration (FIV) is widespread problem in energy systems because they rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occurs. Given the importance of reliable components in the nuclear industry, flow-induced vibration has long been a major concern in safety and operation of nuclear reactors. In particular, nuclear fuel rods and steam generators have been known to suffer from flow-induced vibration and related failures. Advanced reactors, such as integral Pressurized Water Reactors (PWRs) considered for Small Modular Reactors (SMR), often rely on innovative component designs to meet cost and safety targets. One component that is the subject of advanced designs is the steam generator, some designs of which forego the usual shell-and-tube architecture in order to fit within the primary vessel. In addition to being more cost- and space-efficient, such steam generators need to be more reliable, since failure of the primary vessel represents a potential loss of coolant and a safety concern. A significant amount of data exists on flow-induced vibration in shell-and-tube heat exchangers, and heuristic methods are available to predict their occurrence based on a set of given assumptions. In contrast, advanced designs have far less data available. Advanced modeling and simulation based on coupled structural and fluid simulations have the potential to predict flow-induced vibration in a variety of designs, reducing the need for expensive experimental programs, especially at the design stage. Over the past seven years, the Reactor Product Line has developed the integrated multi-physics code suite SHARP. The goal of developing such a tool is to perform multi-physics neutronics, thermal/fluid, and structural mechanics modeling of the components inside the full reactor core or portions of it with a user-specified fidelity. In particular SHARP contains high-fidelity single-physics codes Diablo for structural mechanics and Nek5000 for fluid mechanics calculations. Both codes are state-of-the-art, highly scalable tools that have been extensively validated. These tools form a strong basis on which to build a flow-induced vibration modeling capability. In the last two years a significant effort was devoted to 1) develop and validate a fully coupled capability in SHARP for Fluid-Structure calculations (FSI) and 2) validate the use of SHARP for one-way coupled calculations in helical steam generators on legacy data available as well as novel experiments. In this report we discuss an extension of previous one-way coupled calculations performed with Nek5000 and Diablo aimed at simulating available FIV experiments in helical steam generators in the turbulent buffeting regime. It was demonstrated that in the buffeting regime one-way coupling was judged sufficient because the pressure loads do not cause substantial displacements. However, higher speeds led to a reduction in accuracy and the need for a fully coupled capability for velocity or approaching the critical velocity and fluid elastic instability. In this report we demonstrate such fully coupled capability on the Argonne test case discussed in previous reports.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Constructive Quantum Interference in Photochemical Reactions

Interferences emerge when multiple pathways coexist together, leading toward the same result. Here, we report a theoretical study for a reaction scheme that leads to constructive quantum interference in a photoassociation (PA) reaction of a 87 Rb Bose–Einstein condensate where the reactant spin state is prepared in a coherent superposition of multiple bare spin states. This is achieved by changing the reactive scattering channel in the PA reaction. As the origin of coherent control comes from the spin part of the wavefunction, we show that it is sufficient to use radio frequency (RF) coupling to achieve the superposition state. We simulate the RF coupling on a quantum processor (IBMQ Lima), and our results show that interferences can be used as a resource for the coherent control of photochemical reactions. The approach is general and can be employed to study a wide spectrum of chemical reactions in the ultracold regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exactly solvable model of light-scattering errors in quantum simulations with metastable trapped-ion qubits

Here, we analytically solve a model for light scattering in Ising dynamics of metastable atomic qubits, generalizing the approach of Foss-Feig et al. to include leakage outside the qubit manifold. We analyze the influence of these fundamental errors in simulations of proposed experiments with metastable levels of 40 Ca + ions in a Penning trap. We find that “effective magnetic fields” generated by leaked qubits have significant impacts on spin-spin correlation functions for Greenberger-Horne-Zeilinger state preparation or for quantum simulations with strong coupling, while spin squeezing uses a much weaker coupling and is largely insensitive to the simulated leakage errors, even with a few hundred ions. Our theory and results are expected to be useful in modeling a variety of metastable qubit experiments in the future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Failure detection by pilots during automatic landing: Models and experiments

A model of the pilot as a monitor of instrument failures during automatic landing is proposed. The failure detection model consists of two stages: a linear estimator (Kalman Filter) and a decision mechanism which is based on sequential analysis. The filter equations are derived from a simplified version of the linearized dynamics of the airplane and the control loop. The perceptual observation noise is modelled to include the effects of the partition of attention among the several instruments. The final result is a simple model consisting of a high pass filter to produce the observation residuals, and a decision function which is a pure integration of the residuals minus a bias term. The dynamics of a Boeing 707 were used to simulate the fully coupled final approach in a fixed base simulator which also included failures in the airspeed, glideslope, and localizer indicators. Subjects monitored the approaches and detected the failures; their performance was compared with the predictions of the model with good agreement between the experimental data and the model.

Gai, E. G.↗

VHDL simulation with access to transistor models

Hardware description languages such as VHDL have evolved to aid in the design of systems with large numbers of elements and a wide range of electronic and logical abstractions. For high performance circuits, behavioral models may not be able to efficiently include enough detail to give designers confidence in a simulation's accuracy. One option is to provide a link between the VHDL environment and a transistor level simulation environment. The coupling of the Vantage Analysis Systems VHDL simulator and the NOVA simulator provides the combination of VHDL modeling and transistor modeling.

Gibson, J.↗

Simulation of the photo-luminescent coupling transfer function in multijunction nanoheterostructure solar cells

A mathematical model of recombination radiation transfer in multilayer heterostructures for calculating photoluminescent coupling transfer function has been developed. The photoluminescent coupling effectiveness for the GaInP-GaInAs and GaInAs-Ge subcell pairs being a part of multijunction GaInP/GaInAs/Ge SCs has been investigated. The simulation results were compared with experimental data. The electroluminescence effectiveness estimation for a wide band gap subcell in each pair has been performed. An anomalously high value of the photoluminescent coupling for the GaInAs-Ge subcell pair has been observed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Performance of Coupled Physics Solvers for Multidisciplinary Hypersonic Flow Simulations on Several Classes of Computer Architectures

The application of hypersonic flow simulation tools to realistic flight scenarios will require the coupling of multiple physical effects to the baseline fluid dynamics. Such multiphysics effects can include the aerooelastic response of the airframe or engine components, dynamic transport of atmospheric particles, the deformation of solid-fluid interfaces that can ablate, pyrolyze, or erode, as well as a host of other processes, all of which are governed by unique sets of physical equations and models. Coupling multiple (and potentially disparate) physics solvers to a robust compressible flow solver poses additional challenges related to the stability, performance and scalability of the combined solver. The choices made during the software design process can therefore lead to a variation in simulation efficiency across different computer architectures. In this paper, we will consider two representative multiphysics hypersonic flow scenarios: the interaction of solid particulates with the flow field created by a hypersonic lifting body and the aerooelastic deformation of a model airframe under high-Mach-number flow conditions. For these simulations we explore the behavior of several hypersonic simulation tools, including Kestrel, FUN3D, US3D, and JENRE multiphysics framework, on several high performance computing systems containing various CPU and GPU architectures.

architecture↗

Quantifying the agreement between observed and simulated extratropical modes of interannual variability

Using historical simulations of the Coupled Model Intercomparison Project-5 (CMIP5) and multiple observationally-based datasets, we employ skill metrics to analyze the fidelity of the simulated Northern Annular Mode, the North Atlantic Oscillation, the Pacific North America pattern, the Southern Annular Mode, the Pacific Decadal Oscillation, the North Pacific Oscillation, and the North Pacific Gyre Oscillation. We assess the benefits of a unified approach to evaluate these modes of variability, which we call the common basis function (CBF) approach, based on projecting model anomalies onto observed empirical orthogonal functions (EOFs). The CBF approach circumvents issues with conventional EOF analysis, eliminating, for example, corrections of arbitrarily assigned, but inconsistent, signs of the EOF’s/PC’s being compared. It also avoids the problem that sometimes the first observed EOF is more similar to a higher order model EOF, particularly if the simulated EOFs are not well separated. Compared to conventional EOF analysis of models, the CBF approach indicates that models compare significantly better with observations in terms of pattern correlation and root-mean-squared-error (RMSE) than heretofore suggested. In many cases, models are doing a credible job at capturing the observationally-based estimates of patterns; however, errors in simulated amplitudes can be large and more egregious than pattern errors. In the context of the broad distribution of errors in the CMIP5 ensemble, sensitivity tests demonstrate that our results are relatively insensitive to methodological considerations (CBF vs. conventional approach), observational uncertainties in pattern (as determined by using multiple datasets), and internal variability (when multiple realizations from the same model are compared). The skill metrics proposed in this study can provide a useful summary of the ability of models to reproduce the observed EOF patterns and amplitudes. Additionally, the skill metrics can be used as a tool to objectively highlight where potential model improvements might be made. Finally, we advocate more systematic and objective testing of simulated extratropical variability, especially during the non-dominant seasons of each mode, when many models are performing relatively poorly.

58 GEOSCIENCES↗

Simulation of Gas Production from Multilayered Hydrate-Bearing Media with Fully Coupled Flow, Thermal, Chemical and Geomechanical Processes Using TOUGH+Millstone. Part 3: Production Simulation Results

The TOUGH+Millstone simulator has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH 4 - /sub>-hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH+Hydrate simulator, v2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the evolution of hydrate-bearing systems and includes the most recent physical properties relationships, coupled seamlessly with (b) Millstone v1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of the TOUGH+Millstone code are demonstrated in the simulation and analysis of the system flow, thermal, and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the third paper of this series, we apply the simulators described in parts 1 and 2 to a problem of gas production from a complex, multilayered system of hydrate-bearing sediments in an oceanic environment. We perform flow simulations of constant-pressure production via a vertical well and compare those results to a coupled flow-geomechanical simulation of the same process. Finally, the results demonstrate the importance of fully coupled geomechanics when modeling the evolution of reservoir properties during production.

58 GEOSCIENCES↗