Here we extend the conservation of energy in the framework of the thin layer approximation to the asymmetrical case. Four types of interstellar medium are analysed, in which the density follows an inverse square profile, a power law profile, an exponential profile and a toroidal profile. An analytical solution for the radius as a function of time and the polar angle in spherical coordinates is derived in the case of the inverse square profile. The analytical and numerical results are applied to two supernova remnants: SN 1987A and SN 1006. The back reaction due to the radiative losses is evaluated in the case of the inverse square profile for the surrounding medium. Two models for the image formation are presented, which explain the triple ring visible in SN 1987A and the jet feature of SN 1006.

The thin layer approximation applied to the expansion of a supernova remnant assumes that all the swept mass resides in a thin shell. The law of motion in the thin layer approximation is therefore found using the conservation of momentum. Here we instead introduce the conservation of energy in the framework of the thin layer approximation. The first case to be analysed is that of an interstellar medium with constant density and the second case is that of 7 profiles of decreasing density with respect to the centre of the explosion. The analytical and numerical results are applied to 4 supernova remnants: Tycho, Cas A, Cygnus loop, and SN 1006. The back reaction due to the radiative losses for the law of motion is evaluated in the case of constant density of the interstellar medium.

The statistical parameters of five generalizations of the Lindley distribution, such as the average, variance and moments, are reviewed. A new double truncated Lindley distribution with three parameters is derived. The new distributions are applied to model the initial mass function for stars.

The shape of the local bubble is modeled in the framework of the thin layer approximation. The asymmetric shape of the local bubble is simulated by introducing axial profiles for the density of the interstellar medium, such as exponential, Gaussian, inverse square dependence and Navarro--Frenk--White. The availability of some observed asymmetric profiles for the local bubble allows us to match theory and observations via the observational percentage of reliability. The model is compatible with the presence of radioisotopes on Earth.

The gamma function is a good approximation to the luminosity function of astrophysical objects, and a truncated gamma distribution would permit a more rigorous analysis. This paper examines the generalized gamma distribution (GG) and then introduces the scale and the new double truncation. The magnitude version of the truncated GG distribution with scale is adopted in order to fit the luminosity function (LF) for galaxies or quasars. The new truncated GG LF is applied to the five bands of SDSS galaxies, to the 2dF QSO Redshift Survey in the range of redshifts between 0.3 and 0.5, and to the COSMOS QSOs in the range of redshifts between 3.7 and 4.7. The average absolute magnitude versus redshifts for SDSS galaxies and QSOs of 2dF was modeled adopting a redshift dependence for the lower and upper absolute magnitude of the new truncated GG LF.

The conservation of the energy flux in turbulent jets that propagate in the interstellar medium (ISM) allows us to deduce the law of motion when an inverse power law decrease of density is considered. The back-reaction that is caused by the radiative losses for the trajectory is evaluated. The velocity dependence of the jet with time/space is applied to the jet of HH34, for which the astronomical data of velocity versus time/space are available. The introduction of precession and constant velocity for the central star allows us to build a curved trajectory for the superjet connected with HH34. The bow shock that is visible in the superjet is explained in the framework of the theory of the image in the case of an optically thin layer.

The presence of the dark energy allows both the acceleration and the expansion of the universe. In the case of a constant equation of state for dark energy we derived an analytical solution for the Hubble radius in terms of the hypergeometric function. An approximate Taylor expansion of order seven is derived for both the constant and the variable equation of state for dark energy. In the case of the Cardassian cosmology we also derived an analytical solution for the Hubble radius in terms of the hypergeometric function. The astronomical samples of the distance modulus for Supernova (SN) of type Ia allows the derivation of the involved cosmological in the case of constant equation of state, variable equation of state and Cardassian cosmology.

The conservation of the energy flux in turbulent jets that propagate in the intergalactic medium (IGM) allows us to deduce the law of motion in the classical and relativistic cases. Four types of IGM are considered: constant density, hyperbolic decrease of density, inverse power law decrease of density and a Lane--Emden ($n=5$) profile. The conservation of the relativistic flux for the energy allows us to derive, to the first order, an analytical expression for the velocity. It also allows us to numerically determine the trajectory for the four types of medium. In the case of a Lane--Emden ($n=5$) profile, the back-reaction due to the radiative losses for the trajectory is evaluated both in the classical and the relativistic case. Astrophysical applications are made to the centerline intensity of the synchrotron emission and to the evolution of the magnetic field in the case of the radio-galaxy 3C31.

This paper reviews the Lindley distribution and then introduces the scale and the double truncation. The unknown parameters of the truncated Lindley distribution are evaluated with the maximum likelihood estimators. An application of the Lindley distribution with scale is done to the initial mass function for stars. The magnitude version of the Lindley distribution with scale is applied to the luminosity function for the Sloan Digital Sky Survey (SDSS) galaxies and to the photometric maximum of the 2MASS Redshift Survey (2MRS) galaxies. The truncated Lindley luminosity function allows us to model the Malquist bias of the 2MRS galaxies.

The giant arcs in the clusters of galaxies are modeled in the framework of the superbubbles. The density of the intracluster medium is assumed to follow a hyperbolic behavior. The analytical law of motion is function of the elapsed time and the polar angle. As a consequence the flux of kinetic energy in the expanding thin layer decreases with increasing polar angle making the giant arc invisible to the astronomical observations. In order to calibrate the arcsec-parsec conversion three cosmologies are analyzed.

A new analytical solution for the luminosity distance in flat $\Lambda$CDM cosmology is derived in terms of elliptical integrals of first kind with real argument. The consequent derivation of the distance modulus allows evaluating the Hubble constant, $H_0=69.77\pm 0.33$, $\om= 0.295\pm0.008$ and the cosmological constant, $\Lambda= (1.194 \pm 0.017)10^{-52} \frac{1}{m^2}$.

The intersections between a spherical shell and the faces of Voronoi's polyhedrons are numerically evaluated. The nodes of these intersections are the points that share the same distances from three nuclei. The nodes are assumed to be the places where we will find clusters of galaxies. The cosmological environment is given by the coupling between the number of galaxies as function of the redshift and the LCDM cosmology.

We consider a turbulent jet which is moving in a Lane--Emden ($n=5$) medium. The conserved quantity is the energy flux, which allows finding, to first order, an analytical expression for the velocity and an approximate trajectory. The conservation of the relativistic flux for the energy allows deriving, to first order, an analytical expression for the velocity, and numerically determining the trajectory. The back-reaction due to the radiative losses for the trajectory is evaluated both in the classical and the relativistic case.

We model the circumstellar medium with four density profiles: hyperbolic type, power law type, exponential type and Gaussian type. We solve analytically or numerically the four first-order differential equations which arise in the framework of the classical thin layer approximation. The non-cubic dependence of the swept mass with the advancing radius is also considered. We derive the equation of motion for the thin layer approximation in special relativity in two cases. The initial conditions are chosen in order to model the temporal evolution of SN 1987A over 23 years and of SN 1006 over 1000 years. We review the building blocks of the symmetrical and asymmetrical formations of the image.

The evaluation of the physical parameters or effects such as the luminosity function (LF) or the photometric maximum (PM) for galaxies is routinely modeled by their spectroscopic redshift. Here, we model LF and PM for galaxies by the photometric redshift as given by the Spitzer Wide-Area Infrared Extragalactic (SWIRE) catalog in the framework of a lognormal LF. In addition, we compare our model with the Schechter LF for galaxies. The adopted cosmological framework is that of the flat cosmology.

An enhancement in the number of galaxies as function of the redshift is visible on the SDSS Photometric Catalogue DR 12 at z=0.383. This over-density of galaxies is named the Great Wall. This variable number of galaxies as a function of the redshift can be explained in the framework of the luminosity function for galaxies. The differential of the luminosity distance in respect to the redshift is evaluated in the framework of the LCDM cosmology.

We derive some first order differential equations which model the classical and the relativistic thin layer approximations. The circumstellar medium is assumed to follow a density profile which can be exponential, Gaussian, Plummer-like, self-gravitating of Lane--Emden ($n=5$) type, or power law. The first order differential equations are solved analytically, or numerically, or by a series expansion, or by recursion, or by Pad\'e approximation. The initial conditions are chosen in order to model the temporal evolution of SN 1993J over ten years. The Pad\'e approximated equations of motion are applied to four SNRs: Tycho, Cas A, Cygnus loop, and SN~1006.

In order to model the Fermi bubbles we apply the theory of the superbubble (SB). A thermal model and a self-gravitating model are reviewed. We introduce a third model based on the momentum conservation of a thin layer which propagates in a medium with an inverse square dependence for the density. A comparison have been made between the sections of the three models and the section of an observed map of the Fermi bubbles. An analytical law for the SB expansion as function of the time and polar angle is deduced. We derive a new analytical result for the image formation of the Fermi bubbles in an elliptical framework.

Based on well known photometric rules for the number of galaxies as function of the distance in Mpc we model the so called "Great Wall" which is visible on the Cosmicflows-2 catalog. The gravitational field is evaluated at the light of the shell theorem and a finite value for the gravitational field is numerically derived.

The luminosity function (LF) for stars is here fitted by a Schechter function and by a Gamma probability density function. The dependence of the number of stars on the distance, both in the low and high luminosity regions, requires the inclusion of a lower and upper boundary in the Schechter and Gamma LFs. Three astrophysical applications for stars are provided: deduction of the parameters at low distances, behavior of the average absolute magnitude with distance, and the location of the photometric maximum as a function of the selected flux. The use of the truncated LFs allows to model the Malmquist bias.

A classical model based on a power law assumption for the radius-time relationship in the expansion of a Supernova (SN) allows to derive an analytical expression for the flow of mechanical kinetic energy and the time duration of Gamma-ray burst (GRB). A random process based on the ratio of two truncated lognormal distributions for luminosity and luminosity distance allows to derive the statistical distribution for time duration of GRBs. The high velocities involved in the first phase of expansion of a SN requires a relativistic treatment. The circumstellar medium is assumed to follow a density profile of Plummer type with eta=6. A series solution for the relativistic flow of kinetic energy allows to derive in a numerical way the duration time for GRBs. Here we analyse two cosmologies: the standard cosmology and the plasma cosmology.

The spiral structure of the Milky Way can be simulated by adopting percolation theory, where the active zones are produced by the evolution of many supernova (SN). Here we assume conversely that the percolative process is triggered by superbubbles (SB), the result of multiple SNs. A first thermal model takes into account a bursting phase which evolves in a medium with constant density, and a subsequent adiabatic expansion which evolves in a medium with decreasing density along the galactic height. A second cold model follows the evolution of an SB in an auto-gravitating medium in the framework of the momentum conservation in a thin layer. Both the thermal and cold models are compared with the results of numerical hydro-dynamics. A simulation of GW~46.4+5.5, the Gould Belt, and the Galactic Plane is reported. An elementary theory of the image, which allows reproducing the hole visible at the center of the observed SB, is provided.

The initial mass function for the stars is often modeled by a lognormal distribution. This paper is devoted to demonstrating the advantage of introducing a left and right truncated lognormal probability density function, which is characterized by four parameters. Its normalization constant, mean, the variance, second moment about the origin and distribution function are calculated. The chi-square test and the Kolmogorov--Smirnov test are performed on four samples of stars.

The luminosity function for quasars (QSOs) is usually fitted by a Schechter function. The dependence of the number of quasars on the redshift, both in the low and high luminosity regions, requires the inclusion of a lower and upper boundary in the Schechter function. The normalization of the truncated Schechter function is forced to be the same as that for the Schechter function, and an analytical form for the average value is derived. Three astrophysical applications for QSOs are provided: deduction of the parameters at low redshifts, behavior of the average absolute magnitude at high redshifts, and the location (in redshift) of the photometric maximum as a function of the selected apparent magnitude. The truncated Schechter function with the double power law and an improved Schechter function are compared as luminosity functions for QSOs. The chosen cosmological framework is that of the flat cosmology, for which we provided the luminosity distance, the inverse relation for the luminosity distance, and the distance modulus.

In this paper we present formulae for chord length distribution in the framework of Poissonian Voronoi Tessellation (PVT) and non Poissonian Voronoi Tessellation (NPVT). The introduction of the scale parameter in the obtained distributions allows us to model the chord for cosmic voids. A graphical comparison between cosmic voids visible on two catalogs of galaxies, 2dFGRS and VIPERS, and theoretical random chords is reported.

A superbubble which advances in a symmetric Navarro--Frenk--White density profile or in an auto-gravitating density profile generates a thick shell with a radius that can reach 10 kpc. The application of the symmetric and asymmetric image theory to this thick 3D shell produces a ring in the 2D map of intensity and a characteristic `U' shape in the case of 1D cut of the intensity. A comparison of such a ring originating from a superbubble is made with the Einstein's ring. A Taylor approximation of order 10 for the angular diameter distance is derived in order to deal with high values of the redshift.

In this paper we derive three equations of motion for a supernova remnant (SNR) in the framework of the thin layer approximation using the Pad\'e approximant. The circumstellar medium is assumed to follow a density profile of either an exponential type, a Gaussian type, or a Lane--Emden ($n=5$) type. The three equations of motion are applied to four SNRs: Tycho, Cas A, Cygnus loop, and SN~1006. The percentage error of the Pad\'e approximated solution is always less than $10\%$. The theoretical decrease of the velocity over ten years for SNRs is evaluated.

The determination of the luminosity function (LF) in gamma ray bursts (GRBs) depends on the adopted cosmology, each one characterized by its corresponding luminosity distance. Here we analyse three cosmologies: the standard cosmology, the plasma cosmology, and the pseudo-Euclidean universe. The LF of the GRBs is firstly modeled by the lognormal distribution and the four broken power law, and secondly by a truncated lognormal distribution. The truncated lognormal distribution fits acceptably the range in luminosity of GRBs as a function of the redshift.

We present an analytical solution for the luminosity distance in spatially flat cosmology with pressureless matter and the cosmological constant. The complex analytical solution is made of a real part and a negligible imaginary part. The real part of the luminosity distance allows finding the two parameters $H_0$ and $\om$. A simple expression for the distance modulus for SNs of type Ia is reported in the framework of the minimax approximation.

In this paper we derive some first order differential equations which model the classical and the relativistic thin layer approximations in the presence of a circumstellar medium with a density which is decreasing in the distance $z$ from the equatorial plane. The circumstellar medium is assumed to follow a density profile with $z$ of hyperbolic type, power law type, exponential type or Gaussian type. The first order differential equations are solved analytically, or numerically, or by a series expansion, or by Pad\'e approximants. The initial conditions are chosen in order to model the temporal evolution of SN 1987A over 23 years. The free parameters of the theory are found by maximizing the observational reliability which is based on an observed section of SN 1987A.

The luminosity distance in the standard cosmology as given by $\Lambda$CDM and consequently the distance modulus for supernovae can be defined by the Pad\'e approximant. A comparison with a known analytical solution shows that the Pad\'e approximant for the luminosity distance has an error of $4\%$ at redshift $= 10$. A similar procedure for the Taylor expansion of the luminosity distance gives an error of $4\%$ at redshift $=0.7 $; this means that for the luminosity distance, the Pad\'e approximation is superior to the Taylor series. The availability of an analytical expression for the distance modulus allows applying the Levenberg--Marquardt method to derive the fundamental parameters from the available compilations for supernovae. A new luminosity function for galaxies derived from the truncated gamma probability density function models the observed luminosity function for galaxies when the observed range in absolute magnitude is modeled by the Pad\'e approximant. A comparison of $\Lambda$CDM with other cosmologies is done adopting a statistical point of view.

The conservation of the energy flux in turbulent jets which propagate in the intergalactic medium (IGM) allows deducing the law of motion in the classical and relativistic cases. Three types of IGM are considered: constant density, hyperbolic and inverse power law decrease of density. An analytical law for the evolution of the magnetic field along the radio-jets is deduced using a linear relation between the magnetic pressure and the rest density. Astrophysical applications are made to the centerline intensity of synchrotron emission in NGC315 and to the magnetic field of 3C273.

Recently it could be shown ( that the impact crater size-frequency distribution of Pluto (based on an analysis of first images obtained by the recent New Horizons flyby) follows a power law alpha = 2.4926 in the interval of diameter (D) values ranging from 3.75 km to the largest determined value of 37.77 km. A reanalysis of this data set revealed that the whole crater SFD (i.e., with values in the interval of 1.2-37.7 km) can be described by a truncated Pareto distribution.

The number of galaxies at a given flux as a function of the redshift, $z$, is derived when the $z$-distance relation is non-standard. In order to compare different models, the same formalism is also applied to the standard cosmology. The observed luminosity function for galaxies of the zCOSMOS catalog at different redshifts is modelled by a new luminosity function for galaxies, which is derived by the truncated beta probability density function. Three astronomical tests, which are the photometric maximum as a function of the redshift for a fixed flux, the mean value of the redshift for a fixed flux, and the luminosity function for galaxies as a function of the redshift, compare the theoretical values of the standard and non-standard model with the observed value. The tests are performed on the FORS Deep Field (FDF) catalog up to redshift $z=1.5$ and on the zCOSMOS catalog extending beyond $z=4$. These three tests show minimal differences between the standard and the non-standard models.

A new set of laws of motion for turbulent jets propagating in an intergalactic medium characterized by a decreasing density is derived by applying conservation of momentum flux both in the classical and relativistic framework. Two characteristic features of radio-galaxies, such as oscillations and curvature, are modeled by a classical helicoidal jet. A third feature of a radio-galaxy, the appearance of knots, is explained as an effect due to the theory of images.

In order to explain light curve (LC) for Supernova (SN) we derive a classical formula for the conversion of the flux of kinetic energy into radiation. We then introduce a correction for the absorption adopting an optical depth as function of the time. The developed framework allows to fit the LC of type Ia SN 2005cf ( B and V ) and type IIp SN 2004A (B,V,I and R ). A relativistic formula for the flux of kinetic energy is also derived in terms of a Taylor expansion and the application is done to the LC of GRB 050814. The decay of the radioactive isotopes as a driver the LC for SNs is also reviewed and a new formulation is introduced. The Arnett's formula for bolometric luminosity is corrected for the optical depth and applied to SN 2001ay.

In this paper we derive some first order differential equations which model the classical and the relativistic thin layer approximations. The circumstellar medium is assumed to follow a density profile of Plummer type, or of Lane--Emden ($n=5$) type, or a power law. The first order differential equations are solved analytically, or numerically, or by a series expansion, or by recursion. The initial conditions are chosen in order to model the temporal evolution of SN 1993J over ten years and a smaller chi-squared is obtained for the Plummer case with eta=6. The stellar mass ejected by the SN progenitor prior to the explosion, expressed in solar mass, is identified with the total mass associated with the selected density profile and varies from $0.217 $ to $0.402 $ when the central number density is $10^7 $ particles per cubic centimeter. The Full width at half maximum of the three density profiles, which can be identified with the size of the Pre-SN 1993J envelope, varies from 0.0071 pc to 0.0092 pc.

The all-sky maps of the observed variation of the Hubble constant can be reproduced from a theoretical point of view by introducing an intergalactic plasma with a variable number density of electrons. The observed averaged value and variance of the Hubble constant are reproduced by adopting a rim model, an auto-gravitating model, and a Voronoi diagrams model as the backbone for an auto-gravitating medium. We also analyze an astronomer's model based on the 3D spatial distribution of galaxies as given by the 2MASS Redshift Survey and an auto-gravitating Lane--Emden ($n=5$) profile of the electrons. The simulation which involves the Voronoi diagrams is done in a cubic box with sides of 100 Mpc. The simulation which involves the 2MASS covers the range of redshift smaller than 0.05.

The high velocities observed in supernovae require a relativistic treatment for the equation of motion in the presence of gradients in the density of the interstellar medium. The adopted theory is that of the thin layer approximation. The chosen medium is auto-gravitating with respect to an equatorial plane. The differential equation which governs the relativistic conservation of momentum is solved in numerically and by recursion. The asymmetric field of relativistic velocities as well the time dilation are plotted at the age of 1 yr for SN 1987A.

Two recent luminosity function (LF) for galaxies are reviewed and the parameters which characterize the near infrared are fixed. A first LF is a modified Schechter LF with four parameters. The second LF is derived from the generalized gamma and has four parameters. The formulas which give the number of galaxies as function of the redshift are reviewed and a special attention is given to the position of the photometric maximum which is expressed as function of a critical parameter or the flux of radiation or the apparent magnitude. A simulation of the 2MASS Redshift Survey is given in the framework of the non Poissonian Voronoi Tessellation.

We present a new formula which models the rate of decline of supernovae (SN) as given by the light curve in various bands. The physical basis is the conversion of the flux of kinetic energy into radiation. The two main components of the model are a power law dependence for the radius--time relation and a decreasing density with increasing distance from the central point. The new formula is applied to SN 1993J, SN 2005cf, SN 1999ac, and SN 2001el in different bands.

A first new luminosity functions of galaxies can be built starting from a left truncated beta probability density function, which is characterized by four parameters. In the astrophysical conversion, the number of parameters increases by one, due to the addition of the overall density of galaxies. A second new galaxy luminosity function is built starting from a left truncated beta probability for the mass of galaxies once a simple nonlinear relationship between mass and luminosity is assumed; in this case the number of parameters is six because the overall density of galaxies and a parameter that regulates mass and luminosity are added. The two new galaxy luminosity functions with finite boundaries were tested on the Sloan Digital Sky Survey (SDSS) in five different bands; the results produce a "better fit" than the Schechter luminosity function in two of the five bands considered. A modified Schechter luminosity function with four parameters has been also analyzed.

The pair correlation function (PCF) for galaxies presents typical oscillations in the range 20-200 Mpc/h which are named baryon acoustic oscillation (BAO). We first review and test the oscillations of the PCF when the 2D/3D vertexes of the Poissonian Voronoi Tessellation (PVT) are considered. We then model the behavior of the PCF at a small scale in the presence of an auto gravitating medium having a line/plane of symmetry in 2D/3D. The analysis of the PCF in an astrophysical context was split into two, adopting a non-Poissonian Voronoi Tessellation (NPVT). We first analyzed the case of a 2D cut which covers few voids and a 2D cut which covers approximately 50 voids. The obtained PCF in the case of many voids was then discussed in comparison to the bootstrap predictions for a PVT process and the observed PCF for an astronomical catalog. An approximated formula which connects the averaged radius of the cosmic voids to the first minimum of the PCF is given.

The UVOIR bolometric light curves are usually modeled by the radioactive decay. In order to model more precisely the absolute/apparent magnitude versus time relationship the continuous production of radioactive isotopes is introduced. A differential equation of the first order with separable variables is solved.

The gamma density function is usually defined in interval between zero and infinity. This paper introduces an upper and a lower boundary to this distribution. The parameters which characterize the truncated gamma distribution are evaluated. A statistical test is performed on two samples of stars. A comparison with the lognormal and the four power law distribution is made.

We review the model of the Voronoi Diagrams which allows to reproduce the large-scale structures of our universe as given by the astronomical catalogs. The observed number of galaxies in a given solid angle with a chosen flux/magnitude versus the redshift presents a maximum that is a function of the flux/magnitude ; it can be explained by a detailed analysis of the standard luminosity function for galaxies as well by two new luminosity function for galaxies. The current status of the research on the statistics of the Voronoi Diagrams is reviewed.

The introduction of a exponential or power law gradient in the interstellar medium (ISM) allows to produce an asymmetric evolution of the supernova remnant (SNR) when the framework of the thin layer approximation is adopted. Unfortunately both the exponential and power law gradients for the ISM do not have a well defined physical meaning. The physics conversely is well represented by an isothermal self-gravitating disk of particles whose velocity is everywhere Maxwellian. . We derived a law of motion in the framework of the thin layer approximation with a control parameter of the swept mass. The photon's losses ,that are often neglected in the thin layer approximation, are modeled trough a velocity dependence. The developed framework is applied to SNR 1987A and the three observed rings are simulated.

A method is developed to compute the chord length distribution along a line which intersects a cellular Universe. The cellular Universe is here modeled by the Poissonian Voronoi Tessellation (PVT) and by a non-Poissonian Voronoi Tessellation (NPVT). The distribution of the spheres is obtained from common approximations used in modeling the volumes of Voronoi Diagrams. We give analytical formulas for the distributions of the lengths of chords in both the PVT and NPVT. The astrophysical applications are made to the real Eso Slice Project and to an artificial slice of galaxies which simulates the 2dF Galaxy Redshift Survey.

The initial mass function (IMF) for the stars is usually fitted by three straight lines, which means seven parameters. The presence of brown dwarfs (BD) increases to four the straight lines and to nine the parameters. Another common fitting function is the lognormal distribution, which is characterized by two parameters. This paper is devoted to demonstrating the advantage of introducing a left truncated beta probability density function, which is characterized by four parameters. The constant of normalization, the mean, the mode and the distribution function are calculated for the left truncated beta distribution. The normal-beta (NB) distribution which results from convolving independent normally distributed and beta distributed components is also derived. The chi-square test and the K-S test are performed on a first sample of stars and BDs which belongs to the massive young cluster NGC 6611 and on a second sample which represents the star's masses of the cluster NGC 2362.

The release of the 2MASS Redshift Survey (2MRS), with it's 44599 galaxies allows the deduction of the galaxy's masses in nearly complete sample. A cubic box with side of 37 Mpc containing 2429 galaxies is extracted and the Newtonian gravitational field is evaluated both at the center of the box as well in 101 x 101 x 101 grid points of the box. The obtained results are then discussed at the light of the shell theorem which states that at the internal of a sphere the gravitational field is zero.

The Cosmological Principle in it's various versions states that: (i) the Galaxy does not occupy a particular position, (ii) the Universe is homogeneous and isotropic. This statement does not agree with the recent astronomical observations in the range z lower than 0.05 which are in agreement with a cellular structure of the Universe. Here we present a local analysis of the inhomogeneity of the Universe. When z is greater than 0.05 our analysis cannot be applied because the astronomical sample of galaxies here processed is not complete. The two tools of the Poisson Voronoi tessellation (PVT) and the luminosity function for galaxies allow of building a new version of the local Cosmological Principle.

In this chapter first the statistics of the standard and truncated Pareto distributions are derived and used to fit empirical values of asteroids diameters from different families, namely, Koronis, Eos and Themis, and from the Astorb database. A theoretical analysis is then carried out and two possible physical mechanisms are suggested that account for Pareto tails in distributions of asteroids diameter.

Stereology allows shifting from the 3D distribution of the volumes of Poissonian Voronoi Diagrams to their 2D cross-sections. The basic assumption is that the 3D statistics of the volumes of the voids in the local Universe has a distribution function of the gamma-type. The standard rule of conversion from 3D volumes to 2D circles, adopting the standard rules of stereology, produces a new probability density function of the radii which contains the Meijer $G$-function. A non-Poissonian distribution of volumes is also considered. The distribution of the 3D radii of the Sloan Digital Sky Survey Data Release 7 is best fitted by a non-Poissonian distribution in volumes as given by the Kiang function with argument of about two.

The image of planetary nebulae (PN), supernova remnant (SNR) and Eta-Carinae is made by three different physical processes. The first process is the expansion of the shell that can be modeled by the canonical laws of motion in the spherical case and by the momentum conservation when gradients of density are present in the interstellar medium. The quality of the simulations is introduced along one direction as well along many directions. The second process is the diffusion of particles that radiate from the advancing layer. The 3D diffusion from a sphere, the 1D diffusion with drift and 1D random walk are analyzed. The third process is the composition of the image through an integral operation along the line of sight. The developed framework is applied to three PN which are A39, the Ring nebula and the etched hourglass nebula MyCn 18, the hybrid object Eta-Carinae, and to two SNR which are SN 1993J and SN 1006. In all the considered cases a careful comparison between the observed and theoretical profiles in intensity is done.

The expansion of a superbubble is investigated both analytically and numerically. Our model implements the thin layer approximation in a vertical profile of density as given by an isothermal self-gravitating disk. A precise comparison with the results of numerical hydro-dynamics is given. Analogies are drawn with the Kompaneets equation that includes the quadratic hyperbolic-secant law in the list of the plane-parallel stratified media. An astrophysical application is made to the superbubble connected with the two worms 46.4+5.5 and 39.7+5.7. The effects of the rotation of the galaxy on the simulated radius and on the velocity are introduced. The worms with their strong limb-brightening visible on astronomical maps are explained in the framework of image theory.

The time evolution of a Gamma-ray burst (GRB) is associated with the evolution of a supernova remnant (SNR). The time evolution of the flux of a GRB is modeled introducing a law for the density of the medium in the advancing layer. The adopted radiative model for the GRB in the various e.m. frequencies is synchrotron emission. The X-ray ring which appears a few hours after a GRB is simulated.

The surface brightness diameter relationship for supernovae remnants (SNRs) is explained by adopting a model of direct conversion of the flux of kinetic energy into synchrotron luminosity. Two laws of motion are adopted, a power law model for the radius-time relationship, and a model which uses the thin layer approximation. The fluctuations on the log-log surface diameter relationship are modeled by a Monte Carlo simulation. In this model a new probability density function for the density as function of the galactic height is introduced.

A classical and a relativistic law of motion for an advancing shell are deduced applying the thin layer approximation. A new parameter connected with the quantity of absorbed matter in the expansion is introduced; this allows of matching theory and observation.

A new law of motion for supernova remnant (SNR) which introduces the quantity of swept matter in the thin layer approximation is introduced. This new law of motion is tested on 10 years observations of SN1993J. The introduction of an exponential gradient in the surrounding medium allows to model an aspherical expansion. A weakly asymmetric SNR, SN1006, and a strongly asymmetric SNR, SN1987a, are modeled. In the case of SN1987a the three observed rings are simulated.

A classical and a relativistic law of motion for a supernova remnant (SNR) are deduced assuming an inverse power law behavior for the density of the interstellar medium and applying the thin layer approximation. A third equation of motion is found in the framework of relativistic hydrodynamics with pressure, applying momentum conservation. These new formulas are calibrated against a decade of observations of \snr. The existing knowledge of the diffusive processes of ultrarelativistic electrons is reviewed in order to explain the behavior of the `U' shaped profile of intensity versus distance from the center of SN 1993J.

The voids between galaxies are identified with the volumes of the Poisson Voronoi tessellation. Two new survival functions for the apparent radii of voids are derived. The sectional normalized area of the Poisson Voronoi tessellation is modelled by the Kiang function and by the exponential function. Two new survival functions with equivalent sectional radius are therefore derived; they represent an alternative to the survival function of voids between galaxies as given by the self-similar distribution. The spatial appearance of slices of the 2dF Galaxy Redshift Survey is simulated.

Self-similar solutions to the problem of a special relativistic law of motion for thin shells of matter are calculated. These solutions represent the special relativistic generalization of momentum conservation for the thin layer approximation in classical physics. The analytical and numerical results are applied to Supernova Remnant 1987 A.

The asymmetric shape of the nebula around $\eta$-Carinae (Homunculus) can be explained by a spherical expansion in a non-homogeneous medium. Two models are analyzed: an exponential and an inverse power law dependence for the density as a function of distance from the equatorial plane. The presence of a medium with variable density along the polar direction progressively converts the original spherical shell into a bipolar nebula. In the case of the nebula around $\eta$-Carinae, we know the time elapsed since the great outburst in 1840. An exact match between observed radii and velocities can be obtained by fine tuning the parameters involved, such as initial radius, initial velocity and the typical scale that characterizes the gradient in density. The observed radius and velocity of the Homunculus as a function of the polar angle in spherical coordinates can be compared with the corresponding simulated data by introducing the efficiency in a single or multiple directions. Once the 3D spatial structure of the Homunculus is obtained, we can compose the image by integrating along the line of sight. In order to simulate the observed image, we have considered a bipolar nebula with constant thickness and an optically thin emitting layer. Some simulated cuts of the relative intensity are reported and may represent a useful reference for the astronomical cuts.

It is shown that the Hubble constant can be derived from the standard luminosity function of galaxies as well as from a new luminosity function as deduced from the mass-luminosity relationship for galaxies. An analytical expression for the Hubble constant can be found from the maximum number of galaxies (in a given solid angle and flux) as a function of the redshift. A second analytical definition of the Hubble constant can be found from the redshift averaged over a given solid angle and flux. The analysis of two luminosity functions for galaxies brings to four the new definitions of the Hubble constant. The equation that regulates the Malmquist bias for galaxies is derived and as a consequence it is possible to extract a complete sample. The application of these new formulae to the data of the two-degree Field Galaxy Redshift Survey provides a Hubble constant of $( 65.26 \pm 8.22 ) \mathrm{\ km\ s}^{-1}\mathrm{\ Mpc}^{-1}$ for a redshift lower than 0.042. All the results are deduced in a Euclidean universe because the concept of space-time curvature is not necessary as well as in a static universe because two mechanisms for the redshift of galaxies alternative to the Doppler effect are invoked.

Two new luminosity functions of galaxies can be built starting from three and four parameter generalized gamma distributions. In the astrophysical conversion, the number of parameters increases by one, due to the addition of the overall density of galaxies. A third new galaxy luminosity function is built starting from a three parameter generalized gamma distribution for the mass of galaxies once a simple nonlinear relationship between mass and luminosity is assumed; in this case the number of parameters is five because the overall density of galaxies and a parameter that regulates mass and luminosity are added. The three new galaxy luminosity functions were tested on the Sloan Digital Sky Survey (SDSS) in five different bands; the results always produce a "better fit" than the Schechter function. The formalism that has been developed allows to analyze the Schechter function with a transformation of location. A test between theoretical and observed number of galaxies as a function of redshift was done on data extracted from a two-degree field galaxy redshift survey.

The 3D network originated by the faces of irregular Poissonian Voronoi Polyhedrons may represent the backbone on which the galaxies are originated. As a consequence the spatial appearance of the catalogs of galaxies can be reproduced. The selected catalogs to simulate are the 2dF Galaxy Redshift Survey and the Third Reference Catalog of Bright Galaxies. In order to explain the number of observed galaxies for a given flux/magnitude as a function of the redshift, the photometric properties of the galaxies should be carefully examined from both the astronomical and theoretical point of view. The statistics of the Voronoi normalized volume is modeled by two distributions and the Eridanus super-void is identified as the largest volume belonging to the Voronoi Polyhedron. The behavior of the correlation function for galaxies is simulated by adopting the framework of thick faces of Voronoi Polyhedrons on short scales, while adopting standard arguments on large scales.

The publication of new tables of calibration of some fundamental parameters of Galactic B0-B5 supergiants in the two classes $Ia$ and $Ib$ allow to particularize the eight parameters conjecture that model five fundamental parameters. The numerical expression for visual magnitude, radius, mass, luminosity and surface gravity are derived for supergiants in the range of temperature between 29700 and 15200. The presence of accurate tables of calibration allows us to introduce the efficiency of the derived formulas.

The overall properties of the Herbig-Haro objects such as centerline velocity, transversal profile of velocity, flow of mass and energy are explained adopting two models for the turbulent jet. The complex shapes of the Herbig-Haro objects, such as the arc in HH34 can be explained introducing the combination of different kinematic effects such as velocity behavior along the main direction of the jet and the velocity of the star in the interstellar medium. The behavior of the intensity or brightness of the line of emission is explored in three different cases : transversal 1D cut, longitudinal 1D cut and 2D map. An analytical explanation for the enhancement in intensity or brightness such as usually modeled by the bow shock is given by a careful analysis of the geometrical properties of the torus.

The image of planetary nebulae is made by three different physical processes. The first process is the expansion of the shell that can be modeled by the canonical laws of motion in the spherical case and by the momentum conservation when gradients of density are present in the interstellar medium. The second process is the diffusion of particles that radiate from the advancing layer. The 3D diffusion from a sphere as well as the 1D diffusion with drift are analyzed. The third process is the composition of the image through an integral operation along the line of sight. The developed framework is applied to A39, to the Ring nebula and to the etched hourglass nebula MyCn 18.

The turbulent jets are usually described by classical velocities. The relativistic case can be treated starting from the conservation of the relativistic momentum. The two key assumptions which allow to obtain a simple expression for the relativistic trajectory and relativistic velocity are null pressure and constant density.

This paper summarizes the known physics of turbulent jets observed in laboratory experiments. The formula, which gives the power released in turbulence describes the concentration of turbulence/relativistic particles in each point of the astrophysical jets. The same expression is also used to analyze the power released in turbulence in the case of pipe and non Newtonian fluids. Through an integral operation it is possible to deduce the intensity of synchrotron radiation for a profile perpendicular or not to a straight jet, a 2D map for a perpendicular, randomly oriented straight jet as well as a 2D map of complex trajectories such as NCC4061 and 3C31. Presented here is a simulation of the spectral index in brightness of 3C273 as well as a 2D map of the degree of linear polarization. The Sobel operator is applied to the theoretical 2D maps of straight perpendicular jets.

The absolute visual magnitude as function of the observed colour (B-V), also named Hertzsprung-Russell diagram can be described through five equations; that in presence of calibrated stars means eight constants. The developed framework allows to deduce the remaining physical parameters that are mass, radius and luminosity. This new technique is applied to the first 10 pc, the first 50 pc, the Hyades and to the determination of the distance of a cluster. The case of the white dwarfs is analysed assuming the absence of calibrated data: our equation produces a smaller $\chi^2$ in respect to the standard colour-magnitude calibration when applied to the Villanova Catalog of Spectroscopically Identified White Dwarfs. The theoretical basis of the formulae for the colours and the bolometric correction of the stars are clarified through a Taylor expansion in the temperature of the Planck distribution.

A new luminosity function for galaxies can be built starting from the product of two random variables X and Y represented by a gamma variate with argument 2 >. The mean, the standard deviation and the distribution function of this new distribution are computed. This new probability density function is assumed to describe the mass distribution of galaxies. Through a non linear rule of conversion from mass to luminosity a second new luminosity function for galaxies is derived. The test of reliability of these two luminosity functions was made on the Sloan Digital Sky Survey (SDSS) in five different bands. The Schechter function gives a better fit with respect to the two new luminosity functions for galaxies here derived.

The Pareto probability distribution is widely applied in different fields such us finance, physics, hydrology, geology and astronomy. This note deals with an application of the Pareto distribution to astrophysics and more precisely to the statistical analysis of mass of stars and of diameters of asteroids. In particular a comparison between the usual Pareto distribution and its truncated version is presented. Finally a possible physical mechanism that produces Pareto tails for the distribution of the masses of stars is suggested.

The search for a luminosity function for galaxies both alternative or companion to a Schechter function is a key problem in the reduction of data from catalogs of galaxies. Two luminosity functions for galaxies can be built starting from two distributions of mass as given by the fragmentation. A first overall distribution function is the Kiang function that represents a useful description of the area and volume distribution of the Poisson Voronoi diagrams. The second distribution, that covers the case of low mass galaxies, is the truncated Pareto distribution : in this model we have a natural bound due to the minimum mass/luminosity observed and an upper bound (function of the considered environment) represented by the boundary with the observed mass/luminosity overall behaviour. The mass distribution is then converted into a luminosity distribution through a standard mass-luminosity relationship.

The non-thermal X-rays from the SN 1006 NE rim present characteristic scale lengths that are interpreted in the context of diffusion of a relativistic electron. The adopted theoretical framework is the mathematical diffusion in 3D, 1D and 1D with drift as well as the Monte Carlo random walk in 1D with drift. The asymmetric random walk with diffusion from a plane can explain the scale widths of 0.04 pc upstream and 0.2 pc downstream in the non thermal intensity of X-ray emission in SN 1006. A mathematical image of the non thermal X-flux from an supernova remnant as well as profiles function of the distance from the center can be simulated. This model provides a reasonable description of both the limbs and the central region of SN 1006. A new method to deduce the magnetic field in supernova remnant is suggested.

Super-bubbles are shells in the interstellar medium produced by the simultaneous explosions of many supernova remnants. The solutions of the mathematical diffusion and of the Fourier expansion in 1D, 2D and 3D were deduced in order to describe the diffusion of nucleons from such structures. The mean number of visits in the the case of the Levy flights in 1D was computed with a Monte Carlo simulation. The diffusion of cosmic rays has its physical explanation in the relativistic Larmor gyro-radius which is energy dependent. The mathematical solution of the diffusion equation in 1D with variable diffusion coefficient was computed. Variable diffusion coefficient means magnetic field variable with the altitude from the Galactic plane. The analytical solutions allow us to calibrate the code that describes the Monte Carlo diffusion. The maximum energy that can be extracted from the super-bubbles is deduced. The concentration of cosmic rays is a function of the distance from the nearest super-bubble and the selected energy. The interaction of the cosmic rays on the target material allows us to trace the theoretical map of the diffuse Galactic continuum gamma-rays. The streaming of the cosmic rays from the Gould Belts that contains the sun at it's internal was described by a Monte Carlo simulation. Nine new formulas are derived.

In order to find a law characterising the decrease of velocity along a jet, five analytical methods are suggested. The first two simple models examine the variation of velocity in the presence of Newton's or Stoke's resistance. The equation that represents the conservation of the momentum along a pyramidal sector is solved from an analytical point of view (third model). The application of the conservation of the total momentum flux allows us to deduce the velocity of the galaxy as a function of time for classical velocities (fourth model) and relativistic velocities (fifth model). The variation of velocity along the jet combined with an adequate composition of jet precession velocity, rotational velocity of the galaxy, and galaxy dispersion velocity in the cluster allows us to trace the geometrical pattern of the head-tail radio sources. Application of the developed theory/code to the radio galaxies NGC1265,NGC4061,NGC326, and Cygnus A gives the central galaxy's approximate dispersion velocities in the direction perpendicular to the jet. A transition from head--tails to classical double radio galaxies as a function of the increasing jet's mechanical power is introduced.

The size distributions of 2D and 3D Voronoi cells and of cells of $V_p(2,3)$,--2D cut of 3D Voronoi diagram--are explored, with the single-parameter (re-scaled) gamma distribution playing a central role in the analytical fitting. Observational evidence for a cellular universe is briefly reviewed. A simulated $V_p(2,3)$ map with galaxies lying on the cell boundaries is constructed to compare, as regards general appearance, with the observed CFA map of galaxies and voids, the parameters of the simulation being so chosen as to reproduce the largest observed void size.

The dichotomy between a universal mass function (IMF) and a variable IMF which depends on local physical parameters characterises observational and theoretical stellar astronomy. In this contribution the available distributions of probability are briefly reviewed. The physical nature of two of them, gamma variate and lognormal, is then explained once the framework of the fragmentation is introduced. Interpolating techniques are then applied to the sample of the first 10 pc and to the open cluster NGC6649: in both cases lognormal distribution produces the best fit. The three power law function has also been investigated and visual comparison with an artificially generated sample of 100000 stars suggests that the variations in the spectral index are simply due to the small number of stars available in the observational sample. In order to derive the sample of masses, a new formula that allows us to express the mass as a function of the absolute magnitude and (B-V) for MAIN V, GIANTS III and SUPERGIANTS I is derived.

The galactic supershells are cavities in the interstellar medium. These shells can be explained by introducing the concept of superbubbles, the theoretical result of multiple supernova. The superbubbles can be analytically described if the ambient medium has a constant density both in the so-called bursting phase and in the subsequent adiabatic expansion. In order to solve the expansion of superbubbles in the ISM, which is a non-homogeneous medium, a numerical technique is used that divides the sphere into many sectors. By varying the time of the bursting phenomenon and the time over which the phenomenon is followed, elliptical and hour-glass shapes or vertical walls can be obtained. Application of the developed theory/code to the super-shell associated with GW~46.4+5 and with GSH 238 allows us to say that the suggested physical parameters are consistent with our theory. The map of the expanding superbubble's velocity can be tentatively traced by generating random points on the expanding surface. The structure of the galactic plane as a result of the evolution of many super-bubbles was simulated by adopting the percolation theory in order to generate new OB associations.