Comparisons of some difference forms for compressible flow in cylindrical geometry on arbitrary Lagrangian and Eulerian framework

The study of cylindrically symmetric compressible fluid is interesting from both theoretical and numerical points of view. In this paper, the typical spherical symmetry properties of the numerical schemes are discussed, and an area weighted scheme is extended from a Lagrangian method to an arbitrary Lagrangian and Eulerian (ALE) method. Numerical results are presented to compare three discrete configurations, i.e., the control volume scheme, the area weighted scheme, and the plane scheme with the addition of a geometrical source. The fact that the singularity arises from the geometrical source term in the plane scheme is illustrated. A suggestion for choosing the discrete formulation is given when the strong shock wave problems are simulated.     

Dynamic crystal plasticity: An Eulerian approach

In this paper an Eulerian rate-dependent single crystal model that accounts for high-strain rates, large strains and rotations is developed. The viscoplastic law as well as the evolution equations for the lattice are written in terms of vectorial and tensorial quantities associated with the current configuration. The viscoplastic law is obtained from Schmid law using an overstress approach. Such an expression for the viscoplastic law is motivated by the microdynamics of crystal defects. A general analysis of the plane-strain response of the proposed rigid-viscoplastic single crystal model is presented. It is shown that only one differential equation, involving the orientation of one composite in-plane slip system, is necessary to describe the lattice evolution. Several two-dimensional boundary value problems, such as equal-channel die extrusion and channel die compression are selected to illustrate the predictive capabilities of the model. The results show that even at relatively low strain rates the viscosity plays an important role in the development of localized deformation modes. At high crosshead velocity, the plastic properties and crystal anisotropy are less important while inertia effects are dominant. Finally, the grains interaction is investigated by analyzing the compression of a grains multicrystal.     

On the theory of the shear-induced isotropic-to-nematic phase transition of side chain liquid-crystalline polymers

The shear-induced isotropic-to-nematic phase transition of side chain liquid-crystalline polymers is studied theoretically. A modification of the previous models of main-chain liquid crystals to the case of side chain liquid-crystalline polymers is proposed. Orientational and rheological properties of the model are studied in plane-shear flow. Predictions of the present model agree qualitatively with experimental results (Pujolle-Robic, Noirez in Nature 409:167, 2001).     

Apparent viscosity of a mixture of a Newtonian fluid and interacting particles

We investigate the behavior of fluid–particle mixtures subject to shear stress, by mean of direct simulation. This approach is meant to give some hints to explain the influence of interacting red cells on the global behavior of the blood. We concentrate on the apparent viscosity, which we define as a macroscopic quantity which characterizes the resistance of a mixture against externally imposed shear motion. Our main purpose is to explain the non-monotonous variations of this apparent viscosity when a mixture of fluid and interacting particles is submitted to shear stress during a certain time interval. Our analysis of these variations is based on preliminary theoretical remarks, and some computations for some well-chosen static configurations. To cite this article: A. Lefebvre, B. Maury, C. R. Mecanique 333 (2005).     

新型整体-局部高阶剪切变形理论及层合板精化三角形单元

基于整体-局部位移方法,建立了一种高阶剪切变形理论。整体位移部分采用的是Reddy理论的位移模式（1984）,局部位移为LIXY等（1997）建立的1,2-3理论的局部函数。这一理论使满足自由表面条件的Red@理论进一步满足层间位移、应力连续,同时有效减少了1,2—3理论的未知数个数。基于此理论深入开展了有限元法研究,建立了满足C^1连续条件的精化三节点三角形单元（每个节点参数为9个）。计算结果表明：建立的精化单元能准确计算整体位移和层间应力。     

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.     

Linear relationship between the first invariants of the stress and strain tensors in theories of plasticity with strain hardening

Energy-coupled stress and strain measures are defined in Euler coordinates. They are used to analyze the relationship between the first invariants of the stress and strain tensors for linearity and to determine strains at which the plastic component of the first strain invariant can be neglected. It is established that this relationship remains linear within an engineering plastic-strain tolerance of 0.2% irrespective of the value of strain intensity, which depends on the type of material and its stress state __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 3, pp. 60–72, March 2007.     

The Structure Functions of the Velocity and Temperature Fields from the Perspective of Dimensional Scaling

In this paper, dimensional scaling is used to describe the turbulence structure of the velocity and temperature fields in the inertial range and the far dissipation range as well as the intermediate transition range under locally isotropic conditions at sufficiently large Reynolds numbers. This kind of scaling is expressed in a strictly mathematical manner employing dimensional π -invariants analysis. It is shown that in the case of the asymptotic solutions for either the inertial range or the far dissipation range only one π number occurs that has to be considered as a non-dimensional universal constant. This π number may be determined theoretically or/and empirically. In the case of the transition range two π numbers occur. Consequently, a universal function is established that has to be derived theoretically or/and empirically, too. Here, Batchelor's [7] classical interpolation formula for the turbulence structure of the velocity field and the empirical one of Stolovitzky et al. [59], both may serve as universal functions, are compared with the results provided by numerical solutions of Kolmogorov′s [32] structure equation for the velocity field. It is shown that these interpolation formulae match not only the asymptotic solutions of the inertial range and the far dissipation range, respectively, but also these numerical results in an excellent manner. The former may be considered as necessary condition and the latter as sufficient condition. In the case of the temperature field results of the corresponding universal function are predicted using Yaglom's [63] structure equation. These results also match the corresponding asymptotic solutions of both the inertial range and the far dissipation range. However, in contrast to the case of the velocity field, the predicted universal function for the temperature field may notably overshoot its asymptotic solution for the inertial range. This overshooting occurs in the transition range and may be considered as an analogue to the so-called Hill ‘bump’ that usually occurs in the high-wave number portion of the temperature spectrum.     

Study of the aerodynamic characteristics of cone-sphere models in a hypersonic flow

The aerodynamic characteristics of cone-sphere models are studied at Mach numbers M = 6, 8.4, and 12 to 13 over a wide Reynolds number range. Models of a braking device (sphere) were connected with a load (frustum of a cone) by means of shrouds. The dependences of the aerodynamic coefficients C _x and C _y on the angle of attack α were obtained for different relative dimensions of the load and the braking device, shroud lengths, and Mach and Reynolds numbers. The effect of the above-mentioned parameters on the aerodynamic characteristics of the models is analyzed. The C _x (Re_D) dependences of load-parachutemodels in a symmetric flow are determined over the wide Mach and Reynolds number ranges 6 ≤ M ≤ 13 and 3 · 10³ ≤ Re_D ≤ 3 · 10⁶.