794 research outputs found
Polarization Modeling and Predictions for DKIST Part 2: Application of the Berreman Calculus to Spectral Polarization Fringes of Beamsplitters and Crystal Retarders
We outline polarization fringe predictions derived from a new application of
the Berreman calculus for the Daniel K. Inouye Solar Telescope (DKIST) retarder
optics. The DKIST retarder baseline design used 6 crystals, single-layer
anti-reflection coatings, thick cover windows and oil between all optical
interfaces. This new tool estimates polarization fringes and optic Mueller
matrices as functions of all optical design choices. The amplitude and period
of polarized fringes under design changes, manufacturing errors, tolerances and
several physical factors can now be estimated. This tool compares well with
observations of fringes for data collected with the SPINOR spectropolarimeter
at the Dunn Solar Telescope using bi-crystalline achromatic retarders as well
as laboratory tests. With this new tool, we show impacts of design decisions on
polarization fringes as impacted by anti-reflection coatings, oil refractive
indices, cover window presence and part thicknesses. This tool helped DKIST
decide to remove retarder cover windows and also recommends reconsideration of
coating strategies for DKIST. We anticipate this tool to be essential in
designing future retarders for mitigation of polarization and intensity fringe
errors in other high spectral resolution astronomical systems.Comment: Accepted for publication in JATI
Representations of world coordinates in FITS
The initial descriptions of the FITS format provided a simplified method for
describing the physical coordinate values of the image pixels, but deliberately
did not specify any of the detailed conventions required to convey the
complexities of actual image coordinates. Building on conventions in wide use
within astronomy, this paper proposes general extensions to the original
methods for describing the world coordinates of FITS data. In subsequent
papers, we apply these general conventions to the methods by which spherical
coordinates may be projected onto a two-dimensional plane and to
frequency/wavelength/velocity coordinates.Comment: 15 Pages, 1 figure, LaTex with Astronomy & Astrophysics macro
package, submitted to A&A, related papers at
http://www.aoc.nrao.edu/~egreise
An open and parallel multiresolution framework using block-based adaptive grids
A numerical approach for solving evolutionary partial differential equations
in two and three space dimensions on block-based adaptive grids is presented.
The numerical discretization is based on high-order, central finite-differences
and explicit time integration. Grid refinement and coarsening are triggered by
multiresolution analysis, i.e. thresholding of wavelet coefficients, which
allow controlling the precision of the adaptive approximation of the solution
with respect to uniform grid computations. The implementation of the scheme is
fully parallel using MPI with a hybrid data structure. Load balancing relies on
space filling curves techniques. Validation tests for 2D advection equations
allow to assess the precision and performance of the developed code.
Computations of the compressible Navier-Stokes equations for a temporally
developing 2D mixing layer illustrate the properties of the code for nonlinear
multi-scale problems. The code is open source
A rarefaction-tracking method for hyperbolic conservation laws
We present a numerical method for scalar conservation laws in one space
dimension. The solution is approximated by local similarity solutions. While
many commonly used approaches are based on shocks, the presented method uses
rarefaction and compression waves. The solution is represented by particles
that carry function values and move according to the method of characteristics.
Between two neighboring particles, an interpolation is defined by an analytical
similarity solution of the conservation law. An interaction of particles
represents a collision of characteristics. The resulting shock is resolved by
merging particles so that the total area under the function is conserved. The
method is variation diminishing, nevertheless, it has no numerical dissipation
away from shocks. Although shocks are not explicitly tracked, they can be
located accurately. We present numerical examples, and outline specific
applications and extensions of the approach.Comment: 21 pages, 7 figures. Similarity 2008 conference proceeding
The impact of organisational external peer review on colorectal cancer treatment and survival in the Netherlands
Background: Organisational external peer review was introduced in 1994 in the Netherlands to improve multidisciplinary cancer care. We examined the clinical impact of this programme on colorectal cancer care. Methods: Patients with primary colorectal cancer were included from 23 participating hospitals and 7 controls. Hospitals from the intervention group were dichotomised by their implementation proportion (IP) of the recommendations from each peer review (high IP vs low IP). Outcome measures were the introduction of new multidisciplinary therapies and survival. Results: In total, 45 705 patients were included (1990-2010). Patients from intervention hospitals more frequently received adjuvant chemotherapy for stage III colon cancer. T2-3/M0 rectal cancer patients from hospitals with a high IP had a higher chance of receiving preoperative radiotherapy (OR 1.31, 95% CI 1.11-1.55) compared with the controls and low IP group (OR 0.75, 95% CI 0.63-0.88). There were no differences in the use of preoperative chemoradiation for T4/M0 rectal cancer. Survival was slightly higher in colon cancer patients from intervention hospitals but unrelated to the phase of the programme in which the hospital was at the time of diagnosis. Conclusions: Some positive effects of external peer review on cancer care were found, but the results need to be interpreted cautiously due to the ambiguity of the outcomes and possible confounding factors
On the validity of mean-field amplitude equations for counterpropagating wavetrains
We rigorously establish the validity of the equations describing the
evolution of one-dimensional long wavelength modulations of counterpropagating
wavetrains for a hyperbolic model equation, namely the sine-Gordon equation. We
consider both periodic amplitude functions and localized wavepackets. For the
localized case, the wavetrains are completely decoupled at leading order, while
in the periodic case the amplitude equations take the form of mean-field
(nonlocal) Schr\"odinger equations rather than locally coupled partial
differential equations. The origin of this weakened coupling is traced to a
hidden translation symmetry in the linear problem, which is related to the
existence of a characteristic frame traveling at the group velocity of each
wavetrain. It is proved that solutions to the amplitude equations dominate the
dynamics of the governing equations on asymptotically long time scales. While
the details of the discussion are restricted to the class of model equations
having a leading cubic nonlinearity, the results strongly indicate that
mean-field evolution equations are generic for bimodal disturbances in
dispersive systems with \O(1) group velocity.Comment: 16 pages, uuencoded, tar-compressed Postscript fil
Kink instabilities in jets from rotating magnetic fields
We have performed 2.5D and 3D simulations of conical jets driven by the
rotation of an ordered, large-scale magnetic field in a stratified atmosphere.
The simulations cover about three orders of magnitude in distance to capture
the centrifugal acceleration as well as the evolution past the Alfven surface.
We find that the jets develop kink instabilities, the characteristics of which
depend on the velocity profile imposed at the base of the flow. The
instabilities are especially pronounced with a rigid rotation profile, which
induces a shearless magnetic field. The jet's expansion appears to be limiting
the growth of Alfven mode instabilities.Comment: 10 pages, 13 figures, accepted for publication in A&
Relaxation and reconstruction on (111) surfaces of Au, Pt, and Cu
We have theoretically studied the stability and reconstruction of (111)
surfaces of Au, Pt, and Cu. We have calculated the surface energy, surface
stress, interatomic force constants, and other relevant quantities by ab initio
electronic structure calculations using the density functional theory (DFT), in
a slab geometry with periodic boundary conditions. We have estimated the
stability towards a quasi-one-dimensional reconstruction by using the
calculated quantities as parameters in a one-dimensional Frenkel-Kontorova
model. On all surfaces we have found an intrinsic tensile stress. This stress
is large enough on Au and Pt surfaces to lead to a reconstruction in which a
denser surface layer is formed, in agreement with experiment. The
experimentally observed differences between the dense reconstruction pattern on
Au(111) and a sparse structure of stripes on Pt(111) are attributed to the
details of the interaction potential between the first layer of atoms and the
substrate.Comment: 8 pages, 3 figures, submitted to Physical Review
Phase Slips and the Eckhaus Instability
We consider the Ginzburg-Landau equation, , with complex amplitude . We first analyze the phenomenon of
phase slips as a consequence of the {\it local} shape of . We next prove a
{\it global} theorem about evolution from an Eckhaus unstable state, all the
way to the limiting stable finite state, for periodic perturbations of Eckhaus
unstable periodic initial data. Equipped with these results, we proceed to
prove the corresponding phenomena for the fourth order Swift-Hohenberg
equation, of which the Ginzburg-Landau equation is the amplitude approximation.
This sheds light on how one should deal with local and global aspects of phase
slips for this and many other similar systems.Comment: 22 pages, Postscript, A
Scalar field induced oscillations of neutron stars and gravitational collapse
We study the interaction of massless scalar fields with self-gravitating
neutron stars by means of fully dynamic numerical simulations of the
Einstein-Klein-Gordon perfect fluid system. Our investigation is restricted to
spherical symmetry and the neutron stars are approximated by relativistic
polytropes. Studying the nonlinear dynamics of isolated neutron stars is very
effectively performed within the characteristic formulation of general
relativity, in which the spacetime is foliated by a family of outgoing light
cones. We are able to compactify the entire spacetime on a computational grid
and simultaneously impose natural radiative boundary conditions and extract
accurate radiative signals. We study the transfer of energy from the scalar
field to the fluid star. We find, in particular, that depending on the
compactness of the neutron star model, the scalar wave forces the neutron star
either to oscillate in its radial modes of pulsation or to undergo
gravitational collapse to a black hole on a dynamical timescale. The radiative
signal, read off at future null infinity, shows quasi-normal oscillations
before the setting of a late time power-law tail.Comment: 12 pages, 13 figures, submitted to Phys. Rev.
- …