Abstract: In this paper, we present a simple effective numerical geometric method to estimate the divergence of a vector field over a curved surface. The conservation law is an important principle in physics and mathematics. However, many well-known numerical methods for solving diffusion equations do not obey conservation laws. Our presented method in this paper combines the divergence theorem with a generalized finite difference method and obeys the conservation law on discrete closed surfaces. We use the similar method to solve the Cahn-Hilliard equations on evolving spherical surfaces and observe stability results in our numerical simulations.
Abstract: The diffusion-reaction equations are important Partial Differential Equations in mathematical biology, material science, physics, and so on. However, finding efficient numerical methods for diffusion-reaction systems on curved surfaces is still an important and difficult problem. The purpose of this paper is to present a convergent geometric method for solving the reaction-diffusion equations on closed surfaces by an O(r)-LTL configuration method. The O(r)-LTL configuration method combining the local tangential lifting technique and configuration equations is an effective method to estimate differential quantities on curved surfaces. Since estimating the Laplace-Beltrami operator is an important task for solving the reaction-diffusion equations on surfaces, we use the local tangential lifting method and a generalized finite difference method to approximate the Laplace-Beltrami operators and we solve this reaction-diffusion system on closed surfaces. Our method is not only conceptually simple, but also easy to implement.
Abstract: In this paper, a backward semi-Lagrangian scheme
combined with the second-order backward difference formula
is designed to calculate the numerical solutions of nonlinear
advection-diffusion equations. The primary aims of this paper are
to remove any iteration process and to get an efficient algorithm
with the convergence order of accuracy 2 in time. In order to achieve
these objects, we use the second-order central finite difference and the
B-spline approximations of degree 2 and 3 in order to approximate
the diffusion term and the spatial discretization, respectively. For the
temporal discretization, the second order backward difference formula
is applied. To calculate the numerical solution of the starting point
of the characteristic curves, we use the error correction methodology
developed by the authors recently. The proposed algorithm turns out
to be completely iteration free, which resolves the main weakness
of the conventional backward semi-Lagrangian method. Also, the
adaptability of the proposed method is indicated by numerical
simulations for Burgers’ equations. Throughout these numerical
simulations, it is shown that the numerical results is in good
agreement with the analytic solution and the present scheme offer
better accuracy in comparison with other existing numerical schemes.
Abstract: In this paper, the difference between the Alternating
Direction Method (ADM) and the Non-Splitting Method (NSM) is
investigated, while both methods applied to the simulations for 2-D
multimaterial radiation diffusion issues. Although the ADM have the
same accuracy orders with the NSM on the uniform meshes, the
accuracy of ADM will decrease on the distorted meshes or the
boundary of domain. Numerical experiments are carried out to
confirm the theoretical predication.
Abstract: In this paper, we develop an accurate and efficient Haar wavelet method for well-known FitzHugh-Nagumo equation. The proposed scheme can be used to a wide class of nonlinear reaction-diffusion equations. The power of this manageable method is confirmed. Moreover the use of Haar wavelets is found to be accurate, simple, fast, flexible, convenient, small computation costs and computationally attractive.
Abstract: In this paper we present modeling and simulation for
physical vapor deposition for metallic bipolar plates. In the models
we discuss the application of different models to simulate the
transport of chemical reactions of the gas species in the gas chamber.
The so called sputter process is an extremely sensitive process to
deposit thin layers to metallic plates. We have taken into account
lower order models to obtain first results with respect to the gas
fluxes and the kinetics in the chamber.
The model equations can be treated analytically in some
circumstances and complicated multi-dimensional models are solved
numerically with a software-package (UG unstructed grids, see [1]).
Because of multi-scaling and multi-physical behavior of the models,
we discuss adapted schemes to solve more accurate in the different
domains and scales. The results are discussed with physical
experiments to give a valid model for the assumed growth of thin
layers.
Abstract: Calcium is a vital second messenger used in signal transduction. Calcium controls secretion, cell movement, muscular contraction, cell differentiation, ciliary beating and so on. Two theories have been used to simplify the system of reaction-diffusion equations of calcium into a single equation. One is excess buffer approximation (EBA) which assumes that mobile buffer is present in excess and cannot be saturated. The other is rapid buffer approximation (RBA), which assumes that calcium binding to buffer is rapid compared to calcium diffusion rate. In the present work, attempt has been made to develop a model for calcium diffusion under excess buffer approximation in neuron cells. This model incorporates the effect of [Na+] influx on [Ca2+] diffusion,variable calcium and sodium sources, sodium-calcium exchange protein, Sarcolemmal Calcium ATPase pump, sodium and calcium channels. The proposed mathematical model leads to a system of partial differential equations which have been solved numerically using Forward Time Centered Space (FTCS) approach. The numerical results have been used to study the relationships among different types of parameters such as buffer concentration, association rate, calcium permeability.