Evaporation heat transfer of R-134a inside a microfin tube with different tube inclinations

An experimental investigation has been carried out to study the heat transfer characteristics during evaporation of R-134a inside a single helical microfin tube. The microfin tube has been provided with different tube inclination angles of the direction of fluid flow from horizontal, α. The experiments were performed for seven different tube inclinations, α, in a range of −90° to +90° and four mass velocities of 53, 80, 107 and 136 kg/m² s for each tube inclination angle during evaporation of R-134a. The results demonstrate that the tube inclination angle, α, affects the boiling heat-transfer coefficient in a significant manner. For all refrigerant mass velocities, the best performing tube is that having inclination angle of α = +90°. The effect of tube inclination angle, α, on heat-transfer coefficient, h, is more prominent at low vapor quality and mass velocity. An empirical correlation has also been developed to predict the heat-transfer coefficient during flow boiling inside a microfin tube with different tube inclinations.     

Numerical simulation of the flow field in a reactor for titanium dioxide production

A physicochemical and fluid dynamic model is formulated for the numerical simulation of the flow field in a reactor for titanium dioxide production, the turbulence motion is described by theKε equation, the governing equations are solved by the SIMPLER algorithm devised by Partankar and Spalding. The velocity, tmperature and concentration fields are obtained for three cases: A) with chemical reaction and thermal insulation on the walls; B) with chemical reaction and wall temperature is 450K; C) without chemical reaction and thermal insulation on the walls, and the physicochemical numerical simulation for the titanium dioxide production has been done. The results of the paper can be used as a theoretical guide for the engineers in the design of such reactors.     

A Numerical Study of Laminar Reciprocating Flow in a Pipe of Finite Length

A numerical investigation has been carried out for a laminar incompressible reciprocating flow in a circular pipe with a finite length. An examination of the governing equations and boundary conditions indicates that a sinusoidally reciprocating flow is governed by three similarity parameters: the kinetic Reynolds number Re, the dimensionless oscillation amplitude A_o, and the length to diameter ratio L/D. The numerical solution for the velocity profiles of a developing reciprocating flow shows that at any instant of time, there exist three flow regimes in the pipe, namely, an entrance regime, a fully developed regime and an exit regime. The numerical results for the fully developed region are shown to be in excellent agreement with the analytical solution. Based on the numerical results, a correlation equation of the space-cycle averaged friction coefficient for a laminar developing reciprocating pipe flow has been obtained in terms of the three similarity parameters.     

Investigation of the structure of a strong electric discharge in a high-velocity air flow

Many studies devoted to both the theoretical and the experimental analysis of the problem of controlling the aerodynamic characteristics of vehicles using local heat-generating zones organized by means of external fuel combustion or directly as a result of the Joule dissipation of an electric current have recently been published. In this study attention is concentrated on the latter case.Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, 2004, pp. 163–175. Original Russian Text Copyright © 2004 by Alfyorov.     

Brachistochrone on a surface with Coulomb friction

The problem of brachistochronic motion of a particle on a surface with the simultaneous action of gravity and Coulomb friction has been solved. Analytical solutions of the problem in special cases have been found. In the case of a cylindrical surface, it is shown that the solutions found includes, as a special case, all well-known results which refer to the brachistochronic motion of a particle. The results are illustrated with a series of concrete examples.     

The effect of particles on fluid turbulence in a turbulent boundary layer over a cylinder

The turbulent fluid and particle interaction in the turbulent boundary layer for cross flow over a cylinder has been experimentally studied. A phase-Doppler anemometer was used to measure the mean and fluctuating velocities of both phases. Two size ranges of particles (30μm–60μm and 80μm–150μm) at certain concentrations were used for considering the effects of particle sizes on the mean velocity profiles and on the turbulent intensity levels. The measurements clearly demonstrated that the larger particles damped fluid turbulence. For the smaller particles, this damping effect was less noticeable. The measurements further showed a delay in the separation point for two phase turbulent cross flow over a cylinder. The project supported by the National Natural Science Foundation of China     

幂硬化弹塑性-刚性界面Ⅲ型定常扩展裂纹尖端的渐近场

对幂硬化弹塑性材料－刚性材料界面上裂纹以定常方式扩展的Ⅲ型问题进行弹塑性渐近分析，给出裂纹尖端的应力，应变和位移场解。通过数值计算，考察了不同Mach数以及裂纹尖端混合参数对场解的构造以及应力，应变分布的影响，为给出合理的断裂准则提供理论依据。     

Two-phase frictional pressure drop multipliers for SUVA R-134a flowing in a rectangular duct

The adiabatic two-phase frictional multipliers for SUVA, R-134a flowing in a rectangular duct (with D_H = 4.8 mm) have been measured for three nominal system pressures (0.9 MPa, T_sat = 35.5 °C; 1.38 MPa, T_sat = 51.8 °C; and 2.41 MPa, T_sat = 75.9 °C) and three nominal mass fluxes (510, 1020 and 2040 kg/m²/s). The data is compared with several classical correlations to assess their predictive capabilities. The Lockhart–Martinelli model gives reasonable results at the lowest pressure and mass flux, near the operating range of most refrigeration systems, but gives increasingly poor comparisons as the pressure and mass flux are increased. The Chisholm B-coefficient model is found to best predict the data over the entire range of test conditions; however, there is significant disagreement at the highest pressure tested (with the model over predicting the data upwards of 100% for some cases). The data shows an increased tendency toward homogeneous flow as the pressure and flow rate are increased, and in fact the homogeneous model best predicts the bulk of the data at the highest pressure tested.     

Pressure field around a well in a stratified inhomogeneous bed

The steady axisymmetric flow of an incompressible fluid into a vertical well hydrodynamically perfectly drilled into a stratified inhomogeneous half-space consisting of three layers with different permeabilities is considered. The boundaries of the layers are assumed to be horizontal planes and the roof of the upper layer is assumed to be impermeable. The flow obeys a linear Darcy’s law. The pressure distribution on the well is assumed to be given, which is the main obstacle to finding an exact solution of the problem. Beginning with the classical studies of Muskat and Charnyi [1, 2], approximate solutions of such problems have been constructed as a superposition of flows generated by point sources with given intensities, distributed along the well axis in accordance with a fairly simple law. In the present study, this approach is used to obtain an integral equation for the source density distribution, which is then solved numerically. Comparison with the known exact solution for a thin elongated ellipsoid (“needle”) shows that this approach makes it possible to ensure an accuracy which at any rate is sufficient for applications.     

Slow gas motions and the negative drag of a strongly heated spherical particle

In the slow flows of a strongly and nonuniformly heated gas, in the continuum regime (Kn → 0) thermal stresses may be present. The theory of slow nonisothermal continuum gas flows with account for thermal stresses was developed in 1969–1974. The action of the thermal stresses on the gas results in certain paradoxical effects, including the reversal of the direction of the force exerted on a spherical particle in Stokes flow. The propulsion force effect is manifested at large but finite temperature differences between the particle and the gas. This study is devoted to the thermal-stress effect on the drag of a strongly heated spherical particle traveling slowly in a gas for small Knudsen numbers (M ～ Kn → 0), small but finite Reynolds numbers (Re ≤ 1), a linear temperature dependence of the transport coefficients µ ∝ T, and large but finite temperature differences ((T _w ? T _∞ )/T _M8 ～ 1). Two different systems of equations are solved numerically: the simplified Navier-Stokes equations and the modified Navier-Stokes equations with account for the thermal stresses.     

Repair of a Plate with a Circular Hole by Applying a Patch

The stressed state of a thin elastic infinite plate with a circular hole covered by a circular patch of a greater radius is considered. The center of the hole coincides with the center of the patch. The patch is attached to the plate along its entire boundary. Stresses are prescribed at infinity on the plate and at the hole boundary. Complex Muskhelishvili potentials are found by the method of power series, and the behavior of stresses on the patch–plate interface and at the hole boundary is studied.     

Flow of a Viscoplastic Fluid on a Rotating Disk

ＦＬＯＷＯＦＡＶＩＳＣＯＰＬＡＳＴＩＣＦＬＵＩＤＯＮＡＲＯＴＡＴＩＮＧＤＩＳＫＦａｎＣｈｕｎ（范椿）（ＩｎｓｔｉｕｉｅｏｆＭｅｃｈａｎｉｃｓ，ＡｃａｄｅｍｉａＳｉｎｉｃａ，Ｂｅｉｊｉｎｇ）（ＲｅｃｅｉｖｅｄＮｏｖ．２０，１９９２；Ｃｏｍｍｕｎｉｃａｔ...     

The Explicit Determination of the Logarithm of a Tensor and Its Derivatives

The logarithm of a tensor is often used in nonlinear constitutive relations of elastic materials. Here we show how the logarithm of an arbitrary tensor can be explicitly evaluated for any underlying space dimension n. We also present a method for the explicit evaluation of the derivatives of the logarithm of a tensor.      

Dynamics of a Flexible Cantilever Beam Carrying a Moving Mass

The motion of a flexible cantilever beam carrying a moving spring-mass system is investigated. The beam is assumed to be an Euler–Bernouli beam. The motion of the system is described by a set of two nonlinear coupled partial differential equations where the coupling terms have to be evaluated at the position of the mass. The nonlinearities arise due to the coupling between the mass and the beam. Due to the nonlinearities the system exhibits internal resonance which is investigated in this work. The equations of motion are solved numerically using the Rayleigh–Ritz method and an automatic ODE solver. An approximate solution using the perturbation method of multiple scales is also obtained.     

Inertial motion of a body in an ideal fluid from the state at rest

The motion of a body in an ideal incompressible fluid flow without vortices in the absence of external forces is considered. It is demonstrated that the body can move inertially from the state at rest if its shape satisfies certain conditions. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 214–219, July–August, 2008.     

A least square extrapolation method for the a posteriori error estimate of the incompressible Navier Stokes problem

A posteriori error estimators are fundamental tools for providing confidence in the numerical computation of PDEs. To date, the main theories of a posteriori estimators have been developed largely in the finite element framework, for either linear elliptic operators or non‐linear PDEs in the absence of disparate length scales. On the other hand, there is a strong interest in using grid refinement combined with Richardson extrapolation to produce CFD solutions with improved accuracy and, therefore, a posteriori error estimates. But in practice, the effective order of a numerical method often depends on space location and is not uniform, rendering the Richardson extrapolation method unreliable. We have recently introduced (Garbey, 13th International Conference on Domain Decomposition, Barcelona, 2002; 379–386; Garbey and Shyy, J. Comput. Phys. 2003; 186 :1–23) a new method which estimates the order of convergence of a computation as the solution of a least square minimization problem on the residual. This method, called least square extrapolation, introduces a framework facilitating multi‐level extrapolation, improves accuracy and provides a posteriori error estimate. This method can accommodate different grid arrangements. The goal of this paper is to investigate the power and limits of this method via incompressible Navier Stokes flow computations. Copyright © 2005 John Wiley & Sons, Ltd.     

Stability and Transition on a Swept Cylinder in a Supersonic Flow

Results of experimental investigations of the evolution of natural disturbances and laminar–turbulent transition in a supersonic boundary layer on the attachment line of a circular cylinder with a sweep angle of 68° and a freestream Mach number M = 2 are presented. The experimental studies are supplemented by calculations of the mean flow and stability characteristics. Flow regimes in the boundary layer on the attachment line are determined by a hotwire technique as functions of the Reynolds number and height of twodimensional roughness elements. The results are compared with NASA (Ames) experiments.     

Negative wake generation of FENE-CR fluids in uniform and Poiseuille flows past a cylinder

The effect of flow conditions on the negative wake generation (longitudinal velocity overshoot behind a cylinder in the viscoelastic fluid flow along the centerline) has been investigated. FENE-CR model that predicts constant shear viscosity and controlled extensional viscosity was considered as a constitutive equation. The discrete elastic viscous split stress-G/streamline upwind Petrov–Galerkin (DEVSS-G/SUPG) formulation was employed and the high-resolution solutions were obtained with an efficient iterative solver based on the incomplete LU(0)-type preconditioner and BiCGSTAB. We found that the negative wake generation was more obvious in uniform flow conditions than in Poiseuille flow, which suggests that the experimentally unrevealed negative wake generation of Boger fluids could be partially attributed to the geometrical effect of Poiseuille flow. The negative wake generation was more discernable at low extensibility and high value of viscosity ratio, which agrees well with the previous studies. In addition, we could observe an undershoot phenomenon in Poisseuille flow condition, which has never been reported.     

Flow of viscoelastic fluids through a porous channel—I

The flow of viscoelastic fluids through a porous channel with one impermeable wall is computed. The flow is characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. Three solutions are developed: (i) an exact numerical solution, (ii) a perturbation solution for small R, the cross-flow Reynold's number and (iii) an asymptotic solution for large R. The results from exact numerical integration reveal that the solutions for a non-Newtonian fluid are possible only up to a critical value of the viscoelastic fluid parameter, which decreases with an increase in R. It is further demonstrated that the perturbation solution gives acceptable results only if the viscoelastic fluid parameter is also small. Two more related problems are considered: fluid dynamics of a long porous slider, and injection of fluid through one side of a long vertical porous channel. For both the problems, exact numerical and other solutions are derived and appropriate conclusions drawn.     

Bending of a circular cylinder containing a crack

The study of bending of cracked circular cylinders is of more significance. The bending of cylinders containing radical crack or cracks was discussed by refs. [1]–[4] and that of concentrically craked circular cylinders was studied by [5]. Continuing [6] and using complex variable methods in elasticity, this paper deals with the bending problems of a circular cylinder, containing an internal linear crack at any position under an acting force perpendicular to the crack. The general forms of displacements, stresses, and stressintensity factors, expressed in terms of series, are obtained and to this bending problems with small Ah are presented good approximate formulas for the stress-intensity factors whose variations with the center of the crack are analysed. Finally, the twist angle per unit length and the center of bending for the radically cracked circular cylinder, one of whose crack-tips is located at the origin, have been computed and the results are almost the same as that calculated in [1].