Dynamical and radiative properties of astrophysical supersonic jets I. Cocoon morphologies

We present the results of a numerical analysis of the propagation and interaction of a supersonic jet with the external medium. We discuss the motion of the head of the jet into the ambient in different physical conditions, carrying out calculations with different Mach numbers and density ratios of the jet to the exteriors. Performing the calculation in a reference frame in motion with the jet head, we can follow in detail its long term dynamics. This numerical scheme allows us also to study the morphology of the cocoon for different physical parameters. We find that the propagation velocity of the jet head into the ambient medium strongly influences the morphology of the cocoon, and this result can be relevant in connection to the origin and structure of lobes in extragalactic radiosources.Comment: 14 pages, TeX. Accepted for A&

Astrophysical fluid simulations of thermally ideal gases with non-constant adiabatic index: numerical implementation

An Equation of State (\textit{EoS}) closes the set of fluid equations. Although an ideal EoS with a constant \textit{adiabatic index} $\Gamma$ is the preferred choice due to its simplistic implementation, many astrophysical fluid simulations may benefit from a more sophisticated treatment that can account for diverse chemical processes. Here, we first review the basic thermodynamic principles of a gas mixture in terms of its thermal and caloric EoS by including effects like ionization, dissociation as well as temperature dependent degrees of freedom such as molecular vibrations and rotations. The formulation is revisited in the context of plasmas that are either in equilibrium conditions (local thermodynamic- or collisional excitation- equilibria) or described by non-equilibrium chemistry coupled to optically thin radiative cooling. We then present a numerical implementation of thermally ideal gases obeying a more general caloric EoS with non-constant adiabatic index in Godunov-type numerical schemes.We discuss the necessary modifications to the Riemann solver and to the conversion between total energy and pressure (or vice-versa) routinely invoked in Godunov-type schemes. We then present two different approaches for computing the EoS.The first one employs root-finder methods and it is best suited for EoS in analytical form. The second one leans on lookup table and interpolation and results in a more computationally efficient approach although care must be taken to ensure thermodynamic consistency. A number of selected benchmarks demonstrate that the employment of a non-ideal EoS can lead to important differences in the solution when the temperature range is $500-10^4$ K where dissociation and ionization occur. The implementation of selected EoS introduces additional computational costs although using lookup table methods can significantly reduce the overhead by a factor $3\sim 4$.Comment: 17 pages, 10 figures, Accepted for publication in A&

Selective large-eddy simulation of hypersonic flows. Procedure to activate the filtering in unresolved regions only (arXiv:1211.1305, submitted to Computer Physics Communications)

A new method for the localization of the regions where the turbulent fluctuations are unresolved is applied to the large-eddy simulation (LES) of a compressible turbulent jet with an initial Mach number equal to 5. The localization method used is called selective LES and is based on the exploitation of a scalar probe function f which represents the magnitude of the stretching-tilting term of the vorticity equation normalized with the enstrophy (Tordella et al. 2007). For a fully developed turbulent field of fluctuations, statistical analysis shows that the probability that f is larger than 2 is almost zero, and, for any given threshold, it is larger if the flow is under-resolved. By computing the spatial field of f in each instantaneous realization of the simulation it is possible to locate the regions where the magnitude of the normalized vortical stretching-tilting is anomalously high. The sub-grid model is then introduced into the governing equations in such regions only. The results of the selective LES simulation are compared with those of a standard LES, where the sub-grid terms are used in the whole domain. The comparison is carried out by assuming as reference field a higher resolution Euler simulation of the same jet. It is shown that the selective LES modifies the dynamic properties of the flow to a lesser extent with respect to the classical LES. In particular, the prediction of the enstrophy distribution and of the energy and density spectra are substantially improve

Making Fanaroff-Riley I radio sources. Numerical Hydrodynamic 3D Simulations of Low Power Jets

Extragalactic radio sources have been classified into two classes, Fanaroff-Riley I and II, which differ in morphology and radio power. Strongly emitting sources belong to the edge-brightened FR II class, and weakly emitting sources to the edge-darkened FR I class. The origin of this dichotomy is not yet fully understood. Numerical simulations are successful in generating FR II morphologies, but they fail to reproduce the diffuse structure of FR Is. By means of hydro-dynamical 3D simulations of supersonic jets, we investigate how the displayed morphologies depend on the jet parameters. Bow shocks and Mach disks at the jet head, which are probably responsible for the hot spots in the FR II sources, disappear for a jet kinetic power L_kin < 10^43 erg/s. This threshold compares favorably with the luminosity at which the FR I/FR II transition is observed. The problem is addressed by numerical means carrying out 3D HD simulations of supersonic jets that propagate in a non-homogeneous medium with the ambient temperature that increases with distance from the jet origin, which maintains constant pressure. The jet energy in the lower power sources, instead of being deposited at the terminal shock, is gradually dissipated by the turbulence. The jets spread out while propagating, and they smoothly decelerate while mixing with the ambient medium and produce the plumes characteristic of FR I objects. Three-dimensionality is an essential ingredient to explore the FR I evolution because the properties of turbulence in two and three dimensions are very different, since there is no energy cascade to small scales in two dimensions, and two-dimensional simulations with the same parameters lead to FRII-like behavior.Comment: 11 pages, 12 figures, to appear on A&

Diffusive shock acceleration in extragalactic jets

We calculate the temporal evolution of distributions of relativistic electrons subject to synchrotron and adiabatic processes and Fermi-like acceleration in shocks. The shocks result from Kelvin-Helmholtz instabilities in the jet. Shock formation and particle acceleration are treated in a self-consistent way by means of a numerical hydrocode. We show that in our model the number of relativistic particles is conserved during the evolution, with no need of further injections of supra-thermal particles after the initial one. From our calculations, we derive predictions for values and trends of quantities like the spectral index and the cutoff frequency that can be compared with observations.Comment: 12 pages containing 7 postscript figures; uses A&A macros. Accepted for publication in Astronomy and Astrophysic

Dynamical and radiative properties of astrophysical supersonic jets

We present the results of a numerical analysis of the propagation and interaction of a supersonic jet with the external medium. We discuss the motion of the head of the jet into the ambient in different physical conditions, carrying out calculations with different Mach numbers and density ratios of the jet to the exteriors. Performing the calculation in a reference frame in motion with the jet head, we can follow in detail its long term dynamics. This numerical scheme allows us also to study the morphology of the cocoon for different physical parameters. We find that the propagation velocity of the jet head into the ambient medium strongly influences the morphology of the cocoon, and this result can be relevant in connection to the origin and structure of lobes in extragalactic radiosources

Consumer attitudes and preference exploration towards fresh-cut salads using best–worst scaling and latent class analysis

This research explored the preferences and buying habits of a sample of 620 consumers of fresh-cut, ready-to-eat salads. A best–worst scaling approach was used to measure the level of preference stated by individuals regarding 12 attributes for quality (intrinsic, extrinsic and credence) of fresh-cut salads. The experiment was carried out through direct interviews at several large-scale retail outlets in the Turin metropolitan area (north-west of Italy). Out of the total number of questioned consumers, 35% said they did not consume fresh-cut salads. On the contrary, the rest of the involved sample expressed the highest degree of preference towards the freshness/appearance attribute, followed by the expiration date and the brand. On the contrary, attributes such as price, organic certification and food safety did not emerge as discriminating factors in consumer choices. Additionally, five clusters of consumers were identified, whose preferences are related both to purchasing styles and socio-demographic variables. In conclusion, this research has highlighted the positive attitude of consumers towards quality products backed by a brand, providing ideas for companies to improve within this sector and implement strategies to answer the needs of a new segment of consumers, by determining market opportunities that aim to strengthen local brands

YSO jets: MHD simulations with radiative cooling

The new High Power Computing facilities available for the scientific community allowed the use of increasingly complex codes for the numerical simulation of physical processes. Recent magnetohydrodynamic (MHD) simulations of astrophysical jets could finally include non-ideal effects, such as the radiative cooling we will discuss in this work in the context of protostellar jets. This makes the simulations more reliable and, with the recent improvements in available observational data, will provide a valuable tool for model discrimination. From 2D adaptivemesh refinement (AMR) simulations, synthetic surface brightness maps for the line emissions are computed, to be compared with observations