A combined analytical–numerical method for treating corner singularities in viscous flow predictions

A combined analytical–numerical method based on a matching asymptotic algorithm is proposed for treating angular (sharp corner or wedge) singularities in the numerical solution of the Navier–Stokes equations. We adopt an asymptotic solution for the local flow around the angular points based on the Stokes flow approximation and a numerical solution for the global flow outside the singular regions using a finite‐volume method. The coefficients involved in the analytical solution are iteratively updated by matching both solutions in a small region where the Stokes flow approximation holds. Moreover, an error analysis is derived for this method, which serves as a guideline for the practical implementation. The present method is applied to treat the leading‐edge singularity of a semi‐infinite plate. The effect of various influencing factors related to the implementation are evaluated with the help of numerical experiments. The investigation showed that the accuracy of the numerical solution for the flow around the leading edge can be significantly improved with the present method. The results of the numerical experiments support the error analysis and show the desired properties of the new algorithm, i.e. accuracy, robustness and efficiency. Based on the numerical results for the leading‐edge singularity, the validity of various classical approximate models for the flow, such as the Stokes approximation, the inviscid flow model and the boundary layer theory of varying orders are examined. Although the methodology proposed was evaluated for the leading‐edge problem, it is generally applicable to all kinds of angular singularities and all kinds of finite‐discretization methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Wave induced thermal boundary layers in a compressible fluid: analysis and numerical simulations

The response of a semi-infinite compressible fluid to a step-wise change in temperature of its boundary is investigated analytically and numerically. Numerical results of the boundary layer structure are compared with Clarke’s analytical solution for a gas with thermal conductivity proportional to temperature. To avoid unwanted numerical dissipation in the numerical analysis, the space-time conservation element and solution element (CESE) method has been adopted to solve the unsteady 1-D Navier-Stokes equations. Good agreement between analytical and numerical results has been found for the development of the thermal boundary layer on a long time scale. Weak shock waves and expansion waves induced by the thermal boundary layer due to its compressibility, are observed in the numerical simulation. Finally, the numerical method has been applied to the reflection of a non-linear expansion wave and to a shock wave from an isothermal wall, thereby illustrating the effect of the boundary layer on the external flow field.     

A hybrid numerical approach to predict the vibrational responses of panels excited by a turbulent boundary layer

In this work, a hybrid numerical approach to predict the vibrational responses of planar structures excited by a turbulent boundary layer is presented. The approach combines an uncorrelated wall plane wave technique with the finite element method. The wall pressure field induced by a turbulent boundary layer is obtained as a set of uncorrelated wall pressure plane waves. The amplitude of these plane waves are determined from the cross spectrum density function of the wall pressure field given either by empirical models from literature or from experimental data. The response of the planar structure subject to a turbulent boundary layer excitation is then obtained from an ensemble average of the different realizations. The numerical technique is computationally efficient as it rapidly converges using a small number of realizations. To demonstrate the method, the vibrational responses of two panels with simply supported or clamped boundary conditions and excited by a turbulent flow are considered. In the case study comprising a plate with simply supported boundary conditions, an analytical solution is employed for verification of the method. For both cases studies, numerical results from the hybrid approach are compared with experimental data measured in two different anechoic wind tunnels.     

MHD boundary layer flow and heat transfer over a stretching sheet with induced magnetic field

In this paper, the problem of steady magnetohydrodynamic boundary layer flow and heat transfer of a viscous and electrically conducting fluid over a stretching sheet is studied. The effect of the induced magnetic field is taken into account. The transformed ordinary differential equations are solved numerically using the finite-difference scheme known as the Keller-box method. Numerical results are obtained for various values of the magnetic parameter, the reciprocal magnetic Prandtl number and the Prandtl number. The effects of these parameters on the flow and heat transfer characteristics are determined and discussed in detail. When the magnetic field is absent, the closed analytical results for the skin friction are compared with the exact numerical results. Also the numerical results for the heat flux from the stretching surface are compared with the results reported by other authors when the magnetic field is absent. It is found that very good agreement exists.     

Falkner-Skan方程的近似解析解

研究了粘性流体绕流楔型物体的Falkner-Skan边界层方程求解问题.利用Adomian拆分方法,通过引入Crocco变量变换将无穷区间的边界值问题转为初值问题并利用Padé逼近技巧确定初值,给出了一种有效的解析分解方法.进一步,本文设计了一种数值解法,将本文得到的近似解析解及数值结果与早期研究者Hartree等人的结果进行了比较,证明了本文提出的解法的有效性和可靠性.     

Thermal Analysis of Flow in a Porous Medium Over a Permeable Stretching Wall

This work presents a similarity solution for boundary layer flow through a porous medium over a stretching porous wall. Two considered wall boundary conditions are power-law distribution of either wall temperature or heat flux which are general enough to cover the isothermal and isoflux cases. In addition to momentum, both first and second laws of thermodynamics analyses of the problem are investigated. Independent numerical simulations are also performed for verification of the proposed analytical solution. The results, from the two independent approaches, are found to be in complete agreement. A comprehensive parametric study is presented and it is shown that heat transfer and entropy generation rates increase with Reynolds number, Prandtl number, and suction to the surface.     

Dual solutions in MHD stagnation-point flow of Prandtl fluid impinging on shrinking sheet

The present article investigates the dual nature of the solution of the magneto- hydrodynamic （MHD） stagnation-point flow of a Prandtl fluid model towards a shrinking surface. The self-similar nonlinear ordinary differential equations are solved numerically by the shooting： method. It is found that the dual solutions of the flow exist for cer- tain values of tile velocity ratio parameter. The special case of the first branch solutions （the classical Newtonian fluid model） is compared with the present numerical results of stretching flow. The results are found to be in good agreement. It is also shown that the boundary layer thickness for the second solution is thicker than that for the first solution.     

Turbulent flow in converging nozzles,part one: boundary layer solution

The boundary layer integral method is used to investigate the development of the turbulent swirling flow at the entrance region of a conical nozzle. The governing equations in the spherical coordinate system are simplified with the boundary layer assumptions and integrated through the boundary layer. The resulting sets of differential equations are then solved by the fourth-order Adams predictor-corrector method. The free vortex and uniform velocity profiles are applied for the tangential and axial velocities at the inlet region, respectively. Due to the lack of experimental data for swirling flows in converging nozzles, the developed model is validated against the numerical simulations. The results of numerical simulations demonstrate the capability of the analytical model in predicting boundary layer parameters such as the boundary layer growth, the shear rate, the boundary layer thickness, and the swirl intensity decay rate for different cone angles. The proposed method introduces a simple and robust procedure to investigate the boundary layer parameters inside the converging geometries.     

Accurate solutions for viscoelastic boundary layer flow and heat transfer over stretching sheet

In this article, we present accurate analytical solutions for boundary layer flow and heat transfer of an incompressible and electrically conducting viscoelastic fluid over a linearly stretching surface subject to a transverse uniform magnetic field using the homotopy analysis method （HAM） for two general types of non-isothermal boundary conditions. In addition, we demonstrate that the previously reported analytical solutions for the temperature field given in terms of Kummer＇s function do not converge at the boundary. We provide a graphical and numerical demonstration of the convergence of the HAM solutions and tabulate the effects of various parameters on the skin friction coefficient and wall heat transfer.     

An exact analytical solution of the unsteady magnetohydrodynamics nonlinear dynamics of laminar boundary layer due to an impulsively linear stretching sheet

In this paper, we investigate the theoretical analysis for the unsteady magnetohydrodynamic laminar boundary layer flow due to impulsively stretching sheet. The third-order highly nonlinear partial differential equation modeling the unsteady boundary layer flow brought on by an impulsively stretching flat sheet was solved by applying Adomian decomposition method and Pade approximants. The exact analytical solution so obtained is in terms of rapidly converging power series and each of the variants are easily computable. Variations in parameters such as mass transfer (suction/injection) and Chandrasekhar number on the velocity are observed by plotting the graphs. This particular problem is technically sound and has got applications in expulsion process and related process in fluid dynamics problems.     

Prediction of periodic boundary layers

A relatively simple, yet efficient and accurate finite difference method is developed for the solution of the unsteady boundary layer equations for both laminar and turbulent flows. The numerical procedure is subjected to rigorous validation tests in the laminar case, comparing its predictions with exact analytical solutions, asymptotic solutions, and/or experimental results. Calculations of periodic laminar boundary layers are performed from low to very high oscillation frequencies, for small and large amplitudes, for zero as well as adverse time-mean pressure gradients, and even in the presence of significant flow reversal. The numerical method is then applied to predict a relatively simple experimental periodic turbulent boundary layer, using two well-known quasi-steady closure models. The predictions are shown to be in good agreement with the measurements, thereby demonstrating the suitability of the present numerical scheme for handling periodic turbulent boundary layers. The method is thus a useful tool for the further development of turbulence models for more complex unsteady flows.     

MESHLESS ANALYSIS FOR THREE-DIMENSIONAL ELASTICITY WITH SINGULAR HYBRID BOUNDARY NODE METHOD

The singular hybrid boundary node method (SHBNM) is proposed for solving three-dimensional problems in linear elasticity. The SHBNM represents a coupling between the hybrid displacement variational formulations and moving least squares (MLS) approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the later. The rigid movement method was employed to solve the hyper-singular integrations. The 'boundary layer effect', which is the main drawback of the original Hybrid BNM, was overcome by an adaptive integration scheme. The source points of the fundamental solution were arranged directly on the boundary. Thus the uncertain scale factor taken in the regular hybrid boundary node method (RHBNM) can be avoided. Numerical examples for some 3D elastic problems were given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method were studied through the numerical examples.     

An exact analytical solution of the Falkner-Skan equation with mass transfer and wall stretching

In this paper, an exact analytical solution of the famous Falkner-Skan equation is obtained. The solution involves the boundary layer flow over a moving wall with mass transfer in presence of a free stream with a power-law velocity distribution. Multiple solution branches are observed. The effects of mass transfer and wall stretching are analyzed. Interesting velocity profiles including velocity overshoot and reversal flows are observed in the presence of both mass transfer and wall stretching. These solutions greatly enrich the analytical solution for the celebrated Falkner-Skan equation and the understanding of this important and interesting equation.     

Nonlinear beam formulation incorporating surface energy and size effect: application in nano-bridges

A nonlinear beam formulation is presented based on the Gurtin-Murdoch surface elasticity and the modified couple stress theory. The developed model theoretically takes into account coupled effects of the energy of surface layer and microstructures sizedependency. The mid-plane stretching of a beam is incorporated using von-Karman nonlinear strains. Hamilton’s principle is used to determine the nonlinear governing equation of motion and the corresponding boundary conditions. As a case study, pull-in instability of an electromechanical nano-bridge structure is studied using the proposed formulation. The nonlinear governing equation is solved by the analytical reduced order method (ROM) as well as the numerical solution. Effects of various parameters including surface layer, size dependency, dispersion forces, and structural damping on the pullin parameters of the nano-bridges are discussed. Comparison of the results with the literature reveals capability of the present model in demonstrating the impact of nanoscale phenomena on the pull-in threshold of the nano-bridges.     

Numerical studies for flow and heat transfer of the Powell-Eyring fluid thin film over an unsteady stretching sheet with internal heat generation using the chebyshev finite difference method

An analysis is carried out to study the unsteady two-dimensional Powell-Eyring flow and heat transfer to a laminar liquid film from a horizontal stretching surface in the presence of internal heat generation. The flow of a thin fluid film and subsequent heat transfer from the stretching surface is investigated with the aid of a similarity transformation. The transformation enables to reduce the unsteady boundary layer equations to a system of nonlinear ordinary differential equations. A numerical solution of the resulting nonlinear differential equations is found by using an efficient Chebyshev finite difference method. A comparison of numerical results is made with the earlier published results for limiting cases. The effects of the governing parameters on the flow and thermal fields are thoroughly examined and discussed.     

Extended homotopy perturbation method and the axisymmetric flow past a porous stretching sheet

The extended homotopy perturbation method, which is an extension of the celebrated homotopy perturbation method (HPM), is applied to obtain a solution to the problem of the steady, laminar, axisymmetric flow of a viscous, incompressible fluid past a porous stretching sheet. The solution so obtained is totally analytical and is expressible in terms of the cross‐flow velocity of the fluid past the stretching sheet. Its hallmark is that it does not depend upon computation of any auxiliary parameter for enlarging the convergence region of the solution. Rather, it calculates the solution automatically adjusting the scaling factor of the independent similarity variable normal to the sheet. The results obtained by the extended HPM are in excellent agreement with the exact numerical solution. Also, an asymptotic solution valid for large suction parameter is developed, which matches well with the exact solution even for moderate values of the suction parameter. Copyright © 2011 John Wiley & Sons, Ltd.     

Flow of a thixotropic fluid over an exponentially stretching sheet with heat transfer

This article addresses the boundary layer flow of a thixotropic fluid past an exponentially stretching sheet with heat transfer. The governing partial differential equations are reduced to an ordinary differential equation whose solution is found by the homotopy analysis method. The numerical values of the skin friction coefficient and Nusselt number are compared with available data.     

断裂力学的半解析相似边界元法

本文首先对弹性力学的相似边界元法进行了研究，推导了相应的计算公式。与传统的边界元法相比，相似边界元法由于只需在少数单元上进行数值积分，大大减少了计算量。在此基础上，对断裂力学问题，利用裂纹尖端位移场的解析表达式将裂纹尖端节点未知量转化为几个待定常数，提出了半解析相似边界元法，可大大减少最终形成的线性代数方程组的系数矩阵的阶数，进一步减小计算量。最后给出了算例，说明了本文方法的有效性。     

Time-dependent three-dimensional flow and mass transfer of elastico-viscous fluid over unsteady stretching sheet

This article studies the three-dimensional boundary layer flow of an elasticoviscous luid over a stretching surface. Velocity of the stretching sheet is assumed to be ime-dependent. Effect of mass transfer with higher order chemical reaction is further onsidered. Computations are made by the homptopy analysis method (HAM). Convergence f the obtained series solutions is explicitly analyzed. Variations of embedding arameters on the velocity and concentration are graphically discussed. Numerical computations f surface mass transfer are reported. Comparison of the present results with he numerical solutions is also given.     

Time-dependent three-dimensional flow and mass transfer of elastico-viscous fluid over unsteady stretching sheet

This article studies the three-dimensional boundary layer flow of an elasticoviscous fluid over a stretching surface. Velocity of the stretching sheet is assumed to be time-dependent. Effect of mass transfer with higher order chemical reaction is further considered. Computations are made by the homptopy analysis method (HAM). Convergence of the obtained series solutions is explicitly analyzed. Variations of embedding parameters on the velocity and concentration are graphically discussed. Numerical computations of surface mass transfer are reported. Comparison of the present results with the numerical solutions is also given.