The stokes flow from half-space into semi-infinite circular cylinder

The infinite-series solutions for the creeping motion of a viscous incompressible fluid from half-space into semi-infinite circular cylinder are presented. The results show that inside the cylinder beyond a distance equal to 0.5 times the radius of the tube from the pore opening, the deviation of the velocity profile from Poiseuille flow is less than 1%. The inlet length in this case is comparable to that computed for a finite circular cylinder pore by Dagan et al.^[1]. In the half-space outside the cylinder pore region, the flow is strongly affected by the wall. Beyond one radius of the tube from the orifice, the solutions match almost exactly the flow through an orifice of zero thickness given by Sampson^[2]. The relationship between the pressure drop and the volumetric flow rate is also computed in the present paper for the semi-infinite tube.     

Viscous-inviscid interacting flow theory

In this paper a viscous-inviscid interacting flow theory (IFT) is developed for an incompressible, two—dimensional laminar flow. IFT's main points are as follows. (1) By introducing a concept of interaction layer where the normal momentum exchange is dominating, a new three-layer structure is established. (2) Through the conventional manipulations and by introducing an interaction model, both the streamwise and normal length scales are proved to be functions of a single parameterm, which is related to the streamwise pressure gradient and Reynolds number. (3) The approximate equations governing the flow of each layer as well as the whole interaction flow are derived. The present IFT is applicable to both attached and attached-separation bubble—reattached flows. The classical boundary layer theory^[1] and Triple-deck theory^[2] are shown to be two special cases of the present theory underm=0 and 1/4, respectively. Furthermore IFT provides new distinctions of both the normal and streamwise length scales for flow-field numerical computation and also gives a new approach to developing the simplified Navier-Stokes (SNS) equations. The project is supported by the National Natural Science Foundation of China.     

On the interplay between material flaws and dynamic necking

In this paper we investigate the interplay between material defects and flow localization in elastoplastic bars subjected to dynamic tension. For that task, we have developed a 1D finite difference scheme within a large deformation framework in which the material is modelled using rate-dependent J₂ plasticity. A perturbation of the initial yield stress is introduced in each node of the finite difference mesh to model localized material flaws. Numerical computations are carried out within a wide spectrum of strain rates ranging from 500 s⁻¹ to 2500 s⁻¹. On the one hand, our calculations reveal the effect of the material defects in the necking process. On the other hand, our results show that the necking inception, instead of being a random type process, is the deterministic result of the interplay between the mechanical behaviour of the material and the boundary conditions. This conclusion agrees with the experimental evidence reported by Rittel et al. [1] and Rotbaum et al. [2].     

Drag law for monodisperse gas–solid systems using particle-resolved direct numerical simulation of flow past fixed assemblies of spheres

Gas–solid momentum transfer is a fundamental problem that is characterized by the dependence of normalized average fluid–particle force F on solid volume fraction ? and the Reynolds number based on the mean slip velocity Re_m. In this work we report particle-resolved direct numerical simulation (DNS) results of interphase momentum transfer in flow past fixed random assemblies of monodisperse spheres with finite fluid inertia using a continuum Navier–Stokes solver. This solver is based on a new formulation we refer to as the Particle-resolved Uncontaminated-fluid Reconcilable Immersed Boundary Method (PUReIBM). The principal advantage of this formulation is that the fluid stress at the particle surface is calculated directly from the flow solution (velocity and pressure fields), which when integrated over the surfaces of all particles yields the average fluid–particle force. We demonstrate that PUReIBM is a consistent numerical method to study gas–solid flow because it results in a force density on particle surfaces that is reconcilable with the averaged two-fluid theory. The numerical convergence and accuracy of PUReIBM are established through a comprehensive suite of validation tests. The normalized average fluid–particle force F is obtained as a function of solid volume fraction ? (0.1 ? ? ? 0.5) and mean flow Reynolds number Re_m (0.01 ? Re_m ? 300) for random assemblies of monodisperse spheres. These results extend previously reported results of  and  to a wider range of ?, Re_m, and are more accurate than those reported by Beetstra et al. (2007). Differences between the drag values obtained from PUReIBM and the drag correlation of Beetstra et al. (2007) are as high as 30% for Re_m in the range 100–300. We take advantage of PUReIBM’s ability to directly calculate the relative contributions of pressure and viscous stress to the total fluid–particle force, which is useful in developing drag correlations. Using a scaling argument, Hill et al. (2001b) proposed that the viscous contribution is independent of Re_m but the pressure contribution is linear in Re_m (for Re_m > 50). However, from PUReIBM simulations we find that the viscous contribution is not independent of the mean flow Reynolds number, although the pressure contribution does indeed vary linearly with Re_m in accord with the analysis of Hill et al. (2001b). An improved correlation for F in terms of ? and Re_m is proposed that corrects the existing correlations in Re_m range 100–300. Since this drag correlation has been inferred from simulations of fixed particle assemblies, it does not include the effect of mobility of the particles. However, the fixed-bed simulation approach is a good approximation for high Stokes number particles, which are encountered in most gas–solid flows. This improved drag correlation can be used in CFD simulations of fluidized beds that solve the average two-fluid equations where the accuracy of the drag law affects the prediction of overall flow behavior.     

An experimental study of the regimes of motion of spheres falling or ascending freely in a Newtonian fluid

This paper presents the results of an experimental investigation aimed at verifying some of the interesting conclusions of the numerical study by Jenny et al. concerning the instability and the transition of the motion of solid spheres falling or ascending freely in a Newtonian fluid. The phenomenon is governed by two dimensionsless parameters: the Galileo number G, and the ratio of the density of the spheres to that of the surrounding fluid ρ_s/ρ. Jenny et al. showed that the (G, ρ_s/ρ) parameter space may be divided into regions with distinct features of the trajectories followed eventually by the spheres after their release from rest. The characteristics of these ‘regimes of motion’ as described by Jenny et al., agree well with what was observed in our experiments. However, flow visualizations of the wakes of the spheres using a Schlieren optics technique raise doubts about another conclusion of Jenny et al., namely the absence of a bifid wake structure.     

The equilibrium range cascades of wind-generated waves

It was established theoretically by Zakharov and Filonenko (1966) that the direct cascade region of the equilibrium range of the spectrum should follow an f⁻⁴ variation. This has since been verified from experimental data, by Toba (1973), Donclan et al (1985) and others. In this study we present a numerical verification of this f⁻⁴ variation, assuming physically realistic parameterizations for nonlinear wave-wave interactions, S_nl, for energy input to waves by the wind, S_in and removed by wave - breaking dissipation, S_ds.     

The application of the finite element method to solving Biot's consolidation equation

Biot's theory of consolidation of saturated soil, regards the consolidation process as a coupling problem between stress of elastic body and flow of fluid existing in pores^[1]. It can more precisely reflect the mechanism of consolidation than Terzhigi's theory^[2]. In this article, we obtain the general Biot's finite element equations of consolidation with classical variational principles. The equations have clear physical meaning and have been applied to analysing the consolidation of Bajiazui earth dam. The computational results are in accord with engineering practice.     

A similarity solution for the flow and heat transfer over a moving permeable flat plate in an external free stream: case of strong injection

The previous work of Bachok et al. (Heat Mass Transf. 47:1643–1649, 2011) on the forced convection heat transfer on an isothermal moving surface in an external free stream is extended to the case when fluid injection through the surface, characterized by the parameter γ, is large. The asymptotic solution derived in this limit shows that the boundary layer has a double region structure, with an inviscid inner region of thickness O(γ) and an outer shear layer. Some further aspects of the original problem not treated in Bachok et al. (Heat Mass Transf. 47:1643–1649, 2011) are discussed as well as the analogous problem for a constant surface heat flux, where relatively small injection rates are seen to give rise to large increases in the surface temperature.     

A strength reliability model by Markov process of unidirectional composites with fibers placed in hexagonal arrays

This paper proposes a strength reliability model based on a Markov process for unidirectional composites with fibers in a hexagonal array. The model assumes that a group of fiber breaking points, a so-called cluster, evolves with increased stress. The cluster evolution process branches because of various fiber-breakage paths. Load-sharing structure of intact fibers around clusters was estimated from geometric and mechanical local load-sharing rules. Composites fracture if a cluster achieves a critical size, so the model expresses a fracture criterion by setting an absorbing state. Next, the author constituted a state transition diagram concerning cluster evolutions of 1-fiber to 7-fiber breaks and analytically solved simultaneous differential equations obtained from the diagram. Results showed that, as critical cluster size increases, slope of the fracture probability distribution is given in a Weibull probability scale as follows: m_c=i×m_f (i, the number of broken fibers in a cluster; m_c and m_f, Weibull shape parameters for fracture probabilities of a critical cluster and fiber strength, respectively). This relation between m_c and m_f had been shown by Smith et al. [Proc. R. Soc. London, A 388 (1983) 353–391], but the present study demonstrated it analytically without any lower tail of the Weibull distribution used in that paper. In addition, the present model can be approximated by a one-state birth model.     

Nonlinear MHD Kelvin-Helmholtz instability in a pipe

Nonlinear MHD Kelvin-Helmholtz (K-H) instability in a pipe is treated with the derivative expansion method in the present paper. The linear stability problem was discussed in the past by Chandrasekhar (1961)^[1] and Xu et al. (1981).^[6]Nagano (1979)^[3] discussed the nonlinear MHD K-H instability with infinite depth. He used the singular perturbation method and extrapolated the obtained second order modifier of amplitude vs. frequency to seek the nonlinear effect on the instability growth rate γ. However, in our view, such an extrapolation is inappropriate. Because when the instability sets in, the growth rates of higher order terms on the right hand side of equations will exceed the corresponding secular producing terms, so the expansion will still become meaningless even if the secular producing terms are eliminated. Mathematically speaking, it's impossible to derive formula (39) when γ ₀ ² is negative in Nagano's paper.^[3]Moreover, even as early as γ ₀ ² → O⁺, the expansion becomes invalid because the 2nd order modifier γ₂ (in his formula (56)) tends to infinity. This weakness is removed in this paper, and the result is extended to the case of a pipe with finite depth. Theproject is supported by the National Natural Science Foundation of China.     

Irreversible thermodynamics of triple junctions during the intergranular void motion under the electromigration forces

A rigorous reformulation of internal entropy production and the rate of entropy flow is developed for multi-component systems consisting of heterophases, interfaces and/or surfaces. The result is a well-posed moving boundary value problem describing the dynamics of curved interfaces and surfaces associated with voids and/or cracks that are intersected by grain boundaries. Extensive computer simulations are performed for void configuration evolution during intergranular motion. In particular we simulate evolution resulting from the action of capillary and electromigration forces in thin film metallic interconnects having a “bamboo” structure, characterized by grain boundaries aligned perpendicular to the free surface of the metallic film interconnects. Analysis of experimental data utilizing previously derived mean time to failure formulas gives consistent values for interface diffusion coefficients and enthalpies of voids. 3.0 × 10⁻⁶ exp(−0.62 eV/kT) m² s⁻¹ is the value obtained for voids that form in the interior of the aluminum interconnects without surface contamination. 6.5 × 10⁻⁶ exp(−0.84 eV/kT) m² s⁻¹ is obtained for those voids that nucleate either at triple junctions or at the grain boundary-technical surface intersections, where the chemical impurities may act as trap centers for hopping vacancies.     

Validation of averaged equations for thermal vibrational convection in near-critical fluidsValidation d'équations moyennées pour l'étude de la convection vibrationnelle dans les fluides supercritiques.

Within an averaging approach, the governing equations and effective boundary conditions describing both the average and pulsation motion of a near-critical fluid subjected to high-frequency vibrations are obtained. Vibrations induce the non-homogeneities in average temperature. Owing to these non-homogeneities, the average flows can be generated even in isothermal cavity under weightlessness. These flows are examined for 1D and 2D configurations. The direct numerical simulations fulfilled earlier confirm the averaged model, we obtain the same flow structures by essentially smaller requirements for computational time. To cite this article: A.Vorobev et al., C. R. Mecanique 332 (2004).     

An adaptation of the low Mach number approximation for supercritical fluid buoyant flows

This Note describes an acoustic filtering of the equations governing the supercritical fluid buoyant flow driven by a weak heating. The resulting low Mach number approximation takes into account the compressibility of the fluid with respect to the hydrostatic pressure. Using the direct numerical simulation of a supercritical fluid flow in the Rayleigh–Bénard configuration, we show that the density stratification may be taken into account without further numerical effort and is fundamental for the prediction of the convective instability threshold induced by a weak heating. To cite this article: G. Accary et al., C. R. Mecanique 333 (2005).     

Theoretical analysis of a class of mixed,C0 continuity formulations for general dipolar Gradient Elasticity boundary value problems

Mixed weak formulations, with two or three main (tensor) variables, are stated and theoretically analyzed for general multi-dimensional dipolar Gradient Elasticity (biharmonic) boundary value problems. The general structure of constitutive equations is considered (with and without coupling terms). The mixed formulations are based on various generalizations of the so-called Ciarlet–Raviart technique. Hence, C⁰ continuity conforming basis functions may be employed in the finite element approximations (or even, C⁻¹ basis functions for the Cauchy stress variable). All the complicated boundary conditions, especially in the multi-dimensional scenario, are naturally considered. The main variables are the displacement vector, the double stress tensor and the Cauchy stress tensor. The latter variable may be eliminated in some of the formulations, depending on the structure of the constitutive equations. The standard continuous and discrete Babuška–Brezzi inf–sup conditions for the constraint equation, as well as, solution uniqueness for both the continuous statements and discrete approximations, are established in all cases. For the purpose of completeness, two one-dimensional mixed formulations are also analyzed. The respective constitutive equations possess general structure (with coupling terms). For the 1-D formulations, all the inf–sup conditions are satisfied, for both the continuous and discrete statements (assuming proper selection of the polynomial spaces for the main variables). Hence, the general Babuška–Brezzi theory results in quasi-optimality and stability. For multi-dimensional problems, the difficulty of deducing the inf–sup condition on the kernel is examined. Certain aspects of methodologies employed to theoretically by-pass this problem, are also discussed.     

Exact and numerical solutions of a fully developed generalized second-grade incompressible fluid with power-law temperature-dependent viscosity

We compute exact and numerical solutions of a fully developed flow of a generalized second-grade fluid, with power-law temperature-dependent viscosity (μ=θ^-M), down an inclined plane. Analytical solutions are found for the case when M=m+1, m≠1, m being a constant that models shear thinning (m<0) or shear thickening (m>0). The exact solutions are given in terms of Bessel functions. The numerical solutions indicate that both the velocity and the temperature increase with decreasing Froude number and that there is a critical value of Fr below which temperature “overshoots” its free surface value of unity. This phenomena is not reported in the work of Massoudi and Phuoc [Fully developed flow of a modified second grade fluid with temperature dependent viscosity, Acta Mech. 150 (2001) 23-37.] for viscosity that depends exponentially on temperature.     

The bending of elastic circular ring of non-homogeneous and variable cross section under the actions of arbitrary loads

On the basis of the stepped reduction method suggested in [1], we investigate the problem of the bending of elastic circular ring of non-homogeneous and variable cross section under the actions of arbitrary loads. The general solution of this problem is obtained so that it can be used for the calculations of strength and rigidity of practical problems such as arch, tunnel etc. In order to examine results of this paper and explain the application of this new method, an example is brought out at the end of this paper. Circular ring and arch are commonly used structures in engineering. Timoshenko, S.^[2], Barber, J. R.^[3], Tsumura Rimitsu^[4] et al. have studied these problems of bending, but, so far as we know, it has been solely restricted to the general solution of homogeneous uniform cross section ring. The only known solution for the problems with variable cross section ones has been solely restricted to the solution of special case of flexural rigidity in linear function of coordinates. On account of fundamental equations of the non-homogeneous variable cross section problem being variable coefficients, it is very difficult to solve them. In this paper, we use the stepped reduction method suggested in [1] to transform the variable coefficient differential equation into equivalent constant coefficient one. After introducing virtual internal forces, we obtain general solution of an elastic circular ring with non-homogeneity and variable cross section under the actions of arbitrary loads.     

A new procedure for the prediction of forced convection heat transfer at near- and supercritical pressure

Based on heat transfer experiments with tubes of different inside geometry, correlations for forced convection heat transfer at near- and supercritical pressure are presented. They are of the Dittus-Boelter typeNu=CRe ^m Pr ⁿ and over the whole enthalpy range. The simple form is made possible by a new method of determining a representative specific heat capacity. At five reference temperatures heat capacities are computed with semi-empirical equations. The two highest values closest to the (pseudo-)critical state are then sorted out, because they are not relevant for turbulent flow. The average of the remaining three values represents a reasonable characteristic heat capacity. Comparisons to measurements show the very good accuracy of the new correlations. The procedure can be recommended whenever strong dependencies of fluid properties on temperature or even singularities occur.     

On Boussinesq's paradigm in nonlinear wave propagation

Boussinesq's original derivation of his celebrated equation for surface waves on a fluid layer opened up new horizons that were to yield the concept of the soliton. The present contribution concerns the set of Boussinesq-like equations under the general title of ‘Boussinesq's paradigm’. These are true bi-directional wave equations occurring in many physical instances and sharing analogous properties. The emphasis is placed: (i) on generalized Boussinesq systems that involve higher-order linear dispersion through either additional space derivatives or additional wave operators (so-called double-dispersion equations); and (ii) on the ‘mechanics’ of the most representative localized nonlinear wave solutions. Dissipative cases and two-dimensional generalizations are also considered. To cite this article: C.I. Christov et al., C. R. Mecanique 335 (2007).     

A new method for obtaining the exact solutions of the Navier-Stokes (NS) equations and its applications

In this paper we propose a new method for obtaining the exact solutions of the Mavier-Stokes (NS) equations for incompressible viscous fluid in the light of the theory of simplified Navier-Stokes (SNS) equations developed by the first author^[1,2], Using the present method we can find some new exact solutions as well as the well-known exact solutions of the NS equations. In illustration of its applications, we give a variety of exact solutions of incompressible viscous fluid flows for which NS equations of fluid motion are written in Cartesian coordinates, or in cylindrical polar coordinates, or in spherical coordinates. The project supported by National Natural Science Foundation of China.