Rings, Shells, and Arc Structures Around B[e] Supergiants. I. Classical Tools of Nonlinear Hydrodynamics

Structures in circumstellar matter reflect both fast processes and quasi-equilibrium states. A geometrical diversity of emitting circumstellar matter is observed around evolved massive stars, in particular around B[e] supergiants. We recapitulate classical analytical tools of linear and nonlinear potential theory, such as Cole–Hopf transformation and Grad–Shafranov theory, and develop them further to explain the occurrence of the circumstellar matter structures and their dynamics. We use potential theory to formulate the nonlinear hydrodynamical equations and test dilatations of the quasi-equilibrium initial conditions. We find that a wide range of flow patterns can basically be generated and the timescales can switch, based on initial conditions, and lead to eruptive processes, reinforcing that the nonlinear fluid environment includes both quasi-stationary structures and fast processes like finite-time singularities. Some constraints and imposed symmetries can lead to Keplerian orbits, while other constraints can deliver quasi-Keplerian ones. The threshold is given by a characteristic density at the stellar surface.


INTRODUCTION
During their post-main sequence evolution, massive stars (M > 8 M ⊙ ) can undergo phases of enhanced mass loss and material ejections, which may lead to the formation of circumstellar shells and disks.One group of objects is particularly peculiar.These are the B[e] supergiants.The early findings by Zickgraf et al. (1985) of a hybrid character of the UV and optical spectra have led to the assumption of a two-component wind emanating from these objects with a classical line-driven wind in polar direction along with a much slower, cool and dense equatorial (disk-forming) outflow in equatorial direction.
Intense infrared excess emission points to the presence of significant amounts of hot circumstellar dust (Zickgraf et al. 1986), which imprints its specific emission features on the infrared spectra of these objects (Kastner et al. 2010).The existence of a disk-shaped dusty structure has been inferred from optical linear polarization observations (Magalhaes 1992;Melgarejo et al. 2001) and has later-on been reinforced for a few close-by Galactic objects by optical interferometric observations (Domiciano de Souza et al. 2007;Millour et al. 2011;Borges Fernandes et al. 2011;Cidale et al. 2012;Wheelwright et al. 2012a).
A warm and dense circumstellar disk provides an ideal environment for the formation of molecules, and rovibrational emission from hot molecular gas has been detected in the near-infrared from CO (McGregor et al. 1988;Morris et al. 1996;Miroshnichenko et al. 2005;Kraus et al. 2013Kraus et al. , 2014Kraus et al. , 2016Kraus et al. , 2020;;Wheelwright et al. 2012b;Oksala et al. 2012Oksala et al. , 2013;;Muratore et al. 2015;Sholukhova et al. 2015;Kourniotis et al. 2018) as well as from SiO (Kraus et al. 2015), and in the optical possibly from TiO (Zickgraf et al. 1989;Kraus et al. 2016;Torres et al. 2012Torres et al. , 2018) ) in a large number of objects (for a detailed description we refer to the review by Kraus 2019).The detection of significant enrichment of the disk material in the isotope 13 C, traced by intense emission of the molecule 13 CO, reinforces that the circumstellar matter of B[e] supergiants must have been released from the stellar surface and cannot be a remnant from star formation (Kraus 2009;Liermann et al. 2010).
Several theoretical approaches have been presented to explain the formation of dense outflowing disks from B[e] supergiants.Bjorkman & Cassinelli (1993) proposed that the winds emanating from the polar regions of rapidly rotating massive stars are bend towards the equator regions where they collide forming the so-called wind-compressed disk.This analysis considered a spherical shape of the star, but the inclusion of the nonradial forces occurring on the rotationally distorted stellar surface seems to prevent the formation of a windcompressed disk (Owocki et al. 1996).
Rapid stellar rotation facilitates another mechanism that might be considered, the rotation-induced bistability (Pelupessy et al. 2000).Due to the decrease in surface temperature from the pole to the equator caused by the rotation of the star (known as gravity darkening), the threshold temperature of ∼ 25 000 K for recombination of Fe iv into Fe iii might be crossed.Because Fe iii has significantly more lines suitable to drive the wind, a substantial increase in mass flux can be expected at lower temperature, that is, towards equatorial regions (Vink et al. 1999).However, the density enhancement that can be achieved by this scenario remains a factor 10-100 below the expectations from observations.When the bi-stability mechanism is combined with the slowwind solutions discovered by Curé (2004), the situation improves and a higher density contrast can be generated in the vicinity of the star (Curé et al. 2005), but for the price of an equatorial wind velocity that is 10-20 times higher than what has been inferred from observations.The effects of gravity darkening have also been utilized in computations of the latitude-dependent ionization structure in the winds of B[e] supergiants, and Kraus (2006) has shown that, despite the fact that rapid rotation alone leads to a lower wind density in equatorial directions as was shown by Owocki et al. (1996) and Maeder & Desjacques (2001), the wind material (in particular hydrogen and elements with a similar ionization potential) can still recombine in the equatorial wind region of these luminous objects leading to a zone of neutral gas confined to the equatorial plane, in which also molecules and dust can form.The existence of such neutral (in hydrogen) zones (or disk-like structures) has been concluded from the analysis of the observed line luminosity of the [O i] lines that arise in these disks (Kraus et al. 2007).Follow-up 2D models revealed that such recombination scenarios require very high equatorial stellar mass-loss rates (Zsargó et al. 2008).In a dif-ferent approach of 2D dense, viscous outflowing disks it was found that viscous heating dominates the innermost disk regions leading to extremely high temperatures within the disk mid-plane and to instabilities with significant waves or bumps in density and temperature (Kurfürst et al. 2018).
Many new observations have been carried out in the past years providing clearer insight with respect to the density distribution and the dynamics within the circumstellar (disk) matter.In particular, it has been found that the circumstellar material is confined in series of rings, arcs, or spiral-arm like structures revolving the central object on (quasi-)Keplerian orbits, rather than being spread over a disk in the classical picture (e.g., Kraus et al. 2010Kraus et al. , 2016Kraus et al. , 2023;;Millour et al. 2011;Aret et al. 2012;Maravelias et al. 2018;Torres et al. 2018).The arrangement of these rings is thereby unique for each object (Maravelias et al. 2018), and each ring can have a different density that follows not necessarily the usual radial distribution expected in an outflow.Moreover, these rings can have gaps or inhomogeneities, and they can be either stable in time (Kraus et al. 2016(Kraus et al. , 2023) ) or display temporal variabilities (Maravelias et al. 2018;Torres et al. 2018) including fading (Kraus et al. 2020), complete disappearance (Liermann et al. 2014), but occasionally also a sudden appearance of a new structure (Oksala et al. 2012) possibly caused by a pile-up of matter in a steadily decelerating outflow (Kraus et al. 2010).These findings indicate that the mass-loss from these stars is not a smooth process, but could be related to ejection phases, possibly triggered by instabilities acting in the strongly inflated envelopes of such massive and luminous objects (Glatzel & Kiriakidis 1993;Kiriakidis et al. 1993;Glatzel & Mehren 1996).
Motivated by this great diversity of circumstellar environments of B[e] supergiants ranging from stationary density distributions in the form of rotating rings with sometimes alternating densities, or arc-like features, to decelerating equatorial outflows with sudden pileup of matter, and the deficiency of existing models to describe them, we develop in this paper new perspectives with different hydrodynamical scenarios that might help understanding their formation1 and stationary structure.In particular, the formerly discussed scenarios and phenomena should be interconnected with basic, generic properties of fluid dynamics.We want to consider the following open questions from an abstract perspective: • How can a stationary, ideal fluid representation (without dissipative effects) of persistent matter (stable, quasi-stationary rings, arcs, complete or incomplete spiral structures) be constructed?
• How can quasi-stationary mass distributions and probably time-dependent flows (e.g. in the form of episodic mass loss of the central star) appear together?
• How can simplified single-fluid time-dependent and time-independent (stationary) velocity fields be constructed for such above mentioned cases, if not many detailed physical parameters are known?
• Do such time-dependent velocity fields exist at all for a stationary density distribution?What is their nature?
The paper is structured in the following way: In Section 2 we present the stationary 2D solution techniques but also the transition to non-stationary and 3D problems within incompressible ideal hydrodynamics (HD).In Section 3 we derive a non-linear Schrödinger-type equation as equivalent formulation of the compressible HD equations.In Section 4 we analyze a general 3D compressible flow on the basis of stationary stellar wind solutions.Our results are discussed and our conclusions are summarized in Section 5.

Basic equations of ideal hydrodynamics
The scenarios described in the introduction will be treated by using the basic equations of HD.These are given by the mass continuity equation (Equation (1)) and the Euler equation (Equation (2)), in which ρ is the mass density, v the gas velocity, p is the gas pressure which may include also the radiation pressure given that it can be assumed to be isotropic2 , and g = −∇φ is the gravitational acceleration of the star with the gravitational potential φ.Self-gravitation effects of the circumstellar matter are neglected.We focus on systems in which viscosity can be neglected, meaning that the length scales of the pressure force are much larger than the so-called deflection length, and the flows are super-sonic but not highly super-sonic (i.e.no shocks involved, see Frank et al. 1992).
To allow for a profound investigation of the physical aspects, we will split our analysis into two distinct physical extremes: incompressible velocity fields (Section 2.2) and irrotational potential velocity fields (Section 3).One should be aware of the fact that the incompressibility condition ∇•v = 0 is neither valid for classical viscous disks models, nor for classical wind solutions, as these usually assume isothermy.Our incompressibility model is not isothermal, and it makes shocks less likely to occur.

Stationary 3D incompressible flows
An interesting case that has emerged from the observations are the (presumably) Keplerian rotating rings detected around the B[e]SG star LHA 120-S 73 (Kraus et al. 2016).These rings did not show any measurable radial motion within the observation period spanning 16 years, justifying the assumption of a quasistationary model, i.e., time-dependent changes of physical quantities are small and can be neglected.Moreover, due to the preferential rotational and relaxed motion of the gas along closed orbits, the compressibility of the gas is negligible meaning that ∇•v = 0.In this case, the mass conservation Equation (1) simplifies to implying that the flow is always perpendicular to the gradient of the density (v ⊥ ∇ρ), respectively the density is constant along streamlines.We now introduce the streaming vector (or auxiliary flow vector, Nickeler & Fahr 2005;Nickeler et al. 2006) The total pressure, or as we will call this in the following, the Bernoulli-pressure Π, is defined by Applying the identity to the stationary Euler Equation (2) and the definitions for the streaming vector in Equations ( 4)-( 6), the resulting momentum equation in 3D can be written in the The Equations ( 5) and ( 8) form the set of incompressible HD equations that need to be solved.

Solutions in 2D
Due to flat or disk-like structuring and concentration of the matter around the star, we assume that the material is confined in a thin sheet around the equatorial plane.In general, circumstellar disks can be assumed to be symmetric around the mid-plane, having their maximum density and pressure along z = 0 so that ∂/∂z = 0. Our analysis is restricted to the xy-plane and fringe effects are neglected.The investigated scenario is not intended to generate an outflowing disk in the classical sense, but to find representations for revolving rings or arcs in a quasi-stationary state.In the limit, we assume a purely azimuthal flow and use a different geometry (basically 2D in Cartesian components, i.e. (x, y) = (R, Φ)), where in contrast most works about disks use the (R, z) coordinate system (i.e.∂/∂Φ = 0).
To satisfy the mass conservation Equation ( 5) we define the stream function ψ = ψ(x, y), the streaming vector w = ∇ψ×e z , and the mass density ρ = ρ(ψ).In the general case, ∇ψ and ∇φ are not parallel to each other almost everywhere, such that the stream function ψ and the gravitational potential φ can be regarded as coordinates replacing x and y, and Π is consequently an explicit function of ψ and φ.Inserting these relations into the momentum Equation ( 8) and expanding with respect to the basis vectors ∇ψ and ∇φ delivers From a comparison of coefficients we obtain one nonlinear Poisson-like partial differential equation for the stream function and one equation for the external gravitational potential This set of relations is an analogy to results in magnetohydrostatic (MHS) theory with gravity (Schindler et al. 1983), while the first quasi-or non-linear elliptic type equation has already been found and described in HD by Stokes (1848).Formal integration of the right-hand equation of Equation ( 11) results in a relation for the Bernoulli-pressure however, this formal integration does not solve the system completely.Inserting this Bernoulli-pressure into the left-hand equation of Equation ( 11) delivers where the prime denotes the derivative with respect to ψ.
To be able to solve this equation, it is necessary to know the function ρ(ψ), i.e. the density as a function of the stream function.While the stream function ψ is not known a priori, the density along each stream-line label (i.e. the value of the stream function) must be constant, but can vary across stream lines.With this knowledge, we are able to calculate the 2D density structure, as we can basically choose the density function ρ = ρ(ψ) arbitrarily.The choice of the density function implies a nonlinear feedback on the solution of the non-linear Laplace Equation ( 13).The solution of this Laplace equation finally delivers a stream function ψ = ψ(x, y), based on which the spatial density distribution ρ = ρ(x, y) can be computed3 .
Taking the most simple, nontrivial approach given by and assuming that Π 0 is constant, results in The relation between the density and the stream function, Equation ( 14), is chosen in such a way that the monotonicity of the density function, and therefore a unique relation, i.e. the bijective character of the density function, is guaranteed.We will recognize in the following that λ controls the influence of the gravitation on the geometry and dynamics of the flow: the larger λ, the larger the deviation from ordinary potential flows (see Equations ( 17) and ( 18) below).
To derive this influence we first of all have to facilitate the solution of Equation ( 15), by switching to Wirtinger calculus (e.g.Remmert 1991) and using the coordinates (or variables) u := x + iy and v := x − iy with i 2 = −1.Then the Poisson equation, Equation ( 15) can be written as where we inserted the definition of the gravitational potential given by φ = −GM * R , in which G is the gravitational constant, M * is the mass of the central star, and where R defined as R 2 = uv is the radial distance from the center of the star within the x-y-plane, with R ≥ R * , where R * denotes the radius of the star.
Integration of the solution of the Poisson Equation ( 16) with respect to the u-v coordinates delivers the stream function where ψ 1 and ψ 2 are free functions.With this stream function we can compute the streaming vector and the Bernoulli-pressure, where the thermal pressure p can then be calculated by subtracting the kinetic pressure.
The general solution (Equation ( 17)) shows an inhomogeneous part depending on the coupling constant λ.This simplified approach leads in a mathematical limit to a quasi-Kepler rotation.We assume that the homogeneous part of the general solution (Equation ( 17)), i.e. the meromorphic part, is zero.Therefore the solution for ψ and the derived streaming vector can be used to calculate the (rotational) velocity of the gas around the star where we have inserted for the density function the solution of Equation ( 14), i.e., ρ(ψ) = −λψ + ρ 0 = λ 2 GM * R + ρ 0 .For the limit λ, ρ 0 → 0 and ρ 0 /λ 2 → 0, or for ρ 0 = 0, the velocity is identical to the Kepler rotation.For a vanishing integration constant ρ 0 , the density at the stellar surface is given by ρ(R * ) = λ 2 GM * R * .While the density of the circumstellar matter in the vicinity of the star can take small values, depending on the choice of λ, it is not trivial to generate a Keplerian rotating disk or ring detached from the stellar surface.
Other cases, e.g., for large values of ρ 0 /λ 2 , lead to strong deviation from Kepler rotation close to the star.But, even in the latter case, the velocity can approach Keplerian behavior for sufficiently large values of R. It should be noted that our approach only represents the close-by circumstellar environment, i.e., regions between the stellar surface and the first ring or onset of the disc.The linear function chosen for the density ρ(ψ) might be regarded as the lowest order part of a Taylor expansion in the region of our interest of a complex density function, and is only meant to have a prescription from a low density, e.g.close to the star ('gap'), to high values (ring/disk).The lack of information of the real radial density distribution from observations hampers a more precise theoretical description.
Inserting the streaming vector, Equation ( 18), and the density function ρ(ψ) into the equation for the pressure function, Equation ( 12), Equation ( 6) delivers If ρ 0 would be zero, the density would be given by λ 2 GM * R, the pressure gradient (and therefore the pressure force) would vanish, and the motion of the gas would be purely ballistic (Keplerian).
2.2.2.On the existence of Π(ψ, φ) for one-dimensional equilibria In Section 2.2.1 we showed for two-dimensional equilibria that the pressure Π can be expressed as a function of the scalar fields ψ and φ with where Ω w is the vortex strength (or vorticity in 2D) associated with the flow vector w , which is defined by Ω w := ∇×w and reduces in our case to Ω w = Ω w e z .
In the one-dimensional case, e.g.∂/∂Φ = 0, where (x, y) ≡ (R, Φ) (with R and Φ being the radial coordinate and the azimuthal angle, respectively), it is not obvious at first that there is an incompressible stationary solution, characterized by Π(R), ψ(R) and φ(R), the vorticity Ω w (R), and the mass density ρ(R).
A possible proof has been proposed by Hornig (1996, private communication) and presented by Fleischer (1996), which we recapitulate for clarification: We consider the two-dimensional space IR 2 with the coordinates (ψ, φ), in which the solution ψ(R), φ(R) is a curve L, parameterized by R. On L the vector field We are looking for a potential Π on IR 2 such that V = ∇Π| L .We assume that there is a ε 0 -hose around L that does not overlap anywhere.Then there are local coordinates (s, ε), where s is the arc length along L and ε is the distance from L. These are orthogonal to L, e s •e ε | L = 0.There is then a decomposition of V (s) = V s (s)e s + V ε (s)e ε on L, and the potential Π, with Π ≡ 0 outside the tube, and Thus there is a pressure function Π(ψ, φ) whose partial derivatives according to the Equations ( 22) and ( 23) agree with the vorticity Ω w and the mass density ρ as functions of ψ and φ.Even if a great variety of geometrical structures with regard to flow patterns can be generated in 2D with the method described in Section 2.2.1, we know that even quasi-stationary structures have to undergo a formation phase that can be quite eruptive and therefore strongly time-dependent.This leads to the question whether quasi-stationary structures can be preserved in case of a non-linear time-dependent change of the system or relax respectively blow up.Or, equivalently, is it at least possible to extend the stationary models via 'slight amplitude modulation', i.e. the assumption that stationary fields can be multiplied with a specific purely time-dependent function, to deliver a time-dependent solution including the stationary solution?To answer this question we focus on the dynamical behavior in the non-linear case.We additionally assume that the system reacts instantaneously on non-local perturbations.

2.2.3.
The slight amplitude modulation is equivalent to posing the question if a time-separability from the spatial components guarantees regular solutions within a small, finite time interval.The regularity would at least make it plausible that the solutions are non-linearly stable.The separability of solutions of evolutionary equations with respect to time is possible for many examples in continuum and quantum physics (see, e.g., Galaktionov & Svirshchevskii 2006).The general solutions of such kind of equations can then be constructed by linear superposition of multiple summands constructed by such a separation ansatz, or by nonlinear superposition of a similar set of summands.
We assume a time-separability and take the general incompressible 3D representation into account.Considering the Euler equation for vanishing time derivative ∂ρ/∂t ≈ 0 and external force, and using the definition of the stream vector (Equations ( 4) and ( 5)), we obtain With w = w 1 (t) w 0 (x ), and analogously for the total pressure Π, we can rewrite the second form of the Euler Equation ( 27) as where the dot in ẇ1 denotes the derivative with respect to time.Assuming that ∃ (w 0 , Π 0 , ρ) with with the integration constant w 10 .For w 10 > 0, t 0 > 0 and t 0 > t this solution constitutes a blow-up solution or finite-time singularity for t → t 0 , whereas in the limit t → −∞ the velocity and pressure decay.We would like to emphasize that Equation ( 29) does not represent an equilibrium solution, but poses an additional constraint for the stationary part of the flow, respectively the steady-state flow pattern.
If we take also gravity into account for the above calculations then the time-independent form of the momentum equation is where e r is the unit vector directing radially outwards from the center of gravity, and M * (t) is the stellar mass changing over time.For t 0 < 0 and t = 0 the relation ẇ1 = w 10 /t 2 0 = 1 fulfills the condition that M * (t = 0) = M * .Then it follows that for t > 0 and t → ∞ the mass, pressure and velocity decrease, whereas they increase for t 0 > 0, t 0 > t and t → t 0 .While the former can be interpreted with mass loss, the latter would imply mass accretion.It should be noted that both scenarios only hold for a meaningful (small) time interval, as the temporal variation of the velocity field is given by a (pure) dilatation.
The assumptions leading to Equations ( 27)-( 31) reflect a nonlinear perturbation of a system, whose original stationary (geometrical) structure should be preserved, whereas the amplitudes of the stationary field components should vary only slowly.The calculations show that the time-independent field components can in fact no longer satisfy any stationary equation, and that even the time-dependent amplitudes diverge, i.e. develop finite-time singularities, or decay completely.Although the dynamics is limited to incompressible flows in these considerations, the behavior of the time-dependent amplitudes implies an eruptive mass loss (or mass gain) of the star.This means that within a short period of time the stellar mass loss (or gain) would be considerably enhanced when the system is forced by non-linear variations (perturbations) which can change the character of the system, driving it away from a relaxed incompressible flow.Our analysis shows that time-dependent incompressible flows around a stationary state cannot exist but they lead either to a completely decaying flow or to a blow up.

POTENTIAL FLOWS IN 3D
In contrast to the studies in Section 2.2, the timedependent change in velocity field is now included.Also, we drop the restriction to incompressible flows.Instead, we assume that the flow is irrotational.Moreover, we replace the stream function with a scalar potential, which allows us to easily implement the time-dependence of the velocity field.
In contrast to the stream function model used in Section 2.2, reflecting the isocontours of the density distribution in a relaxed flow, we now utilize a scalar velocity potential.The surfaces of constant potential represent a family of surfaces, which has an affinity to the family of surfaces of constant radial coordinate.This would favor a radial outflow, instead of the more azimuthal flows described by the stream function model in the previous section.In case of a time dependency, the potential describes expanding (or eventually shrinking) surfaces, to which the velocity vector is perpendicular, thus pointing in a certain sense 'radially' outwards (or inwards).The density, decreasing outwards (eventually inwards), should be therefore, at least locally, diffeomorphic with respect to the potential, as will be introduced and explained in the following.
First of all, maintaining the condition of a stationary density distribution, the equation of mass conservation reduces again to The Euler Equation ( 2) can be re-written, using the identity Equation ( 7), as Introducing a 3D velocity potential v := ∇ϕ, leading to ∇×v = 0, and utilizing a barotropic law p = p(ρ), Equations ( 32) and ( 33) take the following form The latter, Equation ( 35), is the classical Bernoulli equation for compressible, unsteady flows without viscosity.Without loss of generality, we have incorporated a possible time-dependent integration constant into the function ϕ.This equation needs further integration, as the mass conservation Equation ( 34) must also be fulfilled.
The study of possible temporal evolutions of the Euler Equation ( 33), together with the mass continuity equation Equation ( 34) is subject of Section (3.1).
In analogy to the incompressibility representation (ρ = ρ(ψ)), we introduce the affinity between isosurfaces of the velocity potential and density, namely ρ = ρ(ϕ), where ρ ′ = ρ ′ (ϕ) = dρ/dϕ.Adding an explicit time dependence to the mass continuity Equation ( 34) results in the time-depending continuity equation Using this equation, we substitute either the quadratic derivative of the potential ϕ or its time derivative in the Euler Equation ( 35) and obtain The function G can be identified as the specific enthalpy.Equation ( 37) is a non-linear diffusion equation, similar to the non-linear Schrödinger equation, and Equation ( 38) is a convection-diffusion-type equation.
In the following we first discuss a special case with a specific time dependency where only ϕ is time dependent and ρ is stationary and depends only on the spatial parts of ϕ (Section 3.1), whereas in Section 3.2 we return to the general case of a density depending explicitly on the time-dependent velocity potential.

Separable time dependence and neglection of gravity
We want to search for the general time dependence of (circumstellar) flows in order to find out whether instabilities or collapse processes necessarily occur, even if the system stipulates a stationary density at far distances from the central star as observations suggest.In particular, we are interested in an expansion of a solution around an equilibrium state, and this can be naturally achieved by the separation of the temporal and spatial parts of the corresponding fields.Therefore, we assume a separability of time dependency for the pressure and velocity potential, consider a stationary density, depending on the spatial part of the velocity potential, and neglect gravity4 , i.e. p = p 0 (t) p 1 (ϕ 1 ), ϕ = ϕ 0 (t) ϕ 1 (x, y, z) and ρ = ρ(ϕ 1 ).The chosen time-separability for the pressure enables a barotropic law at any time.Inserting this separability ansatz into the conservation of mass equation for stationary density, Equation ( 34), we obtain the relation With the 3D velocity potential v = ∇ϕ, and neglecting gravity, Equation ( 33) reduces to Inserting the separation ansatz for the pressure and the velocity potential leads to ) .
(42) The right-hand side can be written as With this, Equation ( 42) can be formally integrated, and if we replace the term (∇ϕ 1 ) 2 by Equation ( 40), we obtain As one can see by inspection of Equation ( 39), F is up to a factor p 0 (t) identical to the specific enthalpy G.
Inserting Equation (45) into Equation ( 44), we get Equation ( 47) is only valid for φ0 = 0. To achieve the separability of this equation into spatial and timedependent parts, we factorize the first terms in Equation ( 47) and recognize that it is indispensable that the term p0 ϕ 2 0 must be constant.We label this constant as c 1 and obtain Equation ( 48) can only be solved if delivering a finite-time singularity for the timedependent part of the velocity potential.The constant ϕ 00 describes the amplitude of the velocity potential imprinted onto the spatial part of the velocity potential for the time t = 0.
For the spatial part we find Equation ( 50) is basically a nonlinear elliptic partial differential equation whose solution reflects the spatial part of the scalar velocity potential ϕ 1 and thus determines the structure of the entire time-dependent flow.To solve this equation requires the knowledge of the specific enthalpy F respectively G which results from the specification of the density function ρ(ϕ 1 ) along with pressure dependency p 1 (ρ) or explicitly p 1 (ϕ 1 ).
A detailed analysis of this equation is beyond the scope of the current paper, however, we wish to draw the attention to the fact that, if the time separability is requested along with a time independent density distribution (stationarity), then according to Equation (49) for ϕ 00 > 0, t 0 > 0 and t 0 > t the system inevitably develops into a finite-time singularity for t → t 0 , whereas in the limit t → −∞ the velocity and pressure decay.
For a simplified case we can find a subspace.For example, for the case φ0 = 0, Equation ( 46) takes the following form This equation can only be fulfilled for p 0 = const.In this case, the equation takes the form of a quasi-linear partial differential equation which can also be derived directly from the non-linear Schrödinger Equation (37).Equation ( 52) is a non-linear Poisson equation for which exact analytical solutions are known in 2D.

General time dependent approach and generalized Cole-Hopf transformation
We return now to the system of non-linear diffusion equations (Equations ( 36)-( 39)) which we want to reformulate in a more compact form.For this, we use a generalized form of the Cole-Hopf transformation.The basic form of this transformation type has been invented by Hopf (1950) and Cole (1951).
In the case of non-linear diffusion equations the spatial part of the differential operators can have the form Let us assume ϕ = F (Λ), we can rewrite Equation ( 53) as where the primes at F denote derivatives with respect to Λ.The transformation ϕ = F (Λ) should be now done in such a way that terms with (∇Λ) 2 are eliminated.This demand can be formulated as To solve Equation ( 55 Inserting the general solution (Equation ( 56)) along with the functions a, b, c into Equation (57) leads to The last equation of Equation ( 59) is a nonlinear diffusion equation and is the Cole-Hopf transformed mass continuity equation.For this kind of non-linear diffusion or Schrödinger equations (Equation ( 59) and Equation (37)) it is known that they can have blow-up solutions (see, e.g., Galaktionov & Vázquez 2002).
Next, we transform Equation ( 38).Due to the mathematical similarity of Equation ( 38) and Equation ( 53) we can use the transformation Equation ( 56) for This leads to The advantage of carrying out the transformation is that Equation (61) does not contain time derivatives and quadratic terms anymore, and that it is again a nonlinear Poisson equation which can be solved with standard procedures.The only problematic and non-trivial issue is that ρ(ϕ) must be chosen adequately in order to find solutions of the equation, based on which it will be possible to obtain the velocity field, pressure and density distribution.

SOLUTIONS OF EULER EQUATION FOR SPECIAL, PERSISTENT GEOMETRICAL FLOW PATTERNS
Now we drop the condition for incompressible or irrotational flows and treat the hydrodynamical problem of non-linear and non-local instabilities from a more general point of view.It is possible to implement the gravitation of the star into the pressure in case that either p can be considered as a function of ρ (barotropic fluid), or ρ as a function of φ, although in the former case it then be necessary to slightly modify the equation of motion.Without significant loss of generality, we therefore use only one pressure gradient force and write the Euler equation for the stellar wind in the following form For this, the equilibrium solution, which one may also consider as the 'ground state' of the system, is If we assume the separability for the pressure p = p 1 (t) p 0 (x ) and analogously also for the density ρ = ρ 1 (t)ρ 0 (x ) and for the velocity field v = v 1 (t) v 0 (x ), then we will show in the following that for some choices of constraints there can exist regular solutions, and other choices can lead to a blow-up of the stationary nonlinear spatial solution.
The separation implies that locally, i.e. inside some small subset of a time interval, the fluid variables are not to be considered as fast evolving, around a known stationary wind solution.The mass continuity equation can then be written as and with the conditions of Equation ( 63) it follows that ρ1 = 0 without loss of generality.The Euler equation ( 62) together with the initial condition Equation ( 63) can be expressed in the form given in Equation (66).Taking first the divergence operator on both sides (∇•), and second the curl operator (∇×), one receives ⇒ ∃ ϕ 0 , ∆ϕ 0 = 0 , v 00 = const., with The condition v 00 < 0 is subject to the physical approach that the direction of the velocity should be antiparallel to the pressure gradient in the relation between the mass-current vector ρ 0 v 0 , and the pressure gradient Equation ( 68).The assumption that v 00 < 0 reflects the decrease of the pressure from upwind to downwind direction, which usually drives a stellar wind or emulates 'decretion'.Equation ( 69) represents an ordinary differential equation for the time dependent dynamics of the flow, based on the equilibrium solution Equation ( 63).Already under simplified assumptions, like p 1 = const., we can find solutions with finite-time singularities (blow-up solutions), but these solutions originate from a finite-time singularity and converge towards the equilibrium solution.In contrast to the decretion solution, the accretion solution would originate from the equilibrium solution and converge to a finite-time singularity.In the case p 1 = 0 it can clearly be seen and easily calculated that v 1 ∝ 1/(t 0 − t).This solution must be excluded, as this would imply that the time perturbation reduces the complete pressure to zero for all times, implying a 'pressure-less fluid', i.e., a completely force-free flow.
For p 1 = const.= 0 we calculate the formal solution: leading to (71) For t > t 0 and v 1 , v 10 > √ p 1 Equation ( 71) can be reformulated as Then, the solution of Equation ( 73) implies that the velocity decays with time and relaxes into an equilibrium state.
If we assume that p 1 is not constant anymore but p 1 = p 1 (t) can be imposed, then Equation (69) turns into a ordinary differential equation of Riccati type.
As the problem is nonlinear anyway, the function p 1 can be regarded, especially in this situation of a dependency of only one variable, namely t, as function of v 1 , i.e. p 1 = p 1 (v 1 ).For certain constraints we are able to show that equilibria could converge to blow-up solutions.The assumption that p 1 (v 1 ) ≥ 0 guarantees that the pressure is always positive.For a continuous outflow we must also guarantee that v1 is always larger than zero, if we exclude oscillatory phases of acceleration and deceleration.These conditions are fulfilled by p 1 − v 2 1 ≥ 0 (see Equation ( 69)).Thus, it is further necessary that p 1 is not only positive, and a sufficient criterion to guarantee the validity of p 1 − v 2 1 ≥ 0 is that there exists a p 10 > 1 and p 1 = p 10 v 2 1 such that Integration of this equation results in with v 1 (t 0 ) = v 10 .For t → t crit the system develops a blow-up, where To summarize, the choice of p 1 (v 1 ) determines the temporal course of the flow and we have shown two contrasting examples, namely a decaying perturbation (for p 1 = const.)and a perturbation that causes a blow-up of the system (for p 1 = p 10 v 2 1 ).More complex choices of p 1 (v 1 ), e.g. with higher orders of the polynomial, can lead to correspondingly diverse results such as multiple blow-ups and decays.Other time courses are quite possible and will occur, but one should keep in mind that the blow-up solutions can be damped by normal and anomalous dissipative processes, leading to a regular time-dependent behavior.However, the blow-up solutions indicate that non-linear instabilities can develop in quasi-ideal systems and hence can lead to abrupt changes of the typical temporal scales.

DISCUSSION AND CONCLUSIONS
The environments of certain types of evolved massive stars, such as the B[e] supergiants, display indications for disks, or multiple ring-like structures of yet unknown origin, which cannot be reconciled with the classical theory of a (viscous) outflowing disk (e.g., Lee et al. 1991;Okazaki 2001;Kurfürst et al. 2018).These configurations can be either steady (Kraus et al. 2016(Kraus et al. , 2023) ) or display temporal variability (Maravelias et al. 2018;Torres et al. 2018).Therefore we have two problems, namely (i) either long time scales (equilibria), or systems close to or converging to equilibrium states, or (ii) we need short time scales to explain sudden appearance of new structures (such as detected by Oksala et al. 2012) resulting from non-linear instabilities (occurring either within the circumstellar matter or already in the stellar atmosphere) which might be connected to collapse processes (finite-time singularities).To study the high diversity of geometrical structures and their possible formation mechanism we recourse to elementary mathematical tools of classical hydrodynamics such as linear and non-linear potential theory.
For the first case of long time scales, complex geometries of streamlines are known and can be constructed from potential representation of equilibrium fields, e.g.potential fields (Laplace equation), or by non-linear solutions of Grad-Shafranov-type equations, where closed and open streamline configurations (arcs, rings, radial structures) are possible (see, e.g., Nickeler et al. 2013Nickeler et al. , 2014)).In our analysis we include flows and gravity and find that the equations are of similar structure, which means that they can be solved in an analogous manner.We could prove that for a non-pressureless gas we can find quasi-Keplerian or even Keplerian rotation of the circumstellar matter.
To bridge the gap from these equilibrium structures to short time scales a 'homotopic' modulation, i.e. a dilatation of the quasi-equilibria is performed, leading eventually to restrictions for the equilibrium values and to ordinary differential equations as constraints for the time dependencies of the fluid equations.These ordinary differential equations in time lead either to regular (decay or growth, or relaxation) or non-regular time dependencies where the latter includes finite-time singularities, which are known from (magneto-)hydrodynamical systems (e.g., Klapper et al. 1996;Nickeler & Karlický 2008), and a broad overview about mathematical solution techniques for different physical problems can be found in Galaktionov & Svirshchevskii (2006, and references therein).
For future investigations, non-separable time dependencies can be taken into account to get a better understanding of the variability and changes of the circumstellar matter.Moreover, for studies of the conditions under which rings and arcs can form within the circumstellar material it will be important to analyze in more detail the influence of different conformal mappings or multipole components, being able to create arcs, rings, and spiral-arm like structures.Such an analysis will then allow to derive the properties of the gas (density, temperature, and emissivity) which can be compared to observed quantities.
Preliminary discussion: Do time-dependent incompressible flows exist in 2D or in 3D ?
), it is necessary that the function a/b only explicitly depends on Λ.The transformation F can be derived by integration of the condition Equation (55) which has the general solution of coefficients between Equation (36) and Equation (53), we can derive the functions a, b and c