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

Thermodynamic and magnetic properties of a three-state Potts model on a triangular lattice with next-neighbor interactions

Using Monte Carlo calculations, we study the thermodynamic and magnetic properties of 2D structures that can be described by a three-state Potts model on a triangular lattice with nearest- and secondnearest- neighbor interactions characterized by the coupling constants J{sub 1} and J{sub 2}, respectively. Analyzing the thermodynamic parameters of heat capacity, the order parameter, the susceptibility, and fourth-order Binder cumulants, we show that the three-state Potts model with coupling constants J{sub 1} > 0 and J{sub 2} < 0 predicts a phase transition of the second kind for 0 ≤ |r| ≤ 1/3, where r = J{sub 2}/J{sub 1}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Coupling DAEM and CFD for simulating biomass fast pyrolysis in fluidized beds

We report results from computational simulations of an experimental, lab-scale bubbling bed biomass pyrolysis reactor that include a distributed activation energy model (DAEM) for the kinetics. In this study, we utilized multiphase computational fluid dynamics (CFD) to account for the turbulent hydrodynamics, and this was combined with the DAEM kinetics in a multi-component, multi-step reaction network. Our results indicate that it is possible to numerically integrate the coupled CFD–DAEM system without significantly increasing computational overhead. It is also clear, however, that reactor operating conditions, reaction kinetics, and multiphase flow dynamics all have major impacts on the pyrolysis products exiting the reactor. We find that, with the same pre-exponential factors and mean activation energies, inclusion of distributed activation energies in the kinetics can shift the predicted average value of the exit vapor-phase tar flux and its statistical distribution, compared to single-valued activation-energy kinetics. Perhaps the most interesting observed trend is that increasing the diversity of the DAEM activation energies appears to increase the mean tar yield, all else being equal. As a result, these findings imply that accurate resolution of the reaction activation energy distributions will be important for optimizing biomass pyrolysis processes.

09 BIOMASS FUELS↗
DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗
DOE OSTI2018

Verification and Validation of Simulation Framework for Analysis of Traumatic Brain Injury

Many U.S. troops deployed in Iraq and Afghanistan have suffered blast-related, closed-head injuries as a result of being within non-lethal distance of a detonated improvised explosive device (IED). Consequently, we are conducting a modeling & simulation investigation, coupled with clinical diagnostics of brain injury, into the mechanisms associated with blast to the human head that give rise to traumatic brain injury (TBI). In particular, we are attempting to identify the specific conditions of focused stress wave energy within the brain, resulting from blast loading to the head, which define a threshold for persistent brain injury. As these threshold injury conditions are identified, one can then employ a modeling & simulation toolset to assess the efficacy of various helmet designs to mitigate blast exposure that would otherwise lead to traumatic brain injury. In this report, we describe our efforts to establish the role of stress wave interactions in the genesis of traumatic brain injury (TBI) from exposure to explosive blast. We discuss our effort to correlate simulation predictions with clinical assessments of TBI subjects whose brain injuries have been diagnosed by advanced magnetic resonance imaging techniques. This work has been funded through the U.S. Naval Health Research Center, Office of Naval Research, Mr. James Mackiewicz, project funding manager.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2016

Impact of Soil Moisture Assimilation on Land Surface Model Spin-Up and Coupled LandAtmosphere Prediction

Advances in satellite monitoring of the terrestrial water cycle have led to a concerted effort to assimilate soil moisture observations from various platforms into offline land surface models (LSMs). One principal but still open question is that of the ability of land data assimilation (LDA) to improve LSM initial conditions for coupled short-term weather prediction. In this study, the impact of assimilating Advanced Microwave Scanning Radiometer for EOS (AMSR-E) soil moisture retrievals on coupled WRF Model forecasts is examined during the summers of dry (2006) and wet (2007) surface conditions in the southern Great Plains. LDA is carried out using NASAs Land Information System (LIS) and the Noah LSM through an ensemble Kalman filter (EnKF) approach. The impacts of LDA on the 1) soil moisture and soil temperature initial conditions for WRF, 2) land-atmosphere coupling characteristics, and 3) ambient weather of the coupled LIS-WRF simulations are then assessed. Results show that impacts of soil moisture LDA during the spin-up can significantly modify LSM states and fluxes, depending on regime and season. Results also indicate that the use of seasonal cumulative distribution functions (CDFs) is more advantageous compared to the traditional annual CDF bias correction strategies. LDA performs consistently regardless of atmospheric forcing applied, with greater improvements seen when using coarser, global forcing products. Downstream impacts on coupled simulations vary according to the strength of the LDA impact at the initialization, where significant modifications to the soil moisture flux- PBL-ambient weather process chain are observed. Overall, this study demonstrates potential for future, higher-resolution soil moisture assimilation applications in weather and climate research.

soil moisture↗
DOE OSTI2018

A method for measuring coupling coefficients between cores and corrections to mode propagation constants in multicore fibres

A new experimental method is proposed for reconstructing supermode profiles of a multicore fibre (MCF), small corrections to the propagation constants of its cores and coupling coefficients between the cores without phase-sensitive measurements. The method is based on intensity measurements in the cores at the output of several MCF segments differing in length, using several configurations of controlled light launch into the cores at the fibre input (in the simplest case, by sequentially exciting each core), followed by numerical treatment using an iterative algorithm. The use of the method is demonstrated experimentally for a seven-core fibre. Using numerical simulation, we assess the stability of the method and the effect of the number of measurements on the accuracy of fibre parameter reconstruction. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The mechanism for acetate formation in electrochemical CO (2) reduction on Cu: selectivity with potential, pH, and nanostructuring

Nanostructured Cu catalysts have increased the selectivities and geometric activities for high value C–C coupled (C 2 ) products (ethylene, ethanol, and acetate) in the electrochemical CO (2) reduction reaction (CO (2) RR). The selectivity among the high-value C 2 products is also altered, where for instance the yield of acetate increases with alkalinity and is dependent on the catalyst morphology. The reaction mechanisms behind the selectivity towards acetate vs. other C 2 products remain controversial. In this work, we elucidate the reaction mechanism for acetate formation by using ab initio simulations, a coupled kinetic-transport model, and loading dependent experiments. We find that trends in acetate selectivity can be rationalized from variations in electrolyte pH and the local mass transport properties of the catalyst and not from changes in Cu's intrinsic activity. The selectivity mechanism originates from the transport of ketene, a stable (closed shell) intermediate, away from the catalyst surface into solution where it reacts to form acetate. While this type of mechanism has not yet been discussed in the CO (2) RR, variants of it may explain similar selectivity fluctuations observed for other stable intermediates like CO and acetaldehyde. Our proposed mechanism suggests that acetate selectivity increases with increasing pH, decreasing catalyst roughness and significantly varies with the applied potential.

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

The Hybrid Diabatization Method and Its Application to the CaH Quasi-Molecule

A procedure for diabatization of quantum-chemical data is proposed that allows one to surmount the problem of -shaped peaks of nonadiabatic matrix elements by passing to a hybrid diabatic representation of a system of coupled channel equations of nuclear dynamics with the use of model nonadiabatic matrix elements that approach initially specified matrix elements. Results of the application of this method within the framework of the Landau–Zener model are presented for the example of the CaH collision system.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

What can simulation test beds teach us about social science? Results of the ground truth program

The ground truth program used simulations as test beds for social science research methods. The simulations had known ground truth and were capable of producing large amounts of data. This allowed research teams to run experiments and ask questions of these simulations similar to social scientists studying real-world systems, and enabled robust evaluation of their causal inference, prediction, and prescription capabilities. We tested three hypotheses about research effectiveness using data from the ground truth program, specifically looking at the influence of complexity, causal understanding, and data collection on performance. We found some evidence that system complexity and causal understanding influenced research performance, but no evidence that data availability contributed. The ground truth program may be the first robust coupling of simulation test beds with an experimental framework capable of teasing out factors that determine the success of social science research.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Design of 4-electrode optical device for application of vector electric fields to self-assembled quantum dot complexes

Self-assembled InAs quantum dots (QDs) are of great interest as components of optoelectronic devices that can operate at the quantum limit. The charge configuration, interdot coupling, and symmetry of complexes containing multiple QDs can all be tuned with applied electric fields, but the magnitude and angle of the electric field required to control each of these parameters depend on the orientation of the QD complex. We present a 4-electrode device compatible with optical excitation and emission that allows application of electric fields with arbitrary magnitudes and angles relative to isolated QD complexes. We demonstrate the electric field tunability of this device with numerical simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Explicit environmental and vibronic effects in simulations of linear and nonlinear optical spectroscopy

Accurately simulating the linear and nonlinear electronic spectra of condensed phase systems and accounting for all physical phenomena contributing to spectral line shapes presents a significant challenge. Vibronic transitions can be captured through a harmonic model generated from the normal modes of a chromophore, but it is challenging to also include the effects of specific chromophore–environment interactions within such a model. We work to overcome this limitation by combining approaches to account for both explicit environment interactions and vibronic couplings for simulating both linear and nonlinear optical spectra. We present and show results for three approaches of varying computational cost for combining ensemble sampling of chromophore–environment configurations with Franck–Condon line shapes for simulating linear spectra. We present two analogous approaches for nonlinear spectra. Simulated absorption spectra and two-dimensional electronic spectra (2DES) are presented for the Nile red chromophore in different solvent environments. Furthermore, employing an average Franck–Condon or 2DES line shape appears to be a promising method for simulating linear and nonlinear spectroscopy for a chromophore in the condensed phase.

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

Simulating the Effect of Fracture Connectivity on Repository Performance with GDSA Framework - 18589

This work examines the possible migration of radionuclides from a deep geologic nuclear waste repository sited in fractured crystalline host rock. The key safety concern is the potential effect on waste isolation of the inter-connectivity of the fracture network, which is primarily established by the temporal evolution of the temperature and stress fields at the time of original rock intrusion and cooling. Two end members are considered in these simulations, one with a high degree of connectivity to the biosphere, such that advective transport through the fracture and fault network controls radionuclide migration, and the other with a low degree of connectivity, such that slow diffusive transport through the crystalline rock matrix is the controlling process for radionuclide migration to the biosphere. Both end members have several areally extensive, high-transmissivity deformation zones, whose distance to the repository is an important factor for waste isolation capability. Uncertainties in fracture properties (transmissivity, orientation, radius) and fracture distribution also give rise to uncertainty in waste isolation capability (i.e., the overall connectivity of the repository with the surface biosphere). These types of natural system uncertainties are generally present in any geologic site-characterization program and thus represent an important factor in assessing repository performance in hard rock environments. The waste considered in these simulations are spent fuel rod assemblies from the U.S. commercial nuclear reactor fleet. Thermal output of the waste must be considered in the simulations and can influence the rate and timing of waste package failure, waste form degradation, and fluid flux due to thermal expansion around the repository horizon. To achieve a high-fidelity representation of radionuclide transport in fractures and rock matrix, combined with thermal energy transport and fluid flow in fractures and matrix, the mathematical model is solved numerically in a parallel high-performance computing (HPC) environment on a finite volume unstructured grid consisting of approximately 5 million cells, using Geologic Disposal Safety Assessment (GDSA) Framework (https://pa.sandia.gov), an open-source performance assessment tool for deep underground disposal of nuclear waste. GDSA Framework uses PFLOTRAN to solve the balance equations on a three-dimensional grid with heterogeneous properties, using multiple processors in a parallel configuration based on domain decomposition. The fracture networks in these simulations are originally generated as discrete fracture networks (DFNs), which are sets of two-dimensional planes distributed in a three-dimensional domain. The method used in GDSA Framework maps the stochastically generated DFN to an equivalent continuous porous medium (ECPM) domain that allows for the simulation of coupled heat flow, fluid flow, and radionuclide transport, including heat conduction through the matrix of the fractured rock, which is a process not easily modeled in a DFN representation. Computational efficiency is also greatly enhanced using the ECPM method, allowing for a realistic representation and analysis of uncertainties in a multi-realization performance assessment of a deep geologic repository. The effect of fracture connectivity on the waste isolation safety function, as brought to light by these GDSA Framework simulations, points to the importance of including a realistic representation of uncertainties in fracture properties and distribution (effectively, uncertainty in spatial heterogeneity) in repository safety assessment simulations. (authors)

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

Dark matter relics and the expansion rate in scalar-tensor theories

Here we study the impact of a modified expansion rate on the dark matter relic abundance in a class of scalar-tensor theories. The scalar-tensor theories we consider are motivated from string theory constructions, which have conformal as well as disformally coupled matter to the scalar. We investigate the effects of such a conformal coupling to the dark matter relic abundance for a wide range of initial conditions, masses and cross-sections. We find that exploiting all possible initial conditions, the annihilation cross-section required to satisfy the dark matter content can differ from the thermal average cross-section in the standard case. We also study the expansion rate in the disformal case and find that physically relevant solutions require a nontrivial relation between the conformal and disformal functions. We study the effects of the disformal coupling in an explicit example where the disformal function is quadratic.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Online transfer learning strategy for enhancing the scalability and deployment of deep reinforcement learning control in smart buildings

In recent years, advanced control strategies based on Deep Reinforcement Learning (DRL) proved to be effective in optimizing the management of integrated energy systems in buildings, reducing energy costs and improving indoor comfort conditions when compared to traditional reactive controllers. However, the scalability and implementation of DRL controllers are still limited since they require a considerable amount of time before converging to a near-optimal solution. This issue is currently addressed in literature through the offline pre-training of the DRL agent. However this solution results in two main critical issues: (1) the need to develop a building surrogate model to perform the training task, and (2) the need to perform a fine-tuning process over several training episodes to obtain a near-optimal control policy. In this context, this paper introduces an Online Transfer Learning (OTL) strategy that exploits two knowledge-sharing techniques, weight-initialization and imitation learning, to transfer a DRL control policy from a source office building to various target buildings in a simulation environment coupling EnergyPlus and Python. A DRL controller based on discrete Soft Actor–Critic (SAC) is trained on the source building to manage the operation of a cooling system consisting of a chiller and a thermal storage. Several target buildings are defined to benchmark the performance of the OTL strategy with that of a Rule-Based Controller (RBC) and two DRL-based control strategies, deployed in offline and online fashion. The strategy adopted for OTL emulates the real world implementation with a simulation process by implementing the transferred DRL agent for a single episode in the target buildings. Target buildings have the same geometrical features and are served by the same energy system as the source building, but differ in terms of weather conditions, electricity price schedules, occupancy patterns, and building envelope efficiency levels. The results show that the OTL strategy can reduce the cumulated sum of temperature violations on average by 50% and 80% respectively when compared to RBC and online DRL while enhancing the energy system operation with electricity cost savings ranging between 20% and 40%. Furthermore, the OTL agent performs slightly worse than the offline DRL controller but it does not require any modeling effort and can be implemented directly on target buildings emulating a real-world implementation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Bayesian calibration of strength parameters using hydrocode simulations of symmetric impact shock experiments of Al-5083

We report that predictive modeling of materials requires accurately parameterized constitutive models. Parameterizing models that describe dynamic strength and plasticity require experimentally probing materials in a variety of strain rate regimes. Some experimental protocols (e.g., plate impact) probe the constitutive response of a material using indirect measures such as free surface velocimetry. Manual efforts to parameterize constitutive models using indirect experimental measures often lead to non-unique optimizations without quantification of parameter uncertainty. This study uses a Bayesian statistical approach to find model parameters and to quantify the uncertainty of the resulting parameters. The technique is demonstrated by parameterizing the Johnson-Cook strength model for aluminum alloy 5083 by coupling hydrocode simulations and velocimetry measurements of a series of plate impact experiments. Simulation inputs and outputs are used to calibrate an emulator that mimics the outputs of the computationally intensive simulations. Varying the amount of experimental data available for emulator calibration showed clear differences in the degree of uncertainty and uniqueness of the resulting optimized Johnson-Cook parameters for Al-5083. The results of the optimization provided a numerical evaluation of the degree of confidence in model parameters and model performance. Lastly, given an understanding of the physical effects of certain model parameters, individual parameter uncertainty can be leveraged to quickly identify gaps in the physical domains covered by completed experiments.

36 MATERIALS SCIENCE↗