BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

John C. B. Papaloizou

John C. B. Papaloizou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-phastro-ph.SRastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2019arXiv

Low-degree mixed modes in red giant stars with moderate core magnetic fields

Observations of pressure-gravity mixed modes, combined with a theoretical framework for understanding mode formation, can yield a wealth of information about deep stellar interiors. In this paper, we seek to develop a formalism for treating the effects of deeply buried core magnetic fields on mixed modes in evolved stars, where the fields are moderate, i.e. not strong enough to disrupt wave propagation, but where they may be too strong for non-degenerate first-order perturbation theory to be applied. The magnetic field is incorporated in a way that avoids having to use this. Inclusion of the Lorentz force term is shown to yield a system of differential equations that allows for the magnetically-affected eigenfunctions to be computed from scratch, rather than following the approach of first-order perturbation theory. For sufficiently weak fields, coupling between different spherical harmonics can be neglected, allowing for reduction to a second-order system of ordinary differential equations akin to the usual oscillation equations that can be solved analogously. We derive expressions for (i) the mixed-mode quantisation condition in the presence of a field and (ii) the frequency shift associated with the magnetic field. In addition, for modes of low degree we uncover an extra offset term in the quantisation condition that is sensitive to properties of the evanescent zone. These expressions may be inverted to extract information about the stellar structure and magnetic field from observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the energy dissipation rate at the inner edge of circumbinary discs

We study, by means of numerical simulations and analysis, the details of the accretion process from a disc onto a binary system. We show that energy is dissipated at the edge of a circumbinary disc and this is associated with the tidal torque that maintains the cavity: angular momentum is transferred from the binary to the disc through the action of compressional shocks and viscous friction. These shocks can be viewed as being produced by fluid elements which drift into the cavity and, before being accreted, are accelerated onto trajectories that send them back to impact the disc. The rate of energy dissipation is approximately equal to the product of potential energy per unit mass at the disc's inner edge and the accretion rate, estimated from the disc parameters just beyond the cavity edge, that would occur without the binary. For very thin discs, the actual accretion rate onto the binary may be significantly less. We calculate the energy emitted by a circumbinary disc taking into account energy dissipation at the inner edge and also irradiation arising there from reprocessing of light from the stars. We find that, for tight PMS binaries, the SED is dominated by emission from the inner edge at wavelengths between 1-4 and 10 $μ$m. This may apply to systems like CoRoT 223992193 and V1481 Ori.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Three body resonances in close orbiting planetary systems: Tidal dissipation and orbital evolution

We study the orbital evolution of a three planet system with masses in the super-Earth regime resulting from the action of tides on the planets induced by the central star which cause orbital circularization. We consider systems either in or near to a three body commensurability for which adjacent pairs of planets are in a first order commensurability. We develop a simple analytic solution, derived from a time averaged set of equations, that describes the expansion of the system away from strict commensurability as a function of time, once a state where relevant resonant angles undergo small amplitude librations has been attained. We perform numerical simulations that show the attainment of such resonant states focusing on the Kepler 60 system. The results of the simulations confirm many of the scalings predicted by the appropriate analytic solution. We go on to indicate how the results can be applied to put constraints on the amount of tidal dissipation that has occurred in the system. For example, if the system has been in a librating state since its formation, we find that its present period ratios imply an upper limit on the time average of 1/Q', with Q' being the tidal dissipation parameter. On the other hand if a librating state has not been attained, a lower upper bound applies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Disk-planets interactions and the diversity of period ratios in Kepler's multi-planetary systems

The Kepler mission is dramatically increasing the number of planets known in multi-planetary systems. Many adjacent planets have orbital period ratios near resonant values, with a tendency to be larger than required for exact first-order mean-motion resonances. This intriguing feature has been shown to be a natural outcome of orbital circularization of resonant planetary pairs due to star-planet tidal interactions. However, this feature holds in multi-planetary systems with periods longer than ten days, for which tidal circularization is unlikely to provide efficient divergent evolution of the planets orbits. Gravitational interactions between planets and their parent protoplanetary disk may instead provide efficient divergent evolution. For a planet pair embedded in a disk, we show that interactions between a planet and the wake of its companion can reverse convergent migration, and significantly increase the period ratio from a near-resonant value. Divergent evolution due to wake-planet interactions is particularly efficient when at least one of the planets opens a partial gap around its orbit. This mechanism could help account for the diversity of period ratios in Kepler's multiple systems comprising super-Earth to sub-jovian planets with periods greater than about ten days. Diversity is also expected for planet pairs massive enough to merge their gap. The efficiency of wake-planet interactions is then much reduced, but convergent migration may stall with a variety of period ratios depending on the density structure in the common gap. This is illustrated for the Kepler-46 system, for which we reproduce the period ratio of Kepler-46b and c.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interaction between massive planets on inclined orbits and circumstellar discs

We study the interaction between massive planets and a gas disc with a mass in the range expected for protoplanetary discs. We use SPH simulations to study the orbital evolution of a massive planet as well as the dynamical response of the disc for planet masses between 1 and $6\ \rmn{M_J}$ and the full range of initial relative orbital inclinations. Gap formation can occur for planets in inclined orbits. For given planet mass, a threshold relative orbital inclination exists under which a gap forms. At high relative inclinations, the inclination decay rate increases for increasing planet mass and decreasing initial relative inclination. For an initial semi-major axis of 5 AU and relative inclination of $i_0=80^\circ,$ the times required for the inclination to decay by $10^\circ$ is $\sim10^{6}\ \rmn{yr}$ and $\sim10^{5}\ \rmn{yr}$ for $1\ \rmn{M_J}$ and $6\ \rmn{M_J}$. Planets on inclined orbits warp the disc by an extent that is negligible for $1\ \rmn{M_J}$ but increases with increasing mass becoming quite significant for a planet of mass $6\ \rmn{M_J}$. We also find a solid body precession of both the total disc angular momentum vector and the planet orbital momentum vector about the total angular momentum vector. Our results illustrate that the influence of an inclined massive planet on a protoplanetary disc can lead to significant changes of the disc structure and orientation which can in turn affect the orbital evolution of the planet significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Planet-Disc Interactions and Early Evolution of Planetary Systems

The great diversity of extrasolar planetary systems has challenged our understanding of how planets form, and how their orbits evolve as they form. Among the various processes that may account for this diversity, the gravitational interaction between planets and their parent protoplanetary disc plays a prominent role in shaping young planetary systems. Planet-disc forces are large, and the characteristic times for the evolution of planets orbital elements are much shorter than the lifetime of protoplanetary discs. The determination of such forces is challenging, because it involves many physical mechanisms and it requires a detailed knowledge of the disc structure. Yet, the intense research of the past few years, with the exploration of many new avenues, represents a very significant improvement on the state of the discipline. This chapter reviews current understanding of planet-disc interactions, and highlights their role in setting the properties and architecture of observed planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Type I planet migration in weakly magnetised laminar discs

[Abridged] The migration of low mass planets has been studied in hydrodynamical disc models for more than three decades, but the impact of a magnetic field in the protoplanetary disc is less known. When the disc's magnetic field is strong enough to prevent horseshoe motion, the corotation torque is replaced by a torque arising from magnetic resonances. For weak enough magnetic fields, horseshoe motion and a corotation torque exist, and recent turbulent MHD simulations have reported the existence of a new component of the corotation torque in the presence of a mean toroidal field. The aim of this paper is to investigate the physical origin and the properties of this new corotation torque. We performed MHD simulations of a low mass planet embedded in a 2D laminar disc threaded by a weak toroidal magnetic field, with the effects of turbulence modelled by a viscosity and a resistivity. We confirm that the interaction between the magnetic field and the horseshoe motion results in an additional corotation torque on the planet, which we dub the MHD torque excess. It is caused by the accumulation of the magnetic field along the downstream separatrices of the planet's horseshoe region, which gives rise to an azimuthally asymmetric underdense region at that location. The properties of the MHD torque excess are characterised by varying the slope of the density, temperature and magnetic field profiles, as well as the diffusion coefficients and the strength of the magnetic field. The sign of the MHD torque excess depends on the density and temperature gradients only, and is positive for profiles expected in protoplanetary discs. Its magnitude is in turn mainly determined by the strength of the magnetic field and the turbulent resistivity. The MHD torque excess can be strong enough to reverse migration, even when the magnetic pressure is less than one percent of the thermal pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hysteresis and thermal limit cycles in MRI simulations of accretion discs

The recurrent outbursts that characterise low-mass binary systems reflect thermal state changes in their associated accretion discs. The observed outbursts are connected to the strong variation in disc opacity as hydrogen ionises near 5000 K. This physics leads to accretion disc models that exhibit bistability and thermal limit cycles, whereby the disc jumps between a family of cool and low accreting states and a family of hot and efficiently accreting states. Previous models have parametrised the disc turbulence via an alpha (or `eddy') viscosity. In this paper we treat the turbulence more realistically via a suite of numerical simulations of the magnetorotational instability (MRI) in local geometry. Radiative cooling is included via a simple but physically motivated prescription. We show the existence of bistable equilibria and thus the prospect of thermal limit cycles, and in so doing demonstrate that MRI-induced turbulence is compatible with the classical theory. Our simulations also show that the turbulent stress and pressure perturbations are only weakly dependent on each other on orbital times; as a consequence, thermal instability connected to variations in turbulent heating (as opposed to radiative cooling) is unlikely to operate, in agreement with previous numerical results. Our work presents a first step towards unifying simulations of full MHD turbulence with the correct thermal and radiative physics of the outbursting discs associated with dwarf novae, low-mass X-ray binaries, and possibly young stellar objects.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Orbital evolution of a planet on an inclined orbit interacting with a disc

We study the dynamics of a planet on an orbit inclined with respect to a disc. If the initial inclination of the orbit is larger than some critical value, the gravitational force exerted by the disc on the planet leads to a Kozai cycle in which the eccentricity of the orbit is pumped up to large values and oscillates with time in antiphase with the inclination. On the other hand, both the inclination and the eccentricity are damped by the frictional force that the planet is subject to when it crosses the disc. We show that, by maintaining either the inclination or the eccentricity at large values, the Kozai effect provides a way of delaying alignment with the disc and circularization of the orbit. We find the critical value to be characteristically as small as about 20 degrees. Typically, Neptune or lower mass planets would remain on inclined and eccentric orbits over the disc lifetime, whereas orbits of Jupiter or higher mass planets would align and circularize. This could play a significant role in planet formation scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Weakly nonlinear theory for spiral density waves excited by accretion disc turbulence

We develop an analytic theory to describe spiral density waves propagating in a shearing disc in the weakly nonlinear regime. Such waves are generically found to be excited in simulations of turbulent accretion disks, in particular if said turbulence arises from the magneto-rotational instability (MRI). We derive a modified Burgers equation governing their dynamics, which includes the effects of nonlinear steepening, dispersion, and a bulk viscosity to support shocks. We solve this equation approximately to obtain nonlinear sawtooth solutions that are asymptotically valid at late times. In this limit, the presence of shocks is found to cause the wave amplitude to decrease with time as 1/t^2. The validity of the analytic description is confirmed by direct numerical solution of the full nonlinear equations of motion. The asymptotic forms of the wave profiles of the state variables are also found to occur in MRI simulations indicating that dissipation due to shocks plays a significant role apart from any effects arising from direct coupling to the turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the dynamics of multiple systems of hot super-Earths and Neptunes: Tidal circularization, resonance and the HD 40307 system

[Abridged] We consider the dynamics of a system of hot super-Earths or Neptunes such as HD 40307. We show that, as tidal interaction with the central star leads to small eccentricities, the planets in this system could be undergoing resonant coupling even though the period ratios depart significantly from very precise commensurability. In a three planet system, this is indicated by the fact that resonant angles librate or are associated with long term changes to the orbital elements. We propose that the planets in HD 40307 were in a strict Laplace resonance while they migrated through the disc. After entering the disc inner cavity, tidal interaction would cause the period ratios to increase from two but with the inner pair deviating less than the outer pair, counter to what occurs in HD 40307. However, the relationship between these pairs that occurs in HD 40307 might be produced if the resonance is impulsively modified by an event like a close encounter shortly after the planetary system decouples from the disc. We find this to be in principle possible for a small relative perturbation on the order of a few 1.d-3 but then we find that evolution to the present system in a reasonable time is possible only if the masses are significantly larger than the minimum masses and tidal dissipation is very effective. On the other hand we found that a system like HD 40307 with minimum masses and more realistic tidal dissipation could be produced if the eccentricity of the outermost planet was impulsively increased to about 0.15. The form of resonantly coupled tidal evolution we consider here is quite general and could be of greater significance for systems with inner planets on significantly shorter orbital periods characteristic of for example CoRoT 7 b.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stochastic orbital migration of small bodies in Saturn's rings

Many small moonlets, creating propeller structures, have been found in Saturn's rings by the Cassini spacecraft. We study the dynamical evolution of such 20-50m sized bodies which are embedded in Saturn's rings. We estimate the importance of various interaction processes with the ring particles on the moonlet's eccentricity and semi-major axis analytically. For low ring surface densities, the main effects on the evolution of the eccentricity and the semi-major axis are found to be due to collisions and the gravitational interaction with particles in the vicinity of the moonlet. For large surface densities, the gravitational interaction with self-gravitating wakes becomes important. We also perform realistic three dimensional, collisional N-body simulations with up to a quarter of a million particles. A new set of pseudo shear periodic boundary conditions is used which reduces the computational costs by an order of magnitude compared to previous studies. Our analytic estimates are confirmed to within a factor of two. On short timescales the evolution is always dominated by stochastic effects caused by collisions and gravitational interaction with self-gravitating ring particles. These result in a random walk of the moonlet's semi-major axis. The eccentricity of the moonlet quickly reaches an equilibrium value due to collisional damping. The average change in semi-major axis of the moonlet after 100 orbital periods is 10-100m. This translates to an offset in the azimuthal direction of several hundred kilometres. We expect that such a shift is easily observable.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Type III migration in a low viscosity disc

We study the type III migration of a Saturn mass planet in low viscosity discs. The planet is found to experience cyclic episodes of rapid decay in orbital radius, each amounting to a few Hill radii. We find this to be due to the scattering of large- scale vortices present in the disc. The origin and role of vortices in the context of type III migration is explored. It is shown through numerical simulations and semi- analytical modelling that spiral shocks induced by a sufficiently massive planet will extend close to the planet orbital radius. The production of vortensity across shock tips results in thin high vortensity rings with a characteristic width of the local scale height. For planets with masses equal to and above that of Saturn, the rings are co-orbital features extending the entire azimuth. Linear stability analysis show there exists unstable modes that are localised about local vortensity minima which coincide with gap edges. Simulations show that vortices are non-linear a outcome. We used hydrodynamic simulations to examine vortex-planet interactions. Their effect is present in discs with kinematic viscosity less than about an order of magnitude smaller than the typically adopted value of ν= 10^{-5}Ω_pr_p(0)^2, where r_p(0) and Ω_p are the initial orbital radius and angular velocity of the planet respectively. We find that the magnitude of viscosity affects the nature of type III migration but not the extent of the orbital decay. The role of vortices as a function of initial disc mass is also explored and it is found that the amount of orbital decay during one episode of vortex-planet interaction is independent of initial disc mass. We incorporate the concept of the co-orbital mass deficit in the analysis of our results and link it to the presence of vortices at gap edges.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Conditions for the occurrence of mean-motion resonances in a low mass planetary system

The dynamical interactions that occur in newly formed planetary systems may reflect the conditions occurring in the protoplanetary disk out of which they formed. With this in mind, we explore the attainment and maintenance of orbital resonances by migrating planets in the terrestrial mass range. Migration time scales varying between millions of years and thousands of years are considered. In the former case, for which the migration time is comparable to the lifetime of the protoplanetary gas disk, a 2:1 resonance may be formed. In the latter, relatively rapid migration regime commensurabilities of high degree such as 8:7 or 11:10 may be formed. However, in any one large-scale migration several different commensurabilities may be formed sequentially, each being associated with significant orbital evolution. We also use a simple analytic theory to develop conditions for first order commensurabilities to be formed. These depend on the degree of the commensurability, the imposed migration and circularization rates, and the planet mass ratios. These conditions are found to be consistent with the results of our simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dynamical architectures of planetary systems induced by orbital migration

The aim of this talk is to present the most recent advances in establishing plausible planetary system architectures determined by the gravitational tidal interactions between the planets and the disc in which they are embedded during the early epoch of planetary system formation. We concentrate on a very well defined and intensively studied process of the disc-planet interaction leading to the planet migration. We focus on the dynamics of the systems in which low-mass planets are present. Particular attention is devoted to investigation of the role of resonant configurations. Our studies, apart from being complementary to the fast progress occurring just now in observing the whole variety of planetary systems and uncovering their structure and origin, can also constitute a valuable contribution in support of the missions planned to enhance the number of detected multiple systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Mean-Motion Resonances of High Order in Extrasolar Planetary Systems

Many multi-planet systems have been discovered in recent years. Some of them are in mean-motion resonances (MMR). Planet formation theory was successful in explaining the formation of 2:1, 3:1 and other low resonances as a result of convergent migration. However, higher order resonances require high initial orbital eccentricities in order to be formed by this process and these are in general unexpected in a dissipative disk. We present a way of generating large initial eccentricities using additional planets. This procedure allows us to form high order MMRs and predict new planets using a genetic N-body code.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Dynamical Origin of the Multi-Planetary System HD45364

The recently discovered planetary system HD45364 which consists of a Jupiter and Saturn mass planet is very likely in a 3:2 mean motion resonance. The standard scenario to form planetary commensurabilities is convergent migration of two planets embedded in a protoplanetary disc. When the planets are initially separated by a period ratio larger than two, convergent migration will most likely lead to a very stable 2:1 resonance for moderate migration rates. To avoid this fate, formation of the planets close enough to prevent this resonance may be proposed. However, such a simultaneous formation of the planets within a small annulus, seems to be very unlikely. Rapid type III migration of the outer planet crossing the 2:1 resonance is one possible way around this problem. In this paper, we investigate this idea in detail. We present an estimate for the required convergent migration rate and confirm this with N-body and hydrodynamical simulations. If the dynamical history of the planetary system had a phase of rapid inward migration that forms a resonant configuration, we predict that the orbital parameters of the two planets are always very similar and hence should show evidence of that. We use the orbital parameters from our simulation to calculate a radial velocity curve and compare it to observations. Our model can explain the observational data as good as the previously reported fit. The eccentricities of both planets are considerably smaller and the libration pattern is different. Within a few years, it will be possible to observe the planet-planet interaction directly and thus distinguish between these different dynamical states.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Protoplanet Magnetosphere Interactions

In this paper, we study a simple model of an orbiting protoplanet in a central magnetospheric cavity, the entry into such a cavity having been proposed as a mechanism for halting inward orbital migration. We have calculated the gravitational interaction of the protoplanet with the magnetosphere using a local model and determined the rate of evolution of the orbit. The interaction is found to be determined by the outward flux of MHD waves and thus the possibility of the existence of such waves in the cavity is significant. The estimated orbital evolution rates due to gravitational and other interactions with the magnetosphere are unlikely to be significant during protoplanetary disk lifetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

3D-MHD simulations of an accretion disk with star-disk boundary layer

We present global 3D MHD simulations of geometrically thin but unstratified accretion disks in which a near Keplerian disk rotates between two bounding regions with initial rotation profiles that are stable to the MRI. The inner region models the boundary layer between the disk and an assumed more slowly rotating central, non magnetic star. We investigate the dynamical evolution of this system in response to initial vertical and toroidal fields imposed in a variety of domains contained within the near Keplerian disk. Cases with both non zero and zero net magnetic flux are considered and sustained dynamo activity found in runs for up to fifty orbital periods at the outer boundary of the near Keplerian disk. Simulations starting from fields with small radial scale and with zero net flux lead to the lowest levels of turbulence and smoothest variation of disk mean state variables. For our computational set up, average values of the Shakura & Sunyaev (1973) $α$ parameter in the Keplerian disk are typically $0.004\pm 0.002.$ Magnetic field eventually always diffuses into the boundary layer resulting in the build up of toroidal field inward angular momentum transport and the accretion of disk material. The mean radial velocity, while exhibiting large temporal fluctuations is always subsonic. Simulations starting with net toroidal flux may yield an average $α\sim 0.04.$ While being characterized by one order of magnitude larger average $α$, simulations starting from vertical fields with large radial scale and net flux may lead to the formation of persistent non-homogeneous, non-axisymmetric magnetically dominated regions of very low density.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The Migration and Growth of Protoplanets in Protostellar Discs

We investigate the gravitational interaction of a Jovian mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions have been investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Nonlinear hydrodynamic simulations are employed using three independent numerical codes. A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near circular orbit after a characteristic viscous time scale of approximately 10,000 initial orbital periods. This was found to be independent of whether the protoplanet was allowed to accrete mass or not. Inward migration is helped through the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion onto the central star.Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is approximately 5 Jupiter masses, because of the reduced accretion rates obtained for planets of increasing mass. Assuming some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets (Marcy, Cochran, & Mayor 1999; Marcy & Butler 1998) as well as their inward orbital migration can be accounted for by consideration of disc--protoplanet interactions during the late stages of giant planet formation. Maximally accreting protoplanets reached about four Jovian masses on reaching the neighbourhood of the central star.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The hydrodynamical response of a tilted circumbinary disc: linear theory and non-linear numerical simulations

In this paper we present an analytical and numerical study of the response of a circumbinary disc subject to the tidal-forcing of a binary with a fixed circular orbit. We consider fluid discs with a range of thicknesses and binaries with a range of mass ratios, orbital separations and inclination angles. Our numerical simulations are implemented using a SPH code. For our unperturbed disc models, we write a scaling relation for the shear viscosity and deduce that the disc thickness cannot be varied without affecting the viscosity in these kinds of SPH disc models. It is found that maintainance of an inner cavity owing to the tidal truncation of the disc is effective for non-zero orbital inclinations. Also we show that our model discs may precess approximately like rigid bodies, provided that the disc is able to communicate on a length scale comparable to the inner boundary radius by either sonic or viscous effects, in a sufficiently small fraction of the local precession period. Furthermore, the disc precession period may tend to infinity if the disc outer edge is allowed to become arbitrarily large, the disc suffering only a modest quasi-steady warp near the inner boundary. When the disc response is linear, or weakly non-linear, the precession periods and the forms of warping that we measure yield reasonable quantitative agreement with the analytical expressions that we derive from a linear response calculation. For a stronger disc response the results can agree poorly with our linear analysis, although some qualitative features of the response remain intact. We show that the response of a disc of non-interacting particles is qualitatively different from this. The work presented in this paper is of relevance to a number of astrophysical phenomena of current interest in star and planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Variable Smoothing Lengths and Energy Conservation in Smoothed Particle Hydrodynamics

We present a new formulation of the equations of motion used in smoothed particle hydrodynamics (SPH). The spatial resolution in SPH is determined by the smoothing length, $h$, and it has become common practice for each particle to be given its own adaptive smoothing length, $h_i$. Consequently, the dynamic range that may be spatially resolved is greatly increased, but additional ($\nabla h$) terms, which account for the variability of the smoothing lengths, should appear in the equations of motion in order to satisy conservation requirements. Previous implementations of SPH have neglected these additional terms, whereas we have included them. This is achieved by defining a functional form for the $h_i$s, that depends on inter-particle distances only, and then deriving the equations of motion using a Hamiltonian formalism. A number of test calculations are presented. We find that including the $\nabla h$ terms has no detrimental effect on the ability of SPH to model known problems with reasonable accuracy. For problems where the energy conservation was rather poor, including the $\nabla h$ terms results in a dramatic improvement. In particluar non-conservation of energy during a collision between two polytropes can occur at the level of 10% when the $\nabla h$ terms are neglected. When the $\nabla h$ terms are included, then this error reduces to 0.8%.

0 citations0 reviewsNo code linkedOpen access