Abstract
In this paper we consider a hydrodynamic model for the matter density distribution in a self gravitating, isentropic 2-d disk of gas where the isentropy coefficient is allowed to be a function of position. For this model we prove analytically the existence of steady state and time dependent solutions in which the matter density in the disk is oscillatory and pattern forming. This research is motivated in part by recent astronomical observations and Laplace conjecture (made in 1796) that planetary systems evolve from a family of isolated rings that are formed within a primitive interstellar gas cloud.
Similar content being viewed by others
1 Introduction
In the last few decades many researcher explored different kinds of systems which self-organize and lead to pattern formation. Physical systems which display such phenomena include magnetization, Benard convective cells, cloud formations in stripes or rolls and many others. The challenge in all these instances is to find an appropriate mathematical model for the emergence of patterns in these self-organizing systems (Chandrasekhar 1946; Kiguchi et al. 1987; Lissauer 1993; Matsumoto and Hanawa 1999; Petigura et al. 2013; Spitzer 1968; Von Weizsacker 1944; Woolfson 2000).
One of the first conjecture about pattern formation under the action of gravity was formulated by Laplace (2010). He put forward the hypothesis that a primitive interstellar gas cloud evolve under gravity to form a system of isolated rings which may in turn lead to the formation of planetary systems. More recently Belrage (1968) and Prentice (1978) followed on Laplace hypothesis and explored possible mechanisms for the formation of ring structure in the solar nebula.
The issue of ring formation has become especially relevant and timely since the discovery of a remarkable pattern of bright and dark rings in the circumstellar disc surrounding the young solar analogue star HL Tau (ALMA Partnership 2015). It is also a relevant issue due to the large number of exoplanetary systems that have now been discovered. The orbital architecture of these systems is consistent with their condensation from dense rings of gas in proto-stellar disks, just as is believed to be the case for our solar system (Prentice 1978). The new astronomical data leads us to believe that there is a fundamental physical process at work, not yet fully understood, in which ring pattern formation is crucial to our understanding of the origin of planetary systems.
Many theories were put forward in the past about the origin of solar system. Currently the leading theory about the formation of planetary systems is the “Nebula Theory” whereby a cloud of interstellar gas accreted under it own gravitation, to form in stages, the protostar and the planets. Many of the results related to this theory were obtained through elaborate modeling and large scale numerical simulations. These involve, in general, thermodynamic considerations, magnetohydrodynamics modeling and turbulence (Ya Marov et al. 2013).
In a series of paper we addressed the issue of pattern formation under gravity from idealized analytical point of view (Humi 2006, 2009, 2016). In particular we showed that a self-gravitating incompressible rotating 2-dimensional gas admits steady states with ring structure in the matter distribution. It is our objective in this paper to examine this pattern formation problem under the assumption that the gas is isentropic. This is motivated by the desire to take into account (at least partially) the irreversible thermodynamic processes taking place in the disk which lead to the emission of radiation and heat. Under this idealized scenario the entropy produced within the disk due to the irreversible thermodynamic and turbulent processes taking place is removed by heat and radiation and the gas maintains a constant entropy. Obviously this assumption implies that this is a prototype “barebone” model that neglects many other important aspects of the processes taking place in an actual primitive cloud of interstellar gas.
The plan of the paper is as follows: In Sect. 2 we present the basic model equations. In Sect. 3 we consider the radial steady state of these equations and demonstrate that they admit solutions with patterns of density oscillations. In Sect. 4 we discuss radial time dependent states and establish again the existence of pattern forming solutions. We end up in Sect. 5 with some conclusions and observations.
2 Basic Equations
The basic equations that govern the time dependent state of non-relativistic self gravitating, compressible, inviscid two dimensional gas (in the \(x-y\) plane) are Prentice (1978), Humi (2006), Humi (2009) and Li et al. (2006)
where subscripts indicate differentiation with respect to the indicated variable, \(\mathbf{u}=(u,v)\) is the gas velocity, \(\rho\) is its density, p is the pressure, \(\phi\) is the gravitational potential, G is the gravitational constant, t is time and \(\omega\) represents a solid body rotation of the gas.
In this set of Eq. (2.1) is the continuity equation, and (2.2–2.3) are the momentum equations. Finally (2.4) is the equation for the gravitational potential.
We can nondimensionalize these equations by introducing the following scalings
where \(L,U_0,\rho _0\) are some characteristic length, velocity and mass density respectively that characterize the gas under consideration. Substituting these scalings in Eqs. (2.1–2.4) and dropping the tildes these equations remain unchanged (but the quantities that appear in these equations become non-dimensional).
3 Radial Steady States
In this section we consider radial steady states of the isentropic gas. In polar coordinates \((r,\theta )\) when all variables depend only on r. Eqs. (2.1–2.4) become;
Solving (3.2) algebraically for \(\frac{d\rho }{dr}\) and substituting in (3.1) it follows that
and therefore
Similarly one can solve (3.1) for \(\rho u\) and obtain that
(where \(c_1\), \(c_2\) are integration constants).
The classical relationship between pressure and density for isentropic gas is
where A is a constant which can be normalized to one and \(\alpha\) is the isentropy index. A more general functional relationship between the pressure and density will allow the isentropy coefficient A to be a function of r. That is
This generalization is introduced to account for the possibility that the thermodynamic processes taking place within the gas are position dependent. An additional generalization can be obtained if we let \(\alpha\) to be a function of r. However in this paper we shall consider only the case where \(\alpha\) is a constant.
Substituting (3.6, 3.7 and 3.9) in (3.3) and assuming that \(\alpha\) is a constant we obtain (after differentiation to eliminate \(\phi\)) the following equation for \(\rho\)
In the following we present analytic solutions of (3.10) with density pattern.
3.1 Steady State Analytic Solutions with Density Patterns
When \(\alpha =2\) and A is a constant (3.10) becomes
When \(c_1=c_2=0\) (i.e. \(u=v=0\)) this equation has an exact analytic solution,
where \(w=\sqrt{\frac{2\pi G}{A}}\) and \(J_0\), \(Y_0\) are zero order Bessel functions of the first and second kind. Since \(Y_0\) is singular at the origin we let \(C_2=0\). Assuming that the disk has a “central body” of radius 1 and \(\rho (1)=0\) the solution becomes
For \(\omega =0.36\), \(A=0.02\) and \(G=0.05\) this solution is positive and have oscillatory pattern beyond the boundary of the central body viz. \(r >1\) (See Fig. 1).
Analytic solution of (3.11) in terms of hypergeometric functions can be obtained also when \(c_1=0\) and \(c_2 \ne 0\). To show the impact of this change on the solution we display the difference between this solution with \(c_2=0.02\) and the one above in Fig. 2 (Same values for A, G and \(\omega\) and same strategy to determine the integration constants). We see that the major difference between the two solutions happens in the region close to \(r=0\).
Another analytic oscillatory solution is obtained for \(A(r)=\beta r\), \(c_1=c_2=0\) where \(\beta\) is a constant. The general solution in this case can be written as
To avoid the singularity in this expression at \(r=0\) we let \(C_1=0\) and \(C_2=\frac{\beta \omega ^2}{8}\). The solution is then
For \(G=1\), \(\omega =1\), \(\beta =0.1\) this solution displays some density fluctuations whose amplitude decays as r increases.
One can obtain similar (but more complicated) analytic expressions for \(\rho\) when \(c_2 \ne 0\).
These solutions demonstrate that steady state density patterns do arise in isontropic gas under gravity.
For \(\alpha =1\) and \(\omega =c_1=c_2=0\) we obtain the following analytic solution for \(\rho\)
where \(d=\sqrt{\frac{C_1}{A}}\) and \(b=C_2 d\) and \(C_1,\,C_2\) are integration constants.
4 Time Dependent Radial Solutions
When all the variables in (2.1–2.4) depend only on r and t we obtain the following system of equations:
To simplify (somewhat) these equations we solve (algebraically) (4.2) for \(\frac{\partial \rho }{\partial t}\) and substitute in (4.1) to obtain
Similarly solving (4.3) algebraically for \(\frac{\partial \rho }{\partial t}\) and substituting in (4.1) yields,
4.1 Further Reductions for \(\boldsymbol{\alpha} \boldsymbol{=}{\bf 2}\)
Assuming that v(r, t) is independent of t it follows from (4.5) that
where \(c_1\) is a constant. Equation (4.4) can be used to express \(\phi\) in terms of \(\rho\).
where \(f_1(t)\) is some function of t.
Substituting (4.7) in (4.6) and using the fact that \(\alpha =2\) (viz. \(p=A\rho ^2\)) leads to
Solving this equation for \(\rho (r,t)+\phi (r,t)\) yields,
where \(f_2(t)\) is arbitrary. Using (4.8) to substitute for \(\phi (r,t)\), differentiating the result with respect to r, and then multiplying by r and differentiating again we obtain the following differential equation for \(\rho\)
where
The system (4.1), (4.11) for \(\rho\) and u can solved numerically. However it can be decoupled also by solving (4.11) for \(\rho\)
where \(a=\sqrt{\frac{2\pi G}{A}}\) and \(f_3,\,f_4\) are arbitrary. However since \(Y_0\) is singular at the origin we let \(f_4=0\). Substituting this result in (4.1) yields a (complicated) intgro-differential equation for u(r, t) . However the key point here is that \(\rho (r,t)\) in (4.13) is expressed in terms of Bessel functions whose values are modulated by H(r, t). Since Bessel functions are oscillatory it is plausible to assume that for proper functional values of u(r, t) this may lead to time dependent pattern formations in the density.
4.2 Time Dependent Perturbations to a Steady State
The system (4.1–4.4) admits a steady state \(\rho =\rho _0(r)\) and \(u=v=0\). We consider here a perturbation from this steady state with \(\alpha =2\) in the form
Since \(\alpha =2\) the relevant equations to satisfy are (4.1) and (4.11). Substituting the expressions in (4.14) in these equations we obtain to first order in \(\epsilon\) the following equations
Introducing the unsatz
we find that the time dependence in (4.15), (4.16) factors out and we obtain the following two equations for \(g_1\) and \(g_2\)
It is possible to reduce the system (4.18), (4.19) to a single equation by solving (4.18) algebraically for \(g_1\) and substituting the result in (4.19). For the steady state \(\rho _0=\frac{1}{2}\frac{\omega ^2}{G\pi }\) this procedure leads to the equation for \(g_2\),
The solutions for \(g_2\) and \(g_1\) (which have no singularity at \(r=0\)) are
Where \(J_0\) and \(J_1\) are Bessel functions of the first kind and
This solution demonstrates that perturbations from the steady state \(\rho =\frac{1}{2}\frac{\omega ^2}{G\pi }\) yield patterns of density oscillations within the cloud.
4.3 Explicit Solutions
To provide some explicit time dependent pattern forming solutions we consider the system (4.1), (4.5), (4.6) and (4.4) under the assumption that \(p(r,t)=A\rho (r,t)^{\alpha (r)}\), viz. we let the isentropy index to vary with r. To derive these solutions we use separation of variables on this system of equations. That is we look for solutions in the form
From (4.1) we have
(Here and in the following \(\mu _i\) are separation of variables constants). Similarly (4.5) leads to
Equation (4.4) yields that \(g_3(t)=g_4(t)\) (up to a multiplicative constant which can be normalized to 1) and \(f_4(r),\,f_3(r)\) must satisfy
Finally substituting the ansatz (4.22) in (4.6) leads to
To perform separation of variables on this equation it is obvious that additional assumptions are needed. We explore here two of these possibilities.
-
1.
Assume that \(g_5=g_3g_1^2\), \(g_3=g_1^2\), \(\omega =0\) and \(g_1=g_2\). (Here and in the following we ignore possible multiplicative constants which can be normalized to 1).
Using these assumptions we can perform separation of variables on (4.28) and obtain;
To solve for \(g_1,\,g_2\) and \(g_3\) we use (4.23) and (4.25). Since we assumed that \(g_1=g_2\) we can solve (4.25) to obtain
Substituting this result in (4.23) leads to
where \(C_1\), \(C_2\) are an integration constants. However since we assumed that \(g_3=g_1^2\) we must have \(C_2=1\) and \(\mu _1=2\mu _2\) These results will be consistent with (4.29) if we let \(\mu _2=\mu _3\).
This completes the determination of all the functions \(g_i(t)\,\,i=1,\ldots ,5\) (since \(g_4=g_3\) and \(g_5=g_3g_1^2\)).
Equations (4.24), (4.26), (4.28) and (4.30) form a redundant system of four equations for five functions \(f_i(r),\,i=1,\ldots ,5\). Therefore we can choose one of these functions as a “parameter function”. For example if we let
then from (4.24) and (4.26) we have
where \(C_4\), \(C_5\) are constants. The expressions for \(g_3\) and \(f_3\) imply that \(\rho\) is oscillatory in r but the amplitude of these oscillations decays in time.
Substituting \(f_3\) in (4.27) we can solve for \(f_4\). Using all these results (4.29) becomes a first order differential equation for \(f_5\) which can be solved numerically. A plot of the isentropic index (at \(t=0\)) in this case with \(a=2\), \(k=5\pi\), \(C_4=C_5=1\), \(\mu _2=\mu _3=-1\) \(\mu _1=-2\) is shown in Fig. 3. In this figure the isentropic index has an approximate value of 2 at \(r=0\) but then rises rapidly to a constant value of 2.6.
-
2.
Assume that \(g_5=g_3g_1^2\), \(g_3=g_1^2\), \(\omega =0\) and \(f_2^2=rf_1\).
Under present assumptions, separation of variables for (4.28) yields,
To determine the functions \(g_i(t)\) \(i=1...5\) in this case we use first (4.23) and (4.25) to obtain
Hence
However since we assumed that \(g_3=g_1^2\) it follows that
Substituting this result in (4.25) and solving for \(g_2\) yields
where \(C_1,\, C_2\) are integration constants. Therefore
These results will be consistent with (4.33) if we impose the constraints,
To solve for the functions \(f_i(r)\), \(i=1...5\) we use the assumption \(f_2^2=rf_1\) in (4.26) and solve for \(f_1\). This yields
where \(E_1\) is a constant. However to avoid the singularity at the origin we let \(E_1=0\) and hence \(f_1=-\frac{\mu _2}{2}r\). Substituting this in (4.24) and solving for \(f_3\) we find that
where \(E_2\) is a constant. However \(\mu _1=2\mu _2\), hence,
From (4.27) we obtain then (after applying the boundary conditions at \(r=0\))
Similarly (4.34) yields a polynomial expression for \(f_5\).
5 Conclusions
In this paper we explored the solutions of a simple two dimensional hydrodynamic model for a self gravitating and rotating isentropic gas cloud. Our main objective was to establish the existence of solutions with density oscillations (viz. pattern formation).
It might be argued that this paper does not present a general theory for pattern formation under gravity. However it should be pointed out that the existence and classification of periodic solutions for a system of nonlinear differential equations is an open problem (Wayne 1997). Nevertheless our results establish analytically the existence of such solutions for the steady and time dependent states of a rotating, isentropic self-gravitating gas in two dimensions.
For the radial steady state case we were able to reduce the original set of four partial differential equations to one differential equation for the density of the gas. The steady state solutions of this equation were explored with constant isentropic index and isentropic coefficient that is a function of r.
For the time dependent case we showed, perturbations from a steady state lead to oscillation in the the density. Furthermore using separation of variables, we showed explicitly that analytic oscillatory closed form solutions of the equations exist.
References
H.P. Belrage, The Origin of the Solar System (Pergamon Press, Oxford, 1968)
S. Chandrasekhar, On a new theory of Weizsacker on the origin of the solar system. Rev. Mod. Phys. 18, 94–102 (1946)
M. Humi, Steady states of self gravitating incompressible fluid. J. Math. Phys. 47, 093101 (2006)
M. Humi, Steady states of self gravitating incompressible fluid with axial symmetry. Int. J. Mod. Phys. A 24(23), 4287–4303 (2009)
M. Humi, A model for pattern formation under gravity. Appl. Math. Model. 40(2016), 41–49 (2016)
M. Kiguchi, S. Narita, S.M. Miyama, C. Hayashi, The equilibria of rotating isothermal clouds. Astrophys. J. 317, 830–845 (1987)
P.S. Laplace, Exposition du Systeme du Monde (Reprinted by Cambridge University Press, Paris, 2010). ISBN 0511693338
J. Li, T. Zhang, Y. Zheng, Simple waves and a characteristic decomposition of the two dimensional compressible Euler equations. Commun. Math. Phys. 267, 1–12 (2006)
J.J. Lissauer, Planet formation. Annu. Rev. Astron. Astrophys. 31, 120–174 (1993)
T. Matsumoto, T. Hanawa, Bar and disk formation in gravitationally collapsing clouds. Astrophys. J. 521(2), 659–670 (1999)
A.L.M.A. Partnership, Astrophys. J. Lett. 808, L3 (2015). https://doi.org/10.1088/2041-8205/808/1/L3
E.A. Petigura, A.W. Howard, G.W. Marcy, Prevalence of Earth-size planets orbiting Sun-like stars. Proc. Nat. Acad. Sci. 110(48), 19273 (2013)
A.J.R. Prentice, Origin of the solar system. Earth Moon Planet 19, 341–398 (1978)
L. Spitzer Jr., Diffuse Matter in Space (Interscience Publishers, New-York, 1968)
C.F. Von Weizsacker, Z. Astrophys. 22, 319–355 (1944)
C.E. Wayne, Periodic solutions of nonlinear partial differential equations. Not. AMS 44, 895–902 (1997)
M.M. Woolfson, Astron. Geophys. 41(1), 1.12–1.19 (2000). https://doi.org/10.1046/j.1468-4004.2000.00012.x
M. Ya Marov, A.V. Kolesnichenko, Turbulence and Self-Organization (Modeling Astrophysical Objects) (Springer, New York, 2013)
Acknowledgements
The author is indebted to Prof. A. Prentice whose input and comments improved the quality of this paper.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Humi, M. Patterns Formation in a Self-Gravitating Isentropic Gas. Earth Moon Planets 121, 1–12 (2018). https://doi.org/10.1007/s11038-017-9512-y
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11038-017-9512-y