首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The structure of supercritical western boundary currents is investigated using a quasi-geostrophic numerical model. The basic flow is of meridional Munk balance, and the input boundary is perturbed by the most unstable wave solution obtained from linear spatial instability calculations. Self-preserving (or equilibrium) solutions are obtained for the model runs at Re=30, 60, 90, and 120, and their energy and vorticity budgets are analyzed. In an analogy with the laboratory turbulence of wall boundary layers, the western boundary layer is divided into inner and outer layers. In the inner layer, the mean energy is dissipated via direct viscous dissipation, while in the outer layer it is converted to the eddy energy via turbulence production. The main scenario is that the mean energy is produced in the inner layer via ageostrophic pressure work divergence, and it is partly removed due to viscous action within a narrow region near the wall, defined here as viscous sub-layer. The remaining portion is converted to the eddy energy via turbulence production in the outer layer, which is in turn transported to the inner layer, then again to the viscous sub-layer where it is ultimately dissipated. In the near-wall side, the vorticity balance of the mean flow is maintained by viscous effect and Reynolds flux divergence, while in the offshore side it is maintained by beta effect and Reynolds flux divergence. The length scale of the supercritical boundary current is roughly , where LM is the Munk length, as observed from a dimensional analysis.  相似文献   

2.
Two-dimensional asymmetric merger of two like-signed vorticity monopoles with different sizes and vorticities is examined by combining simplified analytical models and contour dynamics experiments. The model results can capture the key dynamics and hence allow the prediction of the critical merger distance in a number of the situations. The models ignore deformation of one of the two vortices, replacing it with a point vortex, and employ a corotating frame of reference with a rotation rate estimated by point vortices. Thus, the two vortex problem becomes two separate problems of a single vortex in a background shear flow. Vortex merger is found to happen when the vortex cannot resist the background shear flow. Vortex merger and merging processes depend on the centroid distance d, the circulation ratio, (qi and ri are the vorticity and radius, respectively) and initial conditions. In the lowest order, the background flow is approximated by a uniform shear field, and the behavior of an elliptical vortex can be described by the Kida (1981) equation supplemented with one describing the time evolution of the centroid distance. This model reveals that merger takes place because the natural rotation of an elliptical vortex is overcome by the background uniform shear flow; the ellipse inversely rotates and is drawn out by the background straining field. The vortex deformation in a background flow field induces an inward flow at the position of the other vortex; as a result, the centroid distance decreases and two vortices merge. The critical merger distance from this model agrees quite well with the results from contour dynamics experiments for two vortices. Inclusion of higher order non-uniform shear in the background flow extends the critical merger distance, which gives almost perfect estimates for the experiment. In the non-uniform shear flow, partial merger occurs, where the vortex sheds off a filament, but the remaining part of the vortex resumes its natural rotation.  相似文献   

3.
  相似文献   

4.
A Field Study of the Mean Pressure About a Windbreak   总被引：3，自引：0，他引：3  
To provide additional field data for assessingwindbreak flow models, mean ground-level pressurehas been measured upstream and downstream from along porous fence (height H = 1.25 m, resistancecoefficient k r = 2.4). Measurements were madeduring periods of near-neutral stability and near-normallyincident flow, with the fence standing on bare soil(roughness length, z 0 0.8 cm;H/z 0 160), or within a plant canopy. The mean pressure field,measured far from the ends of the fence, was foundto be quite insensitive to mean wind direction( , zero for perpendicular flow), for| | less than about 25°.In the absence of a canopy, during each measurementperiod the minimum pressure occurred at the closestsampling location to leeward of the windbreak, thepressure-gradient in most cases beingmaximally-adverse in the immediate lee, and decayingwith increasing downwind distance (x). On one day ofmeasurements, however, the pressure gradient over2 x/H 6 (H = windbreak height) resembled theleeward plateau identified by Wang and Taklein their numerical studies. Perhaps thisoccasional feature was only due to instrumenterror. Nevertheless a plateau of sorts wasindicated in similar measurements by Judd andPrendergast (with H = 1.92 m, z 0 1.2 cm;H/z 0 160, k r 3). Therefore,existence of a leeward pressure plateau behind athin fence cannot be definitely ruled out.When the windbreak was placed in a canopy, minimumsurface pressure was displaced downwind. Thisagrees with the wind-tunnel study of Judd, Raupach and Finnigan,and is consistent with a simple simulation reported here.  相似文献   

5.
This study attempts to determine the scales of turbulence in a high Reynolds number shear flow near which transition to isotropy occurs and the scales for which Taylor's hypothesis is applicable. The flow studied was the wind near height x 3 = 2 m above a flat land surface. Four hot-wire anemometers were mounted in a three-dimensional array with equal separations between 1.8 m and 2 cm in three different directions. Theoretical cross-spectra were computed from the observed spectra of downwind velocity fluctuations assuming isotropy and Taylor's hypothesis. Comparison between these and the observed cross-spectra revealed that the turbulence in the flow studied was consistent with both assumptions provided k 1x3&> 20, where k 1 is the radian wavenumber; this was the lower bound to which no departure from isotropy could be detected by the experiment. For 4 k 1x3 20, the observations are consistent with symmetry of the turbulence about the downstream direction. That part of Taylor's hypothesis relating observed frequency at a stationary sensor to the downstream wavenumber component appears to be justified within experimental error for k 1x3& > 3.  相似文献   

6.
This article outlines a critical gap in the assessment methodology used to estimate the macroeconomic costs and benefits of climate and energy policy, which could lead to misleading information being used for policy-making. We show that the Computable General Equilibrium (CGE) models that are typically used for assessing climate policy use assumptions about the financial system that sit at odds with the observed reality. These assumptions lead to ‘crowding out’ of capital and, because of the way the models are constructed, negative economic impacts (in terms of gross domestic product (GDP) and welfare) from climate policy in virtually all cases.

In contrast, macro-econometric models, which follow non-equilibrium economic theory and adopt a more empirical approach, apply a treatment of the financial system that is more consistent with reality. Although these models also have major limitations, they show that green investment need not crowd out investment in other parts of the economy – and may therefore offer an economic stimulus. Our conclusion is that improvements in both modelling approaches should be sought with some urgency – both to provide a better assessment of potential climate and energy policy and to improve understanding of the dynamics of the global financial system more generally.

POLICY RELEVANCE

This article discusses the treatment of the financial system in the macroeconomic models that are used in assessments of climate and energy policy. It shows major limitations in approach that could result in misleading information being provided to policy-makers.  相似文献   


7.
The note presents a rational approach to modelling the source/sink due to vegetation or buoyancy effects that appear in the turbulent kinetic energy, E, equation and a supplementary equation for a length-scale determining variable, φ, when two-equation closure is applied to canopy and atmospheric boundary-layer flows. The approach implements only standard model coefficients C φ1 and C φ2 in the production and destruction terms of the φ equation, respectively. Numerical tests illustrate the practical applicability of the method, where, for example, simulations with the E–ω model (where is the specific dissipation and is the dissipation rate of E) properly reproduce both the surface-layer wind profile estimated from the Monin-Obukhov similarity theory and the mixing-height evolution observed above forested terrain in Southern Finland.  相似文献   

8.
The relation between the turbulence Reynolds numberR and a Reynolds numberz* based on the friction velocity and height from the ground is established using direct measurements of the r.m.s. longitudinal velocity and turbulent energy dissipation in the atmospheric surface layer. Measurements of the relative magnitude of components of the turbulent kinetic energy budget in the stability range 0 >z/L 0.4 indicate that local balance between production and dissipation is maintained. Approximate expressions, in terms of readily measured micrometeorological quantities, are proposed for the Taylor microscale and the Kolmogorov length scale .  相似文献   

9.
Summary A zonally averaged global energy balance model with feedback mechanisms was constructed to simulate (i) the poleward limits of ITCZ over the continent and over the ocean and (ii) a simple monsoon system as a result of differential heating between the continent and the ocean. Three numerical experiments were performed with lower boundary as (1) global continent, (2) global ocean and (3) continent-ocean, with freezing latitudes near the poles. Over the continent, midlatitude deserts were found and the ITCZ migrates 25° north and south with seasons. Over a global swamp ocean results do not show migration of ITCZ with time but once the ocean currents are introduced the ITCZ migrates 5° north and south with seasons. It was found that the seasonal migration of ITCZ strongly depends on the meridional distribution of the surface temperature. It was also found that continent influences the location of the oceanic ITCZ. In the tropics northward progression of quasi-periodic oscillations called events are found during the pre- and post-monsoon periods with a period of 8 to 15 days. This result is consistent with the observed quasi-periodic oscillations in the tropical region. Northward propagation of the surface temperature perturbation appears to cause changes in the sensible heat flux which in turn causes perturbations in vertical velocity and latent heat flux fields.List of Symbols vertical average - 0 zonal average - vertical mean of the zonal average - 0s zonal average at the surface - 0a zonal average at 500 mb level - latitude We now define the various symbols used in the model rate of atmospheric heating due to convective cloud formation (K/sec) - dp/dt (N/m2/sec) - density - potential temperature (K) - rate of rotation of the earth (rad/sec) - empirical constant - humidity mixing ratio - * saturated humidity mixing ratio - opacity of the atmosphere - 1,2 factors for downward and upward effective black body long wave radiation from the atmosphere - Stefan-Boltzmann constant - emissivity of the surface - D subsurface temperature (K) - a specific volume - 0xs ,0ys eastward and northward components of surface frictional stress - * vertical velocity at the top of the boundary layer (N/m2/sec) - P Thickness of the boundary layer (mb) - nondimensional function of pressure - P pressure - P a pressure of the model atmosphere (N/m2) - P s pressure at the surface (N/m2) - t time (sec) - U eastward wind speed (m/sec) - V northward wind speed (m/sec) - surface water availability - T absolute temperature (K) - heat addition due to water phase changes - g acceleration due to gravity (m2/sec) - a radius of the earth (m) - R gas constant for dry air (J/Kg/K) - C p specific heat of air at constant pressure (J/Kg/K) - k R/C p - L latent heat of condensation (J/Kg) - f coriolis parameter (rad/sec) - H s H 0s (1) +H 0s (2) +H 0s (3) +H 0s (4) +H 0s (5) (J/m2/Sec)=sum of the rates of vertical heat fluxes per unit surface area, directed toward the surface - H a H 0a (1) +H 0a (2) +H 0a (3) +H 0a (4) (J/m2/Sec)=sum of the rates of heat additions to the atmospheric column per unit horizontal area by all processes - H 0s (1) ,H 0a (1) heat flux due to short wave radiation - H 0s (2) ,H 0a (2) heat flux due to long wave radiation - H 0s (3) ,H 0a (3) heat flux due to small scale convection - H 0s (4) heat flux due to evaporation - H 0a (4) heat flux due to condensation - H 0s (5) heat flux due to subsurface conduction and convection - e * saturation vapor pressure - R solar constant (W/m2) - r a albedo of the atmosphere - r s albedo of the surface - b 2 empirical constant (J/m2/sec) - c 2 empirical constant (J/m2/sec) - e 2 nondimensional empirical constant - f 2 empirical constant (J/m2/sec) - factor proportional to the conductive capacity of the surface medium - a s constant used in Sellers model - b s positive constant of proportionality used in the Sellers model (kg m2/J/sec2) - K HT coefficient for eddy diffusivity of heat (m2/sec) - K HE exchange coefficient for water vapor (m2/sec) - h depth of the water column (m) - z height (m) - V 0ws meridional component of surface current (m/sec) - n cloud amount - G 0,n long wave radiation form the atmosphere for cloud amount n (W/m2) - B 0 long wave radiation from the surface (W/m2) - S 0,n short wave radiation from the atmosphere for cloud amount n (W/m2) - A n albedo factor for a cloud amount n - R f1 large scale rainfall (mm/day) - R f2 small scale rainfall (mm/day) With 22 Figures  相似文献   

10.
In a recent paper, the author introduced a new viscous boundary layer, called the mesolayer, in turbulent shear flow. Its importance stems from its location between the inner and outer regions which are controlled by the law of the wall and Reynolds number similarity, respectively. This intrusion prevents the classical overlap assumption which appears to be fundamental in the derivation of the classical logarithmic behavior. The mesolayer has a thickness proportional to Taylor's microscale . This, and the analogy between the energy equation for the spectrum function of isotropic turbulence and the momentum equation for shear flow, suggest the existence of a similar region in wavenumber space with wavenumber k ~ -1. This mesoregion separates the inner region k ~ k s(where k s-1 and is the Kolmogorov length) and the outer region k k e(where k e -1 and l is the energy-containing eddy size) and again invalidates the overlap assumption which appears to be fundamental in the derivation of the classical k -5/3-behavior of the inertial subrange.Incorporation of the mesoregion into the argument leads to a new theory with k -5/3-behavior in two regions (-1 k k s) and (k e k -1) although with two different coefficients of proportionality (Kolmogorov constants). This leads to a wandering of the spectrum curve about the classical k -5/3 line similar to a wandering in turbulent shear flow about the logarithmic curve. This is clearly indicated by the data for the variation of the Kolmogorov constant.Other data support the new theory. In particular, the location of the point k mwhere the curve of the nonlinear energy-transfer function goes through zero shows agreement with the theory, i.e., k m-1.  相似文献   

11.
It is well known that if turbulent mass convection is modelled as diffusion, errors result unless trajectories from the source (ath) to the point of observation (z p ) comprise many statistically-independent segments (Taylor, 1921). We show that this is not guaranteed merely by the Lagrangian timescale () at the source being small (e.g., source at ground), but that a better criterion istmax[(h), (z p )], wheret is a typical travel time toz p .  相似文献   

12.
Water vapour is the largest and, radiatively, the most important trace gas in the Earth's atmosphere. Cloud amount and cloud optical depth feedbacks are, as yet, poorly understood and improvements in model parameterization schemes await an adequate observational data base. Satellite retrievals do not, and will not for some time, provide more than snapshot (a few months to a few years) records. Conventional surface-based observations of cloudiness could complement the global coverage offered by satellite retrievals if a sound relationship between the two observational measurements could be found. Observations underline the importance of the vertical dimension of clouds which affects the observational geometries from satellites and the surface. A new basis for the relationship between the (vertical) earthview of cloud amount and the (whole dome) skycover of cloud amount has been sought. Over four and a half thousand all-sky camera photographs, representing a considerable range of seasonal and climatological conditions, have been analyzed to give rise to a database from which predictive relationships for earthview, E, and skycover, S, have been established.Cubic functions are the most soundly based both physically and empirically. We find: S=0.647+2.192E–0.461E 2+0.037E 3 and E=-0.001+0.082S+0.269S 2–0.019S 3 for the prediction of skycover from earthview and earthview from skycover, respectively. If earthview is required from skycover observations then ES could be used with little additional error. Hence, conventional surface observations of skycover could be compared directly with satellite-derived earthview values. More importantly, these results do not support the widespread assumption that conventional (surface) observations of cloud amount always exceed the earthview could retrieval. Furthermore, climate model predictions of total cloud amount may also be interpreted via these relationships. GCM-predicted layer cloud amounts can be synthesized into modelled E values using the random overlap formula and hence it is possible to construct modelled S values which are directly comparable with conventionally observed cloud climatologies. The baseline observation of skycover by clouds therefore provides a valuable validation tool for both satellite programmes and climate models.  相似文献   

13.
On mountain wave drag over complex terrain   总被引：1，自引：0，他引：1  
Summary Mountain wave drag is calculated for rotating, stratified, nonhydrostatic Boussinesq flow over a mountain ridge using linear theory for a variety of mountain profiles representing complex/irregular terrain. The inclusion of a sinusoidal corrugation to the familiar witch-of-Agnesi profile creates a stegosaurus profile. The associated drag is greatly enhanced for mesoscale mountains when the corrugation wave-number matches that for the dominant inertia-gravity wave contribution to the cross-mountain surface pressure gradient. Similarly, increasing the jaggedness (by decreasing the exponentb) increases the drag for mesoscale mountains whose topographic spectral intensity,M(k), has the form of a power law:M(k)=mk –b wherek is the zonal wavenumber.Spectral analysis of one-kilometer resolution topographic data for the Appalachian Mountains suggests that a power law profile withb=1.7 accurately represents the topographic spectral intensity and that it yields good estimates of the drag.The application of these results to the parameterization of mountain wave drag in general circulation models is discussed.With 7 Figures  相似文献   

14.
A modification of the most popular two-equation (E–φ) models, taking into account the plant drag, is proposed. Here E is the turbulent kinetic energy (TKE) and φ is any of the following variables: El (product of E and the mixing length l), (dissipation rate of TKE), and ω (specific dissipation of TKE, ). The proposed modification is due to the fact that the model constants estimated experimentally for ‘free-air’ flow do not allow for adequate reconstruction of the ratio between the production and dissipation rates of TKE in the vegetation canopy and have to be adjusted. The modification is universal, i.e. of the same type for all E–φ models considered. The numerical experiments carried out for both homogeneous and heterogeneous plant canopies with E–φ models (and with the E–l model taken as a kind of reference) show that the modification performs well. They also suggest that E– and E–ω schemes are more promising than the E–El scheme for canopy flow simulation since they are not limited by the need to use a wall function.In addition, a new parameterization for enhanced dissipation within the plant canopy is derived. It minimizes the model sensitivity to C μ, the key parameter for two-equation schemes, and whose estimates unfortunately vary considerably from experiment to experiment. The comparison of results of new modified E– and E –ω models with observations from both field and wind-tunnel experiments shows that the proposed parameterization is quite robust. However, because of uncertainties with the turbulence Prandtl and Schmidt numbers for the E– model within the canopy, the E–ω model is recommended for future implementation, with the suggested modifications.  相似文献   

15.
Numerical solutions are obtained for a two-dimensional eddy diffusion equation describing aerosol particles in the atmospheric surface layer. Particles are assumed to be from either an elevated line source or a surface area source. Characteristics of the solutions are exhibited as a function of source type, lower boundary condition, and the ratio of terminal velocity () to friction velocity (u *). This theory indicates that suspended particles tend to escape the near-surface layer, and are thus subject to long-range transport, when the ratio /u * is less than von Karman's constant. This prediction is consistent with available data.  相似文献   

16.
A comprehensive model for the prediction of concentration fluctuations in plumes dispersing in the complex and highly disturbed wind flows in an urban environment is formulated. The mean flow and turbulence fields in the urban area are obtained using a Reynolds-averaged Navier-Stokes (RANS) flow model, while the standard k-ϵ turbulence model (k is the turbulence kinetic energy and ϵ is the viscous dissipation rate) is used to close the model. The RANS model provides a specification of the velocity statistics of the highly disturbed wind flow in the urban area, required for the solution of the transport equations for the mean concentration and concentration variance (both of which are formulated in the Eulerian framework). A physically-based formulation for the scalar dissipation time scale t d , required for the closure of the transport equation for , is presented. This formulation relates t d to an inner time scale corresponding to “internal” concentration fluctuation associated with relative dispersion, rather than an outer time scale associated with the entire portion of the fluctuation spectrum. The two lowest-order moments of concentration ( and ) are used to determine the parameters of a pre-chosen functional form for the concentration probability density function (clipped-gamma distribution). Results of detailed comparisons between a water-channel experiment of flow and dispersion in an idealized obstacle array and the model predictions for mean flow, turbulence kinetic energy, mean concentration, concentration variance, and concentration probability density function are presented.  相似文献   

17.
We propose a simple model for esimating the average number of occurrences per unit time (c o) that a threshold concentration c o is exceeded. It is based on the joint probability density of the observed concentration c(t) and its time derivative (t) under the assumption that c(t) is a stationary time series; this assumption leads to the hypothesis that c(t) and (t) are statistically independent. Adopting plausible forms of the frequency distributions of c and , we apply the model to diffusion from an infinite area source and from an elevated point source, both in the neutral boundary layer, and obtain simple results for (c o) and the average duration of one excursion above c o as functions of c o, the mean and the standard deviation of the concentration, and surface-layer variables.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

18.
The energy balance was measured for the dry canopy of narrow-leaved snow tussock (Chionochloa rigida), and measurements of transpiration were obtained from a large weighing lysimeter.Typical maximum summer transpiration rates of 0.21–0.43 mmhr-1 (140–290 W m-2) were recorded. The latent heat flux accounted for less than 40% of net radiation. The estimated value of the bulk stomatal resistance (r ST) for 29 days was 158 s m-1, and the decoupling parameter () was 0.17. Transpiration rates were not driven directly by net radiation, but were closely linked to the size of the regional saturation deficit imposed at the level of the canopy by efficient overhead mixing, and were constrained by a large bulk stomatal resistance. A linear relationship between r ST and the saturation deficit is proposed as a realistic method for estimating transpiration for water yield studies of tussock catchments.  相似文献   

19.
An easily-evaluated expression for the dimensionless concentration profile (z/z 0,/z 0, z 0/L) = = cu */kQ (or z 0cu*/kQ) downwind of a continuous ground-level area (or line) source in the stable surface layer is obtained by integrating the diffusion equation using the Shwetz approximation method (c = concentration, Q = source strength, k = von Kárman's constant). The analytical solution compares closely with concentration profiles obtained using a trajectory-simulation model over a useful range of heights, the important discrepancies occurring at the upper edge of the plume. The analytical solution is used to generate predictions of ground-level concentration for the Project Prairie Grass experiments; good agreement with the observations is obtained at all downwind distances (50 to 800 m).  相似文献   

20.
The turbulent flow inside dense canopies is characterized by wake production and short-circuiting of the energy cascade. How these processes affect passive scalar concentration variability in general and their spectral properties in particular remains a vexing problem. Progress on this problem is frustrated by the shortage of high resolution spatial concentration measurements, and by the lack of simplified analytical models that connect spectral modulations in the turbulent kinetic energy (TKE) cascade to scalar spectra. Here, we report the first planar two-dimensional scalar concentration spectra (ϕ cc ) inside tall canopies derived from flow visualization experiments. These experiments were conducted within the deeper layers of a model canopy composed of densely arrayed cylinders welded to the bottom of a large recirculating water channel. We found that in the spectral region experiencing wake production, the ϕ cc exhibits directional scaling power laws. In the longitudinal direction (x), or the direction experiencing the largest drag force, the ϕ cc (k x ) was steeper than and followed an approximate at wavenumbers larger than the injection scale of wake energy, where k x is the longitudinal wavenumber. In the lateral direction (y), the spectra scaled as up to the injection scale, and then decayed at an approximate power law. This departure from the classical inertial subrange scaling (i.e., k −5/3) was reproduced using a newly proposed analytical solution to a simplified scalar spectral budget equation. Near the velocity viscous dissipation range, the scalar spectra appear to approach an approximate k −3, a tantalizing result consistent with dimensional analysis used in the inertial-diffusive range. Implications to subgrid modelling for large-eddy simulations (LES) inside canopies are briefly discussed.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号