Evolution of governing mass and momentum balances following an abrupt pressure impact in a porous medium

A mathematical model is developed of an abrupt pressure impact applied to a compressible fluid flowing through a porous medium domain. Nondimensional forms of the macroscopic fluid mass and momentum balance equations yield two new scalar numbers relating storage change to pressure rise. A sequence of four reduced forms of mass and momentum balance equations are shown to be associated with a sequence of four time periods following the onset of a pressure change. At the very first time period, pressure is proven to be distributed uniformly within the affected domain. During the second time interval, the momentum balance equation conforms to a wave form. The behavior during the third time period is governed by the averaged Navier-Stokes equation. After a long time, the fourth period is dominated by a momentum balance similar to Brinkman's equation which may convert to Darcy's equation when friction at the solid-fluid interface dominates.     

Turbulent sediment transport using momentum balance — The strong turbulence approximation

The turbulent flow of a sediment-fluid mixture near a flat bottom is studied. An equation is derived expressing the balance between the gravitational force on the sediment and the turbulent sediment momentum flux. A mixing length theory is used to model the turbulent momentum flux, with the mixing lengths assumed proportional to the distance from the bottom. The ratio of the sediment momentum mixing length to the mixture momentum mixing length is assumed to be a universal constant. An approximate solution for the mean sediment flux profile is derived for situations where the turbulent strength near the bottom dominates the gravitational effects there. The resulting concentration profile is constant near the bottom and drops off rapidly above a certain predicted depth.     

On canonical equations of continuum thermomechanics

It is shown that the canonical balance of momentum of continuum mechanics can be formulated in a general way, but not independently of the usual balance of linear momentum, even in the absence of specified constitutive equations. A parallel construct is made of necessity for the accompanying time-like canonical energy equation. On specifying the energy, previous particular cases can be deduced including pure elasticity, inhomogeneous thermoelasticity of conductors, and the case of dissipative solid-like materials described by means of a diffusive internal variable (such as in damage or weakly non-local plasticity). A redefinition of the entropy flux is necessarily accompanied by a redefinition of the Eshelby stress tensor.     

Wave momentum and scattering of elastic waves by two-dimensional thin objects

The elastic wave momentum equation is applied to scattering of dilatational and shear waves by two-dimensional thin objects. It is shown that the sources of wave momentum are located at the edges of these objects. For a stress-zero crack or for a rigid inclusion there are two sources at each edge, for a fluid-filled crack there is just one. The scattered wave is expressed in terms of these sources. This reduces the number of independent variables by a factor two. An application to inverse scattering problems is also given.     

A momentum coupling method for the unsteady incompressible Navier–Stokes equations on the staggered grid

A new numerical method is developed to efficiently solve the unsteady incompressible Navier–Stokes equations with second-order accuracy in time and space. In contrast to the SIMPLE algorithms, the present formulation directly solves the discrete x- and y-momentum equations in a coupled form. It is found that the present implicit formulation retrieves some cross convection terms overlooked by the conventional iterative methods, which contribute to accuracy and fast convergence. The finite volume method is applied on the fully staggered grid to solve the vector-form momentum equations. The preconditioned conjugate gradient squared method (PCGS) has proved very efficient in solving the associate linearized large, sparse block-matrix system. Comparison with the SIMPLE algorithm has indicated that the present momentum coupling method is fast and robust in solving unsteady as well as steady viscous flow problems. © 1998 John Wiley & Sons, Ltd.     

NUMERICAL SIMULATION OF SHALLOW WATER WAVE PROPAGATION USING A BOUNDARY ELEMENT WAVE EQUATION MODEL

The present paper makes use of a wave equation formulation of the primitive shallow water equations to simulate one-dimensional free surface flow. A numerical formulation of the boundary element method is then developed to solve the wave continuity equation using a time-dependent fundamental solution, while an explicit finite difference scheme is used to derive velocities from the primitive momentum equation. One-dimensional free surface flows in open channels are treated and the results compared with analytical and numerical solutions. © 1997 John Wiley & Sons, Ltd.     

Consistent discretization for simulations of flows with moving generalized curvilinear coordinates

We develop a consistent discretization of conservative momentum and scalar transport for the numerical simulation of flow using a generalized moving curvilinear coordinate system. The formulation guarantees consistency between the discrete transport equation and the discrete mass conservation equation due to grid motion. This enables simulation of conservative transport using generalized curvilinear grids that move arbitrarily in three dimensions while maintaining the desired properties of the discrete transport equation on a stationary grid, such as constancy, conservation, and monotonicity. In addition to guaranteeing consistency for momentum and scalar transport, the formulation ensures geometric conservation and maintains the desired high‐order time accuracy of the discretization on a moving grid. Through numerical examples we show that, when the computation is carried out on a moving grid, consistency between the discretized scalar advection equation and the discretized equation for flow mass conservation due to grid motion is required in order to obtain stable and accurate results. We also demonstrate that significant errors can result when non‐consistent discretizations are employed. Copyright © 2009 John Wiley & Sons, Ltd.     

Similarity and decay laws of wakes with zero momentum and angular momentum

In [1–4] the laws of decay of the average and fluctuating velocities in momentumless turbulent wakes were experimentally investigated with and without swirl. In [5, 6] unswirled momentumless wakes and in [7] wakes with a nonzero angular momentum were theoretically investigated. However, turbulent wakes with zero momentum and angular momentum were not covered by these investigations. This class of flows is the subject of the present study.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 35–41, September–October, 1993.     

Construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body

We consider the problem of construction of optimal laws of variation in the angular momentum vector of a dynamically symmetric rigid body so as to ensure the transition of the rigid body from an arbitrary initial angular position to the required final angular position. For the functionals to be minimized, we use combined performance functionals, one of which characterizes the expenditure of time and of the squared modulus of the angular momentum vector in a given proportion, while the other characterizes the expenditure of time and momentum of the modulus of the angular momentum vector necessary to change the rigid body orientation. The control (the vector of the rigid body angular momentum) is assumed to be bounded in the modulus. The problem is solved by using Pontryagin’s maximum principle and the quaternion differential equation [1, 2] relating the vector of the dynamically symmetric rigid body angular momentum to the quaternion of orientation of the coordinate system rotating with respect to the rigid body about its dynamical symmetry axis at an angular velocity proportional to the angular momentum vector projection on the axis. The use of such a model of rotational motion leads to the problem of optimal control with the moving right end of the trajectory and significantly simplifies the analytic study of the problem of construction of optimal laws of variation in the angular momentum vector, because this model explicitly exploits the body angular momentum quaternion (control) instead of the rigid body absolute angular velocity quaternion. We construct general analytic solutions of the differential equations for the boundary-value problems which form systems of nine nonlinear differential equations. It is shown that the process of solving the differential boundary-value problems is reduced to solving two scalar algebraic transcendental equations.     

非平整运动海底上n层流动中波浪传播的Hamilton逼近

在海洋水域，界面波对大尺度变化流的作用是一种典型的分层流动现象．考虑一不可压缩、无黏的分层势流运动，建立了一个在非平整运动海底上的n层流体演化系统，并对其进行了Hamilton描述．每层流体具有各自的常密度、均匀流水平速度，其厚度由未扰动和扰动部分构成．相对于顶层流体的自由表面，刚性、运动的海底具有一般地形变化特征．在明确指出n层流体运动的控制方程和各层交界面上的运动学、动力学边界条件(包含各层交界面上张力效应)后，对该分层流动力系统进行了Hamilton构造，即给出其正则方程和其下述的正则变量：各交界面位移和各交界面上的动量势密度差。     

On the momentum of Eulerian phonons

For a one-dimensional dispersive medium the linear momentum of a phonon is discussed in both the Lagrangian and the Eulerian picture, i.e. with the use of substantial (material) and local coordinates, respectively. As phonons are usually considered as solutions of the linearized equations of motion, in the Eulerian picture the linear momentum of a phonon is only defined up to linear terms in the fields. To obtain results relevant towards higher order in the fields, one has to solve the nonlinear equations of motion. This is done to obtain expressions for the linear momentum up to terms quadratic in the fields.     

Flow distribution in parallel and reverse flow manifolds

Numerical predictions of the flow distribution in parallel and reverse flow manifold systems have been obtained by a novel finite-difference procedure. A one-dimensional elliptic formulation is proposed. An iterative numerical scheme solves the differential equations for momentum in the longitudinal direction and continuity, while in the lateral direction an integral equation for momentum is solved. The predicted pressure distributions compare well with the available experimental results. A parametric computation has been performed to demonstrate the effect of the governing non-dimensional parameters on the flow distribution.     

Modeling and simulation of the coupled mechanical-electrical response of soft solids

Dielectric elastomer (DE) is one type of electro-active polymers (EAP) that responds to electrical stimulation with a significant shape and size change. As EAPs, dielectric elastomers are lightweight, inexpensive, pliable and can be fabricated into various shapes, all of which are attractive properties to justify the intense research in the field. This paper presents a nonlinear, electrical and mechanical coupled, large deformation finite element formulation for DEAs. Maxwell’s equations for the electroquasistatic fields were solved simultaneously with equation of linear momentum. The hyperelastic Ogden model and total Maxwell stress method were combined to describe the material. The formulation was based on the weak forms of Maxwell’s equation and linear momentum expressed in the reference configuration. The closed form consistent tangent moduli for dielectric elastomers were derived. The results of the simulation compared with the experiments have demonstrated the validity of the method from the computational aspect.     

Mechanical wave momentum from the first principles

Axial momentum carried by waves in a uniform waveguide is considered based on the conservation laws and a kind of the causality principle. Specifically, we examine (without resorting to constitutive data) steady-state waves of an arbitrary shape, periodic waves which speed differs from the speed of its form and binary waves carrying self-equilibrated momentum. The approach allows us to represent momentum as a product of the wave mass and the wave speed. The propagating wave mass, positive or negative, is the excess of that in the wave over its initial value. This general representation is valid for mechanical waves of arbitrary nature and intensity. The finite-amplitude longitudinal and periodic transverse waves are examined in more detail. It is shown in particular, that the transverse excitation of a string or an elastic beam results in the binary wave. The closed-form expressions for the drift in these waves functionally reduce to the Stokes’ drift in surface water waves (a half the latter by the value). Besides, based on the general representation an energy–momentum relation is discussed and the physical meaning of the so-called “wave momentum” is clarified.     

Momentum/continuity coupling with large non‐isotropic momentum source terms

Pressure‐based methods such as the SIMPLE algorithm are frequently used to determine a coupled solution between the component momentum equations and the continuity equation. This paper presents a colocated variable pressure correction algorithm for control volumes of polyhedral/polygonal cell topologies. The correction method is presented independent of spatial approximation. The presence of non‐isotropic momentum source terms is included in the proposed algorithm to ensure its applicability to multi‐physics applications such as gas and particulate flows. Two classic validation test cases are included along with a newly proposed test case specific to multiphase flows. The classic validation test cases demonstrate the application of the proposed algorithm on truly arbitrary polygonal/polyhedral cell meshes. A comparison between the current algorithm and commercially available software is made to demonstrate that the proposed algorithm is competitively efficient. The newly proposed test case demonstrates the benefits of the current algorithm when applied to a multiphase flow situation. The numerical results from this case show that the proposed algorithm is more robust than other methods previously proposed. Copyright © 2009 John Wiley & Sons, Ltd.     

Accurate numerical estimation of interphase momentum transfer in Lagrangian–Eulerian simulations of dispersed two-phase flows

The Lagrangian–Eulerian (LE) approach is used in many computational methods to simulate two-way coupled dispersed two-phase flows. These include averaged equation solvers, as well as direct numerical simulations (DNS) and large-eddy simulations (LES) that approximate the dispersed-phase particles (or droplets or bubbles) as point sources. Accurate calculation of the interphase momentum transfer term in LE simulations is crucial for predicting qualitatively correct physical behavior, as well as for quantitative comparison with experiments. Numerical error in the interphase momentum transfer calculation arises from both forward interpolation/approximation of fluid velocity at grid nodes to particle locations, and from backward estimation of the interphase momentum transfer term at particle locations to grid nodes. A novel test that admits an analytical form for the interphase momentum transfer term is devised to test the accuracy of the following numerical schemes: (1) fourth-order Lagrange Polynomial Interpolation (LPI-4), (3) Piecewise Cubic Approximation (PCA), (3) second-order Lagrange Polynomial Interpolation (LPI-2) which is basically linear interpolation, and (4) a Two-Stage Estimation algorithm (TSE). A number of tests are performed to systematically characterize the effects of varying the particle velocity variance, the distribution of particle positions, and fluid velocity field spectrum on estimation of the mean interphase momentum transfer term. Numerical error resulting from backward estimation is decomposed into statistical and deterministic (bias and discretization) components, and their convergence with number of particles and grid resolution is characterized. It is found that when the interphase momentum transfer is computed using values for these numerical parameters typically encountered in the literature, it can incur errors as high as 80% for the LPI-4 scheme, whereas TSE incurs a maximum error of 20%. The tests reveal that using multiple independent simulations and higher number of particles per cell are required for accurate estimation using current algorithms. The study motivates further testing of LE numerical methods, and the development of better algorithms for computing interphase transfer terms.     

Solution of the incompressible mass and momentum equations by application of a coupled equation line solver

This paper describes an iterative technique for solving the coupled algebraic equations for mass and momentum conservation for an incompressible fluid flow. The technique is based on the simultaneous solution for pressure and velocity along lines. In a manner similar to ADI methods for a single variable, the solution domain is entirely swept line-by-line in each co-ordinate direction successively until a converged solution is obtained. The tight coupling between the equations that is guaranteed by the method results in an economical solution of the equation set.     

Impact dynamics of multibody systems with frictional contact using joint coordinates and canonical equations of motion

This paper presents a methodology in computational dynamics for the analysis of mechanical systems that undergo intermittent motion. A canonical form of the equations of motion is derived with a minimal set of coordinates. These equations are used in a procedure for balancing the momenta of the system over the period of impact, calculating the jump in the body momentum, velocity discontinuities and rebounds. The effect of dry friction is discussed and a contact law is proposed. The present formulation is extended to open and closed-loop mechanical systems where the jumps in the constraints' momenta are also solved. The application of this methodology is illustrated with the study of impact of open-loop and closed-loop examples.     

Hamilton体系下介电弹性体圆形薄膜的动力学建模与辛求解

采用辛算法研究了Hamilton体系下介电弹性体圆形薄膜的动力学响应。首先,将该问题引入Hamilton对偶变量体系,借助Legendre变换,给出系统的广义动量和Hamilton函数,通过对Hamilton函数作用量的变分,得到Hamilton体系下的正则方程。其次,对于得到的正则方程给出了辛Runge-Kutta的计算格式。最后,采用二级四阶辛Runge-Kutta算法对动力学系统进行了数值求解,和四级四阶经典Runge-Kutta算法进行对比,结果表明,二级四阶辛Runge-Kutta算法具有保能量以及长时间数值稳定的优势,同时说明四级四阶经典Runge-Kutta算法对于步长依赖的局限性。