Multi-Level quasi-Newton methods for the partitioned simulation of fluid-structure interaction

In previous work of the authors, Fourier stability analyses have been performed of Gauss-Seidel iterations between the flow solver and the structural solver in a partitioned fluid-structure interaction simulation. These analyses of the flow in an elastic tube demonstrated that only a number of Fourier modes in the error on the interface displacement are unstable. Moreover, the modes with a low wave number are most unstable and these modes can be resolved on a coarser grid. Therefore, a new class of quasi-Newton methods with more than one grid level is introduced. Numerical experiments show a significant reduction in run time

Patient-specific CFD simulation of intraventricular haemodynamics based on 3D ultrasound imaging

Background: The goal of this paper is to present a computational fluid dynamic (CFD) model with moving boundaries to study the intraventricular flows in a patient-specific framework. Starting from the segmentation of real-time transesophageal echocardiographic images, a CFD model including the complete left ventricle and the moving 3D mitral valve was realized. Their motion, known as a function of time from the segmented ultrasound images, was imposed as a boundary condition in an Arbitrary Lagrangian-Eulerian framework. Results: The model allowed for a realistic description of the displacement of the structures of interest and for an effective analysis of the intraventricular flows throughout the cardiac cycle. The model provides detailed intraventricular flow features, and highlights the importance of the 3D valve apparatus for the vortex dynamics and apical flow. Conclusions: The proposed method could describe the haemodynamics of the left ventricle during the cardiac cycle. The methodology might therefore be of particular importance in patient treatment planning to assess the impact of mitral valve treatment on intraventricular flow dynamics

The Potential of the Timing Method to Detect Evolved Planetary Systems

The timing method, using either stellar pulsations or eclipse timing of close binaries as a clock, is proving to be an efficient way to detect planets around stars that have evolved beyond the red giant branch. In this article we present a short review of the recent discoveries and we investigate the potential of the timing method using data both from ground-based facilities as well as from the Kepler and CoRoT space missions.Comment: Part of PlanetsbeyondMS/2010 proceedings http://arxiv.org/html/1011.6606v1, Proc. of the workshop on "Planetary Systems beyond the Main Sequence" (Bamberg, 11-14 August 2010), AIPC in press (eds. S. Schuh, H. Drechsel and U. Heber), 15 pages, 5 figure

The blue-edge problem of the V1093 Her instability strip revisited using evolutionary models with atomic diffusion

We have computed a new grid of evolutionary subdwarf B star (sdB) models from the start of central He burning, taking into account atomic diffusion due to radiative levitation, gravitational settling, concentration diffusion, and thermal diffusion. We have computed the non-adiabatic pulsation properties of the models and present the predicted p-mode and g-mode instability strips. In previous studies of the sdB instability strips, artificial abundance enhancements of Fe and Ni were introduced in the pulsation driving layers. In our models, the abundance enhancements of Fe and Ni occur naturally, eradicating the need to use artificial enhancements. We find that the abundance increases of Fe and Ni were previously underestimated and show that the instability strip predicted by our simulations solves the so-called blue edge problem of the subdwarf B star g-mode instability strip. The hottest known g-mode pulsator, KIC 10139564, now resides well within the instability strip {even when only modes with low spherical degrees (l<=2) are considered.Comment: 7 pages, 7 figures. Accepted for publication in Astronomy & Astrophysic

IRAS 19135+3937: An SRd variable as interacting binary surrounded by a circumbinary disc

Semi-regular (SR) variables are not a homogeneous class and their variability is often explained due to pulsations and/or binarity. This study focuses on IRAS 19135+3937, an SRd variable with an infra-red excess indicative of a dusty disc. A time-series of high-resolution spectra, UBV photometry as well as a very accurate light curve obtained by the Kepler satellite, allowed us to study the object in unprecedented detail. We discovered it to be a binary with a period of 127 days. The primary has a low surface gravity and an atmosphere depleted in refractory elements. This combination of properties unambiguously places IRAS 19135+3937 in the subclass of post-Asymptotic Giant Branch stars with dusty discs. We show that the light variations in this object can not be due to pulsations, but are likely caused by the obscuration of the primary by the circumbinary disc during orbital motion. Furthermore, we argue that the double-peaked Fe emission lines provide evidence for the existence of a gaseous circumbinary Keplerian disc inside the dusty disc. A secondary set of absorption lines has been detected near light minimum, which we attribute to the reflected spectrum of the primary on the disc wall, which segregates due to the different Doppler shift. This corroborates the recent finding that reflection in the optical by this type of discs is very efficient. The system also shows a variable Halpha profile indicating a collimated outflow originating around the companion. IRAS 19135+3937 thus encompasses all the major emergent trends about evolved disc systems, that will eventually help to place these objects in the evolutionary context.Comment: Accepted to MNRA

Testing the asymptotic relation for period spacings from mixed modes of red giants observed with the Kepler mission

Dipole mixed pulsation modes of consecutive radial order have been detected for thousands of low-mass red-giant stars with the NASA space telescope Kepler. Such modes have the potential to reveal information on the physics of the deep stellar interior. Different methods have been proposed to derive an observed value for the gravity-mode period spacing, the most prominent one relying on a relation derived from asymptotic pulsation theory applied to the gravity-mode character of the mixed modes. Our aim is to compare results based on this asymptotic relation with those derived from an empirical approach for three pulsating red-giant stars. We developed a data-driven method to perform frequency extraction and mode identification. Next, we used the identified dipole mixed modes to determine the gravity-mode period spacing by means of an empirical method and by means of the asymptotic relation. In our methodology, we consider the phase offset, $\epsilon_{\mathrm{g}}$, of the asymptotic relation as a free parameter. Using the frequencies of the identified dipole mixed modes for each star in the sample, we derived a value for the gravity-mode period spacing using the two different methods. These differ by less than 5%. The average precision we achieved for the period spacing derived from the asymptotic relation is better than 1%, while that of our data-driven approach is 3%. Good agreement is found between values for the period spacing derived from the asymptotic relation and from the empirical method. Full abstract in PDF file.Comment: 14 pages, 13 figures, accepted for publication in A&

Novel methods for spatial prediction of soil functions within landscapes (SP0531)

Previous studies showed that soil patterns could be predicted in agriculturally managed landscapes by modelling and extrapolating from extensive existing but related integrated datasets. Based on these results we proposed to develop and apply predictive models of the relationships between environmental data and known soil patterns to predict capacity for key soil functions within diverse landscapes for which there is little detailed underpinning soil information available. Objectives were: To develop a high-level framework in which the non-specialist user-community could explore questions. To generate digital soil maps for three selected catchments at a target resolution of 1:50000 to provide the base information for soil function prediction. To use a modelling approach to predict the performance of key soil functions in catchments undergoing change but where only sparse or low resolution soil survey data are available. To use a modelling approach to assess the impact of different management scenarios and/or environmental conditions on the delivery of multiple soil functions within a catchment. To create a detailed outline of the requirements for ground-truthing to test the predicted model outputs at a catchment scale. To contribute to the development of a high-level framework for decision makers

Phoebe 2.0 – Triple and multiple systems

Some close binary formation theories require the presence of a third body so that the binary orbit can shrink over time. Tidal friction and Kozai cycles transfer energy from the binary to its companion, resulting in a close inner binary and a wide third body orbit. Spectroscopy and imaging studies have found 40% of binaries with periods less than 10 days, and 96% with periods less than 3 days, have a wide tertiary companion. With recent advancements in large photometric surveys, we are now beginning to detect many of these triple systems by observing tertiary eclipses or through the effect they have on the eclipse timing variations (ETVs) of the inner-binary. In the sample of 2600 Kepler EBs, we have detected the possible presence of a third body in ∼20%, including several circumbinary planets. Some multiple systems are quite dynamical and feature disappearing and reappearing eclipses, apsidal motion, and large disruptions to the inner-binary. phoebe is a freely available binary modeling code which can dynamically model all of these systems, allowing us to better test formation theories and probe the physics of eclipsing binaries

Mass ratio from Doppler beaming and R{\o}mer delay versus ellipsoidal modulation in the Kepler data of KOI-74

We present a light curve analysis and radial velocity study of KOI-74, an eclipsing A star + white dwarf binary with a 5.2 day orbit. Aside from new spectroscopy covering the orbit of the system, we used 212 days of publicly available Kepler observations and present the first complete light curve fitting to these data, modelling the eclipses and transits, ellipsoidal modulation, reflection, and Doppler beaming. Markov Chain Monte Carlo simulations are used to determine the system parameters and uncertainty estimates. Our results are in agreement with earlier studies, except that we find an inclination of 87.0 \pm 0.4\degree, which is significantly lower than the previously published value. We find that the mass ratio derived from the radial velocity amplitude (q=0.104 \pm 0.004) disagrees with that derived from the ellipsoidal modulation (q=0.052 \pm 0.004} assuming corotation). This was found before, but with our smaller inclination, the discrepancy is even larger than previously reported. Accounting for the rapid rotation of the A-star is found to increase the discrepancy even further by lowering the mass ratio to q=0.047 \pm 0.004. These results indicate that one has to be extremely careful in using the amplitude of an ellipsoidal modulation signal in a close binary to determine the mass ratio, when a proof of corotation is not firmly established. The radial velocities that can be inferred from the detected Doppler beaming in the light curve are found to be in agreement with our spectroscopic radial velocity determination. We also report the first measurement of R{\o}mer delay in a light curve of a compact binary. This delay amounts to -56 \pm 17 s and is consistent with the mass ratio derived from the radial velocity amplitude. The firm establishment of this mass ratio at q=0.104 \pm 0.004 leaves little doubt that the companion of KOI-74 is a low mass white dwarf.Comment: 9 pages, 7 figures, 2 tables; accepted for publication in MNRA