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

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

Identification of Changes Needed in Supermarket Design for Energy Demand Reduction

Supermarkets use 3 percent of UK energy. To satisfy building regulations supermarket buildings are modeled in considerable detail. Lighting, occupancy, and small electrical energy impacts are included in this modeling. However, refrigeration energy is not, as it is classified as process energy rather than building related. Refrigeration energy, which can be very significant, is therefore currently unregulated and as a result, heat transfers related to refrigeration cabinets are typically not incorporated in modeling of the building at design stage. This paper explores the comparative energy demands of supermarket stores modeled, using a simple first order dynamic model, executed on Excel, and optimized firstly with, and secondly without, the cooling effect of refrigeration cabinets included in the model. A recently built supermarket is modeled. Results suggest that the energy demand of a new store could be reduced by 15 to 25 percent by improvement of the building envelope design with process energy included in the modeling

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.

Vanishing viscosity limit of navier-stokes equations in gevrey class

In this paper we consider the inviscid limit for the periodic solutions to Navier-Stokes equation in the framework of Gevrey class. It is shown that the lifespan for the solutions to Navier-Stokes equation is independent of viscosity, and that the solutions of the Navier-Stokes equation converge to that of Euler equation in Gevrey class as the viscosity tends to zero. Moreover the convergence rate in Gevrey class is presented

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

A faster algorithm for smoothed particle hydrodynamics with radiative transfer in the flux-limited diffusion approximation

We describe a new, faster implicit algorithm for solving the radiation hydrodynamics equations in the flux-limited diffusion approximation for smoothed particle hydrodynamics. This improves on the method elucidated in Whitehouse & Bate by using a Gauss-Seidel iterative method rather than iterating over the exchange of energy between pairs of particles. The new algorithm is typically many thousands of times faster than the old one, which will enable more complex problems to be solved. The new algorithm is tested using the same tests performed by Turner & Stone for ZEUS-2D, and repeated by Whitehouse & Bate.Comment: 12 pages, 7 figures. Accepted for publication in MNRA

Optimal prediction for moment models: Crescendo diffusion and reordered equations

A direct numerical solution of the radiative transfer equation or any kinetic equation is typically expensive, since the radiative intensity depends on time, space and direction. An expansion in the direction variables yields an equivalent system of infinitely many moments. A fundamental problem is how to truncate the system. Various closures have been presented in the literature. We want to study moment closure generally within the framework of optimal prediction, a strategy to approximate the mean solution of a large system by a smaller system, for radiation moment systems. We apply this strategy to radiative transfer and show that several closures can be re-derived within this framework, e.g. $P_N$, diffusion, and diffusion correction closures. In addition, the formalism gives rise to new parabolic systems, the reordered $P_N$ equations, that are similar to the simplified $P_N$ equations. Furthermore, we propose a modification to existing closures. Although simple and with no extra cost, this newly derived crescendo diffusion yields better approximations in numerical tests.Comment: Revised version: 17 pages, 6 figures, presented at Workshop on Moment Methods in Kinetic Gas Theory, ETH Zurich, 2008 2 figures added, minor correction

Analytic solutions and Singularity formation for the Peakon b--Family equations

Using the Abstract Cauchy-Kowalewski Theorem we prove that the $b$-family equation admits, locally in time, a unique analytic solution. Moreover, if the initial data is real analytic and it belongs to $H^s$ with $s > 3/2$, and the momentum density $u_0 - u_{0,{xx}}$ does not change sign, we prove that the solution stays analytic globally in time, for $b\geq 1$. Using pseudospectral numerical methods, we study, also, the singularity formation for the $b$-family equations with the singularity tracking method. This method allows us to follow the process of the singularity formation in the complex plane as the singularity approaches the real axis, estimating the rate of decay of the Fourier spectrum