Search results
1 – 10 of 910Sanjay Komala Sheshachala and Ramon Codina
This paper aims to present a finite element formulation to approximate systems of reaction–diffusion–advection equations, focusing on cases with nonlinear reaction. The…
Abstract
Purpose
This paper aims to present a finite element formulation to approximate systems of reaction–diffusion–advection equations, focusing on cases with nonlinear reaction. The formulation is based on the orthogonal sub-grid scale approach, with some simplifications that allow one to stabilize only the convective term, which is the source of potential instabilities. The space approximation is combined with finite difference time integration and a Newton–Raphson linearization of the reactive term. Some numerical examples show the accuracy of the resulting formulation. Applications using classical nonlinear reaction models in population dynamics are also provided, showing the robustness of the approach proposed.
Design/methodology/approach
A stabilized finite element method for advection–diffusion–reaction equations to the problem on nonlinear reaction is adapted. The formulation designed has been implemented in a computer code. Numerical examples are run to show the accuracy and robustness of the formulation.
Findings
The stabilized finite element method from which the authors depart can be adapted to problems with nonlinear reaction. The resulting method is very robust and accurate. The framework developed is applicable to several problems of interest by themselves, such as the predator–prey model.
Originality/value
A stabilized finite element method to problems with nonlinear reaction has been extended. Original contributions are the design of the stabilization parameters and the linearization of the problem. The application examples, apart from demonstrating the validity of the numerical model, help to get insight in the system of nonlinear equations being solved.
Details
Keywords
The purpose of this paper is to develop a new finite-volume method of lines for one-dimensional reaction-diffusion equations that provides piece-wise analytical solutions in space…
Abstract
Purpose
The purpose of this paper is to develop a new finite-volume method of lines for one-dimensional reaction-diffusion equations that provides piece-wise analytical solutions in space and is conservative, compare it with other finite-difference discretizations and assess the effects of the nonlinear diffusion coefficient on wave propagation.
Design/methodology/approach
A conservative, finite-volume method of lines based on piecewise integration of the diffusion operator that provides a globally continuous approximate solution and is second-order accurate is presented. Numerical experiments that assess the accuracy of the method and the time required to achieve steady state, and the effects of the nonlinear diffusion coefficients on wave propagation and boundary values are reported.
Findings
The finite-volume method of lines presented here involves the nodal values and their first-order time derivatives at three adjacent grid points, is linearly stable for a first-order accurate Euler’s backward discretization of the time derivative and has a smaller amplification factor than a second-order accurate three-point centered discretization of the second-order spatial derivative. For a system of two nonlinearly-coupled, one-dimensional reaction-diffusion equations, the amplitude, speed and separation of wave fronts are found to be strong functions of the dependence of the nonlinear diffusion coefficients on the concentration and temperature.
Originality/value
A new finite-volume method of lines for one-dimensional reaction-diffusion equations based on piecewise analytical integration of the diffusion operator and the continuity of the dependent variables and their fluxes at the cell boundaries is presented. The method may be used to study heat and mass transfer in layered media.
Details
Keywords
The purpose of this paper is to develop a new transversal method of lines for one-dimensional reaction–diffusion equations that is conservative and provides piecewise–analytical…
Abstract
Purpose
The purpose of this paper is to develop a new transversal method of lines for one-dimensional reaction–diffusion equations that is conservative and provides piecewise–analytical solutions in space, analyze its truncation errors and linear stability, compare it with other finite-difference discretizations and assess the effects of the nonlinear diffusion coefficients, reaction rate terms and initial conditions on wave propagation and merging.
Design/methodology/approach
A conservative, transversal method of lines based on the discretization of time and piecewise analytical integration of the resulting two-point boundary-value problems subject to the continuity of the dependent variables and their fluxes at the control-volume boundaries, is presented. The method provides three-point finite difference expressions for the nodal values and continuous solutions in space, and its accuracy has been determined first analytically and then assessed in numerical experiments of reaction-diffusion problems, which exhibit interior and/or boundary layers.
Findings
The transversal method of lines presented here results in three-point finite difference equations for the nodal values, treats the diffusion terms implicitly and is unconditionally stable if the reaction terms are treated implicitly. The method is very accurate for problems with the interior and/or boundary layers. For a system of two nonlinearly-coupled, one-dimensional reaction–diffusion equations, the formation, propagation and merging of reactive fronts have been found to be strong function of the diffusion coefficients and reaction rates. For asymmetric ignition, it has been found that, after front merging, the temperature and concentration profiles are almost independent of the ignition conditions.
Originality/value
A new, conservative, transversal method of lines that treats the diffusion terms implicitly and provides piecewise exponential solutions in space without the need for interpolation is presented and applied to someone.
Details
Keywords
Ali Şahin, İdris Dağ and Bülent Saka
This paper seeks to develop an efficient B‐spline Galerkin scheme for solving the Fisher's equation, which is a nonlinear reaction diffusion equation describing the relation…
Abstract
Purpose
This paper seeks to develop an efficient B‐spline Galerkin scheme for solving the Fisher's equation, which is a nonlinear reaction diffusion equation describing the relation between the diffusion and nonlinear multiplication of a species.
Design/methodology/approach
The solution domain is partitioned into uniform mesh and, using the quartic B‐spline functions, the Galerkin method is applied to the Fisher's equation.
Findings
The method yields stable accurate solutions. Obtained results are acceptable and in unison with some earlier studies.
Originality/value
Using the uniform mesh, quartic B‐spline Galerkin method is employed for finding the numerical solutions of Fisher's equation.
Details
Keywords
In the nonlinear model of reaction–diffusion, the Fitzhugh–Nagumo equation plays a very significant role. This paper aims to generate innovative solitary solutions of the…
Abstract
Purpose
In the nonlinear model of reaction–diffusion, the Fitzhugh–Nagumo equation plays a very significant role. This paper aims to generate innovative solitary solutions of the Fitzhugh–Nagumo equation through the use of variational formulation.
Design/methodology/approach
The partial differential equation of Fitzhugh–Nagumo is modified by the appropriate wave transforms into a dimensionless nonlinear ordinary differential equation, which is solved by a semi-inverse variational method.
Findings
This paper uses a variational approach to the Fitzhugh–Nagumo equation developing new solitary solutions. The condition for the continuation of new solitary solutions has been met. In addition, this paper sets out the Fitzhugh–Nagumo equation fractal model and its variational principle. The findings of the solitary solutions have shown that the suggested method is very reliable and efficient. The suggested algorithm is very effective and is almost ideal for use in such problems.
Originality/value
The Fitzhugh–Nagumo equation is an important nonlinear equation for reaction–diffusion and is typically used for modeling nerve impulses transmission. The Fitzhugh–Nagumo equation is reduced to the real Newell–Whitehead equation if β = −1. This study provides researchers with an extremely useful source of information in this area.
Details
Keywords
J.I. Ramos and Carmen María García López
The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the…
Abstract
Purpose
The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.
Design/methodology/approach
An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.
Findings
The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.
Originality/value
The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.
Details
Keywords
Mehdi Dehghan, Jalil Manafian Heris and Abbas Saadatmandi
The purpose of this paper is to use He's Exp‐function method (EFM) to construct solitary and soliton solutions of the nonlinear evolution equation.
Abstract
Purpose
The purpose of this paper is to use He's Exp‐function method (EFM) to construct solitary and soliton solutions of the nonlinear evolution equation.
Design/methodology/approach
This technique is straightforward and simple to use and is a powerful method to overcome some difficulties in the nonlinear problems.
Findings
This method is developed for searching exact traveling wave solutions of the nonlinear partial differential equations. The EFM presents a wider applicability for handling nonlinear wave equations.
Originality/value
The paper shows that EFM, with the help of symbolic computation, provides a straightforward and powerful mathematical tool for solving nonlinear evolution equations. Application of EFM to Fitzhugh‐Nagumo equation illustrates its effectiveness.
Details
Keywords
The purpose of this study is to introduce the reproducing kernel algorithm for treating classes of time-fractional partial differential equations subject to Robin boundary…
Abstract
Purpose
The purpose of this study is to introduce the reproducing kernel algorithm for treating classes of time-fractional partial differential equations subject to Robin boundary conditions with parameters derivative arising in fluid flows, fluid dynamics, groundwater hydrology, conservation of energy, heat conduction and electric circuit.
Design/methodology/approach
The method provides appropriate representation of the solutions in convergent series formula with accurately computable components. This representation is given in the W(Ω) and H(Ω) inner product spaces, while the computation of the required grid points relies on the R(y,s) (x, t) and r(y,s) (x, t) reproducing kernel functions.
Findings
Numerical simulation with different order derivatives degree is done including linear and nonlinear terms that are acquired by interrupting the n-term of the exact solutions. Computational results showed that the proposed algorithm is competitive in terms of the quality of the solutions found and is very valid for solving such time-fractional models.
Research limitations/implications
Future work includes the application of the reproducing kernel algorithm to highly nonlinear time-fractional partial differential equations such as those arising in single and multiphase flows. The results will be published in forthcoming papers.
Practical implications
The study included a description of fundamental reproducing kernel algorithm and the concepts of convergence, and error behavior for the reproducing kernel algorithm solvers. Results obtained by the proposed algorithm are found to outperform in terms of accuracy, generality and applicability.
Social implications
Developing analytical and numerical methods for the solutions of time-fractional partial differential equations is a very important task owing to their practical interest.
Originality/value
This study, for the first time, presents reproducing kernel algorithm for obtaining the numerical solutions of some certain classes of Robin time-fractional partial differential equations. An efficient construction is provided to obtain the numerical solutions for the equations, along with an existence proof of the exact solutions based upon the reproducing kernel theory.
Details
Keywords
Mostafa Abbaszadeh, AliReza Bagheri Salec and Shurooq Kamel Abd Al-Khafaji
The space fractional PDEs (SFPDEs) play an important role in the fractional calculus field. Proposing a high-order, stable and flexible numerical procedure for solving SFPDEs is…
Abstract
Purpose
The space fractional PDEs (SFPDEs) play an important role in the fractional calculus field. Proposing a high-order, stable and flexible numerical procedure for solving SFPDEs is the main aim of most researchers. This paper devotes to developing a novel spectral algorithm to solve the FitzHugh–Nagumo models with space fractional derivatives.
Design/methodology/approach
The fractional derivative is defined based upon the Riesz derivative. First, a second-order finite difference formulation is used to approximate the time derivative. Then, the Jacobi spectral collocation method is employed to discrete the spatial variables. On the other hand, authors assume that the approximate solution is a linear combination of special polynomials which are obtained from the Jacobi polynomials, and also there exists Riesz fractional derivative based on the Jacobi polynomials. Also, a reduced order plan, such as proper orthogonal decomposition (POD) method, has been utilized.
Findings
A fast high-order numerical method to decrease the elapsed CPU time has been constructed for solving systems of space fractional PDEs.
Originality/value
The spectral collocation method is combined with the POD idea to solve the system of space-fractional PDEs. The numerical results are acceptable and efficient for the main mathematical model.
Details
Keywords
Nizhnik–Novikov–Veselov system (NNVS) is a well-known isotropic extension of the Lax (1 + 1) dimensional Korteweg-deVries equation that is also used as a paradigm for an…
Abstract
Purpose
Nizhnik–Novikov–Veselov system (NNVS) is a well-known isotropic extension of the Lax (1 + 1) dimensional Korteweg-deVries equation that is also used as a paradigm for an incompressible fluid. The purpose of this paper is to present a fractal model of the NNVS based on the Hausdorff fractal derivative fundamental concept.
Design/methodology/approach
A two-scale transformation is used to convert the proposed fractal model into regular NNVS. The variational strategy of well-known Chinese scientist Prof. Ji Huan He is used to generate bright and exponential soliton solutions for the proposed fractal system.
Findings
The NNV fractal model and its variational principle are introduced in this paper. Solitons are created with a variety of restriction interactions that must all be applied equally. Finally, the three-dimensional diagrams are displayed using an appropriate range of physical parameters. The results of the solitary solutions demonstrated that the suggested method is very accurate and effective. The proposed methodology is extremely useful and nearly preferable for use in such problems.
Practical implications
The research study of the soliton theory has already played a pioneering role in modern nonlinear science. It is widely used in many natural sciences, including communication, biology, chemistry and mathematics, as well as almost all branches of physics, including nonlinear optics, plasma physics, fluid dynamics, condensed matter physics and field theory, among others. As a result, while constructing possible soliton solutions to a nonlinear NNV model arising from the field of an incompressible fluid is a popular topic, solving nonlinear fluid mechanics problems is significantly more difficult than solving linear ones.
Originality/value
To the best of the authors’ knowledge, for the first time in the literature, this study presents Prof. Ji Huan He's variational algorithm for finding and studying solitary solutions of the fractal NNV model. The reported solutions are novel and present a valuable addition to the literature in soliton theory.
Details