2,927 research outputs found
Critical Collapse of an Ultrarelativistic Fluid in the Limit
In this paper we investigate the critical collapse of an ultrarelativistic
perfect fluid with the equation of state in the limit of
. 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 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
- …