The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincaré variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator Mean Exponential Growth factor of Nearby Orbits (MEGNO), a measure of the maximal Lyapunov exponent, that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best-fitting solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best-fitting solutions are found in dynamically active region of the phase space. The long-term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myr of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behaviour.     

Potentials and limitations of modelling vertical distributions of root water uptake of an Austrian pine forest on a sandy soil

Root water uptake patterns are often studied with simulation models of the unsaturated soil water flow, as they are difficult to measure directly. Calibration of these models is not straightforward and causes uncertainties in simulated uptake distributions. In this paper we study how uncertainties in the calibration of the SWIF model affect uncertainty intervals in simulated uptake patterns of an Austrian pine stand (Pinus nigra var. nigra) on a sandy soil. After calibrating and validating SWIF with a large data set of more than 125 000 measured soil water contents over a three year period, uncertainty ranges in simulated soil water dynamics and root water uptake distributions were estimated with a Monte Carlo analysis. In general, uncertainties in root uptake patterns were small (typically <2 10⁻⁴ m³ m⁻³ day⁻¹) and were higher for trees with a shallow rooting system (0·8 m) than for trees with a deep rooting system (2·5 m). Uncertainties arose mainly from uncertainties in simulated soil water fluxes and from variations in the reduction of uptake during periods of drought. Uncertainties in soil water contents were far higher (typically 0·01 m³ m⁻³) than uncertainties in uptake, illustrating that uncertainties in uptake parameters and those in the distribution of water uptake hardly affect the modelling of soil water dynamics. Root water uptake models should therefore be validated against measured uptake distributions, which can be determined on sandy soils during dry periods with a high water use when soil fluxes are negligible to uptake. Copyright © 2000 John Wiley & Sons, Ltd.     

聚焦波作用下平面网衣结构的水动力特性研究

由于沿海区域的限制以及愈加严重的环境污染,渔业养殖正从近海走向深远海。深远海海域的海况更加恶劣,给养殖装备的设计与性能评估带来新的挑战。为解决该问题,对极端波浪与养殖装备网衣结构的相互作用开展研究。基于waves2Foam建立数值波浪水池,极端波浪模拟采用基于NewWave理论的聚焦波模型,网衣结构模拟采用多孔介质模型,并通过与Morison模型计算的网衣受力等效分析,获得多孔介质模拟网衣结构阻力系数的直接估计方法。然后将多孔介质模型嵌入waves2Foam中,开展聚焦波与网衣结构相互作用的数值模拟,同时开展水槽试验,验证数值模拟的准确性。基于数值模拟结果,系统地分析了不同网衣密实度及不同波浪参数下网衣结构的升阻力特性以及网衣结构对波浪场的扰动规律。研究表明：聚焦波波峰幅值和网衣密实度对网衣结构的升阻力影响较大,且升力峰值出现在阻力为0的时刻;网衣结构对聚焦波的时空演化特性有影响,改变了聚焦波波形。     

非静力中β尺度模式对渤海地区水平能见度的模拟

运用各方面的观测资料,和非流体静力学中尺度模式第5版本对2002年12月2日发生在环渤海地区的一次低能见度天气事件进行了分析研究.利用MM5模式对本次事件进行了模拟,并计算了大气的水平能见度分布.模拟结果显示当能见度低于某一低值即1km,同时大面积爆发时,模拟结果与实况比较吻合;当能见度处于轻雾并零星分布状态时,模式还不能非常灵敏地将其有效模拟出来.     

Use of Large-Eddy Simulation for the bed shear stress estimation over a dune

Environmental flows are generally characterized by complex bed morphology and high current speeds. Such configurations favor the formation of vortex structures that strongly affect hydrody-namics and sediment transport. Large-Eddy Simulation (LES) enables investigation of the dynam-ics of the largest turbulence scales and, thanks to enhanced calculation resources, has now become applicable for simulating environmental flows. In this paper, a LES approach is developed in a CFD code (TELEMAC-3D), which was originally developed to simulate free surface flows using RANS methods. The present developments involve implementing subgrid models, boundary con-ditions and numerical schemes suitable for LES. The LES version of TELEMAC-3D was validated by comparing results on the model with experimental data for flow past a cylinder. Then, the model was applied to a test case representing flow over dunes. After validating the hydrodynamics, the model was used to assess the bottom shear stress, using both a RANS and a LES approach. Com-parison highlighted the potential contribution of LES to investigating the hydrodynamic forces acting on the bottom.     

Effects of a Building's Density on Flow in Urban Areas

The effects of a building's density on urban flows are investigated using a CFD model with the RNG k-ε turbulence closure scheme. Twenty-seven cases with different building's density parameters (e.g., building and street-canyon aspect ratios) are numerically simulated. As the building's density parameters vary, different flow regimes appear. When the street canyon is relatively narrow and high, two counter-rotating vortices in the vertical direction are generated. The wind speed along streets is mainly affected by the building's length. However, it is very difficult to find or generalize the characteristics of the street-canyon flows in terms of a single building's density parameter. This is because the complicated flow patterns appear due to the variation of the vortex structure and vortex number. Volume-averaged vorticity magnitude is a very good indicator to reflect the flow characteristics despite the strong dependency of flows on the variation of the building's density parameters. Multi-linear regression shows that the volume-averaged vorticity magnitude is a strong function of the building's length and the street-canyon width. The increase in the building's length decreases the vorticity of the street-canyon flow, while, the increase in the street-canyon width increases the vorticity.     

21世纪重庆中雨以上天数的预估分析

利用用于IPCC-AR4的全球气候模式产品,验证其对重庆地区极端降水指数中雨以上天数(dR10)模拟能力的基础上,对模拟能力较好的模式进行组合,预估高(A2)、中(A1B)、低(B1)三种排放情景下未来21世纪重庆地区dR10的变化。不同排放情景下未来重庆dR10的变化不太一致。与目前气候(1980—1999年)相比,不同情景下未来21世纪重庆地区dR10在多数时期将可能减少。21世纪的后90a(2011—2100年),A2情景下重庆dR10减少最多,平均减少1.3d;3种情景平均将减少0.5d。21世纪初期(2011—2040年)、中期(2041—2070年)和后期(2071—2100年),A2情景下重庆dR10减少都最多,分别平均减少1.6d、1.6d和0.7d;3种情景平均分别减少0.8d、0.6d和0.1d。     

Study of near-surface models for large-eddy simulations of a neutrally stratified atmospheric boundary layer

In large-eddy simulations (LES) of the atmospheric boundary layer (ABL), near-surface models are often used to supplement subgrid-scale (SGS) turbulent stresses when a major fraction of the energetic scales within the surface layer cannot be resolved with the temporal and spatial resolution at hand. In this study, we investigate the performance of both dynamic and non-dynamic eddy viscosity models coupled with near-surface models in simulations of a neutrally stratified ABL. Two near-surface models that are commonly used in LES of the atmospheric boundary layer are considered. Additionally, a hybrid Reynolds- averaged/LES eddy viscosity model is presented, which uses Prandtl’s mixing length model in the vicinity of the surface, and blends in with the dynamic Smagorinsky model away from the surface. Present simulations show that significant portions of the modelled turbulent stresses are generated by the near-surface models, and they play a dominant role in capturing the expected logarithmic wind profile. Visualizations of the instantaneous vorticity field reveal that flow structures in the vicinity of the surface depend on the choice of the near-surface model. Among the three near-surface models studied, the hybrid eddy viscosity model gives the closest agreement with the logarithmic wind profile in the surface layer. It is also observed that high levels of resolved turbulence stresses can be maintained with the so-called canopy stress model while producing good agreement with the logarithmic wind profile.     

基于CFD的地震液化研究新进展

综述了近年来关于液化土体流体动力学特征的试验发展状况,以及基于计算流体动力学（简称CFD）的地震液化数值模拟现状,重点介绍了目前比较活跃的可以较高精度模拟液化土体流动状态的三次伪质点数值方法（简称CIP法）。通过对CFD和传统固体力学在地震液化研究中的应用比较,指出了应用CFD的三大优势,即土体大变形问题、液化土体参数分析以及液化土体中结构物的变形应用CFD分析,均可获得较好的结果。进一步提出,在地震液化应用中,未来CFD的发展应该考虑整合液化前的土体性状研究和地震液化中桩-土-结构物的综合分析。     

Zebra and boxwork fabrics in hydrothermal dolomites of northern Canada: Indicators for dilational fracturing,dissolution or in situ replacement?

Zebra and boxwork fabrics in Manetoe and Presqu'ile hydrothermal dolomites provide critical data for discrimination between hypotheses for their origin. In a unique occurrence, sets of weakly imbricated, decimetre‐sized curvilinear lenticles of white dolospar are developed within sub‐metre sized tongues of greyish‐white dolomitized crinoidal packstone and grainstone. This nascent zebra fabric has developed along wavy to nodular intra‐bed sedimentary partings. Primary bedding can be traced through dolomitized masses containing these fabrics, into the surrounding undolomitized limestone. All zebra and boxwork fabrics are confined within single beds, consistent with an early burial time of origin. No change in bed thickness occurs across zebra or boxwork fabrics within beds; this indicates a dissolutional, rather than a dilational, origin for the creation of pore space partly cemented with white dolospar. The dolomitized groundmass displays an abruptly gradational or sharp transition to centripetal saddle dolomite cement, which partially to nearly completely, occludes vugs. Circulation of geothermally heated hypersaline Devonian Elk Point basinal sea water brines led to dissolution and replacive dolomitization of limestone adjacent to vugs, and near‐contemporaneous precipitation of white dolospar within dissolutional vugs, consistent with geochemical simulations. The presence of downward‐extending galleries of white dolospar‐cemented solution‐collapse breccia provides further support for an interpretation of regional thermally driven convection of hydrothermal evaporative sea water brines across a broad area of northern Canada. The absence of gases under pressure in vacuoles within fluid inclusions from Manetoe and Presqu'ile dolospars is also more consistent with a dissolutional and contemporaneous dolomitization origin for these fabrics, rather than an origin involving dilational fracturing for space creation and dolospar precipitation. The ubiquitous presence of zebra and boxwork fabrics in hydrothermal dolomite reservoirs indicates that they are not confined fault zones and instead occur wherever precursor limestones had relatively greater porosity and permeability prior to dolomitization.