首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 36 毫秒
1.
The transition to “multi-α” solutions of the upper bounding problem for thermal convection is discussed. For convection in a fluid contained between parallel stress-free perfectly conducting boundaries, the “single-α” solutions of Straus (1973) are used to determine the Rayleigh number R at which the first transition occurs. Two upper bounding problems are treated: one valid for all values of the Prandtl number and one valid only for large Prandtl numbers. A significant difference between the two problems is noted. The former has a transition at R ～ 28200; the latter has no transition within the range of Rayleigh numbers treated here: R ? 2.3 · 105.  相似文献   

2.
Multifractal analysis can provide parameters associated with different scales of rainfall, which may be useful for setting up parsimonious downscaling models of rainfall, or for revealing climate-specific properties. Time series of rain rate with 1-min resolution collected from ten stations over a monsoon watershed in eastern China were used to study the multifractal properties. The power spectra estimated by fast Fourier transform (FFT) and discrete Haar wavelet transform (DWT) showed three scaling regimes: the sub-hourly scaling regime with β?≈?1.2, the scaling regime from 1 h to 1 day with β close to 0.6, and the low-frequency spectra plateau with β?≈?0.1. From the hyperbolic tails of exceeding probability distributions, the estimated values of parameter q c are in 2–2.5, which were consistent with the critical order of K(q) curves. The statistical moments display two main scaling regimes: the high-frequency regime from 3 min to 5 days and the scaling regime beyond 5 days. The scales of 5–10 days seem a transitional regime. The reason that the regimes, revealed by the power spectra, disagree with the statistical moments may be that both FFT and DWT power spectra have limited abilities of analyzing low-frequency scaling but are sensitive to the properties in high-frequency scales. The H values estimated for the regime of sub-hourly scales are larger than 0.4, and the values for the regime 1 h–1 day are close to 0.1. For the low-frequency scales beyond 1 day, negative H is obtained by DWT power spectra. The parameters of universal multifractal models were also estimated. The values of α for the scaling range of 1 min–5 days are 0.486?±?0.047, and for the low-frequency scaling range, its values are 0.808?±?0.323. For the high- and low-frequency scaling ranges, the values of C 1 are 0.5 and 0.169, respectively, which is different from the values for daily rainfall series collected at the same rain gages.  相似文献   

3.
The steady hydrostatic flow through a channel of rectangular cross section connecting reservoirs of infinite width and depth and containing inviscid fluids of different densities and levels is studied. The main goal is the determination of the discharges of the lighter and denser fluids in terms of the external conditions (reservoir levels, fluid densities and variation of width and depth along a channel). It is shown that the key parameter is δ, which is the ratio of relative reservoir level difference, γ, to relative density difference, ε. If δ<0 then the denser fluid plunges under the stationary lighter layer. If δ>δ* (1<δ*<1.5) then the lighter fluid runs up on a wedge of stationary heavier fluid. Here δ* depends on the geometry of the constriction. The solutions describing these regimes are stated. If 0<δ<δ* then both layers are in motion. A qualitative analysis of the solution for arbitrary bottom shape and channel width and arbitrary ε is presented and the problem is reduced to a system of two equations which can be easily solved numerically for any particular channel profile. We give detailed analyses for the following two cases: 1) the narrowest width of the channel is on the side of the heavier fluid and the top of the sill is on the side of lighter fluid; 2) the minima in channel depth and width coincide. In the second case the discharges for one class of geometries in the Boussinesq approximation are calculated and discussed.  相似文献   

4.
The dynamics of oceanic bottom currents are examined both theoretically and in the laboratory. A class of similarity solutions for steady flow indicates that the geostrophic current is drained by Ekman flux at its downslope edge and ultimately extinguished at a downstream distance of order (fQs2gr)12 magnified by E?12.Laboratory source flows are found to be consistently wave-like. Nevertheless, certain gross features of the steady Ekman flux mechanism are observed. Instabilities are classified according to the magnitudes of Rossby (?) and Ekman (E) numbers for the steady flow scaling. Characteristic were forms include: (1) a meandering jet (E < ?, 10?2); (2) transverse waves on broad contour current (? < E, 10?2); and (3) transverse waves on a viscous flow (?, E > 10?2). The dispersion relation for source flow waves resembles that of baroclinic instabilities for a uniform two-layer channel flow, and an empirically determined stability boundary is in rough agreement with the inviscid channel flow criterion.An interpretation of field measurements from the Denmark Strait Overflow in terms of the laboratory results is presented.  相似文献   

5.
The changes with timet of a temperature deviation δT(t,α) and of a vertical velocityW i(t,α) of an isolated dry thermal have been investigated theoretically. Solutions for the functionW i(t, α) have been derived for stable and unstable environmental stratifications. Comparing these solutions with the corresponding ones for the rise of an adiabatic thermal yield some interesting conclusions. Firstly, there is the evident relation between the rate of entrainment of environmental air (expressed by the parameter α=(1/M i) dM i/dz whereM i is the mass of the thermal) and the vertical velocity of the thermal: an increase in α decreases the velocity. Two similar thermals in stably stratified surroundings, one of them moving adiabatically (α=0) the other nonadiabatically (α>0), would rise for the same length of timet 2=π/N, whereN is a typical Brunt-Väisälä frequency, but with different velocities and to different heights: the ascent timet 2 depends only on environmental parameters. In an unstable stratification, the vertical non-adiabatic velocity of the thermal, instead of increasing without limit, tends towards a finite asymptotic velocity $$W_t (\infty ) = (\sqrt { - \mathcal{N}^2 } )/\alpha $$ the value of which depends upon both the stratification of the surroundings and upon the entrainment rate α. In a real atmosphere, where additional retarding forces exist, the motion will certainly be damped.  相似文献   

6.
Measurements of atmospheric turbulence made over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to determine the limits of applicability of Monin–Obukhov similarity theory (in the local scaling formulation) in the stable atmospheric boundary layer. Based on the spectral analysis of wind velocity and air temperature fluctuations, it is shown that, when both the gradient Richardson number, Ri, and the flux Richardson number, Rf, exceed a ‘critical value’ of about 0.20–0.25, the inertial subrange associated with the Richardson–Kolmogorov cascade dies out and vertical turbulent fluxes become small. Some small-scale turbulence survives even in this supercritical regime, but this is non-Kolmogorov turbulence, and it decays rapidly with further increasing stability. Similarity theory is based on the turbulent fluxes in the high-frequency part of the spectra that are associated with energy-containing/flux-carrying eddies. Spectral densities in this high-frequency band diminish as the Richardson–Kolmogorov energy cascade weakens; therefore, the applicability of local Monin–Obukhov similarity theory in stable conditions is limited by the inequalities Ri < Ri cr and Rf < Rf cr. However, it is found that Rf cr  =  0.20–0.25 is a primary threshold for applicability. Applying this prerequisite shows that the data follow classical Monin–Obukhov local z-less predictions after the irrelevant cases (turbulence without the Richardson–Kolmogorov cascade) have been filtered out.  相似文献   

7.
考虑了一类具p-Laplacian算子型带积分边界条件的三点奇异方程组边值问题正解的存在性.通过使用锥上的不动点定理,在适当的条件下,建立了这类方程组边值问题存在一个或多个正解的充分条件,并给出两个例子来验证主要结果.  相似文献   

8.
We investigate the momentum and energy exchange across the wave boundary layer (WBL). Directly at the air–sea interface, we test three wave-growth parametrizations by comparing estimates of the wave-induced momentum flux derived from wave spectra with direct covariance estimates of the momentum flux. An exponential decay is used to describe the vertical structure of the wave-induced momentum in the atmospheric WBL through use of a decay rate, a function of the dimensionless decay rate and wavenumber (A?=?α k). The decay rate is varied to minimize the difference between the energy extracted from the WBL and the energy flux computed from wave spectra using our preferred wave-growth parametrization. For wave ages (i.e. the peak phase speed to atmospheric friction velocity ratio) in the range \( 15 < c_{p}/u_{*} < 35 \) we are able to balance these two estimates to within 10%. The decay rate is used to approximate the WBL height as the height to which the wave-induced flux is 0.1 of its surface value and the WBL height determined this way is found to be between 1–3 m. Finally, we define an effective phase speed with which to parametrize the energy flux for comparison with earlier work, which we ultimately attempt to parametrize as a function of wind forcing.  相似文献   

9.
Flow in the urban boundary layer is strongly influenced by the surface roughness, which is composed principally of isolated buildings or groups of buildings. Previous research has shown that the flow regime depends on the characteristic height of these obstacles (H), and the spacing between them (W). In reality, the urban boundary layer contains roughness elements with a wide range of length scales; in many practical situations these can be classified into large-scale roughness—buildings, or groups of buildings—and small-scale roughness, such as street furniture and elements on the façades and roofs. It is important to understand how the small-scale roughness might modify mass and momentum transfer in the urban boundary layer, but relatively little information is available concerning the potential interaction between large- and small-scale roughness elements in the different flow regimes. This problem has been studied using wind-tunnel experiments, by measuring vertical velocity profiles over a two-dimensional obstacle array, adding small-scale roughness elements to the top of larger parallel square bars. The experiments were performed for different cavity aspect ratios: the results show that the small-scale roughness increases the turbulence intensities and the momentum transfer when the large-scale obstacles are closely packed (H/W > 1) but it has very little effect for more widely-spaced obstacles (H/W < 1).  相似文献   

10.
By use of the small parameter expansion method, the nonlinear planetary boundary layer (PBL) is studied in this paper. The PBL is divided into the surface layer and the Ekman layer, which is divided into several sublayers. In the surface-layer, the eddy coefficient K is taken as a linear function of height; in the Ekman layer, different constant K values are taken within different sublayers: these values are determined from O'Brien's formula (O'Brien, 1970) approximately. Under the upper and lower boundary conditions and the continuity conditions of the wind velocities and turbulent stresses at each boundary between sublayers, analytical expressions for wind velocity in all sublayers and the vertical velocity at the top of the PBL are obtained. A specific example of steady axisymmetrical circular high and low pressure areas is analysed, and some new conclusions are obtained. The results are in better agreement with reality than previous results. This example also shows that the vertical velocity at the top of the PBL caused by friction approaches zero near the center of a high or low pressure system for this model, but attains its maximum absolute values near the center of the high or low pressure area for Wu's (1984) model. This is due to the fact that in our model, the geostrophic wind speed near the center of this specific vortex approaches zero, which causes the wind shear and the friction effect to be very weak. Therefore the wind distribution in the PBL is very sensitive to the type of eddy coefficient.  相似文献   

11.
The stationary, Ekman-layer equations have been solved in closed form for two expressions of the eddy viscosity as a function of height, z: v τ=cu*z(1?z/h)and v τ=cu*z(1?z/h) 2, where u* is the friction velocity, h the boundary-layer height and c a constant. The main difference between both solutions is that the quadratic K-profile leads to a velocity discontinuity at the top of the boundary layer, while the solution for the cubic profile approaches the geostrophic wind at z=h smoothly. We discuss the characteristics of the solutions in terms of a dimensionless parameter C=fh/cu*, where f is the Coriolis parameter. The dependence on C can be interpreted in terms of a varying boundary-layer height or in terms of stability. The results for C ~ 1 are related to a neutral boundary layer. They agree well with results of a second-order model. The limit C → 0 is investigated in detail. We find that the stress profile becomes linear. The velocity profile shows different characteristics depending on whether we consider a shallow or a very unstable boundary layer. The results agree with observations. Finally we consider the influence of baroclinicity on the wind and stress profiles.  相似文献   

12.
Solutions to the sheared Fickian advection–diffusion equation in a half-space with arbitrary surface source are given using a ‘transfer function’ method. The method uses Fourier transforms in two horizontal coordinates and time, along with complex Airy functions in the vertical coordinate. Surface deposition and tracer decay are included in the formulation. ‘Puff’ and steady ‘plume’ solutions are compared with Saffman’s moment formulae. The inclusion of a decay rate factor (α) allows the average tracer age to be computed from steady state solutions for concentration C(x, y, z) according to Age = − dln C/dα. A comparison between the puff centroid formula of Saffman and plume Age computations confirms that shear causes tracer puffs to accelerate horizontally as they diffuse upward into a different wind regime. In forward shear, tracer ages are younger than in unsheared flow but the range of ages is greater due to the existence of a high fast pathway and a low slow pathway. In reverse shear, concentrations, ages and the range of ages all rise markedly near the source. Large tracer age suggests that some tracer has taken a very distant path involving a low-level outbound trip and a high-level return. The effect of surface deposition is to reduce the influence of the distant path. In the case of reverse shear, deposition makes the tracer younger. In a turning wind, the time needed to reach a given radius increases due to the curved path of the plume.  相似文献   

13.
Data collected during the SHEBA and CASES-99 field programs are employed to examine the flux–gradient relationship for wind speed and temperature in the stably stratified boundary layer. The gradient-based and flux-based similarity functions are assessed in terms of the Richardson number Ri and the stability parameter z/Λ*, z being height and Λ* the local Obukhov length. The resulting functions are expressed in an analytical form, which is essentially unaffected by self-correlation, when thermal stratification is strong. Turbulence within the stably stratified boundary layer is classified into four regimes: “nearly-neutral” (0 < z/Λ* < 0.02), “weakly-stable” (0.02 < z/Λ* < 0.6), “very-stable” (0.6 < z/Λ* < 50), and “extremely-stable” (z/Λ* > 50). The flux-based similarity functions for gradients are constant in “nearly-neutral” conditions. In the “very-stable” regime, the dimensionless gradients are exponential, and proportional to (z/Λ*)3/5. The existence of scaling laws in “extremely-stable” conditions is doubtful. The Prandtl number Pr decreases from 0.9 in nearly-neutral conditions and to about 0.7 in the very-stable regime. The necessary condition for the presence of steady-state turbulence is Ri < 0.7.  相似文献   

14.
The geostrophic Ekman boundary layer for large Rossby number (Ro) has been investigated by exploring the role played by the mesolayer (intermediate layer) lying between the traditional inner and outer layers. It is shown that the velocity and Reynolds shear stress components in the inner layer (including the overlap region) are universal relations, explicitly independent of surface roughness. This universality of predictions has been supported by observations from experiment, field and direct numerical simulation (DNS) data for fully smooth, transitionally rough and fully rough surfaces. The maxima of Reynolds shear stresses have been shown to be located in the mesolayer of the Ekman boundary layer, whose scale corresponds to the inverse square root of the friction Rossby number. The composite wall-wake universal relations for geostrophic velocity profiles have been proposed, and the two wake functions of the outer layer have been estimated by an eddy viscosity closure model. The geostrophic drag and cross-isobaric angle predictions yield universal relations, which are also supported by extensive field, laboratory and DNS data. The proposed predictions for the geostrophic drag and the cross-isobaric angle compare well with data for Rossby number Ro ≥ 105. The data show low Rossby number effects for Ro < 105 and higher-order effects due to the mesolayer compare well with the data for Ro ≥ 103.  相似文献   

15.
The heterogeneous chemistry of nitryl chloride and nitryl bromide by salt containing solutions was studied as a function of temperature in the range from 275 to 293 K with the wetted-wall flowtube combined with FTIR and mass spectrometry detection. Uptake coefficients and values of the product Hk1/2 on these saline solutions have been determined. For nitryl halides interacting with NaI and NaBr solutions, the values of the product Hk1/2 are respectively 4384.7±326.7 and 103.1±18.7 M atm-1 s-1/2 for nitryl chloride at 275 K and 544.2±94.7 and 47.7±15.2 M atm-1 s-1/2 for nitryl bromide at 278 K. When reacting with NaI or NaBr solutions, these heterogeneous reactions release, as major products, the molecular forms of the halogen i.e., respectively I2 and Br2. A simplified reaction scheme explaining the formation of these products is presented and is inserted into a model simulating the chemistry in the marine boundary layer. The modelling effort showed Cl and BrO atoms concentrations up to 5×104 and 1.8×106 molecules cm-3 respectively, which are comparable to values actually measured in field campaigns.  相似文献   

16.
Atmospheric surface layer meteorological observations obtained from 20-m-high meteorological tower at Mangalore, situated along the west coast of India are used to estimate the surface layer scaling parameters of roughness length (z o) and drag coefficient (C D), surface layer fluxes of sensible heat and momentum. These parameters are computed using the simple flux–profile relationships under the framework of Monin–Obukhov (M–O) similarity theory. The estimated values of z o are higher (1.35–1.54 m) than the values reported in the literature (>0.4–0.9 m) probably due to the undulating topography surrounding the location. The magnitude of C D is high for low wind speed (<1.5 m s?1) and found to be in the range 0.005–0.03. The variations of sensible heat fluxes (SHF) and momentum fluxes are also discussed. Relatively high fluxes of heat and momentum are observed during typical days on 26–27 February 2004 and 10–11 April 2004 due to the daytime unstable atmospheric conditions. Stable or near neutral conditions prevail after 1700 h IST with negative SHF. A mesoscale model PSU/NCAR MM5 is run using a high-resolution (1 km) grid over the study region to examine the influence of complex topography on the surface layer parameters and the simulated fluxes are compared with estimated values. Spatial variations of the frictional velocity (u *), C D, surface fluxes, planetary boundary layer (PBL) height and surface winds are noticed according to the topographic variations in the simulation.  相似文献   

17.
The boundary currents over the Western Australian continental shelf and slope consist of the poleward flowing Leeuwin Current (LC) and the equatorward flowing Leeuwin Undercurrent (LUC). Key properties of the LC are its poleward strengthening, deepening to the south, and shelfbreak intensification. The alongshore flow reverses direction below about 300 m, forming the LUC at greater depths. To investigate the processes that cause these features, we obtain solutions to an idealized, regional ocean model of the South Indian Ocean. Solutions are forced by relaxing surface density to a prescribed, meridionally varying density profile ρ*(y) with a timescale of δt. In addition, vertical diffusion is intensified near the ocean surface. This diffusion establishes the minimum thickness over which density is well-mixed. We define this thickness as the “upper layer”. Solutions are obtained with and without a continental shelf and slope off Western Australia and for a range of values of δt and mixing parameters. Within this upper layer, there is a meridional density gradient that balances a near-surface, eastward geostrophic flow. The eastward current downwells near the eastern boundary, leading to westward flow at depth. The upper layer's meridional structure and zonal currents crucially depend on coastal processes, including the presence of topography near the eastern boundary. Kelvin waves inhibit the upper layer from deepening at the coast. Rossby waves propagate the coastal density structure offshore, hence modifying the interior currents. A comparison of the solutions with or without a continental shelf and slope demonstrate that topographic trapping of Rossby waves is a necessary process for maintaining realistic eastern boundary current speeds. Significant poleward speeds occur only onshore of where the upper layer intersects the slope, that is, at a grounding line. Its poleward transport increases when surface-enhanced vertical mixing is applied over a greater depth. When the timescale δt is sufficiently short, the poleward current is nearly barotropic. The current's spatial structure over the shelf is controlled by horizontal mixing, having the structure of a Munk layer. Increasing vertical diffusion deepens the upper layer thickness and strengthens the alongshore current speed. Bottom drag leads to an offshore flow along the bottom, reducing the net onshore transport and weakening the current's poleward acceleration. When δt is long, poleward advection of buoyancy forms a density front near the shelf break, intensifying poleward speeds near the surface. With bottom drag, a bottom Ekman flow advects density offshore, shifting the jet core offshore of the shelf break. The resulting cross-shelf density gradient reverses the meridional current's direction at depth, leading to an equatorward undercurrent.  相似文献   

18.
The performance of the general bulk formulation in estimating sensible heat flux at Nigerian Micrometeorological Experimental site was assessed. Reliable sensible heat flux was estimated with the use of accurate diurnal values of transfer coefficient of sensible heat. The performances of one α, two β and a modified α formulations in the estimation of latent heat flux were also assessed at the station. The Lee and Pielke (β), modified Kondo (α), Jacquemin and Noilhan (α) and Noilhan and Planton (β) parameterizations gave good estimation of latent heat flux. The coefficient of determination (R 2) of the models between measured and estimated values were greater than 0.7. Low diurnal mean absolute error and root mean squared error values were found between measured and estimated fluxes. All the parameterizations gave reliable latent heat flux when diurnal values of transfer coefficients of moisture were used.  相似文献   

19.
A Lagrangian statistical-trajectory model based on a Markov chain relation is used to investigate vertical dispersion from elevated sources into the neutral planetary boundary layer. The model is fully two-dimensional, in that both vertical and longitudinal velocity fluctuations, and their correlation, are simulated explicitly. The best observational information currently available is used to characterize the mean and turbulent structure of the neutral boundary layer. In particular, a realistic vertical profile of the Lagrangian integral time scale is proposed, based partly on a review of direct measurements and partly on a comparison of the model predictions with published diffusion data. The model predictions are shown to agree well with a variety of dispersion observations. The model is used to study vertical diffusion as a function of release height H, friction velocity u* and surface roughness z 0 for downwind distances up to 10 km from the source. The equivalent Gaussian dispersion parameter Σ z is shown to decrease slightly with an increase in H, and to increase with increases in z 0 or u*. It is demonstrated that relationships valid in a field of homogeneous turbulence can be applied to vertical dispersion in the atmosphere if the release occurs above the region of strongest gradients in the mean and turbulent parameters. Scaling in terms of the standard deviation in elevation angle of the wind at the release point leads to a universal curve which provides accurate estimates of Σ z over a wide range of values of H, z 0 and the meteorological parameters.  相似文献   

20.
Integral Turbulence Statistics Over a Central European City Centre   总被引：2，自引：2，他引：0  
Atmospheric measurements over 5 years (2005–2010) at two sites in ?ód?, central Poland have been analyzed to develop a better understanding of turbulence in urban areas. Fast response wind velocity, temperature, humidity and CO2 concentration were measured using sonic anemometers and gas analyzers, placed on narrow masts at 37 and 42 m above the ground. The measurements were used to calculate standard deviations of each parameter, and were then normalized according to local Monin–Obukhov similarity theory and plotted as a function of stability parameter ζ = z′/L. Results for the wind components show typical scaling with a power law with exponent ±1/3 in the free convection limit, and that approaches a constant value close to neutral stratification. For stable conditions, the constant value in the neutral limit remains the same for stability parameters lower than 0.1–0.2, then increases. The normalized standard deviation of temperature fits the ?1/3 law in the free convection limit, approaching a constant value within a stable limit. However, it exhibits hyperbolic characteristics for close to neutral stratification. The normalized standard deviations for humidity and CO2 concentration exhibit scaling similar to the wind components in the unstable regime and remain constant in the stable domain. The results for the wind components and for temperature are in the range of various functions found in other studies. The absolute values for humidity and CO2 concentration seem to be slightly higher, but only single examples of such investigations can be found in the literature.  相似文献   

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

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