Low thrust propulsion in a coplanar circular restricted four body problem

This paper formulates a circular restricted four body problem (CRFBP), where the three primaries are set in the stable Lagrangian equilateral triangle configuration and the fourth body is massless. The analysis of this autonomous coplanar CRFBP is undertaken, which identies eight natural equilibria; four of which are close to the smaller body, two stable and two unstable, when considering the primaries to be the Sun and two smaller bodies of the solar system. Following this, the model incorporates `near term' low-thrust propulsion capabilities to generate surfaces of articial equilibrium points close to the smaller primary, both in and out of the plane containing the celestial bodies. A stability analysis of these points is carried out and a stable subset of them is identied. Throughout the analysis the Sun-Jupiter-Asteroid-Spacecraft system is used, for conceivable masses of a hypothetical asteroid set at the libration point L4. It is shown that eight bounded orbits exist, which can be maintained with a constant thrust less than 1:5 10􀀀4N for a 1000kg spacecraft. This illustrates that, by exploiting low-thrust technologies, it would be possible to maintain an observation point more than 66% closer to the asteroid than that of a stable natural equilibrium point. The analysis then focusses on a major Jupiter Trojan: the 624-Hektor asteroid. The thrust required to enable close asteroid observation is determined in the simplied CRFBP model. Finally, a numerical simulation of the real Sun-Jupiter-624 Hektor-Spacecraft is undertaken, which tests the validity of the stability analysis of the simplied model

Natural and sail-displaced doubly-symmetric Lagrange point orbits for polar coverage

This paper proposes the use of doubly-symmetric, eight-shaped orbits in the circular restricted three-body problem for continuous coverage of the high-latitude regions of the Earth. These orbits, for a range of amplitudes, spend a large fraction of their period above either pole of the Earth. It is shown that they complement Sun-synchronous polar and highly eccentric Molniya orbits, and present a possible alternative to low thrust pole-sitter orbits. Both natural and solar-sail displaced orbits are considered. Continuation methods are described and used to generate families of these orbits. Starting from ballistic orbits, other families are created either by increasing the sail lightness number, varying the period or changing the sail attitude. Some representative orbits are then chosen to demonstrate the visibility of high-latitude regions throughout the year. A stability analysis is also performed, revealing that the orbits are unstable: it is found that for particular orbits, a solar sail can reduce their instability. A preliminary design of a linear quadratic regulator is presented as a solution to stabilize the system by using the solar sail only. Finally, invariant manifolds are exploited to identify orbits that present the opportunity of a ballistic transfer directly from low Earth orbit

Reconstructing phylogenetic level-1 networks from nondense binet and trinet sets

Binets and trinets are phylogenetic networks with two and three leaves, respectively. Here we consider the problem of deciding if there exists a binary level-1 phylogenetic network displaying a given set T of binary binets or trinets over a taxon set X, and constructing such a network whenever it exists. We show that this is NP-hard for trinets but polynomial-time solvable for binets. Moreover, we show that the problem is still polynomial-time solvable for inputs consisting of binets and trinets as long as the cycles in the trinets have size three. Finally, we present an O(3^{|X|} poly(|X|)) time algorithm for general sets of binets and trinets. The latter two algorithms generalise to instances containing level-1 networks with arbitrarily many leaves, and thus provide some of the first supernetwork algorithms for computing networks from a set of rooted 1 phylogenetic networks

Neighborhoods of trees in circular orderings

In phylogenetics, a common strategy used to construct an evolutionary tree for a set of species X is to search in the space of all such trees for one that optimizes some given score function (such as the minimum evolution, parsimony or likelihood score). As this can be computationally intensive, it was recently proposed to restrict such searches to the set of all those trees that are compatible with some circular ordering of the set X. To inform the design of efficient algorithms to perform such searches, it is therefore of interest to find bounds for the number of trees compatible with a fixed ordering in the neighborhood of a tree that is determined by certain tree operations commonly used to search for trees: the nearest neighbor interchange (nni), the subtree prune and regraft (spr) and the tree bisection and reconnection (tbr) operations. We show that the size of such a neighborhood of a binary tree associated with the nni operation is independent of the tree’s topology, but that this is not the case for the spr and tbr operations. We also give tight upper and lower bounds for the size of the neighborhood of a binary tree for the spr and tbr operations and characterize those trees for which these bounds are attained

Triangleland. I. Classical dynamics with exchange of relative angular momentum

In Euclidean relational particle mechanics, only relative times, relative angles and relative separations are meaningful. Barbour--Bertotti (1982) theory is of this form and can be viewed as a recovery of (a portion of) Newtonian mechanics from relational premises. This is of interest in the absolute versus relative motion debate and also shares a number of features with the geometrodynamical formulation of general relativity, making it suitable for some modelling of the problem of time in quantum gravity. I also study similarity relational particle mechanics (`dynamics of pure shape'), in which only relative times, relative angles and {\sl ratios of} relative separations are meaningful. This I consider firstly as it is simpler, particularly in 1 and 2 d, for which the configuration space geometry turns out to be well-known, e.g. S^2 for the `triangleland' (3-particle) case that I consider in detail. Secondly, the similarity model occurs as a sub-model within the Euclidean model: that admits a shape--scale split. For harmonic oscillator like potentials, similarity triangleland model turns out to have the same mathematics as a family of rigid rotor problems, while the Euclidean case turns out to have parallels with the Kepler--Coulomb problem in spherical and parabolic coordinates. Previous work on relational mechanics covered cases where the constituent subsystems do not exchange relative angular momentum, which is a simplifying (but in some ways undesirable) feature paralleling centrality in ordinary mechanics. In this paper I lift this restriction. In each case I reduce the relational problem to a standard one, thus obtain various exact, asymptotic and numerical solutions, and then recast these into the original mechanical variables for physical interpretation.Comment: Journal Reference added, minor updates to References and Figure

Statistical analysis of the ion flux to the JET outer wall

Statistical analysis of the ion flux to the JET outer-wall is conducted in outer-wall limiter mounted Langmuir probe (OLP) time-series across a wide range of plasma current and line-averaged density during Ohmically heated horizontal target L-mode plasmas. The mean, μ, and the standard deviation, σ, of the ion-saturation current measured by the OLP show systematic variation with plasma current and density. Both increase as either plasma current decreases and/or density increases. Upon renormalization, achieved by subtraction of μ and rescaling by σ, the probability distribution functions (PDFs) of each signal collapse approximately onto a single curve. The shape of the curve deviates from a Γ distribution in the tail of the PDF and is better described by a log-normal distribution. The invariance in the shape of the PDF, which occurs over approximately four decades of the ordinate, is shown to be the result of a balance between the duration time of the average burst wave-form, τd and the waiting time between bursts, τw. This implies that the intermittency parameter, τd/τw , can be considered constant at the JET outer wall during horizontal target Ohmic L-mode operation. This result may be important both for model validation and prediction