N-body simulations using customized potential–density pair basis sets

Potential–density pair basis sets can be used for highly efficient N -body simulation codes, but they suffer from a lack of versatility, i.e. a basis set has to be constructed for each different class of stellar system. We present numerical techniques for generating a biorthonormal potential–density pair basis set that has a general specified pair as its lowest-order member. We go on to demonstrate how the set can be used to construct N -body equilibria, which we then evolve using an N -body code that calculates forces using the basis set.     

Magnetic field bending in accretion discs with dynamos

We consider the problem of poloidal magnetic field advection and bending of an initially vertical field owing to radial inflow in thin accretion discs. For a ratio of kinematic viscosity to magnetic diffusivity of order unity, significant bending of an externally applied vertical field cannot occur in a disc with no internal dynamo. However, we show that if poloidal field is generated by a dynamo operating near its critical state, then significant field bending may be possible. Our results are of particular relevance to wind launching from accretion discs.     

An efficient methodology for simulating earthquake‐induced 3D pounding of buildings

The current paper presents an efficient methodology for numerically simulating in three dimensions adjacent buildings that may experience pounding during strong earthquakes. In particular, a new approach to the numerical problem of spatial impact modeling that does not require the ‘a priori’ determination of the contact points is presented, taking also into account the geometry at the vicinity of an impact. In the current study, the buildings are simulated as linear multi‐degree‐of‐freedom‐systems, but the methodology can be easily extended to consider nonlinear behavior as well. A software application has been specifically developed to implement the proposed methodology, using modern object‐oriented design and programming. The developed software is utilized in a simple example, and the computed results are compared with the corresponding analysis results obtained from a commercial general‐purpose software application that uses typical contact elements for the simulation of impacts. A discussion follows on the advantages and capabilities of the proposed methodology and the developed software. Copyright © 2013 John Wiley & Sons, Ltd.     

Consequences of tidal interaction between disks and orbiting protoplanets for the evolution of multi-planet systems with architecture resembling that of Kepler 444

We study orbital evolution of multi-planet systems with masses in the terrestrial planet regime induced through tidal interaction with a protoplanetary disk assuming that this is the dominant mechanism for producing orbital migration and circularization. We develop a simple analytic model for a system that maintains consecutive pairs in resonance while undergoing orbital circularization and migration. This model enables migration times for each planet to be estimated once planet masses, circularization times and the migration time for the innermost planet are specified. We applied it to a system with the current architecture of Kepler 444 adopting a simple protoplanetary disk model and planet masses that yield migration times inversely proportional to the planet mass, as expected if they result from torques due to tidal interaction with the protoplanetary disk. Furthermore the evolution time for the system as a whole is comparable to current protoplanetary disk lifetimes. In addition we have performed a number of numerical simulations with input data obtained from this model. These indicate that although the analytic model is inexact, relatively small corrections to the estimated migration rates yield systems for which period ratios vary by a minimal extent. Because of relatively large deviations from exact resonance in the observed system of up to 2 %, the migration times obtained in this way indicate only weak convergent migration such that a system for which the planets did not interact would contract by only \({\sim }1\,\%\) although undergoing significant inward migration as a whole. We have also performed additional simulations to investigate conditions under which the system could undergo significant convergent migration before reaching its final state. These indicate that migration times have to be significantly shorter and resonances between planet pairs significantly closer during such an evolutionary phase. Relative migration rates would then have to decrease allowing period ratios to increase to become more distant from resonances as the system approached its final state in the inner regions of the protoplanetary disk.     

The convective stability of a thin viscous disc

Abstract

We prove that the presence of viscosity does not affect stability to axisymmetric convective modes of a thin differentially rotating disc with no thermal conduction but in which viscosity is taken fully into account. In such a case the Schwarzschild criterion is necessary and sufficient for convective stability to local perturbations. In the proof we use a general formulation of local stability analysis, which allows a rigorous demonstration. Restricted particular forms of the viscous stress tensor introduced in the modelling of thin accretion discs may lead to viscous overstabilities. The additional instability found by Elstner et al. (1989) and described by the authors as a correction to the Schwarzschild criterion is a manifestation of these. However, when viscosity is taken fully into account, such instabilities cannot be discussed within the framework of a local analysis, a fully global treatment being required.