TY - JOUR
T1 - Proper Generalized Decomposition using Taylor expansion for non-linear diffusion equations
AU - Deeb, Ahmad
AU - Kalaoun, Omar
AU - Belarbi, Rafik
N1 - Publisher Copyright:
© 2023 International Association for Mathematics and Computers in Simulation (IMACS)
PY - 2023/6
Y1 - 2023/6
N2 - For a physical problem described by a parameterized mathematical model, different configurations of the problem require computing the solution over a range of parameters in order to study the phenomenon when parameters change. In other words, it is a process of looking for a continuum of solutions of the equation, relative to these parameters, in order to find the ones that fit the experimental data. However, solving a direct problem for each parametric configuration will generate a cascade of direct problems, which will cost a huge amount of time, especially when we deal with non-linear equations. Therefore, the parametric solution is a suitable alternative strategy to compute the solution of the equation. In this paper, we will use the Proper Generalized Decomposition (PGD) method to solve non-linear diffusion equations and produce parametric solutions. To treat the non-linear functions, we will not use the Discrete Empirical Interpolation Methods (DEIM), which has proven their utility, but the non-linear terms will be replaced by their Taylor series expansion up to an order m. This will produce a new model, which we call here the ”developed equation” and therefore the PGD is applied on. Polynomial equations appear for each tensor element computation. While space and time tensor elements’ equations are to be solved using Finite Elements Methods (FEM) and Borel–Padé–Laplace (BPL) integrator respectively, Newton solver is used for tensors relative to the parameters’ equations. Here, rational polynomial functions arise for parametric tensor elements, which are known to extrapolate solutions. Numerical simulations are done for a non-linear diffusion equation with exponential diffusion coefficient as first trial, and with a magnetic diffusion coefficient as a second one.
AB - For a physical problem described by a parameterized mathematical model, different configurations of the problem require computing the solution over a range of parameters in order to study the phenomenon when parameters change. In other words, it is a process of looking for a continuum of solutions of the equation, relative to these parameters, in order to find the ones that fit the experimental data. However, solving a direct problem for each parametric configuration will generate a cascade of direct problems, which will cost a huge amount of time, especially when we deal with non-linear equations. Therefore, the parametric solution is a suitable alternative strategy to compute the solution of the equation. In this paper, we will use the Proper Generalized Decomposition (PGD) method to solve non-linear diffusion equations and produce parametric solutions. To treat the non-linear functions, we will not use the Discrete Empirical Interpolation Methods (DEIM), which has proven their utility, but the non-linear terms will be replaced by their Taylor series expansion up to an order m. This will produce a new model, which we call here the ”developed equation” and therefore the PGD is applied on. Polynomial equations appear for each tensor element computation. While space and time tensor elements’ equations are to be solved using Finite Elements Methods (FEM) and Borel–Padé–Laplace (BPL) integrator respectively, Newton solver is used for tensors relative to the parameters’ equations. Here, rational polynomial functions arise for parametric tensor elements, which are known to extrapolate solutions. Numerical simulations are done for a non-linear diffusion equation with exponential diffusion coefficient as first trial, and with a magnetic diffusion coefficient as a second one.
KW - Diffusion equation
KW - Non-linearity
KW - Proper Generalized Decomposition
KW - Rational polynomial
KW - Taylor expansion
UR - http://www.scopus.com/inward/record.url?scp=85147191200&partnerID=8YFLogxK
U2 - 10.1016/j.matcom.2023.01.008
DO - 10.1016/j.matcom.2023.01.008
M3 - Article
AN - SCOPUS:85147191200
SN - 0378-4754
VL - 208
SP - 71
EP - 94
JO - Mathematics and Computers in Simulation
JF - Mathematics and Computers in Simulation
ER -