The finite element solution of a generalized Stokes system in terms of the flow variables stream function and vorticity is studied. This system results from a time discretization of the time-dependent Stokes system in...The finite element solution of a generalized Stokes system in terms of the flow variables stream function and vorticity is studied. This system results from a time discretization of the time-dependent Stokes system in stream function-vorticity formulation, or yet by the application of the characteristics method to solve the Navier-Stokes equations in the same representation. Numerical results presented for both cases illustrate the good behaviour of the adopted approach.展开更多
In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinea...In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems （with same stiffness matrix but different right-hand side） on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.展开更多
This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through...This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.展开更多
A nonlinear Galerkin finite element method is presented for the two dimensional incom- pressible Navier-Stokes equations with stream-vorticity form.the scheme is based on two finite ele- ment spaces XH and XH for the ...A nonlinear Galerkin finite element method is presented for the two dimensional incom- pressible Navier-Stokes equations with stream-vorticity form.the scheme is based on two finite ele- ment spaces XH and XH for the approximation of the stream and vorticity function ,defined respec- tively on a coarse grid with grid size H and a fine grid with grid size h<<H.We prove that the difference between the new nonlinear Galerkin method and the standard Galerkin method is of the order H2both in stream function and vorticity.展开更多
Taking the distributing calculation of velocity and concentration as an example, the paper established a series of governing equations by the vorticity stream function method, and dispersed the equations by the finit...Taking the distributing calculation of velocity and concentration as an example, the paper established a series of governing equations by the vorticity stream function method, and dispersed the equations by the finite differencing method. After figuring out the distribution field of velocity, the paper also calculated the concentration distribution in sedimentation tank by using the two dimensional concentration transport equation. The validity and feasibility of the numerical method was verified through comparing with experimental data. Furthermore, the paper carried out a tentative exploration into the application of numerical simulation of sedimentation tanks.展开更多
Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid...Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid.Moreover,the scaling between these two grid sizes is super-linear.Approximation,stability and convergence aspects of a fully discrete scheme are analyzed.At last a numrical example is given whose results show that the algorithm proposed in this paper is effcient.展开更多
The unsteady axisymmetric incompressible flow between two concentric spheres was discussed in this paper. It is useful to most astrophysical, geophysical and engineering applications. In order to get the existence and...The unsteady axisymmetric incompressible flow between two concentric spheres was discussed in this paper. It is useful to most astrophysical, geophysical and engineering applications. In order to get the existence and uniqueness of weak solution of this flow with the stream_velocity form, firstly, the relations among the nonlinear terms in this equation is found; then, the existence is proved by an auxiliary semi_discrete scheme and a compactness argument.展开更多
In this paper, problems of the flow over a fat plate in the large Reynolds numbercase are studied by using the method of multiple scales￣[1,2].We have obtained N-orderuniformly valid asymptotic solutions of the Naver...In this paper, problems of the flow over a fat plate in the large Reynolds numbercase are studied by using the method of multiple scales￣[1,2].We have obtained N-orderuniformly valid asymptotic solutions of the Naver-Stodes equations.展开更多
This paper deals with a new higher order compact difference scheme, which is, O(h4) using coupled approach on the 19-point 3D stencil for the solution of three dimensional nonlinear biharmonic equations. At each inter...This paper deals with a new higher order compact difference scheme, which is, O(h4) using coupled approach on the 19-point 3D stencil for the solution of three dimensional nonlinear biharmonic equations. At each internal grid point, the solution u(x,y,z) and its Laplacian Δ4u are obtained. The resulting stencil algo-rithm is presented and hence this new algorithm can be easily incorporated to solve many problems. The present discretization allows us to use the Dirichlet boundary conditions only and there is no need to discretize the derivative boundary conditions near the boundary. We also show that special treatment is required to handle the boundary conditions. Convergence analysis for a model problem is briefly discussed. The method is tested on three problems and compares very favourably with the corresponding second order approximation which we also discuss using coupled approach.展开更多
In our previous works, we suggest that quantum particles are composite physical objects endowed with the geometric and topological structures of their corresponding differentiable manifolds that would allow them to im...In our previous works, we suggest that quantum particles are composite physical objects endowed with the geometric and topological structures of their corresponding differentiable manifolds that would allow them to imitate and adapt to physical environments. In this work, we show that Dirac equation in fact describes quantum particles as composite structures that are in a fluid state in which the components of the wavefunction can be identified with the stream function and the velocity potential of a potential flow formulated in the theory of classical fluids. We also show that Dirac quantum particles can manifest as standing waves which are the result of the superposition of two fluid flows moving in opposite directions. However, for a steady motion a Dirac quantum particle does not exhibit a wave motion even though it has the potential to establish a wave within its physical structure, therefore, without an external disturbance a Dirac quantum particle may be considered as a classical particle defined in classical physics. And furthermore, from the fact that there are two identical fluid flows in opposite directions within their physical structures, the fluid state model of Dirac quantum particles can be used to explain why fermions are spin-half particles.展开更多
Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level met...Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.展开更多
The authors propose a numerical algorithm for the two-dimensional Navier-Stokes equations written in stream function-vorticity formulation. The total time derivative term is treated with a first order characteristics ...The authors propose a numerical algorithm for the two-dimensional Navier-Stokes equations written in stream function-vorticity formulation. The total time derivative term is treated with a first order characteristics method. The space approximation is based on a piecewise continuous finite element method. The proposed algorithm is used to simulate the mechanical aeration process in lakes. Such process is used to combat the degradation of the water quality due to the eutrophication phenomena. For this application high computing facilities and capacities are required. In order to optimize the computing time and make possible the simulation of real applications, the authors propose a parallel implementation of the numerical algorithm. The parallelization technique is performed using the Message Passing Interface. The efficiency of the proposed numerical algorithm is illustrated by some numerical results.展开更多
文摘The finite element solution of a generalized Stokes system in terms of the flow variables stream function and vorticity is studied. This system results from a time discretization of the time-dependent Stokes system in stream function-vorticity formulation, or yet by the application of the characteristics method to solve the Navier-Stokes equations in the same representation. Numerical results presented for both cases illustrate the good behaviour of the adopted approach.
基金supported by National Foundation of Natural Science under the Grant 11071216
文摘In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems （with same stiffness matrix but different right-hand side） on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.
文摘This article gives a general model using specific periodic special functions, that is, degenerate elliptic Weierstrass P functions composed with the LambertW function, whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of arbitrary r balls of the 3-Torus. The continuity equation is satisfied together with spatially periodic boundary conditions. The yicomponent forcing terms consist of a function F as part of its expression that is arbitrarily small in an r ball where it is associated with a singular forcing expression both for inviscid and viscous cases. As a result, a significant simplification occurs with a v3(vifor all velocity components) only governing PDE resulting. The extension of three restricted subspaces in each of the principal directions in the Cartesian plane is shown as the Cartesian product ℋ=Jx,t×Jy,t×Jz,t. On each of these subspaces vi,i=1,2,3is continuous and there exists a linear independent subspace associated with the argument of the W function. Here the 3-Torus is built up from each compact segment of length 2R on each of the axes on the 3 principal directions x, y, and z. The form of the scaled velocities for non zero scaled δis related to the definition of the W function such that e−W(ξ)=W(ξ)ξwhere ξdepends on t and proportional to δ→0for infinite time t. The ratio Wξis equal to 1, making the limit δ→0finite and well defined. Considering r balls where the function F=(x−ai)2+(y−bi)2+(z−ci)2−ηset equal to −1e+rwhere r>0. is such that the forcing is singular at every distance r of centres of cubes each containing an r-ball. At the centre of the balls, the forcing is infinite. The main idea is that a system of singular initial value problems with infinite forcing is to be solved for where the velocities are shown to be locally Hölder continuous. It is proven that the limit of these singular problems shifts the finite time blowup time ti∗for first and higher derivatives to t=∞thereby indicating that there is no finite time blowup. Results in the literature can provide a systematic approach to study both large space and time behaviour for singular solutions to the Navier Stokes equations. Among the references, it has been shown that mathematical tools can be applied to study the asymptotic properties of solutions.
文摘A nonlinear Galerkin finite element method is presented for the two dimensional incom- pressible Navier-Stokes equations with stream-vorticity form.the scheme is based on two finite ele- ment spaces XH and XH for the approximation of the stream and vorticity function ,defined respec- tively on a coarse grid with grid size H and a fine grid with grid size h<<H.We prove that the difference between the new nonlinear Galerkin method and the standard Galerkin method is of the order H2both in stream function and vorticity.
文摘Taking the distributing calculation of velocity and concentration as an example, the paper established a series of governing equations by the vorticity stream function method, and dispersed the equations by the finite differencing method. After figuring out the distribution field of velocity, the paper also calculated the concentration distribution in sedimentation tank by using the two dimensional concentration transport equation. The validity and feasibility of the numerical method was verified through comparing with experimental data. Furthermore, the paper carried out a tentative exploration into the application of numerical simulation of sedimentation tanks.
文摘Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid.Moreover,the scaling between these two grid sizes is super-linear.Approximation,stability and convergence aspects of a fully discrete scheme are analyzed.At last a numrical example is given whose results show that the algorithm proposed in this paper is effcient.
文摘The unsteady axisymmetric incompressible flow between two concentric spheres was discussed in this paper. It is useful to most astrophysical, geophysical and engineering applications. In order to get the existence and uniqueness of weak solution of this flow with the stream_velocity form, firstly, the relations among the nonlinear terms in this equation is found; then, the existence is proved by an auxiliary semi_discrete scheme and a compactness argument.
文摘In this paper, problems of the flow over a fat plate in the large Reynolds numbercase are studied by using the method of multiple scales￣[1,2].We have obtained N-orderuniformly valid asymptotic solutions of the Naver-Stodes equations.
文摘This paper deals with a new higher order compact difference scheme, which is, O(h4) using coupled approach on the 19-point 3D stencil for the solution of three dimensional nonlinear biharmonic equations. At each internal grid point, the solution u(x,y,z) and its Laplacian Δ4u are obtained. The resulting stencil algo-rithm is presented and hence this new algorithm can be easily incorporated to solve many problems. The present discretization allows us to use the Dirichlet boundary conditions only and there is no need to discretize the derivative boundary conditions near the boundary. We also show that special treatment is required to handle the boundary conditions. Convergence analysis for a model problem is briefly discussed. The method is tested on three problems and compares very favourably with the corresponding second order approximation which we also discuss using coupled approach.
文摘In our previous works, we suggest that quantum particles are composite physical objects endowed with the geometric and topological structures of their corresponding differentiable manifolds that would allow them to imitate and adapt to physical environments. In this work, we show that Dirac equation in fact describes quantum particles as composite structures that are in a fluid state in which the components of the wavefunction can be identified with the stream function and the velocity potential of a potential flow formulated in the theory of classical fluids. We also show that Dirac quantum particles can manifest as standing waves which are the result of the superposition of two fluid flows moving in opposite directions. However, for a steady motion a Dirac quantum particle does not exhibit a wave motion even though it has the potential to establish a wave within its physical structure, therefore, without an external disturbance a Dirac quantum particle may be considered as a classical particle defined in classical physics. And furthermore, from the fact that there are two identical fluid flows in opposite directions within their physical structures, the fluid state model of Dirac quantum particles can be used to explain why fermions are spin-half particles.
文摘Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.
文摘The authors propose a numerical algorithm for the two-dimensional Navier-Stokes equations written in stream function-vorticity formulation. The total time derivative term is treated with a first order characteristics method. The space approximation is based on a piecewise continuous finite element method. The proposed algorithm is used to simulate the mechanical aeration process in lakes. Such process is used to combat the degradation of the water quality due to the eutrophication phenomena. For this application high computing facilities and capacities are required. In order to optimize the computing time and make possible the simulation of real applications, the authors propose a parallel implementation of the numerical algorithm. The parallelization technique is performed using the Message Passing Interface. The efficiency of the proposed numerical algorithm is illustrated by some numerical results.
