CC BY
17The Structure of Induced Diffusion Flows on a Wedge with
Curved Edges
N.F. Dimitrieva1, Yu.D. Chashechkin2
lInstitute of Hydromechanics, National Academy of Sciences of Ukraine, Kiev, Ukraine
e-mail: [email protected] 2A.Yu. Ishlinsky Institute for Problems in Mechanics, Russian Academy of Sciences,
Moscow, Russia e-mail: [email protected]
The mathematical model and the numerical implementation method permitting to study simultaneously all the elements of the internally multi-scale stratified flows are developed. For the first time, the induced diffusion flows on a fixed symmetrical wedge with the curved edges are calculated and visualized using the supercomputer systems. Influence both of the edge effects and the wedge lateral edge curvature upon the flow structure is shown. In a stably stratified environment, the integral force inducing self-motion of a free wedge along the neutral buoyancy horizon towards the top is formed by pressure deficit in the adjacent jet streams.
Keywords: stratified fluid, induced diffusion flows, self-motion.
DOI: 10.22449/1573-160X-2016-3-70-78
© 2016, N.F. Dimitrieva, Yu.D. Chashechkin © 2016, Physical Oceanography
Introduction. Dissolved substances in the fluid and suspended particles under the influence of the mass forces (the gravitational or centrifugal ones) are unevenly distributed and form a stable stratification. Its density field is determined by the temperature, impurity concentration and pressure profiles in accordance with the type of the equation of state. Nonequilibrium medium with molecular flow of the stratifying components is resting only if density gradients are parallel to the gravitational vertical. The molecular flow interruption at the impermeable boundaries of arbitrary shape generates induced diffusion flows. Theoretical studies of these flows started in the early 40s of the last century [1], and a little later they began to be simulated in the laboratory [2]
Interest to the study of gradient flows is maintained due to their prevalence in the environment. Under normal conditions, the perturbations are concentrated in thin layers at the impermeable surfaces and reach storm values forming large temperature gradients in the atmosphere at steep range slopes [3] and nearby glaciers [4]. In recent years, special attention is paid to the study of the influence of induced diffusion flows on the dynamics and structure of processes in the reservoirs, lakes, seas [5] and oceans, taking into account the influence of the global rotation [6]. The complexity of the equation of state of the marine environment and the presence of an additional stratifying component with its own diffusion coefficient («differential diffusion») was considered in some studies [7].
At the end of the last century the models of the process of setting of induced diffusion flows at various obstacles (inclined and horizontal plate [8, 9], cylinder and sphere) started to be developed along with analytical studies of stationary currents [1, 2, 5]. The laboratory tests using highly sensitive shadow devices were also carried out. In the experiments except large eddies the rosettes of dissipative
70
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
non-stationary gravity waves were observed at the obstacle poles. Removing from the obstacles, these rosettes were transformed into extended layers [10].
The integral force action of convergent flows, absent in the symmetric obstacles (sphere, cylinder or horizontal plate), assumes a finite value on the inclined plate and other bodies, non-symmetrical with respect to the line of gravity. The resulting pressure gradients are large enough and capable of inducing self-motion of the neutral buoyancy free bodies («diffusion fish» in parlance of article [11]), which plays an important role in the dynamics of the hydrosphere. Calculation of steady flows on the flat fixed wedge is carried out in [12]. The study of formation mechanism of the propulsive forces on the bodies having various shapes is of practical interest.
In the present paper, for the first time the setting up dynamics and structure of the two-dimensional induced diffusion flows on concave or convex symmetrical wedge was calculated.
Problem formulation and its solving method. Since the velocities of the studied flows are small compared to the sound velocity, the mathematical model of the given physical processes is based on the system of differential balance equations of inhomogeneous multicomponent fluid mechanics in the Boussinesq approximation and neglection of the compressibility effects [13]. This system includes the continuity equation, the momentum balance (Navier - Stokes), the substances and closing algebraic equation of state p(S (z)):
P = Po (exp(- z/A)+s ), (1)
divv = 0, (2)
— + (vV)v =--VP + vAv - sg, (3)
8t po
ds ^ . vz
— + vVs = k< As +—, (4)
dt s A
where S = S0 (z) + s is the total salinity, including the salt compression ratio with initial stationary stratification and perturbation marked S(x, z, t), p0 is the density at zero level (neutral buoyancy level); v = (vx,vz) if the velocity of fluid; P is the pressure excluding the hydrostatic one; v and ks are the constant coefficients of kinematic viscosity and salt diffusion; t is the time; g is the acceleration of gravity; V and A are the operators of Hamilton and Laplace; A = p /1 Vp |, N = sjg/A , Tb = 2n/N are the scale, frequency and period of buoyancy; the z axis is directed vertically upward. Member with the vertical component of the flow velocity in the formula (4) represents the convective transport of salinity, violating the initial stratification.
At the initial time t = 0 in the resting continuously stratified fluid the impenetrable wedge of the selected shape is placed horizontally. On the surface of the body Z the adhesion condition for the substance velocity and impermeability is set, and at an infinite distance - the decay of all disturbances:
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
71
dS 1 dz ---+ A dn ds
dn i dn
= 0, (6)
i
where n is the outward normal to the obstacle.
A set of dimensional parameters (v = 10-6m2/c, ks = 1,41-10-9 m2/s, gz = 9,8 m/s2, N = 1 s )is presented in the problem. They form the characteristic time scales t = Tb, including the linear macroscales (of buoyancy A = 9,8 m, of the
obstacle size L = 0,1 m) and microscales (SN = vN = 10-3 m and Sn = 4Ks IN = 4 -10-5 m). The characteristic velocity scale is given by the ratio of the basic scales UvN =4vN = 10-3 m/s, UKNs = JK~N = 4 -10-5 m/s.
Macroscales A and L characterize the initial stratification and geometry of the current. Dissipation nature microscales SN and SKs determine the lateral dimensions of the fine structure components. A wide range of the values of length scales indicates the complexity of the stratified current internal structure.
Research of the linearized fundamental equations and the results of laboratory modeling showed the large-scale elements of currents (waves and eddies) to be characterized by regular perturbation components of the total solution. Fine structure elements of the flows, visualized in the whole range of parameters of the studied processes, are described by the large variety of singularly perturbed components. All the components coexist and actively interact with each other, forming a fine evolving structure [14].
Scale ratios set the traditional dimensionless complexes, such as the buoyancy scale C = A/L, the Reynolds number Re = |v|L jv , the Peclet number
Pe = |v|Ljks , the Schmidt number Sc = v/ks and the Froude number Fr = |v|/NL,
where |v| is the flow velocity. For the solution of the problem of induced diffusion
flow on a fixed obstacle it is taken to be equal to UKs.
In the nonlinear formulation the system (1) - (6) provides the simultaneous study of all elements of flows within the single description in natural physical variables without additional constants and ratios. Due to multi-scale phenomena and nonlinear equations, one of the main tools for analysis of such evolutionary processes is the numerical modeling.
Computational domain sampling was carried out on the basis of the SALOME open-source generic platform (www.salome-platform.org), including blockMesh, topoSet and refmeMesh generating utilities supplied by OpenFOAM open source software (www.openfoam.com), taking into consideration their individual scales. The computational domain is a rectangle with the horizontally arranged symmetrical wedge with L = 10 cm length and h = 2 cm base height. The sides of the body are arcs of circles and radii. The positions of their centers in a Cartesian coordinate system centered at the wedge top are shown in Tab. 1.
72
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
Coordinates of the center,
CM
(-7.6; 126.7) (17.6; 125.7) (9.2; 41.4) (5.8; 126.7)
Wedge Geometry
Curvature radius, cm
126.9 125.7 43.2 10.0
T a b l e 1
Deviation from the face of the chord AR, cm
-0.1 0.1 0.3 1.4
Taking into account the shapes of the body the hexahedral block-structured grids with registration of the lines on the block boundaries were reconstructed. For
resolution of the fine structure components of flows to a minimum microscale SNf several cells of the computational grid in domains with large gradients were nested. The fulfilled test calculations on a coarser grid showed the local non-physical oscillations of the solution near the obstacle. These oscillations lead to the accumulation of errors and stop of the counting during continuous computations.
Numerical solution was based finite volume method in the open source software OpenFOAM. To take into account the stratification and diffusion effects the standard numerical model icoFoam, realizing the non-stationary flows of the homogeneous fluid, was supplemented with new variables (p and s) and the appropriate equations (1) and (4), as well as with new additional parameters (N, A , K, g, etc. [15]).
The boundary condition for salinity perturbation (6) was implemented by the extended utility funkySetBoundaryField, which permits to set the analytical expressions for the physical variables. To conduct the further analysis applying such utilities, as vorticity, funkySetFields, probes, forces and forceCoeffs, the informative physical quantities (vorticity and its baroclinic generation, rate of dissipation of mechanical energy) lacking in the standard OpenFOAM software package, were calculated.
Results and discussion. The complex multi-scale structure of the induced diffusion flows is illustrated by the field of the horizontal component of salinity perturbation gradient (Fig. 1). Positive values of the visualized value are indicated in light gray, negative ones - in dark gray. There is the same range of value change between the isolines. The image structure is consistent with the data of shadow flow visualization in a laboratory pool [11].
Nearby the corner points of the wedge the additional fine structure elements -the dissipative gravitational waves - are formed. Their maxima of the longitudinal component of salinity gradient perturbation reached |ds/dc|max = 4 -10-2 (salinity
perturbations s had order of 10-5). With the removal of obstacles, the salinity perturbation gradient sharply decreases and reaches about 10-6 at 5 cm distance horizontally and 0.5 cm vertically from the corner point.
PHYSICAL OCEANOGRAPHY NO. 3 (2016) 73
Fig. 1. Field of the horizontal component of salinity perturbation gradient ds/dx on the wedge with the curved edges (L = 10 cm, h = 2 cm, Tb = 6.28 c, t = t/Tb = 16): a - AR = -0.1 cm; b - AR = 0.3 cm; c - AR = 1.4 cm
Field Structure of the horizontal component of salinity gradient perturbation essentially depends on the sign of the curvature of the lateral edges: the more oblique the tops at the base of the wedge, the more pronounced the visualized beams of alternating bands are (Fig. 1, a). At the same time the maximum values |3s/dx|max weakly depend on the curvature of the wedge faces.
For a convex wedge (Fig. 1, b) with the angle between the base and the side close to 90°, the beam of fine structure elements spreads.
For the wedge with maximum height shifted relatively to the base, the pattern of flows is reconstructed (Fig. 1, c) and becomes similar to the one observed in the cylinder [8].
Impenetrable obstacle immersed in the fixed stably stratified fluid forms a complex system of flows, including the main thin jets from the adjacent compensation counterflows along the sloping sides (Fig. 2). When the obstacle is removed, the width of layers increases, and the maximum velocity drops.
In the distribution of the horizontal velocity component vx the main jet thickness of the adjacent flow increases along with a decrease of the local tangential tilt angle to the neutral buoyancy horizon. Accordingly, for the concave wedge the thickness is reduced from 0.8 cm nearby the sharp top to 0.2 cm near the base (Fig. 2, a), and as for the convex wedge, it, on the contrary, increases (Fig. 2, b).
74
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
Fig. 2. Components of the flow velocity on the wedge (L = 10 cm, h = 2 cm, Tb = 6.28 s, t = T/Tb = 16): a, b) - longitudinal v ; c, d) - vertical v ; a, c) - AR = -0.1 cm; b, d) - AR = 0.3 cm
In the field of the vertical velocity component vz (Fig. 2, c, d) the main jets (layers having maximum velocity value) detach from the surface at the corner points at the wedge base and form a system of dissipative gravitational waves. There is an unperturbed area in the middle of the wedge base, as the vertical surface does not violate the diffusion impurity transfer.
Distributions of the flow velocity components near the wedge tops are shown in Fig. 3. Near the sharp top, the flow bifurcation point, there are precipitous changes in the vertical component of the velocity (curves 1 and 2 in Fig. 3, b). The values of the velocity components depend on the angle of wedge inclination edge to the x axis.
Around the corner points of the wedge base the complex alternating structure is formed. There the vertical component of the velocity vz is comparable to the absolute value of the horizontal one vx (curves 3 and 4 in Fig. 3, b).
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
75
Fig. 3. The profiles of velocity component at top (curves 1 and 2) and lateral base (curves 3 and 4) of the wedge (L = 10 cm, h = 2 cm, Tb = 6.28 s, t = T/Tb = 16): a) - longitudinal vx; 6) - vertical v z ; curves 1, 3 - Ar = -0.1 cm; 2, 4 - AR = 0.3 cm
In the sharp top of the wedge, and also in the small flow along its lateral sides the pressure perturbation decreases (Fig. 4). Its value depends on the curvature of the sides. Thus, the thickness of the pressure deficit extended area on the concave wedge under AR = -0.1 cm is about 3 cm (Fig. 4, a). It is significantly more than on the convex wedge: under AR = 0.3 cm the thickness of this area is 2.0 cm, under AR = 1.4 - 1.5 cm (Fig. 4, c). For a convex wedge in Fig. 4, there is a deficit of pressure in all the tops, due to change of the edge inclination angle.
The pressure difference creates a horizontal component of an integral force Fx , pushing the horizontal wedge towards the top, which is consistent with the observations in [11, 12]. The calculated values of Fx are shown in Tab. 2 when t = t/Tb = 16.
T a b l e 2
Dependence of the Horizontal Component of the Composite Force Acting on a Wedge and the Curvature of its Surfaces (t = t/Tb = 16)
AR, cm -0.1 0.0 0.1 0.3 1.4
Fx, H 8.4-10-6 8.3-10-6 8.0-10-6 -6.1-10-6 -2.2-10-6
76
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
The vertical component of the composite force Fz, calculated separately for the top and bottom edges, is in excess of Fx in two orders. Due to the uncollinearity, they compensate each other's action.
Conclusion. The mathematical model and the numerical implementation method permitting to study simultaneously all the elements of the internally multi-scale stratified currents without additional hypotheses and links are developed on the basis of the fundamental system of equations. Criteria for the reconstruction and optimization of the grids are formulated based on the large-scale analysis of systems of equations. For the first time, the induced diffusion flows on a fixed symmetrical wedge with the curved edges are calculated and visualized using the supercomputer systems.
Influence both of the edge effects and the wedge lateral edge curvature upon the current structure is shown. In a stably stratified environment, the integral force inducing self-motion of a free wedge along the neutral buoyancy horizon towards the top is formed by pressure deficit in the adjacent jet streams. In accordance with the results of the experiments, the body in a stably stratified medium moves toward its top. The pressure perturbation value depends on the curvature of the lateral sides of the wedge.
Thin fluid jet streams, formed along each side of the wedge, generate the internal waves in the areas of convergence with the corner points. These waves form thin layers.
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
77
The results of calculations of the induced diffusion flows can be used as initial conditions for solving problems of flow of the stratified marine environment past obstacles.
Acknowledgements. The computations were carried out in UniHUB web laboratory (www.unihub.ru) and visualized by supercomputer in CC MSU Lomonosov (www.parallel.ru). The research was conducted in A. Ishlinsky Institute for Problems in Mechanics RAS under the financial support of RFBR (Project No. 15-01-09235) and Department of Energy, Engineering, Mechanics and Control Processes, RAS (Program IV-4-12 Dynamics of formation and interaction of waves and vortices in continuous media).
REFERENCES
1. Prandtl', L., 1949, "Gidroaeromekhanika [Hydroaeromechanics]", Moscow, Izd-vo inostr. lit., 488 p. (in Russian).
2. Phillips, O.M., 1970, "On flows induced by diffusion in a stably stratified fluid", Deep-Sea Res, vol. 17, pp. 435-443.
3. Shapiro, A., Fedorovich, E., 2014, "A boundary-layer scaling for turbulent katabatic flow", Bound.-Lay. Meteorol., vol. 153, iss. 1, pp. 1-17.
4. Oerlemans, J., Van Pelt, W.J.J., 2015, "A model study of Abrahamsenbreen, a surging glacier in northern Spitsbergen", The Cryosph., vol. 9, no 2, pp. 767-779.
5. Zyryanov, V.N., Lapina, L.E., 2012, "Sklonovye techeniya v moryakh, ozerakh i vodokhranilishchakh, obuslovlennye diffuzionnymi effektami [Slope currents in the seas, lakes and reservoirs due to diffusion effects]", Vodnye resursy, vol. 39, no. 3, pp. 292-303 (in Russian).
6. Garrett, C., MacCready, P. & Rhines, P.B., 1993, "Boundary mixing and arrested Ekman layers: rotating, stratified flow near a sloping boundary", Ann. Rev. Fluid Mech., vol. 25, pp. 291-323.
7. Gargett, A.E., 2003, "Differential diffusion: an oceanographic primer", Progr. Oceanogr., vol. 56, pp. 559-570.
8. Linden, P.F., Weber, J.E., 1977, "The formation of layers in a double-diffusive system with a sloping boundary", J. Fluid Mech., vol. 81, pp. 757-773.
9. Chashechkin, Yu.D., Zagumennyy, Ya.V., 2012, "Techeniya nepreryvno stratifitsirovannoy zhidkosti, indutsirovannye preryvaniem diffuzionnogo perenosa nepodvizhnoy plastinoy [Continuously stratified fluid flows induced by disruption of the diffusion transfer by fixed plate]", Morskoy gidrofizicheskiy zhurnal, no. 5, pp. 3-22 (in Russian).
10. Zagumennyi, Ia.V., Chashechkin, Yu.D., 2013, "The structure of convective flows driven by density variations in a continuously stratified fluid", Physica Scripta, vol. 155, article id. 014034.
11. Mercier, M.J., Ardekani, F.M. & Allshouse, M.R. [et al.], 2014, "Self-propulsion of immersed object via natural convection", Phys. Rev. Lett., vol. 112, no. 20, p. 204501(5).
12. Page, M.A., 2010, "Propelled by diffusion", Nature Physics, vol. 6, pp. 486-487.
13. Landau, L.D., Lifshits, E.M., 1986, "Teoreticheskaya fizika. T. VI. Gidrodinamika [Theoretical physics. Vol. VI. Hydrodynamics]", Moscow, Nauka, 736 p. (in Russian).
14. Chashechkin, Yu.D., 2014, "Differentsial'naya mekhanika zhidkostey: soglasovannye analiticheskie, chislennye i laboratornye modeli stratifitsirovannykh techeniy [Differential fluid mechanics: coordinated analytical, numerical and laboratory model of stratified currents]", Vestnik MGTU im. N.E. Baumana. Ser. Estestvennye nauki, no. 6, pp. 67-95 (in Russian).
15. Dimitrieva, N.F., Zagumennyy, Ya.V., 2014, "Chislennoe modelirovanie stratifitsirovannykh techeniy s ispol'zovaniem OpenFOAM [Numerical modeling of stratified flows using OpenFOAM]", Tr. ISP RAN, vol. 26, no. 5, pp. 187-200 (in Russian).
78
PHYSICAL OCEANOGRAPHY NO. 3 (2016)
