A non-overlapping optimized Schwarz method for the heat equation with non linear boundary conditions and with applications to de-icing

Abstract When simulating complex physical phenomena such as aircraft icing or de-icing, several dedicated solvers often need to be strongly coupled. In this work, a non-overlapping Schwarz method is constructed with the unsteady simulation of de-icing as the targeted application. To do so, optimized coupling coefficients are first derived for the one dimensional unsteady heat equation with linear boundary conditions and for the steady heat equation with non-linear boundary conditions. The choice of these coefficients is shown to guarantee the convergence of the method. Using a linearization of the boundary conditions, the method is then extended to the case of a general unsteady heat conduction problem. The method is tested on simple cases and the convergence properties are assessed theoretically and numerically. Finally the method is applied to the simulation of an aircraft electrothermal de-icing problem in two dimensions.

[1]  P. Trontin,et al.  Description and assessment of the new ONERA 2D icing suite IGLOO2D , 2017 .

[2]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .

[3]  Victorita Dolean,et al.  An introduction to domain decomposition methods - algorithms, theory, and parallel implementation , 2015 .

[4]  Martin J. Gander,et al.  Schwarz Methods over the Course of Time , 2008 .

[5]  Martin J. Gander,et al.  Overlapping Schwarz Waveform Relaxation for the Heat Equation in N Dimensions , 2002 .

[6]  Martin J. Gander,et al.  Optimized Schwarz Methods , 2006, SIAM J. Numer. Anal..

[7]  M. Salaün,et al.  A mixed adhesion–brittle fracture model and its application to the numerical study of ice shedding mechanisms , 2016 .

[8]  Martin J. Gander,et al.  Optimized Schwarz methods for a diffusion problem with discontinuous coefficient , 2014, Numerical Algorithms.

[9]  Olof B. Widlund,et al.  An Overlapping Schwarz Algorithm for Almost Incompressible Elasticity , 2009, SIAM J. Numer. Anal..

[10]  Alain Dervieux,et al.  INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2012) Published online in Wiley Online Library (wileyonlinelibrary.com/journal/nmf). DOI: 10.1002/fld.3733 Efficiency and scalability of a two-level Schwarz algorithm for i , 2022 .

[11]  Philippe Villedieu,et al.  Numerical simulation and modeling of ice shedding: Process initiation , 2014 .

[12]  Philippe Villedieu,et al.  An implicit time marching Galerkin method for the simulation of icing phenomena with a triple layer model , 2018, Finite Elements in Analysis and Design.

[13]  Marc-Paul Errera,et al.  Optimal solutions of numerical interface conditions in fluid-structure thermal analysis , 2013, J. Comput. Phys..

[14]  H. Brezis.  Functional Analysis, Sobolev Spaces and Partial Differential Equations , 2010 .

[15]  A. V. Luikov,et al.  Analytical Heat Diffusion Theory , 1969 .

[16]  Alfio Quarteroni,et al.  Domain Decomposition Methods for Partial Differential Equations , 1999 .

[17]  Jeffrey W. Banks,et al.  A stable and accurate partitioned algorithm for conjugate heat transfer , 2017, J. Comput. Phys..

[18]  P. Gosselet,et al.  Non-overlapping domain decomposition methods in structural mechanics , 2006, 1208.4209.

[19]  Martin J. Gander,et al.  Optimized Schwarz Methods without Overlap for the Helmholtz Equation , 2002, SIAM J. Sci. Comput..

[20]  Martin J. Gander,et al.  Optimized Schwarz waveform relaxation for advection reaction diffusion equations in two dimensions , 2014, Numerische Mathematik.

[21]  M. Salaün,et al.  Two Dimensional Model of an Electro-Thermal Ice Protection System , 2013 .

[22]  P. Trontin,et al.  A robust coupling algorithm applied to thermal ice protection system unsteady modeling , 2014 .