Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. ...Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.展开更多
A numerical study of a non-Darcy mixed convective heat and mass transfer flow over a vertical surface embedded in a dispersion, melting, and thermal radiation is porous medium under the effects of double investigated....A numerical study of a non-Darcy mixed convective heat and mass transfer flow over a vertical surface embedded in a dispersion, melting, and thermal radiation is porous medium under the effects of double investigated. The set of governing boundary layer equations and the boundary conditions is transformed into a set of coupled nonlinear ordinary differential equations with the relevant boundary conditions. The transformed equations are solved numerically by using the Chebyshev pseudospectral method. Comparisons of the present results with the existing results in the literature are made, and good agreement is found. Numerical results for the velocity, temperature, concentration profiles, and local Nusselt and Sherwood numbers are discussed for various values of physical parameters.展开更多
A mathematical model is elaborated for the laminar flow of an Eyring-Powell fluid over a stretching sheet.The considered non-Newtonian fluid has Prandtl number larger than one.The effects of variable fluid properties ...A mathematical model is elaborated for the laminar flow of an Eyring-Powell fluid over a stretching sheet.The considered non-Newtonian fluid has Prandtl number larger than one.The effects of variable fluid properties and heat generation/absorption are also discussed.The balance equations for fluid flow are reduced to a set of ordinary differential equations through a similarity transformation and solved numerically using a Chebyshev spectral scheme.The effect of various parameters on the rate of heat transfer in the thermal boundary regime is investigated,i.e.,thermal conductivity,the heat generation/absorption ratio and the mixed convection parameter.Good agreement appears to exist between theoretical predictions and the existing published results.展开更多
In this paper, we consider solving the Helmholtz equation in the Cartesian domain , subject to homogeneous Dirichlet boundary condition, discretized with the Chebyshev pseudo-spectral method. The main purpose of this ...In this paper, we consider solving the Helmholtz equation in the Cartesian domain , subject to homogeneous Dirichlet boundary condition, discretized with the Chebyshev pseudo-spectral method. The main purpose of this paper is to present the formulation of a two-level decomposition scheme for decoupling the linear system obtained from the discretization into independent subsystems. This scheme takes advantage of the homogeneity property of the physical problem along one direction to reduce a 2D problem to several 1D problems via a block diagonalization approach and the reflexivity property along the second direction to decompose each of the 1D problems to two independent subproblems using a reflexive decomposition, effectively doubling the number of subproblems. Based on the special structure of the coefficient matrix of the linear system derived from the discretization and a reflexivity property of the second-order Chebyshev differentiation matrix, we show that the decomposed submatrices exhibits a similar property, enabling the system to be decomposed using reflexive decompositions. Explicit forms of the decomposed submatrices are derived. The decomposition not only yields more efficient algorithm but introduces coarse-grain parallelism. Furthermore, it preserves all eigenvalues of the original matrix.展开更多
The evolution of the three-dimensional time-developing mixing layer was simulated numerically using the pseudo-spectral method. The initial perturbations consisted of the two-dimensional fundamental wave and the stre...The evolution of the three-dimensional time-developing mixing layer was simulated numerically using the pseudo-spectral method. The initial perturbations consisted of the two-dimensional fundamental wave and the streamwise-invariant three-dimensional disturbance. A comparison of the formations of the streamwise vortices with different amplitude functions for three-dimensional disturbances was made. In one case the results are similar to that of Rogers and Moser (1992), whereas a different way in which the quadrupole forms and sudden expansion of the rib were observed in another case. The simulation also confirms that stretching by the forming roller rather than Rayleigh centrifugal instability is responsible for the formation of the rib. Finally, numerical flow visualization results were presented. (Edited author abstract) 9 Refs.展开更多
The evolution of the coherent structures in a two-dimensional time-developing mixing layer of the FENE-P fluids is examined numerically. By the means of an appropriate filtering for the polymer stress, some characteri...The evolution of the coherent structures in a two-dimensional time-developing mixing layer of the FENE-P fluids is examined numerically. By the means of an appropriate filtering for the polymer stress, some characteristics of the coherent structures at high b were obtained, which Azaiez and Homsy did not address. The results indicate that adding polymer to the Newtonian fluids will cause stronger vorticity diffusion, accompanied with weaker fundamental and subharmonical perturbations and slower rotational motion of neighbouring vortices during pairing. This effect decreases with the Weissenberg number, but increases with b. In addition, the time when the consecutive rollers are completely coalesced into one delays in the viscoelastic mixing layer compared with the Newtonian one of the same total viscosity.展开更多
Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high ...Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high Reynolds numbers and low Deborah numbers. The effect of viscoelastics on the evolution of the large coherent structure was shown by making a comparison between the second-order and Newtonian fluids at the same Reynolds numbers.展开更多
In this paper we describe a new pseudo-spectral method to solve numerically two and three-dimensional nonlinear diffusion equations over unbounded domains,taking Hermite functions,sinc functions,and rational Chebyshev...In this paper we describe a new pseudo-spectral method to solve numerically two and three-dimensional nonlinear diffusion equations over unbounded domains,taking Hermite functions,sinc functions,and rational Chebyshev polynomials as basis functions.The idea is to discretize the equations by means of differentiation matrices and to relate them to Sylvester-type equations by means of a fourth-order implicit-explicit scheme,being of particular interest the treatment of three-dimensional Sylvester equations that we make.The resulting method is easy to understand and express,and can be implemented in a transparent way by means of a few lines of code.We test numerically the three choices of basis functions,showing the convenience of this new approach,especially when rational Chebyshev polynomials are considered.展开更多
文摘Fractional differential equations have recently been applied in various areas of engineering, science, finance, applied mathematics, bio-engineering and others. However, many researchers remain unaware of this field. In this paper, an efficient numerical method for solving the fractional Advection-dispersion equation (ADE) is considered. The fractional derivative is described in the Caputo sense. The method is based on Chebyshev approximations. The properties of Chebyshev polynomials are used to reduce ADE to a system of ordinary differential equations, which are solved using the finite difference method (FDM). Moreover, the convergence analysis and an upper bound of the error for the derived formula are given. Numerical solutions of ADE are presented and the results are compared with the exact solution.
文摘A numerical study of a non-Darcy mixed convective heat and mass transfer flow over a vertical surface embedded in a dispersion, melting, and thermal radiation is porous medium under the effects of double investigated. The set of governing boundary layer equations and the boundary conditions is transformed into a set of coupled nonlinear ordinary differential equations with the relevant boundary conditions. The transformed equations are solved numerically by using the Chebyshev pseudospectral method. Comparisons of the present results with the existing results in the literature are made, and good agreement is found. Numerical results for the velocity, temperature, concentration profiles, and local Nusselt and Sherwood numbers are discussed for various values of physical parameters.
文摘A mathematical model is elaborated for the laminar flow of an Eyring-Powell fluid over a stretching sheet.The considered non-Newtonian fluid has Prandtl number larger than one.The effects of variable fluid properties and heat generation/absorption are also discussed.The balance equations for fluid flow are reduced to a set of ordinary differential equations through a similarity transformation and solved numerically using a Chebyshev spectral scheme.The effect of various parameters on the rate of heat transfer in the thermal boundary regime is investigated,i.e.,thermal conductivity,the heat generation/absorption ratio and the mixed convection parameter.Good agreement appears to exist between theoretical predictions and the existing published results.
文摘In this paper, we consider solving the Helmholtz equation in the Cartesian domain , subject to homogeneous Dirichlet boundary condition, discretized with the Chebyshev pseudo-spectral method. The main purpose of this paper is to present the formulation of a two-level decomposition scheme for decoupling the linear system obtained from the discretization into independent subsystems. This scheme takes advantage of the homogeneity property of the physical problem along one direction to reduce a 2D problem to several 1D problems via a block diagonalization approach and the reflexivity property along the second direction to decompose each of the 1D problems to two independent subproblems using a reflexive decomposition, effectively doubling the number of subproblems. Based on the special structure of the coefficient matrix of the linear system derived from the discretization and a reflexivity property of the second-order Chebyshev differentiation matrix, we show that the decomposed submatrices exhibits a similar property, enabling the system to be decomposed using reflexive decompositions. Explicit forms of the decomposed submatrices are derived. The decomposition not only yields more efficient algorithm but introduces coarse-grain parallelism. Furthermore, it preserves all eigenvalues of the original matrix.
基金The project supported by the Zhejiang Province Natural Science Special Fund for Youth Scientists' Cultivation.
文摘The evolution of the three-dimensional time-developing mixing layer was simulated numerically using the pseudo-spectral method. The initial perturbations consisted of the two-dimensional fundamental wave and the streamwise-invariant three-dimensional disturbance. A comparison of the formations of the streamwise vortices with different amplitude functions for three-dimensional disturbances was made. In one case the results are similar to that of Rogers and Moser (1992), whereas a different way in which the quadrupole forms and sudden expansion of the rib were observed in another case. The simulation also confirms that stretching by the forming roller rather than Rayleigh centrifugal instability is responsible for the formation of the rib. Finally, numerical flow visualization results were presented. (Edited author abstract) 9 Refs.
文摘The evolution of the coherent structures in a two-dimensional time-developing mixing layer of the FENE-P fluids is examined numerically. By the means of an appropriate filtering for the polymer stress, some characteristics of the coherent structures at high b were obtained, which Azaiez and Homsy did not address. The results indicate that adding polymer to the Newtonian fluids will cause stronger vorticity diffusion, accompanied with weaker fundamental and subharmonical perturbations and slower rotational motion of neighbouring vortices during pairing. This effect decreases with the Weissenberg number, but increases with b. In addition, the time when the consecutive rollers are completely coalesced into one delays in the viscoelastic mixing layer compared with the Newtonian one of the same total viscosity.
文摘Numerical simulations have been performed in time-developing plane mixing layers of the viscoelastic second-order fluids with pseudo-spectral method. Roll-up, pairing and merging of large eddies were examined at high Reynolds numbers and low Deborah numbers. The effect of viscoelastics on the evolution of the large coherent structure was shown by making a comparison between the second-order and Newtonian fluids at the same Reynolds numbers.
文摘In this paper we describe a new pseudo-spectral method to solve numerically two and three-dimensional nonlinear diffusion equations over unbounded domains,taking Hermite functions,sinc functions,and rational Chebyshev polynomials as basis functions.The idea is to discretize the equations by means of differentiation matrices and to relate them to Sylvester-type equations by means of a fourth-order implicit-explicit scheme,being of particular interest the treatment of three-dimensional Sylvester equations that we make.The resulting method is easy to understand and express,and can be implemented in a transparent way by means of a few lines of code.We test numerically the three choices of basis functions,showing the convenience of this new approach,especially when rational Chebyshev polynomials are considered.
