Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.     

Effect of bed surface roughness on longitudinal dispersion in artificial open channels

In several empirical and modelling studies on river hydraulics, dispersion was negatively correlated to surface roughness. In this study, it was aimed to investigate the influence of surface roughness on longitudinal dispersion under controlled conditions. In artificial flow channels with a length of 104 m, tracer experiments with variations in channel bed material were performed. By use of measured tracer breakthrough curves, average flow velocity, mean longitudinal dispersion, and mean longitudinal dispersivity were calculated. Longitudinal dispersion coefficients ranged from 0·018 m² s^?1 in channels with smooth bed surface up to 0·209 m² s^?1 in channels with coarse gravel as bed material. Longitudinal dispersion was linearly related to mean flow velocity. Accordingly, longitudinal dispersivities ranged between 0·152 ± 0·017 m in channels with smooth bed surface and 0·584 ± 0·015 m in identical channels with a coarse gravel substrate. Grain size and surface roughness of the channel bed were found to correlate positively to longitudinal dispersion. This finding contradicts several existing relations between surface roughness and dispersion. Future studies should include further variation in surface roughness to derive a better‐founded empirical equation forecasting longitudinal dispersion from surface roughness. Copyright © 2011 John Wiley & Sons, Ltd.     

分数阶混合随机泛函微分方程的能控性

主要研究了分数阶混合随机泛函微分方程的能控性.在无限维空间下,假设所考虑方程线性部分生成半群不是紧的,使用非紧性测度技术和Mönch不动点定理,给出了方程能控性充分条件,并通过一个例子说明了结论的有效性.     

Generation of oceanic internal gravity waves by a cyclonic surface stress disturbance

Atmospheric cyclones with strong winds significantly impact ocean circulation, regional sea surface temperature, and deep water formation across the global oceans. Thus they are expected to play a key role in a variety of energy transport mechanisms. Even though wind-generated internal gravity waves are thought to contribute significantly to the energy balance of the deep ocean, their excitation mechanisms are only partly understood.The present study investigates the generation of internal gravity waves during a geostrophic adjustment process in a Boussinesq model with axisymmetric geometry. The atmospheric disturbance is set by an idealized pulse of cyclonic wind stress with a Rankine vortex structure. Strength, radius and duration of the forcing are varied. The effect upon wave generation of stratification with variable mixed-layer depth is also examined.Results indicate that internal gravity waves are generated after approximately one inertial period. The outward radial energy flux is dominated by waves having structure close to vertical mode-1 and with frequency close to the inertial frequency. Less energetic higher mode waves are observed to be generated close to the sea floor underneath the storm. The total radiated energy corresponds to approximately 0.02% of the wind input. Deeper mixed-layer conditions as well as weaker stratification reduce this fraction.The low energy transfer rates suggest that other processes that drive vertical motion like surface heat fluxes, turbulent motion, mixed region collapse and storm translation are essential for significant energy extraction by internal gravity waves to occur.     

One‐dimensional dynamics of saturated incompressible porous media: analytical solutions and influence of inertia terms

The complexity of formulations for the hydromechanical coupled mechanics of porous media is typically minimised by simplifying assumptions such as neglecting the effect of inertia terms. For example, three formulations commonly employed to model practical problems are classified as fully dynamic, simplified dynamic and quasi‐static. Thus, depending on the porous media conditions, each formulation will have advantages and limitations. This paper presents a comprehensive analysis of these limitations when solving one‐dimensional fully saturated porous media problems in addition to a new solution that considers a more general loading situation. A phase diagram is developed to assist on the selection of which formulation is more appropriate and convenient regarding particular cases of porosity and hydraulic conductivity values. Non‐dimensional formulations are proposed to achieve this goal. Results using the analytical solutions are compared against numerical values obtained with the finite element method, and the effect of porosity is investigated. Copyright © 2016 John Wiley & Sons, Ltd.     

A survey of numerical methods for stochastic differential equations

The development of numerical methods for stochastic differential equations has intensified over the past decade. The earliest methods were usually heuristic adaptations of deterministic methods, but were found to have limited accuracy regardless of the order of the original scheme. A stochastic counterpart of the Taylor formula now provides a framework for the systematic investigation of numerical methods for stochastic differential equations. It suggests numerical schemes, which involve multiple stochastic integrals, of higher order of convergence. We shall survey the literature on these and on the earlier schemes in this paper. Our discussion will focus on diffusion processes, but we shall also indicate the extensions needed to handle processes with jump components. In particular, we shall classify the schemes according to strong or weak convergence criteria, depending on whether the approximation of the sample paths or of the probability distribution is of main interest.     

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

无结构网格上二维浅水流动的数值模拟

提出在无结构网格上建立有限体积高性能格式族的统一框架,通过引入跨单元界面法向数值通量的逆风分解,将一维Osher、TVD两种通量分裂格式自然地推广至二维浅水方程组。给出了各种情况下浅水方程组有限体积法边界处理的计算公式。最后,利用该格式对穹包溢流、陆地动边界和河口潮流三个问题进行计算。结果表明,这类守恒格式具有高精度、无振荡性以及处理复杂流态过渡、自动捕俘间断和模拟陆地动边界的功能。     

IMPLICIT NONLINEAR NORMAL MODE INITIALIZATION:A MULTIGRID APPROACH

The implicit nonlinear normal mode initialization (INMI) is applied to a tropical limited area shallow water modelin spherical coordinates.The boundary condition for the INMI scheme is based on the boundary formulation of themodel.The INMI scheme is found to be very efficient in suppressing spurious gravity wave oscillation and providing awell balanced initial data set for the model.The INMI scheme involves solving a number of elliptic type equations withvarying complexity.and hence an efficient numerical technique is required for solving such equations.In order to makeINMI computationally more attractive,we are employing the multigrid method for solving all the elliptic type equationsin the INMI scheme.The numerical procedures for the development of such multigrid solvers are briefly described.     

The SOMS model and its application to Lake Neuchâtel

A three-dimensional numerical circulation model (SOMS) based on primitive equations is described. The algorithm, by which Coriolis and vertical diffusion terms are treated implicitly while mass is still conserved exactly (algebraically), is discussed in detail. The model is applied to Lake Neuchâtel (Switzerland), to determine the general circulation under influence of the most prevailing wind.