Next Article in Journal
Laser Doppler Velocimetry Test of Flow Characteristics in Draft Tube of Model Pump Turbine
Next Article in Special Issue
Study on the Microscopic Mechanism of Grouting in Saturated Water-Bearing Sand Stratum Based on VOF-DEM Method
Previous Article in Journal
Research on Micro-Pore Structure and 3D Visual Characterization of Inter-Salt Shale Based on X-CT Imaging Digital Core Technology
Previous Article in Special Issue
Numerical and Experiment Investigation on Novel Guide Vane Structures of Turbo Air Classifier
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simple Particle Relaxation Modeling in One-Dimensional Oscillating Flows

1
Boysen-TU Dresden-Research Training Group, Technische Universität Dresden, Chemnitzer Straße 48b, 01062 Dresden, Germany
2
Chair of Energy Process Engineering, Institute of Process Engineering and Environmental Technology, Technische Universität Dresden, George-Bähr-Straße 3b, 01069 Dresden, Germany
*
Author to whom correspondence should be addressed.
Processes 2022, 10(7), 1322; https://doi.org/10.3390/pr10071322
Submission received: 1 June 2022 / Revised: 28 June 2022 / Accepted: 2 July 2022 / Published: 5 July 2022
(This article belongs to the Special Issue Multiphase Flows and Particle Technology)

Abstract

:
The relaxation of a rigid particle suspended in a one-dimensional oscillating flow is calculated according to different drag models and the results are compared. Conditions are derived under which relaxation can be neglected or drag models can be substituted by simpler ones. This investigation is conducted analytically and graphically via the plane defined by the Reynolds number and amplitude parameter. This work matches various, mostly analytic drag models together to consider simple particle relaxation with a few, broad range input parameters and cover large parts of the plane spanned by Reynolds number and amplitude parameter.

1. Introduction

A rigid, spherical particle suspended in a harmonically oscillating one-dimensional flow is considered. This configuration can be found, for example, with ultrasonic levitators [1,2] or pulsation reactors [3,4,5]. In case the particle is fixed in position (or the particle executes harmonic oscillations in a fluid at rest), the interaction between the fluid and the particle as well as the resulting flow state is defined by two dimensionless numbers: the particle Reynolds number R e and the amplitude parameter ϵ . The particle Reynolds number R e = U ρ f d / η expresses the ratio between inertial and viscous forces with the slip velocity amplitude U, the density of the fluid ρ f , the particle diameter d, and the dynamic viscosity of the fluid η . The amplitude parameter ϵ = A / d sets the particle displacement A in relation to the particle diameter, while the particle displacement A = U / ω itself is the ratio of slip velocity amplitude to the angular frequency ω of the oscillating fluid. Reynolds number and amplitude parameter span a plane in which all flow states can be pinpointed, as shown in Figure 1. Note that the origin of Figure 1 has the coordinates ϵ = 1 , R e = 1 . In case the particle is fixed in place, the slip velocity amplitude equals the far-field fluid velocity amplitude U f . Within this work, the far-field fluid velocity will be referred to as fluid velocity and can be described as u f ( t ) = U f cos ( ω t ) . In case the particle is not considered fixed in place, the slip velocity u ( t ) differs from the flow velocity due to inertia and drag forces acting on the particle (relaxation). The particle will perform oscillations with the same frequency as the fluid, which is also true for the slip velocity. If harmonic behavior is assumed, the only unknown quantities, which are necessary to describe the slip velocity, are the slip velocity amplitude and the phase shift ϕ with respect to the fluid velocity [6]. The basic geometrical relations of these central quantities are displayed in Figure 2. In Section 2.1, several drag models are presented with their respective ranges of validity to calculate relaxation in large parts of the ϵ - R e plane. Afterward, in Section 2.2, the slip velocity is calculated and discussed. Distinguishing between cases where the particle can be considered fixed and where this is not the case is discussed in Section 2.3. In the last section (Section 3), criteria are derived under which the simple Stokes model can be applied in the validity ranges of the other models.

2. Method

2.1. Drag Models

If all quantities ( ρ p , ρ f , U f , ω , η , and d ) but U are known and constant, U can be determined by equating the inertia of the particle (here treated as a mass point) with its drag when omitting all other forces. The slip velocity is not only central in working with the ϵ - R e plane (since it is part of ϵ and R e ), but it is also crucial for further considerations like flow patterns and heat and mass transfer, which are not subjects of this work. While the validity of the Navier-Stokes equations (NSE) is assumed, the complex process of solving them can be circumvented by applying various simplified drag models in their respected ranges of validity. Several drag models are necessary to cover large parts of the ϵ - R e plane, while the relations in Figure 2 holds true, independent of the applied drag model. It will now be investigated how these models behave and relate to each other. The expressions for the drag force of each drag model can be found in Table 1.
While the described ranges of validity are listed on the right side of Figure 1, they are plotted in the ϵ - R e plane on the left side. In case two or three models are valid in the same area, the simplest one has a solid color, while the others have a striped pattern.
The first step is to distinguish between a steady and an unsteady flow. The main criterion for differentiating is whether the boundary layer thickness δ = η ρ f / ω is larger or smaller than the particle itself. This phenomenon is expressed by the Womersley number W o 2 = d 2 / δ 2 = R e / ϵ [7]. The Womersley number is sometimes also referred to as Stokes number [8] or frequency parameter M 2 [9] in literature. While the Womersley number is one bisector of the ϵ - R e plane, the other one is the Streaming Reynolds number R e S = R e ϵ , which is central in the consideration of Steady Streaming [10]. If W o 2 1 , the flow can be considered steady [11]. In case the Reynolds number is also smaller than unity R e < 1 , the viscous forces dominate and creeping flow can be assumed. This assumption enables the application of the simple Stokes drag model [12], which considers pure drag only and leads to closed analytic solutions. If W o 2 1 , while the Reynolds number exceeds the validity limit of the Stokes model R e > 1 , the flow can still be considered steady, but a drag model with a larger range of validity regarding the Reynolds number has to be applied. In this case, the Schiller & Naumann model [13] (SN) was chosen, but for this case, many other drag models can be applied as well [14]. The SN model differs in nature from the Stokes model in two important aspects. On the one hand, it is an empirical model derived via the correlation of experimental data, rather than being derived from the NSE by omitting inertial forces. On the other hand, applying the SN drag model leads to a non-linear motion ODE, which can only be solved numerically. In the opposite case of W o 2 1 , the drag model by Landau & Lifshitz [11] (LL) can be applied, which linearizes the NSE around small ϵ < 1 to ensure no detachment of the flow from the particle. In case neither W o 2 1 nor W o 2 1 , but R e < 1 , the drag model by Basset [15] is valid, which is explicitly derived for unsteady flow conditions. It is an extension of the Stokes model that, besides pure drag, also considers added mass and the pressure gradient. Additionally, the “history” of the fluid is factored in by a theoretically derived term. In this work, only the stationary state is considered, and therefore the lower boundary of the history integral t 0 is set to . A discussion of the transient case can be found with Coimbra & Rangel [16]. An uncolored area in the ϵ - R e plane is recognizable where none of the simplified models are valid. In this area, the NSEs need to be solved numerically to calculate the slip velocity. The introduced drag models are well established in the literature and have been validated via experiments (and DNS) several times. The SN model was derived from experimental data in the first place. Therefore, in this work, the drag models themselves and their minimal ranges of validity are considered valid.
Since the particle oscillates with the same frequency as the fluid, the particle velocity and, in extension, the slip velocity are bound by this condition and change only slightly with respect to the applied drag model. This is in contrast to a steady flow, where the differences between models tend to keep increasing over time and often grow exponentially. This means that the simple Stokes drag model is adequate over a wider parameter range in an oscillating flow than its perceived limit of validity of small Reynolds numbers R e < 1 in a steady flow. This phenomenon encourages the substitution of more complex models of drag with Stokes. This simplification is addressed in Section 3.

2.2. Slip Velocity Amplitude

The slip velocity amplitude can be calculated with each respected drag model under the assumption of a harmonic oscillating flow that can be described with u f ( t ) = U f cos ( ω t ) , the results of which are displayed in Table 1. The derivation of the slip velocity amplitude with Stokes, SN, and the LL model is given in Appendix A, Appendix B and Appendix C, respectively, while the solution for the Basset model is taken from Hjelmfelt & Mockros [6].
Each solution in Table 1 contains only the Womersley number W o 2 and the density ratio γ (In the case of the SN model also the Reynolds number). The density ratio γ = ρ p / ρ f is the ratio between the particle density ρ p and the fluid density. This implies that, in cases where relaxation needs to be considered, the above introduced ϵ - R e plane can still be utilized with γ as a parameter. The normalized slip velocity amplitude is displayed in Figure 3 for γ = 25,000, corresponding to iron particles in hot air (left), and γ = 8, corresponding to iron particles in water at STP (right). First, neglecting the particle relaxation will be discussed, followed by the comparison of drag models, both by means of Figure 3.

2.3. Particle Relaxation

The ϵ - R e plane can be loosely classified into two areas with respect to the slip velocity amplitude. In one part, related to large W o 2 , U / U f 1 holds true. In this area, relaxation can be neglected. The size of this area depends only on γ , while its border can be represented by a value of the Womersley number W o 2 . This means each value of W o 2 can be related to a density ratio value set (and vice versa) for which relaxation can be neglected. The lowest value of this set is now defined as the critical density ratio γ c r i t ( W o 2 ) . Since relaxation is always present with this problem a cut-off criterion needs to be defined for neglecting particle relaxation. Depending on the specific applied case, the criterion can, for example, be set to U / U f > 0.9 , U / U f > 0.95 , or U / U f > 0.99 . The critical density ratio γ c r i t with respect to W o 2 is displayed in Figure 4 for the suggested cut-off criteria. If the density ratio of a specific problem is above the line for the respective Womersley number, relaxation can be neglected, while the choice of the cut-off criteria becomes more important for high Womersley numbers, as can be derived from Figure 4. Graphically, it is even suggested that this finding is also valid for the white area of the ϵ - R e plane, where the NSEs need to be solved numerically.
In the other part of the ϵ - R e plane, U / U f 1 and relaxation does need to be considered via adequate drag models.

3. Result and Discussion

For large γ , the quasi-steady approximation is valid for the entire relaxation part of the ϵ - R e plane, and the slip velocity amplitude can be calculated by Stokes or SN. This case is displayed on the left side of Figure 3, where γ is large enough for Stokes to reach U / U f 1, within its limit of validity R e < 1 . For small γ (right side of Figure 3) this is not the case. The same approach can be applied to the comparison between Stokes and SN. Here it is not γ , but rather the ratio between R e and S t k that is the determining factor. This relates to the distinction of U / U f 1 being reached before or after R e becomes unity. Even though the limit W o 2 < < 1 was set here to W o 2 = 0.01 , the slip velocity amplitude calculated with the Basset model still deviates substantially from the Stokes model. (Note here the logarithmic scale exaggerating the deviation between Basset and Stokes). Similar to the comparison between Stokes and Basset, on the right-hand side of Figure 3, the values at the model limits between Stokes and SN and even between Basset and LL do not match well, although this effect is not as pronounced. This is contrary to the high-density ratio case displayed on the left side of Figure 3, where the values at the model limits do match well. This raises the question of under which conditions the models connect to each other well or can even be applied interchangeably. The relation between Stokes and other models can be calculated with the expressions presented in Table 2, while harmonic particle velocities and slip velocities are assumed. The derivation of the deviation between Stokes and the SN model is given in detail in Appendix B. The expressions in Table 2 are derived from the expressions in Table 1 and also depend only on γ and W o 2 . This time, however, the oscillation Stokes number S t k = ω τ p = γ W o 2 / 18 is utilized [17]. It sets the characteristic time of the particle τ p = ρ p d 2 / 18 η in relation to the characteristic time of the flow τ f = 1 / ω . It works as a criterion for how well the particle can adapt to the changing flow conditions. When the normalized slip velocity amplitude in an oscillating flow is calculated with the Stokes model, it is only dependent on the oscillation Stokes number U / U f = 1 / 1 + 1 / S t k 2 . Therefore, the oscillation Stokes number is utilized in comparison of the Stokes model with the other models. The criterion for an acceptable deviation between the models was defined as 5%, and the parameter areas where the deviation is above and below this criterion are displayed in Figure 5 for the various models. The density ratio γ is the decisive parameter for the comparison between Stokes and Basset as well as Stokes and the LL model, while for the comparison between Stokes and the SN model it is the Reynolds number. Additionally, at least one criterion is presented in Table 2 for the deviation between Stokes and the other drag models to stay below 5%, while multiple criteria are set in relation to each other via logical operators. In this case, the Stokes model can be applied in the validity range of the respective model. These criteria are also plotted in Figure 5. Simple criteria were chosen for the sake of convenience, as they can not match the 5% error interface perfectly. This means that the deviation will stay below 5% if the criterion is met, but it can also stay below 5% if the criterion is violated, as can be seen in Figure 5. Stokes can be used in order to calculate the slip velocity amplitude in the entire colored area of the ϵ - R e plane when all criteria in Table 2 are met, hence ( γ > 1000 ) ( γ S t k 6 / 5 > 370 ) ( S t k < 3 R e ) .

4. Conclusions

The special case of a rigid particle in a one-dimensional oscillating flow was considered. The ϵ - R e plane was introduced, in which all flow states of this problem are defined. The conditions under which particle relaxation can be neglected were derived. Subsequently, various drag models from the literature were considered to model particle relaxation as simply as possible in large parts of the ϵ - R e plane. It was pointed out under which conditions the simple Stokes model can be applied in the validity range of various other drag models. This was conducted graphically and analytically.
It was shown that the Stokes drag model can be applied validly for the calculation of particle relaxation for large sets of parameter combinations. This can even be true for R e 1 , as long as other criteria, connected to high-density ratios between fluid and particle and/or high oscillating Stokes numbers, are met. These parameter combinations are found with many gas-solid multi-phase flows, which makes this work especially interesting for these kinds of systems.

Author Contributions

Conceptualization, methodology, writing—original draft preparation, writing—review and editing and visualization, S.H.; supervision, project administration and funding acquisition, S.U. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank the Boysen-TU Dresden-Research Training Group for the financial support that has made this publication possible. The Research Training Group is co-financed by TU Dresden and the Friedrich and Elisabeth Boysen Foundation. Grant number: BOY-135.

Data Availability Statement

Data are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Adisplacement amplitude
aparticle acceleration
C D drag coefficient
dparticle diameter
Fforce
R e Reynolds number
S t k oscillation Stokes number
ttime
uvelocity
Uvelocity amplitude
W o 2 Womersley number (frequency parameter)
γ density ratio
ϵ amplitude parameter
η dynamic viscosity
ρ density
τ relaxation time
ϕ phase shift
ω angular frequency
Abbreviations
L L Landau & Lifshitz
N S E Navier-Stokes equations
O D E ordinary differential equation
SStokes
S N Schiller & Naumann
S T P standard temperature and pressure
Indices
Ddrag
Iinertia
pparticle
ffluid
0initial state

Appendix A. Derivation of the Slip Velocity Amplitude Calculated with the Stokes Drag Model

Starting point is the ever-present force equilibrium on the particle between inertia and drag. The drag force F D , S predicted by the Stokes model, as listed in Table 1, is F D , S ( t ) = 3 π η d u ( t ) . A spherical particle is considered with an inertia of F I ( t ) = ρ p π 6 d 3 u ˙ p ( t ) . Combining this with an oscillating flow around a zero mean that can be described with u f ( t ) = U f cos ( ω t ) leads to the ordinary differential equation (ODE)
u ˙ p + 1 τ p u p = 1 τ p U f cos ( ω t ) ,
Here, u ( t ) = u f ( t ) u p ( t ) is the slip velocity, u ˙ p ( t ) the instantaneous acceleration, and τ p = the characteristic time of the particle. The solution to ODE (A1) is
u ( t ) = u f ( t ) u p ( t ) = u 0 + U f 1 + S t k 2 cos ( ϕ 0 ) + S t k sin ( ϕ 0 ) e t τ p transient + U f S t k 1 + S t k 2 S t k cos ( ω t + ϕ 0 ) sin ( ω t + ϕ 0 ) resident
Here, Stk is the Stokes number as introduced in Section 3. Neglecting the transient part of Equation (A2) ( t > 3 τ p ) and finding the extreme value delivers the slip velocity amplitude, which can be expressed normalized as
U U f = 1 1 + 1 S t k 2

Appendix B. Derivation of Slip Velocity Amplitude Calculated with SN and Deviation from the Stokes Model

Considering a spherical rigid particle and starting from the basic equilibrium between the drag force F D = C D π 2 ρ f d 2 u 2 and the inertial force F I = ρ p 4 π 3 d 3 a p for Stokes (index S) and Schiller & Naumann (index SN):
C D , S π 2 ρ f d 2 u S 2 ( t ) = ρ p 4 π 3 d 3 a p , S ( t )
C D , S N π 2 ρ f d 2 u S N 2 ( t ) = ρ p 4 π 3 d 3 a p , S N ( t )
Dividing (A4) by (A5) and rearranging leads to
u S N 2 ( t ) u S 2 ( t ) = C D , S ( t ) C D , S N ( t ) a p , S N ( t ) a p , S ( t )
For harmonic relations, the following statements with particle acceleration a p and particle velocity u p hold true:
a p ( t ) = u ˙ p ( t )
U p = ω A p
Here, U p and A p denote the particle velocity amplitude and the particle acceleration amplitude, respectively. u ˙ p is the particle velocity derivation in respect of time. Relation (A6) also holds true for the extreme values and inserting (A8) leads to
U S N 2 U S 2 = C ^ D , S C ^ D , S N U p , S N U p , S
C ^ D is the maximum value of the drag coefficient when the slip velocity amplitude U is reached with the maximum Reynolds number
R e = U ρ f d η
Using the relations U S = S t k U p , S = τ p A p , S gives the expression
U S N U S = S t k C ^ D , S C ^ D , S N U p , S N U S N
When the Pythagorean theorem
U f 2 = U 2 + U p 2
is applied, which can be derived from Figure 2, Equation (A11) can be further simplified:
U S N U S = S t k C ^ D , S C ^ D , S N U f U S N 2 1
In order to set the slip velocity amplitude in relation to the fluid velocity amplitude, each side of Equation (A13) can be multiplied with the right side of Equation (A3):
U S N U f = S t k 1 + 1 S t k 2 C ^ D , S C ^ D , S N U f U S N 2 1
Inserting the drag coefficients C ^ D , S = 24 / R e and C ^ D , S N = ( 24 / R e ) S N with the amplitude Reynolds number given in (A10), leads to:
U S N U f = 1 1 + S N 2 S t k 2
U S N U S = S t k 2 + 1 S t k 2 + S N 2
Since only now are the different drag coefficients C D , i inserted, Equation (A14) can also be used in order to calculate the deviation between the Stokes drag model and any other drag model, which follows Newton’s drag formulation and can be expressed by a drag coefficient.

Appendix C. Derivation and Solution of Particle Motion with the Landau & Lifshitz Model

The point of origin is the constant equilibrium of the drag force and the inertia of the particle:
F D = F I
3 π η d 1 + d 2 δ u + 3 4 π d 2 2 η ρ ω 1 + d 9 δ d u d t = ρ p π 6 d 3 u ˙ p
In combination with an oscillating flow around a zero mean that can be described with u f ( t ) = U f cos ( ω t ) , Equation (A17) leads to the ODE
u ˙ p + 18 η 1 + d 2 δ 9 d η ρ f 2 ω 1 + d 9 δ + ρ p d 2 u p = 18 η 1 + d 2 δ 9 d η ρ f 2 ω 1 + d 9 δ + ρ p d 2 U f cos ( ω t ) 9 d η ρ f 2 ω 1 + d 9 δ 9 d η ρ f 2 ω 1 + d 9 δ + ρ p d 2 ω U f sin ( ω t )
Solving ODE (A19) and recovering the slip velocity is achieved by applying the relation u ( t ) = u f ( t ) u p ( t ) , while the extreme value is discovered by detecting zeros in the first differential, which leads to the expression for the normalized slip velocity amplitude
U L L U f = 1 18 γ W o 2 + 9 γ W o 2 + 9 2 γ W o + 1 2 γ + 1 2

References

  1. Rothlisberger, M.; Schmidli, G.; Schuck, M.; Kolar, J.W. Multi-Frequency Acoustic Levitation and Trapping of Particles in All Degrees of Freedom. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2022, 69, 1572–1575. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, Y.; Wu, L.; Wang, Y. Study on Particle Manipulation in a Metal Internal Channel under Acoustic Levitation. Micromachines 2021, 13, 18. [Google Scholar] [CrossRef] [PubMed]
  3. Heidinger, S.; Spranger, F.; Dostál, J.; Zhang, C.; Klaus, C. Material Treatment in the Pulsation Reactor—From Flame Spray Pyrolysis to Industrial Scale. Sustainability 2022, 14, 3232. [Google Scholar] [CrossRef]
  4. Klaus, C.; Wegner, K.; Rammelt, T.; Ommer, M. New Challenges in Thermal Processing. Interceram Int. Ceram. Rev. 2021, 70, 22–25. [Google Scholar] [CrossRef]
  5. Hoffmann, C.; Ommer, M. Reactors for fluid-solid reactions: Pulsation reactors. In Handbuch Chemische Reaktoren; Reschetilowski, W., Ed.; Springer Reference Naturwissenschaften; Springer: Berlin/Heidelberg, Germany, 2019; pp. 1–19. [Google Scholar] [CrossRef]
  6. Hjelmfelt, A.T.; Mockros, L.F. Motion of discrete particles in a turbulent fluid. Appl. Sci. Res. 1966, 16, 149–161. [Google Scholar] [CrossRef]
  7. Womersley, J.R. Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known. J. Physiol. 1955, 127, 553–563. [Google Scholar] [CrossRef] [PubMed]
  8. Clift, R.; Grace, J.R.; Weber, M.E. Bubbles, Drops, and Particles; Dover Publications: Mineola, NY, USA, 2013. [Google Scholar]
  9. Riley, N. On a Sphere Oscillating in a Viscous Fluid. Q. J. Mech. Appl. Math. 1966, 19, 461–472. [Google Scholar] [CrossRef]
  10. Riley, N. Steady Streaming. Annu. Rev. Fluid Mech. 2001, 33, 43–65. [Google Scholar] [CrossRef]
  11. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics, 2nd ed.; Course of Theoretical Physics; Pergamon: Oxford, UK, 1987; Volume 6. [Google Scholar]
  12. Stokes, G.G. On the Effect of the Internal Friction of Fluids on the Motion of Pendulums. In Mathematical and Physical Papers; Cambridge University Press: Cambridge, UK, 2009; pp. 1–10. [Google Scholar] [CrossRef]
  13. Schiller, L.; Naumann, A. VDI Zeitung; VDI: Düsseldorf, Germany, 1935. [Google Scholar]
  14. Goossens, W.R. Review of the empirical correlations for the drag coefficient of rigid spheres. Powder Technol. 2019, 352, 350–359. [Google Scholar] [CrossRef]
  15. Basset, A.B. On the motion of a sphere in a viscous liquid. Philos. Trans. R. Soc. Lond. (A) 1888, 179, 43–63. [Google Scholar] [CrossRef] [Green Version]
  16. Coimbra, C.F.M.; Rangel, R.H. Spherical Particle Motion in Harmonic Stokes Flows. AIAA J. 2001, 39, 1673–1682. [Google Scholar] [CrossRef]
  17. Sazhin, S.; Shakked, T.; Sobolev, V.; Katoshevski, D. Particle grouping in oscillating flows. Eur. J. Mech. B/Fluids 2008, 27, 131–149. [Google Scholar] [CrossRef]
Figure 1. Drag models considered for various combinations of Reynolds number R e , amplitude parameter ϵ , and Womersley number W o 2 . (left) Various drag models displayed in the ϵ - R e plane. The striped patterns indicate the areas of validity, while the solid colors mark the preferable model in the area. Note that the origin has the coordinates ϵ = 1 , R e = 1 . (right) Various drag models with their limits of validity.
Figure 1. Drag models considered for various combinations of Reynolds number R e , amplitude parameter ϵ , and Womersley number W o 2 . (left) Various drag models displayed in the ϵ - R e plane. The striped patterns indicate the areas of validity, while the solid colors mark the preferable model in the area. Note that the origin has the coordinates ϵ = 1 , R e = 1 . (right) Various drag models with their limits of validity.
Processes 10 01322 g001
Figure 2. Vector representation of fluid velocity amplitude U f , particle velocity amplitude U p , and slip velocity amplitude U, including phase shift angles in reference to the fluid oscillation of the particle velocity ϕ p and slip velocity ϕ , respectively.
Figure 2. Vector representation of fluid velocity amplitude U f , particle velocity amplitude U p , and slip velocity amplitude U, including phase shift angles in reference to the fluid oscillation of the particle velocity ϕ p and slip velocity ϕ , respectively.
Processes 10 01322 g002
Figure 3. Normalized slip velocity amplitude, calculated with various drag models. (left) γ = 25,000 (corresponds to iron particles in hot air); (right) γ = 8 (corresponds to iron particles in water at STP).
Figure 3. Normalized slip velocity amplitude, calculated with various drag models. (left) γ = 25,000 (corresponds to iron particles in hot air); (right) γ = 8 (corresponds to iron particles in water at STP).
Processes 10 01322 g003
Figure 4. Critical density ratio γ c r i t depending on the Womersley number W o 2 for cut off criteria of U / U f > 90 % , U / U f > 95 % , and U / U f > 99 % . Above the respected line, the particle motion is not affected by the flow and the slip velocity can be set equal to the fluid velocity (no relaxation case).
Figure 4. Critical density ratio γ c r i t depending on the Womersley number W o 2 for cut off criteria of U / U f > 90 % , U / U f > 95 % , and U / U f > 99 % . Above the respected line, the particle motion is not affected by the flow and the slip velocity can be set equal to the fluid velocity (no relaxation case).
Processes 10 01322 g004
Figure 5. Areas where the deviation between Stokes and the respective drag model is above (blue—Schiller & Naumann, purple—Basset, green—Landau & Lifshitz) and below (red) 5%. In the red areas, the Stokes model can be applied in the validity range of the respective model.
Figure 5. Areas where the deviation between Stokes and the respective drag model is above (blue—Schiller & Naumann, purple—Basset, green—Landau & Lifshitz) and below (red) 5%. In the red areas, the Stokes model can be applied in the validity range of the respective model.
Processes 10 01322 g005
Table 1. Amplitude ratio of slip to fluid velocity U / U f for several drag models dependent on the Womersley number W o 2 and the density ratio γ .
Table 1. Amplitude ratio of slip to fluid velocity U / U f for several drag models dependent on the Womersley number W o 2 and the density ratio γ .
NameDrag Force F D Slip Velocity Amplitude Ratio U / U f =
Stokes 3 π η d u 1 + 18 W o 2 γ 2 1 / 2
Schiller & Naumann 3 π η d u S N 1 + 18 S N W o 2 γ 2 1 / 2 S N = 1 + 0.158 R e 2 / 3
Basset 3 π η d u π 6 d 3 p x + π 12 d 3 ρ f d u d t + 3 2 d 2 π ρ f η t 0 t d u / d t t t d t 1 ( 1 + f 1 ) 2 f 2 2 1 / 2 f 1 = 1 γ 0.5 + γ + 9 2 W o 18 W o 2 + 9 2 W o 2 + 0.5 + γ + 9 2 W o 2
f 2 = 1 γ 18 W o 2 + 9 2 W o 18 W o 2 + 9 2 W o 2 + 0.5 + γ + 9 2 W o 2
Landau & Lifshitz 3 π η d 1 + d 2 δ u + 3 4 π d 2 2 η ρ w 1 + d 9 δ d u d t 18 γ W o 2 + 9 γ W o 2 + 9 2 γ W o + 1 2 γ + 1 2 1 / 2
Table 2. Relations of the drag models listed in Table 1 to the Stokes model, expressed with the oscillation Stokes number S t k and the density ratio γ . Additionally, the limit is given up to which the deviation stays below 5%, in which case the Stokes model can be applied in the validity range of the respective model. Multiple criteria are set in relation to each other via the logic operator ∧—‘and’.
Table 2. Relations of the drag models listed in Table 1 to the Stokes model, expressed with the oscillation Stokes number S t k and the density ratio γ . Additionally, the limit is given up to which the deviation stays below 5%, in which case the Stokes model can be applied in the validity range of the respective model. Multiple criteria are set in relation to each other via the logic operator ∧—‘and’.
NameRelation to Stokes U i / U Stk | U i U Stk | / U Stk < 5 %
Schiller &
Naumann
S t k 2 + 1 S t k 2 + S N 2 1 / 2 S N = 1 + 0.158 R e 2 / 3 3 R e S t k < 1
Basset 1 + 1 S t k 2 1 ( 1 + f 1 ) 2 f 2 2 1 / 2 f 1 = 2 1 γ 1 + 2 γ + 3 γ S t k 2 γ S t k + 3 γ S t k 2 + 1 + 2 γ + 3 γ S t k 2
f 2 = 2 1 γ 2 γ S t k + 3 γ S t k 2 γ S t k + 3 γ S t k 2 + 1 + 2 γ + 3 γ S t k 2
( γ S t k 6 / 5 > 370 ) γ > 5
Landau &
Lifshitz
S t k 2 + 1 S t k 2 3 2 γ S t k + 1 S t k 2 + 3 2 γ S t k + 1 2 γ + 1 2 1 / 2 ( γ > 1000 )
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Heidinger, S.; Unz, S.; Beckmann, M. Simple Particle Relaxation Modeling in One-Dimensional Oscillating Flows. Processes 2022, 10, 1322. https://doi.org/10.3390/pr10071322

AMA Style

Heidinger S, Unz S, Beckmann M. Simple Particle Relaxation Modeling in One-Dimensional Oscillating Flows. Processes. 2022; 10(7):1322. https://doi.org/10.3390/pr10071322

Chicago/Turabian Style

Heidinger, Stefan, Simon Unz, and Michael Beckmann. 2022. "Simple Particle Relaxation Modeling in One-Dimensional Oscillating Flows" Processes 10, no. 7: 1322. https://doi.org/10.3390/pr10071322

APA Style

Heidinger, S., Unz, S., & Beckmann, M. (2022). Simple Particle Relaxation Modeling in One-Dimensional Oscillating Flows. Processes, 10(7), 1322. https://doi.org/10.3390/pr10071322

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop