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.展开更多
The moist potential vorticity (MPV) equation is derived from complete atmospheric equations including the effect of mass forcing, with which the theory of Up-sliding Slantwise Vorticity Development, (USVD) is proposed...The moist potential vorticity (MPV) equation is derived from complete atmospheric equations including the effect of mass forcing, with which the theory of Up-sliding Slantwise Vorticity Development, (USVD) is proposed based on the theory of Slantwise Vorticity Development (SVD). When an air parcel slides up along a slantwise isentropic surface, its vertical component of relative vorticity will develop, and the steeper the isentropic surface is, the more violent the development will he. From the definition of MPV and the MPV equation produced here in, a complete vorticity equation is then put forward with mass forcing, which explicitly includes the effects of both internal forcings, such as variations of stability, baroclinicity, and vertical shear of horizontal wind, arid external forcings, such as diabatic heating, friction, and mass forcing. When isentropic surfaces are flat, the complete vorticity equation matches its traditional counterpart. The physical interpretations of some of the items which are included in the complete- vorticity equation but not in the traditional one are studied with a simplified model of the Changjiang-Huaihe Meiyu front. A 60-h simulation is then performed to reproduce a torrential rain event in the Changjiang-Huaihe region and the output of the model is studied qualitatively based on the theory of USVD. The result shows that the conditions of the theory of USVD are easily satisfied immediately in front of mesoscale rainstorms in the downwind direction, that is, the theory of USVD is important to the development and movement of these kinds of systems.展开更多
We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditi...We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition (without the incompressibility condition), which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n(n ≥ 3) with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.展开更多
For the potential vorticity (PV) invariant, there is a PV-based complete-form vorticity equation, which we use heuris- tically in the present paper to answer the following question: for the Ertel-Rossby invariant ...For the potential vorticity (PV) invariant, there is a PV-based complete-form vorticity equation, which we use heuris- tically in the present paper to answer the following question: for the Ertel-Rossby invariant (ERI), is there a corresponding vorticity tendency equation? Such an ERI-based thermally-coupled vorticity equation is derived and discussed in detail in this study. From the obtained new vorticity equation, the vertical vorticity change is constrained by the vertical velocity term, the term associated with the slope of the generalized momentum surface, the term related to the horizontal vorticity change, and the baroclinic or solenoid term. It explicitly includes both the dynamical and thermodynamic factors' influence on the vorticity change. For the ERI itself, besides the traditional PV term, the ERI also includes the moisture factor, which is excluded in dry ERI, and the term related to the gradients of pressure, kinetic energy, and potential energy that reflects the fast-manifold property. Therefore, it is more complete to describe the fast motions off the slow manifold for severe weather than the PV term. These advantages are naturally handed on and inherited by the ERI-based thermally-coupled vorticity equation. Then the ERI-based thermally-coupled vorticity equation is further transformed and compared with the traditional vorticity equation. The main difference between the two equations is the term which describes the contribution of the solenoid term to the vertical vorticity development. In a barotropic flow, the solenoid term disappears, then the ERI-based thermally-coupled vorticity equation can regress to the traditional vorticity equation.展开更多
This paper investigates an important high-dimensional model in the atmospheric and oceanic dynamics-(3+1)- dimensional nonlinear baroclinic potential vorticity equation by the classical Lie group method. Its symmet...This paper investigates an important high-dimensional model in the atmospheric and oceanic dynamics-(3+1)- dimensional nonlinear baroclinic potential vorticity equation by the classical Lie group method. Its symmetry algebra, symmetry group and group-invariant solutions are analysed. Otherwise, some exact explicit solutions are obtained from the corresponding (2+1)-dimensional equation, the inviscid barotropic nondivergent vorticy equation. To show the properties and characters of these solutions, some plots as well as their possible physical meanings of the atmospheric circulation are given out.展开更多
The momentum advection vorticity equation in the form of cross multiplication is introduced, in which the divergence term in the classic vorticity equation does not appear explicitly. This equation includes the rotati...The momentum advection vorticity equation in the form of cross multiplication is introduced, in which the divergence term in the classic vorticity equation does not appear explicitly. This equation includes the rotation effect of the horizontal wind advection, which are not explicitly included in the classic vortieity equation. The vorticity and its tendency of Typhoon Fung-Wong (0808) that occurred in July 2008 are analyzed. The computed results show that the rotation effect of the advection of the horizontal wind is a leading factor in determining the change of vertical vorticity for Fung-Wong during its life cycle, especially in the period leading up to landfall. The advection term represents the tendency variation of the vertical vortieity, and the positive-value region of the vertical vorticity tendency is almost in accord with the track of Fung-Wong, which may be taken as a factor to locate the key observational region of Fung-Wong. The equation provides a supplementary diagnostic tool for the systems related with strong advection of horizontal wind.展开更多
In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability...In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.展开更多
A Fourier-Chebyshev pseudospectral scheme is proposed for three-dimensionalvorticily equation with unilaterally periodic boundary condition. The generalized stability and convergence are analysed. The numerical result...A Fourier-Chebyshev pseudospectral scheme is proposed for three-dimensionalvorticily equation with unilaterally periodic boundary condition. The generalized stability and convergence are analysed. The numerical results are presented.展开更多
By means of the reductive perturbation method, three types of generalized (2+l)-dimensional Kadomtsev- Petviashvili (KP) equations are derived from the baroclinic potential vorticity (BPV) equation, including t...By means of the reductive perturbation method, three types of generalized (2+l)-dimensional Kadomtsev- Petviashvili (KP) equations are derived from the baroclinic potential vorticity (BPV) equation, including the modified KP (mKP) equation, standard KP equation and cylindrical KP (cKP) equation. Then some solutions of generalized cKP and KP equations with certain conditions are given directly and a relationship between the generalized mKP equation and the mKP equation is established by the symmetry group direct method proposed by Lou et al. From the relationship and the solutions of the mKP equation, some solutions of the generalized mKP equation can be obtained. Furthermore, some approximate solutions of the baroclinic potential vorticity equation are derived from three types of generalized KP equations.展开更多
Applying the classical Lie symmetry approach to the barotropic and quasi-geostrophic potential vorticity equation without forcing and dissipation on a β-plane channel, we find a new symmetry, which is not presented i...Applying the classical Lie symmetry approach to the barotropic and quasi-geostrophic potential vorticity equation without forcing and dissipation on a β-plane channel, we find a new symmetry, which is not presented in a previous work [F. Huang, Commun. Theor. Phys. (Beijing, China) 42 (2004) 903]. A general finite transformation group is obtained based on the full Lie point symmetry, Furthermore, two new types of similarity reduction solutions are obtained.展开更多
Recently F. Huang [Commun. Theor. Phys. 42 (2004) 903] and X. Tang and P.K. Shukla [Commun. Theor. Phys. 49 (2008) 229] investigated symmetry properties of the barotropic potential vorticity equation without forci...Recently F. Huang [Commun. Theor. Phys. 42 (2004) 903] and X. Tang and P.K. Shukla [Commun. Theor. Phys. 49 (2008) 229] investigated symmetry properties of the barotropic potential vorticity equation without forcing and dissipation on the beta-plane. This equation is governed by two dimensionless parameters, F and β, representing the ratio of the characteristic length scale to the Rossby radius of deformation and the variation of earth' angular rotation, respectively. In the present paper it is shown that in the case F ≠ 0 there exists a well-defined point transformation to set β = 0. The classification of one- and two-dimensional Lie subalgebras of the Lie symmetry algebra of the potential vorticity equation is given for the parameter combination F ≠ 0 and β = 0. Based upon this classification, distinct classes of group-invariant solutions are obtained and extended to the case β ≠0.展开更多
By the Backlund transformation method, an important (2+1)-dimensional nonlinear barotropie and quasigeostrophic potential vorticity (BQGPV) equation is investigated. Some simple special Backlund transformation th...By the Backlund transformation method, an important (2+1)-dimensional nonlinear barotropie and quasigeostrophic potential vorticity (BQGPV) equation is investigated. Some simple special Backlund transformation theorems are proposed and used to get explicit solutions of the BQGPV equation. Furthermore, all solutions of a second order linear ordinary differential equation including an arbitrary function can be used to construct explicit solutions of the (2+1)-dimensional BQGPV equation. Some figures are also given out to describe these solutions.展开更多
A Fourier pseudospectral-finue difference scheme is proposed for three-dimensional vorticity. equalion with unilalerally periodic boundary. condilion. Thegeneralized stability. and conrergence are analyzed. The nunier...A Fourier pseudospectral-finue difference scheme is proposed for three-dimensional vorticity. equalion with unilalerally periodic boundary. condilion. Thegeneralized stability. and conrergence are analyzed. The nunierical results show. theadvantage of this method.展开更多
The incompressible Navier Stokes equations are solved via variables of vorticity and velocity. Firstly, a rigorous variational framework with the equivalence between the velocity pressure and the vorticity velocity fo...The incompressible Navier Stokes equations are solved via variables of vorticity and velocity. Firstly, a rigorous variational framework with the equivalence between the velocity pressure and the vorticity velocity formulations is presented in a Lipschitz domain. Next, a class of Galerkin finite element approximations of the corresponding variational form is introduced, and a convergence analysis is given for the Stokes problem. Finally, an iterative finite element solver for the Navier Stokes problem is proposed.展开更多
To investigate the multiscale interaction characteristics of Landfall Typhoon Lekima(2019),this study analyzed the characteristics of the different scale vortex structure and interactions among different scales based ...To investigate the multiscale interaction characteristics of Landfall Typhoon Lekima(2019),this study analyzed the characteristics of the different scale vortex structure and interactions among different scales based on vorticity equation diagnosis.The analysis is based on the simulation results of the WRF model which has been thoroughly verified.The main results are as follows:the original vorticity dominated by the meso-αscale vorticity increases with height and then decreases,with maximum vorticity distributed at 900 hPa.The meso-βscale vorticity varies significantly with altitude,while the meso-γscale vorticityfield exhibits obvious positive vorticity below 850 hPa.The meso-αscale vorticity tendency primarily maintains negative,contributing significantly to the overall reduction in the original vorticityfield over time.The increase in mid-to-upper-level(above 550 hPa)original vorticity is mainly related to the variations in the meso-βand meso-γscale vorticityfields.The original vorticity dominated by the meso-αscale vorticity increases with height and then decreases,and the whole layer vorticity decreases over time.The meso-βscale vorticity varies significantly with altitude and time,while the meso-γscale vorticityfield consistently exhibits significant positive vorticity below 850 hPa.The vorticity equation diagnosis revealed that the primary source terms of the vorticity tendencies are the twisting and stretching terms,and the main sink terms being horizontal and vertical vorticity transport terms below 900 hPa.The source terms and sink terms exchange above 850 hPa.Scale separation results show that the primary contributions of all impact factors originate from the meso-αand meso-γscalefields(accounting for over 80%of the total),with the contribution of the meso-αscale being less than that of the meso-γscale and a notable contribution over 35.5%of the interactions between different scales.展开更多
The horizontal vorticity equation used in this study was obtained using the equations of motion in the pressure coordinate system without considering friction, to reveal its relationship with vertical shear. By diagno...The horizontal vorticity equation used in this study was obtained using the equations of motion in the pressure coordinate system without considering friction, to reveal its relationship with vertical shear. By diagnostically analyzing each term in the horizontal vorticity equation during a squall line process that occurred on 19 June 2010, we found that the non-thermal wind term had a negative contribution to the local change of upward movement in the low-level atmosphere, and that its impact changed gradually from negative to positive with altitude, which could influence upward movement in the mid-and upper-level atmosphere greatly. The contribution of upward vertical transport to vertical movement was the largest in the low-level atmosphere, but had negative contribution to the upper-level atmosphere. These features were most evident in the development stage of the squall line. Based on analysis of convection cells along a squall line, we found that in the process of cell development diabatic heating caused the subsidence of constant potential temperature surface and non-geostrophic motion, which then triggered strong convergence of horizontal acceleration in the mid-level atmosphere and divergence of horizontal acceleration in the upper-level atmosphere. These changes of horizontal wind field could cause a counterclockwise increment of the horizontal vorticity around the warm cell, which then generated an increase of upward movement. This was the main reason why the non-thermal wind term had the largest contribution to the strengthening of upward movement in the mid-and upper-level atmosphere. The vertical transport of large value of horizontal vorticity was the key to trigger convection in this squall line process.展开更多
The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not...The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi (LBB) condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1(Ω).展开更多
The present work provides a novel method for calculating vertical velocity based on continuity equations in a pressure coordinate system.The method overcomes the disadvantage of accumulation of calculating errors of h...The present work provides a novel method for calculating vertical velocity based on continuity equations in a pressure coordinate system.The method overcomes the disadvantage of accumulation of calculating errors of horizontal divergence in current kinematics methods during the integration for calculating vertical velocity,and consequently avoids its subsequent correction.In addition,through modifications of the continuity equations,it shows that the vorticity of the vertical shear vector(VVSV) is proportional to-ω,the vertical velocity in p coordinates.Furthermore,if the change of ω in the horizontal direction is neglected,the vorticity of the horizontal vorticity vector is proportional to-ω.When ω is under a fluctuating state in the vertical direction,the updraft occurs when the vector of horizontal vorticity rotates counterclockwise;the downdraft occurs when rotating clockwise.The validation result indicates that the present method is generally better than the vertical velocity calculated by the ω equation using the wet Q-vector divergence as a forcing term,and the vertical velocity calculated by utilizing the kinematics method is followed by the O'Brien method for correction.The plus-minus sign of the vertical velocity obtained with this method is not correlated with the intensity of d BZ,but the absolute error increases when d BZ is >=40.This method demonstrates that it is a good reflection of the direction of the vertical velocity.展开更多
Momentum balances are used to derive the Kutta-Joukowsky equation for an infinite cascade of aerofoils and an isolated aerofoil. These derivations are simpler than those based on the Blasius theorem or more complex un...Momentum balances are used to derive the Kutta-Joukowsky equation for an infinite cascade of aerofoils and an isolated aerofoil. These derivations are simpler than those based on the Blasius theorem or more complex unsteady control volumes, and show the close relationship between a single aerofoil and an infinite cascade. The modification of lift due to the presence of another lifting body is similarly derived for a wing in ground effect, a biplane, and tandem aerofoils. The results are identical to those derived from the vector form of the Kutta-Joukowsky equation.展开更多
Taking an extratropical cyclone that produced extreme precipitation as the research object,this paper calculates the contribution of condensation latent heat release(LHR)to relative vorticity tendency based on the com...Taking an extratropical cyclone that produced extreme precipitation as the research object,this paper calculates the contribution of condensation latent heat release(LHR)to relative vorticity tendency based on the complete-form vertical vorticity tendency equation.The results show that the heating rate of convectional condensation LHR can reach up to about 40 times that of stable condensation LHR.Both the stable and convectional heating centers are higher than 700 hPa,which would cause∂Q/∂z>0 and a positive vorticity source in the lower troposphere.The vertical gradient of stable condensation LHR contributes little to the growth of relative vorticity,while the relative vorticity tendency associated with the vertical gradient of convectional condensation LHR can be an order of magnitude higher than the former.The positive vorticity source is always located right below the latent heating center,and its maximum value can always be found in the lower troposphere.Convectional LHR is the primary factor for cyclone development from the perspective of diabatic heating.The horizontal gradient of total condensation LHR can contribute about 65%of the actual vorticity growth,but the effect of the vertical gradient of convectional condensation(LHR)can reach twice as much.The adiabatic heating from LHR can cause vorticity tendency directly.However,it can also change the vertical and horizontal gradient of potential temperature,which can further induce vorticity tendency.展开更多
文摘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.
基金supported by the Chinese Academy of Sciences Program of Well-Known 0verseas Chinese Scholars under Dr.Samuel Shen and the Innovation Project of the Chinese Academy of Sciences under Grant No.KZCX3-SW-213by the National Natural Science foundation of China under Grant Nos.40023001 and 40135020.
文摘The moist potential vorticity (MPV) equation is derived from complete atmospheric equations including the effect of mass forcing, with which the theory of Up-sliding Slantwise Vorticity Development, (USVD) is proposed based on the theory of Slantwise Vorticity Development (SVD). When an air parcel slides up along a slantwise isentropic surface, its vertical component of relative vorticity will develop, and the steeper the isentropic surface is, the more violent the development will he. From the definition of MPV and the MPV equation produced here in, a complete vorticity equation is then put forward with mass forcing, which explicitly includes the effects of both internal forcings, such as variations of stability, baroclinicity, and vertical shear of horizontal wind, arid external forcings, such as diabatic heating, friction, and mass forcing. When isentropic surfaces are flat, the complete vorticity equation matches its traditional counterpart. The physical interpretations of some of the items which are included in the complete- vorticity equation but not in the traditional one are studied with a simplified model of the Changjiang-Huaihe Meiyu front. A 60-h simulation is then performed to reproduce a torrential rain event in the Changjiang-Huaihe region and the output of the model is studied qualitatively based on the theory of USVD. The result shows that the conditions of the theory of USVD are easily satisfied immediately in front of mesoscale rainstorms in the downwind direction, that is, the theory of USVD is important to the development and movement of these kinds of systems.
基金supported in part by the National Science Foundation under Grants DMS-0807551, DMS-0720925, and DMS-0505473the Natural Science Foundationof China (10728101)supported in part by EPSRC grant EP/F029578/1
文摘We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition (without the incompressibility condition), which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n(n ≥ 3) with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.
基金supported by the National Basic Research Program of China(Grant No.2013CB430105)the Key Research Program of the Chinese Academy of Sciences(Grant No.KZZD-EW-05-01)+2 种基金the National Natural Science Foundation of China(Grant Nos.41375054 and 41375052)the Special Scientific Research Fund of the Meteorological Public Welfare of the Ministry of Sciences and Technology,China(Grant No.GYHY201406003)the Opening Foundation of the State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences(Grant No.2012LASW-B02)
文摘For the potential vorticity (PV) invariant, there is a PV-based complete-form vorticity equation, which we use heuris- tically in the present paper to answer the following question: for the Ertel-Rossby invariant (ERI), is there a corresponding vorticity tendency equation? Such an ERI-based thermally-coupled vorticity equation is derived and discussed in detail in this study. From the obtained new vorticity equation, the vertical vorticity change is constrained by the vertical velocity term, the term associated with the slope of the generalized momentum surface, the term related to the horizontal vorticity change, and the baroclinic or solenoid term. It explicitly includes both the dynamical and thermodynamic factors' influence on the vorticity change. For the ERI itself, besides the traditional PV term, the ERI also includes the moisture factor, which is excluded in dry ERI, and the term related to the gradients of pressure, kinetic energy, and potential energy that reflects the fast-manifold property. Therefore, it is more complete to describe the fast motions off the slow manifold for severe weather than the PV term. These advantages are naturally handed on and inherited by the ERI-based thermally-coupled vorticity equation. Then the ERI-based thermally-coupled vorticity equation is further transformed and compared with the traditional vorticity equation. The main difference between the two equations is the term which describes the contribution of the solenoid term to the vertical vorticity development. In a barotropic flow, the solenoid term disappears, then the ERI-based thermally-coupled vorticity equation can regress to the traditional vorticity equation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10735030,90718041 and 40975038)Shanghai Leading Academic Discipline Project(Grant No.B412)Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0734)
文摘This paper investigates an important high-dimensional model in the atmospheric and oceanic dynamics-(3+1)- dimensional nonlinear baroclinic potential vorticity equation by the classical Lie group method. Its symmetry algebra, symmetry group and group-invariant solutions are analysed. Otherwise, some exact explicit solutions are obtained from the corresponding (2+1)-dimensional equation, the inviscid barotropic nondivergent vorticy equation. To show the properties and characters of these solutions, some plots as well as their possible physical meanings of the atmospheric circulation are given out.
基金projects of the Ministry of Sciences and Technology of the People’s Republic of China (GYHY200906004 GYHY200706020)+1 种基金project of the Natural Science Foundation of China (40975034)project of State Key Laboratory of Severe Weather (2008LASW-A01)
文摘The momentum advection vorticity equation in the form of cross multiplication is introduced, in which the divergence term in the classic vorticity equation does not appear explicitly. This equation includes the rotation effect of the horizontal wind advection, which are not explicitly included in the classic vortieity equation. The vorticity and its tendency of Typhoon Fung-Wong (0808) that occurred in July 2008 are analyzed. The computed results show that the rotation effect of the advection of the horizontal wind is a leading factor in determining the change of vertical vorticity for Fung-Wong during its life cycle, especially in the period leading up to landfall. The advection term represents the tendency variation of the vertical vortieity, and the positive-value region of the vertical vorticity tendency is almost in accord with the track of Fung-Wong, which may be taken as a factor to locate the key observational region of Fung-Wong. The equation provides a supplementary diagnostic tool for the systems related with strong advection of horizontal wind.
文摘In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.
文摘A Fourier-Chebyshev pseudospectral scheme is proposed for three-dimensionalvorticily equation with unilaterally periodic boundary condition. The generalized stability and convergence are analysed. The numerical results are presented.
基金supported by National Natural Science Foundation of China (Grant Nos.10735030 and 40775042)Ningbo Natural Science Foundation (Grant No. 2008A610017)+1 种基金National Basic Research Program of China (973 Program) (Grant Nos. 2005CB422301 and 2007CB814800)K.C. Wong Magna Fund in Ningbo University
文摘By means of the reductive perturbation method, three types of generalized (2+l)-dimensional Kadomtsev- Petviashvili (KP) equations are derived from the baroclinic potential vorticity (BPV) equation, including the modified KP (mKP) equation, standard KP equation and cylindrical KP (cKP) equation. Then some solutions of generalized cKP and KP equations with certain conditions are given directly and a relationship between the generalized mKP equation and the mKP equation is established by the symmetry group direct method proposed by Lou et al. From the relationship and the solutions of the mKP equation, some solutions of the generalized mKP equation can be obtained. Furthermore, some approximate solutions of the baroclinic potential vorticity equation are derived from three types of generalized KP equations.
基金The project supported by the Alexander von Humboldt Foundationthe Youth Foundation of Shanghai Jiao Tong UniversityNational Natural Science Foundation of China under Grant No.10475055
文摘Applying the classical Lie symmetry approach to the barotropic and quasi-geostrophic potential vorticity equation without forcing and dissipation on a β-plane channel, we find a new symmetry, which is not presented in a previous work [F. Huang, Commun. Theor. Phys. (Beijing, China) 42 (2004) 903]. A general finite transformation group is obtained based on the full Lie point symmetry, Furthermore, two new types of similarity reduction solutions are obtained.
基金supported by the Austrian Science Fund (FWF),project P20632
文摘Recently F. Huang [Commun. Theor. Phys. 42 (2004) 903] and X. Tang and P.K. Shukla [Commun. Theor. Phys. 49 (2008) 229] investigated symmetry properties of the barotropic potential vorticity equation without forcing and dissipation on the beta-plane. This equation is governed by two dimensionless parameters, F and β, representing the ratio of the characteristic length scale to the Rossby radius of deformation and the variation of earth' angular rotation, respectively. In the present paper it is shown that in the case F ≠ 0 there exists a well-defined point transformation to set β = 0. The classification of one- and two-dimensional Lie subalgebras of the Lie symmetry algebra of the potential vorticity equation is given for the parameter combination F ≠ 0 and β = 0. Based upon this classification, distinct classes of group-invariant solutions are obtained and extended to the case β ≠0.
基金Supported by the National Natural Science Foundation of China under Grant Nos. 10735030, 90718041, and 40975038Shanghai Leading Academic Discipline Project under Grant No. B412Program for Changjiang Scholars and Innovative Research Team in University (IRT0734)
文摘By the Backlund transformation method, an important (2+1)-dimensional nonlinear barotropie and quasigeostrophic potential vorticity (BQGPV) equation is investigated. Some simple special Backlund transformation theorems are proposed and used to get explicit solutions of the BQGPV equation. Furthermore, all solutions of a second order linear ordinary differential equation including an arbitrary function can be used to construct explicit solutions of the (2+1)-dimensional BQGPV equation. Some figures are also given out to describe these solutions.
文摘A Fourier pseudospectral-finue difference scheme is proposed for three-dimensional vorticity. equalion with unilalerally periodic boundary. condilion. Thegeneralized stability. and conrergence are analyzed. The nunierical results show. theadvantage of this method.
文摘The incompressible Navier Stokes equations are solved via variables of vorticity and velocity. Firstly, a rigorous variational framework with the equivalence between the velocity pressure and the vorticity velocity formulations is presented in a Lipschitz domain. Next, a class of Galerkin finite element approximations of the corresponding variational form is introduced, and a convergence analysis is given for the Stokes problem. Finally, an iterative finite element solver for the Navier Stokes problem is proposed.
基金funded by the National Natural Science Foundation of China(grant U2142206).W。
文摘To investigate the multiscale interaction characteristics of Landfall Typhoon Lekima(2019),this study analyzed the characteristics of the different scale vortex structure and interactions among different scales based on vorticity equation diagnosis.The analysis is based on the simulation results of the WRF model which has been thoroughly verified.The main results are as follows:the original vorticity dominated by the meso-αscale vorticity increases with height and then decreases,with maximum vorticity distributed at 900 hPa.The meso-βscale vorticity varies significantly with altitude,while the meso-γscale vorticityfield exhibits obvious positive vorticity below 850 hPa.The meso-αscale vorticity tendency primarily maintains negative,contributing significantly to the overall reduction in the original vorticityfield over time.The increase in mid-to-upper-level(above 550 hPa)original vorticity is mainly related to the variations in the meso-βand meso-γscale vorticityfields.The original vorticity dominated by the meso-αscale vorticity increases with height and then decreases,and the whole layer vorticity decreases over time.The meso-βscale vorticity varies significantly with altitude and time,while the meso-γscale vorticityfield consistently exhibits significant positive vorticity below 850 hPa.The vorticity equation diagnosis revealed that the primary source terms of the vorticity tendencies are the twisting and stretching terms,and the main sink terms being horizontal and vertical vorticity transport terms below 900 hPa.The source terms and sink terms exchange above 850 hPa.Scale separation results show that the primary contributions of all impact factors originate from the meso-αand meso-γscalefields(accounting for over 80%of the total),with the contribution of the meso-αscale being less than that of the meso-γscale and a notable contribution over 35.5%of the interactions between different scales.
基金National Key Basic Research Development Program“973”(2013CB430103)State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences(2015LASW-A07)National Natural Science Funding(41375058,41530427)
文摘The horizontal vorticity equation used in this study was obtained using the equations of motion in the pressure coordinate system without considering friction, to reveal its relationship with vertical shear. By diagnostically analyzing each term in the horizontal vorticity equation during a squall line process that occurred on 19 June 2010, we found that the non-thermal wind term had a negative contribution to the local change of upward movement in the low-level atmosphere, and that its impact changed gradually from negative to positive with altitude, which could influence upward movement in the mid-and upper-level atmosphere greatly. The contribution of upward vertical transport to vertical movement was the largest in the low-level atmosphere, but had negative contribution to the upper-level atmosphere. These features were most evident in the development stage of the squall line. Based on analysis of convection cells along a squall line, we found that in the process of cell development diabatic heating caused the subsidence of constant potential temperature surface and non-geostrophic motion, which then triggered strong convergence of horizontal acceleration in the mid-level atmosphere and divergence of horizontal acceleration in the upper-level atmosphere. These changes of horizontal wind field could cause a counterclockwise increment of the horizontal vorticity around the warm cell, which then generated an increase of upward movement. This was the main reason why the non-thermal wind term had the largest contribution to the strengthening of upward movement in the mid-and upper-level atmosphere. The vertical transport of large value of horizontal vorticity was the key to trigger convection in this squall line process.
基金supported by the National Basic Research Program of China (2005CB321701)NSF of mathematics research special fund of Hebei Province(08M005)
文摘The purpose of this article is to develop and analyze least-squares approximations for the incompressible magnetohydrodynamic equations. The major advantage of the least-squares finite element method is that it is not subjected to the so-called Ladyzhenskaya-Babuska-Brezzi (LBB) condition. The authors employ least-squares functionals which involve a discrete inner product which is related to the inner product in H^-1(Ω).
基金National Key Basic Research Development Program"973"(2013CB430103,2009CB421503)National Natural Science Funding(41375058,41530427)State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences(2015LASW-A07)
文摘The present work provides a novel method for calculating vertical velocity based on continuity equations in a pressure coordinate system.The method overcomes the disadvantage of accumulation of calculating errors of horizontal divergence in current kinematics methods during the integration for calculating vertical velocity,and consequently avoids its subsequent correction.In addition,through modifications of the continuity equations,it shows that the vorticity of the vertical shear vector(VVSV) is proportional to-ω,the vertical velocity in p coordinates.Furthermore,if the change of ω in the horizontal direction is neglected,the vorticity of the horizontal vorticity vector is proportional to-ω.When ω is under a fluctuating state in the vertical direction,the updraft occurs when the vector of horizontal vorticity rotates counterclockwise;the downdraft occurs when rotating clockwise.The validation result indicates that the present method is generally better than the vertical velocity calculated by the ω equation using the wet Q-vector divergence as a forcing term,and the vertical velocity calculated by utilizing the kinematics method is followed by the O'Brien method for correction.The plus-minus sign of the vertical velocity obtained with this method is not correlated with the intensity of d BZ,but the absolute error increases when d BZ is >=40.This method demonstrates that it is a good reflection of the direction of the vertical velocity.
文摘Momentum balances are used to derive the Kutta-Joukowsky equation for an infinite cascade of aerofoils and an isolated aerofoil. These derivations are simpler than those based on the Blasius theorem or more complex unsteady control volumes, and show the close relationship between a single aerofoil and an infinite cascade. The modification of lift due to the presence of another lifting body is similarly derived for a wing in ground effect, a biplane, and tandem aerofoils. The results are identical to those derived from the vector form of the Kutta-Joukowsky equation.
基金This study was supported by the Natural Science Foundation of Jiangsu Province[grant number BK20161603]the National Natural Science Foundation of China[grant numbers 41575010 and 41575070]the China Meteorological Administration[grant number CMAYBY2018-028].
文摘Taking an extratropical cyclone that produced extreme precipitation as the research object,this paper calculates the contribution of condensation latent heat release(LHR)to relative vorticity tendency based on the complete-form vertical vorticity tendency equation.The results show that the heating rate of convectional condensation LHR can reach up to about 40 times that of stable condensation LHR.Both the stable and convectional heating centers are higher than 700 hPa,which would cause∂Q/∂z>0 and a positive vorticity source in the lower troposphere.The vertical gradient of stable condensation LHR contributes little to the growth of relative vorticity,while the relative vorticity tendency associated with the vertical gradient of convectional condensation LHR can be an order of magnitude higher than the former.The positive vorticity source is always located right below the latent heating center,and its maximum value can always be found in the lower troposphere.Convectional LHR is the primary factor for cyclone development from the perspective of diabatic heating.The horizontal gradient of total condensation LHR can contribute about 65%of the actual vorticity growth,but the effect of the vertical gradient of convectional condensation(LHR)can reach twice as much.The adiabatic heating from LHR can cause vorticity tendency directly.However,it can also change the vertical and horizontal gradient of potential temperature,which can further induce vorticity tendency.