Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled models”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be to infer species separation, thermal decoupling, flows or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to underlying dependence on the shape and time evolution of the temperature and density proles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second order analytic expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second order expression shows significant improvement over lower order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling both first and second order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with <10% error the second order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Thermal-Hydraulics Modeling and Simulations of Hot Pool Using the SAS-CFD Coupled Code

The main goal of this activity is to test the dynamic coupling of the SAS4A/SASSYS-1 (SAS) and CFD models, using a recently patched version of the SAS code intended to address an undocumented limitation that hindered the Versatile Test Reactor (VTR) simulation efforts in FY21. As described in previous VTR calculation reports, the undocumented limitation in SAS v5.4 does not allow the user to activate the CFD coupling option during restart calculations. Since the analysts were unaware of this limitation, prior SAS-CFD simulation results for the protected station blackout (PSBO) transient were erroneous. Root-cause analysis was performed to determine the cause of this undocumented limitation in SAS v5.4, the SAS software was updated in a new patch, and the SAS-CFD simulations were repeated with this patched software. The results of the SAS-CFD simulations documented in this report show that the software patch does address the cited issue, and that the patched software indeed supports the activation of the CFD coupling model in restart calculations. The SAS development team will determine the schedule for implementing the patch in an official software release. This report documents updated SAS-CFD simulations of the PSBO transient response in the VTR. The hot pool is modeled with the CFD code STAR-CCM+, which is coupled at the flow boundaries to the SAS model of the primary heat transport system. SAS computes the mass flow rate and temperature at each core subassembly outlet, the thermal insulation cavity bypass, and the IHX inlet windows. CFD in turn computes the absolute pressure and temperature at each of these boundaries. The SAS code will ignore the temperature data at flow boundaries where flow is directed into the hot pool, i.e., at the core subassembly outlets unless flow reversal occurs. Similarly, CFD will ignore temperature data at boundaries where the flow is directed out of the hot pool, i.e., at the IHX inlets except under flow reversal. The focus of this work is to ensure that the SAS software patch addresses the undocumented limitation described in prior VTR calculation reports, rather than the accurate assessment of thermal stratification in the VTR during protected transients. This motivates the development of a new, simplified CFD model with a coarser mesh to accelerate the testing process. The updated model, and simplifying assumptions, are documented in this report. In future work, the thermal stratification assessment should be performed in more detail. The simplified CFD model can be improved by performing grid convergence studies sensitivity studies of turbulence parameters (e.g., Prandtl number, turbulence production and dissipation parameters) on temperature distributions and thermal stratification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A second order yield-temperature relation for accurate inference of burn-averaged quantities in multi-species plasmas

Measured yields and ion temperatures inferred from the fusion product energy spectra can be used as metrics for the performance of an ICF implosion. This can be used to infer species separation, thermal decoupling, flows, or other effects that can cause the inferred ion temperatures to deviate from the true underlying thermal temperature and the yield ratio to deviate from the expected value. Direct inference of the impact of these effects on observed temperatures and yields can be difficult to uncover due to the underlying dependence on the shape and time evolution of the temperature and density profiles of the fusing plasma. Due to differences in the temperature dependence of the reactivities, different fusion products are emitted from different regions and times within the implosion. In order to properly account for this, a second-order analytical expression relating the apparent temperatures and yield ratios is developed. This expression can be coupled to models of yield and/or temperature altering effects to infer their burn-averaged impact on an implosion. The second-order expression shows significant improvement over lower-order expressions in synthetic data studies. Demonstrations of its applications to synthetic data coupled with models of ion thermal decoupling and radial flows are presented. In the case of thermal decoupling, both first and second-order expressions show reasonable levels of accuracy. To consistently infer the amplitude of radial flow with a <10% error, the second-order equation is required.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Toward a 3d Ising model with a weakly-coupled string theory dual

It has long been expected that the 3d Ising model can be thought of as a string theory, where one interprets the domain walls that separate up spins from down spins as two-dimensional string worldsheets. The usual Ising Hamiltonian measures the area of these domain walls. This theory has string coupling of unit magnitude. We add new local terms to the Ising Hamiltonian that further weight each spin configuration by a factor depending on the genus of the corresponding domain wall, resulting in a new 3d Ising model that has a tunable bare string coupling g_s g s . We use a combination of analytical and numerical methods to analyze the phase structure of this model as g_s g s is varied. We study statistical properties of the topology of worldsheets and discuss the prospects of using this new deformation at weak string coupling to find a worldsheet description of the 3d Ising transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Design and analysis of CXL performance models for tightly-coupled heterogeneous computing

Truly heterogeneous systems enable partitioned workloads to be mapped to the hardware that nets the best performance. However, current practice requires that inter-device communication between different vendors' hardware use host memory as an intermediary step. To date, there are no widely adopted solutions that allow accelerators to directly transfer data. A new cache-coherent protocol, CXL, aims to facilitate easier, fine-grained sharing between accelerators. In this work we analyze existing methods for designing heterogeneous applications that target GPUs and FPGAs working collaboratively, followed by an exploration to show the benefits of a CXL-enabled system. Specifically, we develop a test application that utilizes both an NVIDIA P100 GPU and a Xilinx U250 FPGA to show current communication limitations. From this application, we capture overall execution time and throughput measurements on the FPGA and GPU. We use these measurements as inputs to novel CXL performance models to show that using CXL caching instead of host memory results in a 1.31X speedup, while a more tightly-coupled pipelined implementation using CXL-enabled hardware would result in a speedup of 1.45X.

Cabrera, Anthony↗
DOE OSTI2020

Simulating the Trajectory and Biomass Growth of Free-Floating Macroalgal Cultivation Platforms along the U.S. West Coast

Trajectory tracking and macroalgal growth models were coupled to support a novel macroalgae-harvesting concept known as the Nautical Off-shore Macroalgal Autonomous Device (NOMAD). The NOMAD consists of 5 km long carbon-fiber longlines that are seeded and free float southward along the U.S. West Coast for approximately 3 months before harvesting off the California coast, taking advantage of favorable environmental conditions. The trajectory and macroalgal growth models were applied to answer planning questions pertinent to the techno-economic analysis such as identifying the preferred release location, approximate pathway, timing until harvest, and estimated growth. Trajectories were determined with the General NOAA Operational Modeling Environment (GNOME) model, using 11 years of current and wind data, determining probabilities by running nearly 40,000 Monte Carlo simulations varying the start time and location. An accompanying macroalgal growth model was used to estimate the growth of macroalgae based on the trajectory tracks and environmental forcing products, including light, temperature and nutrients. Model results show that NOMAD lines transit south in the months of April to September due to seasonal currents, taking approximately 3 months to reach Southern California. During transit, NOMAD lines are dispersed but typically avoid beaching or passing through marine sanctuaries. NOMAD lines can yield up to 30 kg wet weight per meter of cultivation line.

09 BIOMASS FUELS↗
DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗
DOE OSTI2024

Kicking the can down the road: understanding the effects of delaying the deployment of stratospheric aerosol injection

Abstract Climate change is a prevalent threat, and it is unlikely that current mitigation efforts will be enough to avoid unwanted impacts. One potential option to reduce climate change impacts is the use of stratospheric aerosol injection (SAI). Even if SAI is ultimately deployed, it might be initiated only after some temperature target is exceeded. The consequences of such a delay are assessed herein. This study compares two cases, with the same target global mean temperature of ∼1.5° C above preindustrial, but start dates of 2035 or a ‘delayed’ start in 2045. We make use of simulations in the Community Earth System Model version 2 with the Whole Atmosphere Coupled Chemistry Model version 6 (CESM2-WACCM6), using SAI under the SSP2-4.5 emissions pathway. We find that delaying the start of deployment (relative to the target temperature) necessitates lower net radiative forcing (−30%) and thus larger sulfur dioxide injection rates (+20%), even after surface temperatures converge, to compensate for the extra energy absorbed by the Earth system. Southern hemisphere ozone is higher from 2035 to 2050 in the delayed start scenario, but converges to the same value later in the century. However, many of the surface climate differences between the 2035 and 2045 start simulations appear to be small during the 10–25 years following the delayed SAI start, although longer simulations would be needed to assess any longer-term impacts in this model. In addition, irreversibilities and tipping points that might be triggered during the period of increased warming may not be adequately represented in the model but could change this conclusion in the real world.

Brody, Ezra (ORCID:000900030008681X)↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Multifidelity Ensemble Kalman Filtering Using Surrogate Models Defined by Theory-Guided Autoencoders

Data assimilation is a Bayesian inference process that obtains an enhanced understanding of a physical system of interest by fusing information from an inexact physics-based model, and from noisy sparse observations of reality. The multifidelity ensemble Kalman filter (MFEnKF) recently developed by the authors combines a full-order physical model and a hierarchy of reduced order surrogate models in order to increase the computational efficiency of data assimilation. The standard MFEnKF uses linear couplings between models, and is statistically optimal in case of Gaussian probability densities. This work extends the MFEnKF into to make use of a broader class of surrogate model such as those based on machine learning methods such as autoencoders non-linear couplings in between the model hierarchies. We identify the right-invertibility property for autoencoders as being a key predictor of success in the forecasting power of autoencoder-based reduced order models. We propose a methodology that allows us to construct reduced order surrogate models that are more accurate than the ones obtained via conventional linear methods. Numerical experiments with the canonical Lorenz'96 model illustrate that nonlinear surrogates perform better than linear projection-based ones in the context of multifidelity ensemble Kalman filtering. We additionality show a large-scale proof-of-concept result with the quasi-geostrophic equations, showing the competitiveness of the method with a traditional reduced order model-based MFEnKF.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A new approach to model geomaterials with heterogeneous properties in thermo-hydro-mechanical coupled problems

The main objective of this article is to present a new approach to model coupled thermo-hydro-mechanical problems considering geomaterials with heterogeneous properties. This approach has been implemented in the software CODE_BRIGHT and it provides the possibility of considering geomaterials with a spatially correlated heterogeneous field of porosity, following a normal distribution. This spatial correlation can be isotropic or anisotropic. An important feature of this approach is that material properties such as intrinsic permeability, thermal conductivity, diffusivity, retention curve, elastic modulus or cohesion are defined as a function of porosity and, thus, they become heterogeneous with spatial correlation and, eventually, anisotropic. A validation exercise and other basic numerical examples have been carried out to illustrate the possibilities of the proposed approach. The results, which have been compared with a homogeneous case, show that considering heterogeneous fields can be relevant in different modelling problems, especially coupled thermo-hydro-mechanical problems.

42 ENGINEERING↗
DOE OSTI2020

Understanding Microstructure Variability in Vapor-Deposited Energetic Materials by Using Phase-Field Methods

Critical components, such as detonators, in Sandia's stockpile contain heterogeneous materials whose performance and reliability depend on accurate, predictive models of coupled, complex phenomena to predict their synthesis, processing, and operation. Ongoing research in energetic materials has shown that microstructural properties, such as density, pore-size, morphology, and specific surface area are correlated to their initiation threshold and detonation behavior. However, experiments to study these specific characteristics of energetic materials are challenging and time consuming. Therefore, in this work, we turn to mesoscale modeling methods that may be capable of reproducing some observed phenomena to refine and predict outcomes beforehand. Even so, we have no physics-based modeling capability to predict how the microstructure of an energetic material will evolve over near- and long-term time scales. Thus, the goal of this work is to (i) identify any knowledge gaps in how the underlying microstructure forms and evolves during the synthesis process, and (ii) develop and test a mesoscale phase-field model for vapor deposition to capture critical mechanisms of microstructure formation, evolution, and variability in vapor-deposited energetic materials, such as processing conditions, material properties, and substrate interactions. Based on this work, the phase-field method is shown to be a valuable tool for developing the necessary models containing coupled, complex phenomena to investigate and understand the synthesis and processing of energetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of Surface Turbulence Flux Parameterizations on the MJO: The Role of Ocean Surface Waves

This study investigates the sensitivity of the Madden–Julian oscillation (MJO) to changes to the bulk flux parameterization and the role of ocean surface waves in air–sea coupling using a fully coupled ocean–atmosphere–wave model. The atmospheric and ocean model components of the Energy Exascale Earth System Model (E3SM) are coupled to a spectral wave model, WAVEWATCH III (WW3). Two experiments with wind speed–dependent bulk algorithms (NCAR and COARE3.0a) and one experiment with wave-state-dependent flux (COR3.0a-WAV) were conducted. We modify COARE3.0a to include surface roughness calculated within WW3 and also account for the buffering effect of waves on the relative difference between air-side and ocean-side momentum flux. Differences in surface fluxes, primarily caused by discrepancies in drag coefficients, result in significant differences in MJO’s properties. While COARE3.0a has better convection–circulation coupling than NCAR, it exhibits anomalous MJO convection east of the date line. The wave-state-dependent flux (COR3.0-WAV) improves the MJO representation over the default COARE3.0 algorithm. Strong easterlies over the Pacific Ocean in COARE3.0a enhance the latent heat flux (LHFLX). This is responsible for the anomalous MJO propagation after the date line. In COR3.0a-WAV, waves reduce the anomalous easterlies, leading to a decrease in LHFLX and MJO dissipation after the date line. These findings highlight the role of surface fluxes in MJO simulation fidelity. Most importantly, we show that the proper treatment of wave-induced effects in bulk flux parameterization improves the simulation of coupled climate variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Investigating the Impact of Irrigation on Malaria Vector Larval Habitats and Transmission Using a Hydrology‐Based Model

Abstract A combination of accelerated population growth and severe droughts has created pressure on food security and driven the development of irrigation schemes across sub‐Saharan Africa. Irrigation has been associated with increased malaria risk, but risk prediction remains difficult due to the heterogeneity of irrigation and the environment. While investigating transmission dynamics is helpful, malaria models cannot be applied directly in irrigated regions as they typically rely only on rainfall as a source of water to quantify larval habitats. By coupling a hydrologic model with an agent‐based malaria model for a sugarcane plantation site in Arjo, Ethiopia, we demonstrated how incorporating hydrologic processes to estimate larval habitats can affect malaria transmission. Using the coupled model, we then examined the impact of an existing irrigation scheme on malaria transmission dynamics. The inclusion of hydrologic processes increased the variability of larval habitat area by around two‐fold and resulted in reduction in malaria transmission by 60%. In addition, irrigation increased all habitat types in the dry season by up to 7.4 times. It converted temporary and semi‐permanent habitats to permanent habitats during the rainy season, which grew by about 24%. Consequently, malaria transmission was sustained all‐year round and intensified during the main transmission season, with the peak shifted forward by around 1 month. Lastly, we evaluated the spatiotemporal distribution of adult vectors under the effect of irrigation by resolving habitat heterogeneity. These findings could help larval source management by identifying transmission hotspots and prioritizing resources for malaria elimination planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure–property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius–Duhem inequality. The kinetic framework employs a stress–strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model’s results showing good correlations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predictive Capability Maturity Model Demonstration for Cylindrical Cavity Coupling Using Gemma in the Next Generation Workflow

The predictive capability maturity model (PCMM) uses the expert elicitation process to generate credibility evidence for a particular analysis. To ensure Gemma has the capability to efficiently produce this credibility evidence, next generation workflows (NGW) are created for the solution verification, calibration/validation, and input uncertainty quantification portions of the PCMM assessment. These workflows are then used on the Higgins cylinder problem, which is representative of applications involving external-to-internal electromagnetic field coupling through a slot. The uncertainties calculated using these workflows are then used to calculate the validation comparison error and the validation uncertainty for the model following the American Society of Mechanical Engineers (ASME) verification and validation (V&V) 20 standard. These workflows will enable analysts to iterate each element of PCMM more efficiently than if completed without using a NGW workflow. An example of this iterative process is shown in Section 7.2.

42 ENGINEERING↗