Inverse laminar boundary-layer problems with assigned shear. The mechul function revisited

An inverse method is presented which accurately determines the pressure distribution for assigned wall shear in a two-dimensional, laminar, incompressible boundary layer. The method reformulates the mechul function scheme of Cebeci and Keller to produce a stable solution in the marching direction and to increase accuracy in the normal direction. In the reformulation a modified pressure gradient parameter variation in the normal direction is used in conjunction with three-point backward differences for streamwise derivatives and fourth-order accurate splines for normal derivatives. The resulting spline-finite difference equations are solved by Newton-Raphson iteration together with partial pivoting. Numerical solutions are presented for self-similar and non self-similar flows and compared with published results.     

Transient simulation of 2–3D stratified and intermittent two-phase flows. Part II: Applications

In this paper a new type of transient multidimensional two-fluid model has been applied to simulate intermittent or slug flow problems. Three different approaches to modelling interfacial friction, including an interfacial tracking scheme, have been investigated. The numerial method is based on an implicit finite difference scheme, solved directly in two steps applying a separate equation for the pressure. 2D predictions of Taylor bubble propagation in horizontal and inclined channels have been compared with experimental data and analytical solutions. The 2D model has also been applied to investigate a number of special phenomena in slug flow, including slug initiation, bubble turning in downflow and the bubble centring process at large liquid flow rates.     

Effect of Inertial Terms on Fluid–Porous Medium Flow Coupling

The study considers an effect of the nonlinear inertial terms in the Brinkman filtration equation on the characteristics of coupled flows in a pure fluid and porous medium in the frameworks of two independent problems. The first problem is the forced boundary-layer flow overlying the Darcy–Brinkman porous medium. The Prandtl theory is used, and the self-similar equations are built to describe it. It is shown that the inertial terms have a valuable effect on the boundary-layer structure because of the large velocity gradient in the transition zone. The boundary-layer thickness in a porous medium rapidly grows at large Reynolds numbers. The velocity magnitude and gradient at the interface also change. The second independent problem is an analysis of the inertial terms effect on the flow stability. The neutral curves of the full and linearized flow models are built using the shooting method. They have different short-wave asymptotic, but there are no significant changes in the critical Reynolds numbers and corresponding wave numbers.     

Calculations of stratified wavy two-phase flow in pipes

Calculations of fully developed, stratified wavy gas–liquid pipe flow is presented. The wavy interface is represented by an equivalent interfacial roughness obtained from experimental data, which is made non-dimensional following the Charnock formulation. The two-dimensional, steady-state axial momentum equation is solved together with a two-layer turbulence model, which is modified to account for the roughness introduced at the interface. The governing equations are discretized using a finite difference method on a composite, overlapping grid with local grid refinement near the interface and the wall. The immersed interface method is used to make the numerical scheme well-defined across the interface, and a level set function is used to represent the interface. Numerical calculations are found to compare satisfactorily with experimental data.     

MHD non-Darcian flow through a non-isothermal vertical surface embedded in a porous medium with radiation

The effect of MHD on steady two-dimensional laminar mixed flow about a vertical porous surface is numerically analyzed. Also the effects of radiation and heat generation and absorption are considered. A power law variation of temperature along the vertical wall is assumed. The nonlinear boundary-layer equations were transformed and the resulting differential equations were solved by an implicit finite difference scheme (Keller box method). Numerical results for the velocity distribution and the temperature distribution are presented for various values of Prandtl number Pr, magnetic parameter, porous medium parameter and internal heat generation or absorption coefficient. Further validation with previous works is carried out.     

Pseudospectral analysis of the interfacial stability of free jets with different velocity profiles

A double fluid model for a liquid jet surrounded by a coaxial gas stream was constructed. The interfacial stability of the model was studied by Chebyshev pseudospectral method for different basic velocity profiles. The physical variables were mapped into computational space using a nonlinear coordinates transformation. The general eigenvalues of the dispersion relation obtained are solved by QZ method, and the basic characteristics and their dependence on the flow parameters are analyzed.     

Forced convective film condensation inside vertical tubes

A theoretical study of forced convective film condensation inside vertical tubes is presented. We propose a unified procedure for predicting the pressure gradient and condensation heat transfer coefficient of a vapor flowing turbulently in the core and associated with laminar or turbulent film on the tube wall. The analysis for the vapor flows is performed under the condition that the velocity profiles are locally self-similar. The laminar and turbulent film models equate the gravity, pressure and viscous forces, and consider the effect of interfacial shear. The transition from laminar to turbulent film depends not only on the liquid Reynolds number but also on the interfacial shear stress. In this work we also proposed a new eddy viscosity model which is divided into three regions: the inner region in liquid condensate near the wall; the interface region including both liquid and vapor; and the outer region for the vapor core. Comparisons of the theory with some published experimental data showed good agreement.     

Self-similar solutions of three-dimensional laminar magnetohydrodynamic boundary-layer equations

Self-similar solutions of three-dimensional boundary-layer equations of an incompressible fluid in ordinary hydrodynamics were considered in [1–3] et al. The present work looks for self-similar solutions of three-dimensional magnetohydrodynamic boundary-layer equations.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 9, No. 4, pp. 10–17, July–August, 1968.     

Radiative heat transfer in absorbing, emitting, and anisotropically scattering boundary-layer flows with reflecting boundary

The heat transfer in absorbing, emitting, and anisotropically scattering boundary-layer flows with reflecting boundary over a flat plate, over a 90-deg wedge, and in stagnation flow is solved by application of the Galerkin method with the particular solution boundary condition I _p (τ₀,ξ,?μ) of the equation of radiative transfer for an inhomogeneous term and the Box method. The exact integral expressions for the radiation part of this problem are developed. The coupling between convective and radiative heat transfer in boundary-layer flows is described by a set of nonlinear simultaneous equations including differential equations and integrodifferential equations. The Galerkin method and the particular solution boundary condition I _p (τ₀,ξ,?μ) are used to analyze the radiation part of the problem. The nonsimilar boundary-layer equations are solved by the Box method. The present numerical procedure solutions are compared in tables with the other exact treating results, the P-3, and P-1 approximation methods for the case of isotropically scattering boundary-layer flows. The effects of linearly anistropically scattering and reflecting surface are taken into account. It is found that the present method is a reliable and efficient numerical procedure and scattering leads to a reduction in the total heat flux. The influence of the forward-backward scattering parameter on the total heat flux decreases with the increase of the surface reflectivity.     

液体边界层对平板结构振动阻尼效应分析

本文分析了液体边界层对平板振动的阻尼效应，基于液体的Stokes边界层理论，建立了考虑液体边界层阻尼效应的平板振动单自由度系统运动方程，利用Laplace变换方法求解了运动方程，得到了平板结构的位移响应解答，制作了一个单自由度质量-弹簧系统来测量水体粘滞阻尼效应，并将理论与试验结果进行了对比分析，表明Stokes边界层理论能较好模拟水体的粘滞阻尼特性，边界层阻尼明显减小了结构的动力响应，本文最后讨论了所提的分析方法的适用范围．     

An approach to solving mechanics problems for materials with multiscale self-similar microstructure

This article proposes an efficient method for solving mechanics boundary value problems formulated for domains with multiscale self-similar microstructure. In particular, composite materials for which one of the phases has a fractal-like structure with scale cut-offs are considered. The boundary value problems are solved using a finite element procedure with enriched shape functions that incorporate information about the geometric complexity. The use of these shape functions makes possible the definition of a unique, parametrically defined model from which the solution for configurations with an arbitrary number of scales can be derived. The proposed method is primarily useful for structures with a large number of self-similar scales for which using the usual finite element method would be too expensive. In order to exemplify the method, a 2D composite with fractal microstructure is considered and several boundary value problems are solved.     

Simultaneous effect of rotation and natural convection on the flow about a liquid sphere

The problem of mixed convection around a liquid sphere that experiences a rotation about its axis parallel to the free stream is studied numerically using a finite- difference technique. The coupled boundary-layer energy and momentum equations are numerically solved over a wide range of Grashof number that represents the cases of aiding and opposing free convection and for wide range of the spin parameter Ta/Re². The surface of the sphere also rotates as a result of the shear stress exerted from the external flow of air. The effect of both parameters on the velocity components as well as the temperature within the thermal boundary-layer is presented. Results show that increasing the aiding free convection and the spin parameter cause increases in the shear stress and the local heat transfer coefficient.     

Choice of a self-similar solution in boundary-layer theory

If the speed of the outer flow at the edge of the boundary layer does not depend on the time and is specified in the form of a power-law function of the longitudinal coordinate, then a self-similar solution of the boundary-layer equations can be found by integrating a third-order ordinary differential equation (see [1–3]). When the exponent of the power in the outerflow velocity distribution is negative, a self-similar solution satisfying the equations and the usually posed boundary conditions is not uniquely determinable [4], A similar result was obtained in [5] for flows of a conducting fluid in a magnetic field. In the present paper we study the behavior of non-self-similar perturbations of a self-similar solution, enabling us to provide a basis for the choice of a self-similar solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 42–46, July–August, 1974.     

Numerical solution of solidification in a square prism using an algebraic grid generation technique

The solidification of an infinitely long square prism was analyzed numerically. A front fixing technique along with an algebraic grid generation scheme was used, where the finite difference form of the energy equation is solved for the temperature distribution in the solid phase and the solid–liquid interface energy balance is integrated for the new position of the moving solidification front. Results are given for the moving solidification boundary with a circular phase change interface. An algebraic grid generation scheme was developed for two-dimensional domains, which generates grid points separated by equal distances in the physical domain. The current scheme also allows the implementation of a finer grid structure at desired locations in the domain. The method is based on fitting a constant arc length mesh in the two computational directions in the physical domain. The resulting simultaneous, nonlinear algebraic equations for the grid locations are solved using the Newton-Raphson method for a system of equations. The approach is used in a two-dimensional solidification problem, in which the liquid phase is initially at the melting temperature, solved by using a front-fixing approach. The difference of the current study lies in the fact that front fixing is applied to problems, where the solid–liquid interface is curved such that the position of the interface, when expressed in terms of one of the coordinates is a double valued function. This requires a coordinate transformation in both coordinate directions to transform the complex physical solidification domain to a Cartesian, square computational domain. Due to the motion of the solid–liquid interface in time, the computational grid structure is regenerated at every time step.     

Transonic expansion of a penny-shaped crack

The problem of a penny-shaped crack propagating in an unbounded isotropic linear elastic body is solved. The crack expands from a zero radius in a self-similar manner and is assumed to have speed larger than the S-wave speed and less than the P-wave speed. A tensile load is directed normal to the crack at infinity. The method of self-similar potentials and rotational superposition are applied. Attention is given to the stress singularity at the crack border.     

Simultaneous heat and mass transfer inside a vertical tube in evaporating a heated falling alcohols liquid film into a stream of dry air

A numerical study of the evaporation in mixed convection of a pure alcohol liquid film: ethanol and methanol was investigated. It is a turbulent liquid film falling on the internal face of a vertical tube. A laminar flow of dry air enters the vertical tube at constant temperature in the downward direction. The wall of the tube is subjected to a constant and uniform heat flux. The model solves the coupled parabolic governing equations in both phases including turbulent liquid film together with the boundary and interfacial conditions. The systems of equations obtained by using an implicit finite difference method are solved by TDMA method. A Van Driest model is adopted to simulate the turbulent liquid film flow. The influence of the inlet liquid flow, Reynolds number in the gas flow and the wall heat flux on the intensity of heat and mass transfers are examined. A comparison between the results obtained for studied alcohols and water in the same conditions is made.     

Boundary-layer flow of a micropolar fluid due to a stretching wall

Summary  A Theoretical analysis is carried out to study the boundary-layer flow over a continuously moving surface through an otherwise quiescent micropolar fluid. The transformed boundary-layer equations are solved numerically for a power-law surface velocity using the Keller-box method. The effects of the micropolar K and exponent m parameters on the velocity and microrotation field as well as on the skin-friction group are discussed in a detailed manner. It is shown that there is a near-similarity solution of this problem. The accuracy of the present solution is also discussed. Accepted for publication 1 April 1996     

On stagnation point flow of Sisko fluid over a stretching sheet

The steady two-dimensional stagnation-point flow, represented by Sisko fluid constitutive model, over a stretching sheet is investigated theoretically. Using suitable similarity transformations, the governing boundary-layer equations are transformed into the self-similar non-linear ordinary differential equation. The transformed equation is then solved using a very efficient analytic technique namely the homotopy analysis method (HAM) and the HAM solutions are validated by the exact analytic solutions obtain in certain special cases. The influence of the power-law index (n), the material parameter (A) and the velocity ratio parameter (d/c) on the flow characteristics is studied and presented through several graphs. In addition, the local skin friction coefficient for several values of these parameters is tabulated and examined. The similarity solutions for both the Newtonian and the power-law fluids are presented as special cases of the analysis. The results obtained reveal that, in comparison with the Newtonian and the power-law fluids, the velocity profiles of the Sisko fluid are much faster (slower), for d/c<1 (d/c>1), respectively.     

Elasto-plastic buckling analysis of laminated plates including interfacial damage

Elasto-plastic buckling of orthotropic laminated plates, which include interfacial damage, is analyzed in detail. Firstly, a novel mixed hardening yield criterion, as an improvement of Hill’s counterpart, is proposed for the orthotropic materials on the basis of the plastic theory. And differing from Hill’s theory, the present yield criterion is related to the spherical tensor of stress. Then, the incremental elasto-plastic constitutive relations of the mixed hardening orthotropic materials are presented. Secondly, the incremental static equilibrium equations for laminated plates including interfacial damage are established based on Von-Karman type theory and the principle of minimum potential energy. Finally, the elasto-plastic buckling of laminated plates are solved by adopting the Galerkin method and iteration scheme. The numerical results show that buckling of the plate occurs easier due to the existence of interfacial damage, and the critical load trends to constant when the interfacial damage approaches a certain degree. Also, the effect of anisotropy on buckling is obvious and the analysis of elasto-plastic buckling is necessary.