Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Coupling”

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 73 records · Page 4

Documenting numerical experiments in support of the Coupled Model Intercomparison Project Phase 6 (CMIP6)

Abstract. Numerical simulation, and in particular simulation of the earth system, relies on contributions from diverse communities, from those who develop models to those involved in devising, executing, and analysing numerical experiments. Often these people work in different institutions and may be working with significant separation in time (particularly analysts, who may be working on data produced years earlier), and they typically communicate via published information (whether journal papers, technical notes, or websites). The complexity of the models, experiments, and methodologies, along with the diversity (and sometimes inexact nature) of information sources, can easily lead to misinterpretation of what was actually intended or done. In this paper we introduce a taxonomy of terms for more clearly defining numerical experiments, put it in the context of previous work on experimental ontologies, and describe how we have used it to document the experiments of the sixth phase for the Coupled Model Intercomparison Project (CMIP6). We describe how, through iteration with a range of CMIP6 stakeholders, we rationalized multiple sources of information and improved the clarity of experimental definitions. We demonstrate how this process has added value to CMIP6 itself by (a) helping those devising experiments to be clear about their goals and their implementation, (b) making it easier for those executing experiments to know what is intended, (c) exposing interrelationships between experiments, and (d) making it clearer for third parties (data users) to understand the CMIP6 experiments. We conclude with some lessons learnt and how these may be applied to future CMIP phases as well as other modelling campaigns.

Pascoe, Charlotte (ORCID:0000000223616778)↗

Interpreting the nonlinear dielectric response of glass-formers in terms of the coupling model

Nonlinear dielectric measurements at high electric fields of glass-forming glycerol and propylene carbonate initially were carried out to elucidate the dynamic heterogeneous nature of the structural α-relaxation. Recently, the measurements were extended to sufficiently high frequencies to investigate the nonlinear dielectric response of faster processes including the so-called excess wing (EW), appearing as a second power law at high frequencies in the loss spectra of many glass formers without a resolved secondary relaxation. While a strong increase of dielectric constant and loss is found in the nonlinear dielectric response of the α-relaxation, there is a lack of significant change in the EW. A surprise to the experimentalists finding it, this difference in the nonlinear dielectric properties between the EW and the α-relaxation is explained in the framework of the coupling model by identifying the EW investigated with the nearly constant loss (NCL) of caged molecules, originating from the anharmonicity of the intermolecular potential. The NCL is terminated at longer times (lower frequencies) by the onset of the primitive relaxation, which is followed sequentially by relaxation processes involving increasing number of molecules until the terminal Kohlrausch α-relaxation is reached. These intermediate faster relaxations, combined to form the so-called Johari-Goldstein (JG) β-relaxation, are spatially and dynamically heterogeneous, and hence exhibit nonlinear dielectric effects, as found in glycerol and propylene carbonate, where the JG β-relaxation is not resolved and in D-sorbitol where it is resolved. Like the linear susceptibility, χ{sub 1}(f), the frequency dispersion of the third-order dielectric susceptibility, χ{sub 3}(f), was found to depend primarily on the α-relaxation time, and independent of temperature T and pressure P. I show this property of the frequency dispersions of χ{sub 1}(f) and χ{sub 3}(f) is the characteristic of the many-body relaxation dynamics of interacting systems which are governed solely by the intermolecular potential, and thermodynamic condition plays no role in this respect. Although linked to χ{sub 3}(f), dynamic heterogeneity is one of the parallel consequences of the many-body dynamics, and it should not be considered as the principal control parameter for the other dynamic properties of glassforming systems. Results same as χ{sub 3}(f) at elevated pressures had been obtained before by molecular dynamics simulations from the four-points correlation function and the intermediate scattering function. Naturally all properties obtained from the computer experiment, including dynamics heterogeneity, frequency dispersion, the relation between the α- and JG β-relaxation, and the breakdown of the Stokes-Einstein relation, are parallel consequences of the many-body relaxation dynamics governed by the intermolecular potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling Coupled THMC Processes and Brine Migration in Salt at High Temperatures

In this report, we present FY2014 progress by Lawrence Berkeley National Laboratory (LBNL) related to modeling of coupled thermal-hydrological-mechanical-chemical (THMC) processes in salt and their effect on brine migration at high temperatures. LBNL’s work on the modeling of coupled THMC processes in salt was initiated in FY2012, focusing on exploring and demonstrating the capabilities of an existing LBNL modeling tool (TOUGH-FLAC) for simulating temperature-driven coupled flow and geomechanical processes in salt. This work includes development related to, and implementation of, essential capabilities, as well as testing the model against relevant information and published experimental data related to the fate and transport of water. we provide more details on the FY2014 work, first presenting updated tools and improvements made to the TOUGH-FLAC simulator, and the use of this updated tool in a new model simulation of long-term THM behavior within a generic repository in a salt formation. This is followed by the description of current benchmarking and validations efforts, including the TSDE experiment. We then present the current status in the development of constitutive relationships and the dual-continuum model for brine migration. We conclude with an outlook for FY2015, which will be much focused on model validation against field experiments and on the use of the model for the design studies related to a proposed heater experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Coupling methods between finite element-based Boussinesq-type wave and particle-based free-surface flow models [Coupling methods between finite element-based Boussinesq-typewave and particle-based free-surface flow models]

We develop one-way coupling methods between a Boussinesq-type wave model based on the discontinuous Galerkin finite element method and a free-surface flow model based on a mesh-free particle method to strike a balance between accuracy and computational cost. In our proposed model, computation of the wave model in the global domain is conducted first, and the nonconstant velocity profiles in the vertical direction are reproduced by using its results. Computation of the free-surface flow is performed in a local domain included within the global domain with interface boundaries that move along the reproduced velocity field in a Lagrangian fashion. To represent the moving interfaces, we used a polygon wall boundary model for mesh-free particle methods. Verification and validation tests of our proposed model are performed, and results obtained by the model are compared with theoretical values and experimental results to show its accuracy and applicability.

97 MATHEMATICS AND COMPUTING↗

THCM Coupled Model for Hydrate-Bearing Sediments: Data Analysis and Design of New Field Experiments (Marine and Permafrost Settings)

Gas hydrates are solid compounds made of water molecules clustered around low molecular weight gas molecules such as methane, hydrogen, and carbon dioxide. Methane hydrates form under pressure (P) and temperature (T) conditions that are common in sub-permafrost layers and in deep marine sediments. Stability conditions constrain the occurrence of gas hydrates to submarine sediments and permafrost regions. The amount of technically recoverable methane trapped in gas hydrate may exceed 104tcf. Gas hydrates are a potential energy resource, can contribute to climate change, and can cause large-scale seafloor instabilities. In addition, hydrate formation can be used for CO2 sequestration (also through CO2-CH4 replacement), and efficient geological storage seals. The experimental study of hydrate bearing sediments has been hindered by the very low solubility of methane in water (lab testing), and inherent sampling difficulties associated with depressurization and thermal changes during core extraction. This situation has prompted more decisive developments in numerical modeling in order to advance the current understanding of hydrate bearing sediments, and to investigate/optimize production strategies and implications. The goals of this research has been to addresses the complex thermo-hydro-chemo-mechanical THCM coupled phenomena in hydrate-bearing sediments, using a truly coupled numerical model that incorporates sound and proven constitutive relations, satisfies fundamental conservation principles. Analytical solutions aimed at verifying the proposed code have been proposed as well. These tools will allow to better analyze available data and to further enhance the current understanding of hydrate bearing sediments in view of future field experiments and the development of production technology.

58 GEOSCIENCES↗

Coupling WRF Double-Moment 6-Class Microphysics Schemes to RRTMG Radiation Scheme in Weather Research Forecasting Model

A method to explicitly calculate the effective radius of hydrometeors in the Weather Research Forecasting (WRF) double-moment 6-class (WDM6) microphysics scheme is designed to tackle the physical inconsistency in cloud properties between the microphysics and radiation processes. At each model time step, the calculated effective radii of hydrometeors from the WDM6 scheme are linked to the Rapid Radiative Transfer Model for GCMs (RRTMG) scheme to consider the cloud effects in radiative flux calculation. This coupling effect of cloud properties between the WDM6 and RRTMG algorithms is examined for a heavy rainfall event in Korea during 25–27 July 2011, and it is compared to the results from the control simulation in which the effective radius is prescribed as a constant value. It is found that the derived radii of hydrometeors in the WDM6 scheme are generally larger than the prescribed values in the RRTMG scheme. Consequently, shortwave fluxes reaching the ground (SWDOWN) are increased over less cloudy regions, showing a better agreement with a satellite image. The overall distribution of the 24-hour accumulated rainfall is not affected but its amount is changed. A spurious rainfall peak over the Yellow Sea is alleviated, whereas the local maximum in the central part of the peninsula is increased.

Bae, Soo Ya↗

Efficient Reformulation of Solid Phase Diffusion in Electrochemical-Mechanical Coupled Models for Lithium-Ion Batteries: Effect of Intercalation Induced Stresses

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms with an appropriate discretization or approximation for the solid phase diffusion within the electrode particle. One of the major difficulties in simulating Li-ion battery models is the need for simulating solid-phase diffusion in the second radial dimension r within the particle. It increases the complexity of the model as well as the computation time/cost to a great extent. This is Particularly true for the inclusion of pressure induced diffusion inside particles experiencing volume change. A computationally efficient representation for solid-phase diffusion is discussed in this paper. The operating condition has a significant effect on the validity, accuracy, and efficiency of various approximations for the solid-phase transport governed by pressure induced diffusion. This paper introduces efficient methods for solid phase reformulation - (1) parabolic profile approach and (2) a mixed order finite difference method for approximating/representing solid-phase concentration variations within the active materials of porous electrodes for macroscopic models for lithium-ion batteries. (C) 2013 The Electrochemical Society. All rights reserved.

De, S↗

Geometric optics for a coupling model of electromagnetic and gravitational fields

The coupling between the electromagnetic and gravitational fields results in “faster than light” photons, and then the first and third laws of geometric optics are invalid in usual spacetime. By introducing an effective spacetime, we find that the wave vector can be casted into null and then it obeys the geodesic equation, the polarization vector is perpendicular to the rays, and the number of photons is conserved. That is to say, the laws of geometric optics are valid for the modified theory in the effective spacetime. We also show that the focusing theorem of light rays for the modified theory in the effective spacetime can be cast into the usual form.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implementation of a Three-Dimensional Planetary Boundary Layer Parameterization in a Coupled Modeling System and Evaluation of "Gray Zone" Simulations of a Wind-Wave Event Off the US California Coast Using Observations

A three-dimensional (3D) planetary boundary layer (PBL) parameterization was added to the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) model and the first coupled atmosphere/wave "gray zone" simulations were carried out using the 3D PBL parameterization. A case study of a significant wind-wave event from 4-8 November 2020 off the US west coast was examined, with a focus on the impacts of the event within an approximate 280x280 km central California coastal region covering the Bureau of Ocean Energy Management Morro Bay wind energy lease area. Simulations with both one-dimensional (1D; Yonsei University) and 3D PBL parameterizations were examined in this nest with ..delta..x=400 m. Two-way coupling was active, with near-surface winds feeding back to the wave model and bulk wave statistics feeding back to the atmospheric model. Both simulations compared favorably with buoy observations in capturing the timing and magnitude of wind speed, temperature, dewpoint, and significant wave height, as the front associated with the maritime weather system moved southward across coastal California. However, wind speed errors over land were larger. Time series of the vertical profile of winds below z=250 m from the simulations compared favorably with observations from the US Dept. of Energy Morro Bay Doppler lidar buoy. Differences between the coupled 1D and 3D PBL simulations were minor. The most notable difference was that the simulation using the 3D PBL parameterization had approximately 10 % weaker winds at the peak of the event causing 10 % lower significant wave heights. This difference was shown to be mostly due to differences in the vertical mixing treatment between the parameterizations, in particular the nonlocal downward mixing of higher momentum air in the 1D Yonsei University parameterization in convective conditions. Overall, this work demonstrates that the 3D PBL parameterization can be used in a coupled atmosphere/wave modeling framework with similar behaviors as traditional PBL parameterizations that do not provide horizontal turbulent variances and fluxes.

16 TIDAL AND WAVE POWER↗

Modeling Coupled Chemo-Mechanical Behavior of Randomly Oriented NMC811 Polycrystalline Li-Ion Battery Cathodes

This paper develops a three-dimensional, transient, chemo-mechanical model that predicts the performance of single secondary particle Li-ion battery cathodes. The secondary particles are composed of numerous (approximately 60) randomly oriented single-crystal primary particles. The model incorporates concentration-dependent and anisotropic material properties. As much as possible, electrochemical, transport, and structural properties for crystalline NMC811 (Li x Ni 0.8 Mn 0.1 Co 0.1 O 2 ) are taken from prior publications. Weak Van der Waals bonding between primary particles is modeled empirically using a spring analogy, which enables local primary-particle separations (disintegration) and subsequent reattachments. The model fully couples Li diffusion and the mechanical response. Results include predictions of local Li-concentrations and stresses. High stresses are found near grain boundaries, especially when the lattice orientations are greatly misaligned. Particle separations are characterized in terms of a damage parameter. The model is used to predict the effects of design and operating conditions, including charge/discharge rates, cycling scenarios, and particle sizes.

25 ENERGY STORAGE↗

Assessing Impacts of Waves on Hub-Height Winds off the U.S. West Coast Using Lidar Buoys and Coupled Modeling Approaches

Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.

17 WIND ENERGY↗

Geometric phase effects in low-energy dynamics near conical intersections: A study of the multidimensional linear vibronic coupling model

In molecular systems containing conical intersections (CIs), a nontrivial geometric phase (GP) appears in the nuclear and electronic wave functions in the adiabatic representation. We study GP effects in nuclear dynamics of an N-dimensional linear vibronic coupling (LVC) model. The main impact of GP on low-energy nuclear dynamics is reduction of population transfer between the local minima of the LVC lower energy surface. For the LVC model, we proposed an isometric coordinate transformation that confines non-adiabatic effects within a two-dimensional subsystem interacting with an N − 2 dimensional environment. Since environmental modes do not couple electronic states, all GP effects originate from nuclear dynamics within the subsystem. We explored when the GP affects nuclear dynamics of the isolated subsystem, and how the subsystem-environment interaction can interfere with GP effects. Comparing quantum dynamics with and without GP allowed us to devise simple rules to determine significance of the GP for nuclear dynamics in this model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Replica wormhole and information retrieval in the SYK model coupled to Majorana chains

Motivated by recent studies of the information paradox in (1+1)-D anti-de Sitter spacetime with a bath described by a (1+1)-D conformal field theory, we study the dynamics of second Ŕenyi entropy of the Sachdev-Ye-Kitaev (SYK) model (χ) coupled to a Majorana chain bath (ψ). The system is prepared in the thermofield double (TFD) state and then evolved by HL + HR. For small system-bath coupling, we find that the second Rényi entropy $S^(2)_{χL,χR}$ of the SYK model undergoes a first order transition during the evolution. In the sense of holographic duality, the long-time solution corresponds to a “replica wormhole”. The transition time corresponds to the Page time of a black hole coupled to a thermal bath. We further study the information scrambling and retrieval by introducing a classical control bit, which controls whether or not we add a perturbation in the SYK system. The mutual information between the bath and the control bit shows a positive jump at the Page time, indicating that the entanglement wedge of the bath includes an island in the holographic bulk.

1/N Expansion↗

Microscopic core-quasiparticle coupling model for spectroscopy of odd-mass nuclei

Predictions of the spectroscopic properties of low-lying states are critical for nuclear structure studies but are problematic for nuclei with an odd nucleon due to the interplay of the unpaired single particle with nuclear collective degrees of freedom. To predict the spectroscopic properties of odd-mass medium-heavy and heavy nuclei with a model that treats single-particle and collective degrees of freedom within the same microscopic framework. A microscopic core-quasiparticle coupling (CQC) model based on the covariant density functional theory is developed that contains the collective excitations of even-mass cores and spherical single-particle states of the odd nucleon as calculated from a quadrupole collective Hamiltonian combined with a constrained triaxial relativistic Hartree-Bogoliubov model. Additionally, predictions of the new model for excitation energies, kinematic and dynamic moments of inertia, and transition rates are shown to be in good agreement with results of low-lying spectroscopy measurements of the axially deformed odd-proton nucleus 159 Tb and the odd-neutron nucleus 157 Gd. A microscopic CQC model based on covariant density functional theory is developed for odd-mass nuclei and shown to give predictions that agree with measurements of two medium-heavy nuclei. Future studies with additional nuclei are planned.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Reheating in non-minimal derivative coupling model

We consider a model with non-minimal derivative coupling of inflaton to gravity. The reheating process during rapid oscillation of the inflaton is studied and the reheating temperature is obtained. Behaviors of the inflaton and produced radiation in this era are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

High Order Accurate Finite Difference Modeling of Seismo-Acoustic Wave Propagation in a Moving Atmosphere and a Heterogeneous Earth Model Coupled Across a Realistic Topography

Here, we develop a numerical method for simultaneously simulating acoustic waves in a realistic moving atmosphere and seismic waves in a heterogeneous earth model, where the motions are coupled across a realistic topography. We model acoustic wave propagation by solving the linearized Euler equations of compressible fluid mechanics. The seismic waves are modeled by the elastic wave equation in a heterogeneous anisotropic material. The motion is coupled by imposing continuity of normal velocity and normal stresses across the topographic interface. Realistic topography is resolved on a curvilinear grid that follows the interface. The governing equations are discretized using high order accurate finite difference methods that satisfy the principle of summation by parts. We apply the energy method to derive the discrete interface conditions and to show that the coupled discretization is stable. The implementation is verified by numerical experiments, and we demonstrate a simulation of coupled wave propagation in a windy atmosphere and a realistic earth model with non-planar topography.

58 GEOSCIENCES↗

Modeling Coupled Durability and Performance in Polymer-Electrolyte Fuel Cells: Membrane Effects

During operation, proton-exchange-membrane fuel cells (PEMFCs) are subjected to mechanical and chemical stressors that contribute to membrane degradation, performance loss, and eventual failure. Together, synergistic effects between mechanical and chemical degradation mechanisms lead to accelerated degradation. A physics-based model is developed to understand the synergistic effects of chemical and mechanical degradation and the coupled nature of performance and durability in PEMFCs. The model includes pinhole existence and growth in the membrane, which increases crossover of reactant gases as well as subsequent formation of chemical degradation agents that impact both transport and mechanical properties of the membrane. The fuel-cell model is fully coupled with a mechanical model to determine the stresses on the membrane and subsequent growth of pinholes during transient operation. Simulation results demonstrate pinhole growth under relative-humidity cycling and the resultant increased gas-crossover fluxes and decrease in polarization performance. Furthermore, the model results highlight nonlinearities and the importance of coupling mechanical and chemical degradation models in order to explain membrane degradation under various cycles, and serves as a foundation for examining coupled durability and performance.

25 ENERGY STORAGE↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗