# American Institute of Mathematical Sciences

February  2018, 11(1): 21-34. doi: 10.3934/dcdss.2018002

## Optimal number of Schur subdomains: Application to semi-implicit finite volume discretization of semilinear reaction diffusion problem

 1 Université Abdelmalek Essaadi, Faculté des Sciences et Techniques de Tanger, Laboratoire de Mathématiques et Applications, BP 416,90000 Tanger, Maroc 2 Université Moulay Ismail, Faculté des Sciences et Techniques d'Errachidia, Laboratoire de Mathématiques Informatique et Image, Equipe d'Analyse Mathématique et Numérique des EDPs et Applications, BP 509, Boutalamine, 52000 Errachidia, Maroc

* Corresponding author: N. Nagid

Received  September 2016 Revised  April 2017 Published  January 2018

The purpose of this paper is to establish a new numerical approach to solve, in two dimensions, a semilinear reaction diffusion equation combining finite volume method and Schur complement method. We applied our method for q = 2 non-overlapping subdomains and then we generalized in the case of several subdomains (q≥2). A large number of numerical test cases shows the efficiency and the good accuracy of the proposed approach in terms of the CPU time and the order of the error, when increasing the number of subdomains, without using the parallel computing. After several variations of the number of subdomains and the mesh grid, we remark two significant results. On the one hand, the increase related to the number of subdomains does not affect the order of the error, on the other hand, for each mesh grid when we augment the number of subdomains, the CPU time reaches the minimum for a specific number of subdomains. In order to have the minimum CPU time, we resorted to a statistical study between the optimal number of subdomains and the mesh grid.

Citation: Hassan Belhadj, Mohamed Fihri, Samir Khallouq, Nabila Nagid. Optimal number of Schur subdomains: Application to semi-implicit finite volume discretization of semilinear reaction diffusion problem. Discrete & Continuous Dynamical Systems - S, 2018, 11 (1) : 21-34. doi: 10.3934/dcdss.2018002
##### References:

show all references

##### References:
FV structured mesh of domain Ω
Non-overlapping strip decomposition
Domain decomposition and structured conforming mesh of domain Ω
FE solution (left) and FV-1D (right), for h = 0:0138 and t = 0:5 in 3D
FE solution (left) and FV-1D (right), for h = 0:0138 and t = 0:5 in 2D
Point-wise errors between the solution by FV-1Dom and by FE method (with the basis function in P1 (left) and P2 (right)), for $h=0.0138$ and $t=0.5$
Point-wise errors between the solution by FV-SC (2SD) and by FE method (with the basis function in P1 (left) and P2 (right)), for $h=0.0138$ and $t=0.5$
$L^{\infty}(\Omega)$ error history between FE and FV-1Dom (left) (resp. FV-SC (right)), with the basis functions in P1 (FE-P1) and P2 (FE-P2), for $h=0.0138$ and $t=1$
FV-1Dom solutions for different values of $h$, $y=0.0833$ (left) and $y=1.9167$ (right) at $t=0.5$
FV-SC (2 subdomains) solutions for different values of $h$, $y=0.0833$ (left) and $y=1.9167$ (right) at $t=0.5$
CPU time in seconds for FV-1Dom and FV-SC (h = 0:0138 and t = 1)
Graphic representation of the adjustment to the power model
$L^{2}(\Omega)$ errors between FV and FE-P2 ($h_{FE}=0.0069$) for different values of $h$ and of the number of subdomains, at $t=0.5$
 Number of Subdomains 0.1666 h 0.0138 0.0069 1 Domain 0.0243 1.7708E-4 9.8594e-05 2 Subdomains 0.0243 1.7221E-4 8.9848e-05 3 Subdomains 0.0243 1.7430E-4 6.9495e-05 4 Subdomains 0.0243 1.7600E-4 7.0203e-05 6 Subdomains 0.0243 1.7485E-4 7.1316e-05 8 Subdomains - 1.7375E-4 6.922e-05 9 Subdomains - 1.7357E-4 7.0606e-05
 Number of Subdomains 0.1666 h 0.0138 0.0069 1 Domain 0.0243 1.7708E-4 9.8594e-05 2 Subdomains 0.0243 1.7221E-4 8.9848e-05 3 Subdomains 0.0243 1.7430E-4 6.9495e-05 4 Subdomains 0.0243 1.7600E-4 7.0203e-05 6 Subdomains 0.0243 1.7485E-4 7.1316e-05 8 Subdomains - 1.7375E-4 6.922e-05 9 Subdomains - 1.7357E-4 7.0606e-05
CPU time calculation for FV-1Dom and FV-SC at t = 1
 Number of Subdomains 0.1666 h 0.0138 0.0069 1 Domain 0.1139 (second) 0.6285 (hour) 1.9884 (hour) 2 Subdomains 0.4891 (second) 0.4084 (hour) 1.4347 (hour) 3 Subdomains 0.4602 (second) 0.3264 (hour) 1.3845 (hour) 4 Subdomains 0.4823 (second) 0.3198 (hour) 1.4279 (hour) 6 Subdomains 0.5081 (second) 0.3299 (hour) 1.6049 (hour) 8 Subdomains - 0.3555 (hour) 1.1068 (hour) 9 Subdomains - 0.3700 (hour) 1.1577 (hour)
 Number of Subdomains 0.1666 h 0.0138 0.0069 1 Domain 0.1139 (second) 0.6285 (hour) 1.9884 (hour) 2 Subdomains 0.4891 (second) 0.4084 (hour) 1.4347 (hour) 3 Subdomains 0.4602 (second) 0.3264 (hour) 1.3845 (hour) 4 Subdomains 0.4823 (second) 0.3198 (hour) 1.4279 (hour) 6 Subdomains 0.5081 (second) 0.3299 (hour) 1.6049 (hour) 8 Subdomains - 0.3555 (hour) 1.1068 (hour) 9 Subdomains - 0.3700 (hour) 1.1577 (hour)
Optimal number of subdomains for different values of h
 Optimal Number of subdomains h CPU time 1 Domain 0.1666 0.2018 (second) 2 Subdomains 0.083333 1.2547 (second) 3 Subdomains 0.01388 0.1344 (hour) 4 Subdomains 0.012820 0.4736 (hour) 5 Subdomains 0.00952 1.0142 (hour) 6 Subdomains 0.00925 0.9637 (hour) 7 Subdomains 0.00793 0.9124 (hour) 8 Subdomains 0.00694 1.0741 (hour) 9 Subdomains 0.00347 1.1570 (hour) 10 Subdomains 0.002222 1.4655 (hour) 11 Subdomains 0.002164 1.5632 (hour) 12 Subdomains 0.0021367 1.5932 (hour) 13 Subdomains 0.002051 2.1130 (hour) 35 Subdomains 0.00029304 5.4254 (hour)
 Optimal Number of subdomains h CPU time 1 Domain 0.1666 0.2018 (second) 2 Subdomains 0.083333 1.2547 (second) 3 Subdomains 0.01388 0.1344 (hour) 4 Subdomains 0.012820 0.4736 (hour) 5 Subdomains 0.00952 1.0142 (hour) 6 Subdomains 0.00925 0.9637 (hour) 7 Subdomains 0.00793 0.9124 (hour) 8 Subdomains 0.00694 1.0741 (hour) 9 Subdomains 0.00347 1.1570 (hour) 10 Subdomains 0.002222 1.4655 (hour) 11 Subdomains 0.002164 1.5632 (hour) 12 Subdomains 0.0021367 1.5932 (hour) 13 Subdomains 0.002051 2.1130 (hour) 35 Subdomains 0.00029304 5.4254 (hour)
 [1] Qingping Deng. A nonoverlapping domain decomposition method for nonconforming finite element problems. Communications on Pure & Applied Analysis, 2003, 2 (3) : 297-310. doi: 10.3934/cpaa.2003.2.297 [2] Xin Li, Xingfu Zou. On a reaction-diffusion model for sterile insect release method with release on the boundary. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2509-2522. doi: 10.3934/dcdsb.2012.17.2509 [3] Patrick De Kepper, István Szalai. An effective design method to produce stationary chemical reaction-diffusion patterns. Communications on Pure & Applied Analysis, 2012, 11 (1) : 189-207. doi: 10.3934/cpaa.2012.11.189 [4] Mihaela Negreanu, J. Ignacio Tello. On a comparison method to reaction-diffusion systems and its applications to chemotaxis. Discrete & Continuous Dynamical Systems - B, 2013, 18 (10) : 2669-2688. doi: 10.3934/dcdsb.2013.18.2669 [5] Christos V. Nikolopoulos, Georgios E. Zouraris. Numerical solution of a non-local elliptic problem modeling a thermistor with a finite element and a finite volume method. Conference Publications, 2007, 2007 (Special) : 768-778. doi: 10.3934/proc.2007.2007.768 [6] Fang Zeng, Pablo Suarez, Jiguang Sun. A decomposition method for an interior inverse scattering problem. Inverse Problems & Imaging, 2013, 7 (1) : 291-303. doi: 10.3934/ipi.2013.7.291 [7] Jing Xu, Xue-Cheng Tai, Li-Lian Wang. A two-level domain decomposition method for image restoration. Inverse Problems & Imaging, 2010, 4 (3) : 523-545. doi: 10.3934/ipi.2010.4.523 [8] Klemens Fellner, Wolfang Prager, Bao Q. Tang. The entropy method for reaction-diffusion systems without detailed balance: First order chemical reaction networks. Kinetic & Related Models, 2017, 10 (4) : 1055-1087. doi: 10.3934/krm.2017042 [9] Takashi Kajiwara. The sub-supersolution method for the FitzHugh-Nagumo type reaction-diffusion system with heterogeneity. Discrete & Continuous Dynamical Systems - A, 2018, 38 (5) : 2441-2465. doi: 10.3934/dcds.2018101 [10] Shuang Liu, Xinfeng Liu. Krylov implicit integration factor method for a class of stiff reaction-diffusion systems with moving boundaries. Discrete & Continuous Dynamical Systems - B, 2017, 22 (11) : 1-19. doi: 10.3934/dcdsb.2019176 [11] Yuchi Qiu, Weitao Chen, Qing Nie. A hybrid method for stiff reaction–diffusion equations. Discrete & Continuous Dynamical Systems - B, 2019, 24 (12) : 6387-6417. doi: 10.3934/dcdsb.2019144 [12] Martin Burger, José A. Carrillo, Marie-Therese Wolfram. A mixed finite element method for nonlinear diffusion equations. Kinetic & Related Models, 2010, 3 (1) : 59-83. doi: 10.3934/krm.2010.3.59 [13] Zhongyi Huang. Tailored finite point method for the interface problem. Networks & Heterogeneous Media, 2009, 4 (1) : 91-106. doi: 10.3934/nhm.2009.4.91 [14] Mostafa Bendahmane, Mauricio Sepúlveda. Convergence of a finite volume scheme for nonlocal reaction-diffusion systems modelling an epidemic disease. Discrete & Continuous Dynamical Systems - B, 2009, 11 (4) : 823-853. doi: 10.3934/dcdsb.2009.11.823 [15] Feng Ma, Mingfang Ni. A two-phase method for multidimensional number partitioning problem. Numerical Algebra, Control & Optimization, 2013, 3 (2) : 203-206. doi: 10.3934/naco.2013.3.203 [16] Caterina Calgaro, Meriem Ezzoug, Ezzeddine Zahrouni. Stability and convergence of an hybrid finite volume-finite element method for a multiphasic incompressible fluid model. Communications on Pure & Applied Analysis, 2018, 17 (2) : 429-448. doi: 10.3934/cpaa.2018024 [17] Xiaomao Deng, Xiao-Chuan Cai, Jun Zou. A parallel space-time domain decomposition method for unsteady source inversion problems. Inverse Problems & Imaging, 2015, 9 (4) : 1069-1091. doi: 10.3934/ipi.2015.9.1069 [18] Xu Yang, François Golse, Zhongyi Huang, Shi Jin. Numerical study of a domain decomposition method for a two-scale linear transport equation. Networks & Heterogeneous Media, 2006, 1 (1) : 143-166. doi: 10.3934/nhm.2006.1.143 [19] Shi Jin, Xu Yang, Guangwei Yuan. A domain decomposition method for a two-scale transport equation with energy flux conserved at the interface. Kinetic & Related Models, 2008, 1 (1) : 65-84. doi: 10.3934/krm.2008.1.65 [20] Stefan Berres, Ricardo Ruiz-Baier, Hartmut Schwandt, Elmer M. Tory. An adaptive finite-volume method for a model of two-phase pedestrian flow. Networks & Heterogeneous Media, 2011, 6 (3) : 401-423. doi: 10.3934/nhm.2011.6.401

2018 Impact Factor: 0.545

## Metrics

• HTML views (117)
• Cited by (0)

• on AIMS