Well-posedness analysis and numerical implementation of a linearized two-dimensional bottom sediment transport problem

A two-dimensional linearized model of coastal sediment transport due to the action of waves is studied. Up till now, one-dimensional sediment transport models have been used. The model under study makes allowance for complicated bottom relief, the porosity of the bottom sediment, the size and density of sediment particles, gravity, wave-generated shear stress, and other factors. For the corresponding initial–boundary value problem the uniqueness of a solution is proved, and an a priori estimate for the solution norm is obtained depending on integral estimates of the right-hand side, boundary conditions, and the norm of the initial condition. A conservative difference scheme with weights is constructed that approximates the continuous initial–boundary value problem. Sufficient conditions for the stability of the scheme, which impose constraints on its time step, are given. Numerical experiments for test problems of bottom sediment transport and bottom relief transformation are performed. The numerical results agree with actual physical experiments.     

The effects of bed form roughness on total suspended load via the Lattice Boltzmann Method

Bed forms in natural rivers and man-made channels provide the dominant contribution to overall flow resistance and hence significantly affect sediment transport rate. Many laboratory experiments and field observations have been conducted on bed forms, and it was found that theoretical flat-bed assumptions do not give the correct estimation for the total suspended load (TSL). In this study, we present a systematic numerical investigation of turbulent open-channel flows over bed forms using the Lattice Boltzmann Method (LBM). A static Smagorinsky model is incorporated into LBM to account for turbulence, and the dynamic interface between fluid and air is captured by a free-surface model. The time-averaged flow velocity, turbulence intensity and Reynolds shear stress in LBM simulations show an excellent agreement with the available experimental data. In addition, the coherent flow structures induced by the bed forms qualitatively agree with previous numerical results from Large Eddy Simulations based the Navier–Stokes equations. We then proceed to investigate the effects of bed form roughness, quantified by the total friction factor f_T, on sediment transport. It is found that the prediction of the TSL based on the theoretical flat-bed assumptions may lead to an overestimation of up to 30%, depending on the bed form roughness. In addition, the normalized TSL is linearly proportional to f_T and nearly inversely proportional to the ratio of downward settling velocity and upward turbulence induced diffusion. Our work proposes a general law linking these quantities to estimate the TSL, which has the potential for a more efficient and accurate engineering design of man-made channels and improved river management.     

The model of dry friction in the problem of the rolling of rigid bodies

The Contensou model of combined dry friction [1] is considered. The problem of integrating the shear stresses over the contact area is solved in terms of elementary functions, unlike the solution in [1], reduced to elliptic quadratures. The problem of the rolling of a homogeneous sphere over a plane with dry friction is investigated.     

Sediment Transport via Dam-Break Flows over Sloping Erodible Beds

When a semi-infinite body of homogeneous fluid initially at rest behind a vertical retaining wall is suddenly released by the removal of the barrier the resulting flow over either a horizontal or a sloping bed is referred to as a dam-break flow. When resistance to the flow is neglected the exact solution, in the case of a stable horizontal bed with or without "tail water," may be obtained on the basis of shallow-water theory via the method of characteristics and the results are well known. Discrepancies between these shallow-water based solutions and experiments have been partially accounted for by the introduction of flow resistance in the form of basal friction. This added friction significantly modifies the wave speed and flow profile near the head of the wave so that the simple exact solutions no longer apply and various asymptotic or numerical approaches must be implemented to solve these frictionally modified depth-averaged shallow-water equations. When the bed is no longer stable so that solid particles may be exchanged between the bed and the water column the dynamics of the flow becomes highly complex as the buoyancy forces vary in space and time according to the competing rates of erosion and deposition. It is our intention here to study dam-break flows over erodible sloping beds as agents of sediment transport taking into account basal friction as well as the effects of particle concentrations on flow dynamics including both erosion and deposition. We consider shallow flows over initially dry beds and investigate the effects of changes in the depositional and erosional models employed as well as in the nature of the drag acting on the flow. These models include effects hitherto neglected in such studies and offer insights into the transport of sediment in the worst case scenario of the complete and instantaneous collapse of a dam.     

Uniform analysis of a stabilized hybrid finite element method for Reissner-Mindlin plates

This paper presents a low order stabilized hybrid quadrilateral finite element method for ReissnerMindlin plates based on Hellinger-Reissner variational principle,which includes variables of displacements,shear stresses and bending moments.The approach uses continuous piecewise isoparametric bilinear interpolations for the approximations of the transverse displacement and rotation.The stabilization achieved by adding a stabilization term of least-squares to the original hybrid scheme,allows independent approximations of the stresses and moments.The stress approximation adopts a piecewise independent 4-parameter mode satisfying an accuracy-enhanced condition.The approximation of moments employs a piecewise-independent 5-parameter mode.This method can be viewed as a stabilized version of the hybrid finite element scheme proposed in [Carstensen C,Xie X,Yu G,et al.A priori and a posteriori analysis for a locking-free low order quadrilateral hybrid finite element for Reissner-Mindlin plates.Comput Methods Appl Mech Engrg,2011,200:1161-1175],where the approximations of stresses and moments are required to satisfy an equilibrium criterion.A priori error analysis shows that the method is uniform with respect to the plate thickness t.Numerical experiments confirm the theoretical results.     

Velocity and concentration profiles in uniform sediment-laden flow

This paper presents a theoretical model for computing the velocity and sediment concentration profiles in a uniform sediment-laden flow carrying all fine, medium and coarse sediments. The proposed model essentially includes the effect of sediment concentration in total turbulent shear stress and eddy diffusivity in addition to the modified mixing length derived by Umeyama and Gerritsen [J. Hydr. Engng., ASCE, 118 (2) (1992) 229–245] applied to Hunt’s diffusion equation. Numerical solution of coupled differential equations for velocity and sediment concentration is carried out. The theoretical results show quite good agreement with the experimental data.     

Application of the method of partitioning the state of stress to the analysis of thermoelastic shells : PMM vol. 40, n 6, 1976, pp. 1085–1092

By using an asymptotic approach [1], the method of partitioning the state of stress is extended to thermoelastic shells. It is examined in detail in [2] forun-heated shells subjected to the effect of external forces, and consists of representing the total state of stress of the shell as the sum of those simpler states of stress for each of which the simplest methods for their construction can be given.Partitioning of the state of stress was performed in [3] for shells with a constant temperature over the thickness. It was noted in [4] in an analysis of a circular cylindrical shell by bending theory that integrals extended over the whole middle surface, which describe the fundamental state of stress, and integrals which damp out with distance from the edges and represent edge effects are contained in the general solution. In a number of papers, [5] for example, partitioning is performed on the basis of graphic physical representations for simple examples of analyzing circular cylindrical shells.A general approach to the analysis of rigid thermoelastic shells by the partitioning method is described below.     

A Four-node Finite Element for Piezoelectric Shell Structures in Convective Co-ordinates

In this contribution, an isoparametric piezoelectric shell element is presented which is based on convective coordinates and which allows for the analysis of arbitrary shell geometries. A two-field variation formulation [1, 2] is used in which the displacements and the electric potentials serve as independent variables. Especially, for thin-walled structures under certain boundary conditions and load cases, the displacement based element tend to shear and membrane locking. In order to avoid this poor behaviour, the Assumed Natural Strain (ANS) method [3] is introduced into the piezoelectric shell element. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A shallow-water sedimentation model with friction and Coriolis: An existence theorem

We present an existence theorem of a two-dimensional sedimentation model coupling a shallow water system with a sediment transport equation. The shallow water system includes Coriolis and friction terms. A Galerkin method is used to obtain a finite-dimensional problem which is solved using a Brouwer fixed point theorem. We prove that the limits of the resulting solution sequences satisfy the model equations.     

On the numerical modelling of initiation of sediment transport using Smoothed Particle Hydrodynamics

Mobilization of solid particles at the interface between a porous and a free flow domain is a relevant subject in many fields of mechanical, civil and environmental engineering. One example is the initiation of sediment transport as it appears in river beds. To approximate this initiation state, various theoretical models exist. Common approaches use two-domain formulations as in [1] or one-domain formulations as in [6]. The named approaches were compared with Direct Numerical Simulations (DNS) using Smoothed Particle Hydrodynamics (SPH) in [3]. The results of these simulations showed that the theoretical models often underestimate the occurring velocities at the interface and therefore critical velocities to initialize the motion of single grains can be lower than predicted by theoretical approaches. In our numerical simulations, we study creeping flow in a free flow domain coupled to flow in a porous media applying various porous structures. To investigate velocities and shear stresses at the interface more intensively we then compare our numerical results to data from experiments that were performed on equivalent microstructures. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A coupled wave–current–sediment transport model for an estuarine system: Application to the Río de la Plata and Montevideo Bay

Modelling the hydrodynamics and fine-sediment dynamics in estuarine environments is important for coastal engineering design and environmental assessment. This paper presents the application of a coupled wave–current–sediment transport model to a complex estuarine system with high spatial resolution in a harbour area. The case study was the Río de la Plata, with a focus on the Montevideo Bay area. We used a bidimensional depth-averaged approach, and the sediment transport modelling focused on fine cohesive sediments. The model considers realistic forcings, allows for the simulation of complex geometries such as those present in harbour basins and is capable of providing long-term environmental simulations (on the order of several years) within reasonable computational times. The model results are in good agreement with the measured data and satisfactorily represent the main features of the flow and sediment dynamics of the Río de la Plata. The effect of the internal coupling on the hydrodynamic results is analysed, and the computational times with various coupling alternatives are discussed.The dynamics of fine sediment in Montevideo Bay were analysed based on the model results. The current-induced bottom shear stress results are relevant for representing the permanent suspended sediment concentrations, whereas the wave-induced bottom shear stress is fundamental for reproducing the main resuspension events during storm conditions. The suspended fine-sediment dynamics in the estuary are strongly controlled by the sediment exchange between the bed and water column, whereas inside Montevideo Bay, the dynamics are controlled mainly by advection of sediment originating from the nearby coastal area.     

Hysteretic Friction of a Sliding Rubber Element

Friction, especially friction of elastomers, can cause acoustic problems like noise, squeal and comfort drawbacks like vibrations and wear. Therefore, rubber friction affects the function of many products in technical applications, e.g. seals, belts and tires. It can be classified according to different physical phenomena like adhesion, hysteresis, cohesion and viscous friction, see [3]. The topic of this paper is hysteresis friction of rubber that is caused by the energy dissipation due to internal material damping during the process of deformation. The deformation itself occurs during the sliding of a rubber element across the micro‐scaled asperities of a rough surface. In this paper, the sliding process of a rubber element over real surfaces is simulated in time domain and compared to experiments. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Verification of an anisotropic Coulomb adhesion-friction law for contact surfaces with periodic structure

Sliding friction forces and so-called adhesion forces are the main mechanical characteristics to describe contact interaction. Both together are representing 2D surface constitutive laws in analogy to e.g. elasto-plasticity for 3D continua. The classical model to generalize the Coulomb friction law into anisotropic domains is to introduce an anisotropic friction tensor. Michalowski and Mroz in [1] proposed the structure of the friction tensor considering the sliding of a rigid block on an inclined surface. Zmitrowicz in [2] developed the theoretical basis for the structure of the friction tensor on symmetry groups for the tensor. The current contribution is aimed at verification of this modeling process based on a homogenization procedure for a very fine discretization representing the exact structure of the surface. The validation issue with realistic experiments given in [4] is discussed as well. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of Coulomb friction on the fiber/matrix debonding and matrix cracking in unidirectional fiber-reinforced brittle-matrix composites

A model for a macroscopic crack transverse to bridging fibers is developed based upon the Coulomb friction law, instead of the hypothesis of a constant frictional shear stress usually assumed in fiber/matrix debonding and matrix cracking analyses. The Lamé formulation, together with the Coulomb friction law, is adopted to determine the elastic states of fiber/matrix stress transfer through a frictionally constrained interface in the debonded region, and a modified shear lag model is used to evaluate the elastic responses in the bonded region. By treating the debonding process as a particular problem of crack propagation along the interface, the fracture mechanics approach is adopted to formulate a debonding criterion allowing one to determine the debonding length. By using the energy balance approach, the critical stress for propagating a semi-infinite fiber-bridged crack in a unidirectional fiber-reinforced composite is formulated in terms of friction coefficient and debonding toughness. The critical stress for matrix cracking and the corresponding stress distributions calculated by the present Coulomb friction model is compared with those predicted by the models of constant frictional shear stress. The effect of Poisson contraction caused by the stress re distribution between the fiber and matrix on the matrix cracking mechanics is shown and discussed in the present analysis. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 2, pp. 171–190, March–April, 2007.     

Introduction of non-uniformity through linearization of the system governing turbulent, mixing of a scalar quantity in a 2-d mixing layer

The introduction of mathematical non-uniformity in the formulation of the turbulent mixing of a scalar quantity (mass, temperature, etc.) for a 2-d, free shear flow using Goertler's [ZAMM 22 (1942) 244] perturbation argument is discussed. Approximate, i.e. thin shear layer self-similar forms for mass, momentum and the scalar quantity are derived, and then linearized using Goertler's method. Though successful for the mean velocity field, the regular expansion yields inconsistent solutions for the transport of a scalar. Sources of the non-uniformity are identified using appropriate numerical methods for both non-linear and linear formulations. A consistent result is obtained by rescaling the independent variable and equation system and identifying dominant behavior. The results of this corrected formulation are shown to be consistent with the relationships obtained by the author using an approximate matched asymptotic expansion procedure.     

Asymptotic study of turbulent Poiseuille flow with wall transpiration

An incompressible, pressure–driven, fully developed turbulent flow between two parallel walls, with an extra constant transverse velocity component, is considered. A closure condition is formulated, which relates the shear stress to the first and second derivatives of the longitudinal mean velocity. The closure condition is derived without invoking any special hypotheses on the nature of turbulent motion, only taking advantage of the fact that the flow depends on a finite number of governing parameters. By virtue of the closure condition, the momentum equation is reduced to the boundary–value problem for a second–order differential equation, which is solved by the method of matched asymptotic expansions at high values of the logarithm of the Reynolds number based on the friction velocity. A limiting transpiration velocity is obtained, such that the shear stress at the injection wall vanishes, while the maximum point on the velocity profile approaches the suction wall. In this case, a sublayer near the suction wall appears where the mean velocity is proportional to the square root of the distance from the wall. A friction law for Poiseuille flow with transpiration is found, which makes it possible to describe the relation between the wall shear stress, the Reynolds number, and the transpiration velocity by a function of one variable. A velocity defect law, which generalizes the classical law for the core region in a channel with impermeable walls to the case of transpiration, is also established. In similarity variables, the mean velocity profiles across the whole channel width outside viscous sublayers can be described by a one–parameter family of curves. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Moving contact problems involving a rigid punch and a functionally graded coating

Frictional contact mechanics analysis for a rigid moving punch of an arbitrary profile and a functionally graded coating/homogeneous substrate system is carried out. The rigid punch slides over the coating at a constant subsonic speed. Smooth variation of the shear modulus of the graded coating is defined by an exponential function and the variation of the Poisson's ratio is assumed negligible. Coulomb's friction law is adopted. Hence, tangential force is proportional to the normal applied force through the coefficient of friction. An analytical method is developed utilizing the singular integral equation approach. Governing partial differential equations are derived in accordance with the theory of elastodynamics. The mixed boundary value problem is reduced to a singular integral equation of the second kind, which is solved numerically by an expansion-collocation technique. Presented results illustrate the effects of punch speed, coefficient of friction, material inhomogeneity and coating thickness on contact stress distributions and stress intensity factors. Comparisons indicate that the difference between elastodynamic and elastostatic solutions tends to be quite larger especially at higher punch speeds. It is shown that use of the elastodynamic theory provides more realistic results in contact problems involving a moving punch.     

弹性力学平面问题的精确元法

本文在阶梯折算法[1]和精确解析法[2]的基础上,提出构造有限元的新方法——精确元法.该方法不用一般变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法得到弹性力学平面问题的一个非协调任意四边形单元.它具有八个自由度.由于没有采用雅可比变换,该单元可以蜕化为三角形单元,在工程中使用起来较为方便.文中给出收敛性证明.文末给出算例,位移和应力均给出较好的结果,在单元的节点上有较好的数值精度.     

Steady-states for shear flows of a liquid-crystal model: Multiplicity, stability, and hysteresis

In this work, we study shear flows of a fluid layer between two solid blocks via a liquid-crystal type model proposed in [C.H.A. Cheng, L.H. Kellogg, S. Shkoller, D.L. Turcotte, A liquid-crystal model for friction, Proc. Natl. Acad. Sci. USA 21 (2007) 1-5] for an understanding of frictions. A characterization on the existence and multiplicity of steady-states is provided. Stability issue of the steady-states is examined mainly focusing on bifurcations of zero eigenvalues. The stability result suggests that this simple model exhibits hysteresis, and it is supported by a numerical simulation.     

波形钢腹板箱梁的约束扭转剪应力分析

基于周边不变形理论,结合闭口薄壁杆件约束扭转的计算分析,研究了波形钢腹板箱梁在约束扭转时混凝土悬臂板上扭转剪应力的分布,并进行了计算.通过对悬臂板在约束扭转中剪力流计算公式的推导,进一步阐述了其自由扭转剪应力及翘曲扭转剪应力的分布,指出了相关文献在这部分计算中存在的问题.通过一个简支波形钢腹板组合箱梁算例,将该文方法计算结果与ANSYS有限元计算结果进行比较.结果表明:在波形钢腹板箱梁截面中,主要由波形钢腹板承受扭转剪应力,其次是混凝土底板,底板剪应力最大值发生在底板中心处,其数值近似等于腹板剪应力的一半,而混凝土顶板和悬臂部分的扭转剪应力很小;该文计算的扭转剪应力结果在总体上符合有限元得到的扭转剪应力分布规律,在悬臂自由端为0,随着离开悬臂自由端距离的增大,扭转剪应力逐渐增大并达到峰值.