A Continuum Description of Rarefied Gas Dynamics (I)--- Derivation From Kinetic Theory

We describe an asymptotic procedure for deriving continuum equations from the kinetic theory of a simple gas. As in the works of Hilbert, of Chapman and of Enskog, we expand in the mean flight time of the constituent particles of the gas, but we do not adopt the Chapman-Enskog device of simplifying the formulae at each order by using results from previous orders. In this way, we are able to derive a new set of fluid dynamical equations from kinetic theory, as we illustrate here for the relaxation model for monatomic gases. We obtain a stress tensor that contains a dynamical pressure term (or bulk viscosity) that is process-dependent and our heat current depends on the gradients of both temperature and density. On account of these features, the equations apply to a greater range of Knudsen number (the ratio of mean free path to macroscopic scale) than do the Navier-Stokes equations, as we see in the accompanying paper. In the limit of vanishing Knudsen number, our equations reduce to the usual Navier-Stokes equations with no bulk viscosity.Comment: 16 page

Radiative Transfer in Star Formation: Testing FLD and Hybrid Methods

We perform a comparison between two radiative transfer algorithms commonly employed in hydrodynamical calculations of star formation: grey flux limited diffusion and the hybrid scheme, in addition we compare these algorithms to results from the Monte-Carlo radiative transfer code MOCASSIN. In disc like density structures the hybrid scheme performs significantly better than the FLD method in the optically thin regions, with comparable results in optically thick regions. In the case of a forming high mass star we find the FLD method significantly underestimates the radiation pressure by a factor of ~100.Comment: 4 Pages; to appear in the proceedings of 'The Labyrinth of Star Formation', Crete, 18-22 June 201

Vanishing viscosity limits for the degenerate lake equations with Navier boundary conditions

The paper is concerned with the vanishing viscosity limit of the two-dimensional degenerate viscous lake equations when the Navier slip conditions are prescribed on the impermeable boundary of a simply connected bounded regular domain. When the initial vorticity is in the Lebesgue space $L^q$ with $2<q\le\infty$, we show the degenerate viscous lake equations possess a unique global solution and the solution converges to a corresponding weak solution of the inviscid lake equations. In the special case when the vorticity is in $L^\infty$, an explicit convergence rate is obtained

Fluctuation-Response Relations for Multi-Time Correlations

We show that time-correlation functions of arbitrary order for any random variable in a statistical dynamical system can be calculated as higher-order response functions of the mean history of the variable. The response is to a ``control term'' added as a modification to the master equation for statistical distributions. The proof of the relations is based upon a variational characterization of the generating functional of the time-correlations. The same fluctuation-response relations are preserved within moment-closures for the statistical dynamical system, when these are constructed via the variational Rayleigh-Ritz procedure. For the 2-time correlations of the moment-variables themselves, the fluctuation-response relation is equivalent to an ``Onsager regression hypothesis'' for the small fluctuations. For correlations of higher-order, there is a new effect in addition to such linear propagation of fluctuations present instantaneously: the dynamical generation of correlations by nonlinear interaction of fluctuations. In general, we discuss some physical and mathematical aspects of the {\it Ans\"{a}tze} required for an accurate calculation of the time correlations. We also comment briefly upon the computational use of these relations, which is well-suited for automatic differentiation tools. An example will be given of a simple closure for turbulent energy decay, which illustrates the numerical application of the relations.Comment: 28 pages, 1 figure, submitted to Phys. Rev.

Random Walks on a Fluctuating Lattice: A Renormalization Group Approach Applied in One Dimension

We study the problem of a random walk on a lattice in which bonds connecting nearest neighbor sites open and close randomly in time, a situation often encountered in fluctuating media. We present a simple renormalization group technique to solve for the effective diffusive behavior at long times. For one-dimensional lattices we obtain better quantitative agreement with simulation data than earlier effective medium results. Our technique works in principle in any dimension, although the amount of computation required rises with dimensionality of the lattice.Comment: PostScript file including 2 figures, total 15 pages, 8 other figures obtainable by mail from D.L. Stei

Nonequilibrium corrections in the pressure tensor due to an energy flux

The physical interpretation of the nonequilibrium corrections in the pressure tensor for radiation submitted to an energy flux obtained in some previous works is revisited. Such pressure tensor is shown to describe a moving equilibrium system but not a real nonequilibrium situation.Comment: 4 pages, REVTeX, Brief Report to appear in PRE Dec 9

Numerical study of oscillatory regimes in the Kadomtsev-Petviashvili equation

The aim of this paper is the accurate numerical study of the KP equation. In particular we are concerned with the small dispersion limit of this model, where no comprehensive analytical description exists so far. To this end we first study a similar highly oscillatory regime for asymptotically small solutions, which can be described via the Davey-Stewartson system. In a second step we investigate numerically the small dispersion limit of the KP model in the case of large amplitudes. Similarities and differences to the much better studied Korteweg-de Vries situation are discussed as well as the dependence of the limit on the additional transverse coordinate.Comment: 39 pages, 36 figures (high resolution figures at http://www.mis.mpg.de/preprints/index.html

The phase shift of line solitons for the KP-II equation

The KP-II equation was derived by [B. B. Kadomtsev and V. I. Petviashvili,Sov. Phys. Dokl. vol.15 (1970), 539-541] to explain stability of line solitary waves of shallow water. Stability of line solitons has been proved by [T. Mizumachi, Mem. of vol. 238 (2015), no.1125] and [T. Mizumachi, Proc. Roy. Soc. Edinburgh Sect. A. vol.148 (2018), 149--198]. It turns out the local phase shift of modulating line solitons are not uniform in the transverse direction. In this paper, we obtain the $L^\infty$-bound for the local phase shift of modulating line solitons for polynomially localized perturbations

Instabilities in the Envelopes and Winds of Very Massive Stars

The high luminosity of Very Massive Stars (VMS) means that radiative forces play an important, dynamical role both in the structure and stability of their stellar envelope, and in driving strong stellar-wind mass loss. Focusing on the interplay of radiative flux and opacity, with emphasis on key distinctions between continuum vs. line opacity, this chapter reviews instabilities in the envelopes and winds of VMS. Specifically, we discuss how: 1) the iron opacity bump can induce an extensive inflation of the stellar envelope; 2) the density dependence of mean opacity leads to strange mode instabilities in the outer envelope; 3) desaturation of line-opacity by acceleration of near-surface layers initiates and sustains a line-driven stellar wind outflow; 4) an associated line-deshadowing instability leads to extensive small-scale structure in the outer regions of such line-driven winds; 5) a star with super-Eddington luminosity can develop extensive atmospheric structure from photon bubble instabilities, or from stagnation of flow that exceeds the "photon tiring" limit; 6) the associated porosity leads to a reduction in opacity that can regulate the extreme mass loss of such continuum-driven winds. Two overall themes are the potential links of such instabilities to Luminous Blue Variable (LBV) stars, and the potential role of radiation forces in establishing the upper mass limit of VMS.Comment: 44 pages, 13 figures. Chapter to appear in the book "Very Massive Stars in the Local Universe", Springer, J.S. Vink, e

Moment Closure - A Brief Review

Moment closure methods appear in myriad scientific disciplines in the modelling of complex systems. The goal is to achieve a closed form of a large, usually even infinite, set of coupled differential (or difference) equations. Each equation describes the evolution of one "moment", a suitable coarse-grained quantity computable from the full state space. If the system is too large for analytical and/or numerical methods, then one aims to reduce it by finding a moment closure relation expressing "higher-order moments" in terms of "lower-order moments". In this brief review, we focus on highlighting how moment closure methods occur in different contexts. We also conjecture via a geometric explanation why it has been difficult to rigorously justify many moment closure approximations although they work very well in practice.Comment: short survey paper (max 20 pages) for a broad audience in mathematics, physics, chemistry and quantitative biolog