A physically realistic stellar model with a simple expression for the energy density and conformally flat interior is found. The relations between the different physical conditions are used without graphic proofs. It may represnet a real pulsar.

A simple classification is given of the anisotropic relativistic star models, resembling the one of charged isotropic solutions. On the ground of this database, and taking into account the conditions for physically realistic star models, a method is proposed for generating all such solutions. It is based on the energy density and the radial pressure as seeding functions. Numerous relations between the realistic conditions are found and the need for a graphic proof is reduced just to one pair of inequalities. This general formalism is illustrated with an example of a class of solutions with linear equation of state and simple energy density. It is found that the solutions depend on three free constants and concrete examples are given. Some other popular models are studied with the same method.

In order to study the type of collapse, mentioned in the title, we introduce a physically meaningful object, called the horizon function. It directly enters the expressions for many of the stellar characteristics. The main junction equation, which governs the collapse, transforms into a Riccati equation with simple coefficients for the horizon function. We integrate this equation in the geodesic case. The same is done in the general case when one or another of the coefficients vanish. It is shown how to build classes of star models in this formulation of the problem and simple solutions are given.

We introduce a physically important object, called the horizon function, in the study of geodesic collapse. It is closely related to the stellar characteristics and satisfies a simple Riccati equation. This equation is integrated and all of its solutions are found in terms of some generating functions. Previous solutions are regained and further investigated.

A global view is given upon the study of collapsing shear-free perfect fluid spheres with heat flow. We apply a compact formalism, which simplifies the isotropy condition and the condition for conformal flatness. This formalism also presents the simplest possible version of the main junction condition, demonstrated explicitly for conformally flat and geodesic solutions. It gives the right functions to disentangle this condition into well known differential equations like those of Abel, Riccati, Bernoulli and the linear one. It yields an alternative derivation of the general solution with functionally dependent metric components. We bring together the results for static and time- dependent models to describe six generating functions of the general solution to the isotropy equation. Their common features and relations between them are elucidated. A general formula for separable solutions is given, incorporating collapse to a black hole or to a naked singularity.

The fluid models mentioned in the title are classified. All characteristics of the fluid are expressed through a master potential, satisfying an ordinary second order differential equation. Different constraints are imposed on this core of relations, finding new solutions and deriving the classical results for perfect fluids and dust as particular cases. Many uncharged and charged anisotropic solutions, all conformally flat and some uniform density solutions are found. A number of solutions with linear equation among the two pressures are derived, including the case of vanishing tangential pressure.

The fluid models mentioned in the title are studied in a modified approach, based on two formulas for the mass function. All characteristics of the fluid are expressed through a master potential, satisfying an ordinary second order differential equation. Different constraints are imposed on this core of relations, finding new solutions and deriving the classical results for perfect fluids as particular cases. All charged anisotropic solutions, all conformally flat and all uniform density solutions are found. A large class of solutions with linear equation among the two pressures is derived, including the case of vanishing tangential pressure. The mechanism responsible for the appearance of equation of state is elucidated.

It is shown that an effective anisotropic spherically symmetric fluid model with heat flow can absorb the addition to a perfect fluid of pressure anisotropy, heat flow, bulk and shear viscosity, electric field and null fluid. In most cases the induction of effective heat flow can be avoided.

Algebraically special fields with no gravitational radiation are described. Kerr-Schild fields, which include as a concrete case the Kinnersley photon rocket, form an important subclass of them.

It is argued that static electric or magnetic fields induce Weyl-Majumdar-Papapetrou solutions for the metric of spacetime. Their gravitational acceleration includes a term many orders of magnitude stronger than usual perturbative terms. It gives rise to a number of effects, which can be detected experimentally. Four electrostatic and four magnetostatic examples of physical set-ups with simple symmetries are proposed. The different ways in which mass sources enter and complicate the pure electromagnetic picture are described.

It is shown that for realistic anisotropic star models the surface redshift can not exceed the values 3.842 or 5.211 when the tangential pressure satisfies the strong or the dominant energy condition respectively. Both values are higher than 2, the bound in the perfect fluid case.