计算断层自发破裂过程的离散网格划分研究

离散化网格的空间步长选取在各种数值算法中都是一个很受关注的问题,在全空间均匀介质模型和简单离散化方案下,利用边界积分方程方法研究在自发破裂求解过程中动力学参数组合Dc和Te (Dc为临界滑动弱化位移,Te为有效的断层破裂面上的初始应力)对计算网格划分的影响,初步得到了Dc和Te参数空间中的有效计算网格的选取规律,对合理、有效地运用边界积分方程方法计算地震震源的破裂过程具有重要的指导意义.      

Non-point source pollution estimation in the Pingqiao River Basin,China, using a spatial hydrograph-separation approach

For efficient and targeted management, this study demonstrates a recently developed non-point source (NPS) pollution model for a year-long estimation in the Pingqiao River Basin (22.3 km²) in China. This simple but physically reasonable model estimates NPS export in terms of land use by reflecting spatial hydrological features and source runoff measurements under different land-use types. The NPS export was separately analysed by a distributed hydrological model, a spatial hydrograph-separation technique, and an empirical water quality sub-model. Simulation results suggest that 57 890 kg of total nitrogen (TN) and 1148 kg of total phosphorus (TP) were delivered. The results, validated with observed stream concentrations, show relative errors of 23.3% for TN and 47.4% for TP. Countermeasures for urban areas (5.3% of total area) were prioritized because of the high contribution rate to TN (14.1%) and TP (26.2%) which is caused by the high degree of runoff (8.5%) and pollution source.     

Some remarks on cumulus parameterization

A simple parameterization of cumulonimbus convective heating is presented. The model is primarily based on preserving a moisture budget and on detraining cloud air at levels corresponding to the neutral buoyancy of the air converged at low levels. Results are compared with data from the western Pacific and GATE. Agreement is good. Suggestions are offered for improving the model and extending it to other regions.     

Numerical modeling of lock-exchange gravity/turbidity currents by a high-order upwinding combined compact difference scheme

This study presents two-dimensional direct numerical simulations for sediment-laden current with higher density propagating forward through a lighter ambient water.The incompressible NavierStokes equations including the buoyancy force for the density difference between the light and heavy fluids are solved by a finite difference scheme based on a structured mesh.The concentration transport equations are used to explore such rich transport phenomena as gravity and turbidity currents.Within the framework of an Upwinding Combined Compact finite Difference(UCCD)scheme,rigorous determination of weighting coefficients underlies the modified equation analysis and the minimization of the numerical modified wavenumber.This sixth-order UCCD scheme is implemented in a four-point grid stencil to approximate advection and diffusion terms in the concentration transport equations and the first-order derivative terms in the Navier-Stokes equations,which can greatly enhance convective stability and increase dispersive accuracy at the same time.The initial discontinuous concentration field is smoothed by solving a newly proposed Heaviside function to prevent numerical instabilities and unreasonable concentration values.A two-step projection method is then applied to obtain the velocity field.The numerical algorithm shows a satisfying ability to capture the generation,development,and dissipation of the Kelvin-Helmholz instabilities and turbulent billows at the interface between the current and the ambient fluid.The simulation results also are compared with the data in published literatures and good agreements are found to prove that the present numerical model can well reproduce the propagation,particle deposition,and mixing processes of lock-exchange gravity and turbidity currents.     

三维土-结构动力相互作用的一种时域直接分析方法

本文提出了一种分析三维土-结构动力相互作用的时域直接方法。该方法采用集中质量显式有限元和透射人工边界模拟无限域地基,通过编制的FORTRAN程序实现;采用ANSYS软件对上部结构进行建模分析,并通过FORTRAN程序对ANSYS软件的调用,实现了土与结构系统在地震作用下的整体分析。该方法为显隐式相结合的方法,地基和上部结构可采用不同的时间步距进行分析,可大大提高效率。通过两算例,验证了该方法的可行性。     

有限差分格式与物理量守恒性

为了成功地开展大气数值模拟,如何构造相应的离散格式是一项十分重要的工作.本文以大气球面浅水方程为例,从不同形式的微分方程组出发,分别构造了能量守恒格式、准辛格式和多守恒格式,并作相应的数值检验和比较,从中可以清楚地看到三类格式的优劣.这三类格式构造方法不同,所保持的物理守恒性不完全相同,数值计算时所需机时也不同,这些比较结果为实际工作的不同使用需求提供了选择依据.     

地球流体力学的研究与进展

简要介绍中国科学院大气物理研究所七十多年来在理论与计算地球流体力学方面的若干研究及其新的进展.在理论地球流体力学方面,介绍了长波动力学及线性稳定性问题、弱非线性理论及行星波动力学以及用Arnold方法(能量-Casimir方法)研究大气和海洋中各种流体运动的非线性稳定性问题的成果.此外,对扰动演变、扰动和基流相互作用及热带大气动力学中的第二类不稳定条件(CISK)也作了简要的介绍.在计算地球流体力学方面,主要内容包括:用物理观点和数学分析相结合的方法阐述了造成计算紊乱和计算不稳定的机理,论证计算稳定性、算     

Adaptive Multiscale Permeability Estimation

With multiscale permeability estimation one does not select parameterization prior to the estimation. Instead, one performs a hierarchical search for the right parameterization while solving a sequence of estimation problems with an increasing parameterization dimension. In some previous works on the subject, the same refinement is applied all over the porous medium. This may lead to over-parameterization, and subsequently, to unrealistic permeability estimates and excessive computational work. With adaptive multiscale permeability estimation, the new parameterization at an arbitrary stage in the estimation sequence is such that new degrees of freedom are not necessarily introduced all over the porous medium. The aim is to introduce new degrees of freedom only where it is warranted by the data. In this paper, we introduce a novel adaptive multiscale estimation. The approach is used to estimate absolute permeability from two-phase pressure data in several numerical examples.     

单程波方程偏移算子高阶辛格式

针对地震偏移算法中单程波算子中e指数项的近似，提出一种高阶辛格式算法. 与非辛格式的近似、一阶辛格式近似和两步法进行数值对比表明，文中提出的高阶辛格式具有更高的精度. 脉冲响应的计算也表明文中的近似方法具有降低频散效应和适应横向变速的能力. 用Marmousi模型验证了对复杂构造的成像的能力，实际资料计算表明本文方法适合于较大深度的稳定准确成像.     

Understanding of the Common Land Model performance for water and energy fluxes in a farmland during the growing season in Korea

The Common Land Model (CLM) is one of the most widely used land surface models (LSMs) due to the practicality of its simple parameterization scheme and its versatility in embracing a variety of field datasets. The improved assessment of land surface water and energy fluxes using CLM can be an alternative approach for understanding the complex land–atmosphere interactions in data‐limited regions. The understanding of water and energy cycles in a farmland is crucial because it is a dominant land feature in Korea and Asia. However, the applications of CLM to farmland in Korea are in paucity. The simulations of water and energy fluxes by CLM were conducted against those from the tower‐based measurements during the growing season of 2006 at the Haenam site (a farmland site) in Korea without optimization. According to the International Geosphere–Biosphere Programme (IGBP) land cover classification, a homogeneous cropland was selected initially for this study. Although the simulated soil moisture had a similar pattern to that of the observed, the former was relatively drier (at 0·1 m³ m^?3) than the latter. The simulated net radiation showed good agreement with the observed, with a root mean squared error (RMSE) of 41 W m^?2, whereas relatively large discrepancies between the simulation and observation were found in sensible (RMSE of 66 W m^?2) and latent (RMSE of 60 W m^?2) heat fluxes. On the basis of the sensitivity analysis, soil moisture was more receptive to land cover and soil texture parameterizations when compared to soil temperature and turbulent fluxes. Despite the uncertainty in the predictive capability of CLM employed without optimization, the initial performance of CLM suggests usefulness in a data‐limited heterogeneous farmland in Korea. Further studies are required to identify the controls on water and energy fluxes with an improved parameterization. Copyright © 2010 John Wiley & Sons, Ltd.