Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

An HLLC Riemann solver for relativistic flows – II. Magnetohydrodynamics     

A. Mignone  G. Bodo 《Monthly notices of the Royal Astronomical Society》2006,368(3):1040-1054

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The Harten–Lax–van Leer contact wave (HLLC) solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave.
Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the constrained transport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one- and two-dimensional numerical test problems.  相似文献   

Angular-diameter distance estimates from the Sunyaev–Zeldovich effect in hydrodynamical cluster simulations     

S. Ameglio  S. Borgani  A. Diaferio  K. Dolag 《Monthly notices of the Royal Astronomical Society》2006,369(3):1459-1468

  相似文献   

Radiative Transfer in Inhomogeneous Atmospheres. III     

A. G. Nikoghossian 《Astrophysics》2004,47(3):412-421

The approach proposed in the previous parts of this series of papers is used to solve the radiative transfer problem in scattering and absorbing multicomponent atmospheres. Linear recurrence relations are obtained for both the reflectance and transmittance of these kinds of atmospheres, as well as for the emerging intensities when the atmosphere contains energy sources. Spectral line formation in a one-dimensional inhomogeneous atmosphere is examined as an illustration of the possibility of generalizing our approach to the matrix case. It is shown that, in this case as well, the question reduces to solving an initial value problem for linear differential equations. Some numerical calculations are presented.  相似文献   

应用ATA的对数周期馈源方案于FAST的可能性探讨     

吴盛殷  熊继衮  吕善伟  张江林  徐祥  南仁东  金乘进 《天文学进展》2004,22(2):179-187

介绍应用于美国凤凰计划的Allen望远镜阵(ATA)的对数周期馈源(LPA)概况，初步给出此种馈源的基本参数及性能，仿真估算了其驻波比和方向图，并结合正在进行中的我国500m孔径球面射电望远镜(FAST)计划，探讨了应用这种宽带馈源的可能性及其限制。  相似文献   

A qualitative and quantitative evaluation of experimental model catchment evolution     

G. R. Hancock  G. R. Willgoose 《水文研究》2003,17(12):2347-2363

Due to the geological time scales required for observation of catchment evolution, surrogates or analogues of field data are necessary to understand long‐term processes. To investigate long‐term catchment behaviour, two experimental model catchments that developed without rigid boundaries under controlled conditions are examined and a qualitative and quantitative analysis of their evolution is presented. Qualitatively, the experimental catchments have the visual appearance of field scale data. Observation demonstrates that changes in catchment shape and network form are conservative. Quantitative analysis suggests that the catchments reach an equilibrium form while a reduction in the channel network occurs. While the catchments are laboratory scale models, the results provide insights into field scale behaviour. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

Size and shape evolution of pores in a viscoplastic matrix under compression     

M. Berli  M. L. Accorsi  D. Or 《国际地质力学数值与分析法杂志》2006,30(12):1259-1281

The frequent use of soils and earth materials for hydraulic capping and for geo‐environmental waste containment motivated our interest in detailed modelling of changes in size and shape of macro‐pores to establish links between soil mechanical behaviour and concurrent changes in hydraulic and transport properties. The objective of this study was to use finite element analysis (FEA) to test and extend previous analytical solutions proposed by the authors describing deformation of a single macro‐pore embedded in linear viscoplastic soil material subjected to anisotropic remote stress. The FEA enables to consider more complex pore geometries and provides a detailed picture of matrix yield behaviour to explain shortcomings of approximate analytical solutions. Finite element and analytical calculations agreed very well for linear viscous as well as for viscoplastic materials, only limited for the case of isotropic remote stress due to the simplifications of the analytical model related to patterns and onset of matrix‐yielding behaviour. FEA calculations were compared with experimental data obtained from a compaction experiment in which pore deformation within a uniform modelling clay sample was monitored using CAT scanning. FEA predictions based on independently measured material properties and initial pore geometry provided an excellent match with experimentally determined evolution of pore size and shape hence lending credence to the potential use of FEA for more complex pore geometries and eventually connect macro‐pore deformation with hydraulic properties. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

High-resolution calculations of merging neutron stars – II. Neutrino emission     

S. Rosswog  M. Liebendörfer 《Monthly notices of the Royal Astronomical Society》2003,342(3):673-689

The remnant resulting from the merger of two neutron stars produces neutrinos in copious amounts. In this paper we present the neutrino emission results obtained via Newtonian, high-resolution simulations of the coalescence event. These simulations use three-dimensional smoothed particle hydrodynamics together with a nuclear, temperature-dependent equation of state and a multiflavour neutrino leakage scheme. We present the details of our scheme, discuss the neutrino emission results from a neutron star coalescence and compare them with the core-collapse supernova case where neutrino emission has been studied for several decades. The average neutrino energies are similar to those in the supernova case, but contrary to the latter, the luminosities are dominated by electron-type antineutrinos that are produced in the hot, neutron-rich, thick disc of the merger remnant. The cooler parts of this disc contain substantial fractions of heavy nuclei, which, however, do not influence the overall neutrino emission results significantly. Our total neutrino luminosities from the merger event are considerably lower than those found in previous investigations. This imposes constraints on the ability of neutron star mergers to produce a gamma-ray burst via neutrino annihilation. The neutrinos are emitted preferentially along the initial binary rotation axis, an event seen 'pole-on' would appear much brighter in neutrinos than a similar event seen 'edge-on'.  相似文献   

The efficiency of the spiral-in of a black hole to the Galactic Centre     

Piero F. Spinnato  Michael Fellhauer  Simon F. Portegies Zwart 《Monthly notices of the Royal Astronomical Society》2003,344(1):22-32

We study the efficiency at which a black hole or dense star cluster spirals in to the Galactic Centre. This process takes place on a dynamical friction time-scale, which depends on the value of the Coulomb logarithm (ln Λ). We determine the accurate value of this parameter using the direct N -body method, a tree algorithm and a particle-mesh technique with up to two million plus one particles. The three different techniques are in excellent agreement. Our measurement for the Coulomb logarithm appears to be independent of the number of particles. We conclude that  ln Λ= 6.6 ± 0.6  for a massive point particle in the inner few parsec of the Galactic bulge. For an extended object, such as a dense star cluster, ln Λ is smaller, with a value of the logarithm argument Λ inversely proportional to the object size.  相似文献   

BACKGROUND STUDY ON NON-RESIDUAL PHASE OF ELEMENTS IN CHANGJIANG AND HUANGHE ESTUARINE SEDIMENTS     

刘新威  吴景阳  蔡月琴  林荣根  韩桂荣 《中国海洋湖沼学报》1996,14(3):282-290

Analysis of Zn, Cu, Pb, Co, Cr, Li, Ni, K, Al, Fe extracted by 1 mol/L HCl or 0.5 mol/LHCl/H_2O_2, showed concentrations of Zn, Cu, Pb, Co, Cr, Fe, Ni were significantly correlated with Li, Al,K, and clay. Two methods are used to indicate the background value of the non-residual phase of elementsin sediments, and are the same as the methods used to indicate the background value of totalconcentrations in sediments. The first method uses correlograms and regression equations,the second usesthe mean element concentrations normalized with grain size. Li, Al, K can be used as reference elements to determine the background value of Zn, Cu, Pb, Co,Cr, Ni, Fe, while the clay concentration's correlation with some extractable concentrations can be used tocalculate the background value of the non-residual phase of elements as a percentage of clay concentrationin the sediments. Based on this study, the concept of using the background value of the non-residualphase of elements to compare the pollution level in differ  相似文献