Solution of the problem of flow of a non-axisymmetric swirling submerged jets

This paper considers the problem of a non-axisymmetric swirling jet of an incompressible viscous fluid flowing in a space flooded with the same fluid. The far field of the jet is studied under the assumption that the angular momentum vector corresponding to the swirling of the jet is not collinear to the momentum vector of the jet. It is shown that the main terms of the asymptotic expansion of the full solution for the velocity field are determined by the exact integrals of conservation of momentum, mass, and angular momentum. An analytical solution of the problem describing the axisymmetric swirling jet is obtained.     

Local instant formulation of two-phase flow

The local instant formulation of mass, momentum and energy conservations of two-phase flow has been developed. Distribution, an extended notion of a function, has been introduced for this purpose because physical parameters of two-phase flow media change discontinuously at the interface and the Lebesgue measure of an interface is zero. Using a characteristic function of each phase, the physical parameters of two-phase flow have been defined as field quantities. In addition to this, the source terms at the interface are defined in terms of the local instant interfacial area concentration. Based on these field quantities, the local instant field equations of mass, momentum and total energy conservations of two-phase flow have been derived. Modification of these field equations gives the single field representation of the local instant field equations of two-phase flow. Neglecting the interfacial force and energy, this formulation coincides with the field equations of single-phase flow, except in the definition of differentiation. The local instant two-fluid formulation of two-phase flow has also been derived. This formulation consists of six local instant field equations of mass, momentum and total energy conservations of both phases. Interfacial mass, momentum and energy transfer terms appear in these equations, which are expressed in terms of the local instant interfacial area concentration.     

非正交同位网格下压力与速度耦合算法

发展了一种在非正交同位网格下以笛卡儿速度分量作为动量方程的独立变量、压力与速度耦合的S IM-PLER算法。该算法的特点是显式处理界面速度中的压力交叉导数项,得出压力与压力修正方程,使得压力及压力修正值与界面逆变速度直接耦合。通过对分汊通道内的流动问题进行验证计算,结果表明该算法可以有效而准确地模拟复杂区域内的流动与换热问题。     

On the Microscopic Interpretation of Stress and Couple Stress

Exact continuum forms of balance (for mass, linear momentum, and tensor-valued moment of momentum) are established as relations between weighted spatial averages of corpuscular quantities computed at any supra-molecular length scale. Explicit expressions for stress and generalised couple stress in terms of particle interactions are obtained using a theorem due to Noll, and their physical interpretation is discussed for a specific choice of weighting function. Remarks are made on other choices of weighting function, the interpretation of partial stress in mixture theory, a link between couple stress and inhomogeneity, and other forms of moment of momentum balance. Comparison is made with the statistical mechanical viewpoint pioneered by Irving and Kirkwood.     

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.     

Extensions to the Macroscopic Navier–Stokes Equation

A development is provided showing that for any phase, by not neglecting the macroscopic terms of the deviation from the intensive momentum and of the dispersive momentum, we obtain a macroscopic secondary momentum balance equation coupled with a macroscopic dominant momentum balance equation that is valid at a larger spatial scale. The macroscopic secondary momentum balance equation is in the form of a wave equation that propagates the deviation from the intensive momentum while concurrently, in the case of a Newtonian fluid and under certain assumptions, the macroscopic dominant momentum balance equation may be approximated by Darcys equation to address drag dominant flow. We then develop extensions to the dominant macroscopic Navier–Stokes (NS) equation for saturated porous matrices, to account for the pressure gradient at the microscopic solid-fluid interfaces. At the microscopic interfaces we introduce the exchange of inertia between the phases, accounting for the relative fluid square velocities and the rate of these velocities, interpreted as Forchheimer terms. Conditions are provided to approximate the extended dominant NS equation by Forchheimer quadratic momentum law or by Darcys linear momentum law. We also show that the dominant NS equation can conform into a nonlinear wave equation. The one-dimensional numerical solution of this nonlinear wave equation demonstrates good qualitative agreement with experiments for the case of a highly deformable elasto-plastic matrix.     

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.     

On steady motions of a rigid body bearing a two-degree-of-freedom control momentum gyroscope with precession axis arbitrarily oriented in the carrying body

Steady motions of a rigid body with a control momentum gyroscope are studied versus the gimbal axis direction relative to the body and the magnitude of the system angular momentum. The study is based on a formula that gives a parametric representation of the set of the system steady motions in terms of the rotation angle of the gimbal. It is shown that, depending on the values of the parameters, the system has 8, 12, or 16 steady motions and the number of stable motions is 2 or 4.     

Evolution of the balance equations in saturated thermoelastic porous media following abrupt simultaneous changes in pressure and temperature

A mathematical model is developed for saturated flow of a Newtonian fluid in a thermoelastic, homogeneous, isotropic porous medium domain under nonisothermal conditions. The model contains mass, momentum and energy balance equations. Both the momentum and energy balance equations have been developed to include a Forchheimer term which represents the interaction at the solid-fluid interface at high Reynolds numbers. The evolution of these equations, following an abrupt change in both fluid pressure and temperature, is presented. Using a dimensional analysis, four evolution periods are distinguished. At the very first instant, pressure, effective stress, and matrix temperature are found to be disturbed with no attenuation. During this stage, the temporal rate of pressure change is linearly proportional to that of the fluid temperature. In the second time period, nonlinear waves are formed in terms of solid deformation, fluid density, and velocities of phases. The equation describing heat transfer becomes parabolic. During the third evolution stage, the inertial and the dissipative terms are of equal order of magnitude. However, during the fourth time period, the fluid's inertial terms subside, reducing the fluid's momentum balance equation to the form of Darcy's law. During this period, we note that the body and surface forces on the solid phase are balanced, while mechanical work and heat conduction of the phases are reduced.     

Curvature dependent surface energy for a free standing monolayer graphene: Some closed form solutions of the non-linear theory

Continuum modeling of a free-standing graphene monolayer, viewed as a two dimensional 2-lattice, requires specifications of the components of the shift vector that act as an auxiliary variable. The field equations are then the equations ruling the shift vector, together with momentum and moment of momentum equations. We present an analysis of simple loading histories such as axial, biaxial tension/compression and simple shear for a range of problems of increasing difficulty. We start by laying down the equations of a simplified model which can still capture bending effects. Initially, we ignore out of plane deformations. For this case, we solve analytically the equations ruling the auxiliary variables in terms of the shift vector; these equations are algebraic when the loading is specified. As a next step, still working on the simplified model, out-of-plane deformations are taken into account and the equations complicate dramatically. We describe how wrinkling/buckling can be introduced into the model and apply the Cauchy–Kowalevski theorem to get existence and uniqueness in terms of the shift vector for some characteristic cases. Finally, for the treatment of the most general problem, we classify the equations of momentum and give conditions for the Cauchy–Kowalevski theorem to apply.     

Computation of Jump Coefficients for Momentum Transfer Between a Porous Medium and a Fluid Using a Closed Generalized Transfer Equation

The momentum transfer between a homogeneous fluid and a porous medium in a system analogous to the one used by Beavers and Joseph (J Fluid Mech 30:197–207, 1967) is studied using volume averaging techniques. In this article, we present a closed generalized momentum transport equation (GTE) that is valid everywhere and is expressed in terms of position-dependent effective transport coefficients, which are computed from the solution of associated closure problems previously reported. A combination of the velocity profiles from the GTE in the definition of the excess terms that define the jump coefficients allows their computation using numerical techniques. The calculations are in concordance with those resulting from the work of Goyeau et al. (Int J Heat Mass Transf. 46:4071–4081, 2003), showing a strong dependence with the porosity. In addition, the effects of the roughness of the boundary on the computation of the position-dependent permeability tensor in the inter-region are also analyzed.     

Depth‐averaged two‐dimensional curvilinear explicit finite analytic model for open‐channel flows

A depth‐averaged two‐dimensional model has been developed in the curvilinear co‐ordinate system for free‐surface flow problems. The non‐linear convective terms of the momentum equations are discretized based on the explicit–finite–analytic method with second‐order accuracy in space and first‐order accuracy in time. The other terms of the momentum equations, as well as the mass conservation equation, are discretized by the finite difference method. The discretized governing equations are solved in turn, and iteration in each time step is adopted to guarantee the numerical convergence. The new model has been applied to various flow situations, even for the cases with the presence of sub‐critical and supercritical flows simultaneously or sequentially. Comparisons between the numerical results and the experimental data show that the proposed model is robust with satisfactory accuracy. Copyright © 2000 John Wiley & Sons, Ltd.     

Contributions to Theoretical/Experimental Developments in Shock Waves Propagation in Porous Media

Macroscopic balance equations of mass, momentum and energy for compressible Newtonian fluids within a thermoelastic solid matrix are developed as the theoretical basis for wave motion in multiphase deformable porous media. This leads to the rigorous development of the extended Forchheimer terms accounting for the momentum exchange between the phases through the solid-fluid interfaces. An additional relation presenting the deviation (assumed of a lower order of magnitude) from the macroscopic momentum balance equation, is also presented. Nondimensional investigation of the phases' macroscopic balance equations, yield four evolution periods associated with different dominant balance equations which are obtained following an abrupt change in fluid's pressure and temperature. During the second evolution period, the inertial terms are dominant. As a result the momentum balance equations reduce to nonlinear wave equations. Various analytical solutions of these equations are described for the 1-D case. Comparison with literature and verification with shock tube experiments, serve as validation of the developed theory and the computer code.A 1-D TVD-based numerical study of shock wave propagation in saturated porous media, is presented. A parametric investigation using the developed computer code is also given.     

On the Effect of Vortex Grid Density in the Vortex-in-cell Simulation of Mixing Layers

The two-dimensional vortex-in-cell method is used in the simulation of a spatially growing mixing layer at a high Reynolds number. Criteria, showing that the vortex grid density plays a role in obtaining a converged solution, are presented. Numerical experiments are conducted to test the criteria by varying the number of vortices and the grid sizes. The effects on the momentum thickness and on the peak values of the root-mean-square velocity fluctuations, and negative cross-stream correlation are discussed in terms of the vortex grid density.     

A robust algorithm for computing fluid flows on highly non‐smooth staggered grids

A new method for computing the fluid flow in complex geometries using highly non‐smooth and non‐orthogonal staggered grid is presented. In a context of the SIMPLE algorithm, pressure and physical tangential velocity components are used as dependent variables in momentum equations. To reduce the sensitivity of the curvature terms in response to coordinate line orientation change, these terms are exclusively computed using Cartesian velocity components in momentum equations. The method is then used to solve some fairly complicated 2‐D and 3‐D flow field using highly non‐smooth grids. The accuracy of results on rough grids (with sharp grid line orientation change and non‐uniformity) was found to be high and the agreement with previous experimental and numerical results was quite good. Copyright © 2008 John Wiley & Sons, Ltd.     

The Eshelby stress tensor,angular momentum tensor and dilatation flux in gradient elasticity

The Eshelby (static energy momentum) stress tensor, the angular momentum tensor and the dilatation flux are derived for anisotropic linear gradient elasticity in non-homogeneous materials. The divergence of these tensors gives the configurational forces, moments and work terms in gradient elasticity. There are several types of configurational forces, acting on the dislocation density and its gradient, on the inhomogeneities, proportional to the distortion, and linear and quadratic in the distortion gradient, and on the body force.     

Momentum equation for straight electrically charged jet

A one-dimensional momentum conservation equation for a straight jet driven by an electrical field is developed. It is presented in terms of a stress component, which can be applied to any constitutive relation of fluids. The only assumption is that the fluid is incompressible. The results indicate that both the axial and radial constitutive relations are required to close the governing equations of the straight charged jet. However, when the trace of the extra stress tensor is zero, only the axial constitutive relation is required. It is also found that the second normal stress difference for the charged jet is always zero. The comparison with other developed momentum equations is made.     

On the eddy diffusivities for momentum and heat

Summary From a simple consideration, the eddy diffusivities for momentum and heat are shown to be equal in terms of statistical properties of turbulence. A universal semi-empirical expression for the eddy diffusivities in a turbulent shear flow is obtained in terms of the distance from the solid boundary and flow properties.     

Laminar flow and heat transfer in a periodically converging-diverging channel

A finite difference solution for laminar viscous flow through a sinusoidally curved converging-diverging channel is presented. The physical wavy domain is transformed into a rectangular computational domain in order to simplify the application of boundary conditions on the channel walls. The discretized conservation equations for mass, momentum and energy are derived on a control volume basis. The pseudo-diffusive terms that arise from the co-ordinate transformation are treated as source terms, and the resulting system of equations is solved by a semi-implicit procedure based on line relaxation. Results are obtained for both the developing and the fully developed flow for a Prandtl number of 0.72, channel maximum width-to-pitch ratio of 1.0, Reynolds number ranging from 100 to 500 and wall amplitude-to-pitch ratio varying from 0.1 to 0.25. Results are presented here for constant fluid properties and for a prescribed wall enthalpy only.     

Linear dynamic model for porous media saturated by two immiscible fluids

A linear isothermal dynamic model for a porous medium saturated by two immiscible fluids is developed in the paper. In contrast to the mixture theory, phase separation is avoided by introducing one energy for the porous medium. It is an important advantage of the model based on one energy approach that it can account for the couplings between the phases. The volume fraction of each phase is characterized by the porosity of the porous medium and the saturation of the wetting phase. The mass and momentum balance equations are constructed according to the generalized mixture theory. Constitutive relations for the stress, pore pressure are derived from the free energy function. A capillary pressure relaxation model characterizing one attenuation mechanism of the two-fluid saturated porous medium is introduced under the constraint of the entropy inequality. In order to describe the momentum interaction between the fluids and the solid, a frequency independent drag force model is introduced. The details of parameter estimation are discussed in the paper. It is demonstrated that all the material parameters in our model can be calculated by the phenomenological parameters, which are measurable. The equations of motion in the frequency domain are obtained in terms of the Fourier transformation. In terms of the equations of motion in the frequency domain, the wave velocities and the attenuations for three P waves and one S wave are calculated. The influences of the capillary pressure relaxation coefficient and the saturation of the wetting phase on the velocities and attenuation coefficients for the four wave modes are discussed in the numerical examples.