Large-Eddy Simulations of the Stable Boundary Layer Using the Standard Kolmogorov Theory: Range of Applicability

Large-eddy simulations (LES) of the Stable Atmospheric Boundary Layer (SBL) are difficult because the turbulence is not isotropic for strong stratification and the Kolmogorov theory might be no longer valid. This fact compells us to work on modifications to the subgrid turbulence schemes, although there is not any widely accepted theory on anisotropic turbulence. In this work, a LES model is used to see what range of stable stratification can still be simulated with a subgrid turbulence scheme using the Kolmogorov theory for the dissipation. Twenty simulations of increasing stability have been performed using a horizontal resolution of 5 m. The model is able to simulate weakly and moderately stable conditions and experiences runaway cooling for strong stability. The goodness of the successful simulations is inspected through comparison to observations from the experimental campaigns SABLES-98 and CASES-99. Other supplementary tests have been performed on the resolution and the surface boundary condition.     

印度洋、南海和东南沿海海温异常影响江淮流域6~7月降水量的分析及数值模拟

利用统计分析和CCM3数值模式研究了印度洋、南海和东南沿海同期海温异常对江淮梅雨的影响。统计分析结果表明:印度洋、南海和东南沿海的正 (负) 海温异常对应于江淮流域的多 (少) 雨。对OLR资料分析表明:当印度洋OLR为正 (负) 距平时, 江淮流域OLR也为正 (负) 距平, 但南海和东南沿海OLR多为负 (正) 距平。数值模拟结果表明:当印度洋为正海温异常时, 中高纬地区的阻高偏强, 江淮流域为多雨区; 印度洋为负海温异常时, 东亚沿海出现江淮流域旱年的环流形势, 与统计事实相符合。     

Numerical modelling of flow structures over idealized transverse aeolian dunes of varying geometry

A Computational Fluid Dynamics (CFD) model (PHOENICS™ 3.5) previously validated for wind tunnel measurements is used to simulate the streamwise and vertical velocity flow fields over idealized transverse dunes of varying height (h) and stoss slope basal length (L). The model accurately reproduced patterns of: flow deceleration at the dune toe; stoss flow acceleration; vertical lift in the crest region; lee-side flow separation, re-attachment and reversal; and flow recovery distance. Results indicate that the flow field over transverse dunes is particularly sensitive to changes in dune height, with an increase in height resulting in flow deceleration at the toe, streamwise acceleration and vertical lift at the crest, and an increase in the extent of, and strength of reversed flows within, the lee-side separation cell. In general, the length of the separation zone varied from 3 to 15 h from the crest and increased over taller, steeper dunes. Similarly, the flow recovery distance ranged from 45 to >75 h and was more sensitive to changes in dune height. For the range of dune shapes investigated in this study, the differing effects of height and stoss slope length raise questions regarding the applicability of dune aspect ratio as a parameter for explaining airflow over transverse dunes. Evidence is also provided to support existing research on: streamline curvature and the maintenance of sand transport in the toe region; vertical lift in the crest region and its effect on grainfall delivery; relations between the turbulent shear layer and downward forcing of flow re-attachment; and extended flow recovery distances beyond the separation cell. Field validation is required to test these findings in natural settings. Future applications of the model will characterize turbulence and shear stress fields, examine the effects of more complex isolated dune forms and investigate flow over multiple dunes.     

基于均值平移与模拟泛洪的图像分割算法

针对图像处理的特点，综合考虑执行速度与健壮性要求，提出了一种顺序搜索LUV空间候选点的均值平移实现方法。为了克服传统的基于特征空间分析的图像分割方法对像素点空间关系考虑不够充分的缺陷，借鉴Vincent ＆ Soille沉浸模拟分水岭算法思想，构造了一种与均值平移方法相适应的新的模拟泛洪方法执行图像分割。新方法比较符合人的主观感知，且易于容纳更多的特征以满足特定的图像分割需求。对实验用彩色图像的分割结果表明了该方法的有效性。     

Quantifying Uncertainty in Mineral Resources by Use of Classification Schemes and Conditional Simulations

Mineral inventory determination consists of estimating the amount of mineral resources on a block-by-block basis and classifying individual blocks into categories with increasing level of geologic confidence. Such classification is a crucial issue for mining companies, investors, financial institutions, and authorities, but it remains subject to some confusion because of the wide variations in methodologies and the lack of standardized procedures. The first part of this paper considers some of the criteria used to classify resources in practice and their impact through a sensitivity study using data from a Chilean porphyry copper deposit. Five classification criteria are compared and evaluated, namely: Search neighborhoods, absolute and relative kriging variances, absolute and relative conditional simulation variances. It is shown that some classification criteria either favor or penalize the high-grade areas if the grade distribution presents a proportional effect. In the second part of the paper, conditional simulations are used to quantify the uncertainty on the overall mineral resources. This approach is promising for risk analysis and decision-making. Unlike linear kriging, simulations allow inclusion of a cutoff grade in the calculation of the resources and also provide measures of their joint uncertainty over production volumes.     

CFD方法在街巷气象场模拟和预测中的应用

为开发建立米级分辨率的城市街巷气象场模拟预测方案，以PHOENICS为例介绍了计算流体力学（CFD）方法的原理，讨论了该方法应用于街巷气象要素数值计算的可行性。提出了一种基于CFD模型与高分辨率数值天气预报模式嵌套的街渠气象场模拟预测方案，并对边界嵌套问题作了多种敏感性分析，表明嵌套对街渠气象场计算具有必要性和可行性；另外，还讨论了分辨率对计算结果和计算耗时的影响；最后，用实例介绍了本方案的应用，并讨论了本方法的应用前景。     

Morphological development of the retention basin “Hartheim”: A case study

Flow in rivers and on floodplains is complex as it is affected by several interconnected factors such as topography, sediment transport and vegetation characteristics. The resulting processes are explained by the measure “Hartheim” planned for retention purposes at the Upper Rhine river. On the basis of existing formulas and instruments it is demonstrated that a good estimation of the development of the measure is possible. The proposed procedure is a useful tool for estimating morphological developments of restored river sections.     

Block-effective macrodispersion for numerical simulations of sorbing solute transport in heterogeneous porous formations

Heterogeneity is prevalent in aquifers and has an enormous impact on contaminant transport in groundwater. Numerical simulations are an effective way to deal with heterogeneity directly by assigning different hydraulic property values to each numerical grid block. Because hydraulic properties vary on different scales, but they cannot be sampled exhaustively and the number of numerical grid blocks is limited by computational considerations, the dispersive effects of unmodeled heterogeneity need to be accounted for. Dispersion tensors can be used to model the dispersion caused by unmodeled heterogeneity. The concept of block-effective macrodispersion tensors for modeling the effects of small-scale variability on solute transport introduced by Rubin et al. [Rubin Y, Sun A, Maxwell R, Bellin A. The concept of block-effective macrodispersivity and a unified approach for grid-scale- and plume-scale-dependent transport. J Fluid Mech 1999;395:161–80] is extended in this paper for use with reactive solutes. The tensors are derived for reactive solutes with spatially variable retardation factors and for solutes experiencing spatially uniform rate-limited sorption. The longitudinal block-effective macrodispersion coefficient is largest for perfect negative correlation between the log-hydraulic conductivity and the retardation factor. Because dispersion tensors, as they are usually implemented in numerical simulations, produce symmetric spreading, the applicability of the concept depends on the portion of the plume asymmetry caused by small-scale variability. The presented results show that the concept is applicable for rate-limited sorption for block sizes of one and two integral scales.     

Geospatial–temporal dependence among weekly precipitation extremes with applications to observations and climate model simulations in South America

A quantification of the spatio-temporal dependence among precipitation extremes is important for investigating the properties of intense storms as well as flood or flash-flood related hazards. Extreme value theory has been widely applied to the hydrologic sciences and hydraulic engineering. However, rigorous approaches to quantify dependence structures among extreme values in space and time have not been reported in the literature. Previous researchers have quantified the dependence among extreme values through the concept of (pairwise bivariate) tail dependence coefficients. For estimation of the tail dependence coefficients, we apply a recently developed method [Kuhn G. On dependence and extremes. PhD thesis (Advisor: C. Klüppelberg), Munich University of Technology, 2006] which utilized the multivariate tail dependence function of a subclass of elliptical copulas. This study extends the previous approach in the context of space and time by considering pairs of spatial grids in South America and quantifying the dependence among precipitation extremes based on the time series at each spatial grid. In addition, Kendall’s τ is used to estimate the pairwise copula correlation (for an elliptical copula) of precipitation between all grids in South America. The geospatial–temporal dependence measures are applied to precipitation observations from 1940 to 2005 as well as simulations from the Community Climate System Model version 3 (CCSM3) for 1940–2099. New insights are obtained regarding the spatio-temporal dependence structures for precipitation over South America both with regard to correlation as well as tail dependence.     

The decomposition and simulation of the longitudinal circulation in a coastal lagoon

Current measurements during a 32-day study period in late spring, 1977, are used to quantify the magnitude and relative importance of tidal and wind-driven motion in the interior of the Indian River lagoon, on the Atlantic coast of Florida. Harmonic analysis of the total longitudinal flow along the axis of the lagoon isolates the tidal component of the current; non-tidal flow is revealed by subtracting the tidal current from the total current, and making corrections for non-linear relationships between the current and both surface wind stress and bottom friction. A one-layer, one-dimensional model is developed to simulate wind drift. A quadratic bottom friction term with a drag coefficient of 15 × 10^?3 gives results which compare most favourably with observations. Results indicate that tidal forcing explains approximately 45% of the total variance at the study site, 25 km from the nearest inlet. Local wind forcing accounts for 44% of the non-tidal flow. The remainder of the variance is attributed to freshwater outflow through the lagoon and non-local forcing.