THREE-DIMENSIONAL EXTRUDATE SWELL: FORMULATION WITH THE STREAM TUBE METHOD AND NUMERICAL RESULTS FOR A NEWTONIAN FLUID

Stream tube analysis, already applied to two-dimensional extrudate swell problems involving rate and integral constitutive equations for incompressible fluids, is now considered in the problem of free surface determination in a three-dimensional flow situation. The method allows computation of the unknown free surface by considering only a ‘peripheral stream tube’ limited by the wall and the jet surface and an inner stream surface. Those boundary surfaces are determined by considering the conservation equations together with boundary condition equations, solved by the Levenberg/Marquatdt optimization algorithm. The method leads to a considerable reduction in the number of degrees of freedom and the storage area. As in a previous study in the two-dimensional case, singularity problems in the vicinity of the junction points between the wall and the free surface are avoided. However, the numerical method still allows evaluati on of stress peaks due to the singularity at the exit, as may be observed for results obtained with a Newtonian fluid in a duct of square cross-section.     

Numerical study of flows of complex fluids between eccentric cylinders using transformation functions

In this paper, we investigate fluid flows between eccentric cylinders by means of two stream‐tube analyses. The first method considers a one‐to‐one global transformation function that allows the physical domain to be transformed into a mapped domain, used as computational domain, that involves concentric streamlines. The second approach uses local transformations and domain decomposition techniques to deal with mixed flow regimes. Both formulations are particularly adapted for handling time‐dependent constitutive equations, since particle‐tracking problems are avoided. Mass conservation is verified in both formulations and the relevant numerical procedure can be carried out using simple meshes built on the mapped streamlines. Fluids obeying anelastic and viscoelastic constitutive equations are considered in the calculations. The numerical results are consistent with those in the literature for the flow rates tested. Application of the method to the K‐BKZ memory‐integral constitutive equation highlights significant differences between the model predictions and those provided by more simple rheological models. Copyright © 2002 John Wiley & Sons, Ltd.     

Three-dimensional extrudate swell experimental and numerical study of a polyethylene melt obeying a memory-integral equation

The present work deals with the experimental and numerical features of the flow of a linear low-density polyethylene melt (LLDPE) at 160°C at the exit of a die of square cross-section. The rheological properties of the fluid are fitted by a Wagner's memory-integral constitutive equation. The characteristics of the extrudate jet are determined by optical means at different flow rates. The stream-tube analysis, already applied to two-dimensional extrudate swell problems involving rate and integral constitutive equations, is considered to simulate the flow field. The method avoids particle tracking problems related to integral models and allows computation of the unknown free surface by considering only a `peripheral stream tube' limited by the wall and the jet surface and an inner stream surface. Those boundary surfaces are determined by considering the conservation equations together with boundary condition equations, solved by the Levenberg–Marquardt optimization algorithm. The method leads to a considerable reduction in the number of degrees of freedom and the storage area. The numerical results are found to be generally consistent with the experimental data and highlight the growing importance of stress peaks due to the singularity at the exit when the flow rate increases.     

Stress distribution in multilayered composites near loaded edges

The edge effect in layered composite material is studied using the piecewise-homogeneous body model and the exact equations of the theory of elasticity. It is assumed that continuously distributed normal forces act at the edges of the reinforcing layers. A plain strain state is considered and the stresses are expressed in terms of the solutions of a system of dual singular integral equations. The singularity of the stresses is determined by the solution procedure. The concentration of the reinforcing layers is assumed low and the interaction between them is not taken into account. A numerical algorithm is developed and numerical results on the stress distribution are presented Published in Prikladnaya Mekhanika, Vol. 44, No. 4, pp. 134–144, April 2008.     

Structure of the nonisothermal swirling gas-droplet flow behind an abrupt tube expansion

Flow structure and heat and mass transfer in a swirling two-phase stream is numerically modeled using the Reynolds stress transport model. The gas phase is described by the 3DRANS system of equations with account for the inverse influence of particles on the transport processes in the gas. The gas phase turbulence is calculated using the Reynolds stress transport model with account for the presence of disperse particles. The two-phase nonswirling flow behind an abrupt tube expansion contains a secondary corner vortex which is absent from the swirling flow. The disperse phase is redistributed over the tube cross-section. Large particles are concentrated in the wall region of the channel under the action of the centrifugal forces, while the smaller particles are in the central zone of the chamber.     

Some Specific Features of Integral Equations with the Cauchy Kernel on a Closed Contour in Hydrodynamic Problems

Singular integral equations of the first and second kind with the Cauchy kernel on a limiting narrow closed contour are theoretically considered. The initial equations are found to become different on the limiting contour. This singularity of integral equations with the Cauchy kernel does not allow boundary-value problems of the flow around thin airfoils to be solved correctly; therefore, a system consisting of integral equations of the first and second kind is proposed for solving such problems. The results of the present study are tested against an exact solution of the problem of the flow past a flat plate.     

Unsteady motion of a viscous incompressible fluid in a tube with a deformable wall

A flow of a viscous incompressible fluid in a deformable tube is considered. Solutions of unsteady three-dimensional Navier-Stokes equations are obtained for low-Reynolds-number flows in the tube (under the condition of small deformations of the wall): generalized peristaltic flow and flow with elliptical deformations of the vessel walls. At small unsteady deformations of the tube walls, the solutions satisfy the equations and boundary conditions with an error smaller than the tube wall deformation level by an order of magnitude. In the case of elliptical deformations of the vessel, the solution agrees well with experimental data.     

Semiempirical theory of turbulent transfer

A closed system of equations is constructed for flow of a nonisotropic character on the assumption that the mixing path length is not small compared with the characteristic dimension of the stream. It is assumed that the velocity pulsation field can be characterized by a multipoint distribution function, which satisfies the continuity equation. This enables equations to be obtained for the one-point and two-point distribution function. A series of assumptions is made concerning the nature of the forces acting on the turbulent formation (turbule or vortex) in the stream and concerning the correlation time of the random force with the scale and intensity of the turbulence. Assumptions are also made concerning the expression of the integral in the equation for the one-point distribution function and the expression for the correlation tensor in the isotropic case. After the moments are calculated, a system of Reynolds' equations is obtained in which approximations, usually acceptable from dimensional considerations, follow for a series of terms. Here, this is a consequence of the approximations for the forces in the equation for the distribution function. Closing the system of equations for the moments reduces to solving the equation for the distribution function. It turns out that the integral character of the transfer (diffusion of a nongradient type) is connected with taking third-order moments into account. A series of examples of flow is considered, and values of the empirical constants are determined. The system of equations obtained enables us to consider flow with strong anisotropy of turbulent transfer.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 11, No. 2, pp. 62–71, March–April, 1970.     

Lateral migrations of solid spheres in tube flow

Cross-flow displacements of neutrally buoyant solid spheres under Poiseuille-type flow conditions in narrow channels, at low volume fraction, are described by a phenomenological model within the framework of the Fokker-Planck equation. In this model, the effects of particle-particle and particle-wall interactions are taken into account. In addition to the classical wall effect (depleted layer near the wall), numerical solutions of the derived equations display two opposite lateral migrations towards the wall and the tube axis, and an alignment of the particles in adjacent layers parallel to the flow. These numerical results are then compared with experimental results obtained by Goldsmith and Marlow [18] for erythrocyte suspensions. A good qualitative agreement is observed.     

Boundary layer in the vicinity of the stagnation point of a body of axial symmetry in a turbulent stream

The flow in the boundary layer in the vicinity of the stagnation point of a flat plate is examined. The outer stream consists of turbulent flow of the jet type, directed normally to the plate. Assumptions concerning the connection between the pulsations in velocity and temperature in the boundary layer and the average parameters chosen on the basis of experimental data made it possible to obtain an isomorphic solution of the boundary layer equations. Equations are obtained for the friction and heat transfer at the wall in the region of gradient flow taking into account the effect of the turbulence of the impinging stream. It is shown that the friction at the wall is insensitive to the turbulence of the impinging stream, while the heat transfer is significantly increased with an increase in the pulsations of the outer flow. These properties are confirmed by the results of experimental studies [1–4].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 83–87, September–October, 1973.     

Dynamic Response of a Piezoelectric Material with a Conducting Rigid Inclusion

The problem of a conducting rigid inclusion embedded in an infinite piezoelectric matrix is considered under the action of combined electromechanical impact loads. By using integral transform techniques, the mixed initial-boundary value problem for the case of anti-plane shear load and in-plane electric field is transformed into two systems of dual integral equations, the solutions of which give the singularity coefficients of electroelastic field near the inclusion tips in closed-form in the Laplace transform domain. Numerical results for the stress singularity coefficient in the physical space are presented graphically by numerically solving the resulting Fredholm integral equation and carrying out the numerical inversion of Laplace transform for a PZT-5H material with a conducting rigid line inclusion.     

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube

The effect of yield stress on the flow characteristics of a Casson fluid in a homogeneous porous medium bounded by a circular tube is investigated by employing the Brinkman model to account for the Darcy resistance offered by the porous medium. The non-linear coupled implicit system of differential equations governing the flow is first transformed into suitable integral equations and are solved numerically. Analytical solution is obtained for a Newtonian fluid in the case of constant permeability, and the numerical solution is verified with that of the analytic solution. The effect of yield stress of the fluid and permeability of the porous medium on shear stress and velocity distributions, plug flow radius and flow rate are examined. The minimum pressure gradient required to start the flow is found to be independent of the permeability of the porous medium and is equal to the yield stress of the fluid.     

计算动力学中的伪弧长方法研究

动力学问题通常采用微分方程来描绘,但由于工程实际问题的复杂性,微分方程模型常伴随着解的不连续性、刚性或激波间断奇异性特点,传统方法很难求解,奇异性问题是计算动力学难点,同时也是国内外学者研究的热点.伪弧长数值算法是针对计算动力学中的奇异性问题所提出的,其基本思想为通过在解曲线上引入伪弧长参数,并增加一个约束方程,在伪弧长参数作用下,使得原始离散单元发生扭曲形变,从而达到消除或减弱奇异性的目的.本文首先介绍伪弧长方法求解定常对流-扩散方程的奇异性问题,并提出针对双曲守恒定律的局部伪弧长算法,其思想在于首先通过间断解的梯度变换来确定强间断所处位置,进而通过局部网格点重构以及数值修正来达到强间断处奇异性消除与降低的目的.针对高维问题,提出全局伪弧长方法,通过对整个计算区域内的网格点进行重构,使得所有网格点向奇异间断点处移动,从而降低间断点的影响域,达到降低奇异性的目的.重点讨论了三维全局伪弧长算法问题的计算难点,即三维空间网格扭曲大变形导致的数值算法不收敛,并提出在算法设计过程中采用分块重构与整体计算相结合的策略,实现了三维空间中的伪弧长数值算法,最后通过数值实验来验证伪弧长算法对于奇异性问题的有效性.     

3-D elastic stress singularity at polyhedral corner points

The three-dimensional stress singularity at the top of an arbitrary polyhedral corner is considered. Based on the boundary integral equations, the problem is reduced by the Mellin transform to a system of certain one-dimensional integral equations. The orders of stress singularity are spectral points of the integral operators while angular distribution and intensity factors are found as residues at those points. Numerical results are obtained by means of the Galerkin discretization scheme using expansions in terms of orthogonal polynomials with the proper weights. Some of the results illustrating the orders dependence on the elastic properties and corner geometry for a wedge-shaped punch and a crack, for an elastic trihedron and for a surface-breaking crack are given.     

Low Reynolds number hydrodynamic interaction of a solid particle with a planar wall

The slow viscous flow problem of an arbitrary solid particle in motion near a planar wall is recast into a boundary integral formulation. The present formulation employs the Green function appropriate to the planar wall problem and is developed in sufficient generality to allow calculations for arbitrary particles in any base flow which satisfies Stokes equations and no-slip on the wall. The resulting integral equations are easily discretized and solved for the particle surface tractions. Calculations are performed for axisymmetric motions of a variety of ellips?ids near the planar wall. Agreement with existing theory is excellent.     

Numerical simulation of flow of micropolar fluids in a channel with a porous wall

Two‐dimensional steady, laminar, and incompressible flow of a micropolar fluid in a channel with no‐slip at one wall and constant uniform injection through the other wall is considered for different values of the Reynolds number R. The main flow stream is superimposed by constant injection velocity at the porous wall. The micropolar model introduced by Eringen is used to describe the working fluid. An extension of Berman's similarity transformations is used to reduce governing equations to a set of nonlinear coupled ordinary differential equations (ODEs) in dimensionless form. An algorithm based on finite difference method is employed to solve these ODEs and Richardson's extrapolation is used to obtain higher order accuracy. It has been found that the magnitude of shear stress increases strictly at the impermeable wall whereas it decreases steadily at the permeable wall, by increasing the injection velocity. The maximum value of streamwise velocity and that of the microrotation both increase with increasing the magnitude of R. Copyright © 2010 John Wiley & Sons, Ltd.     

Nonstationary hydrodynamic characteristics of a double grid of profiles

The flow about a double grid of solid profiles of arbitrary shape which vibrate in a stream of an ideal incompressible fluid is considered. Behind the grid profiles, the nonstationary vortex traces simulated by the lines of contact velocity discontinuity are taken into account. The problem is reduced to the solution of a system of two integral equations relative to the fluid velocity on the initial profiles of the double grid under the assumption that the vibration amplitudes are small. Formulas for calculating the nonstationary forces and moments are derived. The dependences of these forces on the shape, mutual positions, and laws of vibration of the grid profiles are studied. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 4, pp. 150–155, July–August, 1999.     

In-plane stress analysis of an orthotropic plane containing multiple defects

The solution of Volterra type climb and glide edge dislocations is utilized to formulate integral equations for an orthotropic homogeneous infinite plane weakened by multiple smooth cracks and/or cavities. Cavities are considered as closed curved cracks without singularity. The integral equations are of Cauchy singular type which are converted to hypersingular integral equations. These equations are then solved numerically to determine stress intensity factors for cracks and hoop stress on the cavities. The results for isotropic and orthotropic planes are compared with available solutions in literature and excellent agreement is observed. The formulation allows stress analysis of orthotropic planes with several arbitrarily oriented cracks and cavities.     

计算离心泵叶轮流场的扫描流束有限元方法

本文提出一种求解离心式叶轮流场的数值方法，将流动求解区域离散为有限个由流线构成其边界的单元，采用伽辽金法建立的单元方程在一条流束上集合为方程组，流线上的节点坐标亦作为未知量包含在有限元方程中，通过扫描计算，逐步解得流线位置及流动参数。本文应用叶轮的通流理论流动模型，采用扫描流速有限元方法对离心泵叶轮流场进行了计算，并与有关文献作了比较。