Proper generalized decomposition of a geometrically parametrized heat problem with geophysical applications

The solution of a steady thermal multiphase problem is assumed to be dependent on a set of parameters describing the geometry of the domain, the internal interfaces and the material properties. These parameters are considered as new independent variables. The problem is therefore stated in a multidimensional setup. The proper generalized decomposition (PGD) provides an approximation scheme especially well suited to preclude dramatically increasing the computational complexity with the number of dimensions. The PGD strategy is reviewed for the standard case dealing only with material parameters. Then, the ideas presented in [Ammar et al., Parametric solutions involving geometry: A step towards efficient shape optimization. Comput. Methods Appl. Mech. Eng., 2014; 268:178-193] to deal with parameters describing the domain geometry are adapted to a more general case including parametrization of the location of internal interfaces. Finally, the formulation is extended to combine the two types of parameters. The proposed strategy is used to solve a problem in applied geophysics studying the temperature field in a cross section of the Earth crust subsurface. The resulting problem is in a 10-dimensional space, but the PGD solution provides a fairly accurate approximation (error 1%) using less that 150 terms in the PGD expansion. Copyright (c) 2015John Wiley & Sons, Ltd.

[1]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[2]  Francisco Chinesta,et al.  Recent Advances and New Challenges in the Use of the Proper Generalized Decomposition for Solving Multidimensional Models , 2010 .

[3]  F. Chinesta,et al.  Recent advances on the use of separated representations , 2009 .

[4]  Francisco Chinesta,et al.  A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modeling of complex fluids , 2006 .

[5]  Francisco Chinesta,et al.  Alleviating mesh constraints : Model reduction, parallel time integration and high resolution homogenization , 2008 .

[6]  János Urai,et al.  Multiphase salt tectonic evolution in NW Germany: seismic interpretation and retro-deformation , 2005 .

[7]  G. Rozza,et al.  On the stability of the reduced basis method for Stokes equations in parametrized domains , 2007 .

[8]  T. Belytschko,et al.  Topology optimization with implicit functions and regularization , 2003 .

[9]  Gianluigi Rozza,et al.  A reduced basis hybrid method for the coupling of parametrized domains represented by fluidic networks , 2012 .

[10]  A. Nouy.  A generalized spectral decomposition technique to solve a class of linear stochastic partial differential equations , 2007 .

[11]  Hermann Zeyen,et al.  Modelling of thermal anomalies in the NW border of the Valencia Trough by groundwater convection , 1990 .

[12]  J. Chun,et al.  Shape optimization of electromagnetic devices using immune algorithm , 1997 .

[13]  A. Huerta,et al.  Parametric solutions involving geometry: A step towards efficient shape optimization , 2014 .

[14]  Ole Sigmund,et al.  A monolithic approach for topology optimization of electrostatically actuated devices , 2008 .

[15]  L. H. Olesen,et al.  A high‐level programming‐language implementation of topology optimization applied to steady‐state Navier–Stokes flow , 2004, physics/0410086.

[16]  Francisco Chinesta,et al.  A new family of solvers for some classes of multidimensional partial differential equations encountered in kinetic theory modelling of complex fluids - Part II: Transient simulation using space-time separated representations , 2007 .

[17]  Mathilde Chevreuil,et al.  A multiscale method with patch for the solution of stochastic partial differential equations with localized uncertainties , 2013 .

[18]  A. Ammar,et al.  PGD-Based Computational Vademecum for Efficient Design, Optimization and Control , 2013, Archives of Computational Methods in Engineering.

[19]  Juan Carlos Afonso,et al.  Density structure and buoyancy of the oceanic lithosphere revisited , 2007 .

[20]  S Niroomandi,et al.  Real‐time simulation of biological soft tissues: a PGD approach , 2013, International journal for numerical methods in biomedical engineering.

[21]  Antonio Falcó,et al.  A Proper Generalized Decomposition for the solution of elliptic problems in abstract form by using a functional Eckart–Young approach , 2011 .

[22]  Mathilde Chevreuil,et al.  Fictitious domain method and separated representations for the solution of boundary value problems on uncertain parameterized domains , 2011 .

[23]  Joos Vandewalle,et al.  A Multilinear Singular Value Decomposition , 2000, SIAM J. Matrix Anal. Appl..

[24]  Hermann Zeyen,et al.  Integrated lithospheric modeling combining thermal, gravity, and local isostasy analysis: Application to the NE Spanish Geotransect , 1994 .

[25]  Pedro Díez,et al.  An error estimator for separated representations of highly multidimensional models , 2010 .