Advances in Mathematical Physics
Volume 2017 (2017), Article ID 7042686, 9 pages
https://doi.org/10.1155/2017/7042686

Research Article

Solutions of Navier-Stokes Equation with Coriolis Force

Sunggeun Lee,¹ Shin-Kun Ryi,² and Hankwon Lim¹

¹Department of Advanced Materials and Chemical Engineering, Catholic University of Daegu, Gyeongsan, Gyeongbuk 38430, Republic of Korea
²Korea Institute of Energy Research, Daejeon 305-343, Republic of Korea

Correspondence should be addressed to Shin-Kun Ryi; h2membrane@kier.re.kr and Hankwon Lim; hklim@cu.ac.kr

Received 15 April 2017; Accepted 3 July 2017; Published 14 August 2017

Academic Editor: Eugen Radu

Copyright © 2017 Sunggeun Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We investigate the Navier-Stokes equation in the presence of Coriolis force in this article. First, the vortex equation with the Coriolis effect is discussed. It turns out that the vorticity can be generated due to a rotation coming from the Coriolis effect, . In both steady state and two-dimensional flow, the vorticity vector gets shifted by the amount of . Second, we consider the specific expression of the velocity vector of the Navier-Stokes equation in two dimensions. For the two-dimensional potential flow , the equation satisfied by is independent of . The remaining Navier-Stokes equation reduces to the nonlinear partial differential equations with respect to the velocity and the corresponding exact solution is obtained. Finally, the steady convective diffusion equation is considered for the concentration and can be solved with the help of Navier-Stokes equation for two-dimensional potential flow. The convective diffusion equation can be solved in three dimensions with a simple choice of .

1. Introduction

Our Earth, even though it is hard to feel in daily life, is rotating. However, if we look around to see what is happening, we realize that this rotation leads to many interesting phenomena. These phenomena are unavoidable as long as we live on Earth. On a macroscopic scale, this effect is substantial: the climate is affected and complicated because the rotation changes the flows of the air heated by the sun. The air has a very complicated motion compared to the one without rotation. Typhoons and their motion are another phenomenon resulting from the effect of the rotation of Earth. Furthermore, their direction of rotation is opposite to each other in the northern and southern hemispheres. The flows of the ocean can be affected by the rotation of Earth as well. Changing the motion of moving objects or fluids such as air and ocean illustrated above is caused by a force. This force on all moving bodies is called Coriolis force. The rotating effect may be seen in outer space such as accretion disks of fast rotating black hole [1] as well. In this case, the relativity and the curved space-time should be counted to understand the effect better. It has been applied commercially in the areas of, for instance, Coriolis flow meter, and the solutions to the flow equations of Navier-Stokes equation are very useful. Hence, the Coriolis effect in the Navier-Stokes equation should be included in these cases. In addition, when the long cylinder carrying flows of water or oil is concerned, the Coriolis effect is not negligible and should be considered in order to construct the cylinder safely. In wide aspects, many chemical reactions of the flowing fluids or gases should be understood with the inclusion of this effect.

In summary, in fact, almost all of the fluids flowing on Earth can be affected, so knowing the exact motion of the fluids is important for applications and problem solving. The flow even involves heat flows of the fluids as well. Our work is just a starting point to find such solutions. As a basic exercise, we investigate the vortex equation and the Navier-Stokes equations under the Coriolis force. Specifically, we assume that the solutions are time independent. In [2], the Coriolis effect on the three-dimensional Navier-Stokes equation was considered. That was the case of compressible fluid while our concern here is the case of an incompressible one.

The present work has the following structure. In Section 2, we construct the vortex equation with a term of Coriolis force. The Coriolis effect can generate the vorticity. Especially for two-dimensional flow, the interesting point of view is the fact that the angular velocity coming from Coriolis force comes into play just as a shift of vorticity vector. In Section 3, we try to find the solution to the Navier-Stokes equations with Coriolis force as well. We apply the solutions to convective diffusion for two-dimensional potential flow. The effect of rotation does not appear on the velocity potential and the concentration is independent of its effect. In three dimensions, we take simplified choice of and solve the convective diffusion equation. Finally, Section 4 finishes with conclusions and discussions.

2. Vortex Equation for Two-Dimensional Flow

In this section, we construct the equation satisfied by the vortex [3] arising from the Navier-Stokes equation. When the fluid motion is viscous, the Euler equation becomes replaced by the Navier-Stokes equation. For incompressible flow, the Navier-Stokes equation in the presence of Coriolis force is given by

where is Reynolds number. The last term of the first equation comes from the Coriolis effect with the angular velocity . The vorticity vector is given by

Let us apply curl on both sides of the Navier-Stokes equation and use the vector identity

We will focus on the case of constant . Then, we see that the vorticity vector satisfies the following differential equations for incompressible fluid :

where we changed the notation of Reynolds number as .

Note that even if the vorticity vector in (4) is only time dependent, that is, without spatial dependence, the Coriolis effect gives rise to the vorticity. In this case, the vorticity is given by

Moreover, we see that, by solving the equation

the vorticity can be generated even for a constant (for independent of both time and space in (4)). This means that the phenomenon of vortex can arise everywhere on Earth even though the effect of Coriolis on Earth is too small to be seen.

If we define such that

and consider stationary solution such that , then we have the following equation:

Thus, the angular velocity coming from Coriolis effect only appears as a shift of vorticity vector. Furthermore, if we consider two-dimensional flows such that

where and are unit vectors on and , respectively, and the angular velocity is a two-dimensional one, then the above equations are reduced to the following:

Note that if has one component, say , affecting the vorticity vector, say , in two-dimensional flow (we express the velocity vector as ), hence we can discard vector notation and the equation becomes

Considering the constant pressure case, the Navier-Stokes equation becomes

Even though the above equation is reduced to a simplified form, it is still hard to be solved. By restricting two dimensions, we may simplify the equation and have a chance to get a solution. Alternatively, changing the above second-order differential equations into a first-order differential equation like BPS solution [4, 5] in field theory is an open and challenging problem.

Solving the Navier-Stokes equation directly is a straightforward way to get a vorticity though the exact solutions are quite restricted. In the next section, we try to solve the steady two-dimensional Navier-Stokes equation.

3. Solution of Navier-Stokes Equations and Its Applications

Let us find the solution of Navier-Stokes equations of which velocities have the following expression [6]:

where is the velocity potential for the flow in the plane and is out-of-plane velocity component. The vorticity vector field is given by

and from the expression in (13) we have the vorticity vector

The contours of are known as the vortex lines of the flow. The zero vorticity comes out if the third component of the velocity is equal to zero in which it turns out to be a two-dimensional flow.

In components, after inserting , , and into the Navier-Stokes equation and focusing on the time independent velocities such that , we have the following set of equations:

The incompressible conditions give the continuity equation:

For the expression of the velocity as in (13) in plane, the independence of the third direction gives the following relation:

Therefore, the continuity equation can be reexpressed as

Finally, we assume the pressure is also independent of the third direction:

Since the specific choice of the velocities reduces the three-dimensional problems to two-dimensional problems and working in the complex coordinates might have a number of advantages thanks to the numerous available techniques in complex analysis (e.g., conformal mapping), we now introduce the complex coordinates such that

So, by applying chain rule, we get the following relations:

Now, let us rewrite the Navier-Stokes equation in complex coordinates by replacing derivatives in real coordinates with complex coordinates:

The addition and subtraction of the first two equations in (23) give

where and and as already noticed , , and are assumed to be constants. The third equation becomes

It is interesting to note that the Coriolis effect does not appear in the equation which is related to two-dimensional velocity. Its effect appears on the equations for and the pressure.

Now suppose, for simplicity, that one of ’s is nonzero such that (in other words, ) and . In this case, we have

The solution would be

However, we try to make it simpler by taking . Then, the special solution of this equation becomes

So, we have the Navier-Stokes equation in complex coordinates under the Coriolis force with only one nonvanishing component of :

The first equation means that is a holomorphic function. In this case, can be written as . By choosing negative sign and putting this holomorphic function with into the above equations, we get

So, when only one of ’s, say , is nonvanishing, the Coriolis effect appears only in the pressure equation.

When , the solution in complex coordinates is given by [6]

When we consider two-dimensional flow where the velocities are given by ,  , and and , we get the result similar to the case obtained when only one component of is nonvanishing. In Figure 1, we plotted and velocity of (31) for a given value of the viscosity 100 when the Coriolis effect does not exist. We see that the velocity vanishes along the -axis.

(a)

(b)

Figure 1: Plots of profiles of absolute value of and velocity when ’s vanish.

Now, let us consider the case where all three components of ’s are nonvanishing. When all components of ’s are not zero, the Navier-Stokes equations become

To solve these equations, we first find the solution to the equation that comes from the incompressible condition of the fluid. As a simple solution, we put as follows:

Then, the above equations become

By differentiating with and on the first and the second equations above, in turn we get equations like

As a result of further manipulation, we get the following equations:

which turns out to be

The further application of partial derivatives gives the relation

The remaining equation

becomes

With the help of , we have the second-order differential equations. Let us now solve the second-order partial differential equation:

First, let us set . Then, we have

or

Introducing the integrating factor ,

we have the solution to the differential equation

Introducing , we get

By dividing by obtained above, we have the following :

Integrating further leads to the following:

where the pressure can be determined by solving (36).

In Figure 2, we plotted and velocity for a nonvanishing value of the rotation. Compared to Figure 1, we see that the zero value axis rotates as value increases.

(a)

(b)

(c)

Figure 2: Plots of profiles of absolute value of (a) and velocities when (b) and (c).

In summary, we have solved the Navier-Stokes equation with the assumption of in the presence of Coriolis force:

Now, let us consider the convective diffusion which is an application of the solution of Navier-Stokes equation. The equation of convective diffusion is given by

where is a concentration and is a diffusion coefficient. For two-dimensional flow where and assuming that is independent of direction, the convective diffusion equation becomes

which can be rewritten as

In this case, we assumed the potential flow in which the velocity is given by . Following the same procedure in [7], we get

where is introduced such that

satisfies Laplace’s equation

where we introduced and such that

Solving the Navier-Stokes equation with the convective diffusion simultaneously, we can obtain the exact form of the concentration in our specific two-dimensional potential flow. If we choose and , then we have

Here, and are given by

In three dimensions, we see that the convective diffusion equation becomes

If we assume that the three-dimensional velocity is represented by and the last term vanishes, we can put

In this case, we can solve the above equation by setting concentration as . The solution of can be obtained from the last equation and can be obtained as was done above. In order to know the concentration in this simplified example, we need to use solution (49). Here, we have

In Figure 3, we plotted the concentration for vanishing , , and , respectively. The concentration is high along one of the axes when the rotation is zero but reduces along that direction as the rotation is included. Our further study, however, shows that the concentration grows along the diagonal line as the absolute value of the rotation becomes large.

(a)  

(b)

(c)

Figure 3: Plots of profiles of absolute value of concentration when (a), (b), and (c).

Using the curl operator, we are able to read the equation of vorticity. It is interesting to note that the above results show that there is no dependence on . The consideration of the independence on the third () direction leads to the result of independence of that direction. The situation would be different if we take another expression of the velocity depending on the coordinates, say . Let us suppose time independent and vanishing energy dissipation for the heat transfer in fluids. Then, we have the following equation:

where is thermal conductivity. Following the same procedure as was done for convective diffusion above, we will get the same formula for heat.

The solutions of vortices have a lot of applications in many areas. For example, let us consider the mixing of chemical species in a reactor. The molecular mixing and chemical reaction in the vortex pair were studied in [8]. The vortices enhance the diffusion rate and the reaction rate by increasing the interfacial area. In this case, the concentration depends on the vortex circulation which is defined as and Reynolds number called vortex strength . From the constant shift of the vorticity vector under the constant Coriolis effect, we draw the result that the circulation reduces by the term coming from . In [9], it was studied that vorticity accelerates the diffusion. It is closely related to the enhancement by the turbulence. It is interesting to find the relation between Coriolis effect and turbulence. We can study further these vorticity effects around the world. Those are about air of the atmosphere and the water at the sea. There are lots of chemical species and their reactions every moment including pollution as a kind of liquids or gases. The presence and the production of the vortices may help clean the environment by accelerating reactions. The solutions of vorticity and the solutions to the Navier-Stokes equation will give a better understanding of our nature.

4. Conclusions and Discussions

The Navier-Stokes equations under the addition of the Coriolis force are discussed. The angular speed of the rotation is assumed to be a constant. We have seen that the vorticity can arise either for vortices depending on time only or for constant vortices as a result of Coriolis effect. It is interesting to note that, in two-dimensional flows, its effect on the vorticity vector results in a shift of the vorticity vector by . It would be interesting to see what happens for the pair of vorticities or a number of vorticities.

By taking the specific expression for the velocity , we can solve exactly the nonlinear Navier-Stokes equations in the presence of Coriolis effect. The specific expression of the velocity leads to the second-order partial differential equations of the velocity in two dimensions. With the independence of the velocity on the third direction, say , we can analyze the velocity in two dimensions with complex coordinates.

Even though the Coriolis effect might be negligible in small scale, we may find some crucial role of its effect in the macroscopic regime. Particularly, to describe the air in the atmosphere, we need to include pressure and couple to temperature. The application of complex analysis such as conformal mapping (or transformation) and Joukowski transformation may help study more mathematically (by these transformations, many complicated flow boundaries may be transformed into regular flow boundaries and it is possible to obtain solutions for the flow around ellipses) [6, 10] and find a more general solution of the partial differential equations. If we can find the transformation between nonrotating and rotating frame, we may apply it to generate all known solutions easily.

The convective diffusion equation is solved for two-dimensional potential flow. With the help of the solution of the Navier-Stokes equation, we get the velocity and find the concentration. The effect of the Coriolis rotation is absent in that case. Perhaps it can be deduced that this result is related to the chance of the rotation vector. The only possibility of the rotation vector is through the plane.

The application of the two-dimensional flows gives the exact motion of the air in our atmosphere and may predict the drifting particle motion that is harmful to creatures breathing air including humans. Moreover, the motion of the reactants inside a rotating reactor is also one of the applications yet with the coupling of temperature.

In fluid (gas or liquid) motions related to the application to chemical reactions, the vorticity would be important. The behavior of the gas might affect the reactions. In [8], the effect of vortex pair was discussed in molecular mixing. In [9, 11], the enhancement of reaction by vortex was discussed. In [12], CFD method was applied to study the vortex effect on dispersion rate, pressure drop, and turbulent kinetic energy of mixing particles. It would be interesting to extend our work to these studies as well. In addition, it is interesting to study the stabilities of our work and other various flows under the Coriolis effect (see [13] for Couette flow with the relation to Taylor-Görtler vortices; they discussed that the Taylor-Görtler vortices are induced by destabilizing effect of accelerating Couette flow).

The understanding of the vortex and the solutions to the Navier-Stokes equation in various situations is very important. For instance, high Reynolds number turbulence is related to the vortex [14]. We can obtain the vortex from the solutions of Navier-Stokes equation. The multivortices solution is an interesting application to the reactions extending the effect of vortex pair. It might be necessary to use computers for the multivortices. It is also interesting to find the solution of Navier-Stokes equation giving multivortices.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This project is supported by the “R&D Center for Reduction of Non-CO₂ Greenhouse Gases (RE2015001690003)” funded by the Korean Ministry of Environment (MOE) as “Global Top Environment R&D Program.”

References

A. R. Prasanna and B. Mukhopadhyay, “Effect of Coriolis Force on Accretion Flows around Rotating Compact Object,” International Journal of Modern Physics D, vol. 12, p. 157, 2003.

A. Korshunova, “Exact solutions to the compressible Navier-Stokes equations with the Coriolis and friction terms,” http://arxiv.org/abs/0807.4311.

P. G. Saffman, Vortex Dynamics, Cambridge University Press, Cambridge, UK, 1992.

E. B. Bogomolny, “Stability of classical solutions,” Soviet Journal of Nuclear Physics, vol. 24, pp. 449–870, 1976.

M. K. Prasad and C. M. Sommerfeld, “Exact Classical Solution for the 't Hooft Monopole and the Julia-Zee Dyon,” Physical Review Letters, vol. 35, p. 760, 1975.

M. Z. Bazant and H. K. Moffatt, “Exact solutions of the Navier–Stokes equations having steady vortex structures,” Journal of Fluid Mechanics, vol. 541, no. 55, 2005.

M. Sen and K. T. Yang, “Laplace's equation for convective scalar transport in potential flow,” Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, vol. 456, p. 3041, 2004.

B. M. Categen and J. P. Aguirre, “Analysis of molecular mixing and chemical reaction in a vortex pair,” Physics of Fluids, vol. 2, no. 12, p. 2211, 1990.

P. Flohr and J. C. Vassilicos, “Accelerated scalar dissipation in a vortex,” Journal of Fluid Mechanics, vol. 348, pp. 295–317, 1997.

M. Z. Bazant, “Conformal mapping of some non-harmonic functions in transport theory,” Proceedings of the Royal Society A, vol. 460, p. 1433, 2004.

J. P. Crimaldi, J. R. Cadwell, and J. B. Weiss, “Reaction enhancement of isolated scalars by vortex stirring,” Physics of Fluids, vol. 20, no. 7, Article ID 073605, 2008.

J. D. Chung, J. W. Kim, and Y. M. Park, “A study on vortex generators to improve the mixing rate in the dry sorbent injection process of the flue gas desulfurization system,” Korean Journal of Chemical Engineering, vol. 27, no. 1, p. 83, 2010.

Y. H. Kim and M. C. Kim, “Relative energy stability analysis on the onset of Taylor-Görtler vortices in impulsively accelerating Couette flow,” Korean Journal of Chemical Engineering, vol. 31, no. 12, p. 2145, 2014.

D. I. Pullin and P. G. Saffman, “Vortex dynamics in turbulence,” Annual Review of Fluid Mechanics, vol. 30, pp. 31–51, 1998.