An improved SPH method for saturated soils and its application to investigate the mechanisms of embankment failure: Case of hydrostatic pore‐water pressure     

Ha H. Bui  Ryoichi Fukagawa 《国际地质力学数值与分析法杂志》2013,37(1):31-50

The method of smoothed particle hydrodynamics (SPH) has recently been applied to computational geomechanics and has been shown to be a powerful alternative to the standard numerical method, that is, the finite element method, for handling large deformation and post‐failure of geomaterials. However, very few studies apply the SPH method to model saturated or submerged soil problems. Our recent studies of this matter revealed that significant errors may be made if the gradient of the pore‐water pressure is handled using the standard SPH formulation. To overcome this problem and to enhance the SPH applications to computational geomechanics, this article proposes a general SPH formulation, which can be applied straightforwardly to dry and saturated soils. For simplicity, the current work assumes hydrostatic pore‐water pressure. It is shown that the proposed formulation can remove the numerical error mentioned earlier. Moreover, this formulation automatically satisfies the dynamic boundary conditions at a submerged ground surface, thereby saving computational cost. Discussions on the applications of the standard and new SPH formulations are also given through some numerical tests. Furthermore, techniques to obtain the correct SPH solution are also proposed and discussed throughout. As an application of the proposed method, the effect of the dilatancy angle on the failure mechanism of a two‐sided embankment subjected to a high groundwater table is presented and compared with that of other solutions. Finally, the proposed formulation can be considered a basic formulation for further developments of SPH for saturated soils. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Pressure and temperature dependence of cation distribution in Mg-Mn olivine     

Tadashi Akamatsu  Kiyoshi Fujino  Mineo Kumazawa  Akio Fujimura  Manabu Kato  Hiroshi Sawamoto  Takamitsu Yamanaka 《Physics and Chemistry of Minerals》1988,16(2):105-113

Synthetic (Mg_0.51, Mn_0.49)₂SiO₄ olivine samples are heat-treated at three different pressures; 0, 8 and 12 GPa, all at the same temperature (～500° C). X-ray structure analyses on these single crystals are made in order to see the pressure effect on cation distribution. The intersite distribution coefficient of Mg and Mn in M1 and M2 sites, K _D = (Mn/Mg) _M1/(Mn/Mg) _M2, of these samples are 0.192 (0 GPa), 0.246 (8 GPa) and 0.281 (12 GPa), indicating cationic disordering with pressure. The small differences of cell dimensions between these samples are determined by powder X-ray diffraction. Cell dimensions b and c decrease, whereas a increases with pressure of equilibration. Cell volume decreases with pressure as a result of a large contraction of the b cell dimension. The effect of pressure on the free energy of the cation exchange reaction is evaluated by the observed relation between the cell volume and the site occupancy numbers. The magnitude of the pressure effect on cation distribution is only a fifth of that predicted from the observed change in volume combined with thermodynamic theory. This phenomenon is attributed to nonideality in this solid solution, and nonideal parameters are required to describe cation distribution determined in the present and previous experiments. We use a five-parameter equation to specify the cationic equilibrium on the basic of thermodynamic theory. It includes one energy parameter of ideal mixing, two parameters for nonideal effects, one volume parameter, and one thermal parameter originated from the lattice vibrational energy. The present data combined with some of the existing data are used to determine the five parameters, and the cation distribution in Mg-Mn olivine is described as a function of temperature, pressure, and composition. The basic framework of describing the cationic behavior in olivine-type mineral is worked out, although the result is preliminary: each of the determined parameters is not accurate enough to enable us to make a reliable prediction.  相似文献   

Structural change of GeO2 under pressure     

S. Kawasaki  O. Ohtaka  T. Yamanaka 《Physics and Chemistry of Minerals》1994,20(8):531-535

Pressure-induced amorphization of α-quartz type GeO₂ was studied with a newly developed X-ray diffraction system which consists of a 4-circle goniometer and a curved position sensitive detector. Single-crystal diffraction was measured under pressurs up to 7.3 GPa at room temperature in order to investigate pretransitional phenomena. Diffraction intensity and line width of the diffraction profiles showed no remarkable change up to 5.9 GPa. However, no sharp diffraction line was observed at pressures over 6.5 GPa. The bulk modulus at 0.1 MPa and its pressure derivative of α-quartz type GeO₂ were determined to be K _T =32.8(3.3) GPa and K′ _T =6.0(2.0), respectively. In situ microscopic observations of the amorphization transformation was also performed. The large volume change due to amorphization was observed and estimated to be about 10%.  相似文献   

The antiferromagnetic-paramagnetic transition in αFe2O3 in the single prism approximation of the cluster variation method     

Benjamin Burton  Ryoichi Kikuchi 《Physics and Chemistry of Minerals》1984,11(3):125-131

The single prism approximation SPA of the cluster variation method has been used to model the antiferromagnetic-paramagnetic transition in hematite. This calculation yields insight into the accuracy of the SPA and other approximate methods for modeling order-disorder phenomena. Published values of the magnetic coupling constants were used to calculate the Néel temperature, sublattice magnetization, and magnetic specific heat. The calculated Néel temperature is found to be 1.21 times the observed value, as compared to 1.36 times observed for a mean field theory approximation, an improvement that reflects the superior treatment of configurational entropy in the SPA. Qualitative to semiquantitative agreement is obtained between observed and calculated values for sublattice magnetization and magnetic specific heat; however we find that previously published values for the magnetic specific heat are too large by a factor of two.  相似文献   

Electric conductivity of Fe2SiO4–Fe3O4 spinel solid solutions     

T. Yamanaka  H. Shimazu  K. Ota 《Physics and Chemistry of Minerals》2001,28(2):110-118

 The spinel solid solution was found to exist in the whole range between Fe₃O₄ and γ-Fe₂SiO₄ at over 10 GPa. The resistivity of Fe_{3− x} Si _x O₄ (0.0<x<0.288) was measured in the temperature range of 80∼300 K by the AC impedance method. Electron hopping between Fe³⁺ and Fe²⁺ in the octahedral site of iron-rich phases gives a large electric conductivity at room temperature. The activation energy of the electron hopping becomes larger with increasing γ-Fe₂SiO₄ component. A nonlinear change in electric conductivity is not simply caused by the statistical probability of Fe³⁺–Fe²⁺ electron hopping with increasing the total Si content. This is probably because a large number of Si⁴⁺ ions occupies the octahedral site and the adjacent Fe²⁺ keeping the local electric neutrality around Si⁴⁺ makes a cluster, which generates a local deformation by Si substitution. The temperature dependence of the conductivity of solid solutions indicates the Verwey transition temperature, which decreases from 124(±2) K at x=0 (Fe₃O₄) to 102(±5) K at x=0.288, and the electric conductivity gap at the transition temperature decreases with Si⁴⁺ substitution. Received: 15 March 2000 / Accepted: 4 September 2000  相似文献   

Spatial distribution of deuterium in atmospheric water vapor: Diagnosing sources and the mixing of atmospheric moisture   总被引：1，自引：0，他引：1  

Tsutomu Yamanaka  Ryosuke Shimizu 《Geochimica et cosmochimica acta》2007,71(13):3162-3169

We measured the stable isotopic composition of hydrogen (δD) within atmospheric water vapor collected simultaneously at six sites in the vicinity of a lake (Lake Kasumigaura, Eastern Japan) to determine its spatial distribution characteristics and thereby diagnose sources and mixing of atmospheric moisture. The measured spatial distribution of δD showed no relation to distance from the lake, although it showed a correlation with the distribution of the water-vapor mixing ratio Q. For two of the three sampling days, we found a simple two-component (i.e., water vapor transpiring from local land surfaces and pre-existing vapor in the background atmosphere) mixing line in a Keeling plot (i.e., δ − 1/Q diagram). On a third day, however, contributions from lake evaporation were detected in addition to the above components. On this day, lake-derived vapor accounted for approximately 10-20% of atmospheric water vapor at the sites located leeward of the lake. The observed differences in mixing patterns among sampling days can be explained by a simple atmospheric moisture budget. Thus, it is likely that simultaneous isotopic measurements of atmospheric water vapor at multiple locations with aid of Keeling plot are capable of giving us useful information in diagnosing the sources and mixing pattern of the vapor.  相似文献   

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model     

Ha H. Bui  Ryoichi Fukagawa  Kazunari Sako  Shintaro Ohno 《国际地质力学数值与分析法杂志》2008,32(12):1537-1570

Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

Application of multi-isotope ratios to study the source and quality of urban groundwater in Metro Manila,Philippines   总被引：1，自引：0，他引：1  

Takahiro Hosono  Fernando Siringan  Tsutomu Yamanaka  Yu Umezawa  Shin-ichi Onodera  Takanori Nakano  Makoto Taniguchi 《Applied Geochemistry》2010

  相似文献   

Bayesian and Neural Network Approaches to Estimate Deep Temperature Distribution for Assessing a Supercritical Geothermal System: Evaluation Using a Numerical Model     

Ishitsuka  Kazuya  Kobayashi  Yosuke  Watanabe  Norihiro  Yamaya  Yusuke  Bjarkason  Elvar  Suzuki  Anna  Mogi  Toru  Asanuma  Hiroshi  Kajiwara  Tatsuya  Sugimoto  Takeshi  Saito  Ryoichi 《Natural Resources Research》2021,30(5):3289-3314

The temperature distribution at depth is a key variable when assessing the potential of a supercritical geothermal resource as well as a conventional geothermal resource. Data-driven estimation by a machine-learning approach is a promising way to estimate temperature distributions at depth in geothermal fields. In this study, we developed two methodologies—one based on Bayesian estimation and the other on neural networks—to estimate temperature distributions in geothermal fields. These methodologies can be used to supplement existing temperature logs, by estimating temperature distributions in unexplored regions of the subsurface, based on electrical resistivity data, observed geological/mineralogical boundaries, and microseismic observations. We evaluated the accuracy and characteristics of these methodologies using a numerical model of the Kakkonda geothermal field, Japan, where a temperature above 500 °C was observed below a depth of about 3.7 km. When using geological and geophysical knowledge as prior information for the machine learning methods, the results demonstrate that the approaches can provide subsurface temperature estimates that are consistent with the temperature distribution given by the numerical model. Using a numerical model as a benchmark helps to understand the characteristics of the machine learning approaches and may help to identify ways of improving these methods.

  相似文献   

Polymetallic Mineralization at the Shin-Ohtoyo Deposit, Harukayama District, Hokkaido, Japan     

Akira IMAI  Hiroharu MATSUEDA    Ryoichi YAMADA  Kenzo MASUTA 《Resource Geology》1999,49(2):75-88

Abstract: The Shin-Ohtoyo Cu–Au deposit is located in the Harukayama district, 20 km west of Sapporo, Hokkaido, Japan. Both acid-type disseminated and adularia–quartz–type vein Au mineralizations have been recognized within a small distance of less than 500 m in the district. Mineralogical characteristics of sulfide ores from the Shin-Ohtoyo deposit have been proved to be polymetallic. Ore minerals containing Sn, V, Bi and Te are recognized. Nine ore types are recognized in terms of characteristic mineral assemblage; (1) chalcedonic quartz veinlets in silicified zone around the deposit, (2) bismuthinite, emplectite, friedrichite and tetrahedrite, (3) an unnamed Cu–Sn–Fe–Zn sulfide, colusite-series minerals, stannoidite, emplectite and tetrahedrite, (4) bournonite, Se-bearing galena and tetrahedrite, (5) luzonite/famatinite and Ag-bearing tetrahedrite, (6) colusite-series minerals, emplectite, aikinite and tetrahedrite/goldfieldite, (7) luzonite/famatinite, colusite-series minerals, mawsonite and tetra–hedrite/goldfieldite, (8) enargite, luzonite/famatinite and tetrahedrite, and (9) colusite-series minerals and tetrahedrite. The first occurrence of friedrichite and stibiocolusite from Japan are reported. The chemical formula of the unnamed phase corresponds to Cu₆(Cu, Fe, Zn)Sn₃S₁₀. Sulfur isotopic ratios (δ³⁴S) of sulfides from the stockpile range from –0. 5% to +1. 9%, and those from drill cores recovered by Metal Mining Agency of Japan (MMAJ) vary from –2. 7% to +0. 8%. Sulfur isotopic ratio of barite in a cavity in the silicified tuff breccia collected from the stock pile yields +27. 1%, while that of barite collected from MMAJ core is +21. 7%. Sulfur isotopic thermometry applied for a pair of barite (+21. 7%) and associated pyrite (+1. 8%) indicates about 300°C. High–Te tetrahedrite composition from both the chalcedonic quartz vein in the silicified zone around the Shin-Ohtoyo deposit and the polymetallic sulfide ores from the adit of the deposit, suggests that the Au mineralization in the former is attributed to a hydrothermal system marginal to the polymetallic mineralization.  相似文献