Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 937 records · Page 52

Analysis of vibronic coupling in a 4f molecular magnet with FIRMS

Vibronic coupling, the interaction between molecular vibrations and electronic states, is a fundamental effect that profoundly affects chemical processes. In the case of molecular magnetic materials, vibronic, or spin-phonon, coupling leads to magnetic relaxation, which equates to loss of magnetic memory and loss of phase coherence in molecular magnets and qubits, respectively. The study of vibronic coupling is challenging, and most experimental evidence is indirect. Here we employ far-infrared magnetospectroscopy to directly probe vibronic transitions in [Yb(trensal)] (where H 3 trensal = 2,2,2-tris(salicylideneimino)trimethylamine). We find intense signals near electronic states, which we show arise due to an “envelope effect” in the vibronic coupling Hamiltonian, which we calculate fully ab initio to simulate the spectra. We subsequently show that vibronic coupling is strongest for vibrational modes that simultaneously distort the first coordination sphere and break the C 3 symmetry of the molecule. With this knowledge, vibrational modes could be identified and engineered to shift their energy towards or away from particular electronic states to alter their impact. Hence, these findings provide new insights towards developing general guidelines for the control of vibronic coupling in molecules.

36 MATERIALS SCIENCE↗

Urbanization-induced land and aerosol impacts on storm propagation and hail characteristics

Changes in land surface and aerosol characteristics from urbanization can affect dynamic and microphysical properties of severe storms, thus affecting hazardous weather events resulting from such storms such as hail and tornado. We examine the joint and individual effects of urban land and anthropogenic aerosols of Kansas City on a severe convective storm observed during the 2015 Plains Elevated Convection At Night (PECAN) field campaign, focusing on storm evolution, convective intensity, and hail characteristics. The simulations are carried out at the cloud-resolving scale (1 km) using a version of WRF-Chem in which the spectral-bin microphysics (SBM) is coupled with the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC). It is found that the urban land effect of Kansas City initiated a much stronger convective cell and the storm got further intensified when interacting with stronger turbulence induced by the urban land. The urban land effect also changed the storm path by diverting the storm toward the city mainly resulting from enhanced urban land-induced convergence in the urban area and around the urban-rural boundaries. The joint effect of urban land and anthropogenic aerosols enhances occurrences of both severe hail and significant severe hail by ~ 20% by enhancing hail formation and growth from riming. Overall the urban land effect on convective intensity and hail is relatively larger than the anthropogenic aerosol effect, but the joint effect is much notable than either of the individual effects, emphasizing the importance to consider both effects in evaluating urbanization effects.

Lin, Yun↗

Effects of Coulomb coupling on stopping power and a link to macroscopic transport

Molecular dynamics simulations are used to assess the influence of Coulomb coupling on the energy evolution of charged projectiles in the classical one-component plasma. The average projectile kinetic energy is found to decrease linearly with time when ν α /ω p ≲ 10 –2 , where να is the Coulomb collision frequency between the projectile and the medium, and ω p is the plasma frequency. Stopping power is obtained from the slope of this curve. In comparison to the weakly coupled limit, strong Coulomb coupling causes the magnitude of the dimensionless stopping power, ( a/k B T ) dE / dx , to increase, the Bragg peak to shift to several times the plasma thermal speed, and for the stopping power curve to broaden substantially. The rate of change of the total projectile kinetic energy averaged over many independent simulations is shown to consist of two measurable components: a component associated with a one-dimensional friction force and a thermal energy exchange rate. In the limit of a slow and massive projectile, these can be related to the macroscopic transport rates of self-diffusion and temperature relaxation in the plasma. Simulation results are compared with available theoretical models for stopping power, self-diffusion, and temperature relaxation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial Neutronics and Thermal-Hydraulic Coupling for Spent Nuclear Fuel Canister

This report documents work performed supporting the US Department of Energy (DOE) Nuclear Energy Spent Fuel and Waste Disposition, Spent Fuel and Waste Science and Technology, under work breakdown structure element 1.08.01.03.05, “Direct Disposal of Dual Purpose Canisters.” In particular, this appendix fulfills the M3 milestone, M3SF-19OR010305015, “Multiphysics criticality consequence analysis capability development status report,” as Revision 1 to M3SF-19OR010305016, “Initial neutronic and thermal hydraulic coupling for waste package” within work package SF-19OR01030501, “Direct Disposal of Dual Purpose Canisters–ORNL.” This report presents the initial development status of a multiphysics criticality consequence simulation framework, Terrenus. Terrenus currently couples two physics codes: the Monte Carlo radiation transport code Shift and the thermal-hydraulics code COBRA-SFS. The coupling capability has been demonstrated using a simplified 3 x 3 fuel pin lattice. Terrenus will be developed further to include (1) general geometry package support, allowing modeling of a spent nuclear fuel (SNF) cask with all structural details, (2) a depletion solver to determine the change in nuclide composition at the end of each time step, and (3) a mechanics code to determine any structural impact due to a criticality event. Future research and development works include identifying/developing a modern two-phase thermal-hydraulics code and developing an approach to solve fast neutronic transients. The goal of this multiphysics framework is to determine the feasibility of direct disposal of currently loaded SNF canisters by including or excluding a criticality event from a repository performance analysis framework in terms of consequences. The revision 1 to this report includes extending Terrenus to a 17 x 17 pressurized water reactor (PWR) fuel assembly. We also investigated the phenomena necessary to quench a critical reaction, and investigated the water level necessary to achieve criticality in a typical canister. Finally, we describe the initial implementation of a critical power search capability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Report of the Atmosphere to Electrons Mesoscale-to-Microscale Coupling Project (FY2019)

The overall goal of the Mesoscale-to-Microscale Coupling (MMC) project is to improve coupling between mesoscale and microscale simulations via improved guidance and new strategies for setting up simulations and the development of new tools that can be used across the community. Including the mesoscale forcing is critical to modeling the full energy transfer across scales in the atmosphere. The project-specific objectives include: Apply rigorous verification and validation (V&V) techniques to the new modeling tools that are developed as part of the project to ensure the accuracy of our codes and results and develop estimates of the relative uncertainty. Improve computational performance of the coupled MMC models through the development of methods that can be used to reduce turbulence spin-up time and hence the size of computational domains. Improve representation of the surface layer in microscale models to enhance simulations of hub-height wind speed. Develop guidance for the community describing the best ways to couple mesoscale and microscale models, including specific spatial scales at which the handoff to the microscale model should occur. Prepare documentation and a suite of software tools that can be used across the community. Transition MMC research to the offshore environment.

17 WIND ENERGY↗

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗

BlueCRAB Domain Overlapping Coupling: Theory & Verification

Coupling low- and high-fidelity codes is a useful way to model complex engineering systems. Thanks to the high-fidelity code, complex phenomena can be resolved in areas of the system where this is required, yet the efficiency of the lower-fidelity code is still retained in modeling the rest of the system. This document details the theory and implementation of coupling two different thermal-hydraulic codes: the system thermal-hydraulics (STH) code System Analysis Module (SAM) and the coarse-mesh computational fluid dynamics (CFD) code Pronghorn. Both applications are included in the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite and are based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. The domain overlapping (DO) coupling approach was adopted, as it offers proven advantages over more conventional domain decomposition methods. In the DO coupling, SAM provides Pronghorn with boundary conditions that depend on the system-level simulation of the entire plant. In return, the overlapping coupled SAM components–termed “surrogate components”–receive friction factors and source terms computed online based on the Pronghorn simulation. The framework is designed to be generic and enable coupling regardless of geometry and the number of inlet/outlet boundaries in the DO coupled domain. The developed method leads to consistent pressure drops, enthalpies, and scalar concentrations when comparing the coupled SAM and Pronghorn simulations. This document presents the DO coupling approach, along with two verification and two demonstration cases. The proposed problems explore different physical aspects relevant to nuclear reactor analysis, including buoyancy-driven flows, complex flow patterns, and multiple inlets and outlets. Periodically, new versions of this “BlueCRAB Domain Overlapping Coupling: Theory & Verification” report will be issued to reflect future developments and verification tests.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Higher-form symmetry and chiral transport in real-time Abelian lattice gauge theory

We study classical lattice simulations of theories of electrodynamics coupled to charged matter at finite temperature, interpreting them using the higher-form symmetry formulation of magnetohydrodynamics (MHD). We compute transport coefficients using classical Kubo formulas on the lattice and show that the properties of the simulated plasma are in complete agreement with the predictions from effective field theories. In particular, the higher-form formulation allows us to understand from hydrodynamic considerations the relaxation rate of axial charge in the chiral plasma observed in previous simulations. A key point is that the resistivity of the plasma – defined in terms of Kubo formulas for the electric field in the 1-form formulation of MHD – remains a well-defined and predictive quantity at strong electromagnetic coupling. However, the Kubo formulas used to define the conventional conductivity vanish at low frequencies due to electrodynamic fluctuations, and thus the concept of the conductivity of a gauged electric current must be interpreted with care.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Global variability in atmospheric new particle formation mechanisms

A key challenge in aerosol pollution studies and climate change assessment is to understand how atmospheric aerosol particles are initially formed. Although new particle formation (NPF) mechanisms have been described at specific sites, in most regions, such mechanisms remain uncertain to a large extent because of the limited ability of atmospheric models to simulate critical NPF processes. Here we synthesize molecular-level experiments to develop comprehensive representations of 11 NPF mechanisms and the complex chemical transformation of precursor gases in a fully coupled global climate model. Combined simulations and observations show that the dominant NPF mechanisms are distinct worldwide and vary with region and altitude. Previously neglected or underrepresented mechanisms involving organics, amines, iodine oxoacids and HNO 3 probably dominate NPF in most regions with high concentrations of aerosols or large aerosol radiative forcing; such regions include oceanic and human-polluted continental boundary layers, as well as the upper troposphere over rainforests and Asian monsoon regions. These underrepresented mechanisms also play notable roles in other areas, such as the upper troposphere of the Pacific and Atlantic oceans. Accordingly, NPF accounts for different fractions (10–80%) of the nuclei on which cloud forms at 0.5% supersaturation over various regions in the lower troposphere. The comprehensive simulation of global NPF mechanisms can help improve estimation and source attribution of the climate effects of aerosols.

54 ENVIRONMENTAL SCIENCES↗

Modeling DOE Standard Canister Configurations with Updated Surface Chemistry

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation rate as the total dose applied in increased. A previous study created a small-scale chemical model was built to replicate a mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. The previous iteration of the model utilized step function for its G-values for the hydrogen generation over a 200-year period. This model makes an update to the surface chemistry to account for theorized surface chemistry reactions allowing for oxygen to remain bounded to the oxyhydroxide layer. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. Three canister configurations are modeled – the base 18-inch, 15-foot DOR standard canister with 30 ATR fuel elements, an 18-inch, 10-foot DOE standard canister with 32 ATR fuel elements and a 24-inch, 10-foot DOE standard canister with 40 ATR fuel elements. In previous reports, the primary sensitivity of the canister conditions was identified as the decay heat of the fuel and the dried condition of the fuel. Only these parameters are studied in the report. For the nominal case with dried fuel, the hydrogen generation is only about 2%, so it is even less than the flammability limit if it were exposed to oxygen. For the densely packed 18-inch case the total hydrogen concentration is only around 4% with the nominal decay heat. For the densely packed 24-inch case the total hydrogen concentration is only 2.5%, roughly 20% higher than the original packing design, and still under flammability conditions if exposed to oxygen.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Guidelines for Publicly Archiving Terrestrial Model Data to Enhance Usability, Intercomparison, and Synthesis

Scientific communities are increasingly publishing data to evaluate, accredit, and build on published research. However, guidelines for curating data for publication are sparse for model-related research, limiting the usability of archived simulation data. In particular, there are no established guidelines for archiving data related to terrestrial models that simulate land processes and their coupled interactions with climate. Terrestrial modelers have a unique set of challenges when publishing data due to the diversity of scientific domains, research questions, and the types and scales of simulations. Researchers in the U.S. Department of Energy’s (DOE) projects use a variety of multiscale models to advance robust predictions of terrestrial and subsurface ecosystem processes. Here, we synthesize archiving needs for data associated with different DOE models, and provide guidelines for publishing terrestrial model data components following FAIR (Findable, Accessible, Interoperable, Reusable) principles. The guidelines recommend archiving model inputs and testing data used in final simulation runs along with associated codes, workflow scripts, and metadata in public repositories. Researchers should consider archiving model outputs if they are within the storage limits of the repository. We also provide considerations for how to bundle files into different data publications with citable digital object identifiers. Finally, we identify repository features and tools that would enable storage and reuse of model data. Given the diversity of DOE terrestrial models, these guidelines are transferable to other model types and will enable efficient reuse of simulation data for purposes such as model intercomparisons, initialization, benchmarking, synthesis, and comparisons with field observations.

58 GEOSCIENCES↗

First-principles based plasma profile predictions for optimized stellarators

In the present Letter, first-of-its-kind computer simulations predicting plasma profiles for modern optimized stellarators—while self-consistently retaining neoclassical transport, turbulent transport with 3D effects, and external physical sources—are presented. These simulations exploit a newly developed coupling framework involving the global gyrokinetic turbulence code GENE-3D, the neoclassical transport code KNOSOS, and the 1D transport solver TANGO. This framework is used to analyze the recently observed degradation of energy confinement in electron-heated plasmas in the Wendelstein 7-X stellarator, where the central ion temperature was 'clamped' to $T_\textrm i \approx 1.5$ keV regardless of the external heating power. By performing first-principles based simulations, we provide key evidence to understand this effect, namely the inefficient thermal coupling between electrons and ions in a turbulence-dominated regime, which is exacerbated by the large $T_\textrm e/T_\textrm i$ ratios, and show that a more efficient ion heat source, such as direct ion heating, will increase the on-axis ion temperature. This work paves the way towards the use of high-fidelity models for the development of the next generation of stellarators, in which neoclassical and turbulent transport are optimized simultaneously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Systematic Benchmarking of Climate Models: Methodologies, Applications, and New Directions

As climate models become increasingly complex, there is a growing need to comprehensively and systematically assess model performance with respect to observations. Given the increasing number and diversity of climate model simulations in use, the community has moved beyond simple model intercomparison and toward developing methods capable of benchmarking a large number of simulations against a suite of climate metrics. Here, we present a detailed review of evaluation and benchmarking methods and approaches developed in the last decade, focusing primarily on scientific implications for Coupled Model Intercomparison Project (CMIP) simulations and CMIP6 results that contributed to the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report (AR6). Based on this review, we explain the resulting contemporary philosophy of model benchmarking, and provide clear distinctions and definitions of the terms model verification, process validation, evaluation, and benchmarking. While significant progress has been made in model development based on systematic evaluation and benchmarking efforts, some climate system biases still remain. The development of open‐source community software packages has played a fundamental role in identifying areas of significant model improvement and bias reduction. We review the key features of several software packages that have been commonly used over the past decade to evaluate and benchmark global and regional climate models. Additionally, we discuss best practices for the selection of evaluation and benchmarking metrics and for interpreting the obtained results, the importance of selecting suitable sources of reference data and accurate uncertainty quantification.

Environmental sciences↗

Efficient generation and amplification of intense vortex and vector laser pulses via strongly-coupled stimulated Brillouin scattering in plasmas

The past decade has seen tremendous progress in the production and utilization of vortex and vector laser pulses. Although both are considered as structured light beams, the vortex lasers have helical phase fronts and phase singularities, while the vector lasers have spatially variable polarization states and polarization singularities. In contrast to the vortex pulses that carry orbital angular momentum (OAM), the vector laser pulses have a complex spin angular momentum (SAM) and OAM coupling. Despite many potential applications enabled by such pulses, the generation of high-power/-intensity vortex and vector beams remains challenging. Here, we demonstrate using theory and three-dimensional simulations that the strongly-coupled stimulated Brillouin scattering (SC-SBS) process in plasmas can be used as a promising amplification technique with up to 65% energy transfer efficiency from the pump beam to the seed beam for both vortex and vector pulses. We also show that SC-SBS is strongly polarization-dependent in plasmas, enabling an all-optical polarization control of the amplified seed beam. Additionally, the interaction of such structured lasers with plasmas leads to various angular momentum couplings and decouplings that produce intense new light structures with controllable OAM and SAM. This scheme paves the way for novel optical devices such as plasma-based amplifiers and light field manipulators.

Physics↗

Representation of atmosphere-induced heterogeneity in land–atmosphere interactions in E3SM–MMFv2

In the Energy Exascale Earth System Model (E3SM) Multi-scale Modeling Framework (MMF), where parameterizations of convection and turbulence are replaced by a 2-D cloud-resolving model (CRM), there are multiple options to represent land–atmosphere interactions. Here, we propose three different coupling strategies, namely the (1) coupling of a single land surface model to the global grid (MMF), (2) coupling a single land copy directly to the embedded CRM (SFLX2CRM), and (3) coupling a single copy of land model to each column of the CRM grid (MAML). In the MAML (Multi-Atmosphere Multi-Land) framework, a land model is coupled to CRM at the CRM-grid scale by coupling an individual copy of a land model to each CRM grid. Therefore, we can represent intra-CRM heterogeneity in the land–atmosphere interaction processes. There are 5-year global simulations run using these three coupling strategies, and we find some regional differences but overall small changes with respect to whether a land model is coupled to CRM or a global atmosphere. In MAML, the spatial heterogeneity within CRM induces stronger turbulence, which leads to the changes in soil moisture, surface heat fluxes, and precipitation. However, the differences in the MAML from the other two cases are rather weak, suggesting that the impact of using MAML does not justify the increase in cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Intercalation of Lithium into Graphite: Insights from First-Principles Simulations

Understanding ion intercalation at electrode–electrolyte interfaces is key to the development of energy storage and water desalination. Here, we investigate Li + kinetics at a prototypical interface between graphite anodes and an organic electrolyte, and we elucidate key factors that determine ion transport, using first-principles methodology coupling ab initio molecular dynamics simulations with a solvation model. We show that surface chemical composition significantly influences the kinetics of ion intercalation from the liquid into graphite. We find that this is partly related to the ion desolvation process, which varies notably for different graphite surface chemical terminations. In addition, interfacial polarization is found to play an important role in determining energy barriers for ion transfer. We also discuss the impact of electrode potentials, which is often neglected in conventional first-principles calculations despite being a key factor in device configurations. Our study provides insights into the coupling of electronic and ionic effects of interfacial chemistry on ion transport at complex electrode–electrolyte interfaces.

25 ENERGY STORAGE↗

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗

General overview of multiphysics modeling and simulation at CEA/DES/IRESNE supporting safety, operating and development of nuclear reactors and facilities

This paper presents a general overview of the innovative multiphysics research/development activities carried on by the CEA/DES/IRESNE institute. Eight main different multiphysics fields addressed by IRESNE have been listed (among others: fuel fabrication, behavior of fuels, of reactors, of corium, radionuclide transfer to environment,...) in relation with the involved coupled physics: neutron transport, fluid mechanics, electro-magnetism, heat transfer, solid mechanics and (physico-)chemistry. The major mission of IRESNE institute is to study current and future nuclear systems integrated into a low carbon energy system. To fulfill this mission, IRESNE's activities rely on the development and the implementation of multiphysics modeling, advanced coupling scheme, data transfer techniques and dedicated multiphysics experimentations. Indeed, such multiphysics simulation approach of nuclear systems is today required to deal with the complex intrinsic features of studied phenomena. This approach aims at taking into account physical phenomena interactions, at each interest modeling scale, in order to reach an accurate and predictive behavior representation of a system component or of the whole system. The final objective is to contribute to the development of digital reactor, core and fuel at CEA, by progressing on applied mathematics, high-performance computing and physical mechanism understanding.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗