后台阶流动的数值模拟

访述了大涡模拟的基本思想，指出大涡模拟的效率主要取决于四个因素，即流动中须有大尺度涡存在、合理的计算格式、合适的滤波器和亚格子应力模型。在深入考虑粘性不可压缩流Navier—Stokes方程各个子项作用的基础上，提出二阶全展开Euler—Taylor—Galerkin有限元方法作为大涡模拟的离散格式，并采用Gauss滤波器，对典型算例——后台阶处的流动进行大涡模拟，计算结果与相关文献符合的很好。从计算结果还可以看出大涡模拟与二阶全展开ETG有限元方法的结合在捕捉涡系及反映涡动时变过程方面具有明显的优势，说明大涡模拟适合于边界几何形状复杂区域流动的模拟。同时应用二阶全展开ETG有限元方法对低雷诺数粘性不可压缩后台阶流动进行了计算，得到与相关文献符合良好的计算结果，即该方法也可独立用于对低雷诺数粘性不可压缩流动的计算。     

Numerical solution to the shearing flow of granular materials between two plates

We study the shearing flow of granular materials between two horizontal flat plates where the top plate is moving with a constant speed. The constitutive relation used for the stress is based on the continuum model proposed by Rajagopal and Massoudi (DOE Report, DOE/PETC/TR-90/3, 1990). The material coefficients such as viscosity and normal stress coefficients are based on the model of Boyle and Massoudi (Int. J. Eng. Sci 28 (1990) 1261). The governing equations are non-dimensionalized and the resulting system of non-linear differential equations is solved numerically using finite difference technique.     

Couette flow of granular materials

The flow of granular materials between rotating cylinders is studied using a continuum model proposed by Rajagopal and Massoudi (A method for measuring material moduli for granular materials: flow in an orthogonal rheometer, DOE/PETC/TR90/3, 1990). For a steady, fully developed condition, the governing equations are reduced to a system of coupled non-linear ordinary differential equations. The resulting boundary value problem is non-dimensionalized and is then solved numerically. The effect of material parameters, i.e., dimensionless numbers on the volume fraction and the velocity fields are studied.     

Relationship between mass discharge rate and granular temperature of silo flow with variance of outlets

A lower granular temperature zone near the orifice in silos was found by speckle visibility spectroscopy technology. Surprisingly, the values of granular temperature in this area are inversely proportional to the mass discharge rate (MDR). This finding that links the macroscopic MDR with mesoscopic dynamics temperature of the particle system, is obtained by looking at the effect that different outlet sizes of silos on MDR and granular temperature field. Apart from evidencing the relevance of a parameter that has been traditionally overlooked in silos, this relationship supposes a benchmark with which to explore the influence of small orifice sizes on MDR. Among these, we found that the frequency of intermittent clogging and their dependence on the size of the outlet.     

Numerical simulation study on particle breakage behavior of granular materials in confined compression tests

In order to study the fragmentation law, the confined compression experiment of granular assemblies has been conducted to explore the particle breakage characteristic by DEM approach in this work. It is shown that contact and contact force during the loading process gradually transform from anisotropy to isotropy. Meanwhile, two particle failure modes caused by different contact force states are analyzed, which are single-through-crack failure and multi-short-crack failure. Considering the vertical distribution of the number of cracks and the four characteristic stress distributions (the stress related to the maximum contact force, the major principal stress, the deviatoric stress and the mean stress), it is pointed out that the stress based on the maximum contact force and the major principal stress can reflect the distribution of cracks accurately. In addition, the size effect of particle crushing indicates that small size particles are prone to break. The lateral pressure coefficient of four size particles during the loading process is analyzed to explain the reason for the size effect of particle breakage.     

Numerical simulation of silencers

The problem of attenuating the noise from weapons firing is studied experimentally and numerically. As a possible method of attenuating the noise significantly, a silencer with no internal baffles is attached to the M242 cannon. The internal pressures inside the muffler are measured. The near-field overpressures outside the muffler at various polar angles are also measured. A numerical simulation of the flow through the muffler is performed, using Harten's shock-capturing method to solve the Euler equations of ideal compressible flow. The numerical simulation yields a detailed picture of the flow field as displayed by the pressure and Mach contours. Pressure–time curves at selected locations are obtained and compared with experimental data. There is good agreement, except that the numerical simulation generates more vigorous oscillations.     

Attenuation of shock waves by granular filters

Attenuation of shock waves in granular filters has been studied. Both pressurized air and solid explosives have been used for generating shock waves in a shock tube. The shock tube had a total length of m, and an internal diameter of 355 mm. Two large scale experiments have also been carried out in a tunnel with a cross-sectional area of 6.5 m². The results are compared with results found in the literature (Zloch, 1976; Medvedev et al., 1990; Britan et al., 2001) and previous experiments in a smaller scale by Slungaard (2002). A simple correlation based on the work of Zloch (1976) for shock wave attenuation in tube bundles and an extensive amount of experiments, is proposed. The correlation can be used to estimate the attenuation of the shock wave through a granular filter with filter characteristic . Setting B=6 will give a conservative estimate of the attenuation, while setting B=3 will give the best fit to all the results from this study and the results found in the literature. The correlation is independent of the type of driver (pressurised air or solid explosives) and upstream shock strength.Received: 30 September 2002, Accepted: 19 December 2002, Published online: 6 March 2003     

Numerical simulation of expansion/compression effect on shock induced turbulent flow

The interaction of homogeneous and isotropic turbulence with a shock wave is observed by solving the Reynolds-averaged Navier–Stokes equations with the k–? turbulence model. All turbulent fluctuations are measured at the period of expansion in the turbulent field and during compression by the reflected shock on turbulent field, and it is observed that the longitudinal turbulent velocity fluctuation is enhanced more at the period of expansion due to incident shock wave movement far from the turbulent field. The amplification of the turbulent kinetic energy (TKE) level in the shock/turbulence interaction depends on the shock wave strength and the longitudinal velocity difference across the shock wave. On decreasing the longitudinal velocity difference across the shock, the turbulent kinetic energy (TKE) level is less amplified. The TKE level is amplified by the factor of 1.5–1.8 in the shock/turbulence interaction where the dissipation rate of TKE decreases in all cases of shock/turbulence interaction. After the shock/turbulence interaction, the turbulent dissipative-length scale is amplified slightly and the amplification of the length scales decreases when increasing the shock strength. To cite this article: M.A. Jinnah, K. Takayama, C. R. Mecanique 333 (2005).     

Assessing continuum postulates in simulations of granular flow

Continuum mechanics relies on the fundamental notion of a mesoscopic volume “element” in which properties averaged over discrete particles obey deterministic relationships. Recent work on granular materials suggests that a continuum law may be inapplicable, revealing inhomogeneities at the particle level, such as force chains and slow cage breaking. Here, we analyze large-scale three-dimensional discrete-element method (DEM) simulations of different granular flows and show that an approximate “granular element” defined at the scale of observed dynamical correlations (roughly three to five particle diameters) has a reasonable continuum interpretation. By viewing all the simulations as an ensemble of granular elements which deform and move with the flow, we can track material evolution at a local level. Our results confirm some of the hypotheses of classical plasticity theory while contradicting others and suggest a subtle physical picture of granular failure, combining liquid-like dependence on deformation rate and solid-like dependence on strain. Our computational methods and results can be used to guide the development of more realistic continuum models, based on observed local relationships between average variables.     

Numerical investigation of transient nozzle flow

Abstract. The starting process of two-dimensional and axisymmetric nozzle flows has been investigated numerically. Special attention has been paid to the early phase of the starting process and to the appearance of a strong secondary shock wave. For both cases, shock intensities and velocities are obtained and discussed. The flow evolution in the axisymmetric case is proved to be more complex and the transient starting process is slower than in the plane case. Finally, the effects of changing the nozzle angle and the incident shock wave Mach number on the transient flow are addressed. It is shown that a faster start-up can be induced either by decreasing the nozzle angle or increasing the Mach number of the incident shock wave. Received 16 November 2001 / Accepted 24 September 2002 / Published online 4 December 2002 Correspondence to:A.-S. Mouronval (e-mail: mouronv@coria.fr)     

垂直上升气液环状流流动参数的数值计算

提出了一个新的气核-液膜耦合模型来求解垂直上升气液环状流在充分发展段的流动参数.本模型考虑了液膜、气核以及它们之间的相互影响和作用.模型中基本的气核区域和液膜区域的质量和动量方程由Fluent6.3.26进行求解,而液滴方程以及相界面上的夹带和沉积作用通过用户自定义接口函数UDF(User Defined Functi...     

Combining X-ray microtomography with computer simulation for analysis of granular and porous materials

The use of X-ray microtomographic (XMT) methods in analysing particulate systems has expanded rapidly in recent years with the availability of affordable desk-top apparatus. This review presents a summary of the major applications in which computer simulations are explicitly coupled with XMT in the area of granular and porous materials. We envisage two main ways of establishing the coupling between both techniques, based on the transference or exchange of information by using physical or geometrical paramet...     

NUMERICAL SIMULATION OF WING-BODY JUNCTION TURBULENCE FLOW

Wing-body junction turbulence flow is simulated by using RANS equation and boundary fitted coordinate technique.Three order differential scheme is used in the computation of convection term and two layers turbulence model are employed in the calculation.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured.The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

DEM simulation of particle mixing in a sheared granular flow

Mixing behaviors of particles are simulated in a sheared granular flow using differently colored but otherwise identical glass spheres, with five different bottom wall velocities. By DEM simulation, the solid fractions, velocities, velocity fluctuations and granular temperatures are measured. The mixing layer thicknesses are compared with the calculations from a simple diffusion equation using the data of apparent self-diffusion coefficients obtained from the current simulation measurements. The calculations and simulation results showed good agreements, demonstrating that the mixing process of granular materials occurred through the diffusion mechanism.     

SPH-based simulation of granular collapse on an inclined bed

The paper presents a numerical model for simulating a granular flow and its deposition on an inclined bed. A granular material is described as an elastic–plastic continuum and its constitutive law, namely Hooke's law, is discretized on the basis of the Smoothed Particle Hydrodynamics (SPH) method. In the equation of motion, however, the artificial viscosity, which is widely used in SPH, is not applied. The diffusive term derived from Hooke's law is introduced with a diffusion coefficient that varies depending on the stress and strain rate based on the Drucker–Prager yield function. The model is verified and validated through two numerical tests. It is shown that the basic elastic–perfectly plastic characteristics are reproduced with a simple shearing test. The effects of the diffusion coefficient and spatial resolution are investigated to show the validity of the model. In the simulation of the gravitational collapse of a granular column on an inclined bed, the performance of the model from the final deposition profile, the time history of the front position of the granular flow, the maximum runout distance, and the velocity profile are investigated for several cases of basal inclinations. The calculated results show good agreement with the experimental results.     

Computation of compaction in compressible granular material

Equations modeling compaction in a mixture of granular high explosive and interstitial gas are solved numerically. Both phases are modeled as compressible, viscous fluids. This overcomes well known difficulties associated with computing shock jumps in the inviscid version of the equations, which cannot be posed in a fully conservative form. One-dimensional shock tube and piston-driven compaction solutions compare favorably with experiment and known analytic solutions. A simple two-dimensional extension is presented.     

NUMERICAL SIMULATION OF COMPLICATED PIPE FLOW WITH k-εMODEL

In this paper, turbulence in a complicated pipe is simulated by using the k-ε model. The ladder-like mesh approximation is used to solve the problem of complicated boundary with the result of numerical simulation favorable. Two computational examples are given to validate the strong adaptability and stability of k-ε model.     

The influence of the drag force due to the interstitial gas on granular flows down a chute

Fully-developed steady flow of granular material down an inclined chute has been a subject of much research interest, but the effect of the interstitial gas has usually been ignored. In this paper, new expressions for the drag force and energy dissipation caused by the interstitial gas (ignoring the turbulent fluctuations of the gas phase) are derived and used to modify the governing equations derived from the kinetic theory approach for granular–gas mixture flows, where particles are relatively massive so that velocity fluctuations are caused by collisions rather than the gas flow. This new model is applied to fully-developed, steady mixture flows down an inclined chute and the results are compared with other simulations. Our results show that the effect of the interstitial gas plays a significant role in modifying the characteristics of fully developed flow. Although the effect of the interstitial gas is less pronounced for large particles than small ones, the flowfields with large particles are still very different from granular flows which do not incorporate any interactions with the interstitial gas.