Final Report on LDRD Project: Coupling Strategies for Multi-Physics Applications

Many current and future modeling applications at Sandia including ASC milestones will critically depend on the simultaneous solution of vastly different physical phenomena. Issues due to code coupling are often not addressed, understood, or even recognized. The objectives of the LDRD has been both in theory and in code development. We will show that we have provided a fundamental analysis of coupling, i.e., when strong coupling vs. a successive substitution strategy is needed. We have enabled the implementation of tighter coupling strategies through additions to the NOX and Sierra code suites to make coupling strategies available now. We have leveraged existing functionality to do this. Specifically, we have built into NOX the capability to handle fully coupled simulations from multiple codes, and we have also built into NOX the capability to handle Jacobi Free Newton Krylov simulations that link multiple applications. We show how this capability may be accessed from within the Sierra Framework as well as from outside of Sierra. The critical impact from this LDRD is that we have shown how and have delivered strategies for enabling strong Newton-based coupling while respecting the modularity of existing codes. This will facilitate the use of these codes in a coupled manner to solve multi-physic applications.

[1]  J. Shadid,et al.  Studies of the Accuracy of Time Integration Methods for Reaction-Diffusion Equations ∗ , 2005 .

[2]  A. Sangiovanni-Vincentelli,et al.  A multilevel Newton algorithm with macromodeling and latency for the analysis of large-scale nonlinear circuits in the time domain , 1979 .

[3]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[4]  Rainald Löhner,et al.  Edge-based finite element scheme for the Euler equations , 1994 .

[5]  J. Shadid,et al.  Studies on the accuracy of time-integration methods for the radiation-diffusion equations , 2004 .

[6]  Ljiljana Trajkovic,et al.  Artificial parameter homotopy methods for the DC operating point problem , 1993, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[7]  John N. Shadid,et al.  Comparison of Operators for Newton-Krylov Method for Solving Compressible Flows on Unstructured Meshes , 2004 .

[8]  James Wilkes,et al.  A New Method for Solving Systems of Nonlinear Equations in Circuit Simulation , 1994 .

[9]  Thomas M. Smith,et al.  Intelligent Nonlinear Solvers for Computational Fluid Dynamics , 2006 .

[10]  Kartikeya Mayaram,et al.  Coupling algorithms for mixed-level circuit and device simulation , 1992, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[11]  Arun Pandy,et al.  Fixed‐point convergence of modular, steady‐state heat transfer models coupling multiple scales and phenomena for melt–crystal growth , 2006 .

[12]  Roman Weber,et al.  Radiation in Enclosures , 2000 .

[13]  Harry K. Moffat,et al.  CADS:Cantera Aerosol Dynamics Simulator. , 2007 .

[14]  Naji Qatanani,et al.  Analytical and numerical investigation of the Fredholm integral equation for the heat radiation problem , 2006, Appl. Math. Comput..

[15]  Jiri Blazek,et al.  Accurate and Efficient Discretization of Navier-Stokes Equations on Mixed Grids , 2000 .

[16]  C. L. Hilmes,et al.  Nonlinear Elimination in Aerodynamic Analysis and Design Optimization , 2003 .

[17]  Paul J. Lanzkron,et al.  An Analysis of Approximate Nonlinear Elimination , 1996, SIAM J. Sci. Comput..

[18]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[19]  Eric R. Keiter,et al.  Xyce parallel electronic simulator design : mathematical formulation, version 2.0. , 2004 .

[20]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[21]  Louis A. Romero,et al.  Bifurcation Tracking Algorithms and Software for Large Scale Applications , 2005, Int. J. Bifurc. Chaos.

[22]  M. Pinar Mengüç,et al.  Thermal Radiation Heat Transfer , 2020 .