Free boundary effects on stability of two phase planar fluid motion in annulus. Part I: Migration of the stable mode

We study the linearized stability of a planar dynamical model describing two-phase perfect fluid circulating around a circle with a sufficiently large radius within a central gravitational field. The model is associated with the spatial and temporal structure of the zonally averaged global-scale atmospheric longitudinal circulation around the Earth. Two cases are studied separately; in the first one, the simulations were carried out using the rigid lid approximation at the upper boundary of the outer atmospheric layer. In the second one, the free boundary nonlinear conditions (kinematic and dynamic) were assumed on the outer atmospheric layer. For the both cases, a certain family of steady, explicit solutions which have circular streamlines was considered. The governing equations were linearized at these solutions to find the typical wave numbers of the interfacial wave perturbation to the basic state at which the destabilizing effect of shear, which overcomes the stabilizing effect of stratification, occurs. It is shown that for the both cases, the model always have the same two potentially unstable wave modes while there always exist two wave modes which are stable for any wavelengths. The behavior of the stable and unstable modes were compared for the both cases to investigate the effects of the free boundary on the mixing process at the interface.     

Asymptotic determination of the formation process of nonlinear distortion of one-dimensional pulses in a layered medium : PMM vol.40, n 6, 1976, pp. 1093–1103

Nonlinear effects in the propagation, reflection, and refraction of one-dimensional pulses in a medium consisting of two layers lying on a half-space are considered and analyzed. Properties of layers and of the half-space are different, and stresses are defined by an expansion in powers of strains. The initial pulse of finite duration is specified in the form of boundary condition at the surface of the external layer either for the deformation or for the dislocation rate, and the problem of wave pattern when the initial pulse amplitude tends to zero,i.e. in the case of small nonlinear effects, is solved.Problem is solved by the method of successive integration of nonhomogeneous linear wave equations, in which the solution of the linear problem is taken as the first approximation and the subsequent approximations are derived by approximating the nonlinear terms with the use of the preceding approximation.     

Stability of reflection and refraction of shocks on interface

In this paper we study the stability of the nonlinear wave structure caused by the attack of an incident shock on an interface of two different kinds of media. The attack will produce a reflected wave and a refracted wave, and also let the interface deflected. In this paper we will mainly study the case, when the reflected wave is a shock, and the flow between the reflected wave and the refracted shock is relatively subsonic. Our result indicates that the wave structure and the flow field for the reflection-refraction problem in this case is conditionally stable.To describe the motion of the fluid we use the inviscid Euler system as the mathematical model. The reflection-refraction problem can be reduced to a free boundary value problem, where the unknown reflected shock and refracted shock are free boundaries, and the deflected interface is also to be determined. In the proof of the existence and the stability of the corresponding wave structure we apply the Lagrange transformation to fix the interface and the decoupling technique to decouple the elliptic-hyperbolic composite system in its principal part. Meanwhile, some efficient weighted Sobolev estimates are established to derive the existence for corresponding nonlinear problems.     

On the use of homotopy analysis method for solving unsteady MHD flow of Maxwellian fluids above impulsively stretching sheets

In the present work, unsteady MHD flow of a Maxwellian fluid above an impulsively stretched sheet is studied under the assumption that boundary layer approximation is applicable. The objective is to find an analytical solution which can be used to check the performance of computational codes in cases where such an analytical solution does not exist. A convenient similarity transformation has been found to reduce the equations into a single highly nonlinear PDE. Homotopy analysis method (HAM) will be used to find an explicit analytical solution for the PDE so obtained. The effects of magnetic parameter, elasticity number, and the time elapsed are studied on the flow characteristics.     

Analytical approximate solution of the cooling problem by Adomian decomposition method

The Adomian decomposition method (ADM) can provide analytical approximation or approximated solution to a rather wide class of nonlinear (and stochastic) equations without linearization, perturbation, closure approximation, or discretization methods. In the present work, ADM is employed to solve the momentum and energy equations for laminar boundary layer flow over flat plate at zero incidences with neglecting the frictional heating. A trial and error strategy has been used to obtain the constant coefficient in the approximated solution. ADM provides an analytical solution in the form of an infinite power series. The effect of Adomian polynomial terms is considered and shows that the accuracy of results is increased with the increasing of Adomian polynomial terms. The velocity and thermal profiles on the boundary layer are calculated. Also the effect of the Prandtl number on the thermal boundary layer is obtained. Results show ADM can solve the nonlinear differential equations with negligible error compared to the exact solution.     

Rossby Solitary Waves in the Presence of a Critical Layer

This study considers the evolution of weakly nonlinear long Rossby waves in a horizontally sheared zonal current. We consider a stable flow so that the nonlinear time scale is long. These assumptions enable the flow to organize itself into a large‐scale coherent structure in the régime where a competition sets in between weak nonlinearity and weak dispersion. This balance is often described by a Korteweg‐de‐Vries equation. The traditional assumption of a weak amplitude breaks down when the wave speed equals the mean flow velocity at a certain latitude, due to the appearance of a singularity in the leading‐order equation, which strongly modifies the flow in a critical layer. Here, nonlinear effects are invoked to resolve this singularity, because the relevant geophysical flows have high Reynolds numbers. Viscosity is introduced in order to render the nonlinear‐critical‐layer solution unique, but the inviscid limit is eventually taken. By the method of matched asymptotic expansions, this inner flow is matched at the edges of the critical layer with the outer flow. We will show that the critical‐layer–induced flow leads to a strong rearrangement of the related streamlines and consequently of the potential‐vorticity contours, particularly in the neighborhood of the separatrices between the open and closed streamlines. The symmetry of the critical layer vis‐à‐vis the critical level is also broken. This theory is relevant for the phenomenon of Rossby wave breaking and eventual saturation into a nonlinear wave. Spatially localized solutions are described by a Korteweg‐de‐Vries equation, modified by new nonlinear terms; depending on the critical‐layer shape, this leads to depression or elevation waves. The additional terms are made necessary at a certain order of the asymptotic expansion while matching the inner flow on the dividing streamlines. The new evolution equation supports a family of solitary waves. In this paper we describe in detail the case of a depression wave, and postpone for further discussion the more complex case of an elevation wave.     

弱粘性流体中垂直激励表面波的阻尼效应

在竖直振动的圆柱形容器中,将Navier-Stokes方程线性化,利用两时间尺度奇异摄动展开法研究了弱粘性流体的单一自由面驻波运动．整个流场被分为外部势流区和内部边界层区两部分,对两部分区域分别求解,得到包含阻尼项和外驱动影响的线性振幅方程．利用稳定性分析,得到形成稳定表面波的条件,给出了临界曲线．此外,还获得了阻尼系数的解析表达式．最后,将线性阻尼加到理想流体条件下所得到的色散关系中对其进行修正,理论结果证明修正后的驱动频率更加接近实验的结果．通过计算发现,当驱动的频率较低时,流体的粘性对表面波模式选择有重要影响,而表面张力的影响不明显;但当驱动频率较高时,流体的表面张力起主要作用,而流体的粘性影响甚小．     

A numerical method to compute the dissolution of second phases in ternary alloys

Dissolution of stoichiometric multi-component particles in ternary alloys is an important process occurring during the heat treatment of as-cast aluminium alloys prior to hot extrusion. A mathematical model is proposed to describe such a process. In this model an equation is given to determine the position of the particle interface in time, using two diffusion equations which are coupled by nonlinear boundary conditions at the interface. Some results concerning existence, uniqueness, and monotonicity are given. Furthermore, for an unbounded domain an analytical approximation is derived. The main part of this work is the development of a numerical solution method. Finite differences are used on a grid which changes in time. The discretization of the boundary conditions is important to obtain an accurate solution. The resulting nonlinear algebraic system is solved by the Newton-Raphson method. Numerical experiments illustrate the accuracy of the numerical method. The numerical solution is compared with the analytical approximation.     

Diffusion approximation for nonlinear evolutionary equations with large interaction and fast boundary fluctuation

Approximations are derived for both nonlinear heat equations and singularly perturbed nonlinear wave equations with highly oscillating random force on boundary and strong interaction. By a diffusion approximation method, if the interaction is large and the singular perturbation is small enough, the approximation of the nonlinear wave equation is an one dimensional stochastic ordinary differential equation with white noise from the boundary which is exactly the same as that of the nonlinear heat equation.     

Numerical treatment of singularly perturbed two point boundary value problems exhibiting boundary layers

A boundary value method for solving a class of nonlinear singularly perturbed two point boundary value problems with a boundary layer at one end is proposed. Using singular perturbation analysis the method consists of solving two problems; namely, a reduced problem and a boundary layer correction problem. We use Pade’ approximation to obtain the solution of the latter problem and to satisfy the condition at infinity. Numerical examples will be given to illustrate the method.     

Continuum models in problems of the hypersonic flow of a rarefied gas around blunt bodiest

All possible continuum (hydrodynamic) models in the case of two-dimensional problems of supersonic and hypersonic flows around blunt bodies in the two-layer model (a viscous shock layer and shock-wave structure) over the whole range of Reynolds numbers, Re, from low values (free molecular and transitional flow conditions) up to high values (flow conditions with a thin leading shock wave, a boundary layer and an external inviscid flow in the shock layer) are obtained from the Navier-Stokes equations using an asymptotic analysis. In the case of low Reynolds numbers, the shock layer is considered but the structure of the shock wave is ignored. Together with the well-known models (a boundary layer, a viscous shock layer, a thin viscous shock layer, parabolized Navier-Stokes equations (the single-layer model) for high, moderate and low Re numbers, respectively), a new hydrodynamic model, which follows from the Navier-Stokes equations and reduces to the solution of the simplified (“local”) Stokes equations in a shock layer with vanishing inertial and pressure forces and boundary conditions on the unspecified free boundary (the shock wave) is found at Reynolds numbers, and a density ratio, k, up to and immediately after the leading shock wave, which tend to zero subject to the condition that (k/Re)^1/2 → 0. Unlike in all the models which have been mentioned above, the solution of the problem of the flow around a body in this model gives the free molecular limit for the coefficients of friction, heat transfer and pressure. In particular, the Newtonian limit for the drag is thereby rigorously obtained from the Navier-Stokes equations. At the same time, the Knudsen number, which is governed by the thickness of the shock layer, which vanishes in this model, tends to zero, that is, the conditions for a continuum treatment are satisfied. The structure of the shock wave can be determined both using continuum as well as kinetic models after obtaining the solution in the viscous shock layer for the weak physicochemical processes in the shock wave structure itself. Otherwise, the problem of the shock wave structure and the equations of the viscous shock layer must be jointly solved. The equations for all the continuum models are written in Dorodnitsyn--Lees boundary layer variables, which enables one, prior to solving the problem, to obtain an approximate estimate of second-order effects in boundary-layer theory as a function of Re and the parameter k and to represent all the aerodynamic and thermal characteristic; in the form of a single dependence on Re over the whole range of its variation from zero to infinity.

An efficient numerical method of global iterations, previously developed for solving viscous shock-layer equations, can be used to solve problems of supersonic and hypersonic flows around the windward side of blunt bodies using a single hydrodynamic model of a viscous shock layer for all Re numbers, subject to the condition that the limit (k/Re)^1/2 → 0 is satisfied in the case of small Re numbers. An aerodynamic and thermal calculation using different hydrodynamic models, corresponding to different ranges of variation Re (different types of flow) can thereby, in fact, be replaced by a single calculation using one model for the whole of the trajectory for the descent (entry) of space vehicles and natural cosmic bodies (meteoroids) into the atmosphere.     

Method of successive approximations for three-dimensional laminar boundary layer problems (locally self-similar case) : PMM vol. 37, n6, 1973, pp. 974–983

We present an analytical method for the computation of problems of incompressible boundary layer theory based on an application of the method of successive approximations. The system of equations is reduced to a form suitable for integration. Parameters characterizing the external flow and the body geometry are contained only in the coefficients of the system and do not enter into the boundary conditions. The transformed momentum equations are integrated across the boundary layer from a current value to infinity with the boundary conditions taken into account. If the integration is made from zero to infinity, then the equations pass over into the Kármán relations. Integrating the system of equations a second time, using the boundary conditions at the wall, we obtain a system of nonlinear integro-differential equations. To solve this system of equations we apply the method of successive approximations. To satisfy the boundary Conditions at infinity we introduce, at each step of the iterations, unknown “governing” functions. From the conditions at the outer side of the boundary layer we obtain additional equations for their determination. With the iterational algorithm formulated in this way, the boundary conditions, both on the body and at the outer side of the boundary layer; are satisfied automatically.We consider a locally self-similar approximation. In this case, relative to the “governing” functions, we obtain an algebraic system of equations. We write out the solution in the first approximation. The results obtained in the first approximation are compared with the results of finite-difference computations for a wide range of problems. The results obtained in this paper are compared with those obtained in [1] for the flow in the neighborhood of a stagnation point. An indication is given of the nonuniqueness of the solutions of the three-dimensional boundary layer equations.     

Parabolic flow over a flat plate with wave disturbance in the main stream

In this paper the effects of unsteady flow past a flat plate, which is at rest, are studied when the steady parabolic flow outside the boundary layer is in a wave disturbance. Solutions are developed for large and small times. The behaviour of the amplitudes and phase leads of skin friction components along chordwise and spanwise directions are studied in detail. Two cases, one when the main stream lines are straight and yawed, and the other when the main stream lines are parabolic and concave with respect to the points on the chordwise direction are studied separately and the results are interpreted graphically.     

Effects of roughness on nonlinear stationary vortices in rotating disk flows

Primary instability of rotating disk boundary layer flow over a rough surface for stationary modes was investigated by using the weakly nonlinear theory where the Reynolds number R is close to its critical value R_c as determined by linear theory. Both the single mode case, where the wave vector K equals its critical K_c at the onset of stationary primary instability, and the bimodal case, where the wave vectors K_n (n = 1, 2) are close to K_c for the primary instability of the flow, are considered. The analysis leads to stable solutions for particular roughness forms and magnitude, and particular wave vectors ˜K_n (n = 1, 2) of the surface roughness.     

Wall effects on the slow steady motion of a particle in a viscous incompressible fluid

In this paper, asymptotic expansions with respect to small Reynolds numbers are proved for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid past a single plane wall. The flow problem is modelled by a certain boundary value problem for the stationary, nonlinear Navier-Stokes equations. The coefficients of these expansions are the solutions of various, linear Stokes problems which can be constructed by single layer potentials and corresponding boundary integral equations on the boundary surface of the particle. Furthermore, some asymptotic estimates at small Reynolds numbers are presented for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid between two parallel, plane walls and in an infinitely long, rectilinear cylinder of arbitrary cross section. In the case of the flow problem with a single plane wall, the paper is based on a thesis, which the author has written under the guidance of Professor Dr. Wolfgang L. Wendland.     

Modeling of flow in a very small surface separation

When a fluid flows in a very small surface separation, the very thin boundary layer physically adhering to the solid surface will participate in the flow, while between the two boundary layers is a continuum fluid flow. An analysis is here presented for this multiscale flow. The continuum fluid is treated as Newtonian. The physical adsorbed boundary layer is treated as non-continuum across the layer thickness. The interfacial slippage can occur on the adsorbed layer-solid surface interface, while it is absent on the adsorbed layer-fluid interface. Three flow equations are derived respectively for the two adsorbed layers and the intermediate continuum fluid. They together govern the multiscale flow in such a small surface separation.     

Spatial Stability for Mixed Convection Boundary Layer over a Heated Horizontal Plate

The spatial stability properties of a mixed convection boundary layer developing over a heated horizontal plate is studied here under linear and quasi-parallel flow assumption. The main aim of the present work is to find out if there is a critical buoyancy parameter that would indicate the importance of heat transfer in destabilizing mixed convection boundary layers, when the buoyancy effect is given by Boussinesq approximation. The undisturbed flow used here is that given by the similarity solution of [ 1 ] that implies the wall temperature to vary as the inverse square root of the distance from the leading edge of the plate. The stability of this flow has been investigated by using the compound matrix method (CMM)—that allows finding all the modes in the chosen range in the complex wave number plane for spatial stability analysis. Presented neutral curves for mixed convection boundary layer show the existence of two types of disturbances present simultaneously, for large buoyancy parameter. One notices very unstable high-frequency mode when the buoyancy parameter exceeds the above-mentioned critical value. This unstable thermal mode is in addition to the hydrodynamic mode of isothermal flow given by corresponding similarity profile. The calculated critical buoyancy parameter is shown to qualitatively match with experimental results.     

Group method analysis of melting effect on MHD mixed convection flow from radiate vertical plate embedded in a saturated porous media

The linear transformation group approach is developed to simulate problem of hydromagnetic heat transfer by mixed convection along vertical plate in a liquid saturated porous medium in the presence of melting and thermal radiation effects for opposing external flow. The application of a one-parameter transformation group reduces the number of independent variables by one so that the governing partial differential equations with the boundary conditions reduce to an ordinary differential equations with appropriate corresponding conditions. The Runge-Kutta shooting method is used to solve the determining equations of the set of nonlinear ordinary differential equations. are presented in the form of the temperature and flow fields in the melting region within the boundary layer for different parameters entering into the analysis. Also the effects of the pertinent parameters on the rate of the heat transfer in terms of the local Nusselt number at the solid–liquid interface are also discussed.