We study the evolution of an isolated, spherical halo of self-interacting dark matter (SIDM) in the gravothermal fluid formalism. We show that the thermal relaxation time, $t_r$, of a SIDM halo with a central density and velocity dispersion of a typical dwarf galaxy is significantly shorter than its age. We find a self-similar solution for the evolution of a SIDM halo in the limit where the mean free path between collisions, $\lambda$, is everywhere longer than the gravitational scale height, $H$. Typical halos formed in this long mean free path regime relax to a quasistationary gravothermal density profile characterized by a nearly homogeneous core and a power-law halo where $\rho \propto r^{-2.19}$. We solve the more general time-dependent problem and show that the contracting core evolves to sufficiently high density that $\lambda$ inevitably becomes smaller than $H$ in the innermost region. The core undergoes secular collapse to a singular state (the ``gravothermal catastrophe'') in a time $t_{coll} \approx 290 t_r$, which is longer than the Hubble time for a typical dark matter-dominated galaxy core at the present epoch. Our model calculations are consistent with previous, more detailed, N-body simulations for SIDM, providing a simple physical interpretation of their results and extending them to higher spatial resolution and longer evolution times. At late times, mass loss from the contracting, dense inner core to the ambient halo is significantly moderated, so that the final mass of the inner core may be appreciable when it becomes relativistic and radially unstable to dynamical collapse to a black hole.

We have studied the distribution of forces in gravitational systems through numerical experiments. Data were taken from an N-body simulation in an expanding universe. Before clustering, the distribution of random forces was represented as a Holtsmark distribution; the nearest-neighbor distribution is also shown as a comparison. The analytical and simulation distributions are in good agreement. When clustering becomes strong, the simulation result showed that the contribution of the forces acting on each galaxy, which is generated from all other galaxies, is almost entirely due to the gravitational attraction of its nearest neighbor. This implies that nearest-neighbor galactic encounters may play the main role in the dynamics of galaxy clustering.

The evolution of tidally truncated globular clusters is investigated by integrating two-dimensional Fokker-Planck equation that allows the development of velocity anisotropy. We start from the isotropic Plummer model with tidal cut off and followed the evolution through the corecollapse. The heating by three-binary is included to obtain the evolution past the corecollapse. The anisotropy in velocity dispersion develops during the precollapse evolution. However, the anisotropy becomes highly depressed during the post-collapse evolution because of rapid loss of radial orbits. Maximum radial anisotropy appears just after the beginning of the expansion, and degree of anisotropy decreases slowly as the total mass of the cluster decreases. Thus it may be possible to determine the evolutionary status of a cluster if the velocity anisotropy can be measured in the sense that the postcollapse clusters always have very little degree of anisotropy. The structure of the post-collapse cluster can be well fitted to King models because the degree of anisotropy is rather small.

Violent relaxation process of spherical stellar systems is examined by numerical simulations of shell model. The collapse of uniform density sphere both with and without external force is investigated. It is found that time variation of mean gravitational potential field induced by external force makes the system closer to Lynden-Bell distribution. Our result imply that systems which are gravitationary affected by other systems may approach the Lynden-Bell distribution.

The following principle of minimum energy may be a powerful substitute to the dynamical perturbation method, when the latter is hard to apply. Fluid elements of self-gravitating barotropic flows, whose vortex lines extend to the boundary of the fluid, are labelled in such a way that any change of trial configurations automatically preserves mass and circulation. The velocity field is given by a mass conserving Clebsch representation. With three independent Lagrangian functions, the total energy is stationary for all small variations about a flow with fixed linear and angular momenta provided Euler's equations for steady motion are satisfied. Thus, steady flows are stable if their energy is minimum. Since energy is here minimized subject to having local and global contants of the motion fixed, stability limits obtained that way are expected to be close to limits given by dynamical perturbation methods. Moreover, the stability limits are with respect to arbitrary, not necessary small, perturbations. A weaker form of the energy principle is also given which may be easier to apply. The Lagrangian functional, with the same three Lagrange variables is stationary for the fully time dependent Euler equations. It follows that the principle of minimum energy gives stability conditions that are both necessary and sufficient if terms linear in time derivatives (gyroscopic terms) are absent from the Lagrangian. The gyroscopic term for small deviations around steady flows is given explicitly. Key words: Energy variational principle; Self-gravitating systems; Stability of fluids.

This paper describes the numerical evolution of an $N$-body system with a slight ``temperature inversion''; i.e. the maximum velocity dispersion occurs not at the centre but further out. Fluid models predict that the core of such a system expands on a time-scale of thousands of central relaxation times, and here this behaviour is qualitatively confirmed for an $N$-body system of over 3000 bodies. With certain qualifications, this demonstrates the existence in N-body systems of one of the fundamental mechanisms which, in fluid models, drive the gravothermal oscillations discovered by Bettwieser & Sugimoto.

We analyze stability conditions of "Maclaurin flows" (self-gravitating, barotropic, two dimensional, stationary streams moving in closed loops around a point) by minimizing their energy, subject to fixing all the constants of the motion including mass and circulations. Necessary and sufficient conditions of stability are obtained when gyroscopic terms in the perturbed Lagrangian are zero. To illustrate and check the properties of this new energy principle, we have calculated the stability limits of an ordinary Maclaurin disk whose dynamical stability limits are known. Perturbations are in the plane of the disk. We find all necessary and sufficient conditions of stability for single mode symmetrical or antisymmetrical perturbations. The limits of stability are identical with those given by a dynamical analysis. Regarding mixed types of perturbations the maximally constrained energy principle give for some the necessary and sufficient condition of stability, for others only sufficient conditions of stability. The application of the new energy principle to Maclaurin disks shows the method to be as powerful as the method of dynamical perturbations.