We describe a new numerical 3D relativistic hydrodynamical code, the results of validation tests, and a comparison with earlier, 2D studies. The 3D code has been used to study the deflection and precession of relativistic flows. We find that even quite fast jets (gamma~10) can be significantly influenced by impinging on an oblique density gradient, exhibiting a rotation of the Mach disk in the jet's head. The flow is bent via a potentially strong, oblique internal shock that arises due to asymmetric perturbation of the flow by its cocoon. In extreme cases this cocoon can form a marginally relativistic flow orthogonal to the jet, leading to large scale dynamics quite unlike that normally associated with astrophysical jets. Exploration of a gamma=5 flow subject to a large amplitude precession (semi-angle 11.25dg) shows that it retains its integrity, with modest reduction in Lorentz factor and momentum flux, for almost 50 jet-radii, but thereafter, the collimated flow is disrupted. The flow is approximately ballistic, with velocity vectors not aligned with the local jet `wall'. We consider simple estimators of the flow emissivity in each case and conclude that a) while the oblique internal shocks which mediate a small change in the direction of the deflected flows have little impact on the global dynamics, significantly enhanced flow emission (by a factor of 2-3) may be associated with such regions; and b) the convolution of rest frame emissivity and Doppler boost in the case of the precessed jet invariably leads to a core-jet-like structure, but that intensity fluctuations in the jet cannot be uniquely associated with either change in internal conditions or Doppler boost alone, but in general are a combination of both factors.

The amplitude of jet distortions and accompanying pressure and velocity fluctuations resulting from Kelvin-Helmholtz instability of three dimensional relativistic jets are explored. The effect of instability on jets as they accelerate from sub- to super-Alfv\'enic speeds is explored and a relativistic stabilization mechanism for trans-Alfv\'enic jets is proposed. The level to which asymmetric instabilities on supermagnetosonic relativistic jets will grow is predicted theoretically and a Doppler boosted ``apparent'' emissivity is computed. Effects due to helically twisted filamentary structure produced by asymmetric modes of instability should be readily observable on relativistic jets.

We compare results from a relativistic and a nonrelativistic set of 2D axisymmetric jet simulations. For a set of five relativistic simulations that either increase the Lorentz factor or decrease the adiabatic index we compute nonrelativistic simulations with equal useful power or thrust. We examine these simulations for morphological and dynamical differences, focusing on the velocity field, the width of the cocoon, the age of the jets, and the internal structure of the jet itself. The primary result of these comparisons is that the velocity field of nonrelativistic jet simulations cannot be scaled up to give the spatial distribution of Lorentz factors seen in relativistic simulations. Since the local Lorentz factor plays a major role in determining the total intensity for parsec scale extragalactic jets, this suggests that a nonrelativistic simulation cannot yield the proper intensity distribution for a relativistic jet. Another general result is that each relativistic jet and its nonrelativistic equivalents have similar ages (in dynamical time units, = R/a_a, where R is the initial radius of a cylindrical jet and a_a is the sound speed in the ambient medium). In addition to these comparisons, we have completed four new relativistic simulations to investigate the effect of varying thermal pressure on relativistic jets. The simulations generally confirm that faster (larger Lorentz factor) and colder jets are more stable, with smaller amplitude and longer wavelength internal variations. The apparent stability of these jets does not follow from linear normal mode analysis, which suggests that there are available growing Kelvin-Helmholtz modes. (Abridged.)

A series of simulated maps showing the appearance in total intensity of flows computed using a recently developed relativistic hydrodynamic code (Duncan \& Hughes 1994: ApJ, 436, L119) are presented. The radiation transfer calculations were performed by assuming the flow is permeated by a magnetic field and fast particle distribution in energy equipartition, with energy density proportional to the hydrodynamic energy density (i.e., pressure). We find that relativistic flows subject to strong perturbations exhibit a density structure consisting of a series of nested bow shocks, and that this structure is evident in the intensity maps for large viewing angles. However, for viewing angles $<30^{\circ}$, differential Doppler boosting leads to a series of axial knots of emission, similar to the pattern exhibited by many VLBI sources. The appearance of VLBI knots is determined primarily by the Doppler boosting of parts of a more extended flow. To study the evolution of a perturbed jet, a time series of maps was produced and an integrated flux light curve created. The light curve shows features characteristic of a radio loud AGN: small amplitude variations and a large outburst. We find that in the absence of perturbations, jets with a modest Lorentz factor ($\sim 5$) exhibit complex intensity maps, while faster jets (Lorentz factor $\sim 10$) are largely featureless. We also study the appearance of kiloparsec jet-counterjet pairs by producing simulated maps at relatively large viewing angles; we conclude that observed hot spot emission is more likely to be associated with the Mach disk than with the outer, bow shock.

We describe a method for the numerical solution of the relativistic Euler equations which we have found to be both robust and efficient, and which has enabled us to simulate relativistic jets. The technique employs a solver of the Godunov-type, with approximate solution of the local Riemann problems, applied to laboratory frame variables. Lorentz transformations provide the rest frame quantities needed for the estimation of wave speeds, etc. This is applied within the framework of an adaptive mesh refinement algorithm, allowing us to perform high-resolution, 2-D simulations with modest computing resources. We present the results of nonrelativistic, and relativistic ($\gamma=5$ and $10$) runs, for adiabatic indices of $5/3$ and $4/3$. We find the same gross morphology in all cases, but the relativistic runs exhibit little instability and less well-defined structure internal to the jet: this might explain the difference between (relatively slow) BL~Lacs and (faster) QSOs. We find that the choice of adiabatic index makes a small but discernible difference to the structure of the shocked jet and ambient media.