Modeling of the thermoerosive action of two-phase media

The action of a supersonic two-phase flow on the surface of a solid is considered. Various methods of modeling the particle concentration, the dynamic head, and the temperature and velocity disequilibria of the phases in connection with the investigation of the thermoerosive action of atmospheric formations in gas dynamic installations are analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 92–98, July–August, 1986.     

Modeling the evolution of droplet size distribution in two-phase flows

A theoretical model is developed in the present study to simulate droplet motion and the evolution of droplet size distribution (DSD) in two-phase air/dispersed water spray flows. The model takes into account several processes which influence DSD and droplet trajectory: droplet collision and coalescence, evaporation and cooling, gravitational settling, and turbulent dispersion of dispersed phase. The DSDs determined by the model at different locations in a two-phase flow are evaluated by comparing them to experimental observations obtained in an icing wind tunnel. The satisfactory coincidence between simulation and experimental results proves that the model is reliable when modeling two-phase flows under icing conditions. The model is applied for two particular examples in which the modification of DSD is calculated in two-phase flows under conditions describing in-cloud icing and freezing drizzle.     

Existence of the detonation cellular structure in two-phase hybrid mixtures

The cellular detonation structure has been recorded for hybrid hydrogen/air/aluminium mixtures on 1.0 m 0.110 m soot plates. Addition of aluminium particles to the gaseous mixture changes its detonation velocity. For very fine particles and flakes, the detonation velocity is augmented and, in the same time, the cell width diminishes as compared with the characteristic cell size of the mixture without particles. On the contrary, for large particles, the detonation velocity decreases and the cell size becomes larger than . It appears that the correlation law between the cell size and the detonation velocity in the hybrid mixture is similar to the correlation between the cell size and the rate of detonation overdrive displayed for homogeneous gaseous mixtures. Moreover, this correlation law remains valid in hybrid mixtures for detonation velocities smaller than the value D of the mixture without particles. Received 10 May 2001 / Accepted 12 August 2002 Published online 19 December 2002 Correspondence to: B. Veyssiere (e-mail: veyssiere@lcd.ensma.fr)     

Numerical modelling of detonation structure in two-phase flows

A one-dimensional physical model and a numerical method for the simulation of heterogeneous detonation were proposed based on an Eulerian approach for heterogeneous flows. The combination of modern shock-capturing schemes in combination with a dynamically moving, adaptive grid ensure the properresolution of both reaction zones and flow discontinuities. Numerical examples illustrate the effect of the heat release due to heterogeneous combustion. Received August 4, 1995 / Accepted December 12, 1995     

Stagnation region of a blunt body in two-phase hypersonic flow

The fluid-mechanics equations of a two-velocity, two-temperature medium are used to investigate flow near the stagnation point of a blunt body washed by a hypersonic stream of gas containing solid or liquid deformed particles. The effect of particles of the gasdynamic flow parameters is analyzed. A relaxation layer was found to occur near the body, with marked changes in the gas parameters. It is shown that the presence of particles in the flow reduces the shock stand-off distance. The results of computations on the dynamics and heating of particles in the shock layer are discussed. A solution in finite form is obtained in the limiting case of fine particles by the method of asymptotic expansions. The motion of solid or liquid particles in hypersonic shock layers has been the subject of several papers [1–6], in which particle dynamics was examined, assuming that the particles have a negligible influence on the gasdynamic flow parameters. The solutions obtained are therefore limited to the case of low mass particle concentration in the incident flow. A numerical solution not subject to this limitation was obtained in [7] for supersonic two-phase flow over a wedge.     

Modeling of two-phase flow in porous media with heat generation

The main purpose of this work is investigation of coolability of a boiling debris bed. The main governing equations are derived using volume averaging technique. From this technique some specific interfacial areas between phases are appeared and proper relations for modeling these areas are proposed. Using these specific areas, a modification for the Tung/Dhir model in the annular flow regime is proposed. The proposed modification is validated and the agreements with experimental data are good. Finally, governing equations and relations are implemented in the THERMOUS program to model two-phase flow in the debris bed in the axisymmetric cylindrical coordinate. Two typical configurations including flat and mounted beds are considered and the main physical phenomena during boiling of water in the debris bed are studied. Comparing the results with the one-dimensional analysis shows higher specific power of the bed.     

二维细胞在剪切流中的运动特性

主要模拟二维细胞在剪切流中的运动特性.计算过程采用浸入边界法,将细胞模化成Navier-Stokes方程中的力源,而不是真实物体.假设细胞的初始形状为椭圆,细胞内外流体粘性相同,细胞膜的弹性力模型选用E-S模型.本文模拟四种不同真圆度情况下细胞的形变情况,观测到初始阶段细胞沿着长轴方向做拉伸和旋转运动,达到稳定状态后细胞作类坦克履带式运动;并且发现细胞达到稳定状态所需要的时间随真圆度的增加而增加,而细胞的稳态倾角随真圆度的增加而减少.     

Modeling the macroscopic behavior of two-phase nonlinear composites by infinite-rank laminates

A new approach is proposed for estimating the macroscopic behavior of two-phase nonlinear composites with random, particulate microstructures. The central idea is to model composites by sequentially laminated constructions of infinite rank whose macroscopic behavior can be determined exactly. The resulting estimates incorporate microstructural information up to the two-point correlation functions, and require the solution to a Hamilton–Jacobi equation with the inclusion concentration and the macroscopic fields playing the role of ‘time’ and ‘spatial’ variables, respectively. Because they are realizable, by construction, these estimates are guaranteed to be convex, to satisfy all pertinent bounds, to exhibit no duality gap, and to be exact to second order in the heterogeneity contrast. Sample results are provided for two- and three-dimensional power-law composites, and are compared with other homogenization estimates, as well as with numerical simulations available from the literature. The estimates are found to give physically sensible predictions for all the cases considered, even for extreme values of the nonlinearity and heterogeneity contrast. Interestingly, in the case of isotropic porous materials under hydrostatic loadings, the estimates agree exactly with standard Gurson-type models for viscoplastic porous media.     

Modeling of deformation induced anisotropy in free-end torsion

The main purpose of this work is to develop a phenomenological model, which accounts for the evolution of the elastic and plastic properties of fcc polycrystals due to a crystallographic texture development and predicts the axial effects in torsion experiments. The anisotropic portion of the effective elasticity tensor is modeled by a growth law. The flow rule depends on the anisotropic part of the elasticity tensor. The normalized anisotropic part of the effective elasticity tensor is equal to the 4th-order coefficient of a tensorial Fourier expansion of the crystal orientation distribution function. Hence, the evolution of elastic and viscoplastic properties is modeled by an evolution equation for the 4th-order moment tensor of the orientation distribution function of an aggregate of cubic crystals. It is shown that the model is able to predict the plastic anisotropy that leads to the monotonic and cyclic Swift effect. The predictions are compared to those of the Taylor–Lin polycrystal model and to experimental data. In contrast to other phenomenological models proposed in the literature, the present model predicts the axial effects even if the initial state of the material is isotropic.     

Modeling analysis of bubble flow regime in a closed two-phase thermosyphon

Predictions of the operating liquid level in the evaporator of a closed two-phase thermosyphon (gravity heat pipe) are given throughout a simplified analysis which takes the influence of the dimension and condensation heat transfer in the condenser of the heat pipe into account. In order to verify the accuracy of our model comparison of the present study with some published results is made by means of computational examples.     

Modeling the viscoplastic micromechanical response of two-phase materials using Fast Fourier Transforms

A viscoplastic approach using the Fast Fourier Transform (FFT) method for obtaining local mechanical response is utilized to study microstructure-property relationships in composite materials. Specifically, three-dimensional, two-phase digital materials containing isotropically coarsened particles surrounded by a matrix phase, generated through a Kinetic Monte Carlo Potts model for Ostwald ripening, are used as instantiations in order to calculate the stress and strain-rate fields under uniaxial tension. The effects of the morphology of the matrix phase, the volume fraction and the contiguity of particles, and the polycrystallinity of matrix phase, on the stress and strain-rate fields under uniaxial tension are examined. It is found that the first moments of the stress and strain-rate fields have a different dependence on the particle volume fraction and the particle contiguity from their second moments. The average stresses and average strain-rates of both phases and of the overall composite have rather simple relationships with the particle volume fraction whereas their standard deviations vary strongly, especially when the particle volume fraction is high, and the contiguity of particles has a noticeable effect on the mechanical response. It is also found that the shape of stress distribution in the BCC hard particle phase evolves as the volume fraction of particles in the composite varies, such that it agrees with the stress field in the BCC polycrystal as the volume of particles approaches unity. Finally, it is observed that the stress and strain-rate fields in the microstructures with a polycrystalline matrix are less sensitive to changes in volume fraction and contiguity of particles.     

Interfacial debonding of two-phase composite with periodic structure

A mixed analytical-numerical (boundary element method) procedure is presented for estimating the effective elastic moduli of a two-phase periodic composite by application of a unit cell. The two-phase composite consists of a metal/polymer matrix and one/three circular ceramic inclusions with adhesive and partial debonding of the interface. The results are displayed numerically with special attention given to development of plastic zones as debonding occurs. Dependence of load-time history is exhibited.     

Modeling of stratified two-phase flow in pipes,pumps and other devices

A model for stratified two-phase flow in pipes. pumps and other devices is presented. Using the assumption of a hydrostatic distribution of pressure over the cross section of a pipe. the effects of stratification are taken into account by means of specific terms in momentum and energy equations. These terms represent the pressure interactions between phases. Four concrete examples (classical earth gravity stratified flow. stratified flow in a regularly curved pipe. annular stratified flow caused by pre- or postrotation. and stratified flow in the impeller of a pump) and a comparison with an experiment are presented.     

Fluid structure interaction problems in large deformation

The present article deals with the simulation of fluid structure interaction problems in large deformation, and discusses two aspects of their numerical solution: (i) the derivation of energy conserving time integration schemes in presence of fluid structure coupling, moving grids, and nonlinear kinematic constraints such as incompressibility and contact, (ii) the introduction of adequate preconditioners efficiently chaining local fluid and structure solvers. Solutions are proposed, analyzed and tested using nonlinear energy correcting terms, and added mass based Dirichlet Neumann preconditioners. Numerical applications include nonlinear impact problems in elastodynamics and blood flows predictions within flexible arteries. To cite this article: P. Le Tallec et al., C. R. Mecanique 333 (2005).     

Modeling slow deformation of polygonal particles using DEM

We introduce two improvements in the numerical scheme to simulate collision and slow shearing of irregular particles. First, we propose an alternative approach based on simple relations to compute the frictional contact forces. The approach improves efficiency and accuracy of the Discrete Element Method （DEM） when modeling the dynamics of the granular packing. We determine the proper upper limit for the integration step in the standard numerical scheme using a wide range of material parameters. To this end, we study the kinetic energy decay in a stress controlled test between two particles. Second, we show that the usual way of defining the contact plane between two polygonal particles is, in general, not unique which leads to discontinuities in the direction of the contact plane while particles move. To solve this drawback, we introduce an accurate definition for the contact plane based on the shape of the overlap area between touching particles, which evolves continuously in time.     

Modeling slow deformation of polygonal particles using DEM

We introduce two improvements in the numerical scheme to simulate collision and slow shearing of irregular particles. First, we propose an alternative approach based on simple relations to compute the frictional contact forces. The approach improves efficiency and accuracy of the Discrete Element Method (DEM) when modeling the dynamics of the granular packing. We determine the proper upper limit for the integration step in the standard numerical scheme using a wide range of material parameters. To this end, we study the kinetic energy decay in a stress controlled test between two particles. Second, we show that the usual way of defining the contact plane between two polygonal particles is, in general, not unique which leads to discontinuities in the direction of the contact plane while particles move. To solve this drawback, we introduce an accurate definition for the contact plane based on the shape of the overlap area between touching particles, which evolves continuously in time.