A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

关于多项式特征值的条件数

Tisseur(2000)研究了非齐次多项式特征值问题的条件数，所得结果的不足是对无穷特征值必须另外处理；Declieu和Tisseur(2003)运用微分几何方法研究了齐次多项式特征值问题的条件数，所得结果的不足之处是与系数矩阵的尺度化(scaling)有关。这两种条件数都有各自的应用范围。本文基于投影尺度引入齐次多项式特征值问题的条件数，其优点是与系数矩阵的尺度化无关，因而也许会有较广泛的应用。     

南海晚第四纪硅藻遗体及其地质意义：Ⅷ.几个问题的讨论

本文根据分析结果阐述了两种不同样品处理方法所获得的硅藻丰度和微型硅藻的人主要种类和相对百分含量。研究表明洗涤富集法所获得的硅藻丰度存在较大误差；微型硅藻的含量普遍介于１０％－２０％之间，最高可达４２．３１％，最低为５．６６％，并指出了在微型硅藻相对含量大于２０％的样品中 ，其主要种类为具槽直链藻和条纹小环藻等沿岸底格格不入中，并成群体产出。     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

Bowen数和东中国海Bowen数分布

本文讨论Bowen数的意义、功能和计算法。同时,依据多年水文气象实测资料作统计,计算出东中国海的Bo值。其结果绘制成1月至12月的月平均分布图,从而对本海域的Bo分布特点作详细分析介绍。     

玛纳斯河一级电站引水渠沿线提取地下水进行融冰的可行性研究

通过对向斜洼地含水层水文地质特征的勘探及研究,对研究区的水文地质边界条件进行了概化,并用水均衡法估算了地下水资源量,对拟建一级电站引水渠沿线提取地下水进行融冰的可行性进行了研究,并提出了取水方案.     

Simulation of marine boundary layer characteristics using a 1-D PBL model over the Bay of Bengal during BOBMEX-99

The characteristic features of the marine boundary layer (MBL) over the Bay of Bengal during the southwest monsoon and the factors influencing it are investigated. The Bay of Bengal and Monsoon Experiment (BOBMEX) carried out during July–August 1999 is the first observational experiment under the Indian Climate Research Programme (ICRP). A very high-resolution data in the vertical was obtained during this experiment, which was used to study the MBL characteristics off the east coast of India in the north and south Bay of Bengal. Spells of active and suppressed convection over the Bay were observed, of which, three representative convective episodes were considered for the study. For this purpose a one-dimensional multi-level PBL model with a TKE-ε closure scheme was used. The soundings, viz., the vertical profiles of temperature, humidity, zonal and meridional component of wind, obtained onboard ORV Sagar Kanya and from coastal stations along the east coast are used for the study. The temporal evolution of turbulent kinetic energy, marine boundary layer height (MBLH), sensible and latent heat fluxes and drag coefficient of momentum are simulated for different epochs of monsoon and monsoon depressions during BOBMEX-99.The model also generates the vertical profiles of potential temperature, specific humidity, zonal and meridional wind. These simulated values compared reasonably well with the observations available from BOBMEX.     

Magnetotelluric Investigation of the Hydrothermal System and Heat Source in the Muine-Toyoha Geothermal Area, Hokkaido, Japan

Abstract. Magnetotelluric (MT) surveys were carried out around the Muine volcano, Hokkaido, Japan, where it is expected that the heat and metal source forming the polymetallic Ag-Pb-Zn-Cu-In Toyoha deposit is present at depth. Measurements were performed at 20 sites, 18 of which were located along a WSW-ENE profile traversing the north ridge of Mt. Muine. A resistivity model obtained from 2D inversion of the MT data shows subsurface specific conductive and resistive features. Conductive layers are present at the surface of Mt. Muine. The low resistivity is probably due to the clay-rich rocks associated with the hydrothermal alteration. A high resistivity layer, which corresponds to the pre-Tertiary Usubetsu Formation, crops out east of Mt. Muine and dips westward. At the west foot of Mt. Muine, relatively high resistive layers are widely exposed. The resistivity increases with depth and exceeds 1000 ohm-m. This fact indicates that this region is not influenced by the recent hydrothermal activity. An extremely conductive zone about 3–6 km wide and 6–9 km thick exists at a depth of 2 km below Mt. Muine. This zone mostly corresponds to an elastic wave attenuation zone detected by a seismic survey. It is interpreted as a large hydrothermal reservoir or melted magma, which is a heat source of the hydrothermal system in this area.