Critical Collapse of an Ultrarelativistic Fluid in the Γ→1\Gamma\to 1 Limit

In this paper we investigate the critical collapse of an ultrarelativistic perfect fluid with the equation of state $P=(\Gamma-1)\rho$ in the limit of $\Gamma\to 1$. We calculate the limiting continuously self similar (CSS) solution and the limiting scaling exponent by exploiting self-similarity of the solution. We also solve the complete set of equations governing the gravitational collapse numerically for $(\Gamma-1) = 10^{-2},...,10^{-6}$ and compare them with the CSS solutions. We also investigate the supercritical regime and discuss the hypothesis of naked singularity formation in a generic gravitational collapse. The numerical calculations make use of advanced methods such as high resolution shock capturing evolution scheme for the matter evolution, adaptive mesh refinement, and quadruple precision arithmetic. The treatment of vacuum is also non standard. We were able to tune the critical parameter up to 30 significant digits and to calculate the scaling exponents accurately. The numerical results agree very well with those calculated using the CSS ansatz. The analysis of the collapse in the supercritical regime supports the hypothesis of the existence of naked singularities formed during a generic gravitational collapse.Comment: 23 pages, 16 figures, revised version, added new results of investigation of a supercritical collapse and the existence of naked singularities in generic gravitational collaps

Can finite-frequency effects be accounted for in ray theory surface wave tomography?

International audience[ 1] We present a series of synthetic tests showing that regional surface wave tomographies with a dense path coverage of the target region can be safely conducted under ray theory because the shortcomings of ray theory in considering finite-frequency effects can be counterbalanced by a physically-based regularization of the inversion. In particular, we show that with ray theory applied under the above conditions, it is possible to detect heterogeneities with length scales smaller than the wavelength of the data set

Type II critical phenomena of neutron star collapse

We investigate spherically-symmetric, general relativistic systems of collapsing perfect fluid distributions. We consider neutron star models that are driven to collapse by the addition of an initially "in-going" velocity profile to the nominally static star solution. The neutron star models we use are Tolman-Oppenheimer-Volkoff solutions with an initially isentropic, gamma-law equation of state. The initial values of 1) the amplitude of the velocity profile, and 2) the central density of the star, span a parameter space, and we focus only on that region that gives rise to Type II critical behavior, wherein black holes of arbitrarily small mass can be formed. In contrast to previously published work, we find that--for a specific value of the adiabatic index (Gamma = 2)--the observed Type II critical solution has approximately the same scaling exponent as that calculated for an ultrarelativistic fluid of the same index. Further, we find that the critical solution computed using the ideal-gas equations of state asymptotes to the ultrarelativistic critical solution.Comment: 24 pages, 22 figures, RevTeX 4, submitted to Phys. Rev.

Accurate discretization of advection-diffusion equations

We present an exact mathematical transformation which converts a wide class of advection-diffusion equations into a form allowing simple and direct spatial discretization in all dimensions, and thus the construction of accurate and more efficient numerical algorithms. These discretized forms can also be viewed as master equations which provides an alternative mesoscopic interpretation of advection-diffusion processes in terms of diffusion with spatially varying hopping rates

Consistent thermodynamic derivative estimates for tabular equations of state

Numerical simulations of compressible fluid flows require an equation of state (EOS) to relate the thermodynamic variables of density, internal energy, temperature, and pressure. A valid EOS must satisfy the thermodynamic conditions of consistency (derivation from a free energy) and stability (positive sound speed squared). When phase transitions are significant, the EOS is complicated and can only be specified in a table. For tabular EOS's such as SESAME from Los Alamos National Laboratory, the consistency and stability conditions take the form of a differential equation relating the derivatives of pressure and energy as functions of temperature and density, along with positivity constraints. Typical software interfaces to such tables based on polynomial or rational interpolants compute derivatives of pressure and energy and may enforce the stability conditions, but do not enforce the consistency condition and its derivatives. We describe a new type of table interface based on a constrained local least squares regression technique. It is applied to several SESAME EOS's showing how the consistency condition can be satisfied to round-off while computing first and second derivatives with demonstrated second-order convergence. An improvement of 14 orders of magnitude over conventional derivatives is demonstrated, although the new method is apparently two orders of magnitude slower, due to the fact that every evaluation requires solving an 11-dimensional nonlinear system.Comment: 29 pages, 9 figures, 16 references, submitted to Phys Rev

Existence and approximation of probability measure solutions to models of collective behaviors

In this paper we consider first order differential models of collective behaviors of groups of agents based on the mass conservation equation. Models are formulated taking the spatial distribution of the agents as the main unknown, expressed in terms of a probability measure evolving in time. We develop an existence and approximation theory of the solutions to such models and we show that some recently proposed models of crowd and swarm dynamics fit our theoretic paradigm.Comment: 31 pages, 1 figur

Efficient implementation of finite volume methods in Numerical Relativity

Centered finite volume methods are considered in the context of Numerical Relativity. A specific formulation is presented, in which third-order space accuracy is reached by using a piecewise-linear reconstruction. This formulation can be interpreted as an 'adaptive viscosity' modification of centered finite difference algorithms. These points are fully confirmed by 1D black-hole simulations. In the 3D case, evidence is found that the use of a conformal decomposition is a key ingredient for the robustness of black hole numerical codes.Comment: Revised version, 10 pages, 6 figures. To appear in Phys. Rev.

Raman signatures of classical and quantum phases in coupled dots: A theoretical prediction

We study electron molecules in realistic vertically coupled quantum dots in a strong magnetic field. Computing the energy spectrum, pair correlation functions, and dynamical form factor as a function of inter-dot coupling via diagonalization of the many-body Hamiltonian, we identify structural transitions between different phases, some of which do not have a classical counterpart. The calculated Raman cross section shows how such phases can be experimentally singled out.Comment: 9 pages, 2 postscript figures, 1 colour postscript figure, Latex 2e, Europhysics Letters style and epsfig macros. Submitted to Europhysics Letter

Finite difference lattice Boltzmann model with flux limiters for liquid-vapor systems

In this paper we apply a finite difference lattice Boltzmann model to study the phase separation in a two-dimensional liquid-vapor system. Spurious numerical effects in macroscopic equations are discussed and an appropriate numerical scheme involving flux limiter techniques is proposed to minimize them and guarantee a better numerical stability at very low viscosity. The phase separation kinetics is investigated and we find evidence of two different growth regimes depending on the value of the fluid viscosity as well as on the liquid-vapor ratio.Comment: 10 pages, 10 figures, to be published in Phys. Rev.

Are gauge shocks really shocks?

The existence of gauge pathologies associated with the Bona-Masso family of generalized harmonic slicing conditions is proven for the case of simple 1+1 relativity. It is shown that these gauge pathologies are true shocks in the sense that the characteristic lines associated with the propagation of the gauge cross, which implies that the name ``gauge shock'' usually given to such pathologies is indeed correct. These gauge shocks are associated with places where the spatial hypersurfaces that determine the foliation of spacetime become non-smooth.Comment: 7 pages, 5 figures, REVTEX 4. Revised version, including corrections suggested by referee