Search results
1 – 10 of 29Yasir Khan, Ishtiaq Ali, S. Islam and Qing Biao Wu
The purpose of this paper is to present some numerical methods based on different time stepping and space discretization methods for the Allen‐Cahn equation with non‐periodic…
Abstract
Purpose
The purpose of this paper is to present some numerical methods based on different time stepping and space discretization methods for the Allen‐Cahn equation with non‐periodic boundary conditions.
Design/methodology/approach
In space the equation is discretized by the Chebyshev spectral method, while in time the exponential time differencing fourth‐order Runge‐Kutta (ETDRK4) and implicit‐explicit scheme are used. Also, for comparison the finite difference scheme in both space and time is used.
Findings
It is found that the use of implicit‐explicit scheme allows use of a large time‐step, since an explicit method has less order of accuracy as compared to implicit‐explicit method. In time‐stepping the proposed ETDRK4 does not behave well for this special kind of partial differential equation.
Originality/value
The paper presents some numerical methods based on different time stepping and space discretization methods for the Allen‐Cahn equation with non‐periodic boundary conditions.
Details
Keywords
Somnath Santra, Shubhadeep Mandal and Suman Chakraborty
The purpose of this study is to perform a detailed review on the numerical modeling of multiphase and multicomponent flows in microfluidic system using phase-field method. The…
Abstract
Purpose
The purpose of this study is to perform a detailed review on the numerical modeling of multiphase and multicomponent flows in microfluidic system using phase-field method. The phase-field method is of emerging importance in numerical computation of transport phenomena involving multiple phases and/or components. This method is not only used to model interfacial phenomena typical to multiphase flows encountered in engineering and nature but also turns out to be a promising tool in modeling the dynamics of complex fluid-fluid interfaces encountered in physiological systems such as dynamics of vesicles and red blood cells). Intrinsically, a priori unknown topological evolution of interfaces offers to be the most concerning challenge toward accurate modeling of moving boundary problems. However, the numerical difficulties can be tackled simultaneously with numerical convenience and thermodynamic rigor in the paradigm of the phase field method.
Design/methodology/approach
The phase-field method replaces the macroscopically sharp interfaces separating the fluids by a diffuse transition layer where the interfacial forces are smoothly distributed. As against the moving mesh methods (Lagrangian) for the explicit tracking of interfaces, the phase-field method implicitly captures the same through the evolution of a phase-field function (Eulerian). In contrast to the deployment of an artificially smoothing function for the interface as used in the volume of a fluid or level set method, however, the phase-field method uses mixing free energy for describing the interface. This needs the consideration of an additional equation for an order parameter. The dynamic evolution of the system (equation for order parameter) can be described by Allen–Cahn or Cahn–Hilliard formulation, which couples with the Navier–Stokes equation with the aid of a forcing function that depends on the chemical potential and the gradient of the order parameter.
Findings
In this review, first, the authors discuss the broad motivation and the fundamental theoretical foundation associated with phase-field modeling from the perspective of computational microfluidics. They subsequently pinpoint the outstanding numerical challenges, including estimations of the model-free parameters. They outline some numerical examples, including electrohydrodynamic flows, to demonstrate the efficacy of the method. Finally, they pinpoint various emerging issues and futuristic perspectives connecting the phase-field method and computational microfluidics.
Originality/value
This paper gives unique perspectives to future directions of research on this topic.
Details
Keywords
Yunxia Sun, Xufeng Xiao, Zhiming Gao and Xinlong Feng
The purpose of this paper is to propose an efficient space-time operator-splitting method for the high-dimensional vector-valued Allen–Cahn (AC) equations. The key of the…
Abstract
Purpose
The purpose of this paper is to propose an efficient space-time operator-splitting method for the high-dimensional vector-valued Allen–Cahn (AC) equations. The key of the space-time operator-splitting is to devide the complex partial differential equations into simple heat equations and nolinear ordinary differential equations.
Design/methodology/approach
Each component of high-dimensional heat equations is split into a series of one-dimensional heat equations in different spatial directions. The nonlinear ordinary differential equations are solved by a stabilized semi-implicit scheme to preserve the upper bound of the solution. The algorithm greatly reduces the computational complexity and storage requirement.
Findings
The theoretical analyses of stability in terms of upper bound preservation and mass conservation are shown. The numerical results of phase separation, evolution of the total free energy and total mass conservation show the effectiveness and accuracy of the space-time operator-splitting method.
Practical implications
Extensive 2D/3D numerical tests demonstrated the efficacy and accuracy of the proposed method.
Originality/value
The space-time operator-splitting method reduces the complexity of the problem and reduces the storage space by turning the high-dimensional problem into a series of 1D problems. We give the theoretical analyses of upper bound preservation and mass conservation for the proposed method.
Details
Keywords
Xufeng Xiao, Dongwei Gui and Xinlong Feng
This study aims to present a highly efficient operator-splitting finite element method for the nonlinear two-dimensional/three-dimensional (2D/3D) Allen–Cahn (AC) equation which…
Abstract
Purpose
This study aims to present a highly efficient operator-splitting finite element method for the nonlinear two-dimensional/three-dimensional (2D/3D) Allen–Cahn (AC) equation which describes the anti-phase domain coarsening in a binary alloy. This method is presented to overcome the higher storage requirements, computational complexity and the nonlinear term in numerical computation for the 2D/3D AC equation.
Design/methodology/approach
In each time interval, the authors first split the original equation into a heat equation and a nonlinear equation. Then, they split the high-dimensional heat equation into a series of one-dimensional (1D) heat equations. By solving each 1D subproblem, the authors obtain a numerical solution for heat equation and take it as an initial for the nonlinear equation, which is solved analytically.
Findings
The authors show that the proposed method is unconditionally stable. Finally, various numerical experiments are presented to confirm the high accuracy and efficiency of this method.
Originality/value
A new operator-splitting method is presented for solving the 2D/3D parabolic equation. The 2D/3D parabolic equation is split into a sequence of 1D parabolic equations. In comparison with standard finite element method, the present method can save much central processing unit time. Stability analysis and error estimates are derived and numerical results are presented to support the theoretical analysis.
Details
Keywords
Jae Seok Choi, Takayuki Yamada, Kazuhiro Izui, Shinji Nishiwaki, Heeseung Lim and Jeonghoon Yoo
The purpose of this paper is to present an optimization method for flux barrier designs in interior permanent magnet (IPM) synchronous motors that aims to produce an advantageous…
Abstract
Purpose
The purpose of this paper is to present an optimization method for flux barrier designs in interior permanent magnet (IPM) synchronous motors that aims to produce an advantageous sinusoidal flux density distribution in the air-gap.
Design/methodology/approach
The optimization is based on the phase field method using an Allen-Cahn equation. This approach is a numerical technique for tracking diffuse interfaces like the level set method based on the Hamilton-Jacobi equation.
Findings
The optimization results of IPM motor designs are highly dependent on the initial flux barrier shapes. The authors solve the optimization problem using two different initial shapes, and the optimized models show considerable reductions in torque pulsation and the higher harmonics of back-electromotive force.
Originality/value
This paper presents the optimization method based on the phase field for the design of rotor flux barriers, and proposes a novel interpolation scheme of the magnetic reluctivity.
Details
Keywords
Przemysław Smakulski, Sławomir Pietrowicz and Jun Ishimoto
This paper aims to describe and investigate the mathematical models and numerical modeling of how a cell membrane is affected by a transient ice freezing front combined with the…
Abstract
Purpose
This paper aims to describe and investigate the mathematical models and numerical modeling of how a cell membrane is affected by a transient ice freezing front combined with the influence of thermal fluctuations and anisotropy.
Design/methodology/approach
The study consists of mathematical modeling, validation with an analytical solution, and shows the influence of thermal noises on phase front dynamics and how it influences the freezing process of a single red blood cell. The numerical calculation has been modeled in the framework of the phase field method with a Cahn–Hilliard formulation of a free energy functional.
Findings
The results show an influence scale on directional phase front propagation dynamics and how significant are stochastic thermal noises in micro-scale freezing.
Originality/value
The numerical calculation has modeled in the framework of the phase field method with a Cahn–Hilliard formulation of a free energy functional.
Details
Keywords
Gang Peng, Zhiming Gao, Wenjing Yan and Xinlong Feng
This paper aims to consider numerical simulation for radionuclide transport calculations in geological radioactive waste repository.
Abstract
Purpose
This paper aims to consider numerical simulation for radionuclide transport calculations in geological radioactive waste repository.
Design/methodology/approach
The nonlinear two-point flux approximation is used to discretize the diffusion flux and has a fixed stencil. The cell-vertex unknowns are applied to define the auxiliary unknowns and can be interpolated by the cell-centered unknowns. The approximation of convection flux is based on the second-order upwind method with a slope limiter.
Findings
Numerical results illustrate that the positivity-preserving is satisfied in solving this convection-diffusion system and has a second-order convergence rate on the distorted meshes.
Originality/value
A new positivity-preserving nonlinear finite volume scheme is proposed to simulate the far-field model used in the geological radioactive waste repository. Numerical results illustrate that the positivity-preserving is satisfied in solving this convection-diffusion system and has a second-order convergence rate on the distorted meshes.
Details
Keywords
A. A. Alanazi, Sultan Z. Alamri, S. Shafie and Shazirawati Mohd Puzi
The purpose of this paper is to obtain the nonlinear Schrodinger equation (NLSE) numerical solutions in the presence of the first-order chromatic dispersion using a second-order…
Abstract
Purpose
The purpose of this paper is to obtain the nonlinear Schrodinger equation (NLSE) numerical solutions in the presence of the first-order chromatic dispersion using a second-order, unconditionally stable, implicit finite difference method. In addition, stability and accuracy are proved for the resulting scheme.
Design/methodology/approach
The conserved quantities such as mass, momentum and energy are calculated for the system governed by the NLSE. Moreover, the robustness of the scheme is confirmed by conducting various numerical tests using the Crank-Nicolson method on different cases of solitons to discuss the effects of the factor considered on solitons properties and on conserved quantities.
Findings
The Crank-Nicolson scheme has been derived to solve the NLSE for optical fibers in the presence of the wave packet drift effects. It has been founded that the numerical scheme is second-order in time and space and unconditionally stable by using von-Neumann stability analysis. The effect of the parameters considered in the study is displayed in the case of one, two and three solitons. It was noted that the reliance of NLSE numeric solutions properties on coefficients of wave packets drift, dispersions and Kerr nonlinearity play an important control not only the stable and unstable regime but also the energy, momentum conservation laws. Accordingly, by comparing our numerical results in this study with the previous work, it was recognized that the obtained results are the generalized formularization of these work. Also, it was distinguished that our new data are regarding to the new communications modes that depend on the dispersion, wave packets drift and nonlinearity coefficients.
Originality/value
The present study uses the first-order chromatic. Also, it highlights the relationship between the parameters of dispersion, nonlinearity and optical wave properties. The study further reports the effect of wave packet drift, dispersions and Kerr nonlinearity play an important control not only the stable and unstable regime but also the energy, momentum conservation laws.
Details
Keywords
A.A. Aderogba, M. Chapwanya and J.K. Djoko
For a partial differential equation with a fourth-order derivative such as the Cahn-Hilliard equation, it is always a challenge to design numerical schemes that can handle the…
Abstract
Purpose
For a partial differential equation with a fourth-order derivative such as the Cahn-Hilliard equation, it is always a challenge to design numerical schemes that can handle the restrictive time step introduced by this higher order term. The purpose of this paper is to employ a fractional splitting method to isolate the convective, the nonlinear second-order and the fourth-order differential terms.
Design/methodology/approach
The full equation is then solved by consistent schemes for each differential term independently. In addition to validating the second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.
Findings
The scheme is second-order accuracy, the authors will demonstrate the efficiency of the proposed method by validating the dissipation of the Ginzberg-Lindau energy and the coarsening properties of the solution.
Originality/value
The authors believe that this is the first time the equation is handled numerically using the fractional step method. Apart from the fact that the fractional step method substantially reduces computational time, it has the advantage of simplifying a complex process efficiently. This method permits the treatment of each segment of the original equation separately and piece them together, in a way that will be explained shortly, without destroying the properties of the equation.
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