A tightlycoupled domaindecomposition approach for highly nonlinear stochastic multiphysics systems
Abstract
Multiphysics simulations often involve nonlinear components that are driven by internally generated or externally imposed random fluctuations. When used with a domaindecomposition (DD) algorithm, such components have to be coupled in a way that both accurately propagates the noise between the subdomains and lends itself to a stable and costeffective temporal integration. We develop a conservative DD approach in which tight coupling is obtained by using a Jacobianfree Newton–Krylov (JfNK) method with a generalized minimum residual iterative linear solver. This strategy is tested on a coupled nonlinear diffusion system forced by a truncated Gaussian noise at the boundary. Enforcement of pathwise continuity of the state variable and its flux, as opposed to continuity in the mean, at interfaces between subdomains enables the DD algorithm to correctly propagate boundary fluctuations throughout the computational domain. Reliance on a single Newton iteration (explicit coupling), rather than on the fully converged JfNK (implicit) coupling, may increase the solution error by an order of magnitude. Increase in communication frequency between the DD components reduces the explicit coupling's error, but makes it less efficient than the implicit coupling at comparable error levels for all noise strengths considered. Finally, the DD algorithm with the implicit JfNK couplingmore »
 Authors:
 Publication Date:
 OSTI Identifier:
 22622251
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 330; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; CORRELATIONS; COUPLING; DECOMPOSITION; DIFFUSION; ENFORCEMENT; ERRORS; FLUCTUATIONS; HYDROGEN; INTERFACES; ITERATIVE METHODS; MATHEMATICAL SOLUTIONS; NOISE; NONLINEAR PROBLEMS; RANDOMNESS; STOCHASTIC PROCESSES
Citation Formats
Taverniers, Søren, and Tartakovsky, Daniel M., Email: dmt@ucsd.edu. A tightlycoupled domaindecomposition approach for highly nonlinear stochastic multiphysics systems. United States: N. p., 2017.
Web. doi:10.1016/J.JCP.2016.10.052.
Taverniers, Søren, & Tartakovsky, Daniel M., Email: dmt@ucsd.edu. A tightlycoupled domaindecomposition approach for highly nonlinear stochastic multiphysics systems. United States. doi:10.1016/J.JCP.2016.10.052.
Taverniers, Søren, and Tartakovsky, Daniel M., Email: dmt@ucsd.edu. Wed .
"A tightlycoupled domaindecomposition approach for highly nonlinear stochastic multiphysics systems". United States.
doi:10.1016/J.JCP.2016.10.052.
@article{osti_22622251,
title = {A tightlycoupled domaindecomposition approach for highly nonlinear stochastic multiphysics systems},
author = {Taverniers, Søren and Tartakovsky, Daniel M., Email: dmt@ucsd.edu},
abstractNote = {Multiphysics simulations often involve nonlinear components that are driven by internally generated or externally imposed random fluctuations. When used with a domaindecomposition (DD) algorithm, such components have to be coupled in a way that both accurately propagates the noise between the subdomains and lends itself to a stable and costeffective temporal integration. We develop a conservative DD approach in which tight coupling is obtained by using a Jacobianfree Newton–Krylov (JfNK) method with a generalized minimum residual iterative linear solver. This strategy is tested on a coupled nonlinear diffusion system forced by a truncated Gaussian noise at the boundary. Enforcement of pathwise continuity of the state variable and its flux, as opposed to continuity in the mean, at interfaces between subdomains enables the DD algorithm to correctly propagate boundary fluctuations throughout the computational domain. Reliance on a single Newton iteration (explicit coupling), rather than on the fully converged JfNK (implicit) coupling, may increase the solution error by an order of magnitude. Increase in communication frequency between the DD components reduces the explicit coupling's error, but makes it less efficient than the implicit coupling at comparable error levels for all noise strengths considered. Finally, the DD algorithm with the implicit JfNK coupling resolves temporallycorrelated fluctuations of the boundary noise when the correlation time of the latter exceeds some multiple of an appropriately defined characteristic diffusion time.},
doi = {10.1016/J.JCP.2016.10.052},
journal = {Journal of Computational Physics},
number = ,
volume = 330,
place = {United States},
year = {Wed Feb 01 00:00:00 EST 2017},
month = {Wed Feb 01 00:00:00 EST 2017}
}

Multiphysics problems often involve components whose macroscopic dynamics is driven by microscopic random fluctuations. The fidelity of simulations of such systems depends on their ability to propagate these random fluctuations throughout a computational domain, including subdomains represented by deterministic solvers. When the constituent processes take place in nonoverlapping subdomains, system behavior can be modeled via a domaindecomposition approach that couples separate components at the interfaces between these subdomains. Its coupling algorithm has to maintain a stable and efficient numerical time integration even at high noise strength. We propose a conservative domaindecomposition algorithm in which tight coupling is achieved by employingmore »

Tightly Coupled Multiphysics Algorithm for Pebble Bed Reactors
We have developed a tightly coupled multiphysics simulation tool for the pebblebed reactor (PBR) concept, a type of Very HighTemperature gascooled Reactor (VHTR). The simulation tool, PRONGHORN, takes advantages of the Multiphysics ObjectOriented Simulation Environment library, and is capable of solving multidimensional thermalfluid and neutronics problems implicitly with a Newtonbased approach. Expensive Jacobian matrix formation is alleviated via the Jacobianfree NewtonKrylov method, and physicsbased preconditioning is applied to minimize Krylov iterations. Motivation for the work is provided via analysis and numerical experiments on simpler multiphysics reactor models. We then provide detail of the physical models and numerical methods in PRONGHORN.more » 
The Role of Data Transfer on the Selection of a Single vs Multiple Mesh Architecture for Tightly Coupled Multiphysics Applications
Data transfer from one distinct mesh to another may be necessary in any number of applications, including prolongation operations supporting multigrid solution methods, spatial adaptation, remeshing, and arbitrary LagrangianEulerian (ALE) and multiphysics simulation. This data transfer process is also referred to as remapping, rezoning and interpolation. Intermesh data transfer has the potential to introduce error into a simulation; the magnitude and importance of which depends on the transfer scenario and the algorithm used to perform the transfer. For a transient analysis, data transfer may occur many times during a simulation, with possible error accumulation at each transfer. The present studymore »