过冷水滴撞击螺旋桨桨叶表面的数值模拟

本文介绍了某型螺旋桨在不同飞行状态下.求解桨叶水滴撞击特性的数值方法.该方法对桨叶运动模型进行简化,并在对绕桨叶运动的气流场计算的基础上,采用拉格朗日方法求解气流场中水滴运动方程,得到水滴运动轨迹.进而,确定了水滴对桨叶的撞击特性参数,为桨叶防冰系统设计提供条件.主要结论如下:(1)在巡航状态下,桨叶沿展向方向上总收集系数Em和局部收集系数β不断增大;(2)在爬升状态下,随着爬升高度H不断增大,飞行速度V0不断增大,水滴撞击在桨叶表面的范围有所增加,而且β随之增大;(3)随着水滴平均有效直径(MVD)的增大,水滴撞击在桨叶表面的范围明显增加同时,β在桨叶表面同一位置的值也随之增大.     

A finite element method study of Eulerian droplets impingement models

To compute droplet impingement on airfoils, an Eulerian model for air flows containing water droplets is proposed as an alternative to the traditional Lagrangian particle tracking approach. Appropriate boundary conditions are presented for the droplets equations, with a stability analysis of the solution near the airfoil surface. Several finite element formulations are proposed to solve the droplets equations, based on conservative and non‐conservative forms and using different stabilization terms. Numerical results on single and multi‐elements airfoils for droplets of mean volume diameter, as well as for a Langmuir distribution of diameters, are presented and validated against measured values. Copyright © 1999 John Wiley & Sons, Ltd.     

A dynamic eddy-viscosity closure model for large eddy simulations of two-dimensional decaying turbulence

This paper puts forth a simplified dynamic eddy-viscosity subgrid-scale model for the vorticity transport equation which is employed in a large eddy simulation study of freely evolving isotropic two-dimensional turbulent flows. The dynamic parameter is averaged in space, thereby retrieving a spatially constant value which only varies in time. The proposed dynamic model is applied to a two-dimensional decaying turbulence problem in a square periodic box, which is a standard prototype of more realistic turbulent flows in the atmosphere and oceans, in order to eliminate any possible errors associated with the boundary conditions or mesh non-uniformities. Compared with high-resolution direct numerical simulations, the performance of the dynamic model is systematically investigated considering various filtering strategies by means of test filters. The effects of the computational resolution and the filtering ratio between the test and the grid filters are also studied by using a huge set of parameters.     

Modelling impingement of hollow metal droplets onto a flat surface

Despite many theoretical and experimental works dealing with the impact of dense melt droplets on the substrate during the process of thermal spray coating, the dynamics of the impingement of hollow melt droplet and the subsequent splat formation are not well addressed. In this paper a model study for the dynamic impingement of hollow droplet is presented. The hollow droplet is modelled such that it consists of a liquid shell enclosing a gas cavity. The impingement model considers the transient flow dynamics during impact, spreading and solidification of the droplet using the volume of fluid surface tracking method (VOF) coupled with a solidification model within a one-domain continuum formulation. The results for spreading, solidification and formation of splats clearly show that the impingement process of hollow droplet is distinctly different from the dense droplet. Study with different droplet void fractions and void distribution indicates that void fraction and void distribution have a significant influence on the flow dynamics during impact and on the final splat shape. The results are likely to provide insights for the less-explored behaviour of hollow melt droplets in thermal spray coating processes.     

A two-dimensional model for linear elastic thick shells

In this paper, a two-dimensional model for linear elastic thick shells is deduced from the three-dimensional problem of a shell thickness 2ε, ε > 0. From different scalings on the tangent and normal components of the displacement u^ε as widely used in recent works, the limit displacement appears to be Kirchhoff–Love displacement of a different type. It contains additional terms to those found in the Reissner–Mindlin model and satisfies more general equations containing the classical terms found in the literature and some new terms related to the third fundamental form. Such terms could not be well handled in the usual framework. Shear stresses across the thickness are also computed. This model appears to be appropriate to handle stiffened shells which, in fact, cannot be considered uniformly as shallow shells. As a by-product it also lays the mathematical background to justify the Reissner–Mindlin model for plates and will probably contribute to a better understanding of the locking phenomenon encountered in computational mechanics.     

A two-equation turbulence model for jet flows laden with vaporizing droplets

A two-equation turbulence model for steady incompressible two-phase flows including phase change has been recently developed by Mostafa & Elghobashi (1984). This model is tested for the flow of a turbulent axisymmetric gaseous jet laden with evaporating liquid droplets. To avoid the problem of density fluctuations of the carrier phase at this stage, only isothermal flow is considered and vaporization is assumed to be due to the vapor concentration gradient. The continuous size distribution of the droplets is approximated by finite size groups. Each group is considered as a continuous phase interpenetrating and interacting with the carrier phase. Two test cases have been predicted by the model. The first is for a Freon-11 spray issuing from a round nozzle, where experimental data are available at distances equal to or greater than 170 nozzle diameters. Good agreement between the data and the predictions was achieved. The second is for a methanol spray where no experiments are available yet and the predictions consider the flow region close to the nozzle ($\text{z/D < 40}$). The results of the methanol spray include distributions of the mean velocity, volume fractions of the different phases, concentration of the evaporated material in the carrier phase, turbulence intensity and shear stress of the carrier phase, droplet diameter distribution, and the jet spreading rate. In this case the results are analyzed based on a qualitative comparison with the corresponding single phase jet flow.     

Experimental investigation on splashing and nonlinear fingerlike instability of large water drops

The fluid physics of the splashing and spreading of a large-scale water drop is experimentally observed and investigated. New phenomena of drop impact that differ from the conventional Rayleigh–Taylor instability theory are reported. Our experimental data shows good agreement with previous work at low Weber number but the number of fingers or instabilities begins to deviate from the R–T equation of Allen at high Weber numbers. Also observed were multiple waves (or rings) on the spreading liquid surface induced from pressure bouncing (or pulsation) within the impacting liquid. The first ring is transformed into a radially ejecting spray whose initial speed is accelerated to a velocity of 4–5 times that of the impacting drop. This first ring is said to be “splashing,” and its structure is somewhat chaotic and turbulent, similar to a columnar liquid jet surrounded by neighboring gas jets at relatively high impact speed. At lower impact speeds, splashing occurs as a crown-shaped cylindrical sheet. A second spreading ring is observed that transforms into fingers in the circumferential direction during spreading. At higher Weber number, the spreading of a third ring follows that of the second. This third ring, induced by the pressure pulsation, overruns and has fewer fingers than the second, which is still in a transitional spreading stage. Several important relationships between the drop impact speed, the spray ejection speed of the first ring, and the number of fingers of the second and third rings are presented, based on data acquired during a set of drop impact experiments. Issues related to the traditional use of the R–T instability are also addressed.     

A non-local two-dimensional foundation model

Classical foundation models such as the Pasternak and the Reissner models have been recently reformulated within the framework of non-local mechanics, by using the gradient theory of elasticity. To contribute to the research effort in this field, this paper presents a two-dimensional foundation model built by using a mechanically based non-local elasticity theory, recently proposed by the authors. The foundation is thought of as an ensemble of soil column elements resting on an elastic base. It is assumed that each column element is acted upon by a local Winkler-like reaction force exerted by the elastic base, by contact shear forces and volume forces due, respectively, to adjacent and non-adjacent column elements. As in the Pasternak model, the contact shear forces involve the second-order derivative of the column element displacement. The volume forces are non-local forces assumed to depend (1) on the relative displacement between the interacting column elements through power-law distance-decaying attenuation functions and (2) on the product between the volumes of the interacting column elements. As a result, the equilibrium equations are fractional differential equations, for which a numerical solution can be readily found based on the finite difference method. Solutions are built for different foundation shapes and loading conditions.     

A consistent finite element model for the two-dimensional continuum

Summary The finite element approximation to the continuum problem is examined from the viewpoint of the principle of virtual work. It is shown that the usual nodal equilibrium equations for triangular elements are a consistent consequence of a piecewise constant strain field, thus guaranteeing that many results of general continuum theory can be directly applied to the finite element model, and also clarifying the relation between the two models.

---

Übersicht Das Verfahren, ein Kontinuum durch finite Elemente anzunähern wird vom Standpunkt des Prinzips der virtuellen Arbeiten untersucht. Es wird gezeigt, daß die üblichen Knotenpunktsgleichungen für dreieckförmige Elemente eine Folge des stückweise konstanten Verformungsfeldes sind. Auf diese Weise wird sichergestellt, daß viele Ergebnisse der allgemeinen Kontinuumstheorie unmittelbar auf das aus endlichen Elementen aufgebaute Modell übertragen werden können. Gleichzeitig werden die Beziehungen zwischen beiden Modellen geklärt.

---

---

Dedicated to Professor Dr. H. Ziegler on the occasion of his 60th birthday.

---

This research was sponsored by the National Science Foundation, Grant GK 10549.     

A moment method for splashing and evaporation processes of polydisperse sprays

A new moment method for the modelling of polydisperse sprays is proposed that simultaneously takes into account the dispersion in droplet size and droplet velocity. For the derivation of this Eulerian method the kinetic spray equation is used which constitutes a partial differential equation for the probability density function of droplets. To reduce the complex kinetic spray equation to a form that can be managed with the available numerical procedures, moment transforms with respect to the droplet velocity and the droplet size are conducted. The resulting moment equations are closed by choosing an approximate probability density function which applies to polydisperse sprays. The method is successfully tested for configurations in which a polydisperse spray is either splashed, evaporated or effected by a Stokes drag force. The tests are organised in such a way that crossing of two spray distributions is always included. The new method is able to capture the polydisperse nature of sprays as well as the bi-(or multi-) modal character of the droplet velocity distribution function, for example, when droplets cross each other.     

A rheological model with constant approximate volume for immiscible blends of ellipsoidal droplets

We develop a set of evolution equations describing the effects of a general deformation field on the shape, size, and orientation of constant-volume droplets suspended in a Newtonian fluid. The rheological characteristic functions of this incompressible and immiscible polymer blend model are also derived in tandem with the abovementioned set. The constant-volume constraint is implemented using a recent methodology (Edwards BJ, Dressler M, Grmela M, Ait-Kadi A (2002) Rheological models with microstructural constraints. Rheo Acta (in press)) and discussed relative to the similarly volume-constrained model of Almusallam et al. (Almusallam AS, Larson RG, Solomon MJ (2000) A constitutive model for the prediction of ellipsoidal droplet shapes and stresses in immiscible blends. J Rheol 44:1055–1083). Sample results are reported for step-strain and shear relaxation profiles.     

A rate-dependent two-dimensional free energy model for ferroelectric single crystals

The one-dimensional free energy model for ferroelectric materials developed by Smith et al. [29–31] is generalized to two dimensions. The two-dimensional free energy potential proposed in this paper consists of four energy wells that correspond to four variants of the material. The wells are separated by four saddle points, representing the barriers for 90°-switching processes, and a local maximum, across which 180°-switching processes take place. The free energy potential is combined with evolution equations for the variant fractions based on the theory of thermally activated processes. The model is compared to recent measurements on BaTiO₃ single crystals by Burcsu et al. [8], and predicitions are made concerning the response to the application of in-plane multi-axial electric fields at various frequencies and loading directions. The kinetics of the 90°- and 180°-switching processes are discussed in detail.     

Lattice Boltzmann simulation of multiple droplets impingement and coalescence in an inkjet-printed line patterning process

Three-dimensional computations on the basis of the index-function lattice Boltzmann method are performed to simulate the process of multiple droplets impinging and coalescing into a line pattern on a solid substrate. The employed calculation model is validated by theoretical calculated values and experimental data from the literature. The influences of the equilibrium contact angle, droplet spacing and impinging velocity on the droplets impingement and coalescence behaviours are investigated. Numerical results demonstrate the width of the formed line depends significantly on the equilibrium contact angle and droplet spacing. The droplet spacing plays a significant role in controlling the coalescence moment of multiple droplets. The resolution of the printed pattern can be slightly increased with increase in impinging velocity.     

Nearly singular two-dimensional Kolmogorov flows for large Reynolds numbers

We study the convergence of two-dimensional stationary Kolmogorov flows as the Reynolds number increases to infinity. Since the flows considered are stationary solutions of Navier-Stokes equations, they are smooth whatever the Reynolds number may be. However, in the limit of an infinite Reynolds number, they can, at least theoretically, converge to a nonsmooth function. Through numerical experiments, we show that, under a certain condition, some smooth solutions of the Navier-Stokes equations converge to a nonsmooth solution of the Euler equations and develop internal layers. Therefore the Navier-Stokes flows are nearly singular for large Reynolds numbers. In view of this nearly singular solution, we propose a possible scenario of turbulence, which is of an intermediate nature between Leray's and Ruelle-Taken's scenarios.     

Flow characteristics of slot jet impingement on a wedge

An experimental investigation is made to study the flow characteristics of slot jet impingement on a wedge whose included angle is 90 degrees. The aim of this investigation is to study the characteristics of the flow field near the wedge surface for various parameters. The different parameters like, jet velocity, slot width, distance of wedge vertex from the jet exit and the inclination of the wedge to jet axis are systematically varied to see their effect on the flow field. The flow field near the wedge vertex is similar to stagnation point flow. Far away from the vertex, the flow field is like that of wall jet. Near the vertex, very large variations of static pressure are observed in streamwise and transverse directions. This is due to large streamwise curvature and stagnation of flow. The transverse pressure gradient slowly decays in the streamwise direction, as a result, the velocity profiles are different from the similarity profiles of stagnation point flow and wall jet in the respective regions. Experiments are conducted for slot widths of 10 mm, 15 mm and 30 mm each for the distance between slot and wedge vertex of 80 mm, 120 mm and 240 mm. The static pressure and velocity profiles are measured by calibrated disk type static pressure probe and pitot tube respectively at various streamwise locations.Es wurde eine experimentelle Studie über die Strömmungscharakteristiken eines auf einem rechtwinkligen Keil auftretenden Spaltstrahls durchgeführt. Das Ziel dieser Untersuchungen ist die Studie des Strömungsverhaltens in Keilnähe, in Abhängigkeit verschiedener Parameter wie Strahlgeschwindigkeit, Spaltbreite, Entfernung der Keilspitze und Winkel zwischen Strahlachse und Keil. Diese Parameter werden systematisch variiert um deren Einfluß auf das Strömungsverhalten zu bestimmen. Die Strömung an der Keilspitze ist ähnlich der Staupunktströmung. Weit hinter der Spitze gleicht das Strömungsfeld dem eines auf einer Wand auftreffenden Strahls. In Nähe der Keilspitze sind die Änderungen des statischen Drucks in Strömungsrichtung und quer zur Strömung groß. Grund dafür sind die starke Krümmung der Strömung und die Stagnation der Strömung. Der Druckgradient in Querrichtung nimmt langsam in Strömungsrichtung ab, daher unterscheiden sich die Geschwindigkeitsprofile von den Ähnlichkeitsprofilen der Staupunktströmung und des Wandstrahls in den jeweiligen Bereichen. Die Versuche wurden für Spaltbreiten von 10, 15 und 30 mm und Keilentfernungen von 80, 120 und 240 mm durchgeführt. Der statische Druck und die Geschwindigkeitsprofile wurden mit kalibrierten scheibenförmigen Drucksonden bzw. einer Pitot-Sonde an verschiedenen Orten gemessen.     

A new two-dimensional model for electro-mechanical response of thick laminated piezoelectric actuator

This paper investigates the electro-mechanical behaviour of a thick, laminated actuator with piezoelectric and isotropic lamina under externally applied electric loading using a new two-dimensional computational model. The elastic core is relatively thick and thus it is modelled by Timoshenko thick-beam theory. Although the piezoelectric lamina is a beam-like layer, it is formulated via a two-dimensional model because of not only the strong electro-mechanical coupling, but also of the presence of a two-dimensional electric field. It is shown in this paper that a one-dimensional model for the piezoelectric beam-like layer is inadequate. The piezoelectric model is constructed within the scope of linear piezoelectricity. The actuation response is induced through the application of external electric voltage. Under the strong coupling of elasticity and electricity, the strain energy and work of electric potential are presented. The electro-mechanical response of the laminated Timoshenko beam is formulated and determined via a variational energy principle. Numerical examples presented illustrate convincing comparison with finite element solutions and existing published data. New numerical solutions are also presented to investigate the geometric effect on the electro-mechanical bending behaviour.     

A constitutive model for orthotropic elasto-plasticity at large strains

Summary The paper presents a thermodynamically consistent constitutive model for elasto-plastic analysis of orthotropic materials at large strain. The elastic and plastic anisotropies are assumed to be persistent in the material but the anisotropy axes can undergo a rigid rotation due to large plastic deformations. The orthotropic yield function is formulated in terms of the generally nonsymmetric Mandel stress tensor such that its skew-symmetric part is additionally taken into account. Special attention is focused on the convexity of the yield surface resulting in the nine-dimensional stress space. Of particular interest are new convexity conditions which do not appear in the classical theory of anisotropic plasticity. They impose additional constraints on the material constants governing the plastic spin. The role of the plastic spin is further studied in simple shear accompanied by large elastic and large plastic deformations. If the plastic spin is neglected, the shear stress response is characterized by oscillations with an amplitude strictly dependent on the degree of the plastic anisotropy.accepted for publication 2 March 2004     

A hybrid dynamic Smagorinsky model for large eddy simulation

A hybrid dynamic subgrid-scale model (HDSM) pertaining to Large-eddy simulation (LES) has been developed. The coefficient obtained by German dynamic Smagorinsky model (DSM) was integrated with a new dynamic coefficient, based on the dynamic subgrid characteristic length and controlled by the subgrid-scale (SGS) motions. In HDSM, the characteristic wave number determining the characteristic length of the dynamic subgrid is calculated from a new energy weighted mean method when the subgrid scale turbulent kinetic energy and the dissipation wave number are known. The dissipation wave-number is derived from the SGS turbulent kinetic energy spectrum equation. The total dissipation rate spectrum equation is based on the Pao energy spectrum and local equilibrium assumption. The dynamic subgrid characteristic length could take into account the rapidly fluctuating small scale behaviours and the spatial variation of turbulent characteristics. HDSM was used to simulate the fully developed channel and turbulent flow past a circular cylinder, and to determine the impact of the dam-break flow on downstream structure. The HDSM is robust in respect to anisotropic mesh and is less sensitive to grid resolution, and would accurately describe the energy transfer from large-scale to SGS fluctuations and capture more fluctuations of turbulence with same meshes compared to the DSM.     

Flow characteristics of spray impingement in PFI injection systems

The present paper addresses an experimental study of the dynamic exchanges between the impact of an intermittent spray and the liquid film formed over the target, based on detailed phase Doppler anemometry (PDA) measurements of droplet size, velocity and volume flux in the vicinity of the impact. The flow configuration is that of a pulsed injector spraying gasoline onto a flat disc to simulate the port fuel injection (PFI) of an internal combustion engine operating at cold-start conditions. The measurements evidence that the outcome of impact cannot be accurately predicted based on the characteristics of the free spray, but requires precise knowledge of the flow structure, induced by the target. The implications for spray–wall interaction modelling are then discussed based on the application of conservation equations to the mass, momentum and energy exchanged between the impinging droplets and the liquid film. The results show that the liquid film starts to form in the vicinity of the stagnation region at early stages of injection and a non-negligible proportion of droplets impinging at outer regions splash after interaction with the film. Film disruption is mainly driven by the intermittent axial momentum of impinging droplets, which enhances the vertical oscillations. The radial momentum imparted to the liquid film at the stagnation region is fed back onto secondary droplets emerging later during the injection cycle at outwards locations, where momentum of impacting droplets is much smaller. As a consequence, although the number of splashed droplets is enhanced by normal momentum, their size and ejection velocity depends more on the radial spread induced onto the liquid film and, hence, on the radial momentum at impact. The analysis further shows that existing spray–wall interaction models can be improved if the dynamic exchanges between the impacting spray and the liquid film are accounted.     

A two-dimensional,higher-order,elasticity-based micromechanics model

A new, two-dimensional (2D) homogenization theory is proposed. The theory utilizes a higher-order, elasticity-based cell model (ECM) analysis. The material microstructure is modeled as a 2D periodic array of unit cells where each unit cell is discretized into four subregions (or subcells). The analysis utilizes a (truncated) eigenfunction expansion of up to fifth order for the displacement field in each subcell. The governing equations for the theory are developed by satisfying the pointwise governing equations of geometrically linear continuum mechanics exactly up through an order consistent with the order of the subcell displacement field. The formulation is carried out independently of any specified constitutive models for the behavior of the individual phases (in the sense that the general governing equations hold for any constitutive model). The fifth order theory is subsequently specialized to a third order theory. Additionally, the higher order analyzes reduce to a theory equivalent to the original 2D method of cells (MOC) theory when all higher order terms are eliminated. The proposed 2D theory is the companion theory to an equivalent 3D theory [T.O. Williams, A three-dimensional, higher-order, elasticity-based micromechanics model, Int. J. Solids Struc., in press].Comparison of the predicted bulk and local responses with published results indicates that the theory accurately predicts both types of responses. The high degree of agreement between the current theory results and published results is due to the correct incorporation of the coupling effects between the local fields.The proposed theory represents the necessary theoretical foundations for the development of exact homogenization solutions of generalized, two-dimensional microstructures.