Feng Shi-De

Wind-driven ocean circulation in shallow water lattice Boltzmann model

A lattice Boltzmann (LB) model with overall second-order accuracy is applied to the 1.5-layer shallow water equation for a wind-driven double-gyre ocean circulation. By introducing the second-order integral approximation for the collision operator, the model becomes fully explicit. In this case, any iterative technique is not needed. The Coriolis force and other external forces are included in the model with second-order accuracy, which is consistent with the discretized accuracy of the LB equation. The numerical results show correct physics of the ocean circulation driven by the double-gyre wind stress with different Reynolds numbers and different spatial resolutions. An intrinsic low-frequency variability of the shallow water model is also found. The wind-driven ocean circulation exhibits subannual and interannual oscillations, which are comparable to those of models in which the conventional numerical methods are used.

Lattice Boltzmann equation model in the Coriolis field

In a large-scale field of rotational fluid, various unintelligible and surprising dynamic phenomena are produced due to the effect of the Coriolis force. The lattice Boltzmann equation (LBE) model in the Coriolis field is developed based on previous works.[1-4] Geophysical fluid dynamics equations are derived from the model. Numerical simulations have been made on an ideal atmospheric circulation of the Northern Hemisphere by using the model and they reproduce the Rossby wave motion well. Hence the applicability of the model is verified in both theory and experiment.

Numerical analysis of fluid flow through a cylinder array using a lattice Boltzmann model

In this paper we present a detailed computational study of an incompressible Newtonian fluid flow across a periodic array of two-dimensional cylinders which is a simplest non-trivial representation of a porous media. A two-dimensional Lattice Boltzmann Method is used to solve the governing Navier–Stokes equation taking into account of viscous dissipation effects and influence of nonlinear fluid drag. Both the flow fields and the Darcy–Forchheimer drag coefficient as a function of the solid volume fraction are calculated for a wide range of flow Reynolds numbers. The predictions were compared with the results from conventional numerical and empirical models for verification. Apart from confirming that inertial effects can cause a significant deviation from Darcys law for large velocities the results also show that the characteristics of the vorticity field vary considerably as the Reynolds number increases, which will have major implications to the transport of passive particulate substances within the pores and their removal rate.

SOME PROGRESS IN THE LATTICE BOLTZMANN MODEL

A lattice Boltzmann equation model has been developed by using the equilibrium distribution function of the Maxwell-Boltzmann-like form, which is third order in fluid velocity uα. The criteria of energy conservation between the macroscopic physical quantities and the microscopic particles are introduced into the model, thus the thermal hydrodynamic equations containing the effect of buoyancy force can be recovered in terms of the Taylor and Chapman-Enskog asymptotic expansion methods. The two-dimensional thermal convection phenomena in a square cavity and between two concentric cylinders have been calculated by implementing a heat flux boundary condition. Both numerical results are in good agreement with the conventional numerical results.

A Climatology of the Southwest Vortex during 1979-2008

Abstract Using a new vortex detection and tracing method, a dataset of the Southwest Vortex (SWV) is established based on Japanese 25-year Reanalysis (JRA-25) reanalysis data during 1979–2008. The spatiotemporal features of the SWV are derived from the dataset. In comparison to other seasons, summer yields the least SWVs, but with the highest probability that they will migrate from their region of origin. SWVs mostly emerge in the southwest of the Sichuan Basin and the southeast of the Tibetan Plateau. Migratory SWVs mainly move along either an eastward or southeastward path. Detailed composite analysis of warm-season SWVs shows that the subtropical high is a key factor in determining the direction of migratory SWVs. Furthermore, the steering wind at 700 hPa dominates the moving direction of migratory SWVs. Potential stability diagnosed by pseudo-equivalent potential temperature θse is of certain significance for the evolution and movement of SWVs. On the other hand, migratory SWVs possess relatively greater strength than stationary SWVs, due to a stronger low-level jet with enhanced baroclinicity and moisture transport providing more energy to support the growth of SWVs along their paths of movement.

A Lattice Boltzmann Model for Two-Dimensional Magnetohydrodynamics

A lattice Boltzmann model (LBM) has been developed for simulating magnetohydrodynamics (MHD) along the line of Dellar [J. Comput. Phys. 179 (2002) 95]. In this model the magnetic field is presented by a vector valued magnetic distribution function which obeys a vector Boltzmann equation. The truncated error of the equilibrium distribution in the present model is up to order O(u4) in velocity u rather than the usual O(u3). For verification, the model is applied to solve the shock tube problem and the main features of the flow predicted by the model are found to compare well with the corresponding results obtained with high-order semi-discrete schemes [J. Comput. Phys. 201 (2004) 261]. The numerical experiments have also shown that the present LBM model with the equilibrium distribution truncated at O(u4) performs much better in terms of numerical stability than those truncated at O(u3).

Corrections to the collision term in the BGK Boltzmann equation

With the discrete method of the hexagonal cell and three different velocities of particle population in each cell, a two-dimensional lattice Boltzmann model is developed in this paper. The collision operator in the Boltzmann equation is expanded to fourth order using the Taylor expansion. With this model, good results have been obtained from the numerical simulation of the reflection phenomenon of the shock wave on the surface of an obstacle and the numerical stability is also good. Thus the applicability of the D2Q19 model is verified.

A STATISTICAL ANALYSIS ON THE DIFFERENT FEATURES OF TC EXTRATROPICAL TRANSITION OVER LAND AND OCEAN

The extratropical transitions(ETs)of tropical cyclones(TCs)over China and the ocean east to 150°E are investigated by the use of best-track data and JRA-25 reanalysis spanning 1979–2008.The ET events occurring north of 25°N and in the warm season(from May to October)are extracted from the reanalysis to emphasize the interaction between TC and midlatitude circulation.Statistical analysis shows that 18.5%of the warm-season TCs go through land ETs north of 25°N in the western North Pacific.And 20.5%of the ET events occur over the ocean east of 150°E.Most(62.2%)ET TCs over China gradually die out after ET,but more(70.7%)ocean ET cases have post-ET reintensification.The evolutions in cyclone phase space and the composite fields for land and ocean ETs,as well as the ET cases with and without post-ET reintensification,are further analyzed.It is found that most TCs with ET over China and those without post-ET reintensification evolve along the typical ET phase path as follows:emergence of thermal asymmetry→losing upper-level warm core→losing lower-level cold core→evolving as extratropical cyclone.The TCs undergoing ETs over ocean and those with post-ET reintensification form a high-level cold core before the ET onset.The TCs with land ET have long distance between the landing TC and a high-level trough.That makes the TC maintain more tropical features and isolates the TC flow from the upstream and downstream jets of the midlatitude trough.The structure of circulation leads to weak development of baroclinicity in land ET.On the contrary,shorter distance between ocean TC and high-level trough makes the high-level trough absorb the TC absolutely.Under that baroclinicity-favorable environment,strong cold advection makes the TC lose its high-level warm core before ET onset.The composite fields confirm that the TC with ocean ET has stronger baroclinic features.Generally,the TC at land ET onset is located to the south of the ridge of the subtropical high,which tends to prevent the TCs from interacting with midlatitude circulation.But for the ocean ET,the situation is just the opposite.Similar analyses are also carried out for the TCs with and without post-ET reintensification over both land and ocean east of 150°E.The results further prove that the TC with stronger baroclinic characteristics,especially in the circumstance favorable to its interaction with high-level midlatitude systems,has more opportunity to reintensify as an extratropical cyclone after ET.