IFOSMONDI Co-simulation Algorithm with Jacobian-Free Methods in PETSc

by   Yohan Eguillon, et al.

IFOSMONDI iterative algorithm for implicit co-simulation of coupled physical systems (introduced by the authors in july 2019 during the Simultech conference, p.176-186) enables us to solve the nonlinear coupling function while keeping the smoothness of interfaces without introducing a delay. Moreover, it automatically adapts the size of the steps between data exchanges among the systems according to the difficulty of the solving of the coupling constraint. The latter was solved by a fixed-point algorithm in the original implementation whereas this paper introduces the JFM version (standing for Jacobian-Free Methods). Most implementations of Newton-like methods require a jacobian matrix which can be difficult to compute in the co-simulation context, except in the case where the interfaces are represented by a Zero-Order-Hold (ZOH). As far as IFOSMONDI coupling algorithm uses Hermite interpolation for smoothness enhancement (up to Third-Order-Hold), we propose hereafter a new formulation of the non-linear coupling function including both the values and the time-derivatives of the coupling variables. This formulation is well designed for solving the coupling through jacobian-free Newton type methods. Consequently, successive function evaluations consist in multiple simulations of the systems on a co-simulation time step using rollback. The orchestrator-workers structure of the algorithm enables us to combine the PETSc framework on the orchestrator side for the non-linear Newton-type solvers with the parallel integrations of the systems on the workers side thanks to MPI processes. Different nonlinear methods will be compared to one another and to the original fixed-point implementation on a newly proposed 2-systems academic test-case (mass-spring-damper type) with direct feedthrough on both sides.


page 1

page 2

page 3

page 4


Quasi-Newton Waveform Iteration for Partitioned Fluid-Structure Interaction

We present novel coupling schemes for partitioned multi-physics simulati...

Treating Smoothness and Balance during Data Exchange in Explicit Simulator Coupling or Cosimulation

Cosimulation methods allow combination of simulation tools of physical s...

Newton Type Methods for solving a Hasegawa-Mima Plasma Model

In [1], the non-linear space-time Hasegawa-Mima plasma equation is formu...

F3ORNITS: A Flexible Variable Step Size Non-Iterative Co-simulation Method handling Subsystems with Hybrid Advanced Capabilities

This paper introduces the F3ORNITS non-iterative co-simulation algorithm...

SymX: Energy-based Simulation from Symbolic Expressions

Optimization time integrators have proven to be effective at solving com...

HPC compact quasi-Newton algorithm for interface problems

In this work we present a robust interface coupling algorithm called Com...

Please sign up or login with your details

Forgot password? Click here to reset