Cauchyperturbative matching reexamined: Tests in spherical symmetry
Abstract
During the last few years progress has been made on several fronts making it possible to revisit Cauchyperturbative matching (CPM) in numerical relativity in a more robust and accurate way. This paper is the first in a series where we plan to analyze CPM in the light of these new results. One of the new developments is an understanding of how to impose constraintpreserving boundary conditions (CPBC); though most of the related research has been driven by outer boundaries, one can use them for matching interface boundaries as well. Another front is related to numerically stable evolutions using multiple patches, which in the context of CPM allows the matching to be performed on a spherical surface, thus avoiding interpolations between Cartesian and spherical grids. One way of achieving stability for such schemes of arbitrary high order is through the use of penalty techniques and discrete derivatives satisfying summation by parts (SBP). Recently, new, very efficient and highorder accurate derivatives satisfying SBP and associated dissipation operators have been constructed. Here we start by testing all these techniques applied to CPM in a setting that is simple enough to study all the ingredients in great detail: Einstein's equations in spherical symmetry, describingmore »
 Authors:
 MaxPlanckInstitut fuer Astrophysik, KarlSchwarzschildStr. 1, 85741 Garching (Germany)
 Department of Physics and Astronomy, 202 Nicholson Hall, Louisana State University, Baton Rouge, Louisiana 70803 (United States)
 (United States)
 (Guatemala)
 Publication Date:
 OSTI Identifier:
 20782910
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Physical Review. D, Particles Fields; Journal Volume: 73; Journal Issue: 8; Other Information: DOI: 10.1103/PhysRevD.73.084011; (c) 2006 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; BLACK HOLES; BOUNDARY CONDITIONS; COSMIC GAMMA BURSTS; COSMOLOGY; EINSTEIN FIELD EQUATIONS; INTERPOLATION; SCALAR FIELDS; SCHWARZSCHILD METRIC; SPATIAL RESOLUTION; SPHERICAL CONFIGURATION; SYMMETRY
Citation Formats
Zink, Burkhard, Pazos, Enrique, Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803, Departamento de Matematica, Universidad de San Carlos de Guatemala, Edificio T4, Facultad de Ingenieria, Ciudad Universitaria Z. 12, Diener, Peter, Tiglio, Manuel, and Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803. Cauchyperturbative matching reexamined: Tests in spherical symmetry. United States: N. p., 2006.
Web. doi:10.1103/PHYSREVD.73.084011.
Zink, Burkhard, Pazos, Enrique, Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803, Departamento de Matematica, Universidad de San Carlos de Guatemala, Edificio T4, Facultad de Ingenieria, Ciudad Universitaria Z. 12, Diener, Peter, Tiglio, Manuel, & Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803. Cauchyperturbative matching reexamined: Tests in spherical symmetry. United States. doi:10.1103/PHYSREVD.73.084011.
Zink, Burkhard, Pazos, Enrique, Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803, Departamento de Matematica, Universidad de San Carlos de Guatemala, Edificio T4, Facultad de Ingenieria, Ciudad Universitaria Z. 12, Diener, Peter, Tiglio, Manuel, and Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803. Sat .
"Cauchyperturbative matching reexamined: Tests in spherical symmetry". United States.
doi:10.1103/PHYSREVD.73.084011.
@article{osti_20782910,
title = {Cauchyperturbative matching reexamined: Tests in spherical symmetry},
author = {Zink, Burkhard and Pazos, Enrique and Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803 and Departamento de Matematica, Universidad de San Carlos de Guatemala, Edificio T4, Facultad de Ingenieria, Ciudad Universitaria Z. 12 and Diener, Peter and Tiglio, Manuel and Center for Computation and Technology, 302 Johnston Hall, Louisana State University, Baton Rouge, Louisiana 70803},
abstractNote = {During the last few years progress has been made on several fronts making it possible to revisit Cauchyperturbative matching (CPM) in numerical relativity in a more robust and accurate way. This paper is the first in a series where we plan to analyze CPM in the light of these new results. One of the new developments is an understanding of how to impose constraintpreserving boundary conditions (CPBC); though most of the related research has been driven by outer boundaries, one can use them for matching interface boundaries as well. Another front is related to numerically stable evolutions using multiple patches, which in the context of CPM allows the matching to be performed on a spherical surface, thus avoiding interpolations between Cartesian and spherical grids. One way of achieving stability for such schemes of arbitrary high order is through the use of penalty techniques and discrete derivatives satisfying summation by parts (SBP). Recently, new, very efficient and highorder accurate derivatives satisfying SBP and associated dissipation operators have been constructed. Here we start by testing all these techniques applied to CPM in a setting that is simple enough to study all the ingredients in great detail: Einstein's equations in spherical symmetry, describing a black hole coupled to a massless scalar field. We show that with the techniques described above, the errors introduced by Cauchyperturbative matching are very small, and that very longterm and accurate CPM evolutions can be achieved. Our tests include the accretion and ringdown phase of a Schwarzschild black hole with CPM, where we find that the discrete evolution introduces, with a low spatial resolution of {delta}r=M/10, an error of 0.3% after an evolution time of 1,000,000M. For a black hole of solar mass, this corresponds to approximately 5s, and is therefore at the lower end of timescales discussed e.g. in the collapsar model of gammaray burst engines.},
doi = {10.1103/PHYSREVD.73.084011},
journal = {Physical Review. D, Particles Fields},
number = 8,
volume = 73,
place = {United States},
year = {Sat Apr 15 00:00:00 EDT 2006},
month = {Sat Apr 15 00:00:00 EDT 2006}
}

Matching conditions relate couplings to particle masses. We discuss the importance of oneloop matching conditions in the Higgs boson and top quark sector as well as the choice of the matching scale. We argue for matching scales {mu}{sub 0,t}{approx_equal}m{sub t} and {mu}{sub 0,H}{approx_equal}max{l_brace}m{sub t},M{sub H}{r_brace}. Using these results, the twoloop Higgs boson mass upper bounds are reanalyzed. Previous results for {Lambda}{approx} few TeV are found to be too stringent, and a recent update is found to be wrong. For {Lambda}=10{sup 19} GeV we find M{sub H}{lt}180{plus_minus}4{plus_minus}5 GeV, the first error indicating the theoretical uncertainty and the second error reflecting themore »

A NEW MULTIDIMENSIONAL GENERAL RELATIVISTIC NEUTRINO HYDRODYNAMIC CODE FOR CORECOLLAPSE SUPERNOVAE. I. METHOD AND CODE TESTS IN SPHERICAL SYMMETRY
We present a new general relativistic code for hydrodynamical supernova simulations with neutrino transport in spherical and azimuthal symmetry (one dimension and two dimensions, respectively). The code is a combination of the COCONUT hydro module, which is a Riemannsolverbased, highresolution shockcapturing method, and the threeflavor, fully energydependent VERTEX scheme for the transport of massless neutrinos. VERTEX integrates the coupled neutrino energy and momentum equations with a variable Eddington factor closure computed from a model Boltzmann equation and uses the 'raybyray plus' approximation in two dimensions, assuming the neutrino distribution to be axially symmetric around the radial direction at every pointmore » 
Ensembleaverage versus suspensionscale Cauchy continuummechanical definitions of stress in polarized suspensions: Global homogenization of a dilute suspension of dipolar spherical particles
The macroscale rheological properties of a dilute suspension exposed to a uniform external field and composed of identical, rigid, inhomogeneous, dipolar, spherical particles dispersed in an incompressible Newtonian fluid and possessing the same mean density as the latter fluid are derived from knowledge of its microscale properties by applying a global ensembleaveraging technique. Each dipole, which is permanently embedded in the particle, is assumed to be generated by the presence of an inhomogeneous external bodyforce field in the particle interior resulting from the action of the uniform external field on an inhomogeneous distribution of interior matter. It is shown thatmore » 
Generalized density matrix reexamined: Microscopic approach to collective dynamics in soft spherical nuclei
The generalized density matrix method is used to calculate microscopically the parameters of the collective Hamiltonian. Higherorder anharmonicities are obtained consistently with the lowestorder results, the mean field (HartreeFockBogoliubov equation), and the harmonic potential (quasiparticle randomphase approximation). The method is applied to soft spherical nuclei, where the anharmonicities are essential for restoring the stability of the system, as the harmonic potential becomes small or negative. The approach is tested in three models of increasing complexity: the Lipkin model, model with factorizable forces, and the quadrupole plus pairing model.