Large-Eddy Simulations of the Baroclinic Mixed Layer

The effects of baroclinicity, imposed on the dry mixed layer by the presenceof large-scale, horizontal temperature gradients, have been investigated basedon a large-eddy simulation model. The purpose of this investigation is to examinesimultaneous impacts of thermal stratification and shear in the atmospheric boundarylayer. For this purpose, five cases are considered – one barotropic, and four baroclinic.Based on the performed simulations, a new parametrization of the interfacial layer isproposed. The parameterization employs new interfacial scaling, which is valid at thetop of the mixed layer. In terms of new scales, dimensionless moments characterizingturbulence at the top of the shearless mixed layer are universal constants. In the shearedcase, dimensionless statistics of turbulence are shown to be functions of the interfacialRichardson number.     

Universal hydrofracturing algorithm for shear-thinning fluids: Particle velocity based simulation

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.     

Ensemble simulation of land evapotranspiration in China based on a multi-forcing and multi-model approach

In order to reduce the uncertainty of offline land surface model(LSM) simulations of land evapotranspiration(ET), we used ensemble simulations based on three meteorological forcing datasets [Princeton, ITPCAS(Institute of Tibetan Plateau Research, Chinese Academy of Sciences), Qian] and four LSMs(BATS, VIC, CLM3.0 and CLM3.5), to explore the trends and spatiotemporal characteristics of ET, as well as the spatiotemporal pattern of ET in response to climate factors over mainland China during 1982–2007. The results showed that various simulations of each member and their arithmetic mean(Ens Mean) could capture the spatial distribution and seasonal pattern of ET sufficiently well, where they exhibited more significant spatial and seasonal variation in the ET compared with observation-based ET estimates(Obs MTE). For the mean annual ET, we found that the BATS forced by Princeton forcing overestimated the annual mean ET compared with Obs MTE for most of the basins in China, whereas the VIC forced by Princeton forcing showed underestimations. By contrast, the Ens Mean was closer to Obs MTE, although the results were underestimated over Southeast China. Furthermore, both the Obs MTE and Ens Mean exhibited a significant increasing trend during 1982–98; whereas after 1998, when the last big EI Ni ?no event occurred, the Ens Mean tended to decrease significantly between 1999 and 2007, although the change was not significant for Obs MTE. Changes in air temperature and shortwave radiation played key roles in the long-term variation in ET over the humid area of China, but precipitation mainly controlled the long-term variation in ET in arid and semi-arid areas of China.     

Probabilistic stability analyses of undrained slopes by 3D random fields and finite element methods

A long slope consisting of spatially random soils is a common geographical feature. This paper examined the necessity of three-dimensional(3 D) analysis when dealing with slope with full randomness in soil properties. Although 3 D random finite element analysis can well reflect the spatial variability of soil properties, it is often time-consuming for probabilistic stability analysis. For this reason, we also examined the least advantageous(or most pessimistic) cross-section of the studied slope. The concept of"most pessimistic" refers to the minimal cross-sectional average of undrained shear strength. The selection of the most pessimistic section is achievable by simulating the undrained shear strength as a 3 D random field. Random finite element analysis results suggest that two-dimensional(2 D) plane strain analysis based the most pessimistic cross-section generally provides a more conservative result than the corresponding full 3 D analysis. The level of conservativeness is around 15% on average. This result may have engineering implications for slope design where computationally tractable 2 D analyses based on the procedure proposed in this study could ensure conservative results.     

Modeling of iceberg drift in the marginal ice zone of the Barents Sea

The paper is focused on the analysis of the drift of tabular iceberg observed in 2009 in the marginal ice zone of the North-West Barents Sea. Momentum balance equations are derived from the Kirchhoff equations describing plane motion of solid body in an ideal fluid. Field works performed on the drift iceberg and on the drift icenear the iceberg are described. Results of the field works and numerical simulations of the iceberg drift and rotation are performed and discussed. It is shown that acceleration of water flow around the iceberg has visible influence on the iceberg drift. Kinetic energy balance of drift iceberg is used to estimate the forces applied to the iceberg by the drift ice.     

Numerical modeling investigation on turbulent oscillatory flow over a plane rough bed composed by randomly arrayed particles

A three-dimensional numerical model is established to simulate the turbulent oscillatory boundary layer over a fixed and rough bed composed by randomly arrayed solid spheres based on the lattice Boltzmann method and the large eddy simulation model.The equivalent roughness height,the location of the theoretical bed and the time variation of the friction velocity are investigated using the log-fit method.The time series of turbulent intensity and Reynolds stress are also investigated.The equivalent roughness height of cases with Reynolds numbers of 1×10~4–6×10~4 is approximately 2.81 d(grain size).The time variation of the friction velocity in an oscillatory cycle exhibits sinusoidal-like behavior.The friction factor depends on the relative roughness in the rough turbulent regime,and the pattern of solid particles arrayed as the rough bed in the numerical simulations has no obvious effect on the friction factor.     

早新生代温室气候及冰期气候转型的模拟研究

早新生代是地质史上最后一个温室气候期,随后南极冰盖形成,地球进入到晚新生代冰期。温室气候的成因和冰期气候转型的机制一直是国际相关学界关注的问题。评述国际上对此开展的古气候模拟,反映了早新生代温室气候受到了海洋和大陆的地理位置、暖海洋温盐环流和海洋热输送、太阳辐射和大气CO2浓度变化的作用和影响。古气候模拟还反映了早新生代温室气候转向冰期气候,受到了大洋通道改变和高原构造隆起、大气成分变化以及海陆生态系相互的作用和反馈。这些古气候模拟试验锁定在气候变化的关键时段和重要驱动因子,对测试地球内外驱动力和地球各圈层反馈作用提供了重要的科学依据;温室气候以及趋向冰期气候的模拟研究对探讨气候变化内在机制、预测未来气候具有重要意义。      

Limit states and failure mechanisms of viscous dampers and the implications for large earthquakes

Fluid viscous dampers are used to control story drifts and member forces in structures during earthquake events. These elements provide satisfactory performance at the design‐level or maximum considered earthquake. However, buildings using fluid viscous dampers have not been subjected to very large earthquakes with intensities greater than the design and maximum considered events. Furthermore, an extensive database of viscous damper performance during large seismic events does not exist. To address these issues, a comprehensive analytical and experimental investigation was conducted to determine the performance of damped structures subjected to large earthquakes. A critical component of this research was the development and verification of a detailed viscous damper mathematical model that incorporates limit states. The development of this model and the laboratory and simulation results conclude good correlation with the new model and the damper limit states and provide superior results compared with the typical damper model when considering near collapse evaluation of structures. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulating the effects of spatial configurations of agricultural ditch drainage networks on surface runoff from agricultural catchments

The study of runoff is a crucial issue because it is closely related to flooding, water quality and erosion. In cultivated catchments, agricultural ditch drainage networks are known to influence runoff. As anthropogenic elements, agricultural ditch drainage networks can therefore be altered to better manage surface runoff in cultivated catchments. However, the relationship between the spatial configuration, i.e. the density and the topology, of agricultural ditch drainage networks and surface runoff in cultivated catchments is not understood. We studied this relationship by using a random network simulator that was coupled to a distributed hydrological model. The simulations explored a large variety of spatial configurations corresponding to a thousand stochastic agricultural ditch drainage networks on a 6.4 km² Mediterranean cultivated catchment. Next, several distributed hydrological functions were used to compute water flow paths and runoff for each simulation. The results showed that (i) denser networks increased the drained volume and the peak discharge and decreased hillslopes runoff, (ii) greater network density did not affect the surface runoff any further above a given network density, (iii) the correlation between network density and runoff was weaker for small subcatchments (< 2 km²) where the variability in the drained area that resulted from changes in agricultural ditch drainage networks increased the variability of runoff and (iv) the actual agricultural ditch drainage network appeared to be well optimized for managing runoff as compared with the simulated networks. Finally, our results highlighted the role of agricultural ditch drainage networks in intercepting and decreasing overland flow on hillslopes and increasing runoff in drainage networks. Copyright © 2011 John Wiley & Sons, Ltd.     

大气重力波触发的Rayleigh－Taylor不稳定性的时空演变

对大气重力波触发的赤道电离层Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性的时空演变进行了数值模拟．结果表明，重力波能在Ｆ区底部触发Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性，这种等离子体不稳定扰动穿过Ｆ峰，到达Ｆ区顶部，形成等离子体泡结构．等离子体泡出现向西倾斜和分岔等特征．当夜Ｆ区较高时，产生等离子体泡的时间仅１９００ｓ．数值模拟结果证实了重力波触发．Ｒａｙｌｅｉｇｈ－Ｔａｙｌｏｒ不稳定性的理论，解释了大量电离层观测现象．