This article gives a general model using specific periodic special functions, which is degenerate elliptic Weierstrass P functions whose presence in the governing equations through the forcing terms simplify the perio...This article gives a general model using specific periodic special functions, which is degenerate elliptic Weierstrass P functions whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of cells of the 3-Torus. Satisfying a divergence-free vector field and periodic boundary conditions respectively with a general spatio-temporal forcing term which is smooth and spatially periodic, the existence of solutions which have finite time singularities can occur starting with the first derivative and higher with respect to time. The existence of a subspace of the solution space where v<sub>3</sub> is continuous and {C, y<sub>1</sub>, y<sub>1</sub><sup>2</sup>}, is linearly independent in the additive argument of the solution in terms of the Lambert W function, (y<sub>1</sub><sup>2</sup>=y<sub>2</sub>, C∈R) together with the condition v<sub>2</sub>=-2y<sub>1</sub>v<sub>1</sub>. On this subspace, the Biot Savart Law holds exactly [see Section 2 (Equation (13))]. Also on this subspace, an expression X (part of PNS equations) vanishes which contains all the expressions in derivatives of v<sub>1</sub> and v<sub>2</sub> and the forcing terms in the plane which are related as with the cancellation of all such terms in governing PDE. The y<sub>3</sub> component forcing term is arbitrarily small in ε ball where Weierstrass P functions touch the center of the ball both for inviscid and viscous cases. As a result, a significant simplification occurs with a v<sub>3 </sub>only governing PDE resulting. With viscosity present as v changes from zero to the fully viscous case at v =1 the solution for v<sub>3</sub> reaches a peak in the third component y<sub>3</sub>. Consequently, there exists a dipole which is not centered at the center of the cell of the Lattice. Hence since the dipole by definition has an equal in magnitude positive and negative peak in y<sub>3</sub>, then the dipole Riemann cut-off surface is covered by a closed surface which is the sphere and where a given cell of dimensions [-1, 1]<sup>3</sup> is circumscribed on a sphere of radius 1. For such a closed surface containing a dipole it necessarily follows that the flux at the surface of the sphere of v<sub>3</sub> wrt to surface normal n is zero including at the points where the surface of sphere touches the cube walls. At the finite time singularity on the sphere a rotation boundary condition is deduced. It is shown that v<sub>3</sub> is spatially finite on the Riemann Sphere and the forcing is oscillatory in y<sub>3</sub> component if the velocity v3</sub> is. It is true that . A boundary condition on the sphere shows the rotation of a sphere of viscous fluid. Finally on the sphere a solution for v3</sub> is obtained which is proven to be Hölder continuous and it is shown that it is possible to extend Hölder continuity on the sphere uniquely to all of the interior of the ball.展开更多
The forward flight of a model butterfly was stud- ied by simulation using the equations of motion coupled with the Navier-Stokes equations. The model butterfly moved under the action of aerodynamic and gravitational f...The forward flight of a model butterfly was stud- ied by simulation using the equations of motion coupled with the Navier-Stokes equations. The model butterfly moved under the action of aerodynamic and gravitational forces, where the aerodynamic forces were generated by flapping wings which moved with the body, allowing the body os- cillations of the model butterfly to be simulated. The main results are as follows: (1) The aerodynamic force produced by the wings is approximately perpendicular to the long-axis of body and is much larger in the downstroke than in the up- stroke. In the downstroke the body pitch angle is small and the large aerodynamic force points up and slightly backward, giving the weight-supporting vertical force and a small neg- ative horizontal force, whilst in the upstroke, the body an- gle is large and the relatively small aerodynamic force points forward and slightly downward, giving a positive horizon- tal force which overcomes the body drag and the negative horizontal force generated in the downstroke. (2) Pitching oscillation of the butterfly body plays an equivalent role of the wing-rotation of many other insects. (3) The body-mass- specific power of the model butterfly is 33.3 W/kg, not very different from that of many other insects, e.g., fruitflies and dragonflies.展开更多
In the present paper, the longitudinal dynamic flight stability properties of two model insects are predicted by an approximate theory and computed by numerical sim- ulation. The theory is based on the averaged model ...In the present paper, the longitudinal dynamic flight stability properties of two model insects are predicted by an approximate theory and computed by numerical sim- ulation. The theory is based on the averaged model (which assumes that the frequency of wingbeat is sufficiently higher than that of the body motion, so that the flapping wings' degrees of freedom relative to the body can be dropped and the wings can be replaced by wingbeat-cycle-average forces and moments); the simulation solves the complete equations of motion coupled with the Navier-Stokes equations. Comparison between the theory and the simulation provides a test to the validity of the assumptions in the theory. One of the insects is a model dronefly which has relatively high wingbeat frequency (164 Hz) and the other is a model hawkmoth which has relatively low wingbeat frequency (26 Hz). The results show that the averaged model is valid for the hawkmoth as well as for the dronefly. Since the wingbeat frequency of the hawkmoth is relatively low (the characteristic times of the natural modes of motion of the body divided by wingbeat period are relatively large) compared with many other insects, that the theory based on the averaged model is valid for the hawkmoth means that it could be valid for many insects.展开更多
A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and...A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.展开更多
文摘This article gives a general model using specific periodic special functions, which is degenerate elliptic Weierstrass P functions whose presence in the governing equations through the forcing terms simplify the periodic Navier Stokes equations (PNS) at the centers of cells of the 3-Torus. Satisfying a divergence-free vector field and periodic boundary conditions respectively with a general spatio-temporal forcing term which is smooth and spatially periodic, the existence of solutions which have finite time singularities can occur starting with the first derivative and higher with respect to time. The existence of a subspace of the solution space where v<sub>3</sub> is continuous and {C, y<sub>1</sub>, y<sub>1</sub><sup>2</sup>}, is linearly independent in the additive argument of the solution in terms of the Lambert W function, (y<sub>1</sub><sup>2</sup>=y<sub>2</sub>, C∈R) together with the condition v<sub>2</sub>=-2y<sub>1</sub>v<sub>1</sub>. On this subspace, the Biot Savart Law holds exactly [see Section 2 (Equation (13))]. Also on this subspace, an expression X (part of PNS equations) vanishes which contains all the expressions in derivatives of v<sub>1</sub> and v<sub>2</sub> and the forcing terms in the plane which are related as with the cancellation of all such terms in governing PDE. The y<sub>3</sub> component forcing term is arbitrarily small in ε ball where Weierstrass P functions touch the center of the ball both for inviscid and viscous cases. As a result, a significant simplification occurs with a v<sub>3 </sub>only governing PDE resulting. With viscosity present as v changes from zero to the fully viscous case at v =1 the solution for v<sub>3</sub> reaches a peak in the third component y<sub>3</sub>. Consequently, there exists a dipole which is not centered at the center of the cell of the Lattice. Hence since the dipole by definition has an equal in magnitude positive and negative peak in y<sub>3</sub>, then the dipole Riemann cut-off surface is covered by a closed surface which is the sphere and where a given cell of dimensions [-1, 1]<sup>3</sup> is circumscribed on a sphere of radius 1. For such a closed surface containing a dipole it necessarily follows that the flux at the surface of the sphere of v<sub>3</sub> wrt to surface normal n is zero including at the points where the surface of sphere touches the cube walls. At the finite time singularity on the sphere a rotation boundary condition is deduced. It is shown that v<sub>3</sub> is spatially finite on the Riemann Sphere and the forcing is oscillatory in y<sub>3</sub> component if the velocity v3</sub> is. It is true that . A boundary condition on the sphere shows the rotation of a sphere of viscous fluid. Finally on the sphere a solution for v3</sub> is obtained which is proven to be Hölder continuous and it is shown that it is possible to extend Hölder continuity on the sphere uniquely to all of the interior of the ball.
基金supported by the National Natural Science Foundation of China(11232002)the Ph.D.Student Foundation of Chinese Ministry of Education(30400002011105001)
文摘The forward flight of a model butterfly was stud- ied by simulation using the equations of motion coupled with the Navier-Stokes equations. The model butterfly moved under the action of aerodynamic and gravitational forces, where the aerodynamic forces were generated by flapping wings which moved with the body, allowing the body os- cillations of the model butterfly to be simulated. The main results are as follows: (1) The aerodynamic force produced by the wings is approximately perpendicular to the long-axis of body and is much larger in the downstroke than in the up- stroke. In the downstroke the body pitch angle is small and the large aerodynamic force points up and slightly backward, giving the weight-supporting vertical force and a small neg- ative horizontal force, whilst in the upstroke, the body an- gle is large and the relatively small aerodynamic force points forward and slightly downward, giving a positive horizon- tal force which overcomes the body drag and the negative horizontal force generated in the downstroke. (2) Pitching oscillation of the butterfly body plays an equivalent role of the wing-rotation of many other insects. (3) The body-mass- specific power of the model butterfly is 33.3 W/kg, not very different from that of many other insects, e.g., fruitflies and dragonflies.
基金supported by the National Natural Science Foundation of China (10732030) and the 111 Project (B07009)
文摘In the present paper, the longitudinal dynamic flight stability properties of two model insects are predicted by an approximate theory and computed by numerical sim- ulation. The theory is based on the averaged model (which assumes that the frequency of wingbeat is sufficiently higher than that of the body motion, so that the flapping wings' degrees of freedom relative to the body can be dropped and the wings can be replaced by wingbeat-cycle-average forces and moments); the simulation solves the complete equations of motion coupled with the Navier-Stokes equations. Comparison between the theory and the simulation provides a test to the validity of the assumptions in the theory. One of the insects is a model dronefly which has relatively high wingbeat frequency (164 Hz) and the other is a model hawkmoth which has relatively low wingbeat frequency (26 Hz). The results show that the averaged model is valid for the hawkmoth as well as for the dronefly. Since the wingbeat frequency of the hawkmoth is relatively low (the characteristic times of the natural modes of motion of the body divided by wingbeat period are relatively large) compared with many other insects, that the theory based on the averaged model is valid for the hawkmoth means that it could be valid for many insects.
文摘A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.