Next Article in Journal
Freestream Turbulence Effects on the Aerodynamics of an Oscillating Square Cylinder at the Resonant Frequency
Previous Article in Journal
Hamiltonian Variational Formulation of Three-Dimensional, Rotational Free-Surface Flows, with a Moving Seabed, in the Eulerian Description
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Gravitational Force on Particle Motion in the Channel Flow Induced by Fluid Injection

by
Konstantin Volkov
Institute of Hydrodynamics and Control Processes, St Petersburg State Marine Technical University, 190121 St Petersburg, Russia
Fluids 2022, 7(10), 328; https://doi.org/10.3390/fluids7100328
Submission received: 8 September 2022 / Revised: 4 October 2022 / Accepted: 9 October 2022 / Published: 14 October 2022
(This article belongs to the Special Issue Multiphase Flow in Pipes with and without Porous Media, Volume II)

Abstract

:
The motion of an individual particle in a circular channel flow induced with fluid injection is considered. Analysis takes into account drag and gravitational forces acting on an individual particle. The change in the radial structure of the flow in a channel with fluid injection on the motion of a particle is studied. A number of simplifying assumptions about the structure of the fluidflow in the channel makes it possible to obtain an analytical solution of the problem for particles. The results of a qualitative analysis of particle trajectories in the channel with fluid injection are compared with numerical simulations. The singular points of the particle trajectory are found in a wide range of characteristic non-dimensional parameters of the problem. The modes of motion of an individual particle in a channel with fluid injection are classified depending on the Stokes and Froude numbers.

1. Introduction

In the products of solid propellant combustion, along with the gas-phase component, there are solid and liquid combustion products of metal powders [1,2]. Simulation of two-phase flows of combustion products plays an important role in the study of intra-chamber processes in solid propellant rocket motors (SRM). The addition of solid particles to the flow complicates the flow pattern, which is associated with a variety of properties (inertia and concentration) of dispersed particles, which lead to the implementation of numerous flow regimes [3]. The concentration of particles is one of the main physical parameters that determine the characteristics of the movement of the dispersed phase, as well as the degree of its influence on the fluidflow [4]. In a dispersed medium, zones free of particles (fragmentation of the phase volume) and regions with intersecting particle trajectories can appear, at the boundaries of which the concentration of the dispersed phase increases sharply [5].
Analytical solutions of the Euler or Navier–Stokes equations contribute to a better understanding of the qualitative features of stationary and non-stationary flows. They allow estimating the domain of applicability of simplified mathematical models (inviscid fluid, slow flows, boundary layers) and are indispensable for validation of numerical and approximate methods. In the approaches developed in [6,7,8], self-similar solutions of the continuity equation and the momentum equations are found under the assumption of a linear dependence of the axial velocity on one of the spatial coordinates and uniform injection from the channel walls. An experimental validation of the analytical solution corresponding to the vortical flow of an inviscid incompressible fluid is carried out in [9,10,11,12], in which measurement data were obtained from the distributions of axial velocity and pressure at various Reynolds numbers.
The influence of compressibility increases at a sufficiently large distance from the inlet section of the channel, which leads to the filling of the axial velocity profile [12,13]. Turbulence has a rather weak effect on the structure of the mean flow, with the exception of the region adjacent to the outlet section of the channel, where there is a significant acceleration of the flow, and the mean velocity profile becomes flatter [11]. The movement of the burning surface of the channel (the effect of non-stationarity of the process) has a comparatively weak effect on the velocity and pressure distributions [14]. The presence of a laminar flow region, which occupies the entire flow region (depending on the Reynolds number), is confirmed by [15,16].
To find a class of analytical solutions describing the flow of a viscous incompressible fluid in channels induced by fluid injection, methods of group analysis [17] are applied. The stability of the flow in a semi-infinite round pipe with a closed left end is considered in [18].
To take into account the influence of particles on the fluid flow, a joint calculation of the motion of the fluid and dispersed phases is required. Under certain restrictions on the particle size (the Stokes law for the particle drag) and the geometry of the computational domain (plane or axisymmetric channel), it is possible to construct analytical or approximate solutions that are convenient for making qualitative estimates and verifying numerical calculations. The mode of single particles is one of the simplest modes of particle-laden flow, which is realized at a low concentration of the dispersed phase. Obtaining new solutions is associated with the use of new dependencies describing the fluid flow and taking into account the influence of additional forces. Analytical or approximate solutions describing the motion of particles make it possible to find the limiting trajectory of particles (separatrix), to study the features of the concentration of the dispersed phase and their deposition on the wall [19,20].
Investigations of flows with particles in combustion chambers are related to the modeling of slagging of sections of the channel [21,22], combustion of aluminum particles [23,24], transport of particles by vortical structures and their interaction with an oscillating flowfield [25,26], and combustion instability [27,28]. The particle motion in channels with fluid injection and their influence on the characteristics of the carrier flow are studied in [29,30,31,32]. The motion of a particle in a propellant channel is considered in [33,34], in which an approach that simplifies the implementation of mathematical models of two-phase flows is proposed.
This study presents analytical solutions of equations describing the motion of an individual particle in a channel with fluid injection from walls, as well as the results of a numerical simulation of particle motion under the influence of non-turbulent factors. The influence of gravitational force on the motion of a particle is discussed and the domains of their applicability are investigated. The combustion surface is modeled by the injection surface of a mixture of fluid and particles, the parameters of which are known and do not depend on the place of injection and the presence of body forces.

2. Flow Induced by Fluid Injection

A flow in a circular channel, formed by injection of a mixture of fluid and dispersed particles with given parameters along the normal to the channel surface is considered. Flow is in a potential field of body forces, the vector of which is normal to the centerline of the channel. This model describes a gas-particle flow in a horizontal channel in the presence of gravity. It is assumed that the gravitational effects do not altered flowfield. The affect particle motion only.
A quasi-developed flow of a viscous incompressible fluid in a sufficiently long cylindrical channel with fluid injection is simulated, when the flow quantities, normalized on the velocity on the channel centerline, change slightly along its length. Fluid injection is uniformly distributed along cylindrical wall. The x axis is aligned with the centerline of the channel, and the y axis is normal to it (Figure 1). The radius of the channel is h, and the injection velocity v w is the same at all points of the channel surface and directed along the normal to it. The spreading of the fluid occurs symmetrically with respect to plane x = 0 (mirror symmetry of the flow).
The channel radius h is chosen as characteristic scale for variables with the dimension of length, and the injection velocity v w is chosen as characteristic scale for variables with the dimension of velocity. The characteristic time is introduced as the ratio of the channel radius to the injection velocity, h / v w . Potential body forces do not affect the fluid velocity in the channel [20], therefore the fluid velocity field is described by the relations that take place in the vortical flow of an inviscid incompressible fluid
u v w = π x h cos π 2 r h 2 , v v w = h r sin π 2 r h 2 .
In this case, the influence of viscous effects manifests itself in rather narrow channels, while the influence of compressibility should be taken into account in long channels (Figure 2). The solid lines correspond to the results of numerical simulation [20], the dotted lines correspond to the vortical flow of an inviscid incompressible fluid, the symbols • correspond to experimental data. Profiles of axial velocity predicted with the developed model (solid lines) are flatter in the core region and steeper in the near wall region than those computed in vortical flow (dashed lines). The distributions of radial velocity computed with turbulence model and model of vortical flow are in a good agreement. Compressibility effects play an important role in a long channel leading to a flatter profile of axial velocity. Impact of Reynolds number on profiles of axial and radial velocities in flow induced by wall injection is small.
The influence of viscous effects on the flow structure in a channel with double-sided injection manifests itself mainly in the axial region, changing the fullness of the axial velocity profile (Figure 3). In a circular channel, as the injection intensity increases, the axial velocity profile becomes less filled and tends to a cosine profile (1) at Re . The radial velocity profiles differ relatively little from each other in a wide range of injection velocities. The maximum radial velocity is located at some distance from the channel wall and exceeds the injection velocity.
Particles are non-deformable spheres of the same diameter; their collisions and their effect on the fluid are neglected. In the model of interaction between a particle and a carrier flow, the drag force and gravity force are taken into account. The drag coefficient of a particle is found from the Stokes law.
The appearance of body forces complicates the problem, since the flow becomes three-dimensional. In this case, the flow parameters in the radial plane of the channel do not depend on the axial coordinate x (Figure 4).

3. Governing Equations

In the discrete-trajectory approach, the equations describing the motion of a particle are written in Lagrangian variables and integrated along the trajectories of individual particles in a known (computed in advance) fluidfield.
The Equation describing the translational motion of a spherical particle has the form
m p d v p d t = 1 2 C D ρ S m v v p 2 + f p ,
where m p is particle mass, S m is particle midsection area, v and v p are fluid velocity and particle velocity, respectively. The term f p in Equation (2) includes forces of different nature, but different from the drag force, the account or disregard of which depends on the specific problem (for example, the gravity force is f p = m p g ). The drag coefficient is represented as
C D = 24 Re p f D ( Re p ) .
The function f D takes into account the correction for the inertia of the particle. The Reynolds number in the relative motion of a particle and fluid is found by the formula
Re p = 2 r p ρ v v p μ ,
where r p is particle radius.
A kinematic relation is added to Equation (2), which makes it possible to calculate the radius vector of the particle’s center of mass
d r p d t = v p .
The temperature of a particle affects its motion through a correction to the drag coefficient. In many flow regimes, such a correction is small and is not taken into account.
Equations (2) and (3) are integrated along the particle trajectory and require the initial conditions (the coordinates and translational velocity of the particle at time t = 0 ). The fluid velocity is calculated at the points lying on the particle trajectory, v = v ( r p ) .

4. Transformation of Coordinates

The motion of the particle is described in the Cartesian coordinate system. The Equations describing the motion of a particle have the form
d u p d t = s ( r ) x p u p Stk ;
d v p d t = p ( r ) y p v p Stk ;
d w p d t = p ( r ) z p w p Stk Fr .
Here, r 2 = y p 2 + z p 2 . The functions s ( r ) and p ( r ) in (1) are found from the relations
s ( r ) = π cos π 2 r 2 , p ( r ) = 1 r 2 sin π 2 r 2 .
Equations (4)–(6) are supplemented with kinematic relations that allow calculating the coordinates of the particle
d x p d t = u p ;
d y p d t = v p ;
d z p d t = w p .
At the initial time t = 0 the particle is on the channel wall, therefore
x p 0 = x 0 , y p 0 = cos α , z p 0 = sin α .
The particle is injected from the wall along the normal to the surface, which gives
u p 0 = 0 , v p 0 = v 0 y p 0 , w p 0 = v 0 z p 0 .
Here, α and v 0 are the initial coordinate and initial velocity of the particle in the cylindrical coordinates (angle and velocity on the injection surface), and α [ 0 , 2 π ] and v 0 [ 0 , 1 ] .
The motion of a particle in a channel is determined by the Stokes and modified Froude numbers
Stk = ρ p d p 2 v w 18 μ h , Fr = h g v w 2 ,
where g is gravitational acceleration.

5. Qualitative Analysis

The radial structure of a two-phase flow in a channel is determined by the system of Equations (4)–(9) and the corresponding initial conditions. Qualitative methods for the analysis of dynamical systems make it possible to systematize and generalize the solutions of the problem under various initial conditions and dimensionless parameters.

5.1. Zone of Gradient Flow

Near the channel centerline ( r 0 , in practice for r < 0.5 ) functions are replaced by linear relations
s ( r ) = π , p ( r ) = π 2 .
For small radii, the radial velocity gradient is equal to π / 2 and does not depend on the radial coordinate.
In the gradient flow zone, the system of Equations (4)–(9) becomes linear
Stk d 2 x p d t 2 + d x p d t + π x p = 0 ;
Stk d 2 y p d t 2 + d y p d t + π 2 y p = 0 ;
Stk d 2 z p d t 2 + d z p d t + π 2 z p + Sf = 0 .
The parameter Sf is the product of the Stokes number and the Froude number ( Sf = Re Fr ).
Equations (10)–(12) are integrated independently of each other. A new variable is introduced
z ^ p = z p + 2 π Sf .
Then, Equations (11) and (12), describing the motion of a particle in the cross section of the channel, take the form
Stk d 2 y p d t 2 + d y p d t + π 2 y p = 0 ;
Stk d 2 z ^ p d t 2 + d z ^ p d t + π 2 z ^ p = 0 .
At the initial time z ^ p 0 = z p 0 + 2 Sf / π .
Equations (13) and (14) have the same form and are independent of each other. The introduction of a new variable makes it possible to eliminate the dependence on the Stokes and Froude numbers (parameter Sf). In the coordinate system ( y p , z ^ p ) , the flow in the gradient zone is axisymmetric and does not depend on body forces.
The characteristic equation for Equations (13) and (14) has the form
Stk λ 2 + λ + π 2 = 0 .
Solving a quadratic equation, find its roots
λ 1 , 2 = 1 ± 1 2 π Stk 1 / 2 2 Stk .
The roots of the characteristic equation show that there is a critical Stokes number Stk * = 1 / 2 π , which determines the change in qualitatively different two-phase flow regimes. For Stk < Stk * , the particle does not cross the channel centerline (particle trajectories do not qualitatively differ from streamlines). For Stk > Stk * there is an intersection of the particle trajectory with the channel centerline.
The form of particle trajectories on the plane ( y p , z ^ p ) with an arbitrary parameter Sf is obtained from the special case corresponding to Sf = 0 . Particle trajectories are straight lines if the initial velocity vector of the particle at the boundary of the gradient zone is directed to a singular point located at the origin of the coordinate system (Figure 5), which requires that the condition
u p 0 v p 0 = z ^ p 0 y p 0 .
Index 0 refers to the coordinates and initial velocity of the particle at the boundary of the gradient zone.
The particle trajectories obtained from the solution of Equations (13) and (14), under the specified condition and any Stokes number, are straight lines directed from the initial point on the injection surface to the singular point at the center of the channel. The oscillatory motion in the neighborhood of the singular point, which occurs at Stk > Stk * , is along this straight line.
From a physical point of view, a singular point is a projection of the corresponding asymptote of the particle trajectories in the longitudinal section of the channel onto the radial plane. The singular point characterizes the limiting state of the particle at an unlimited time of its motion. The type of singular point depends on the Stokes number and its critical value.
At Stk < Stk * the particle trajectories monotonically approach the singular point (stable node). Particle trajectories are similar to streamlines. At Stk > Stk * , an oscillatory regime (with damping) of particle motion along a straight line relative to a singular point (stable focus) is observed. Under these conditions, a specific flow zone is formed in the channel with its characteristic counter-oscillatory motion of particles. As the Stokes number increases, this zone expands. Since the concentration of particles in this zone increases, and the probability of their collision increases, to study the flow region with counter-oscillatory motion of particles, more advanced models are required, including, in particular, the effects of interaction of particles and coupling effects.
In the absence of body forces ( Sf = 0 ), the singular point is at the center of the channel (the flow is axisymmetric). In the presence of body forces, the singular point is displaced in the direction of the vector of these forces. Approximately at Sf < 0.7 , the singular point is in the region of the gradient fluid flow. The particle trajectories in the radial section of the channel, as in the case of Sf = 0 , are still straight lines and are directed to a singular point, the type of which does not change within the gradient zone. The transition to the original coordinate system is graphically carried out by its shift along the z axis by 2 Sf / π (Figure 5).

5.2. General Case

In the range 0.7 < Sf < 1 , Equations (4)–(9) describing particle motion are non-linear, which complicates their analysis and leads to the need to use numerical methods of analysis.
The gradient zone occupies the central, but smaller in terms of cross-sectional area, part of the channel. Outside this zone, the analysis is carried out on the basis of the original model described by Equations (10)–(12).
With the accepted direction of the gravity vector, the coordinates of singular points are found from the conditions
α * = π 2 , 1 r * sin π 2 r * = Sf .
Singular solution points exist only on the z axis in the lower part of the channel.
For small radii (for r 0 ), an approximate relation for the singular point coordinate is obtained, which is valid for the gradient flow zone
r * = 2 π Sf .
The dependence of the radial coordinate of the singular point r * on the parameter Sf is shown in Figure 6. Line 1 corresponds to the numerical solution, and line 2 corresponds to an approximate relation valid for the gradient zone of the flow.
For Sf = 0 , there is one singular point, which is located at the origin of the coordinate system. As Sf increases, the singular point moves along line 1 obtained by the numerical solution. For small r and Sf, line 1 is approximated with good accuracy by line 2. The position of the outer boundaries of the gradient zone is found at r 0.7 , which takes place at Sf 0.7 .
For Sf = 1 , the second singular point appears at r * = 1 (on the injection surface). As Sf increases further, the singular points approach each other (along line 1). For Sf = 1.07 and r * = 0.86 these singular points merge and disappear. For Sf > 1.07 , there are no singular points.
A connection between the presence and position of singular points on the phase plane and the kinematic structure of a two-phase flow in a channel is analyzed. On the phase plane, at the origin of the coordinate system ( r p , v p ) , there is a single singular point of the node type (particles do not cross the channel centerline) or focus (particles cross the channel centerline) depending on the Stokes number.
When body forces are taken into account, the singular point is displaced from the origin of the coordinate system down along the z axis. Changes in the structure of a two-phase flow in the longitudinal section of the channel are schematically shown in Figure 7 (at Stk = 0.1 ). For each Froude number, two trajectories of particles injected from the upper ( z = 1 ) and lower ( z = 1 ) surfaces of the channel are presented.
As the Froude number increases, the asymptote of the axial (along the x axis) particle trajectories, defined as z = r * , shifts downward. At Sf = 1 the second singular point appears (such as a saddle). This state characterizes a new quality of the system, when the drag force acting on the particle from the side of the fluid at the bottom of the channel are balanced by the mass force in opposite direction. This state is dynamically unstable, and with a slight increase in the initial velocity of the particles, they are carried away into the flow. At Sf 1.07 , which corresponds to the extreme value of the radial fluid velocity in the channel at r * 0.86 , both singular points merge. In this case, the particles from the upper half of the channel reach the asymptote of the trajectories ( r * 0.86 ), while the particles on its lower side are held by body forces. The zone free from particles of the corresponding size that appeared at Sf = 1 reaches its maximum size. For Sf > 1.07 there are no singular points. This leads to a new qualitative change in the flow regimes when the asymptote of the trajectories and the zone without particles disappear. Particles in the lower half of the channel cannot leave the blowing surface, and particles leaving the upper half of the channel also fall here (given sufficient time of motion). A further increase in the Stokes and Froude numbers (parameter Sf) does not cause qualitative changes in the flow.

6. Parameter Expansion Method

A modeling of two-phase flows in channels with fluid injection is analyzed with the parameter expansion method.

6.1. Velocity of Non-Equilibrium Motion

The motion of a Stokes particle of variable size, taking into account the body forces, is described by Equation (2), which is written in a dimensionless form
d v p d t = B v g v p + Fr ,
where B = 1 / ( Stk δ 2 ) , δ = r p / r p 0 . Index 0 refers to the start time.
The general solution of Equation (15), which is a linear ordinary differential Equation of the 1st order, is presented as the sum of the general solution of the uniform equation and a partial solution of the non-uniform Equation, v p = v ^ p + v ˜ p . The general solution of the uniform equation has the form
v ^ p = C 1 exp F .
To find a partial solution of non-uniform equation, the method of expansion into a series in terms of the parameter is used
v ˜ p = v g + k = 1 ( 1 ) k Stk k d k v g d t k .
Restricting the expansion to terms of the first order of smallness, the solution of Equation (15) takes the form
v p = C 1 exp F + v g + Stk Fr v g t δ 2 .
Here
F = 1 Stk 0 t d τ δ 2 ( τ ) , C 1 = v p 0 v g Stk Fr v g t 0 .
The relation (16) shows that at Stk 0 the two-phase flow tends to equilibrium ( v p v g ). For Stk 0 and r p 0 , flow reaches equilibrium state, v p v g . At the final stage of evolution, the particles also come into a state of dynamic equilibrium with the fluid. For δ = 1 , a solution for a particle of constant size is obtained ( r p = r p 0 ).
The particle coordinates are found by resolving the kinematic relation
d r p d t = v p .
Integration of Equation (17) over time from 0 to t, for δ 1 gives
r p = C 2 C 1 Stk exp t / Stk + v g t + Stk Fr t v g .
Since the time in the relation (18) enters the term with the Stokes number as a factor, and the role of the exponential term decreases with decreasing Stokes number, the time of the equilibrium phase motion is used as the time.
The integration constants, C 1 and C 2 , are determined from the initial conditions. At time t = 0 , the particle is on the channel wall, therefore
x p ( 0 ) = x p 0 , y p ( 0 ) = 1 , u p ( 0 ) = 0 , v p ( 0 ) = ϕ ,
where ϕ is the parameter of the initial velocity non-equilibrium of the phases.
The obtained solution is linear with respect to the Stokes number and satisfactorily describes the particle trajectories in short channels. For Stk it qualitatively incorrectly reflects the behavior of the particle velocity.
Let us assume that at t t b , where t b is the burning time of a particle, the change in particle size is described by the Equation
r p = r p 0 1 t / t b 1 / q .
For q = 2 , which corresponds to the theoretical model of droplet combustion, the particle velocity is found from the relation
v p = C 1 F t 1 / ( H Stk ) + v g + 1 F t Fr v g t Stk .
The motion of a particle of variable size is determined by the parameter H = t p / t b , which is equal to the ratio of the time the particle stays in the channel t p = L / v w and the burning time of the particle t b = d p 0 q / k . The parameter characterizing the combustion of aluminum is represented as
H = k L v w d p 0 q .
As H increases, the conditions for particle combustion in the channel improve. For H = the particle burns on the channel surface, and for H = 0 the particle does not burn.
The Stokes number characterizes the velocity non-equilibrium of the phases. For Stk 0 the particle is completely entrained by the fluid, while for Stk its motion occurs independently of the fluid (for example, under the action of body forces). In the intermediate region, as the Stokes number increases, the velocity lag of the particle relative to the fluid increases.
At Stk 0 , all parameters, excluding the combustion parameter H, become negligible, and a equilibrium two-phase flow with burning particles is characterized only by the parameter H.

6.2. Domain of Applicability of the Solution

The ratio (16) is transformed to the following form
v g v p = C 1 exp F + Stk Fr v g t δ 2 .
The first term in the right side (19) determines the influence of the initial conditions, the second term is the mass forces, and the third term is the flow gradient. The criterion for the admissibility of the expansion of the particle velocity into a series in terms of the parameter is the condition for the smallness of the degree of velocity non-equilibrium of the phases
v g v p v g 1 .
In this case, the smallness of the Stokes number is not required. For example, for a gradientless flow, the obtained relations are exact for any Stokes number.
The equilibrium conditions correspond to the case C 1 0 . There is no influence of the initial conditions on the further motion of the particle, and the equilibrium initial velocity is found from the relation
v p 0 = v g + Fr v g t 0 Stk .
The intensity of attenuation of the influence of the initial state of the particle is determined by the Stokes number, and when a particle of variable mass moves, the influence of the initial conditions depends on its size.
Body forces affect the rate of non-equilibrium motion of phases through the Froude criterion. The influence of body forces increases with increasing Stokes number, decreases as δ decreases, and manifests itself to the greatest extent in low-velocity regions.
In a gradientless flow d v g / d t = 0 , and there is no influence of the gas-dynamic factor. For d v g / d t 0 the influence of the gas-dynamic factor increases as the Stokes number increases and decreases along with δ . The value of the derivative d v g / d t depends on the geometric shape of the channel and reaches its maximum value in narrow channels and in places where the flow turns.
Taking the injection speed as the characteristic speed, and the ratio of the free volume of the combustion chamber to the fuel combustion area ( u = v w , L = W / S ) as the characteristic size, the Stokes number is written as
Stk = ρ p d p 0 2 18 μ v w S W ,
where W and S are volume and cross-sectional area of combustion chamber. When the propellant burns out, the parameter S / W decreases, which leads to a decrease in the degree of velocity non-equilibrium of the phases due to the weakening of the influence of the gas-dynamic factor.

6.3. Solution for Strong Injection

In channels with strong injection ( v w ), the relations (16) and (18) are presented in the analytical form. The particle coordinates are found from the relations
x = x 0 ch 2 t / 2 , y = 4 π atan exp π 2 t .
The time in the equilibrium motion of the phases is
t = 2 π ln tan π 4 y = 1 π ln 1 cos ( π y / 2 ) 1 + cos ( π y / 2 ) .
The constants of integration are found from the initial conditions
C 1 = π 2 Stk x p 0 , C 2 = 1 ϕ Stk .
Up to the terms O ( Stk 2 ) , the following relations, which give the velocities and coordinates of the particle at an arbitrary moment of time, are valid
u p = u g + π 2 4 Stk x p 0 exp t / Stk π 2 4 Stk x p ; v p = v g + 1 ϕ exp t / Stk Stk π 4 sin π y p ; x p = x π 2 Stk x cos π 2 y ; y p = y ϕ Stk ( 1 ϕ ) Stk exp t / Stk + Stk sin π 2 y .
The velocity distributions of the dispersed phase in the channel are shown in Figure 8 and Figure 9 (results are normalized to the maximum velocity in the cross section, u m , and injection velocity, v w ). Lines 1 correspond to the carrier flow, and lines 2–5 correspond to particles of various sizes. The influence of the initial non-equilibrium flow leads to deformation of the radial velocity profile of the dispersed phase near the injection surface. In general, the distributions of the radial velocity differ relatively little from each other in a wide range of problem parameters.

7. Numerical Results

The results of numerical simulation are shown in Figure 10 and Figure 11 for particles of different sizes. The Froude number is 1, and the Sf parameter is varied by appropriately changing the Stokes number. For each Stokes number, two trajectories for z p 0 = ± 1 (lower and upper walls of the channel) are given, and the position of the singular point (if it exists) is indicated. The prime corresponds to the particles injected from the bottom of the channel. The results are obtained at different initial particle velocity.
It follows from the calculations that for particles injected from the upper part of the channel, a change in the initial velocity has practically no effect on the shape of the phase trajectories (Figure 11a). Particles injected from the lower half of the channel turn out to be more sensitive to changes in the initial conditions (Figure 11b).
For Sf = 0.1 , the singular point shifts slightly down along the z axis, and its type is a stable node. When Sf is increased to 0.5, the singular point type changes from node to focus. Accordingly, the type of phase trajectories also changes, which in the vicinity of the singular point acquire a damped oscillatory character. An increase in the initial particle velocity slightly affects the results.
For Sf = 1 , the flow structure in the channel becomes more complicated. Due to the fact that the singular point is close enough to the lower surface of the channel, and the amplitude of particle oscillations near this point is significant, conditions arise for them to fall out of the flow. Particles located at the bottom of the channel at the initial moment at v p 0 = 0 cannot leave the injection surface. However, an increase in the initial velocity module violates this dynamic equilibrium and leads to separation particles from the injection surface of the channel and their involvement in the flow, followed by asymptotic approximation to the corresponding singular point.
For Sf = 1.5 there are no singular points on the phase plane. Particles leaving the upper half of the channel eventually reach the lower surface, from which they cannot leave at zero initial velocity. An increase in the initial velocity of the particles leads to a temporary removal of particles from the lower part of the channel into the flow, followed by their fallout onto the same surface.
The trajectories of particles emerging from different points on the injection surface of the channel are shown in Figure 12 (for Fr = 1 ).
At Sf = 0.5 , a specific zone of counter-oscillatory motion of particles is formed in the vicinity of the singular point. As the Sf parameter increases to 1, this zone expands, deforms, and captures the lower part of the channel. The change in the structure of this flow zone is accompanied by the precipitation of some particles from the flow onto the channel surface.
Depending on the Stokes number, the singular point is either a node or a focus. Outside the gradient zone, the critical Stokes number, which characterizes the change in the type of singular point, gradually increases as the singular point moves away from the center of the channel.
For Sf = 1 , a second saddle-type singular point appears on the channel surface. The presence of two singular points, their relative proximity to each other and to the channel walls leads to the formation of a complex structure of a two-phase flow. New effects arise, in particular, the precipitation of particles from the flow onto the channel walls and the impossibility of their detachment from the injection surface. This is due to the fact that the oppositely directed drag and mass forces acting on the particles are approximately equal. As a consequence, under these conditions, the influence of secondary factors increases, for example, the initial velocity of particles on the injection surface.

8. Conclusions

A qualitative analysis of the trajectories of an individual particle in a channel with fluid injection from walls has been carried out. The computational domain is divided into a zone of gradient flow, in which the linearized equations of particle motion are applied, and a near-wall zone, in which an approximate solution of the problem is obtained. The results of the qualitative analysis are compared with numerical calculations, on the basis of which the singular points of the individual particle trajectory are found for various Stokes and Froude numbers.
The distributions of the dispersed phase parameters were obtained up to terms linear with respect to the Stokes number. Factors influencing the rate of non-equilibrium motion of phases are identified, and the domains of applicability of the obtained solution are investigated. The velocity and concentration distributions of particles in channels with strong injection are derived.

Funding

The research is partially funded by the Ministry of Science and Higher Education of the Russian Federation as part of World-class Research Center program: Advanced Digital Technologies (contract No. 075-15-2020-903 dated 16 November 2020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Emelyanov, V.N.; Teterina, I.V.; Volkov, K.N.; Garkushev, A.U. Pressure oscillations and instability of working processes in the combustion chambers of solid rocket motors. Acta Astronaut. 2017, 135, 161–171. [Google Scholar] [CrossRef] [Green Version]
  2. Emelyanov, V.N.; Teterina, I.V.; Volkov, K.N. Dynamics and combustion of single aluminium agglomerate in solid propellant environment. Acta Astronaut. 2020, 176, 682–694. [Google Scholar] [CrossRef]
  3. Volkov, K.N.; Emelyanov, V.N. Flows of Gas with Particles; Physical and Mathematical Publishing House: Moscow, Russia, 2008. [Google Scholar]
  4. Varaksin, A.Y. Turbulent Flows of Gas with Solid Particles; Physical and Mathematical Publishing House: Moscow, Russia, 2003. [Google Scholar]
  5. Osiptsov, A.N. Lagrangian modeling of dust admixture in gas flows. Astrophys. Space Sci. 2000, 274, 377–386. [Google Scholar] [CrossRef]
  6. Berman, A.S. Laminar flow in channels with porous walls. J. Appl. Phys. 1953, 24, 1232–1235. [Google Scholar] [CrossRef]
  7. Taylor, G.I. Fluid flow in regions bounded by porous surfaces. Proc. R. Soc. Ser. A 1956, 234, 456–475. [Google Scholar]
  8. Culick, F.E.C. Rotational axisymmetric mean flow and damping of acoustic waves in a solid propellant rocket. AIAA J. 1966, 4, 1462–1464. [Google Scholar] [CrossRef]
  9. Bundy, R.D.; Weissberg, H.L. Experimental study of fully developed laminar flow in a porous pipe with wall injection. Phys. Fluids 1970, 13, 2613–2615. [Google Scholar] [CrossRef]
  10. Olson, R.M.; Eckert, E.R.G. Experimental studies of turbulent flow in a porous circular tube with uniform fluid injection through the tube wall. J. Appl. Mech. 1966, 3, 7–17. [Google Scholar] [CrossRef]
  11. Yamada, K.; Ishikawa, N. Simulative study on the erosive burning of solid rocket motors pages. AIAA J. 1976, 14, 1170–1176. [Google Scholar] [CrossRef]
  12. Dunlap, R.; Blackner, A.M.; Waugh, R.C.; Brown, R.S.; Willoughby, P.G. Internal flow field studies in a simulated cylindrical port rocket chamber. J. Propuls. Power 1990, 6, 690–704. [Google Scholar] [CrossRef]
  13. Traineau, J.C.; Hervat, P.; Kuentzmann, P. Cold-flow simulation of a two-dimensional nozzleless solid-rocket motor. AIAA Pap. 1986, 86-1447. [Google Scholar]
  14. Majdalani, J.; Rienstra, S.W. Two forms of the rotational solution for wave propagation inside porous channels. AIAA Pap. 2022, 2002-2985. [Google Scholar]
  15. Apte, S.V.; Yang, V. Unsteady flow evolution in a porous chamber with surface mass injection, part 1: Free oscillations. AIAA J. 2001, 39, 1577–1586. [Google Scholar] [CrossRef]
  16. Apte, S.V.; Yang, V. Unsteady flow evolution in a porous chamber with surface mass injection, part 2: Acoustic excitations. AIAA J. 2002, 40, 244–253. [Google Scholar] [CrossRef] [Green Version]
  17. Boutros, Y.Z.; Abd-el-Malek, M.B.; Badran, N.A.; Hassan, H.S. Lie-group method for unsteady flows in a semi-infinite expanding or contracting pipe with injection or suction through a porous wall. J. Comput. Appl. Math. 2006, 197, 465–494. [Google Scholar] [CrossRef] [Green Version]
  18. Griffond, J.; Casalis, G.; Pineau, J.-P. Spatial instability of flow in semiinfinite cylinder with fluid injection through its porous walls. Eur. J. Mech. B/Fluids 2000, 19, 69–87. [Google Scholar] [CrossRef]
  19. Denisikhin, S.; Emelyanov, V.; Volkov, K. Fluid dynamics of thrust vectorable submerged nozzle. Fluids 2021, 6, 278. [Google Scholar] [CrossRef]
  20. Emelyanov, V.N.; Volkov, K.N.; Yakovchuk, M.S. Multiparameter optimization of thrust vector control with transverse injection of a supersonic underexpanded gas jet into a convergent-divergent nozzle. Energies 2021, 14, 4359. [Google Scholar] [CrossRef]
  21. Johnston, W.A.; Murdock, J.W.; Koshigoe, S.; Than, P.T. Slag accumulation in the titan SRMU. J. Propuls. Power 1995, 11, 1012–1020. [Google Scholar] [CrossRef]
  22. Salita, M. Defficiencies and requirements in modeling of slag generation in solid rocket motors. J. Propuls. Power 1995, 11, 10–23. [Google Scholar]
  23. Cesco, N.; Dumas, L.; Pevergne, T.; Fabignon, Y. Stochastic models to the investigation of slag accumulation in a large solid rocket motors. AIAA Pap. 1997, 97-2785. [Google Scholar]
  24. Sabnis, J.S. Numerical simulation of distributed combustion in solid rocket motors with metallized propellants. J. Propuls. Power 2003, 19, 48–55. [Google Scholar] [CrossRef]
  25. Dupays, J.; Prevost, M.; Tarrin, P.; Vuillot, F. Effects of particulate phase on vortex shedding driven oscillations in solid rocket motors. AIAA Pap. 1996, 96-3248. [Google Scholar]
  26. Vuillot, F.; Basset, T.; Dupays, J.; Daniel, E.; Lupoglazoff, N. 2D Navier–Stokes stability computations for solid rocket motors: Rotational, combustion, and two-phase flow effects. AIAA Pap. 1997, 97-3326. [Google Scholar]
  27. Cai, W.; Ma, F.; Yang, V. Two-phase vorticoacoustic flow interactions in solid-propellant rocket motors. J. Propuls. Power 2003, 19, 385–396. [Google Scholar] [CrossRef]
  28. Daniel, E.; Thevand, N. Stability of acoustic wave in two-phase dilute flow with mass transfer. AIAA J. 2001, 39, 2121–2130. [Google Scholar] [CrossRef]
  29. Majdalani, J.; van Moorhem, W.K. Laminar cold-flow model for the internal gas dynamics of a slab rocket motor. Aerosp. Sci. Technol. 2011, 5, 193–207. [Google Scholar] [CrossRef]
  30. Feraille, T.; Casalis, G. Channel flow induced by wall injection of fluid and particles. Phys. Fluids 2003, 15, 348. [Google Scholar] [CrossRef]
  31. Kabeel, A.; Elkelawy, M.; Bastawissi, H.; Elbanna, A. Solid particles injection in gas turbulent channel flow. Energy Power Eng. 2016, 8, 367–388. [Google Scholar] [CrossRef] [Green Version]
  32. Al-Malki, M.; Hussain, Z.; Garrett, S.; Calabretto, S. Effects of parietal suction and injection on the stability of the Blasius boundary-layer flow over a permeable, heated plate. Phys. Rev. Fluids 2021, 6, 113902. [Google Scholar] [CrossRef]
  33. Volkov, K.N. Motion of a condensed particle in a channel with injection. J. Eng. Phys. Thermophys. 2006, 79, 83–92. [Google Scholar] [CrossRef]
  34. Volkov, K.N. Application of the parameter-expansion method to calculation of two-phase flows in channels with injection. J. Eng. Phys. Thermophys. 2006, 79, 1169–1178. [Google Scholar] [CrossRef]
Figure 1. Flow induced by fluid injection.
Figure 1. Flow induced by fluid injection.
Fluids 07 00328 g001
Figure 2. Axial velocity profiles in a channel for x / h = 5 (a); 15 (b); 30 (c). Solid lines correspond to the results of numerical simulation [20], dashed lines correspond to vortical flow described by Equation (1), and symbols b u l l e t correspond to experimental data [8].
Figure 2. Axial velocity profiles in a channel for x / h = 5 (a); 15 (b); 30 (c). Solid lines correspond to the results of numerical simulation [20], dashed lines correspond to vortical flow described by Equation (1), and symbols b u l l e t correspond to experimental data [8].
Fluids 07 00328 g002
Figure 3. Profiles of axial velocity (a) and radial velocity (b) in a channel. Lines 1 and 2 correspond to the velocity profiles in channels with weak ( Re 0 ) and strong ( Re ) injection.
Figure 3. Profiles of axial velocity (a) and radial velocity (b) in a channel. Lines 1 and 2 correspond to the velocity profiles in channels with weak ( Re 0 ) and strong ( Re ) injection.
Fluids 07 00328 g003
Figure 4. Frame of reference.
Figure 4. Frame of reference.
Fluids 07 00328 g004
Figure 5. Particle trajectories in the gradient zone without (1) and with (2) of body forces.
Figure 5. Particle trajectories in the gradient zone without (1) and with (2) of body forces.
Fluids 07 00328 g005
Figure 6. Dependence of the radial coordinate of a singular point on the Stokes and Froude numbers.
Figure 6. Dependence of the radial coordinate of a singular point on the Stokes and Froude numbers.
Fluids 07 00328 g006
Figure 7. Change in the structure of two-phase flow in a channel with increasing Froude number.
Figure 7. Change in the structure of two-phase flow in a channel with increasing Froude number.
Fluids 07 00328 g007
Figure 8. Distributions of axial velocity of particle for Stk = 0.025 (2); 0.05 (3); 0.075 (4); 0.1 (5). Line 1 correspond to fluidflow.
Figure 8. Distributions of axial velocity of particle for Stk = 0.025 (2); 0.05 (3); 0.075 (4); 0.1 (5). Line 1 correspond to fluidflow.
Fluids 07 00328 g008
Figure 9. Distributions of radial velocity of particle for ϕ = 1 (a) and ϕ = 0 (b). Lines 2–5 correspond to different Stokes numbers Stk = 0.025 (2); 0.05 (3); 0.075 (4); 0.1 (5). Line 1 correspond to fluidflow.
Figure 9. Distributions of radial velocity of particle for ϕ = 1 (a) and ϕ = 0 (b). Lines 2–5 correspond to different Stokes numbers Stk = 0.025 (2); 0.05 (3); 0.075 (4); 0.1 (5). Line 1 correspond to fluidflow.
Fluids 07 00328 g009
Figure 10. Trajectories of particle for v p 0 = 0 (a) and v p 0 = 1 (b) with different Stokes numbers Stk = 0.1 (1); 0.5 (2); 1.0 (3); 1.5 (4).
Figure 10. Trajectories of particle for v p 0 = 0 (a) and v p 0 = 1 (b) with different Stokes numbers Stk = 0.1 (1); 0.5 (2); 1.0 (3); 1.5 (4).
Fluids 07 00328 g010
Figure 11. Solutions on phase plane ( v p , z p ) for v p 0 = 0 (a) and v p 0 = 1 (b) with different Stokes numbers Stk = 0.1 (1); 0.5 (2); 1.0 (3); 1.5 (4).
Figure 11. Solutions on phase plane ( v p , z p ) for v p 0 = 0 (a) and v p 0 = 1 (b) with different Stokes numbers Stk = 0.1 (1); 0.5 (2); 1.0 (3); 1.5 (4).
Fluids 07 00328 g011
Figure 12. Trajectories of particles in cross section of channel for Sf = 0.5 (a) and Sf = 1 (b).
Figure 12. Trajectories of particles in cross section of channel for Sf = 0.5 (a) and Sf = 1 (b).
Fluids 07 00328 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Volkov, K. Influence of Gravitational Force on Particle Motion in the Channel Flow Induced by Fluid Injection. Fluids 2022, 7, 328. https://doi.org/10.3390/fluids7100328

AMA Style

Volkov K. Influence of Gravitational Force on Particle Motion in the Channel Flow Induced by Fluid Injection. Fluids. 2022; 7(10):328. https://doi.org/10.3390/fluids7100328

Chicago/Turabian Style

Volkov, Konstantin. 2022. "Influence of Gravitational Force on Particle Motion in the Channel Flow Induced by Fluid Injection" Fluids 7, no. 10: 328. https://doi.org/10.3390/fluids7100328

APA Style

Volkov, K. (2022). Influence of Gravitational Force on Particle Motion in the Channel Flow Induced by Fluid Injection. Fluids, 7(10), 328. https://doi.org/10.3390/fluids7100328

Article Metrics

Back to TopTop