Reduced Order Modeling of Coupled Radiation-Conduction Problems Using Learned Embedding Space Dynamics
Explore the source record for details and available documents.
Hajimirza, Shima [Stevens Institute of Technology]↗
SEARCH · Search NASA
Results for “coupled problems”
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.
A higher order approximate static condensation method for multi-material diffusion problems
The paper studies an approximate static condensation method for the diffusion problem with discontinuous diffusion coefficients. The method allows for a general polygonal mesh which is unfitted to the material interfaces. Moreover, the interfaces can be discontinuous across the mesh edges as typical for numerical reconstructions using the volume or moment-of-fluid methods. We apply a mimetic finite difference method to solve local diffusion problems and use P 1 (mortar) edge elements to couple local problems into the global system. The condensation process and the properties of the resulting algebraic system are discussed. It is demonstrated that the method is second order accurate on smooth solutions and performs well for problems with high contrast in diffusion coefficients. Experiments also show the robustness with respect to position of the material interfaces against the underlying mesh.
97 MATHEMATICS AND COMPUTING↗
Projection-Based Model Reduction for Coupled Conduction—Enclosure Radiation Systems
We report a projection-based reduced order model (pROM) methodology has been developed for transient heat transfer problems involving coupled conduction and enclosure radiation. The approach was demonstrated on two test problems of varying complexity. The reduced order models demonstrated substantial speedups (up to 185×) relative to the full order model with good accuracy (less than 3% L ∞ error). An attractive feature of pROMs is that there is a natural error indicator for the ROM solution: the final residual norm at each time-step of the converged ROM solution. Using example test cases, we discuss how to interpret this error indicator to assess the accuracy of the ROM solution. The approach shows promise for many-query applications, such as uncertainty quantification and optimization. The reduced computational cost of the ROM relative to the full-order model (FOM) can enable the analysis of larger and more complex systems as well as the exploration of larger parameter spaces.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media
Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Partitioned exponential methods for coupled multiphysics systems
Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Challenges and Prospects for Reducing Coupled Climate Model SST Biases in the Eastern Tropical Atlantic and Pacific Oceans: The U.S. CLIVAR Eastern Tropical Oceans Synthesis Working Group
Well-known problems trouble coupled general circulation models of the eastern Atlantic and Pacific Ocean basins. Model climates are significantly more symmetric about the equator than is observed. Model sea surface temperatures are biased warm south and southeast of the equator, and the atmosphere is too rainy within a band south of the equator. Near-coastal eastern equatorial SSTs are too warm, producing a zonal SST gradient in the Atlantic opposite in sign to that observed. The U.S. Climate Variability and Predictability Program (CLIVAR) Eastern Tropical Ocean Synthesis Working Group (WG) has pursued here an updated assessment of coupled model SST biases, focusing on the surface energy balance components, on regional error sources from clouds, deep convection, winds, and ocean eddies; on the sensitivity to model resolution; and on remote impacts. Motivated by the assessment, the WG makes the following recommendations: 1) encourage identification of the specific parameterizations contributing to the biases in individual models, as these can be model dependent; 2) restrict multimodel intercomparisons to specific processes; 3) encourage development of high-resolution coupled models with a concurrent emphasis on parameterization development of finer-scale ocean and atmosphere features, including low clouds; 4) encourage further availability of all surface flux components from buoys, for longer continuous time periods, in persistently cloudy regions; and 5) focus on the eastern basin coastal oceanic upwelling regions, where further opportunities for observational–modeling synergism exist.
58 GEOSCIENCES↗
Investigation of the OECD/NEA PWR MOX/UO{sub 2} core transient benchmark using a coupled whole core pin-by-pin route in WIMS
This paper investigates the new CAMELOT coupling route developed in WIMS, which uses prepared tabulated cross-sections combined with the flux solver MERLIN module and the integrated sub-channel thermal hydraulics solver ARTHUR module to solve for the coupled core state. This coupling route is applied to the OECD/NEA pressurized water reactor (PWR) core transient benchmark exercise, involving 3-D neutronics and thermal-hydraulics analysis of a full-sized PWR core featuring MOX fuel at various core states, as well as simulating a control rod ejection transient. Analysis of the results for each part of the benchmark demonstrate that the high-resolution pin-by-pin analysis used in WIMS closely matches the other participants in the benchmark, with differences of between 20-200 pcm for k-effective, 1-2.5% RMS difference for assembly-averaged powers, 2-8% difference in rod worth calculations, < 200 ppm for calculations of critical boron concentration, zero difference in delayed neutron fraction calculation and a very close match for power response to the rod ejection transient. This is the first application of the WIMS-CAMELOT approach for the transient analysis of a full-sized reactor and as such further improvements such as reducing computational cost through parallelization, optimisation and memory reduction, are currently in active development. This paper demonstrates the accuracy and flexibility of the CAMELOT coupling route, providing a straightforward and consistent user image that can be easily used for a variety of modelling problems involving coupled neutronics and thermal-hydraulics. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Preconditioners for multiphase poromechanics with strong capillarity
This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.
42 ENGINEERING↗
Pure electronic metal-insulator transition at the interface of complex oxides
We observed complex materials in electronic phases and transitions between them often involve coupling between many degrees of freedom whose entanglement convolutes understanding of the instigating mechanism. Metal-insulator transitions are one such problem where coupling to the structural, orbital, charge, and magnetic order parameters frequently obscures the underlying physics. We demonstrate a way to unravel this conundrum by heterostructuring a prototypical multi-ordered complex oxide NdNiO3 in ultra thin geometry, which preserves the metal-to-insulator transition and bulk-like magnetic order parameter, but entirely suppresses the symmetry lowering and long-range charge order parameter. Furthermore, these findings illustrate the utility of heterointerfaces as a powerful method for removing competing order parameters to gain greater insight into the nature of the transition, here revealing that the magnetic order generates the transition independently, leading to an exceptionally rare purely electronic metal-insulator transition with no symmetry change.
36 MATERIALS SCIENCE↗
DeepM&Mnet: Inferring the electroconvection multiphysics fields based on operator approximation by neural networks
Electroconvection is a multiphysics problem involving coupling of the flow field with the electric field as well as the cation and anion concentration fields. Here, we use electroconvection as a benchmark problem to put forward a new data assimilation framework, the DeepM&Mnet, for simulating multiphysics and multiscale problems at speeds much faster than standard numerical methods using pre-trained neural networks. We first pre-train DeepONets that can predict independently each field, given general inputs from the rest of the fields of the coupled system. DeepONets can approximate nonlinear operators and are composed of two sub-networks, a branch net for the input fields and a trunk net for the locations of the output field. DeepONets, which are extremely fast, are used as building blocks in the DeepM&Mnet and form constraints for the multiphysics solution along with some sparse available measurements of any of the fields. We demonstrate the new methodology and document the accuracy of each individual DeepONet, and subsequently we present two different DeepM&Mnet architectures that infer accurately and efficiently 2D electroconvection fields for unseen electric potentials. Furthermore, the DeepM&Mnet framework is general and can be applied for building any complex multiphysics and multiscale models based on very few measurements using pre-trained DeepONets in a “plug-and-play” mode.
97 MATHEMATICS AND COMPUTING↗
An energy-based coupling approach to nonlocal interface problems.
Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, settings in which traditional partial differential equation models fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces, such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a physically consistent interface theory which is needed to support numerical developments and, among other features, reduces to classical models in the limit as the extent of nonlocal interactions vanish. In this paper, we use an energy-based approach to develop a formulation of a nonlocal interface problem which provides a physically consistent extension of the classical perfect interface formulation for partial differential equations. Numerical examples in one and two dimensions validate the proposed framework and demonstrate the scope of our theory.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations
While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.
42 ENGINEERING↗
A sharp immersed method for 2D flow-body interactions using the vorticity-velocity Navier-Stokes equations
Immersed methods discretize boundary conditions for complex geometries on background Cartesian grids. Here, this makes such methods especially suitable for two-way coupled flow-body problems, where the body mechanics are partially driven by hydrodynamic forces. However, for the vorticity-velocity form of the Navier-Stokes equations, existing immersed geometry discretizations for two-way coupled problems only achieve first order spatial accuracy near solid boundaries. Here we introduce a sharp-interface approach based on the immersed interface method to handle the one- and two-way coupling between an incompressible flow and one or more rigid bodies using the 2D vorticity-velocity Navier-Stokes equations. Our main contributions are three-fold. First, we develop and analyze a moving boundary treatment for sharp immersed methods that can be applied to PDEs with implicitly defined boundary conditions, such as those commonly imposed on the vorticity field. Second, we develop a two-way coupling methodology for the vorticity-velocity Navier-Stokes equations based on control-volume momentum balance that does not require the pressure field. Third, we show through extensive testing and validation that our resulting flow-body solver reaches second-order accuracy for most practical scenarios, and provides significant efficiency benefits compared to a representative first-order approach.
42 ENGINEERING↗
Development and Verification of a Numerical Library for Solving Global Terrestrial Multiphysics Problems
Abstract Current generation Land Surface Models (LSMs) routinely simulate many nonlinear multiphysics processes. The complexity of future generation LSMs is expected to increase as critical new biophysical and biogeochemical processes are incorporated. Current generation LSMs have several shortcomings including the lack of robust numerical methods to solve discretized equations, monolithic software design that hinders testing of a process representation in isolation from other components, and absence of a flexible coupling framework to solve tightly coupled multiphysics problems. While the LSMs community vigorously evaluates the accuracy of the conceptual model to represent reality (validation), the accuracy of the numerical implementation of the conceptual model (verification) is seldom examined. Method of Manufactured Solutions is a technique for verifying complex codes when analytical solutions are unavailable. In this work, we present a stand‐alone, open source numerical library for solving global terrestrial multiphysics processes that includes a flexible coupling framework. Robust numerical solution for linear and nonlinear equations is provided via the Portable, Extensible Toolkit for Scientific Computation library. We verify the numerical library using Method of Manufactured Solutions for a range of problems comprising single and multiphysics steady state problems that are applied in one or multiple physical domains. This work provides an example of incorporating code verification as an integral part of Land Surface Model development activities.
58 GEOSCIENCES↗
An Asymptotically Compatible Coupling Formulation for Nonlocal Interface Problems with Jumps
Here, we introduce a mathematically rigorous formulation for a nonlocal interface problem with jumps and propose an asymptotically compatible finite element discretization for the weak form of the interface problem. After proving the well-posedness of the weak form, we demonstrate that solutions to the nonlocal interface problem converge to the corresponding local counterpart when the nonlocal data are appropriately prescribed. Several numerical tests in one and two dimensions show the applicability of our technique, its numerical convergence to exact nonlocal solutions, its convergence to the local limit when the horizons vanish, and its robustness with respect to the patch test.
97 MATHEMATICS AND COMPUTING↗
PowerSystems.jl and PowerSimulations.jl [SWR-18-23]
Data structures in Julia to enable power systems analysis. Also, PowerSimulations.jl is a Julia package for power system modeling and simulation found at https://github.com/NREL/PowerSimulations.jl. The objectives of the package are: Provide a flexible modeling framework that can accommodate problems of different complexity and at different time-scales. Streamline the construction of large scale optimization problems to avoid repetition of work when adding/modifying model details. Exploit Julia's capabilities to improve computational performance of large scale power system simulations. The flexible modeling framework is enabled through a modular set of capabilities that enable scalable power system analysis and exploration of new analysis methods. The modularity of PowerSimulations results from the structure of the simulations enabled by the package: Simulations define a set of problems that can be solved using numerical techniques. Problems are generated by expressing model formulations against system data For example, an annual production cost modeling simulation can be created by formulating a unit commitment model against system data to assemble a set of 365 daily time-coupled scheduling problems.
Lara, Jose Daniel↗
Development and benchmarking of transient nodal code SIMULATE5-K neutron kinetics solver
SIMULATE5-K is Studsvik's next generation best estimate transient code. The time dependent diffusion equation is solved with a nodal method consistent with that implemented in the licensed core design code SIMULATE5. Arbitrary number of neutron and delayed neutron precursor groups can be used. For the solution of the spatial problem, the coupling coefficients used to relate the node leakages are found by first converting the time dependent diffusion equation to a static diffusion equation with the use of flux and delayed neutron precursor dynamic frequencies. Once the static-like equations are obtained, the multi-group analytical nodal model is used to obtain the coupling coefficients, expressing the node leakage in terms of adjacent node average fluxes. The coupling coefficients are then inserted into the time dependent nodal balance equation. For the time integration, the time dependent neutron balance equation is solved with the frequency transformation method. The treatment of the temporal dependence yields a fixed source problem which can be solved utilizing the existing fixed-source methodology. The primary purpose of this paper is to describe the neutron kinetics methodology implemented in SIMULATE5-K. The accuracy of the method is demonstrated for a series of well-known, neutronic-only benchmark problems. (author)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past
The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.
42 ENGINEERING↗