We present a code allowing to evolve three-dimensional self-gravitating collisionless systems with a large number of particles N >= 10^7. FLY (Fast Level-based N-bodY code) is a fully parallel code based on a tree algorithm. It adopts periodic boundary conditions implemented by means of the Ewald summation technique. FLY is based on the one-side communication paradigm and was originally developed on a CRAY T3E system using the SHMEM library and it was ported on SGI ORIGIN 2000 and on IBM SP. FLY (http://www.ct.astro.it/fly/) version 1.1 is an open source freely available code. FLY data output can be analysed with AstroMD, an analysis and visualization tool specifically designed to deal with the visualization and analysis of astrophysical data. AstroMD can manage different physical quantities. It can find out structures without well defined shape or symmetries, and perform quantitative calculations on selected regions. AstroMD (http://www.cineca.it/astromd) is a freely available code.

Over the past few years, the role of visualization for scientific purpose has grown up enormously. Astronomy makes an extended use of visualization techniques to analyze data, and scientific visualization has became a fundamental part of modern researches in Astronomy. With the evolution of high performance computers, numerical simulations have assumed a great role in the scientific investigation, allowing the user to run simulation with higher and higher resolution. Data produced in these simulations are often multi-dimensional arrays with several physical quantities. These data are very hard to manage and to analyze efficiently. Consequently the data analysis and visualization tools must follow the new requirements of the research. AstroMD is a tool for data analysis and visualization of astrophysical data and can manage different physical quantities and multi-dimensional data sets. The tool uses virtual reality techniques by which the user has the impression of travelling through a computer-based multi-dimensional model. AstroMD is a freely available tool for the whole astronomical community (http://www.cineca.it/astromd/).

In this paper we present a critical discussion of the algorithms commonly used in N-body simulations of Galaxy Formation to deal with the energy equation governing heating and cooling, to model star formation and the star formation rate, and to account for energy feed-back from stars. First, we propose our technique for solving the energy equation in presence of heating and cooling, which includes some difference with respect to the standard semi-implicit technique. Second, we examine the current criteria for the onset of the star formation activity. We suggest a new approach, in which star formaiton is let depend on the total mass density - baryonic (gas and stars) and dark matter - of the system and on the metal-dependent cooling efficiency. Third, we check and discuss the separate effects of energy (and mass) feed-back from several sources - namely supernovae, stellar winds from massive stars, and UV flux from the same objects. All the simulations are performed in the framework of the formation and evolution of a disk galaxy. We show that the inclusion of these physical phenomena has a signigicant impact on the evolution of the galaxy model.

We present preliminary results on galactic Dark Matter (DM), halo structure, and galaxy evolution. We show how during the first Gyr of the evolution of a $10^{10} M_{\odot}$ dwarf elliptica feed-back from stars (SN\ae~ and stellar winds) leads to an extended constant density isothermal core with radius of 0.15 the virial radius $R_{200}$. We also present first results on galaxy merging as a possibile scenario to form ellipticals, studying in particular how the details of the merging evolution vary as a function of the mass ratio of the interacting galaxies.