Fast nucleon emission as a probe of the isospin momentum dependence

In this article we investigate the structure of the non-local part of the symmetry term, that leads to a splitting of the effective masses of protons and neutrons in asymmetric matter. Based on microscopic transport simulations we suggest some rather sensitive observables in collisions of neutron-rich (unstable) ions at intermediate ($RIA$) energies. In particular we focus the attention on pre-equilibrium nucleon emissions. We discuss interesting correlations between the N/Z content of the fast emitted particles and their rapidity or transverse momentum, that show a nice dependence on the prescription used for the effective mass splitting.Comment: 5 pages, 6 figures, revtex

Robustness of a high-resolution central scheme for hydrodynamic simulations in full general relativity

A recent paper by Lucas-Serrano et al. indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong shock waves. In this paper we apply this scheme in full general relativity (involving {\it dynamical} spacetimes), and assess its suitability by performing test simulations for oscillations of rapidly rotating neutron stars and merger of binary neutron stars. It is demonstrated that this HRC scheme can yield results as accurate as those by the so-called high-resolution shock-capturing (HRSC) schemes based upon Riemann solvers. Furthermore, the adopted HRC scheme has increased computational efficiency as it avoids the costly solution of Riemann problems and has practical advantages in the modeling of neutron star spacetimes. Namely, it allows simulations with stiff equations of state by successfully dealing with very low-density unphysical atmospheres. These facts not only suggest that such a HRC scheme may be a desirable tool for hydrodynamical simulations in general relativity, but also open the possibility to perform accurate magnetohydrodynamical simulations in curved dynamic spacetimes.Comment: 4 pages, to be published in Phys. Rev. D (brief report

Keck Echellette Spectrograph and Imager Observations of Metal-poor Damped Lyα Systems

We present the first results from a survey of SDSS quasars selected for strong H I damped Lyα (DLA) absorption with corresponding low equivalent width absorption from strong low-ion transitions (e.g., C II λ1334 and Si II λ1260). These metal-poor DLA candidates were selected from the SDSS fifth release quasar spectroscopic database, and comprise a large new sample for probing low-metallicity galaxies. Medium-resolution echellette spectra from the Keck Echellette Spectrograph and Imager spectrograph for an initial sample of 35 systems were obtained to explore the metal-poor tail of the DLA distribution and to investigate the nucleosynthetic patterns at these metallicities. We have estimated saturation corrections for the moderately underresolved spectra, and systems with very narrow Doppler parameters (b ≤ 5 km s^(–1)) will likely have underestimated abundances. For those systems with Doppler parameters b > 5 km s^(–1), we have measured low-metallicity DLA gas with [X/H] < –2.4 for at least one of C, O, Si, or Fe. Assuming non-saturated components, we estimate that several DLA systems have [X/H] < –2.8, including five DLA systems with both low equivalent widths and low metallicity in transitions of both C II and O I. All of the measured DLA metallicities, however, exceed or are consistent with a metallicity of at least 1/1000 of solar, regardless of the effects of saturation in our spectra. Our results indicate that the metal-poor tail of galaxies at z ~ 3 drops exponentially at [X/H] ≾ –3. If the distribution of metallicity is Gaussian, the probability of identifying interstellar medium gas with lower abundance is extremely small, and our results suggest that DLA systems with [X/H] < –4.0 are extremely rare, and could comprise only 8 × 10^(–7) of DLA systems. The relative abundances of species within these low-metallicity DLA systems are compared with stellar nucleosynthesis models, and are consistent with stars having masses of 30 M_⊙ < M * < 100 M_⊙. The observed ratio of [C/O] for values of [O/H] < –2.5 exceeds values seen in moderate metallicity DLA systems, and also exceeds theoretical nucleosynthesis predictions for higher mass Population III stars. We also have observed a correlation between the column density N(C IV) with [Si/H] metallicity, suggestive of a trend between mass of the DLA system and its metallicity

Consideration of the Mechanisms for Tidal Bore Formation in an Idealized Planform Geometry

A tidal bore is a positive wave traveling upstream along the estuary of a river, generated by a relatively rapid rise of the tide, often enhanced by the funneling shape of the estuary. The swell produced by the tide grows and its front steepens as the flooding tide advances inland, promoting the formation of a sharp front wave, i.e., the tidal bore. Because of the many mechanisms and conditions involved in the process, it is difficult to formulate an effective criterion to predict the bore formation. In this preliminary analysis, aimed at bringing out the main processes and parameters that control tidal bore formation, the degrees of freedom of the problem are largely reduced by considering a rectangular channel of constant width with uniform flow, forced downstream by rising the water level at a constant rate. The framework used in this study is extremely simple, yet the problem is still complex and the solution is far from being trivial. From the results of numerical simulations, three distinctive behaviors emerged related to conditions in which a tidal bore forms, a tidal bore does not form, and a weak bore forms; the latter has a weakly steep front and after the bore formed it rapidly vanishes. Based on these behaviors, some criteria to predict the bore formation are proposed and discussed. The more effective criterion, suitably rearranged, is checked against data from real estuaries and the predictions are found to compare favorably with the available data

Quiet Sun magnetic fields from space-borne observations: simulating Hinode's case

We examine whether or not it is possible to derive the field strength distribution of quiet Sun internetwork regions from very high spatial resolution polarimetric observations in the visible. In particular, we consider the case of the spectropolarimeter attached to the Solar Optical Telescope aboard Hinode. Radiative magneto-convection simulations are used to synthesize the four Stokes profiles of the \ion{Fe}{1} 630.2 nm lines. Once the profiles are degraded to a spatial resolution of 0\farcs32 and added noise, we infer the atmospheric parameters by means of Milne-Eddington inversions. The comparison of the derived values with the real ones indicates that the visible lines yield correct internetwork field strengths and magnetic fluxes, with uncertainties smaller than $\sim$150 G, when a stray light contamination factor is included in the inversion. Contrary to the results of ground-based observations at 1\arcsec, weak fields are retrieved wherever the field is weak in the simulation.Comment: Accepted for publication in ApJ Letter

Effect of dynamic stall on the aerodynamics of vertical-axis wind turbines

Accurate simulations of the aerodynamic performance of vertical-axis wind turbines pose a significant challenge for computational fluid dynamics methods. The aerodynamic interaction between the blades of the rotor and the wake that is produced by the blades requires a high-fidelity representation of the convection of vorticity within the wake. In addition, the cyclic motion of the blades induces large variations in the angle of attack on the blades that can manifest as dynamic stall. The present paper describes the application of a numerical model that is based on the vorticity transport formulation of the Navier–Stokes equations, to the prediction of the aerodynamics of a verticalaxis wind turbine that consists of three curved rotor blades that are twisted helically around the rotational axis of the rotor. The predicted variation of the power coefficient with tip speed ratio compares very favorably with experimental measurements. It is demonstrated that helical blade twist reduces the oscillation of the power coefficient that is an inherent feature of turbines with non-twisted blade configurations

Binary neutron-star mergers with Whisky and SACRA: First quantitative comparison of results from independent general-relativistic hydrodynamics codes

We present the first quantitative comparison of two independent general-relativistic hydrodynamics codes, the Whisky code and the SACRA code. We compare the output of simulations starting from the same initial data and carried out with the configuration (numerical methods, grid setup, resolution, gauges) which for each code has been found to give consistent and sufficiently accurate results, in particular in terms of cleanness of gravitational waveforms. We focus on the quantities that should be conserved during the evolution (rest mass, total mass energy, and total angular momentum) and on the gravitational-wave amplitude and frequency. We find that the results produced by the two codes agree at a reasonable level, with variations in the different quantities but always at better than about 10%.Comment: Published on Phys. Rev.

Three-dimensional Simulations of Accretion to Stars with Complex Magnetic Fields

Disk accretion to rotating stars with complex magnetic fields is investigated using full three-dimensional magnetohydrodynamic (MHD) simulations. The studied magnetic configurations include superpositions of misaligned dipole and quadrupole fields and off-centre dipoles. The simulations show that when the quadrupole component is comparable to the dipole component, the magnetic field has a complex structure with three major magnetic poles on the surface of the star and three sets of loops of field lines connecting them. A significant amount of matter flows to the quadrupole "belt", forming a ring-like hot spot on the star. If the maximum strength of the magnetic field on the star is fixed, then we observe that the mass accretion rate, the torque on the star, and the area covered by hot spots are several times smaller in the quadrupole-dominant cases than in the pure dipole cases. The influence of the quadrupole component on the shape of the hot spots becomes noticeable when the ratio of the quadrupole and dipole field strengths $B_q/B_d\gtrsim0.5$, and becomes dominant when $B_q/B_d\gtrsim1$. In the case of an off-centre dipole field, most of the matter flows through a one-armed accretion stream, forming a large hot spot on the surface, with a second much smaller secondary spot. The light curves may have simple, sinusoidal shapes, thus mimicking stars with pure dipole fields. Or, they may be complex and unusual. In some cases the light curves may be indicators of a complex field, in particular if the inclination angle is known independently. We also note that in the case of complex fields, magnetospheric gaps are often not empty, and this may be important for the survival of close-in exosolar planets.Comment: 13 pages, 21 figures, accepted for publication in MNRA

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