We present the results of a three dimensional, locally isothermal, non-self-gravitating SPH code which models protoplanetary disks with two fluids: gas and dust. We ran simulations of a 1 Msun star surrounded by a 0.01 Msun disk comprising 99% gas and 1% dust in mass and extending from 0.5 to ~300 AU. The grain size ranges from 0.001 mm to 10 m for the low resolution (~25 000 SPH particles) simulations and from 0.1 mm to 10 cm for the high resolution (~160 000 SPH particles) simulations. Dust grains are slowed down by the sub-Keplerian gas and lose angular momentum, forcing them to migrate towards the central star and settle to the midplane. The gas drag efficiency varies according to the grain size, with the larger bodies being weakly influenced and following marginally perturbed Keplerian orbits, while smaller grains are strongly coupled to the gas. For intermediate sized grains, the drag force decouples the dust and gas, allowing the dust to preferentially migrate radially and efficiently settle to the midplane. The resulting dust distributions for each grain size will indicate, when grain growth is added, the regions when planets are likely to form.

We present the results of three-dimensional numerical simulations that include the effects of hydrodynamical forces and gas drag upon an evolving dusty gas disk. We briefly describe a new parallel, two phase numerical code based upon the smoothed particle hydrodynamics (SPH) technique in which the gas and dust phases are represented by two distinct types of particles. We use the code to follow the dynamical evolution of a population of grains in a gaseous protoplanetary disk in order to understand the distribution of grains of different sizes within the disk. Our ``grains'' range from metre to submillimetre in size.

We describe a new 512-CPU Beowulf cluster with Teraflop performance dedicated to problems in computational astrophysics. The cluster incorporates a cubic network topology based on inexpensive commodity 24-port gigabit switches and point to point connections through the second gigabit port on each Linux server. This configuration has network performance competitive with more expensive cluster configurations and is scaleable to much larger systems using other network topologies. Networking represents only about 9% of our total system cost of USD$561K. The standard Top 500 HPL Linpack benchmark rating is 1.202 Teraflops on 512 CPUs so computing costs by this measure are $0.47/Megaflop. We also describe 4 different astrophysical applications using complex parallel algorithms for studying large-scale structure formation, galaxy dynamics, magnetohydrodynamic flows onto blackholes and planet formation currently running on the cluster and achieving high parallel performance. The MHD code achieved a sustained speed of 2.2 teraflops in single precision or 44% of the theoretical peak.

We describe a parallel, cosmological N-body code based on a hybrid scheme using the particle-mesh (PM) and Barnes-Hut (BH) oct-tree algorithm. We call the algorithm GOTPM for Grid-of-Oct-Trees-Particle-Mesh. The code is parallelized using the Message Passing Interface (MPI) library and is optimized to run on Beowulf clusters as well as symmetric multi-processors. The gravitational potential is determined on a mesh using a standard PM method with particle forces determined through interpolation. The softened PM force is corrected for short range interactions using a grid of localized BH trees throughout the entire simulation volume in a completely analogous way to P$^3$M methods. This method makes no assumptions about the local density for short range force corrections and so is consistent with the results of the P$^3$M method in the limit that the treecode opening angle parameter, $\theta \to 0$. (abridged)

We have developed a new numerical technique for simulating dusty-gas flows. Our unique code incorporates gas hydrodynamics, self-gravity and dust drag to follow the dynamical evolution of a dusty-gas medium. We have incorporated several descriptions for the drag between gas and dust phases and can model flows with submillimetre, centimetre and metre size "dust". We present calculations run on the APAC* supercomputer following the evolution of the dust distribution in the pre-solar nebula. *Australian Partnership for Advanced Computing