Home A modified cubic B-spline differential quadrature method for three-dimensional non-linear diffusion equations
Article Open Access

A modified cubic B-spline differential quadrature method for three-dimensional non-linear diffusion equations

  • Sumita Dahiya EMAIL logo and Ramesh Chandra Mittal
Published/Copyright: July 3, 2017

Abstract

This paper employs a differential quadrature scheme for solving non-linear partial differential equations. Differential quadrature method (DQM), along with modified cubic B-spline basis, has been adopted to deal with three-dimensional non-linear Brusselator system, enzyme kinetics of Michaelis-Menten type problem and Burgers’ equation. The method has been tested efficiently to three-dimensional equations. Simple algorithm and minimal computational efforts are two of the major achievements of the scheme. Moreover, this methodology produces numerical solutions not only at the knot points but also at every point in the domain under consideration. Stability analysis has been done. The scheme provides convergent approximate solutions and handles different cases and is particularly beneficial to higher dimensional non-linear PDEs with irregularities in initial data or initial-boundary conditions that are discontinuous in nature, because of its capability of damping specious oscillations induced by high frequency components of solutions.

1 Introduction

In this work, we have extended an earlier approach [1] to obtain the numerical solutions of three-dimensional non-linear reaction-diffusion equations and three-dimensional Burgers’ equation. Reaction-diffusion models appear very frequently while investigating biological and chemical phenomena and time-dependent PDEs are adopted to model them naturally i.e. reaction-diffusion equations. The non-linear reaction-diffusion system is made up of two different terms: (a) diffusion term, involving the random motion of every particular component because of turbulent nature of the governing field flow and (b) reaction term, which describes the interactions between biological and physical species that are involved [2]. Over the past decades, reaction-diffusion systems have attained enormous interest from researchers because of their complexity and behavioral dynamics such as periodicity in solution, formation of patterns, and chaos [35]. We have discussed some well known and reaction-diffusion models appearing very frequently: the Brusselator model [4] and the Michaelis-Menten model [6]. A number of authors have emphasized the significance of oscillations in biochemical systems, e.g. Turing [7] observed that a stable spatial arrangement can possibly be obtained by coupling the diffusion process with certain reactions leading to the idea of morphogenesis.

The development and analysis of the behavior of a non-linear oscillator related to chemical system

Bin=⇒X1,Ain+X1=⇒X2+D1,2X1+X2=⇒3X1,X1=⇒D2,(1)

was developed and analysed by the Brussels schools [8]. Here Ain, Bin, D1 and D2 denote the input and output chemicals, respectively. X1 and X2 are intermediates. The system (1) consists of only two dependent variables and admits limit cycle oscillations. The trimolecular reaction step in (1) appears in manifold coupling between modes in laser and plasma physics [9, 10], during enzymatic reactions, ozone formation by atomic oxygen by triple collision. The three-dimensional Brusselator system is written as

ut=A+u2v(B+1)u+α2u,(2a)
vt=Buu2v+α2u,(x,t)Ω×(0,T](2b)

together with initial conditions

(u(x,0),v(x,0))=(g1(x),g2(x)),xΩ(3)

and no flux conditions on boundary.

Here, u(x, t) = u(x, y, z, t) and v(x, t) = v(x, y, z, t) represents the concentrations of two reactants in the dimensionless form. A and B are constant concentrations of the two reactants, α is coefficient of diffusivity. For Brusselator system (2), if 1 − A + B2 > 0, then solution reaches the point of equilibria (B, A/B) in steady state, given the values of coefficient of diffusivity α are kept very small.

Michaelis-Menten model is one of the best models known for enzyme kinetics. It involves an enzyme binding to a substrate S that results in the formation of an enzyme-substrate complex that is reduced in a product P and enzyme. The model transforms into an equation that describes the rate of enzymatic reactions by setting a relationship between reaction rate v to the concentration of substrate S. The formula is given by

v=d[P]dt=VmaxSKM+S.

Here, Vmax represents the highest rate of the system at the highest substrate concentrations. KM, known as the Michaelis constant, gives the substrate concentration at which rate of reaction is half of Vmax[12]. Biochemical reactions that involve only one substrate are usually taken to follow kinetics of Michaelis-Menten, with no consideration to underlying assumptions of the model. The three-dimensional counterpart of Michaelis-Menten type one-dimensional problem discussed in [11] is given and reads as:

ut=d(2ux2+2uy2+2uz2)u1u,0<x,y,z<1,t>0,(4)

together with initial condition

u(x,y,z,0)=f1(x,y,z),0x,y,z1.(5)

and homogeneous boundary conditions.

In the present scenario, reaction-diffusion systems have gained significant concern as a model for formation of patterns. The stripes, hexagons, targets, spirals, fronts and dissipative soliton patterns are seen in many types of reaction-diffusion systems despite huge differences e.g. in the local reaction term. The reaction-diffusion processes are found to be significant bases for processes correlated to morphogenesis in biology [13]. Wound healing [14], tumor growth [15], spread of epidemics [16], ecological invasions [17] etc. are some of fields in which reaction-diffusion systems are applicable.

Burgers’ equation has been extensively studied to portray many types of phenomena like mathematical model of turbulence [18], traveling of shock waves in viscous fluid [19], etc. One can take it as Navier-Stokes equation [20, 21] in a simplified version, because of the form of non-linear convective term and the appearance of the viscosity term.

Therefore, study of Burgers’ equation is useful in understanding the convection-diffusion phenomena. Moreover, the study relates to many problems in physics such as acoustic transmission [22], traffic and aerofoil flow theory [23], gravity effect on sedimentation of two types of particles in fluid suspensions [24] as well as a prerequisite to the Navier-Stokes equations [25].

The Burgers’ equation in one dimension is likewise to the Navier-Stokes equations in one dimension barring the stress term. A significant amount of research has been done in the recent years to numerically solve the Burgers’ equation efficiently for varying kinematic viscosity values. Many numerical approaches to solve Burgers’ equation such as Galerkin finite element method [26], finite volume method [27], iterative methods [28], automatic differentiation method [29] have been introduced in literature.

The three-dimensional Burgers’ equation is

ut=uuxuuyuuz+1R(2ux2+2uy2+2zt2),(6)

with prescribed initial and boundary conditions.

2 Differential quadrature method

Differential quadrature method [30] is a numerical approach to approximate numerically the solutions of ODEs and PDEs. The derivatives of the unknown functions are approximated numerically using DQM at some predefined discrete points as a weighted sum of functional values in the complete domain. Computational inexpensiveness, accuracy and ease of application makes the technique more preferable. By employing differential quadrature method, we are capable of approximating the derivatives of functions in space at any of the given knots with the weighted linear sum of function values at all knots in the domain under consideration. To apply DQM in three dimensions, firstly we discretize the domain Ω = {(x,y,z) : 0 ≤ x, y, z ≤ 1} as D1 = {(xi, yj, zk), i = 1, 2,…N; j = 1, 2,…, M; k = 1, 2,…L} by taking step length Δx = xixi − 1 in x-axis direction, Δ y = yjyj − 1 in y-axis direction and Δz = zkzk − 1 in z-axis direction. Approximation of the partial derivative of order one of functional u(x, y, z, t) along x-axis, at point xi, keeping yj and zk fixed, is given as follows

ux(xi,yj,zk,t)=r=1Nair(1)u(xr,yj,zk,t),i=1,2,...,N(7)

where air(1) represent the the weighting coefficients to be determined to approximate the partial derivatives of order one in the direction of x. Various approaches are available in literature that can compute these weighting coefficients, e.g. Shu’s approach [31], Quan and Chang’s approach [32, 33], etc. Korkmaz and Daĝ [34, 35] used DQM based on cosine expansion and sinc DQM for many non-linear PDEs. Lagrange interpolation and sine-cosine expansion based quadrature methods has been very frequently used in literature. A differential quadrature method based on polynomials have been suggested by Mittal et al. [36] for numerical approximation of the the solutions of PDEs both linear and non-linear in nature. The idea of iterative DQM has been introduced by Tomasiello [3740]. In her work, iterative DQM and the stability of DQ solutions have been discussed. Korkmaz [41] in his work employed DQM for solution of KdV equation. Arora and Singh [42] worked on one-dimensional Burgers’ equation, whereas Jiwari and Yuan [43] worked on two-dimensional Brusselator model employing DQM. Bashan et al. [44] used Quintic B-splines DQM to approximate KdVB equation, whereas Korkmaz and Daĝ [45] employed both quartic and quintic B-splines for solving advection-diffusion equations. DQM based on Hermite radial basis function has been used by Krowick [46] to deal with higher order equations and DQM based on quasi-interpolation has been employed by Barrera et al. [47]. Shukla et al. [4850] worked extensively on DQM to solve two and three-dimensional PDEs. In the present work, weighting coefficients of DQM are determined using a methodology that employs modified cubic B-spline functions. In this work, we have used cubic B-splines as basis functions and weighting coefficients are computed accordingly.

2.1 Modified cubic B-spline functions

Computation of weighting coefficients is a key procedure which is of paramount importance. In this work, the authors have computed the weighting coefficients air(1),ajr(1) and akr(1) by adopting modified cubic B-spline functions. These functions can be written as:

φm(x)=1h3(xxm2)3x[xm2,xm1)(xxm2)34(xxm1)3x[xm1,xm)(xm+2x)34(xm+1x)3x[xm,xm+1)(xm+2x)3x[xm+1,xm+2)0elsewhere(8)

Here, {φ0(x), φ1(x),…, φN + 1(x)} defines a basis over the computational domain [a, b]. Table 1 comprises the values of original cubic B-splines and their derivatives at the knots. The B-spline functions are modified in a particular manner in which the diagonal dominance of the resultant matrix can be preserved. Jain [51] defines the modified cubic B-spline basis functions at the knot points as:

ϕ1(x)=φ1(x)+2φ0(x),ϕ2(x)=φ2(x)φ0(x),ϕl(x)=φl(x),l=3,4,....,N2,ϕN1=φN1(x)φN+1(x)ϕN(x)=φN(x)+2φN+1(x)(9)
Table 1

Values of φm(x) and its derivatives at the nodal points

xm − 2xm − 1xmxm+1xm + 2
φm(x)01410
φm(x)03/h03/h0
φm(x)06/h2−12/h26/h20

The function ϕl(x), l = 1, 2,…,N again forms a basis over the interval [a, b].

2.2 Determination of weighting coefficients

Korkmaz and Daĝ [52] used classic cubic B-spline functions and Korkmaz and Akmaz [53] used exponential cubic B-spline functions to compute the weighting coefficients. In their work, the resultant algebraic equation systems has (N + 2) equations in (N + 4) unknowns. Therefore, two more equations has to be added in the system. But, in the present study this is not the case. Determining the weighting coefficients air(1) and using modified B-spline functions ϕm(xi), m = 1, 2,…,N in Eq. (7), gives

ϕm(xi)=k=1Naikϕm(xk),i=1,2,....,N

To compute the weighting coefficients air(1), we fix yaxis and zaxis in Eq. (7) as done in [1]. “Thomas algorithm” is adopted to solve the resultant tri-diagonal system and the weighting coefficients a11(1),a12(1),.....,a1N(1) are computed. Similarly, the weighting coefficients for each knot i = 2, 3,…, N can be computed. With these coefficients, we can approximate partial derivatives of order one in the direction of xaxis. The following recurrence relation is used to compute derivatives of higher order [31].

air(p)=r[air(1)aii(p1)air(p1)xixr],foriri,r=1,2,....,N;p=2,3,....,N1(10)
aii(p)=r=1,riNair(p),fori=r(11)

Here air(p1) represents the weighting coefficients of partial derivatives of order (p − 1) w.r.t. x and air(p) represents the weighting coefficients of derivatives of order (p) w.r.t. x. Following the same methodology, we can determine the weighting coefficients ajr(1) and akr(1) of partial derivatives of order one in the direction of yaxis and zaxis using cubic B-spline functions in modified form. Recurrence formula can be successfully applied to calculate the derivatives of order two and higher.

2.3 Partition and denotations

For time-step discretization, τ is time increment and ∃ two positive integers N1 and n such that τ = T/N1, tn = nt, 0 ≤ nN1. For spatial approximations, let hx, hy and hz be x−, y− and z− dimensional spatial increments, such that hx = (b1a1)/N, hy = (b2a2)/M, hz = (b3a3)/L for positive integers N, M, and L. In this way the set Ωh = {(xi, yj, zk)|1 ≤ iN − 1, 1 ≤ jM − 1, 1 ≤ kL − 1} is composed of the spatial nodes (xi, yj, zk) = (a1 + ihx, a2 + jhy, a3 + khz), which are denoted by xijk.

3 Three-dimensional reaction-diffusion systems

3.1 Brusselator system

For numerical treatment of Brusselator system rewriting Eq. (2) in the first order system, discretizing the derivatives in space using MCB-DQM, we achieve the following system of non-linear ODEs

duijkdt=A+uijk2vijk(B+1)uijk+μ(r=1Nair(2)u(xr,yj,zk,t)+r=1Majr(2)u(xi,yr,zk,t)+r=1Lakr(2)u(xi,yj,zr,t))(12a)
dvijkdt=Buijkuijk2vijk+α(k=1Nair(2)v(xr,yj,zk,t)+r=1Majr(2)v(xi,yr,zk,t)+r=1Lakr(2)v(xi,yj,zr,t)),(12b)
ω(xi,yj,zk)/0<xi,yj,zk<1,t>0(12c)

subject to initial conditions

u(xi,yj,zk,0)=u0(xi,yj,zk)v(xi,yj,zk,0)=v0(xi,yj,zk)(13)

and no flux boundary conditions.

3.2 Michaelis-Menten type reaction-diffusion equation

For numerical treatment of Michaelis-Menten system rewriting Eq. (4) in the first order system, discretizing the spatial derivatives by using MCB-DQM, we obtain the following system of non-linear ODEs

duijkdt=d(r=1Nair(2)u(xr,yj,zk,t)+r=1Majr(2)u(xi,yr,zk,t)r=1Lakr(2)u(xi,yj,zr,t))uijk1uijk,0<x,y,z<1,t>0,(14)

together with homogeneous Dirichlet’s boundary conditions and suitable initial conditions.

4 Discretization of three-dimensional non-linear Burgers’ equation and stability analysis

For numerical treatment of Burgers’ equation rewriting Eq. (6) in the first order system, discreting the spatial derivatives using MCB-DQM, we obtain the following system of non-linear ODEs

duijkdt=uijk(r=1Nair(1)u(xr,yj,zk,t)+r=1Majr(1)u(xi,yr,zk,t)+r=1Lakr(1)u(xi,yj,zr,t))+1R(r=1Nair(2)u(xr,yj,zk,t)+r=1Majr(2)u(xi,yr,zk,t)+r=1Lakr(2)u(xi,yj,zr,t))(15)

with given initial and boundary conditions. Von Neumann stability analysis can be easily done when one is working with the classical finite difference. As with the case here, Von Neumann stability criterion can not be used with differential quadrature method discretized systems. For this purpose, matrix stability or energy stability methods have been extensively discussed in literature [3840].

Let us consider a problem that is time-dependent along with a set of suitable initial and boundary conditions, given as

Vt=l(V)(16)

Here, l is the operator in space. Discretizing Eq. (16), we attain a system of ODEs:

dvdt=[Dc]v+b¯(17)

where v represents unknown functional at knots except at boundaries, Dc denotes the coefficient matrix and b is the non-homogeneous part. Stability conditions of the computational method for solving the Eq. (16) are dependent on the stability of ODEs in Eq. (17). The matrix Dc plays an inevitable part as its eigenvalues directly determine the exact solutions. Therefore, the whole scenario can be well described using the eigenvalues of the coefficient matrix. When the real part of all the eigenvalues of Dc, R(λi) ≤ 0, it shows the stable exact solutions as t → ∞. The matrix Dc can be determined as [Dij]=α(air(1)+ajr(1)+akr(1))+1R(air(2)+ajr(2)+akr(2)), where α = uijk.

The stable solution for t tending to infinity requires

  1. for real eigenvalues −2.78 < Δt. λi < 0,

  2. for eigenvalues having complex part only, −2 2 < Δ t. λi < 2 2,

  3. For complex eigenvalues, Δt. λi should lie in the sector shown in Figure 1.

Figure 1 Stability region
Figure 1

Stability region

For complex eigenvalues, real parts can be small enough positive numbers that lie within a prescribed tolerance limit. [54]. For linearization, we have assumed α = uijk, α constant therefore, the eigenvalues need to be calculated for every problems. As it will induce changes in the entries of the coefficient matrix D. Eigenvalues have been calculated for different grid sizes and have been found in region of stability.

5 Boundary conditions with MCB-DQM

Whenever Dirichlet’s boundary conditions are given, they can be used directly at the boundary without any manipulation. In the case of Neumann or mixed type boundary conditions, they can be discretized adopting MCB-DQM. They are simplified into solutions at boundary points. Dirichlet’s boundary conditions for Michaelis-Menten and Burgers’ equation are given as:

u(a1,j,k)=f1(yj,zk,t),u(b1,j,k)=f2(yj,zk,t)u(i,a2,k)=g1(xi,zk,t),u(i,b2,k)=g2(xi,zk,t)u(i,j,a3)=s1(xi,yj,t),u(i,j,b3)=s1(xi,yj,t),(xi,yj,zk)Ω(18)

Neumann boundary conditions at x = a1 and x = b1 can be approximated as

r=1Na1r(1)u(xr,yj,zk,t)=f1(yj,zk,t),(19)
r=1NaNr(1)u(xr,yj,zk,t)=f2(yj,zk,t),j=1,2,...,M,k=1,2,...,L(20)

Rewritting the above system as:

a11(1)u(1jk)+aN1(1)u(Njk)=f1r=2N1a1r(1)u(xr,yj,zk,t),aN1(1)u(1jk)+aNN(1)u(Njk)=f2r=2N1aNr(1)u(xr,yj,zk,t),j=1,2,...,M,k=1,2,...,L

Solving for u1jk and uNjk we get approximate values at the boundary points x = a1 and x = b1. Similarly, boundary conditions at other boundary points can be calculated. The system of non-linear ODEs can not be solved directly by SSP-RK43 method.

The systems of ordinary differential equations given in system (12), (14) and (15) can be solved using various methods. Using the optimal stability preserving Runge-Kutta method of order three and stage four [SSP-RK43] [55], we solve the system of ODEs. It is not possible to solve this non-linear system of initial-boundary ODEs directly by SSP-RK43 scheme. So, it is dealt by first applying boundary conditions on the system, then a system of non-linear ODEs is obtained. Here, the system in Eqs. (12), (14) and (15) integrated from time t0 to t0 + Δt through the following operations.

v(1)=vm+Δt2L(vm)v(2)=v(1)+Δt2L(v(1))v(3)=23vm+u(2)3+Δt6L(v(2))v(m+1)=v(3)+Δt2L(v(3))(21)

Hence, giving solution v at any specified time level. Here L(vm) is the right hand side of resultant equation after space discretization at the initial time.

6 Why SSP-RK method

SSP methods [55, 56] also sometimes termed as TVD (total variation diminishing) time discretization sometimes are used for solving a system of ODEs. In [57], a general Runge-Kutta method of order m is formulated as:

v(0)=vn,v(i)=k=1i1(αi,kv(k)+Δtβi,kL(v(k))),αi,k0,i=1,...,mv(n+1)=v(m).(22)

For the sake of consistent system k=0i1αi,k=1. The hyperbolic conservation law,

ut=f(u)x,

where the spatial derivative, f(u)x is discretized using TVD finite difference with method of lines approximation results into a system of ODEs

ddtu=L(u),

We invariably assume that the discretization in space −L(u) posses the attribute that when it combines with forward first-order Euler time discretization,

un+1=un+ΔtL(un),

then, for a small enough time step governed by the CFL (Courant-Friedrichs-Lewy) condition,

ΔtcΔtFE,

the total variation in discrete solution does not increase in time, where c=mini,kαi,kβi,k, provided βi,kL is replaced by βi,kL whenever βi,k is negative. The prime objective of multistep time discretization or high-order SSP Runge-Kutta is to achieve temporal accuracy of higher order and to maintain strong stability property. In short, idea behind SSP methods is to start with a method-of-lines semi-discretization, which is strongly stable in a certain norm, semi-norm, or convex functional under forward Euler time stepping, with suitably restricted time step Δ t, and then go for a higher order time discretization maintaining strong stability for the given norm, perhaps under a different time step restriction.

7 Test cases

In this section, the numerical approximations are obtained by implementing MCB-DQM for some examples of both reaction-diffusion system and Burgers’ equation. The illustration of the solutions have been done especially for the same parametric values as taken by researchers in their studies for the sake of valid comparison. The computations are carried out with programming in DEV C++. Availability of analytical solutions aids in analyzing the accuracy and reliability of the proposed scheme.

Error analysis is done by using L and L2 errors for the computed solutions defined below:

L(u)max1iM,1jM,1kL|uijkUijk|

and

L2(u)i=1,j=1,k=1i=N,j=M,k=L|uijkUijk|2

where uijk is the approximate solution computed by the proposed scheme and Uijk is the exact solution at the corresponding mesh points (xi, yj, zk) at time t.

7.1 Reaction-diffusion systems

The following section presents some numerical experiments for three-dimensional problems. The calculated results are compared with exact solutions available in literature or compared our result with the work done by researchers in the same field.

Problem 1. Brusselator model The three-dimensional analogue of present one-dimensional Brusselator system [59]:

ut=A+u2v(B+1)u+α2u,(23a)
vt=Buu2v+α2u,xΩ×(0,T](23b)

with no flux boundary conditions and following initial conditions

u(x,y,z,0)=1+sin(2πx)sin(2πy)sin(2πz),v(x,y,z,0)=3,

To compare our results with those given in [59], the parameters are set accordingly. We have taken A = 2.0, B = 1.0, h = 0.01, k = 0.001, and μ = 0.02. Here Δx = Δy = Δz = h and k is the time step. Table 2 shows that at every grid point, the value of u and v approach 2 and 0.5, respectively as t increases. The plots of approximate values of u and v at the grid point (0.5, 0.5) along z−direction vs. time t are shown in Figure 2. It can be seen from Figure 2 that (u, v) ∈ (2.0, 0.5) along z−direction as time t increases. These observations comply with the theoretical results of [30], which claims that whenever the parameters A and B are chosen such that 1 − A + B2 ≥ 0, then the values of u and v approach B and A/B, respectively.

Figure 2 Plot of u(0.3, 0.3, :) and v(0.3, 0.3, :) vs. time with parameters A = 1.0, B = 2.0, h = 0.1, time step k = 0.001 and α
Figure 2

Plot of u(0.3, 0.3, :) and v(0.3, 0.3, :) vs. time with parameters A = 1.0, B = 2.0, h = 0.1, time step k = 0.001 and α

Table 2

Approximated solution for Problem 1 with A = 1., B = 2., μ = 0.02 at different mesh points and different values of t

t(0.2, 0.2, 0.2)(0.8, 0.9, 0.8)
uvuv
1.02.9517700.3261313.0060900.329108
2.02.2920700.4221342.3033500.419887
3.02.0772100.4739462.0788900.473407
4.02.0173700.4932742.0175800.493198
5.02.0034300.4985652.0034500.498556
6.02.0006200.4999532.0001000.499729
7.02.0001000.4221342.3033500.499953
8.02.000020.4999922.000020.499992
9.02.0000000.4999992.0000000.499999
10.02.0000000.5000002.0000000.500000
2.0000000.5000002.0000000.500000

Figure 3 represents the values of u and v for the same step-length, time-step and μ, but for different values of A and B. Here we have taken A = 3. and B = 1. The values of u and v are approximated until time t = 60. The solutions are found to be oscillatory and not convergent to any fixed concentration. The outcome shows compliance with the observation in [59], that the solution does not converge to any fixed concentration, whenever the parameters A and B are chosen such that 1 − A + B2 ≤ 0

Figure 3 Plot of u(0.3, 0.3, :) and v(0.3, 0.3, :) vs. time with parameters A = 3.0, B = 1.0, h = 0.1, time step k = 0.001 and α
Figure 3

Plot of u(0.3, 0.3, :) and v(0.3, 0.3, :) vs. time with parameters A = 3.0, B = 1.0, h = 0.1, time step k = 0.001 and α

Problem 2. [59] In this problem, the three-dimensional Brusselator system is studied. The exact solution is:

u(x,y,z,t)=exp(xyz0.5t)v(x,y,z,t)=exp(x+y+z+0.5t)(24)

Table 3 gives the L2 and L errors for the concentrations u and v, taking the parametric values as B = 1.0, A = 0, α = 1.0/6.0 at different t values. The errors have been calculated with N = M = L = 21.

Table 3

Errors of Problem 2 with A = 0., B = 1., μ = 1/6 at for varying time

tL2 ErrorL Error
uvuv
0.11.90385E-0054.00995E-0052.19122E-0054.723245E-005
1.02.08566E-0059.03154E-0051.47474E-0057.70677E-005
2.01.27933E-0051.49674E-0048.94614E-0061.12714E-004
3.07.76417E-0062.47037E-0045.42659E-0062.09640E-004
4.04.71003E-0064.07456E-0043.29150E-0063.45643E-004
5.02.85838E-0066.71543E-0041.99741E-0065.69587E-004
10.02.34522E-0078.18585E-0041.62877E-0076.94249E-004

Problem 3. In this example, the three-dimensional Michaelis-Menten system is considered with homogeneous Dirichlet’s boundary and initial condition [59]:

u(x,y,z,0)=1,0x,y,z1.

In the problem, there exists a discontinuity between given initial and boundary conditions, that causes solutions to show spurious oscillations. This problem solved with MCB-DQM shows that this scheme has a capability of damping specious oscillations caused by high frequency components in the solution. The computations are carried out for t = 1, 2 and t = 3. Figure 4 and 5 depicts the results with time-step 0.001 and mesh-size N = M = L = 31. The value of d is taken as 0.2.

Figure 4 Numerical solution of Problem 3 with N = 30 and time step 0.001 for final time t = 0.1 and 1
Figure 4

Numerical solution of Problem 3 with N = 30 and time step 0.001 for final time t = 0.1 and 1

Figure 5 Numerical solution of Problem 3 with N = 30 and time step 0.001 for final time t = 2 and 3
Figure 5

Numerical solution of Problem 3 with N = 30 and time step 0.001 for final time t = 2 and 3

Problem 4. Three-dimensional Burgers’ equation: Here we consider the three-dimensional Burgers’ equation over the domain [−0.5, 0.5] × [−0.5, 0.5] × [−0.5, 0.5], initial and boundary conditions have been taken from the exact solution

u(x,y,t)=1tanh((x+y+z2t)R3)

We have obtained numerical solution for t = 0.1, 1, 3, 5 and 10. Table 4 gives the error estimation for varying values of R.

Table 4

Error estimation for problem 4 for different times and different values of R

Rt = 0.1t = 1t = 3t = 5t = 10
R = 1, h = 1/10, Δt = 0.0012.38653E-0061.00550E-0062.29000E-0061.74900E-0072.23800E-010
R = 10, h = 1/10, Δt = 0.0013.54000E-0043.29200E-0054.46100E-0144.48000E-0144.48400E-014
R = 50, h = 1/20, Δt = 0.0102.30000E-0035.15700E-0069.03100E-0159.03900E-0158.97000E-015
R = 100, h = 1/20, Δt = 0.0102.41000E-0021.32000E-0067.79900E-0158.11500E-0158.08900E-015
R = 500, h = 1/20, Δt = 0.0102.24615E-0014.21000E-0034.13300E-0085.53000E-0148.27170E-015

8 Conclusion

In this work, differential quadrature method based on a modification in cubic B-spline functions have been applied for the first time to non-linear three-dimensional reaction-diffusion equations and Burgers’ equation. The method employed in this work elaborates an approach to numerically approximate solutions for three-dimensional equations. The method efficiently provides accurate approximations for variable sets of parametric values. The efficiency and reliability of the approach have been evaluated for test problems under consideration. Different error norms are calculated for varying times or the results present in literature are compared with the numerical solutions achieved by the used method. The approximated solutions and the corresponding error norms are presented in tables for varying time levels and for varying sets of parametric values. We also note that for large values of R, Burgers’ equation behaves merely as hyperbolic partial differential equation, and the problem becomes very difficult to solve as a steep shock-like wave front is developed. The results show that error decreases with the increase in time. We see that the solutions become better with decreasing values of time step. The method gives convergent solutions and for different cases, handles the equations efficiently even for non-linear ones. The obtained approximations are satisfactory and in compliance with the results available in the literature. The strong point of the method is its ease of implementation and less complexity and calculations. The method is economically cheap and can be relied upon. Therefore, it is suggested that the method can be used as an alternative for the class of three-dimensional reaction-diffusion equation. The possibility of extending present approach may be useful to solve partial differential equations of higher order appearing in various applications of engineering and science.

Acknowledgement

The author Ms. Sumita Dahiya would like to thank Ministry of Human Resource and Development with grant number MHR-02-23-200-429 for their funding.

References

[1] Mittal R.C., Dahiya S., et al., Numerical simulation on hyperbolic diffusion equations using modified cubic B-spline differential quadrature methods, Computers and Mathematics with Applications, 2015, 70.5, 737-749.10.1016/j.camwa.2015.04.022Search in Google Scholar

[2] Zhao S., Ovadia J., Liu X., Zhang Y.T., Nie Q., et al., Operator splitting implicit integration factor methods for stiff reaction–diffusion–advection systems, J. Comput. Phys., 2011, 230(15), 5996-6009.10.1016/j.jcp.2011.04.009Search in Google Scholar

[3] Marek M., and Schreiber I., Chaotic behaviour of deterministic dissipative systems (Vol. 1). Cambridge University Press, 1995.Search in Google Scholar

[4] Nicolis G., and Prigogine I., Self-organization in nonequilibrium systems (Vol. 191977). Wiley, New York, 1977.Search in Google Scholar

[5] Ersoy O., Daĝ I., et al., Numerical solutions of the reaction diffusion system by using exponential cubic B-spline collocation algorithms. Open Phys., 2015, 13: 414-427.10.1515/phys-2015-0047Search in Google Scholar

[6] Cherruault Y., Choubane M., Valleton J.M., Vincent J.C., et al., Stability and asymptotic behavior of a numerical solution corresponding to a diffusion-reaction equation solved by a finite difference scheme (Crank-Nicolson), Computers and Mathematics with Applications, 1990, 20(11), 37-46.10.1016/0898-1221(90)90217-8Search in Google Scholar

[7] Turing A.M., The chemical basis of morphogenesis, Philosophical Transactions of the Royal Society of London B: Biological Sciences, 237(641), 1952, 37-72.10.1098/rstb.1952.0012Search in Google Scholar

[8] Lefever R., Nicolis G., et al., Chemical instabilities and sustained oscillations, Journal of theoretical Biology, 1971, 30(2), 267-284.10.1016/0022-5193(71)90054-3Search in Google Scholar

[9] Twizell E.H., Gumel A.B., Cao Q., et al., A second-order scheme for the Brusselator reaction–diffusion system, J. Math. Chem., 1999, 26(4), 297-316.10.1023/A:1019158500612Search in Google Scholar

[10] Adomian G., The diffusion-Brusselator equation. Computers and Mathematics with Applications, 1995, 29(5), 1-3.10.1016/0898-1221(94)00244-FSearch in Google Scholar

[11] Cherruault Y., Choubane M., Valleton J.M., Vincent J.C., et al., Stability and asymptotic behavior of a numerical solution corresponding to a diffusion–reaction equation solved by a finite difference scheme (Crank–Nicolson Scheme), Comput. Math. Appl., 1990, 20 (11), 37-46.10.1016/0898-1221(90)90217-8Search in Google Scholar

[12] http://www.worthington-biochem.com/introbiochem/substrateconc.htmlSearch in Google Scholar

[13] L.G. Harrison, Kinetic Theory of Living Pattern, Cambridge University Press, 1993.10.1017/CBO9780511529726Search in Google Scholar

[14] Sherratt J.A., Murray J.D., et al., Models of epidermal wound healing, Proc. R. Soc. Lond. B, 1990, 241, 29-36..10.1098/rspb.1990.0061Search in Google Scholar

[15] Chaplain M.A.J., Reaction-diffusion pre-patterning and its potential role in tumour invasion, J. Bio. Systems, 1995, 3, 929-936.10.1142/S0218339095000824Search in Google Scholar

[16] Murray J.D., Stanley E.A., Brown D.L., et al., On the Spatial Spread of Rabies among Foxes, Proc. R. Soc. Lond. B, 1986, 229, 111-150.10.1098/rspb.1986.0078Search in Google Scholar

[17] Holmes E.E., Lewis M.A., Banks J.E., Veit R.R., et al., Spatial Interactions and Population Dynamics, Ecology, 1994, 75(1), 17–29.10.2307/1939378Search in Google Scholar

[18] Burgers J.M., A mathematical model illustrating the theory of turbulence, Advances in applied mechanics, 1948, 1, 171-199.10.1016/S0065-2156(08)70100-5Search in Google Scholar

[19] Khater A.H., Temsah R.S., Hassan M.M., et al., A Chebyshev spectral collocation method for solving Burgers’-type equations, J. Comput. Appl. Math., 2008, 222(2), 333-350.10.1016/j.cam.2007.11.007Search in Google Scholar

[20] Ames W.F. (Ed.). Nonlinear partial differential equations in engineering (Vol. 1). Academic press, 1965.10.1016/S0076-5392(08)60004-1Search in Google Scholar

[21] Karpman V.I. Non-Linear Waves in Dispersive Media: International Series of Monographs in Natural Philosophy (Vol. 71). Elsevier, 2016.Search in Google Scholar

[22] Ali A., Haq S., et al., A computational meshfree technique for the numerical solution of the two-dimensional coupled Burgers’ equations, International Journal for Computational Methods in Engineering Science and Mechanics, 2009, 10(5), 406-422.10.1080/15502280903108016Search in Google Scholar

[23] Perko J., Sarler B., et al., Weight function shape parameter optimization in meshless methods for non-uniform grids, Computer Modeling in Engineering and Sciences, 19(1), 2007, 55.Search in Google Scholar

[24] Nee J., Duan J., et al., Limit set of trajectories of the coupled viscous Burgers’ equations. Applied mathematics letters, 1998, 11(1), 57-61.10.1016/S0893-9659(97)00133-XSearch in Google Scholar

[25] Zhang W, Zhang C.H., Guang X., et al., An explicit Chebyshev pseudospectral multigrid method for incompressible Navier–Stokes equations. Computers and Fluids, 2010, 39, 1, 178-188.10.1016/j.compfluid.2009.08.001Search in Google Scholar

[26] Dogan A., A Galerkin finite element approach to Burgers’ equation. Appl. Math. Comput., 2004, 157(2), 331-346.10.1016/j.amc.2003.08.037Search in Google Scholar

[27] Kumar Pany A., Nataraj, N., Singh, S., et al., A new mixed finite element method for Burgers’ equation, Journal of Applied Mathematics and Computing, 2007, 23(1-2), 43-55.10.1007/BF02831957Search in Google Scholar

[28] Golmankhaneh A.K., Khatuni T., Porghoveh N.A., Baleanu D., et al,. Comparison of iterative methods by solving nonlinear Sturm-Liouville, Burgers and Navier-Stokes equations, Central European Journal of Physics, 2012, 10(4), 966-976.10.2478/s11534-012-0038-7Search in Google Scholar

[29] Asaithambi A., et al., Numerical solution of the Burgers’ equation by automatic differentiation, Appl. Math. Comput., 2010, 216, 2700-2708.10.1016/j.amc.2010.03.115Search in Google Scholar

[30] Bellman R., Kashe B.G., Casti J., et al., Differential quadrature: a technique for the rapid solution of non-linear partial differential equations, J. Comput. Phys., 1972, 10, 40-52.10.1016/0021-9991(72)90089-7Search in Google Scholar

[31] Shu C., Differential Quadrature and its Application in Engineering, Athenaeum Press Ltd., Great Britain, 2000.10.1007/978-1-4471-0407-0Search in Google Scholar

[32] Quan J.R., Chang C.T., et al., Newinsights in solving distributed system equations by the quadrature methods, I. Comput. Chem. Eng., 1989, 13, 779 - 788.10.1016/0098-1354(89)85051-3Search in Google Scholar

[33] Quan J.R., Chang C.T., et al., Newinsights in solving distributed system equations by the quadrature methods-II, Comput. Chem. Eng., 1989, 13, 1017-1024.10.1016/0098-1354(89)87043-7Search in Google Scholar

[34] Korkmaz A., Daĝ I., et al., Shock wave simulations using Sinc DQM, Eng. Comput., 2011, 28(6), 654 - 674.10.1108/02644401111154619Search in Google Scholar

[35] Korkmaz A., Daĝ, I., et al., A differential quadrature algorithm for nonlinear Schrödinger equation. Nonlinear Dynamics, 2009, 56(1-2), 69-83.10.1007/s11071-008-9380-0Search in Google Scholar

[36] Mittal R.C., Jiwari R., Sharma K.K., et al., A numerical scheme based on differential quadrature method to solve time dependent Burgers’ equation. Engineering Computations, 2012, 30(1), 117-131.10.1108/02644401311286071Search in Google Scholar

[37] Tomasiello S., Solving 2D-wave problems by the iterative differential quadrature method, Int. J. Comput. Math., 2011, 88(12), 2550-2566.10.1080/00207160.2010.543133Search in Google Scholar

[38] Tomasiello S., Stability and accuracy of the iterative differential quadrature method. International journal for numerical methods in engineering, 2003, 58(9), 1277-1296.10.1002/nme.815Search in Google Scholar

[39] Tomasiello S., Numerical solutions of the Burgers-Huxley equation by the IDQ method, Int. J. Comput. Math, 2010, 87(1), 129-40.10.1080/00207160801968762Search in Google Scholar

[40] Tomasiello S., Numerical stability of DQ solutions of wave problems, Numerical Algorithms, 2011, 57(3), 289-312.10.1007/s11075-010-9429-2Search in Google Scholar

[41] Korkmaz A., Numerical algorithms for solutions of Korteweg-de Vries equation. Numer. Methods Partial Differential Equations, 2010, 26(6), 1504-1521.10.1002/num.20505Search in Google Scholar

[42] Arora G., Singh B.K., et al., Numerical solution of Burgers’ equation with modified cubic B-spline differential quadrature method, Applied Mathematics and Computation, 2013, 224, 166-177.10.1016/j.amc.2013.08.071Search in Google Scholar

[43] Jiwari R., Yuan J., et al., A computational modeling of two dimensional reaction–diffusion Brusselator system arising in chemical processes, 2014, 52(6), 1535-1551.10.1007/s10910-014-0333-1Search in Google Scholar

[44] Bashan A., Karakoc S.B.G., Geyikli T., et al., Approximation of the KdVB equation by the quintic B-spline differential quadrature method, Kuwait Journal of Science, 2015, 42(2).Search in Google Scholar

[45] Korkmaz A., Daĝ I., et al., Quartic and quintic B-spline methods for advection-diffusion equation, Appl. Math. Comput., 2016, 274, 201-219.10.1016/j.amc.2015.11.004Search in Google Scholar

[46] Krowick A., Hermite type radial basis finction-based differential quadrature method for higher order equations, Applied Mathematical modelling, 2016, 40(3), 2421-2430.10.1016/j.apm.2015.09.069Search in Google Scholar

[47] Barrera D., Gonzàlez P., Ibàñez F., Ibàñez M.J., et al., A general soline differential quadrature method based on quasi-interpolation, J Comput. Appl. Math., 2015, 275, 465-479.10.1016/j.cam.2014.02.006Search in Google Scholar

[48] Shukla H.S., Tamsir M., Srivastava V.K., Kumar J., et al., Numerical solution of two dimensional coupled viscous Burger equation using modified cubic B-spline differential quadrature method, AIP Advances, 2014, 4, 117-134.10.1063/1.4902507Search in Google Scholar

[49] Shukla H.S., Tamsir M., Srivastava V.K., et al., Numerical simulation of two dimensional sine-Gordon solitons using modified cubic B-spline differential quadrature method, AIP Advances, 2015, 5, 017121.10.1063/1.4906256Search in Google Scholar

[50] Shukla H.S., Tamsir M., Srivastava V.K., Rashidi M. M. et al., Modified cubic B-spline differential quadrature method for numerical solution of three-dimensional coupled viscous Burger equation, Modern Physics Letters B 30(11), Modern Physics Letters B, 2016, 30(11), 16501 10-27.10.1142/S0217984916501104Search in Google Scholar

[51] Mittal R.C., Jain R.K., et al., Numerical solutions of nonlinear Burgers’ equation with modified cubic B-splines collocation method, Appl. Math. Comput., 2012, 218, 7839-7855.10.1016/j.amc.2012.01.059Search in Google Scholar

[52] Korkmaz A., Daĝ I., Cubic B-spline differential quadrature method for advection-diffusion equations. International Journal of Numerical Methods for Heat and Fluid flow, 2012, 22(8), 1021-1036.10.1108/09615531211271844Search in Google Scholar

[53] Korkmaz A., Akmaz H.K., et al., Numerical simulations for transport of conservative pollutants, Selcuk Journal of Applied Mathematics, 2015, 16(1).Search in Google Scholar

[54] Jain M.K., Numerical Solution of Differential Equations, 2nd ed., Wiley, New York, NY, 1983.Search in Google Scholar

[55] Ruuth S.J., and Spiteri R.J., et al., A new class of optimal high-order strongstability-preserving time-stepping schemes, SIAM J. Numer. Anal., 2002, 40 (2), 469-491.10.1137/S0036142901389025Search in Google Scholar

[56] Shu C., Total-variation-diminishing time discretizations, SIAM J. Sci. Statist. Comput., 1988, 9, 1073-1084.10.1137/0909073Search in Google Scholar

[57] Ruuth S.J., Spiteri R.J., et al., Two barriers on strong-stability-preserving time discretization methods, J. Sci. Comp., 2002, 17, 211-220.10.1023/A:1015156832269Search in Google Scholar

[58] Stenger F., Numerical methods based on Wittaker cardinal or sinc functions, SIAM Review, 1983, 23, 165-224.10.1137/1023037Search in Google Scholar

[59] Bhatt H.P., Khaliq A.Q.M., et al., The locally extrapolated exponential time differencing LOD scheme for multidimensional reaction–diffusion systems, J. Comput. Appl. Math., 2015, 285, 256 - 278.10.1016/j.cam.2015.02.017Search in Google Scholar

Received: 2016-9-20
Accepted: 2017-2-17
Published Online: 2017-7-3

© 2017 S. Dahiya and R. C. Mittal

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

Articles in the same Issue

  1. Regular Articles
  2. Analysis of a New Fractional Model for Damped Bergers’ Equation
  3. Regular Articles
  4. Optimal homotopy perturbation method for nonlinear differential equations governing MHD Jeffery-Hamel flow with heat transfer problem
  5. Regular Articles
  6. Semi- analytic numerical method for solution of time-space fractional heat and wave type equations with variable coefficients
  7. Regular Articles
  8. Investigation of a curve using Frenet frame in the lightlike cone
  9. Regular Articles
  10. Construction of complex networks from time series based on the cross correlation interval
  11. Regular Articles
  12. Nonlinear Schrödinger approach to European option pricing
  13. Regular Articles
  14. A modified cubic B-spline differential quadrature method for three-dimensional non-linear diffusion equations
  15. Regular Articles
  16. A new miniaturized negative-index meta-atom for tri-band applications
  17. Regular Articles
  18. Seismic stability of the survey areas of potential sites for the deep geological repository of the spent nuclear fuel
  19. Regular Articles
  20. Distributed containment control of heterogeneous fractional-order multi-agent systems with communication delays
  21. Regular Articles
  22. Sensitivity analysis and economic optimization studies of inverted five-spot gas cycling in gas condensate reservoir
  23. Regular Articles
  24. Quantum mechanics with geometric constraints of Friedmann type
  25. Regular Articles
  26. Modeling and Simulation for an 8 kW Three-Phase Grid-Connected Photo-Voltaic Power System
  27. Regular Articles
  28. Application of the optimal homotopy asymptotic method to nonlinear Bingham fluid dampers
  29. Regular Articles
  30. Analysis of Drude model using fractional derivatives without singular kernels
  31. Regular Articles
  32. An unsteady MHD Maxwell nanofluid flow with convective boundary conditions using spectral local linearization method
  33. Regular Articles
  34. New analytical solutions for conformable fractional PDEs arising in mathematical physics by exp-function method
  35. Regular Articles
  36. Quantum mechanical calculation of electron spin
  37. Regular Articles
  38. CO2 capture by polymeric membranes composed of hyper-branched polymers with dense poly(oxyethylene) comb and poly(amidoamine)
  39. Regular Articles
  40. Chain on a cone
  41. Regular Articles
  42. Multi-task feature learning by using trace norm regularization
  43. Regular Articles
  44. Superluminal tunneling of a relativistic half-integer spin particle through a potential barrier
  45. Regular Articles
  46. Neutrosophic triplet normed space
  47. Regular Articles
  48. Lie algebraic discussion for affinity based information diffusion in social networks
  49. Regular Articles
  50. Radiation dose and cancer risk estimates in helical CT for pulmonary tuberculosis infections
  51. Regular Articles
  52. A comparison study of steady-state vibrations with single fractional-order and distributed-order derivatives
  53. Regular Articles
  54. Some new remarks on MHD Jeffery-Hamel fluid flow problem
  55. Regular Articles
  56. Numerical investigation of magnetohydrodynamic slip flow of power-law nanofluid with temperature dependent viscosity and thermal conductivity over a permeable surface
  57. Regular Articles
  58. Charge conservation in a gravitational field in the scalar ether theory
  59. Regular Articles
  60. Measurement problem and local hidden variables with entangled photons
  61. Regular Articles
  62. Compression of hyper-spectral images using an accelerated nonnegative tensor decomposition
  63. Regular Articles
  64. Fabrication and application of coaxial polyvinyl alcohol/chitosan nanofiber membranes
  65. Regular Articles
  66. Calculating degree-based topological indices of dominating David derived networks
  67. Regular Articles
  68. The structure and conductivity of polyelectrolyte based on MEH-PPV and potassium iodide (KI) for dye-sensitized solar cells
  69. Regular Articles
  70. Chiral symmetry restoration and the critical end point in QCD
  71. Regular Articles
  72. Numerical solution for fractional Bratu’s initial value problem
  73. Regular Articles
  74. Structure and optical properties of TiO2 thin films deposited by ALD method
  75. Regular Articles
  76. Quadruple multi-wavelength conversion for access network scalability based on cross-phase modulation in an SOA-MZI
  77. Regular Articles
  78. Application of ANNs approach for wave-like and heat-like equations
  79. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  80. Study on node importance evaluation of the high-speed passenger traffic complex network based on the Structural Hole Theory
  81. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  82. A mathematical/physics model to measure the role of information and communication technology in some economies: the Chinese case
  83. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  84. Numerical modeling of the thermoelectric cooler with a complementary equation for heat circulation in air gaps
  85. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  86. On the libration collinear points in the restricted three – body problem
  87. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  88. Research on Critical Nodes Algorithm in Social Complex Networks
  89. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  90. A simulation based research on chance constrained programming in robust facility location problem
  91. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  92. A mathematical/physics carbon emission reduction strategy for building supply chain network based on carbon tax policy
  93. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  94. Mathematical analysis of the impact mechanism of information platform on agro-product supply chain and agro-product competitiveness
  95. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  96. A real negative selection algorithm with evolutionary preference for anomaly detection
  97. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  98. A privacy-preserving parallel and homomorphic encryption scheme
  99. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  100. Random walk-based similarity measure method for patterns in complex object
  101. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  102. A Mathematical Study of Accessibility and Cohesion Degree in a High-Speed Rail Station Connected to an Urban Bus Transport Network
  103. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  104. Design and Simulation of the Integrated Navigation System based on Extended Kalman Filter
  105. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  106. Oil exploration oriented multi-sensor image fusion algorithm
  107. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  108. Analysis of Product Distribution Strategy in Digital Publishing Industry Based on Game-Theory
  109. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  110. Expanded Study on the accumulation effect of tourism under the constraint of structure
  111. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  112. Unstructured P2P Network Load Balance Strategy Based on Multilevel Partitioning of Hypergraph
  113. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  114. Research on the method of information system risk state estimation based on clustering particle filter
  115. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  116. Demand forecasting and information platform in tourism
  117. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  118. Physical-chemical properties studying of molecular structures via topological index calculating
  119. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  120. Local kernel nonparametric discriminant analysis for adaptive extraction of complex structures
  121. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  122. City traffic flow breakdown prediction based on fuzzy rough set
  123. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  124. Conservation laws for a strongly damped wave equation
  125. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  126. Blending type approximation by Stancu-Kantorovich operators based on Pólya-Eggenberger distribution
  127. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  128. Computing the Ediz eccentric connectivity index of discrete dynamic structures
  129. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  130. A discrete epidemic model for bovine Babesiosis disease and tick populations
  131. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  132. Study on maintaining formations during satellite formation flying based on SDRE and LQR
  133. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  134. Relationship between solitary pulmonary nodule lung cancer and CT image features based on gradual clustering
  135. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  136. A novel fast target tracking method for UAV aerial image
  137. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  138. Fuzzy comprehensive evaluation model of interuniversity collaborative learning based on network
  139. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  140. Conservation laws, classical symmetries and exact solutions of the generalized KdV-Burgers-Kuramoto equation
  141. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  142. After notes on self-similarity exponent for fractal structures
  143. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  144. Excitation probability and effective temperature in the stationary regime of conductivity for Coulomb Glasses
  145. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  146. Comparisons of feature extraction algorithm based on unmanned aerial vehicle image
  147. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  148. Research on identification method of heavy vehicle rollover based on hidden Markov model
  149. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  150. Classifying BCI signals from novice users with extreme learning machine
  151. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  152. Topics on data transmission problem in software definition network
  153. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  154. Statistical inferences with jointly type-II censored samples from two Pareto distributions
  155. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  156. Estimation for coefficient of variation of an extension of the exponential distribution under type-II censoring scheme
  157. Special issue on Nonlinear Dynamics in General and Dynamical Systems in particular
  158. Analysis on trust influencing factors and trust model from multiple perspectives of online Auction
  159. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  160. Coupling of two-phase flow in fractured-vuggy reservoir with filling medium
  161. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  162. Production decline type curves analysis of a finite conductivity fractured well in coalbed methane reservoirs
  163. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  164. Flow Characteristic and Heat Transfer for Non-Newtonian Nanofluid in Rectangular Microchannels with Teardrop Dimples/Protrusions
  165. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  166. The size prediction of potential inclusions embedded in the sub-surface of fused silica by damage morphology
  167. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  168. Research on carbonate reservoir interwell connectivity based on a modified diffusivity filter model
  169. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  170. The method of the spatial locating of macroscopic throats based-on the inversion of dynamic interwell connectivity
  171. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  172. Unsteady mixed convection flow through a permeable stretching flat surface with partial slip effects through MHD nanofluid using spectral relaxation method
  173. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  174. A volumetric ablation model of EPDM considering complex physicochemical process in porous structure of char layer
  175. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  176. Numerical simulation on ferrofluid flow in fractured porous media based on discrete-fracture model
  177. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  178. Macroscopic lattice Boltzmann model for heat and moisture transfer process with phase transformation in unsaturated porous media during freezing process
  179. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  180. Modelling of intermittent microwave convective drying: parameter sensitivity
  181. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  182. Simulating gas-water relative permeabilities for nanoscale porous media with interfacial effects
  183. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  184. Simulation of counter-current imbibition in water-wet fractured reservoirs based on discrete-fracture model
  185. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  186. Investigation effect of wettability and heterogeneity in water flooding and on microscopic residual oil distribution in tight sandstone cores with NMR technique
  187. Special Issue on Advances on Modelling of Flowing and Transport in Porous Media
  188. Analytical modeling of coupled flow and geomechanics for vertical fractured well in tight gas reservoirs
  189. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  190. Special Issue: Ever New "Loopholes" in Bell’s Argument and Experimental Tests
  191. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  192. The ultimate loophole in Bell’s theorem: The inequality is identically satisfied by data sets composed of ±1′s assuming merely that they exist
  193. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  194. Erratum to: The ultimate loophole in Bell’s theorem: The inequality is identically satisfied by data sets composed of ±1′s assuming merely that they exist
  195. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  196. Rhetoric, logic, and experiment in the quantum nonlocality debate
  197. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  198. What If Quantum Theory Violates All Mathematics?
  199. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  200. Relativity, anomalies and objectivity loophole in recent tests of local realism
  201. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  202. The photon identification loophole in EPRB experiments: computer models with single-wing selection
  203. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  204. Bohr against Bell: complementarity versus nonlocality
  205. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  206. Is Einsteinian no-signalling violated in Bell tests?
  207. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  208. Bell’s “Theorem”: loopholes vs. conceptual flaws
  209. Special Issue on Ever-New "Loopholes" in Bell’s Argument and Experimental Tests
  210. Nonrecurrence and Bell-like inequalities
  211. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  212. Three-dimensional computer models of electrospinning systems
  213. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  214. Electric field computation and measurements in the electroporation of inhomogeneous samples
  215. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  216. Modelling of magnetostriction of transformer magnetic core for vibration analysis
  217. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  218. Comparison of the fractional power motor with cores made of various magnetic materials
  219. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  220. Dynamics of the line-start reluctance motor with rotor made of SMC material
  221. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  222. Inhomogeneous dielectrics: conformal mapping and finite-element models
  223. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  224. Topology optimization of induction heating model using sequential linear programming based on move limit with adaptive relaxation
  225. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  226. Detection of inter-turn short-circuit at start-up of induction machine based on torque analysis
  227. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  228. Current superimposition variable flux reluctance motor with 8 salient poles
  229. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  230. Modelling axial vibration in windings of power transformers
  231. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  232. Field analysis & eddy current losses calculation in five-phase tubular actuator
  233. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  234. Hybrid excited claw pole generator with skewed and non-skewed permanent magnets
  235. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  236. Electromagnetic phenomena analysis in brushless DC motor with speed control using PWM method
  237. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  238. Field-circuit analysis and measurements of a single-phase self-excited induction generator
  239. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  240. A comparative analysis between classical and modified approach of description of the electrical machine windings by means of T0 method
  241. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  242. Field-based optimal-design of an electric motor: a new sensitivity formulation
  243. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  244. Application of the parametric proper generalized decomposition to the frequency-dependent calculation of the impedance of an AC line with rectangular conductors
  245. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  246. Virtual reality as a new trend in mechanical and electrical engineering education
  247. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  248. Holonomicity analysis of electromechanical systems
  249. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  250. An accurate reactive power control study in virtual flux droop control
  251. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  252. Localized probability of improvement for kriging based multi-objective optimization
  253. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  254. Research of influence of open-winding faults on properties of brushless permanent magnets motor
  255. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  256. Optimal design of the rotor geometry of line-start permanent magnet synchronous motor using the bat algorithm
  257. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  258. Model of depositing layer on cylindrical surface produced by induction-assisted laser cladding process
  259. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  260. Detection of inter-turn faults in transformer winding using the capacitor discharge method
  261. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  262. A novel hybrid genetic algorithm for optimal design of IPM machines for electric vehicle
  263. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  264. Lamination effects on a 3D model of the magnetic core of power transformers
  265. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  266. Detection of vertical disparity in three-dimensional visualizations
  267. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  268. Calculations of magnetic field in dynamo sheets taking into account their texture
  269. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  270. 3-dimensional computer model of electrospinning multicapillary unit used for electrostatic field analysis
  271. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  272. Optimization of wearable microwave antenna with simplified electromagnetic model of the human body
  273. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  274. Induction heating process of ferromagnetic filled carbon nanotubes based on 3-D model
  275. Special Issue: The 18th International Symposium on Electromagnetic Fields in Mechatronics, Electrical and Electronic Engineering ISEF 2017
  276. Speed control of an induction motor by 6-switched 3-level inverter
Downloaded on 12.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2017-0051/html?lang=en
Scroll to top button