Development of a Single-Column Model in RegCM4 and Its Preliminary Application for Evaluating PBL Schemes in Simulating the Dry Convection Boundary Layer

A single-column model(SCM)is developed in the regional climate model RegCM4.The evolution of a dry convection boundary layer(DCBL)is used to evaluate this SCM.Moreover,four planetary boundary layer(PBL)schemes,namely the Holtslag-Boville scheme(HB),Yonsei University scheme(YSU),and two University of Washington schemes(UW01,Grenier-Bretherton-Mc Caa scheme and UW09,Bretherton-Park scheme),are compared by using the SCM approach.A large-eddy simulation(LES)of the DCBL is performed as a benchmark to examine how well a PBL parameterization scheme reproduces the LES results,and several diagnostic outputs are compared to evaluate the schemes.The results show that the SCM is properly constructed.In general,with the DCBL case,the YSU scheme performs best for reproducing the LES results,which include well-mixed features and vertical sensible heat fluxes;the simulated wind speed,turbulent kinetic energy,entrainment flux,and height of the entrainment zone are all underestimated in the UW09;the UW01 has all those biases of the UW09 but larger,and the simulated potential temperature is not well mixed;the HB is the least skillful scheme,by which the PBL height,entrainment flux,height of the entrainment zone,and the vertical gradients within the mixed layer are all overestimated,and a inversion layer near the top of the surface layer is wrongly simulated.Although more cases and further testing are required,these simulations show encouraging results towards the use of this SCM framework for evaluating the simulated physical processes by the RegCM4.

ALGEBRAIC TURBUL ENCE MODEL WITH MEM ORY FOR COMPUTATION OF 3-D TURBULENT BOUNDARY LAYERS WITH VALIDATION

Additional equations were found based on experiments for an algebraic turbulence model to improve the prediction of the behavior of three dimensional turbulent boundary layers by taking account of the effects of pressure gradient and the historical variation of eddy viscosity, so the model is with memory. Numerical calculation by solving boundary layer equations was carried out for the five pressure driven three dimensional turbulent boundary layers developed on flat plates, swept wing, and prolate spheroid in symmetrical plane. Comparing the computational results with the experimental data, it is obvious that the prediction will be more accurate if the proposed closure equations are used, especially for the turbulent shear stresses.

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers:A Case Study for the Tonga Subducting Slab

For the purpose of investigating the influence of metastable olivine（MO） phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near～410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab：（i） slab pull force exerts above a depth of～230 km;（ii） MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;（iii） as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;（iv） a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.

ATMOSPHERIC BOUNDARY LAYER CONCEPT MODEL OF THE PEARL RIVER DELTA AND ITS APPLICATION

Based on the geographical circumstance, climate and the boundary layer meteorology features of the Pearl River Delta, a boundary layer concept model of the Pearl River Delta was built. The concept model consists of four fundamental factors that affect the boundary layer meteorology of the Pearl River Delta and can convincingly explain the reason of the air quality change in the Pearl River Delta. The model can be used to the diffusion capability analysis, the air pollution potential forecasting or haze forecasting, etc.

Applicability of Temperature Discrete Equation to NMRF Boundary Layer Scheme in GRAPES Model

By deriving the discrete equation of the parameterized equation for the New Medium-Range Forecast(NMRF)boundary layer scheme in the GRAPES model,the adjusted discrete equation for temperature is obviously different from the original equation under the background of hydrostatic equilibrium and adiabatic hypothesis.In the present research,three discrete equations for temperature in the NMRF boundary layer scheme are applied,namely the original(hereafter NMRF),the adjustment(hereafter NMRF-gocp),and the one in the YSU boundary-layer scheme(hereafter NMRF-TZ).The results show that the deviations of height,temperature,U and V wind in the boundary layer in the NMRF-gocp and NMRF-TZ experiments are smaller than those in the NMRF experiment and the deviations in the NMRF-gocp experiment are the smallest.The deviations of humidity are complex for the different forecasting lead time in the three experiments.Moreover,there are obvious diurnal variations of deviations from these variables,where the diurnal variations of deviations from height and temperature are similar and those from U and V wind are also similar.However,the diurnal variation of humidity is relatively complicated.The root means square errors of 2m temperature(T2m)and 10m speed(V10m)from the three experiments show that the error of NMRF-gocp is the smallest and that of NMRF is the biggest.There is also a diurnal variation of T2m and V10m,where T2m has double peaks and V10m has only one peak.Comparison of the discrete equations between NMRF and NMRF-gocp experiments shows that the deviation of temperature is likely to be caused by the calculation of vertical eddy diffusive coefficients of heating,which also leads to the deviations of other elements.

Turbulence Model Investigations on the Boundary Layer Flow with Adverse Pressure Gradients

In this paper, a numerical study of flow in the turbulence boundary layer with adverse and pressure gradients （APGs） is conducted by using Reynolds-averaged Navier-Stokes （RANS） equations. This research chooses six typical turbulence models, which are critical to the computing precision, and to evaluating the issue of APGs. Local frictional resistance coefficient is compared between numerical and experimental results. The same comparisons of dimensionless averaged velocity profiles are also performed. It is found that results generated by Wilcox （2006） k-co are most close to the experimental data. Meanwhile, turbulent quantities such as turbulent kinetic energy and Reynolds-stress are also studied.

LAGRANGIAN MODEL ON THE TURBULENT MOTION OF SMALL SOLID PARTICLE IN TURBULENT BOUNDARY LAYER FLOWS

The Lagrangian equations of motion of small particle in turbulent boundary layer flows, taking into account the effects of the drug force caused by the wall presence, the Saffman and the Magus lift forces et al., is studied. Using the spectral method, analytical expression relating to the Lagrangian power spectra of particle velocity to that of the fluid are developed and the results are used to evaluate various responses statistics. In this paper, the results clearly show that the turbulent diffusivity of the particle may be larger than that of the fluid for a period of long-time.

Wind Structure in an Intermediate Boundary Layer Model Based on Ekman Momentum Approximation

A quasi three–dimensional, intermediate planetary boundary layer (PBL) model is developed by including inertial acceleration with the Ekman momentum approximation, but a nonlinear eddy viscosity based on Blackadar’s scheme was included to improve the theoretical model proposed by Tan and Wu (1993). The model could keep the same complexity as the classical Ekman model in numerical, but extends the conventional Ekman model to include the horizontal accelerated flow with the Ekman momentum approximation. A comparison between this modified Ekman model and other simplified accelerating PBL models is made. Results show that the Ekman model overestimates (underestimates) the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow due to the neglect of the acceleration flow, however, the semi–geostrophic Ekman model overestimates the acceleration effects resulting from the underestimating (overestimating) of the wind speed and pumping velocity in the cyclonic (anticyclonic) shear flow. The Ekman momentum approximation boundary layer model could be applied to the baroclinic atmosphere. The baroclinic Ekman momentum approximation boundary layer solution has both features of classical baroclinic Ekman layer and the Ekman momentum approximate boundary lager.

PRELIMINARY STUDY ON THE FORMATION MECHANISM OF COUNTER WESTERN BOUNDARY UNDERCURRENTS BELOW THE THERMOCLINE——A CONCEPTUAL MODEL

Based on a simple conceptual model of stratified ocean, the criterion of the geostrophic velocity inversion in and below the thermocline was derived as h ′·η′<0 and ρ 1|η′|≤Δρ| h ′|, meaning that the slopes of the thermocline( h ′) and the sea surface(η′) must be opposite to each other, and that h ′ must be strong enough to satisfy the latter inequality. The criterion was applied to discuss the features of the western boundary undercurrents, the counter undercurrents of the western boundary currents below the thermocline, and to discuss the dynamics of their formation finally resulting from the combination of the basin scale circulation and local geostrophic balance. The formation mechanism, multi core structure, and transport variations of the Mindanao Undercurrent and those of other undercurrents, such as the North Equatorial Undercurrent and the Kuroshio undercurrent, can be satisfactorily explained by the above results.

Large Eddy Simulation and Study of the Urban Boundary Layer

Based on a pseudo-spectral large eddy simulation (LES) model, an LES model with an anisotropy turbulent kinetic energy (TKE) closure model and an explicit multi-stage third-order Runge-Kutta scheme is established. The modeling and analysis show that the LES model can simulate the planetary boundary layer (PBL) with a uniform underlying surface under various stratifications very well. Then, similar to the description of a forest canopy, the drag term on momentum and the production term of TKE by subgrid city buildings are introduced into the LES equations to account for the area-averaged effect of the subgrid urban canopy elements and to simulate the meteorological fields of the urban boundary layer (UBL). Numerical experiments and comparison analysis show that: (1) the result from the LES of the UBL with a proposed formula for the drag coefficient is consistent and comparable with that from wind tunnel experiments and an urban subdomain scale model; (2) due to the effect of urban buildings, the wind velocity near the canopy is decreased, turbulence is intensified, TKE, variance, and momentum flux are increased, the momentum and heat flux at the top of the PBL are increased, and the development of the PBL is quickened; (3) the height of the roughness sublayer (RS) of the actual city buildings is the maximum building height (1.5-3 times the mean building height), and a constant flux layer (CFL) exists in the lower part of the UBL.

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers

Kinematic models compute the temperature distribution by prescribing a constant convergent velocity for the subducting slab,resulting in an artificial velocity discontinuity,which may accelerate the heating of subducting slabs.For the purpose to moderate the influence of such artificial discontinuity, we construct a 2D thermal model for subduction zones with a velocity boundary layer,within which the velocities decrease linearly with the distance from the interfaces of slabs.Temperatures are calculated

LONGITUDINAL VORTEX PAIR IMBEDDED IN A TURBULENT BOUNDARY LAYER:COMPARATIVE PERFORMANCES OF SEVERAL TURBULENCE MODELS

Four turbulence models, namely, the basic and nonlinear stress-transport models and the basic and anisotropic k- epsilon models, have been tested in the case of interaction between a longitudinal vortex pair and a flat-plate boundary layer. The results of their predictions were compared with Mehta and Bradshaw's measurements. Part of the results involving those of the nonlinear stress-transport model and anisotropic k- epsilon model are presented and discussed. (Edited author abstract) 13 Refs.

Effect of absorption boundary layer on nonlinear flow in low permeability porous media

Taking low permeability cores of Daqing oilfield for example,the flow characteristics at low velocity were studied with the self-designed micro-flux measuring instrument.Considering the throat distribution and capillary model,the thickness of fluid boundary layer under different pressure gradients was calculated,and the mechanism and influencing factors of nonlinear percolation were discussed.The results show that the percolation curve of ultra-low rocks is nonlinear,and apparent permeability is not a constant which increases with pressure gradient.The absorption boundary layer decreases with the increase of pressure gradient,and changes significantly especially in low pressure gradient,which is the essence of nonlinear percolation.The absorption boundary layer is also found to be impacted by the surface property of rocks.

Footprint Characteristics of Scalar Concentration in the Convective Boundary Layer

Footprint characteristics for passive scalar concentration in the convective boundary layer （CBL） are investigated. A backward Lagrangian stochastic （LS） dispersion model and a large eddy simulation （LES） model are used in the investigation. Typical characteristics of the CBL and their responses to the surface heterogeneity are resolved from the LES. Then the turbulence fields are used to drive the backward LS dispersion. To remedy the spoiled description of the turbulence near the surface, MoninObukhov similarity is applied to the lowest LES level and the surface for the modeling of the backward LS dispersion. Simulation results show that the footprint within approximately 1 km upwind predominates in the total contribution. But influence from farther distances also exists and is even slightly greater than that from closer locations. Surface heterogeneity may change the footprint pattern to a certain degree. A comparison to three analytical models provides a validation of the footprint simulations, which shows the possible influence of along-wind turbulence and the large eddies in the CBL, as well as the surface heterogeneity.

Implementation of Turbulent Mixing over a Stratocumulus-Topped Boundary Layer and Its Impact in a GCM

The effect of a vertical diffusion scheme over a stratocumulus topped boundary layer （STBL） was investigated using the YONU AGCM （Yonsei University Atmospheric General Circulation Model）. To consider the impact of clouds on the turbulence production, the turbulence mixing term, driven by radiative cooling at the cloud top, is implemented as an extended non-local diffusion scheme. In the model with this new scheme, the STBL parameterization significantly influences the lower atmosphere over the tropical and subtropical regions. Consideration of the turbulent mixing within the cloud layer leads to continuous stratocumulus formation. The cloud-top radiative cooling tends to favor more rapid entrainment and produces top-down turbulent mixing. This cooling develops a mixed layer without initiation of deep convection by surface fluxes. Variations in thermodynamical and dynamical features are produced by planetary boundary layer （PBL） cloud development. The simulated stratocumulus induces more mixing of heat and moisture due to the cloud forcing. Over STBL regions, the lower boundary layer becomes warmer and drier. It also weakens vertical motion and zonal trade winds in the eastern Pacific, which indicates that stratocumulus cloud cover plays a role in weakening the Walker circulation; that is, cloud cover damps the tropical circulation.

A numerical study for boundary layer current and sheet flow transport induced by a skewed asymmetric wave

An analytical model with essential parameters given by a two-phase numerical model is utilized to study the net boundary layer current and sediment transport under skewed asymmetric oscillatory sheet flows. The analytical model is the first instantaneous type model that can consider phase-lag and asymmetric boundary layer development. The two-phase model supplies the essential phase-lead, instantaneous erosion depth and boundary layer development for the analytical model to enhance the understanding of velocity skewness and acceleration skewness in sediment flux and transport rate. The sediment transport difference between onshore and offshore stages caused by velocity skewness or acceleration skewness is shown to illustrate the determination of net sediment transport by the analytical model. In previous studies about sediment transport in skewed asymmetric sheet flows, the generation of net sediment transport is mainly concluded to the phase-lag effect.However, the phase-lag effect is shown important but not enough for the net sediment transport, while the skewed asymmetric boundary layer development generated net boundary layer current and mobile bed effect are key important in the transport process.

Mesoscale Vegetation-Breeze Circulations and Their Impact on Boundary Layer Structures at Night

The impact of well watered mesoscale wheat over mid-latitude arid areas on mesoscale boundary layer structures (MBLS) and climate has been investigated in the study .using a mesoscale biophysical, meteorological model (BM) developed in the current study. The BM is composed of six modules:mesoscale atmospheric module, soil module, vegetation module, snow-atmosphere interaction module, underlying surface meteorology module and subgrid scale flux parameterization module. The six modules constitute an interacting system by supplying boundary conditions to each other.The investigation indicates that a horizontal pressure gradient associated with mesoscale perturbations in temperature and humidity is created during the day, which results from more water transpired from the vegetation canopy (VC) and evaporated from underlying wet soil. Non-classical mesoscale circulations (called as vegetation-breeze) are forced by the pressure perturbations with wind speeds about 5 m / s, flowing from the VC to the adjacent bare soil in the low boundary layer.Specific attention in the study is given to study the MBLS at night caused by the interactions between these mesoscale vegetation-breezes.The impacts of the VC's scale and atmospheric background thermal stability (ABTS) on these thermally forced mesoscale vegetation-breezes are also investigated.

Numerical Investigation on Two-dimensional Boundary Layer Flow with Transition

As a basic problem in many engineering applications, transition from laminar to turbulence still remains a difficult problem in computational fluid dynamics （CFD）. A numerical study of one transitional flow in two-dimensional is conducted by Reynolds averaged numerical simulation （RANS） in this paper. Turbulence model plays a significant role in the complex flows＇ simulation, and four advanced turbulence models are evaluated. Numerical solution of frictional resistance coefficient is compared with the measured one in the transitional zone, which indicates that Wilcox （2006） k-ω model with correction is the best candidate. Comparisons of numerical and analytical solutions for dimensionless velocity show that averaged streamwise dimensionless velocity profiles correct the shape rapidly in transitional region. Furthermore, turbulence quantities such as turbulence kinetic energy, eddy viscosity, and Reynolds stress are also studied, which are helpful to learn the transition＇s behavior.

Numerical Experiments of Three Boundary Layer Parameterization Schemes on Simulation Results of Low-vortex Rainstorm Weather in Langfang Region

Using mesoscale numerical model WRFV3.7,simulation tests of a low-vortex rainstorm were carried out in Langfang area,and simulation results of three boundary layer parameterization schemes(MYJ,YSU,ACM2)were contrasted.The results showed that the rainstorm in Langfang area had better dynamic condition.By the influence of northeast cold vortex,cold vortex at high altitude cooperated with the surface wind speed convergence,which provided better dynamic condition for strong convective rainfall.Three boundary layer parameterization schemes all simulated surface wind speed convergence of rainstorm center.Simulation results of 24-h rainfall distribution showed that the simulation of YSU scheme was better than ACM2 and MYJ.The simulation results of flow field and temperature field also verified that YSU scheme was better than other schemes.

A Numerical Research on a Compressibility-correlated Langtry＇s Transition Model for Double Wedge Boundary Layer Flows

A new compressibility correlation is introduced in the Langtry＇s local variable-based transition model to investigate the phe- nomenon on double wedge shock/boundary layer interactions. The cmnputational analysis compared with experimental data has been made to assess the influence of the wall temperature and the leading edge nose radius on a hypersonic double wedge boundary layer. It has been found that the laminar boundary layer separation occurs on the first ramp. Furthermore, the wall temperature and the leading edge nose radius have remarkable influence on the separation characteristics in the kink. Comparison of the calculated pressure coefficient distribution and the boundary layer profile with the experimental data shows that better results can be achieved when using the modified transition model.