The two-dimensional nonlinear shallow water equations in the presence of Coriolis force and bottom topography are solved numerically using the fractional steps method. The fractional steps method consists of splitting...The two-dimensional nonlinear shallow water equations in the presence of Coriolis force and bottom topography are solved numerically using the fractional steps method. The fractional steps method consists of splitting the multi-dimensional matrix inversion problem into an equivalent one dimensional problem which is successively integrated in every direction along the characteristics using the Riemann invariant associated with the cubic spline interpolation. The height and the velocity field of the shallow water equations over irregular bottom are discretized on a fixed Eulerian grid and time-stepped using the fractional steps method. Effects of the Coriolis force and the bottom topography for particular initial flows on the velocity components and the free surface elevation have been studied and the results are plotted.展开更多
A split-step second-order predictor-corrector method for space-fractional reaction-diffusion equations with nonhomogeneous boundary conditions is presented and analyzed for the stability and convergence.The matrix tra...A split-step second-order predictor-corrector method for space-fractional reaction-diffusion equations with nonhomogeneous boundary conditions is presented and analyzed for the stability and convergence.The matrix transfer technique is used for spatial discretization of the problem.The method is shown to be unconditionally stable and second-order convergent.Numerical experiments are performed to confirm the stability and secondorder convergence of the method.The split-step predictor-corrector method is also compared with an IMEX predictor-corrector method which is found to incur oscillatory behavior for some time steps.Our method is seen to produce reliable and oscillatioresults for any time step when implemented on numerical examples with nonsmooth initial data.We also present a priori reliability constraint for the IMEX predictor-corrector method to avoid unwanted oscillations and show its validity numerically.展开更多
We propose a fractional-order improved Fitz Hugh–Nagumo(FHN)neuron model in terms of a generalized Caputo fractional derivative.Following the existence of a unique solution for the proposed model,we derive the numeri...We propose a fractional-order improved Fitz Hugh–Nagumo(FHN)neuron model in terms of a generalized Caputo fractional derivative.Following the existence of a unique solution for the proposed model,we derive the numerical solution using a recently proposed L1 predictor–corrector method.The given method is based on the L1-type discretization algorithm and the spline interpolation scheme.We perform the error and stability analyses for the given method.We perform graphical simulations demonstrating that the proposed FHN neuron model generates rich electrical activities of periodic spiking patterns,chaotic patterns,and quasi-periodic patterns.The motivation behind proposing a fractional-order improved FHN neuron model is that such a system can provide a more nuanced description of the process with better understanding and simulation of the neuronal responses by incorporating memory effects and non-local dynamics,which are inherent to many biological systems.展开更多
An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal str...An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.展开更多
In this paper,two crossover hybrid variable-order derivatives of the cancer model are developed.Grünwald-Letnikov approximation is used to approximate the hybrid fractional and variable-order fractional operators...In this paper,two crossover hybrid variable-order derivatives of the cancer model are developed.Grünwald-Letnikov approximation is used to approximate the hybrid fractional and variable-order fractional operators.The existence,uniqueness,and stability of the proposed model are discussed.Adams Bashfourth’s fifth-step method with a hybrid variable-order fractional operator is developed to study the proposed models.Comparative studies with generalized fifth-order Runge-Kutta method are given.Numerical examples and comparative studies to verify the applicability of the used methods and to demonstrate the simplicity of these approximations are presented.We have showcased the efficiency of the proposed method and garnered robust empirical support for our theoretical findings.展开更多
A fractional-step method of predictor-corrector difference-pseudospectrum with unconditional L2-stability and exponential convergence is presented. The stability and convergence of this method is strictly proved mathe...A fractional-step method of predictor-corrector difference-pseudospectrum with unconditional L2-stability and exponential convergence is presented. The stability and convergence of this method is strictly proved mathematically for a nonlinear convection-dominated flow. The error estimation is given and the superiority of this method is verified by numerical test.展开更多
For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-d...For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L 2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.展开更多
An analytic method of fractional steps, which is unconditionally L_∞-stable, is proposed for the numerical solution to convection-dominated problems. In this paper the stability and convergence of the analytic soluti...An analytic method of fractional steps, which is unconditionally L_∞-stable, is proposed for the numerical solution to convection-dominated problems. In this paper the stability and convergence of the analytic solution with fractional steps to both linear and nonlinear problems are proved, and its error estimates are presented.展开更多
For compressible two-phase displacement problem,the modified upwind finite difference fractionalsteps schemes are put forward.Some techniques,such as calculus of variations,commutative law of multiplicationof differen...For compressible two-phase displacement problem,the modified upwind finite difference fractionalsteps schemes are put forward.Some techniques,such as calculus of variations,commutative law of multiplicationof difference operators,decomposition of high order difference operators,the theory of prior estimates and tech-niques are used.Optimal order estimates in L^2 norm are derived for the error in the approximate solution.Thismethod has already been applied to the numerical simulation of seawater intrusion and migration-accumulationof oil resources.展开更多
For nonlinear coupled system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward, trod...For nonlinear coupled system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward, trod two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources.展开更多
For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fract...For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fractional steps techniques are needed to convert a multi-dimensional problem into a series of successive one-dimensional problems. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators, and the theory of prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources and predicting the consequences of seawater intrusion and protection projects.展开更多
This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the erro...This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.展开更多
A stable finite element method for the time dependent Navier-Stokes equations was used for studying the wind flow and pollutant dispersion within street canyons. A three-step fractional method was used to solve the ve...A stable finite element method for the time dependent Navier-Stokes equations was used for studying the wind flow and pollutant dispersion within street canyons. A three-step fractional method was used to solve the velocity field and the pressure field separately from the governing equations. The Streamline Upwind Petrov-Galerkin(SUPG) method was used to get stable numerical results. Numerical oscillation was minimized and satisfactory results can be obtained for flows at high Reynolds numbers. Simulating the flow over a square cylinder within a wide range of Reynolds numbers validates the wind field model. The Strouhal numbers obtained from the numerical simulation had a good agreement with those obtained from experiment. The wind field model developed in the present study is applied to simulate more complex flow phenomena in street canyons with two different building configurations. The results indicated that the flow at rooftop of buildings might not be assumed parallel to the ground as some numerical modelers did. A counter-clockwise rotating vortex may be found in street canyons with an inflow from the left to right. In addition, increasing building height can increase velocity fluctuations in the street canyon under certain circumstances, which facilitate pollutant dispersion. At high Reynolds numbers, the flow regimes in street canyons do not change with inflow velocity.展开更多
文摘The two-dimensional nonlinear shallow water equations in the presence of Coriolis force and bottom topography are solved numerically using the fractional steps method. The fractional steps method consists of splitting the multi-dimensional matrix inversion problem into an equivalent one dimensional problem which is successively integrated in every direction along the characteristics using the Riemann invariant associated with the cubic spline interpolation. The height and the velocity field of the shallow water equations over irregular bottom are discretized on a fixed Eulerian grid and time-stepped using the fractional steps method. Effects of the Coriolis force and the bottom topography for particular initial flows on the velocity components and the free surface elevation have been studied and the results are plotted.
文摘A split-step second-order predictor-corrector method for space-fractional reaction-diffusion equations with nonhomogeneous boundary conditions is presented and analyzed for the stability and convergence.The matrix transfer technique is used for spatial discretization of the problem.The method is shown to be unconditionally stable and second-order convergent.Numerical experiments are performed to confirm the stability and secondorder convergence of the method.The split-step predictor-corrector method is also compared with an IMEX predictor-corrector method which is found to incur oscillatory behavior for some time steps.Our method is seen to produce reliable and oscillatioresults for any time step when implemented on numerical examples with nonsmooth initial data.We also present a priori reliability constraint for the IMEX predictor-corrector method to avoid unwanted oscillations and show its validity numerically.
文摘We propose a fractional-order improved Fitz Hugh–Nagumo(FHN)neuron model in terms of a generalized Caputo fractional derivative.Following the existence of a unique solution for the proposed model,we derive the numerical solution using a recently proposed L1 predictor–corrector method.The given method is based on the L1-type discretization algorithm and the spline interpolation scheme.We perform the error and stability analyses for the given method.We perform graphical simulations demonstrating that the proposed FHN neuron model generates rich electrical activities of periodic spiking patterns,chaotic patterns,and quasi-periodic patterns.The motivation behind proposing a fractional-order improved FHN neuron model is that such a system can provide a more nuanced description of the process with better understanding and simulation of the neuronal responses by incorporating memory effects and non-local dynamics,which are inherent to many biological systems.
基金the National Metal and Materials Technology Centerthe Thailand Research Fund+1 种基金the Office of Higher Education Commissionthe Chulalongkorn University for supporting the present research
文摘An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.
文摘In this paper,two crossover hybrid variable-order derivatives of the cancer model are developed.Grünwald-Letnikov approximation is used to approximate the hybrid fractional and variable-order fractional operators.The existence,uniqueness,and stability of the proposed model are discussed.Adams Bashfourth’s fifth-step method with a hybrid variable-order fractional operator is developed to study the proposed models.Comparative studies with generalized fifth-order Runge-Kutta method are given.Numerical examples and comparative studies to verify the applicability of the used methods and to demonstrate the simplicity of these approximations are presented.We have showcased the efficiency of the proposed method and garnered robust empirical support for our theoretical findings.
基金the National Natural Science Foundation of China
文摘A fractional-step method of predictor-corrector difference-pseudospectrum with unconditional L2-stability and exponential convergence is presented. The stability and convergence of this method is strictly proved mathematically for a nonlinear convection-dominated flow. The error estimation is given and the superiority of this method is verified by numerical test.
基金This work was supported by the Major State Basic Research Program of China(Grant No. 1990328) the National Tackling Key Problem Program, the National Natural Science Foundation of China (Grant Nos. 19871051 and 19972039) the Doctorate Foundation
文摘For combinatorial system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as implicit-explicit difference scheme, calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L 2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources. Keywords: combinatorial system, multilayer dynamics of fluids in porous media, two-class upwind finite difference fractional steps method, convergence, numerical simulation of energy sources.
文摘An analytic method of fractional steps, which is unconditionally L_∞-stable, is proposed for the numerical solution to convection-dominated problems. In this paper the stability and convergence of the analytic solution with fractional steps to both linear and nonlinear problems are proved, and its error estimates are presented.
基金Supported by the Major State Basic Research Program of China (Grant No.1999032803)the National Natural Science Foundation of China (Grant No.10372052,10271066)the Decorate Foundation of the Ministry Education of China (Grant No.20030422047)
文摘For compressible two-phase displacement problem,the modified upwind finite difference fractionalsteps schemes are put forward.Some techniques,such as calculus of variations,commutative law of multiplicationof difference operators,decomposition of high order difference operators,the theory of prior estimates and tech-niques are used.Optimal order estimates in L^2 norm are derived for the error in the approximate solution.Thismethod has already been applied to the numerical simulation of seawater intrusion and migration-accumulationof oil resources.
文摘For nonlinear coupled system of multilayer dynamics of fluids in porous media, the second order and first order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward, trod two-dimensional and three-dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates, are adopted. Optimal order estimates in L2 norm are derived to determine the error in the second order approximate solution. This method has already been applied to the numerical simulation of migration-accumulation of oil resources.
基金Project supported by the Major State Basic Research Program of China (No.G1999032803)the National Tackling Key Problems Program (No.20050200069)the National Natural Science Foundation of China (Nos.10372052, 10271066)the Doctoral Foundation of Ministry of Education of China (No.20030422047).
文摘For the three-dimensional convection-dominated problem of dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Fractional steps techniques are needed to convert a multi-dimensional problem into a series of successive one-dimensional problems. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators, and the theory of prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources and predicting the consequences of seawater intrusion and protection projects.
文摘This paper deals with the inertial manifold and the approximate inertialmanifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore,we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.
文摘A stable finite element method for the time dependent Navier-Stokes equations was used for studying the wind flow and pollutant dispersion within street canyons. A three-step fractional method was used to solve the velocity field and the pressure field separately from the governing equations. The Streamline Upwind Petrov-Galerkin(SUPG) method was used to get stable numerical results. Numerical oscillation was minimized and satisfactory results can be obtained for flows at high Reynolds numbers. Simulating the flow over a square cylinder within a wide range of Reynolds numbers validates the wind field model. The Strouhal numbers obtained from the numerical simulation had a good agreement with those obtained from experiment. The wind field model developed in the present study is applied to simulate more complex flow phenomena in street canyons with two different building configurations. The results indicated that the flow at rooftop of buildings might not be assumed parallel to the ground as some numerical modelers did. A counter-clockwise rotating vortex may be found in street canyons with an inflow from the left to right. In addition, increasing building height can increase velocity fluctuations in the street canyon under certain circumstances, which facilitate pollutant dispersion. At high Reynolds numbers, the flow regimes in street canyons do not change with inflow velocity.
