We compute multiprecision solutions of the Lane-Emden equation. This differential equation arises when introducing the well-known polytropic model into the equation of hydrostatic equilibrium for a nondistorted star. Since such multiprecision computations are time-consuming, we apply to this problem parallel programming techniques and thus the execution time of the computations is drastically reduced.

The so-called "global polytropic model" is applied to the numerical study of the exoplanet systems HD 10180, Kepler-32, Kepler-33, Kepler-102, and Kepler-186. We compare computed distances of planets from their host stars with corresponding observations and discuss some further orbit predictions made by the model.

According to the so-called "global polytropic model", hydrostatic equilibrium for a planetary system leads to the Lane-Emden differential equation. Solving this equation in the complex plane, we obtain polytropic spherical shells defined by succesive roots of the real part $\mathrm{Re}(\theta)$ of the Lane-Emden function $\theta$. Such shells provide hosting orbits for the members of a planetary system. Within the framework of the global polytropic model, we study the exoplanet systems HD 40307, $\mu$ Ara, Kepler-26, Kepler-62, and Kepler-275. We give emphasis on comparing our results with observations, and with orbit predictions derived from the "generalized Titius-Bode relation".

In this study, we consider the problem of geometry of curved jets. We work with a simplified two-dimensional kinematical model assuming that the jet velocity is varying. We prescribe the velocity variation and study two cases concerning parsec and kilo-parsec scale jets by applying some simplified numerical simulations. For parsec scale jets, we use some general expressions and derive the shape of jets for certain numerical cases. For kilo-parsec scale jets, we study the cases of 4C53.37 and 1159+583 of Abell 2220 and Abell 1446 galaxy clusters, respectively, in an effort to test our model and compare our numerical results with observational data.

In the framework of the so-called "global polytropic model", we assume hydrostatic equilibrium for a planetary system, and solve the resulting Lane--Emden differential equation in the complex plane. We thus obtain polytropic spherical shells defined by succesive roots of the real part $\mathrm{Re}(\theta)$ of the Lane-Emden function $\theta$. These shells seem to provide hosting orbits for the planets of the system(s) under consideration. In the present investigation, we study within this framework the exoplanet systems 55 Cnc, $\upsilon$ And, Kepler-11, Kepler-20, and Kepler-90.

A fundamental assumption in the so-called "global polytropic model" is hydrostatic equilibrium for a system of planets or statellites. By solving the Lane-Emden differential equation for such a system in the complex plane, we find polytropic spherical shells defined by succesive roots of the real part $\mathrm{Re}(\theta)$ of the Lane-Emden function $\theta$. These shells seem to be appropriate places for accomodating planets or satellites. In the present study, we apply the global polytropic model to the systems of satellites of the giant planets: Jupiter, Saturn, Uranus, and Neptune.

The so-called "global polytropic model" is based on the assumption of hydrostatic equilibrium for the solar system, or for a planet's system of statellites (like the jovian system), described by the Lane-Emden differential equation. A polytropic sphere of polytropic index $n$ and radius $R_1$ represents the central component $S_1$ (Sun or planet) of a polytropic configuration with further components the polytropic spherical shells $S_2$, $S_3$, ..., defined by the pairs of radii $(R_1,\,R_2)$, $(R_2,\,R_3)$, ..., respectively. $R_1,\,R_2,\,R_3,\, ...$, are the roots of the real part $\mathrm{Re}(\theta)$ of the complex Lane-Emden function $\theta$. Each polytropic shell is assumed to be an appropriate place for a planet, or a planet's satellite, to be "born" and "Live". This scenario has been studied numerically for the cases of the solar and the jovian systems. In the present paper, the Lane-Emden differential equation is solved numerically in the complex plane by using the Fortran code DCRKF54 (modified Runge-Kutta-Fehlberg code of fourth and fifth order for solving initial value problems in the complex plane along complex paths). We include in our numerical study some trans-Neptunian objects.

We study numerically the so-called ``turn-over scenario'' for rotating magnetic white dwarfs and neutron stars. According to this scenario, the magnetic symmetry axis of the model inclines at a gradually increasing angle (the so-called ``turn-over angle'') relative to the invariant angular momentum axis. Consequently, the model becomes ``perpendicular rotator'' (i.e., its turn-over angle becomes almost 90 degrees) on a ``turn-over timescale'' calculated to be ~ (a few)x10^6 - (a few)x10^7 yr for the examined white dwarf models and ~ 10^1 - 10^3 yr for the examined neutron star models. Furthermore, the initial differential rotation of the model turns to uniform rotation due to angular momentum mixing caused by hydrodynamic Alfven waves propagating along the poloidal magnetic field lines. Our numerical results show that, during the turn-over phase, the spin-down time rate is large, while the spin-down power remains small; so, the turn-over phase is a characteristic case for which an eventually observed large spin-down time rate should not be interpreted as implying a large spin-down power.

When a differentially rotating magnetic star undergoes ``turn-over'', i.e., its magnetic symmetry axis is inclining at a gradually increasing angle with respect to its invariant angular momentum axis up to the perpendicular position (situation known as ``perpendicular rotator''), the poloidal magnetic field causes interchange of angular momentum via hydromagnetic Alfven waves propagating along the poloidal magnetic field lines and, accordingly, the resulting perpendicular rotator obtains uniform rotation. This dissipative mechanism is represented by the so-called ``turn-over viscosity''. For several white dwarf and neutron star models, we compute the energy dissipation due to the turn-over viscosity; we also compare this effect with the corresponding dissipative effects induced by the electron, ion, and turbulent viscosity.

We study a white dwarf with differential rotation and magnetic field, for which the symmetry axis of the toroidal field, the magnetic axis of the poloidal field, and the principal axis I_3 coincide permanently. This common axis, so-called "magnetic symmetry axis", inclines at a small angle relative to the spin axis of the model, the so-called "turn-over angle". Such an almost axisymmetric model undergoes an early evolutionary phase of secular timescale, with moment of inertia along the spin axis, I_{zz}=I_{33}, greater than the moments of inertia along the equatorial axes, I_{11}=I_{22}, since rotation and poloidal field (both deriving oblate configurations) dominate over the toroidal field (deriving prolate configurations). During this phase, the model suffers from secular angular momentum loss due to weak magnetic dipole radiation activated by the poloidal field, leading gradually to "dynamical asymmetry" with I_{11}>I_{33}. However, a dynamically asymmetric configuration tends to turn over spontaneously, i.e., to rotate about axis with moment of inertia greater than I_{33}, with angular momentum remaining invariant. So, the fate of a dynamically asymmetric configuration is to become an oblique rotator and, eventually, a perpendicular rotator. During the so-called "turn-over phase", the turn-over angle increases spontaneously up to 90 degrees on a "turn-over timescale" t_{TOV}. In the present paper, we study numerically the turn-over phase for models of several masses, angular momenta, and magnetic fields.

The white dwarf in AE Aqr is observed to spin down at a steady time rate ~ 5.64*10^{-14} s s^{-1}, while at the same time its UV and X-ray accretion luminosities remain almost unchanged. This is a contradiction, however, since the classicaly estimated spin-down power ~ -10^{34} erg s^{-1} exceeds the accretion luminosity of the primary by a factor ~ 10^2 and, as a dominating power, should lead either to observable luminosity changes or to other detectable effects. This so-called "spin-down problem" can be relaxed under the assumption that the primary is now in a phase of rapidly decreasing its differential rotation under constant angular momentum, undergoing a nonaxisymmetric transition (i.e., turning over its magnetic symmetry axis with respect to its angular momentum axis and eventually becoming a perpendicular rotator) from a "differential rotation state" (DRS) to a "rigid rotation state" (RRS) on a relatively short "nonaxisymmetric DRS-to-RRS transition timescale". If so, then the estimated spin-down power ~ -3*10^{32} erg s^{-1} becomes comparable to the observed luminosities and the spin-down problem is drastically simplified. We present a detailed numerical study of such a "turn-over scenario", which study mainly points out the fact that an observed large spin-down time rate does not always imply a large spin-down power.