Giant HVC’s and the Rotation Curves of the Milky Way and M31

We have modelled, for the cases of Milky Way and M31, the effects on the galactic discs, of the arrival at high velocity (≥150 km s⁻¹) of giant HI clouds, with masses of up to 10⁸M⊙. Predictions are compared with the detailed structure of the observed rotation curves for these two galaxies. The model explains the rises and falls observed at large distances from the centre of each galaxy, distributed with a degree of regularity in radius, in terms of a specific type of perturbations driven by the infall of the high velocity clouds (HVC's) arriving from the intracluster medium of the Local Group. The underlying rotation curve is explained conventionally via the distribution of the baryonic and dark matter components of the galaxy in question. This scenario, though tested here on the two major Local Group objects, is in principle applicable to galaxies undergoing minor mergers with subgalactic mass gas clouds.     

附有限制条件的间接平差秩亏时解法初探

在附有限制条件的间接平差的应用中，若所列观测方程中实际用到的未知参数的个数少于所选参数的总数u时，就会出现法方程系数阵秩亏的问题。为了解决此问题，在理论公式推导的基础上，给出了一种新的解决方法，并结合算例进行了验证。     

A two-dimensional finite volume hydrodynamic model for coastal areas: Model development and validation

A two dimensional implicit finite volume scheme for solving the shallow-water equations is developed. The effects of the Coriolis force, surface wind stress, and waves are included. A non-uniform rectilinear forward staggered grid is used with Cartesian coordinates. The time integration is performed using the Euler implicit technique. The convective flux is treated using the deferred correction method. The viscous terms are discretized using a second order central difference approximation. The SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm is used for coupling the velocity components and the water elevation gradient for the water level correction. The system of equations is solved sequentially using the Strongly Implicit Procedure (SIP). To simulate wave driven current, a phase averaged wave model is used first to simulate wave transformation and calculate radiation stresses. The performance of the developed model is validated for different sources of external forces and different combinations of boundary conditions. The validation cases include tidal circulation in a harbor and wave induced currents behind a breakwater parallel to the coastline. The model is finally applied to simulate the flow pattern in a closed artificial lagoon and along the coastline near Damietta Port located along the Northern coast of Egypt. Results of the developed model agree well with the published results for the considered cases.     

Flows and shocks in active galaxies and their nuclei

Observations of active galactic nuclei imply that shocks must be an essential and important part of their structure. We outline the basic observations, and discuss those features which must be addressed by any physical model of active nuclei. These features, in particular the observed spectrum of strong emission and absorption lines, lead naturally to the conclusion that shocks are present. The velocity widths of these lines, which range from hundreds to many thousands of kilometres per second, are most readily explained by models in which shocks play an important role in the generation of cool gas. The extreme parameters of the shocks in and around active nuclei provide a unique application for the physics discussed in this meeting.     

Parameter estimation for basin-scale ecosystem-linked population models of large pelagic predators: Application to skipjack tuna

A Spatial Ecosystem and Population Dynamic Model (SEAPODYM) is used in a data assimilation study aiming to estimate model parameters that describe dynamics of Pacific skipjack tuna population on ocean-based scale. The model based on advection–diffusion–reaction equations explicitly predicts spatial dynamics of large pelagic predators, while taking into account data on several mid-trophic level components, oceanic primary productivity and physical environment. In order to improve its quantitative ability, the model was parameterized through assimilation with commercial fisheries data, and optimization was carried out using maximum likelihood estimation approach. To address the optimization task we implemented an adjoint technique to obtain an exact, analytical evaluation of the likelihood gradient. We conducted a series of computer experiments in order to (i) determine model sensitivity with respect to variable parameters and, hence, investigate their observability; (ii) estimate observable parameters and their errors; and (iii) justify the reliability of the computed solution. Parameters describing recruitment, movement, habitat preferences, natural and fishing mortality of skipjack population were analysed and estimated. Results of the study suggest that SEAPODYM with achieved parameterization scheme can help to investigate the impact of fishing under various management scenarios, and also conduct forecasts of a given species stock and spatial dynamics in a context of environmental and climate changes.     

Lambert problem solution in the hill model of motion

The goal of this paper is obtaining a solution of the Lambert problem in the restricted three-body problem described by the Hill equations. This solution is based on the use of pre determinate reference orbits of different types giving the first guess and defining the sought-for transfer type. A mathematical procedure giving the Lambert problem solution is described. This procedure provides step-by-step transformation of the reference orbit to the sought-for transfer orbit. Numerical examples of the procedure application to the transfers in the Sun–Earth system are considered. These examples include transfer between two specified positions in a given time, a periodic orbit design, a halo orbit design, halo-to-halo transfers, LEO-to-halo transfer, analysis of a family of the halo-to-halo transfer orbits. The proposed method of the Lambert problem solution can be used for the two-point boundary value problem solution in any model of motion if a set of typical reference orbits can be found.     

Behaviour of the SMF method for the numerical integration of satellite orbits

The SMF algorithms were recently developed by the authors as a multistep generalization of the ScheifeleG-functions one-step method. Like the last, the proposed codes integrate harmonic oscillations without truncation error and the perturbing parameter appears as a factor of that error when integrating perturbed oscillations. Therefore they seemed to be convenient for the accurate integration of orbital problems after the application of linearizing transformations, such as KS or BF. In this paper we present several numerical experiments concerning the propagation of Earth satellite orbits, that illustrate the performance of the the SMF method. In general, it provides greater accuracy than the usual standard algorithms for similar computational cost.     

A Time-Transformed Leapfrog Scheme

We present a time-transformed leapfrog scheme combined with the extrapolation method to construct an integrator for orbits in N-body systems with large mass ratios. The basic idea can be used to transform any second-order differential equation into a form which may allow more efficient numerical integration. When applied to gravitating few-body systems this formulation permits extremely close two-body encounters to be considered without significant loss of accuracy. The new scheme has been implemented in a direct N-body code for simulations of super-massive binaries in galactic nuclei. In this context relativistic effects may also be included.