scholarly journals A note on the numerical approach for the reaction–diffusion problem to model the density of the tumor growth dynamics

2015 ◽  
Vol 69 (12) ◽  
pp. 1504-1517 ◽  
Author(s):  
E. Özuğurlu
2020 ◽  
Vol 26 ◽  
pp. 104
Author(s):  
Carlo Orrieri ◽  
Elisabetta Rocca ◽  
Luca Scarpa

We study a stochastic phase-field model for tumor growth dynamics coupling a stochastic Cahn-Hilliard equation for the tumor phase parameter with a stochastic reaction-diffusion equation governing the nutrient proportion. We prove strong well-posedness of the system in a general framework through monotonicity and stochastic compactness arguments. We introduce then suitable controls representing the concentration of cytotoxic drugs administered in medical treatment and we analyze a related optimal control problem. We derive existence of an optimal strategy and deduce first-order necessary optimality conditions by studying the corresponding linearized system and the backward adjoint system.


2010 ◽  
Vol 51 (3) ◽  
pp. 317-330
Author(s):  
E. ÖZUĞURLU

AbstractThe equation modelling the evolution of a foam (a complex porous medium consisting of a set of gas bubbles surrounded by liquid films) is solved numerically. This model is described by the reaction–diffusion differential equation with a free boundary. Two numerical methods, namely the fixed-point and the averaging in time and forward differences in space (the Crank–Nicolson scheme), both in combination with Newton’s method, are proposed for solving the governing equations. The solution of Burgers’ equation is considered as a special case. We present the Crank–Nicolson scheme combined with Newton’s method for the reaction–diffusion differential equation appearing in a foam breaking phenomenon.


2020 ◽  
Vol 28 (3) ◽  
pp. 147-160
Author(s):  
Andrea Bonito ◽  
Diane Guignard ◽  
Ashley R. Zhang

AbstractWe consider the numerical approximation of the spectral fractional diffusion problem based on the so called Balakrishnan representation. The latter consists of an improper integral approximated via quadratures. At each quadrature point, a reaction–diffusion problem must be approximated and is the method bottle neck. In this work, we propose to reduce the computational cost using a reduced basis strategy allowing for a fast evaluation of the reaction–diffusion problems. The reduced basis does not depend on the fractional power s for 0 < smin ⩽ s ⩽ smax < 1. It is built offline once for all and used online irrespectively of the fractional power. We analyze the reduced basis strategy and show its exponential convergence. The analytical results are illustrated with insightful numerical experiments.


Author(s):  
Oluwaseun Adeyeye ◽  
Ali Aldalbahi ◽  
Jawad Raza ◽  
Zurni Omar ◽  
Mostafizur Rahaman ◽  
...  

AbstractThe processes of diffusion and reaction play essential roles in numerous system dynamics. Consequently, the solutions of reaction–diffusion equations have gained much attention because of not only their occurrence in many fields of science but also the existence of important properties and information in the solutions. However, despite the wide range of numerical methods explored for approximating solutions, the adoption of block methods is yet to be investigated. Hence, this article introduces a new two-step third–fourth-derivative block method as a numerical approach to solve the reaction–diffusion equation. In order to ensure improved accuracy, the method introduces the concept of nonlinearity in the solution of the linear model through the presence of higher derivatives. The method obtained accurate solutions for the model at varying values of the dimensionless diffusion parameter and saturation parameter. Furthermore, the solutions are also in good agreement with previous solutions by existing authors.


2021 ◽  
Vol 15 ◽  
pp. 174830262199958
Author(s):  
Colin L Defreitas ◽  
Steve J Kane

This paper proposes a numerical approach to the solution of the Fisher-KPP reaction-diffusion equation in which the space variable is developed using a purely finite difference scheme and the time development is obtained using a hybrid Laplace Transform Finite Difference Method (LTFDM). The travelling wave solutions usually associated with the Fisher-KPP equation are, in general, not deemed suitable for treatment using Fourier or Laplace transform numerical methods. However, we were able to obtain accurate results when some degree of time discretisation is inbuilt into the process. While this means that the advantage of using the Laplace transform to obtain solutions for any time t is not fully exploited, the method does allow for considerably larger time steps than is otherwise possible for finite-difference methods.


2016 ◽  
Vol 16 (4) ◽  
pp. 609-631 ◽  
Author(s):  
Immanuel Anjam ◽  
Dirk Pauly

AbstractThe results of this contribution are derived in the framework of functional type a posteriori error estimates. The error is measured in a combined norm which takes into account both the primal and dual variables denoted by x and y, respectively. Our first main result is an error equality for all equations of the class ${\mathrm{A}^{*}\mathrm{A}x+x=f}$ or in mixed formulation ${\mathrm{A}^{*}y+x=f}$, ${\mathrm{A}x=y}$, where the exact solution $(x,y)$ is in $D(\mathrm{A})\times D(\mathrm{A}^{*})$. Here ${\mathrm{A}}$ is a linear, densely defined and closed (usually a differential) operator and ${\mathrm{A}^{*}}$ its adjoint. In this paper we deal with very conforming mixed approximations, i.e., we assume that the approximation ${(\tilde{x},\tilde{y})}$ belongs to ${D(\mathrm{A})\times D(\mathrm{A}^{*})}$. In order to obtain the exact global error value of this approximation one only needs the problem data and the mixed approximation itself, i.e., we have the equality$\lvert x-\tilde{x}\rvert^{2}+\lvert\mathrm{A}(x-\tilde{x})\rvert^{2}+\lvert y-% \tilde{y}\rvert^{2}+\lvert\mathrm{A}^{*}(y-\tilde{y})\rvert^{2}=\mathcal{M}(% \tilde{x},\tilde{y}),$where ${\mathcal{M}(\tilde{x},\tilde{y}):=\lvert f-\tilde{x}-\mathrm{A}^{*}\tilde{y}% \rvert^{2}+\lvert\tilde{y}-\mathrm{A}\tilde{x}\rvert^{2}}$ contains only known data. Our second main result is an error estimate for all equations of the class ${\mathrm{A}^{*}\mathrm{A}x+ix=f}$ or in mixed formulation ${\mathrm{A}^{*}y+ix=f}$, ${\mathrm{A}x=y}$, where i is the imaginary unit. For this problem we have the two-sided estimate$\frac{\sqrt{2}}{\sqrt{2}+1}\mathcal{M}_{i}(\tilde{x},\tilde{y})\leq\lvert x-% \tilde{x}\rvert^{2}+\lvert\mathrm{A}(x-\tilde{x})\rvert^{2}+\lvert y-\tilde{y}% \rvert^{2}+\lvert\mathrm{A}^{*}(y-\tilde{y})\rvert^{2}\leq\frac{\sqrt{2}}{% \sqrt{2}-1}\mathcal{M}_{i}(\tilde{x},\tilde{y}),$where ${\mathcal{M}_{i}(\tilde{x},\tilde{y}):=\lvert f-i\tilde{x}-\mathrm{A}^{*}% \tilde{y}\rvert^{2}+\lvert\tilde{y}-\mathrm{A}\tilde{x}\rvert^{2}}$ contains only known data. We will point out a motivation for the study of the latter problems by time discretizations or time-harmonic ansatz of linear partial differential equations and we will present an extensive list of applications including the reaction-diffusion problem and the eddy current problem.


Sign in / Sign up

Export Citation Format

Share Document