Investigations on controlled transition development in a laminar separation bubble by means of LDA and PIV

When a laminar boundary layer separates because of an adverse streamwise pressure gradient, the flow is subject to increased instability with respect to small-amplitude disturbances. Laminar–turbulent transition occurs under a rapid three-dimensional (3D) development within the separated shear layer. When the following turbulent boundary layer reattaches, a laminar separation bubble is formed. To allow controlled measurements, a small-amplitude Tollmien–Schlichting wave (TS wave) was introduced into the boundary layer without (case I) and with (case II) spanwise forcing of steady 3D disturbances. Combined application of laser-Doppler anemometry (LDA) and particle image velocimetry (PIV) demonstrates the suitability of both measurement techniques to capture the development of unsteady, periodic phenomena. The transition mechanism occurring in the flow field under consideration is discussed, and results obtained by controlled measurements are compared to direct numerical simulations (DNS) and predictions from linear stability theory (LST). Flow visualizations and stereoscopic PIV measurements give better insight into the 3D breakdown of the separated shear layer.Nomenclature a amplitude - f₀ fundamental frequency - H₁₂ boundary layer shape factor, H₁₂=₁/₂ - h wavenumber coefficient in time - k wavenumber coefficient in the spanwise direction - l liter - m meter - Re_x Reynolds number based on streamwise distance of the leading edge - Re₁ Reynolds number based on the displacement thickness ₁ - s second - t time - u freestream velocity - u velocity at the boundary layer edge - u,v,w velocities - u,v velocity fluctuations - x streamwise coordinate - y wall-normal coordinate - z spanwise coordinate - boundary layer thickness - ₁ boundary layer displacement thickness - ₂ boundary layer momentum thickness - _z spanwise wavelength - phase angle     

Uniform Convergence for Approximate Traveling Waves in Linear Reaction–Diffusion–Hyperbolic Systems

In this paper we study linear reaction–hyperbolic systems of the form , (i = 1, 2, ..., n) for x > 0, t > 0 coupled to a diffusion equation for p ₀ = p ₀(x, y, θ, t) with “near-equilibrium” initial and boundary data. This problem arises in a model of transport of neurofilaments in axons. The matrix (k _ij ) is assumed to have a unique null vector with positive components summed to 1 and the v _j are arbitrary velocities such that . We prove that as the solution converges to a traveling wave with velocity v and a spreading front, and that the convergence rate in the uniform norm is , for any small positive α.     

Contaminant Transport in Fractured Media with Sources in the Porous Domain

We study contaminant flow with sources in a fractured porous mediumconsisting of a single fracture bounded by a porous matrix. In the fracturewe assume convection, decay, surface adsorption to the interface, and lossto the porous matrix; in the porous matrix we include diffusion, decay,adsorption, and contaminant sources. The model leads to a nonhomogeneous,linear parabolic equation in a quarter-space with a differential equationfor an oblique boundary condition. Ultimately, we study the problemu _t = u _yy – u + f(x,y,t),x,y>0, t>0, u _t = –u _x + u _y – u on y = 0; u(0,0,t) =u₀(t), t>0,with zero initial data. Using Laplace transforms we obtain the Green'sfunction for the problem, and we determine how contaminant sources in theporous media are propagated in time.     

Limiting analytical solutions for complex saturated and unsaturated transport problems

Non-linear diffusion and velocity-dependent dispersion problems are under consideration. The necessary and sufficient conditions allowing the comparison of solutions to the two dimensional convection-dispersion equations with different coefficients are obtained. These conditions provide a framework within which solutions to the complex non-linear problems mentioned above can be estimated by solutions to the problems possessing analytical solvability.Nomenclature c(x, y, t) concentration of solute in solution,ML ^–3 - C(h)=d/dh moisture capacity function - D,D _ij hydrodynamic dispersion coefficient, a second order tensor,L ² T ^–1 - D _L longitudinal hydrodynamic dispersion coefficient,L ² T ^–1 - D _m molecular diffusion coefficient,L ² T ^–1 - D _T transverse hydrodynamic coefficient,L ² T ^–1 - G flow domain for the unsaturated flow problem - G _z , G _w flow domain and complex potential domain, respectively, for the hydrodynamic dispersion problem - h piezometric head,L - I _n given mass flux normal to the boundary,MLT ^–1 - k hydraulic conductivity,LT ^–1 - K(h) unsaturated hydraulic conductivity,LT ^–1 - L continuously differentiable function with respect to all arguments - m porosity - n(x,t) outer normal vector to the boundary - t time,T - V(x, y, t) seepage velocity vector withV=V,LT ^–1 - x Cartesian coordinate system - x horizontal coordinate,L - y vertical coordinate (elevation),L - (x),(x,t) given functions in initial and boundary conditions (3), (4) - ₁(,) angle between vectors ₁c andV - boundary of the flow domain - _L , _T longitudinal and transverse dispersivities, respectively,L - water mass density,ML ^–3 - v _i components of a unit vector in the direction of the outward normal to the boundary - =–kh velocity potential - =/m - stream function defined such thatw=+i is the complex potential - =/m     

The Reissner–Sagoci Problem for a Non-homogeneous Half-space with a Surface Constraint

This paper deals with the problem of twisting a non-homogeneous, isotropic, half-space by rotating a circular part of its boundary surface (0 r < a, z = 0) through a given angle. A ring (a < r < b, z = 0) outside the circle is stress-free and the remaining part (r > b, z = 0) is rigidly clamped. The shear modulus is assumed to vary with the cylindrical coordinates, r, z by the relation (z) = ₁(c + z), c 0 where ₁, c and are real constants. Expressions for some quantities of physical importance, such as torque applied at the surface of the disk and stress intensity factors, are obtained. The effects of non-homogeneity on torque and stress intensity factor are illustrated graphically.     

Singular perturbation of boundary value problem for a vector fourth order nonlinear differential equation

We study the vector boundary value problem with boundary perturbations: ε~2y~((4))=f(x,y,y″,ε, μ) ( μ<χ<1-μ) y(χ,ε,μ)l_(χ-μ)= A_1(ε,μ), y(χ,ε,μ)l_(χ-1-μ)=B_1(ε,μ) y″(χ,ε,μ)l_(χ-μ)=A_2(ε,μ),y″(χ,ε,μ)l_(χ-1-μ)=B_2(ε,μ)where yf, A_j and B_j (j=1,2) are n-dimensional vector functions and ε,μ are two small positive parameters. This vector boundary value problem does not appear to have been studied, although the scalar boundary value problem has been treated. Under appropriate assumptions, using the method of differential inequalities we find a solution of the vector boundary value problem and obtain the uniformly valid asymptotic expansions.     

The Toda System and Clustering Interfaces in the Allen–Cahn equation

We consider the Allen–Cahn equation in a bounded, smooth domain Ω in , under zero Neumann boundary conditions, where is a small parameter. Let Γ₀ be a segment contained in Ω, connecting orthogonally the boundary. Under certain nondegeneracy and nonminimality assumptions for Γ₀, satisfied for instance by the short axis in an ellipse, we construct, for any given N ≥ 1, a solution exhibiting N transition layers whose mutual distances are and which collapse onto Γ₀ as . Asymptotic location of these interfaces is governed by a Toda-type system and yields in the limit broken lines with an angle at a common height and at main order cutting orthogonally the boundary.     

Optimum lifting body shapes in hypersonic flow at high angles of attack

A variational procedure for the determination of lifting body configurations having a maximum lift-to-drag ratio K _max in hypersonic flight at high angles of attack , is proposed. It is based on an analytical solution to the problem for three-dimensional hypersonic flow over small aspect ratio wings using thin shock-layer theory. This reduces the variational problem of finding K _max, and the corresponding optimized wing shape, to the minimization of a linear functional subject to various constraints. The contributions of nonequilibrium thermochemical effects and laminar or turbulent viscous drag effects are also included in the problem formulation. The solution shows that optimized wings have an unbent forward part and a concave lower surface. Due to bifurcation in the optimization process, the planform may have either a sharp apex or a straight nose cut. Corresponding values of K _max() significantly exceed the limiting value K _N=cot for a flat wing. Real thermochemical effects and air viscosity are shown to cause a decrease in K _max and sometimes to influence the optimized wing geometry; however, the relative increment of K _max to K _N is still retained.     

General Nonsymmetric Postbuckling of Orthotropic Annular Plates

ABSTRACT

The general postbuckling behavior of orthotropic annular plates subjected to uniformly distributed forces is studied through the use of the Bubnov-Galerkin method. The transverse diplacement function is assumed to be in the form w = h if /, £ a,,p”+A-’-J cos m8 (=1.1=0 Wheref, = multipliers to be determined, a₁₀= 1, a₁₁(j # 0) = constants depending on the boundary conditions, p, θ = polar coordinates, and m = number of diametral nodal lines. Deflections, bending, and membrane stresses are presented in charts and listed in tables for 12 cases of boundary conditions, 4 loading cases, and various ratios of orthotropy.     

Forced convection heat transfer at the boundary layer of a packed bed

Numerical investigations of the nature of the fluid flow pattern and heat transfer at the boundary layer of a packed bed are reported. A volume averaged Navier-Stokes equation is used to predict the fluid flow and a volume averaged heat balance equation the heat transfer. A variable porosity in the packing is assumed in the region near the wall. Simulations are performed using a modified penalty Galerkin finite element method. The case of fully developed hydrodynamic flow and developing thermal flow is studied. The Nusselt number is found to depend on the Reynolds number, Graetz number and ratio of thermal conductivity of the solid and fluid phases. Comparison is made to some experimental literature values.Nomenclature A constant - [A] Navier-Stokes type matrix - B constant - [B] divergence matrix - C _p constant pressure heat capacity - d characteristic length - D _p particle diameter - D _t tube diameter - {F} solicitation vector - Gz Graetz number, z D _t ^–1 Pr _f Re _p - k permeability term - k _f Thermal conductivity of the fluid phase - k _s Thermal conductivity of the solid phase - [K] coefficient matrix for temperature equation - n normal vector - P pressure - Pr _f Prandtl number for the fluid _f C _p k _f ^-1 - r radial coordinate - R _t tube radius - R residual - R _m residual - Re _p Reynolds number for particle, - t tortuosity factor - T temperature - interstitial velocity - z axial coordinate - effective thermal conductivity - penalty parameter - boundary of solution domain - porosity - viscosity - density - test function - solution domain - test function     

Fluctuating flow in a rotating fluid

Summary Fluctuating flow of a viscous fluid rotating over a disk whose angular velocity oscillates about a nonzero mean is investigated. Initially the disk and the fluid rotate in the same sense with different angular velocities ₁ and ₂ ( ₂> ₁) and at a particular instant of time, the angular velocity of the disk becomes ₁[1+ sin( )]. The problem is solved as an initial boundary value problem and it is found that for small values of the results of analytical and numerical methods are in excellent agreement. The effect of frequency parameter on surface skin frictions has been analysed for various values of angular velocity ratio s and amplitude parameter .

---

Fluktuierende Strömung in einer rotierenden Flüssigkeit

Übersicht Untersucht wird die fluktuierende Strömung einer viskosen Flüssigkeit, die über einer Scheibe, deren Winkelgeschwindigkeit um einen von Null verschiedenen Mittelwert schwankt, rotiert. Anfangs drehen sich die Scheibe und die Flüssigkeit gleichsinnig, aber mit verschiedenen Winkelgeschwindigkeiten ₁ und ₂ ( ₂> ₁). Zu einem Anfangszeitpunkt geht die Winkelgeschwindigkeit der Scheibe über in ₁[1+ sin ( )]. Die Aufgabe wird als Anfangs-/Randwertproblem gelöst. Für kleine Werte stimmen die analytischen und numerischen Ergebnisse hervorragend überein. Für verschiedene Werte des Winkelgeschwindigkeitsverhältnisses und des Amplitudenparameters wurde der Einfluß des Frequenzparameters auf die Reibspannungen an der Scheibe untersucht.

     

Heat transfer by rotational flow in an annulus with porous lining

The flow and heat transfer in an annulus between rotating coaxial cylinders, with non-erodible porous lining, is investigated. The flow in the porous lining is obtained by using Brinkman equation. At the boundary between the porous lining and the free flow (the so called nominal surface), the velocity slip and the temperature slip are used. A quasi-numerical technique developed by the authors is employed in obtaining the solution of the energy equation. The effect of the thickness of the porous lining and the permeability on the velocity and the Nusselt numbers at the walls is studied.

---

Wärmeübergang bei rotierender Strömung in einem Ring mit poröser Wand

Zusammenfassung In dieser Arbeit wird die Strömung und Wärmeübertragung zwischen rotierenden koaxialen Zylindern mit unauswaschbarem porösem Überzug untersucht. Die Strömung innerhalb des porösen Überzugs ist mit Hilfe der Brinkmanschen Gleichung berechnet. An der Grenze (der sogenannten Nominalfläche) zwischen dem Überzug und der freien Strömung wurde die Geschwindigkeitsgleitung und Temperaturgleitung benutzt. Die Energiegleichung ist mit Hilfe eines von den Autoren entwickelten quasi-numerischen Verfahrens gelöst. Der Einfluß der Dicke und der Durchlässigkeit des porösen Überzugs auf die Strömung und die Nusseltschen Zahlen an den Wänden wird untersucht.

---

Nomenclature R ₂ radius of the outer cylinder forming the annulus - ₂ angular velocity of the outer cylinder - T ₂ temperature of the outer cylinder - R _l radius of the inner cylinder forming the annulus - ₁ angular velocity of the inner cylinder - T ₁ temperature of the inner cylinder - h thickness of the porous lining - R radial distance of any point in the annulus - V azimuthal component of velocity in zone 1 (of Fig. 1) - V part of velocity in zone 2 (of Fig. 1) due to transfer of momentum from the main flow - V _p velocity in zone 2 (of Fig. 1) - Q Darcy velocity in the porous medium (zone 2 of Fig. 1) - velocity slip parameter - k absolute permeability of the material used for lining - ₀ nondimensional shearing stress at the outer cylinder - _i nondimensional shearing stress at the inner cylinder - K thermal conductivity in zones 1 and 2 (of Fig. 1) - coefficient of viscosity of the fluid - T temperature in zone 1 (of Fig. 1) - T temperature in zone 2 (of Fig. 1) - temperature slip parameter - (Nu) _o nondimensional Nusselt number at the outer cylinder - (Nu) _i nondimensional Nusselt number at the inner cylinder - radii ratio - nondimensional rotational parameter - nondimensional thickness of the porous lining     

Temperature distribution in a thin rotating disk

The problem of heat conduction in a thin rotating disk with heat input at a fixed point is considered. The disk is cooled by forced convection from its lateral surfaces. By defining a complex temperature, the temperature throughout the disk is presented as a series of Bessel functions of complex argument. Results are given for a range of rotational speeds.Nomenclature R radial coordinate - angular coordinate - a radius of disk - b thickness of disk - T temperature - T ambient temperature - rotational speed of disk - q heat flux into disk - k thermal conductivity of disk - density of disk - c specific heat of disk - h coefficient of convective heat transfer - r dimensionless radial coordinate, R/a - T* characteristic temperature, q ₀ a/ k - t dimensionless temperature, (T–T )/T* - C ₁, C ₂ dimensionless parameters defined in (3)     

Turbulent boundary layer over a smooth wall with widely separated transverse square cavities

Laser-Doppler velocimetry (LDV) measurements and flow visualizations are used to study a turbulent boundary layer over a smooth wall with transverse square cavities at two values of the momentum thickness Reynolds number (R =400 and 1300). The cavities are spaced 20 element widths apart in the streamwise direction. Flow visualizations reveal a significant communication between the cavities and the overlying shear layer, with frequent inflows and ejections of fluid to and from cavities. There is evidence to suggest that quasi-streamwise near-wall vortices are responsible for the ejections of fluid out of the cavities. The wall shear stress, which is measured accurately, increases sharply immediately downstream of the cavity. This increase is followed by a sudden decrease and a slower return to the smooth wall value. Integration of the wall shear stress in the streamwise direction indicates that there is an increase in drag of 3.4% at bothR .Nomenclature C _f skin friction coefficient - C _fsw friction coefficient for a continuous smooth wall - k height of the cavity - k ⁺ ku / - R Reynolds number based on momentum thickness (U ₁ /v) - Rx Reynolds number based on streamwise distance (U ₁ x/) - s streamwise distance between two cavities - t time - t ⁺ tu ₂ / - U ₁ freestream velocity - mean velocity inx direction - u,v,w rms turbulent intensities inx,y andz directions - u local friction velocity - u _sw friction velocity for a continuous smooth wall - w width of the cavity - x streamwise co-ordinate measured from the downstream edge of the cavity - y co-ordinate normal to the wall - z spanwise co-ordinate - y ⁺ yu / - boundary layer thickness - ₀ boundary layer thickness near the upstream edge of the cavity - _i thickness of internal layer - kinematic viscosity of water - ⁺ zu / - momentum thickness     

A note on three-dimensional boundary layer growth

Summary This note is an extension of the work of Görtler²) on two-dimensional boundary layer growth to the three-dimensional case. The solutions of three-dimensional boundary layer equations are obtained by considering the potential flow of the body to be governed by the functions At U ₀(, ) and At U ₀(, ) where is any positive number.     

Calculation of the boundary layer on the electrically conducting wall of a plane channel

There are presently available quite a large number of works devoted to the study of the motion of an electrically conducting fluid in boundary layers formed on electrodes or on the nonconducting walls of various MHD devices. However, the methods of solving the boundary layer equations in these studies are based on various simplifying assumptions which allow the problem to be reduced to the solution of a system of ordinary differential equations. Thus, in [1] there is imposed on the flow the special magnetic fieldH1/x, which enables the problem to be reduced to the self-similar form, while in the studies of other authors [2, 3] either the solution is sought in the form of expansions in x, or it is assumed that the problem is locally self-similar [4]. In the present paper we construct the solution of the MHD boundary layer equations which is obtained by one of the numerical methods which has long been used for solving the boundary layer equations for a nonconducting fluid.     

Magnetohydrodynamic natural convection heat transfer from radiate vertical porous surfaces

Magnetohydrodynamic natural convection heat transfer from radiate vertical surfaces with fluid suction or injection is considered. The nonsimilarity parameter is found to be the conductive fluid injection or suction along the streamwise coordinate = V{4x/² g(T _w – T )}^1/4. Three dimensionless parameters had been found to describe the problem: the magnetic influence number N = B ² _y /V ², the radiation-conduction parameter R _d = k _R /4aT ³ , and the Gebhart number Ge _x = gx/cp to represent the effect of the viscous dissipation. It is found that increasing the magnetic field strength causes the velocity and the heat transfer rates inside the boundary layer to decrease. Its apparent that increasing the radiation-conduction parameter decreases the velocity and enhances the heat transfer rates. The Gebhart number, i.e, the viscous dissipation had no effect on the present problem.Nomenclature a Stefan-Boltzmann constant - B _y Magnetic field flux density Wb/m² - Cf _x Local skin friction factor - c _p Specific heat capacity - f Dimensionless stream function - Ge _x Gebhart number, gx/cp - g Gravitational acceleration - k Thermal Conductivity - L Length of the plate - N Magnetic influence number, B ² _y /V ² - p Pressure - Pr Prandtl number - q _r Radiative heat flux - q _w (x) Local surface heat flux - Q _w (x) Dimensionless Local surface heat flux - R _d Planck number (Radiation-Conduction parameter), k _R /4aT ³ - T Temperature - T Free stream temperature - T _w Wall temperature - u, v Velocity components in x- and y-directions - V Porous wall suction or injection velocity - V _w Porous wall suction or injection velocity - x, y Axial and normal coordinates - Thermal diffusivity Greek symbols _R Roseland mean absorption coefficient, 4/3R _d - Coefficient of thermal expansion - Nonsimilarity parameter, V{4x/² g(T _w – T )}^1/4 - Peseudo-similarity variable - Dimensionless temperature - _w Ratio of surface temperature to the ambient temperature, T _w /T - Dynamice viscosity - Kinemtic viscosity - Fluid density - Electrical conductivity - _w Local wall shear stress - Dimensional stream function     

On the laminar forced convection with axial conduction in a circular tube with exponential wall heat flux

The values of the fully developed Nusselt number for laminar forced convection in a circular tube with axial conduction in the fluid and exponential wall heat flux are determined analytically. Moreover, the distinction between the concepts of bulk temperature and mixing-cup temperature, at low values of the Peclet number, is pointed out. Finally it is shown that, if the Nusselt number is defined with respect to the mixing-cup temperature, then the boundary condition of exponentially varying wall heat flux includes as particular cases the boundary conditions of uniform wall temperature and of convection with an external fluid.

---

Über laminare Zwangskonvektion mit Längswärmeleitung in einem Kreisrohr mit exponentiell veränderlichem Wandwärmefluß

Zusammenfassung Es werden die Endwerte der Nusselt-Zahlen für vollausgebildete laminare Zwangskonvektion in einem Kreisrohr mit Längswärmeleitung und exponentiell veränderlichem Wandwärmefluß analytisch ermittelt. Besondere Betonung liegt auf dem Unterschied zwischen den Konzepten für die Mittel- und die Mischtemperatur bei niedrigen Peclet-Zahlen. Schließlich wird gezeigt, daß bei Definition der Nusselt-Zahl bezüglich der Mischtemperatur die Randbedingung exponentiell veränderlichen Randwärmeflusses die Spezialfälle konstanter Wandtemperatur und konvektiven Wärmeaustausches mit einem umgebenden Fluid einschließt.

---

Nomenclature A _n dimensionless coefficients employed in the Appendix - Bi Biot numberBi=h _e r ₀/ - c _n dimensionless coefficients defined in Eq. (17) - c _p specific heat at constant pressure of the fluid within the tube, [J kg^–1 K^–1] - f solution of Eq. (15) - h ₁,h ₂ specific enthalpies employed in Eqs. (2) and (4), [J kg^–1] - h _e convection coefficient with a fluid outside the tube, [W m^–2 K^–1] - rate of mass flow, [kg s^–1] - Nu bulk Nusselt number,2r ₀ q _w /[(T _w –T _b )] - Nu _H fully developed value of the bulk Nusselt number for the boundary condition of uniform wall heat flux - Nu _T fully developed value of the bulk Nusselt number for the boundary condition of uniform wall temperature - Nu ^* mixing Nusselt number,2r ₀ q _w /[(T _w –T _m )] - Nu _C ^* fully developed value of the mixing Nusselt number for the boundary condition of convection with an external fluid - Nu _H ^* fully developed value of the mixing Nusselt number for the boundary condition of uniform wall heat flux - Nu _T ^* fully developed value of the mixing Nusselt number for the boundary condition of uniform wall temperature - Pe Peclet number, 2r ₀/ - q ₀ wall heat flux atx=0, [W m^–2] - q _w wall heat flux, [W m^–2] - r radial coordinate, [m] - r ₀ radius of the tube, [m] - s dimensionless radius,s=r/r ₀ - T temperature, [K] - T ₀ temperature constant employed in Eq. (14), [K] - T reference temperature of the fluid external to the tube, [K] - T _b bulk temperature, [K] - T _m mixing or mixing-cup temperature, [K] - T _w wall temperature, [K] - u velocity component in the axial direction, [m s^–1] - mean value ofu, [m s^–1] - x axial coordinate, [m] Greek symbols thermal diffusivity of the fluid within the tube, [m² s^–1] - exponent in wall heat flux variation, [m^–1] - dimensionless parameter - dimensionless temperature =(T _w –T)/(T _w –T _b ) - ^* dimensionless temperature ^*=(T _w –T)/(T _w –T _m ) - thermal conductivity of the fluid within the tube, [W m^–1 K^–1] - density of the fluid within the tube, [kg m^–3]     

Laminar boundary layer on blunt bodies with account for vorticity of the external flow

This paper presents the technique for and results from numerical calculations of the hypersonic laminar boundary layer on blunted cones with account for the vorticity of the external flow caused by the curved bow shock wave. It is assumed that the air in the boundary layer is in the equilibrium dissociated state, but the Prandtl number is assumed constant, =0.72. The calculations were made in the range of velocities 3–8 km/sec, cone half-angles _k=0°–20°. With account for the vortical interaction of the boundary layer with the external flow, the distribution of the thermal flux and friction will depend on the freestream Reynolds number (other conditions being the same). In the calculations the Reynolds number R, calculated from the freestream parameters and the radius of the spherical blunting, varies from 2.5·10³ to 5.10⁴. For the smaller Reynolds numbers the boundary layer thickness on the blunting becomes comparable with the shock standoff, and for R<2.5·10³ it is apparent that we must reconsider the calculation scheme. With R>5·10⁴ for cones which are not very long the vortical interaction becomes relatively unimportant. The results of the calculations are processed in accordance with the similarity criteria for hypersonic viscous gas flow past slender blunted cones [1, 2].     

Jet methods of gas injection into fluid boundary layer for drag reduction

Two jet methods for saturating the fluid boundary layer with microbubbles for drag reduction in contrust with gas injection through porous materials are considered. The first method is the gas injection through the slot under a special fluid wall jet. The second method is the saturation of boundary layer by microbubbles via the gas-water mixture injection through the slot. Experimental data, reflecting the skin friction drag reduction on the flat plate and total drag reduction of axisymmetric bodies, are presented. The comparison between a jet methods of gas injection and gas injection through porous materials is made.Nomenclature v free-stream velocity - v _j mean velocity of a water through slot - v _g mean velocity of a gas through slot - h width of slot for realizing water jet - h ₁ width of slot for gas injection - incidence angle - Q volume airflow rate - C _Q airflow rate coefficient (v _g/v ) - C _f skin friction coefficient - v _j/v - C _f0 C _f ifQ=0 andv _j=0 - _f C _f/C _{f 0} - d diameter of an axisymmetric body - L length of body - C _Q 4 · ·Q/d ² v - C _D 4 ·D/1/2v ² ·d ² - C _Q 4 ·Q/d ² v - Q _j volume flow rate of water jet - C 8 ·Q _jv_j/d ² v ² - 1 fluid density of main flow - 2 fluid density of wall jet - B ₁ main stream total pressure - B ₂ wall jet total pressure - v ₁ main stream velocity - Be (B ₂ –B ₁)/1/2₁ v ₁ ² = Bernoulli number - ₂ v ₂/₁ v ₁ - p _st static pressure - p _at atmospheric pressure - p _st/p _at - D hydrodynamic drag of body