印度与欧亚大陆碰撞构造变形数值分析

在全面考虑了大陆构造变形三重非线性特征的基础上，用动态有限元法研究了藏东对印度与欧亚大陆碰撞作用的大变形响应。数值模拟结果显示了明显的纵向收缩变形，变形主轴多具有顺时针转动的趋势，但横向位移在总体上并不显著。其最大值正位于红河和鲜水河断裂之间，这一结果论证了印度板块对欧亚大陆的碰撞可能是川滇菱形块体向南东运动的动力来源。红河、鲜水河、昆仑及阿尔金活断层在碰撞作用下发生了顺时针的旋转变形；而且，计算结果模拟出了红河断裂由左旋变为右旋，另３条断裂始终左旋的运动特征     

Stochastic analysis of immiscible displacement in porous media

The interface of two immiscible fluids flowing in porous media may behave in an unstable fashion. This instability is governed by the pore distribution, differential viscosity and interface tension between the two immiscible fluids. This study investigates the factors that control the interface instability at the wetting front. The development of the flow equation is based on the mass balance principle, with boundary conditions such as the velocity continuity and capillary pressure balance at the interface. By assuming that the two-phase fluids in porous media are saturated, a covariance function of the wetting front position is derived by stochastic theory. According to those results, the unstable interface between two immiscible fluids is governed by the fluid velocity and properties such as viscosity and density. The fluid properties that affect the interface instability are expressed as dimensionless parameters, mobility ratio, capillary number and Bond number. If the fluid flow is driven by gravitational force, whether the interface undergoes upward displacement or downward displacement, the variance of the unstable interface decreases with an increasing mobility ratio, increases with increasing capillary number, and decreases with increasing Bond number. For a circumstance in which fluid flow is horizontal, our results demonstrate that the capillary number does not influence the generation of the unstable interface.     

Shaking table tests on the lateral response of a pile buried in liquefied sand

This paper presents an experimental study on the lateral resistance of a pile subjected to liquefaction-induced lateral flow. To observe the soil surrounding the pile during liquefaction, it was modeled as a buried cylinder that corresponded to a sectional model of the prototype pile at a certain depth in the subsoil. In order to create a realistic stress condition in the model ground, the model was prepared in a sealed container and the overburden pressure was applied to the ground surface by a rubber pressure bag. The model pile was actuated back and forth through rods attached on each side by an electro-hydraulic actuator.This paper focuses on observing the deformation of the liquefied soil surrounding the pile when a large relative displacement between the pile and the soil is induced. The loading rate effect on the lateral resistance of the pile in the liquefied sand and the influence of the relative density are also investigated.Test results show that a larger resistance is mobilized as the loading rate becomes higher. When the loading rate is higher, the cylinder displacement required for the lateral resistance becomes smaller. It has been also observed that as the relative density of the soil increases, dilatancy of the soil in front of the pile also increases.     

The structure of the coastal density front at the outflow of Long Island Sound during spring 2002

South of the eastern end of Long Island (Montauk Point) along the Eastern U.S. coast, a coastal density front forms between the buoyant outflow plume of the Long Island Sound (LIS) and the denser shelf waters offshore. During a 2-day cruise in April 2002, measurements of the density and velocity structure of this front were obtained from high-resolution CTD and ADCP data. Transects show the front intersecting the bottom inshore of the 30 m isobath and shoaling offshore. Variability in the location of the front is small offshore of the 40 m isobath, yet tidal excursions of the front along the bottom are significant (5 km) inshore of this depth.The frontal structure of the LIS plume was similar to observations of bottom-trapped coastal density fronts and shelf break fronts. A coastal jet in the along front direction was the main feature of the mean velocity field and was found to be in thermal wind balance with the mean density field. Stronger than expected offshore velocities near the surface, most likely a result of wind forcing, were the only exception to these similarities. In addition, analysis of temperature and salinity gradients along isopycnals gives evidence of secondary cross-frontal circulation and detachment of the bottom boundary layer. Characteristics of the LIS plume are used to evaluate recent analytical models of bottom-trapped coastal density fronts and bottom-advected plume theory, finding good agreement.     

Intermittency of the solar wind density near the interplanetary shock

The properties of turbulent fluctuations of the solar wind plasma near the interplanetary shock observed at September 12, 2014 by the BMSW instrument are considered. The spectra of the density fluctuations in the solar wind and their statistical characteristics up-and downstream of the shock front are analyzed. They are compared with each other and with characteristics corresponding to different turbulence models. It is shown that the spectral and statistical characteristics of the density fluctuations in the solar wind conserve their basic properties after the arrival of an interplanetary shock. Intermittency is observed both before and after the front, but its level increases on average in the second case. In both regions, the scaling of the structure functions of the density fluctuations in the solar wind differ from the scaling of the classical Kolmogorov model and can be described by the log-Poisson turbulence model. Parameterization of the scaling of the structure functions revealed the presence of filamentary structures in the solar wind plasma, which provide the density intermittency in the studied space regions.     

Nonlinear waves on a coupled density front

Abstract

It is shown that the inclusion of the nonlinear terms in the equations of motion of a coupled density front of zero potential vorticity results in wave solutions which merely propagate with time. The linear theory, on the other hand, predicts an exponential temporal growth. The nonlinear equation admits steady solutions representing standing waves whereas if the nonlinear terms are omitted no steady solutions exist. The general initial value problem is difficult to solve numerically since the linear problem is ill posed.

In addition we prove that the general similarity solution of the nonlinear equation tends to zero for large times, at any point in space, regardless of the initial condition.     

Directions of preferential flow in a hillslope soil, 1. Quasi‐steady flow

Preferred infiltration is mainly perceived as vertically down whereas subsurface storm flow is thought to occur parallel to slopes. The transition from vertical to lateral flow in a layered hillslope soil is the focus of the contribution. Transient flow is assumed to move as a wetting front. Three time‐domain reflectometry (TDR) wave‐guides, each 0·15 m long, were mounted in the shape of a truncated tetrahedron with its peak pointing down. Each wave‐guide focuses the front velocity along its axis. The three front‐velocity vectors are decomposed into their x, y and z components, which are then assembled to the resultant velocity vector. The volume density flux of preferred flow is the product of the front velocity and the mobile water content. The latter is the amplitude of transient soil moisture measured with each wave‐guide. The resultant vector of the volume flux density is computed similarly to the velocity vector. The experimental approach allows for the rapid assessment of transient flows without relying on the variation of water potentials. The experiments indicate that the directions of the resultant vectors of velocity and volume flux density can be estimated if the moisture variations of the three TDR wave‐guides are strongly correlated during the passing of the wetting front. Copyright © 2004 John Wiley & Sons, Ltd.     

Deplacement miscible avec contrastes de densite et de viscosite en milieu poreux. Identification des parametres de deplacement et stabilite en regime de dispersion mecanique

The behaviour and estimation of pollutant transport with water in an aquifer are based on the determination of displacement parameters (velocity field, dispersion tensor, kinematic porosity, etc.). We consider the one-dimensional and isothermal displacement of two miscible fluids through an isotropic and homogeneous porous medium in mechanical dispersion. An in situ measurement technique coupled with an impulse response method is used to verify the theoretical solution in the error function (trace case). With viscosity and density gradients, experimental and numerical (application of instability criterion) results are in good agreement. For a quasi-stable displacement, we have shown an asymptotic dispersion regime. The localization of a pollutant front in the studied configuration, is discussed.

Abstract

Le transfert d'un polluant par l'eau dans un milieu naturel repose sur la détermination des paramètres de déplacement (champ des vitesses tenseur de dispersion, porosité cinematique, ...). On considère un écoulement unidimensionnel et isotherme de deux fluides miscibles dans un milieu poreux homogène et isotrope en régime de dispersion mécanique. Une technique de mesure in situ couplée avec une méthode de réponse impulsionnelle a permis de vérifier la solution analytique en fonction erreur (cas du traceur). Dans le cas d'un déplacement avec contrastes de densité et de viscosité, les résultats expérimentaux et numériques (application d'un critère d'instabilité) concordent. Pour un déplacement quasi-stable, nous avons mis en évidence un régime de dispersion asymptotique. Nous analysons, ensuite, la localisation d'un front de pollution dans la configuration étudiée.     

Space-time model for migration of weak earthquakes along the northern boundary of the Amurian microplate

In this paper, we aimed to investigate the statistical distributions of shallow earthquakes with 2 ≤ М ≤ 4, located in 13 rectangular areas (clusters) bounded by 120°E and 144°E along the northern boundary of the Amurian microplate. As a result of our study, the displacement of seismicity maxima has been determined and three recurrent spatial cycles have been observed. The clusters with similar distribution of earthquakes are suggested to alternate being equally spaced at 7.26° (360–420 km). A comparison of investigation results on the structure of seismicity in various segments of the Amurian microplate reveals the identity between the alternation pattern observed for meridional zones of large earthquakes and a distinguished spatial period. The displacement vector for seismicity in the annual cycles is determined, and the correspondence between its E-W direction and the displacement of the fronts of large earthquakes is established. The elaborated model of seismic and deformation processes is considered, in which subsequent activation of clusters of weak earthquakes (2 ≤ М ≤ 4), tending to extend from the Japanese-Sakhalin island arc to the eastern closure of the Baikal rift zone, is initiated by the displacement of the strain wave front.     

基于数值流形法的反倾层状岩质边坡倾倒破坏分析

边坡破坏是累积性过程,从变形到破坏的过程中会产生永久位移,如果永久位移过大,极有可能产生滑坡.因此根据不同工况下采集到的位移数据,分析地震作用下反倾层状岩质边坡在不同内摩擦角下的破坏特征.利用二维数值流形法(NMM),以青藏高原金沙江流域西藏昌都地区芒康县索多西乡贡扎倾倒滑坡为研究对象,依据实地考察数据及室内力学试验得...     

Hydraulic displacement of dense nonaqueous phase liquids for source zone stabilization

Hydraulic displacement is a mass removal technology suitable for stabilization of a dense, nonaqueous phase liquid (DNAPL) source zone, where stabilization is defined as reducing DNAPL saturations and reducing the risk of future pool mobilization. High resolution three-dimensional multiphase flow simulations incorporating a spatially correlated, heterogeneous porous medium illustrate that hydraulic displacement results in an increase in the amount of residual DNAPL present, which in turn results in increased solute concentrations in groundwater, an increase in the rate of DNAPL dissolution, and an increase in the solute mass flux. A higher percentage of DNAPL recovery is associated with higher initial DNAPL release volumes, lower density DNAPLs, more heterogeneous porous media, and increased drawdown of groundwater at extraction wells. The fact that higher rates of recovery are associated with more heterogeneous porous media stems from the fact that larger contrasts in permeability provide for a higher proportion of capillary barriers upon which DNAPL pooling and lateral migration can occur. Across all scenarios evaluated in this study, the ganglia-to-pool (GTP) ratio generally increased from approximately 0.1 to between approximately 0.3 and 0.7 depending on the type of DNAPL, the degree of heterogeneity, and the imposed hydraulic gradient. The volume of DNAPL recovered as a result of implementing hydraulic displacement ranged from between 9.4% and 45.2% of the initial release volume, with the largest percentage recovery associated with 1,1,1 trichloroethane, the least dense of the three DNAPLs considered.     

Effect of coastal boundary resolution and mixing upon internal wave generation and propagation in coastal regions

A non-linear two-dimensional vertically stratified cross-sectional model of a constant depth basin without rotation is used to investigate the influence of vertical and horizontal diffusion upon the wind-driven circulation in the basin and the associated temperature field. The influence of horizontal grid resolution, in particular the application of an irregular grid with high resolution in the coastal boundary layer is examined. The calculations show that the initial response to a wind impulse is downwelling at the downwind end of the basin with upwelling and convective mixing at the opposite end. Results from a two-layer analytical model show that the initial response is the excitation of an infinite number of internal seiche modes in order to represent the initial response which is confined to a narrow near coastal region. As time progresses, at the downwind end of the basin a density front propagates away from the boundary, with the intensity of its horizontal gradient and associated vertical velocity determined by both horizontal and vertical viscosity values. Calculations demonstrate the importance of high horizontal grid resolution in resolving this density gradient together with an accurate density advection scheme. The application of an irregular grid in the horizontal with high grid resolution in the nearshore region enables the initial response to be accurately reproduced although physically unrealistic short waves appear as the frontal region propagates onto the coarser grid. Parameterization of horizontal viscosity using a Smagorinsky-type formulation acts as a selective grid size-dependent filter, and removes the short-wave problem although enhanced smoothing can occur if the scaling coefficient in the formulation is too large. Calculations clearly show the advantages of using an irregular grid but also the importance of using a grid size-dependent filter to avoid numerical problems.     

General spectral theory for the onset of instabilities in displacement processes in porous media

Summary The problem of penetration of a fluid into a porous medium containing a more viscous liquid is investigated. It is known that the displacement front may become unstable in this case as it may break up into «fingers». The problem of inception of fingers has been treated previously in the literature by describing the displacement front in terms of its Fourier transform. In the present paper, we generalize earlier procedures by making allowance for an arbitrary elemental growth law. Furthermore, we assume that the phenomenon of fingering is not solely governed by the prevailing flow potentials, but also by the spectrum of heterogeneities in the porous medium. This is achieved by introducing a constant characteristic of the frequency of the heterogeneities in the porous medium. It then turns out that the maximum rate of growth as a function of wave length is considerably shifted from that predicted in the literature. At the same time it is also shown that the difficulty encountered by other workers which consists of small wave lengths growing at an infinitely high rate, is being avoided.Paper presented orally at the 1960 Meeting of the A.I.M.E. in Denver (as paper number 1551-G of the program).     

Gravity current down a steeply inclined slope in a rotating fluid

The sinking of dense water down a steep continental slope is studied using laboratory experiments, theoretical analysis and numerical simulation. The experiments were made in a rotating tank containing a solid cone mounted on the tank floor and originally filled with water of constant density. A bottom gravity current was produced by injecting more dense coloured water at the top of the cone. The dense water plume propagated from the source down the inclined cone wall and formed a bottom front separating the dense and light fluids. The location of the bottom front was measured as a function of time for various experimental parameters. In the majority of runs a stable axisymmetric flow was observed. In certain experiments, the bottom layer became unstable and was broken into a system of frontal waves which propagated down the slope. The fluid dynamics theory was developed for a strongly non-linear gravity current forming a near-bottom density front. The theory takes into account both bottom and interfacial friction as well as deviation of pressure from the hydrostatic formula in the case of noticeable vertical velocities. Analytical and numerical solutions were found for the initial (t < 1/f), intermediate (t 1/f), and main (t 1/f) stages, where f is the Coriolis parameter. The model results show that during the initial stage non-linear inertial oscillations are developed. During the main stage, the gravity current is concentrated in the bottom layer which has a thickness of the order of the Ekman scale. The numerical solutions are close to the same analytical one. Stability analysis shows that the instability threshold depends mainly on the Froude number and does not depend on the Ekman number. The results of laboratory experiments confirm the similarity properties of the bottom front propagation and agree well with the theoretical predictions.     

Asymmetry of nonlinear interactions of solar MHD discontinuities with the bow shock

Oblique interaction between the solar fast shock wave, which is a typical nonstationary strong discontinuity in the interplanetary space, and the bow shock front upstream of an Earth-type planetary magnetosphere is studied. Attention has been paid to the qualitative and quantitative (with respect to the proton density distribution) dawn-dusk (or morning-evening) asymmetry of the discontinuities refracted into the magnetosheath, which originates in the ecliptic plane on different sides of the Sun-Earth line. The results under discussion have been corroborated experimentally by the gas-kinetic pattern of the bow-shock front and the WIND and ISEE 3 spacecraft measurements of the plasma density.     

遮帘式板桩码头地基地震液化破坏机理

遮帘式板桩码头作为一种新型的板桩结构型式,其抗震性能研究是设计建造过程中的重要环节。在FEM-FDM水土耦合计算的平台上引入循环弹塑性本构模型,借助FORTRAN编程软件形成饱和砂土动力液化分析的数值方法,可有效模拟饱和砂土在地震动力作用下的非线性及大变形特性,同时也可模拟砂土液化流动对遮帘桩和前墙的动土压力。研究表明:地震作用下可液化土层超孔隙水压力比增长并发生较大的水平流动变形,对前墙的水平破坏大于竖向破坏;前墙剪力最大值位于海床与前墙交界处;遮帘桩剪力最大值位移与前墙底平行的位置;后拉杆拉力逐渐变大,前拉杆拉力逐渐变小。通过对板桩码头地震液化灾害的分析,可为抗震和抗液化设计提供参考依据。     

A laboratory simulation of pyroclastic flows down slopes

Laboratory experiments are described which explore the dynamical consequences of buoyant convective upflow observed above hot pyroclastic flows. In nature, the convection is produced by the hot ash particles exchanging heat with air mixed into the front and top of the pyroclastic flow. This effect on the buoyancy due to the mixing of air and ash has been modelled in the laboratory using mixtures of methanol and ethylene glycol (MEG), which have a nonlinear density behaviour when mixed with water. Intermediate mixtures of these fluids can be denser than either initial component, and so the laboratory experiments were inverted models of the natural situation. We studied MEG flowing up under a sloping roof in a tank filled with water. The experiments were performed both in a narrow channel and on a laterally unconfined slope. The flow patterns were also compared with those of conventional gravity currents formed using fresh and salt water. The presence of the region of reversed buoyancy outside the layer flowing along the slope had two significant effects. First, it periodically protected the flow from direct mixing with the environment, resulting in pulses of relatively undiluted fluid moving out intermittently ahead of the main flow. Second, it produced a lateral inflow towards the axis of the current which kept the current confined to a narrow tongue, even on a wide slope.In pyroclastic flows the basal avalanche portion has a much larger density contrast with its surroundings than the laboratory flows. Calculations show that mixing of air into the dense part of a pyroclastic flow cannot generate a mixture that is buoyant in the atmosphere. However, the overlying dilute ash cloud can behave as a gravity current comparable in density contrast to the laboratory flows and can become buoyant, depending on the temperature and ash content. In the August 7th pyroclastic flow of Mount St. Helens, Hoblitt (1986) describes pulsations in the flow front, which are reminiscent of those observed in the experiments. As proposed by Hoblitt, the pulsations are caused by the ash cloud accelerating away from the front of the dense avalanche as a density current. The ash cloud then mixes with more air, becomes buoyant and lifts off the ground, allowing the avalanche to catch up with and move ahead of the cloud. The pulsing behaviour at the fronts of pyroclastic flows could account for the occurrence of cross-bedded layer 1 deposits which occur beneath layer 2 deposits in many sequences.     

Motion of wetting fronts moving into partially pre-wet soil

We study the motion of wetting fronts for vertical infiltration problems as modeled by Richards’ equation. Parlange and others have shown that wetting fronts in infiltration flows can be described by traveling wave solutions. If the soil layer is not initially dry, but has an initial distribution of water content then the motion of the wetting front will change due to the interaction of the infiltrating flow with the pre-existing soil conditions. Using traveling wave profiles, we construct simple approximate solutions of initial-boundary value problems for Richards’ equation that accurately describe the position and moisture distribution of the wetting front. We show that the influences of surface boundary conditions and initial conditions produce shifts to the position of the wetting front. The shifts can be calculated by examining the cumulative infiltration, and are validated numerically for several problems for Richards’ equation and the linear advection–diffusion equation.     

Tomographic measurements of pore filling at infiltration fronts

For certain porous media and initial conditions, constant flux infiltrations show a saturation profile which exhibits overshoot. This overshoot is the cause of gravity driven fingering, cannot be described by standard models of unsaturated flow, and is likely controlled by the exact nature of the pore filling at the initial front. Here we report synchrotron X-ray micro-tomography measurements of the porous medium and measure which pores are filled by water and air at the initial wetting front as a function of flux. We find that at high fluxes all the pores are filled with water; for intermediate fluxes, the pores along the edge of the column remain unsaturated; and for low fluxes the pores in the bulk of the experimental column remain unsaturated. This suggests that the unsaturated overshoot conditions observed at higher fluxes are primarily an edge effect of the column. The results can help delineate the correct continuum models that can capture overshoot and gravity driven fingering.