共查询到20条相似文献,搜索用时 15 毫秒
1.
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. 相似文献
2.
Randi Moe 《国际流体数值方法杂志》1993,16(11):967-988
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. 相似文献
3.
Kirill Tsiberkin 《Transport in Porous Media》2018,121(1):109-120
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. 相似文献
4.
《International Journal of Multiphase Flow》2005,31(5):571-592
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. 相似文献
5.
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. 相似文献
6.
Pseudospectral analysis of the interfacial stability of free jets with different velocity profiles 总被引:1,自引:0,他引:1
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. 相似文献
7.
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. 相似文献
8.
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. 相似文献
9.
T.-M. Chen 《Heat and Mass Transfer》1997,32(6):411-417
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 (τ0,ξ,?μ) 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 (τ0,ξ,?μ) 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. 相似文献
10.
11.
12.
《International Journal of Solids and Structures》2007,44(24):7877-7890
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. 相似文献
13.
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/Re2. 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. 相似文献
14.
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. 相似文献
15.
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. 相似文献
16.
Y.L. Chung 《Theoretical and Applied Fracture Mechanics》1998,29(3):151-159
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. 相似文献
17.
Simultaneous heat and mass transfer inside a vertical tube in evaporating a heated falling alcohols liquid film into a stream of dry air 总被引:1,自引:0,他引:1
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. 相似文献
18.
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 相似文献
19.
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. 相似文献
20.
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. 相似文献