Maximum dust-to-gas mass flux ratio in spherically expanding dusty-gas flow.
Authors: Zakharov Vladimir , Bykov Nikolay , Rodionov Alexander , Ivanovski Stavro , Della Corte Vincenzo , Rotundi Alessandra 
Published: July 2024

Maximum dust-to-gas mass flux ratio in spherically expanding dusty-gas flow.In the context of an increasing number of complex multiparametric dusty-gas coma models it is convenient to address the analysis of complex physical phenomena in a step-by-step manner. We construct a set of elementary models with a minimum number of parameters selected to represent the key processes acting in a dusty gas coma. Such kind of models enables studies of generic processes occurring in variety of particular cases in order to reveal their characteristic features.In many models of dusty-gas atmosphere of comets the dust-to-gas mass flux ratio is a free parameter. In order to constraint the range of possible values of this parameter, in the present work we propose an estimate of its maximum value. We consider a &#8220;pedagogic&#8221; case which assumes a homogeneous spherical nucleus and spherical expansion of dusty-gas (as an elementary model of the coma). The dust grains are assumed to be homogeneous spheres with size distribution given by a power law n(ad) &#8733; ad-&#946;.After ejection from the surface dust grain is accelerated by the gas flow in the gravitational field of the nucleus. In the present consideration we assume that the radiation pressure due to solar illumination of the grain is negligibly small. With increasing distance to the nucleus the gas drag and the nucleus gravitational attraction decrease and decoupling of gas and dust flows occurs at some distance. In absence of other forces, the dust grains continue motion with constant &#8220;terminal velocity&#8221;. As was shown in [1], the dust grains reach 90% of the terminal velocity at cometocentric distance about 6 comet radii.The numerical integrated terminal velocity of dust grains in a spherically expanding flow is given in [1] for a wide range of conditions. The terminal velocity (vd) can be approximated as:vd(ad)/vmax=Ai IvBi, i=1,2,3,where Iv is a dimensionless parameter characterizing the efficiency of entrainment of the particle within the gas flow (see [1]), vmax is the theoretical maximal velocity of gas expansion, and Ai, Bi are coefficients of approximation defined for three non-overlapping ranges of parameter Iv variation.To escape the gravitational attraction of the nucleus terminal velocity of the dust particles should exceed the escape velocity vesc. This condition combined with the approximation for the terminal velocity allows us to derive the maximum size of the grains that reach escape velocity.The dust mass loss rate qd of grain size ad and the gas production rate qg are related via:qd(ad)=&#967;fmd(ad) qg dad,where &#967;=Qd/qg is the dust-to-gas mass loss rate, Qd is the dust mass loss rate integrated over all grain sizes, and fmd is the normalized mass distribution of dust grains. &#160;Considering the foregoing, for the cometocentric distances greater than 6 comet radii it is easy to derive the energy fluxes of gas and dust flows.The dust flow gets energy from the gas flow. Therefore, the maximum dust-to-gas mass loss rate corresponds a certain fraction of gas flow energy (&#954;