A Multi-Epoch HST Study of the Herbig-Haro Flow from XZ Tauri

We present nine epochs of Hubble Space Telescope optical imaging of the bipolar outflow from the pre-main sequence binary XZ Tauri. Our data monitors the system from 1995-2005 and includes emission line images of the flow. The northern lobe appears to be a succession of bubbles, the outermost of which expanded ballistically from 1995-1999 but in 2000 began to deform and decelerate along its forward edge. It reached an extent of 6" from the binary in 2005. A larger and fainter southern counterbubble was detected for the first time in deep ACS images from 2004. Traces of shocked emission are seen as far as 20" south of the binary. The bubble emission nebulosity has a low excitation overall, as traced by the [S II]/H-alpha line ratio, requiring a nearly comoving surrounding medium that has been accelerated by previous ejections or stellar winds. Within the broad bubbles there are compact emission knots whose alignments and proper motions indicate that collimated jets are ejected from each binary component. The jet from the southern component, XZ Tau A, is aligned with the outflow axis of the bubbles and has tangential knot velocities of 70-200 km/s. Knots in the northern flow are seen to slow and brighten as they approach the forward edge of the outermost bubble. The knots in the jet from the other star, XZ Tau B, have lower velocities of ~100 km/s

Turbulence Model Implementation and Verification in the SENSEI CFD Code

This paper outlines the implementation and verification of the negative Spalart-Allmaras turbulence model into the SENSEI CFD code. The SA-neg turbulence model is implemented in a flexible, object-oriented framework where additional turbulence models can be easily added. In addition to outlining the new turbulence modeling framework in SENSEI, an overview of the other general improvements to SENSEI is provided. The results for four 2D test cases are compared to results from CFL3D and FUN3D to verify that the turbulence models are implemented properly. Several differences in the results from SENSEI, CFL3D, and FUN3D are identified and are attributed to differences in the implementation and discretization order of the boundary conditions as well as the order of discretization of the turbulence model. When a solid surface is located near or intersects an inflow or outflow boundary, higher order boundary conditions should be used to limit their effect on the forces on the surface. When the turbulence equations are discretized using second order spatial accuracy, the edge of the eddy viscosity profile seems to be sharper than when a first order discretization is used. However, the discretization order of the turbulence equation does not have a significant impact on output quantities of interest, such as pressure and viscous drag, for the cases studied

Measuring Gravitational Lensing Flexions in Abell 1689 Using an Analytic Image Model

Measuring dark matter substructure within galaxy cluster haloes is a fundamental probe of the Lambda-CDM model of structure formation. Gravitational lensing is a technique for measuring the total mass distribution which is independent of the nature of the gravitating matter, making it a vital tool for studying these dark-matter dominated objects. We present a new method for measuring weak gravitational lensing flexions, the gradients of the lensing shear field, to measure mass distributions on small angular scales. While previously published methods for measuring flexions focus on measuring derived properties of the lensed images, such as shapelet coefficients or surface brightness moments, our method instead fits a mass-sheet-transformation-invariant Analytic Image Model (AIM) to the each galaxy image. This simple parametric model traces the distortion of lensed image isophotes and constrains the flexion fields. We test the AIM method using simulated data images with realistic noise and a variety of unlensed image properties, and show that it successfully reproduces the input flexion fields. We also apply the AIM method for flexion measurement to Hubble Space Telescope observations of Abell 1689, and detect mass structure in the cluster using flexions measured with the AIM method.Comment: 44 pages, 4 figures, 3 tables. Accepted to ApJ. V2 (published version) has minor changes from V1; ApJ 736 (2011

HST and Spitzer Observations of the HD 207129 Debris Ring

A debris ring around the star HD 207129 (G0V; d = 16.0 pc) has been imaged in scattered visible light with the ACS coronagraph on the Hubble Space Telescope and in thermal emission using MIPS on the Spitzer Space Telescope at 70 microns (resolved) and 160 microns (unresolved). Spitzer IRS (7-35 microns) and MIPS (55-90 microns) spectrographs measured disk emission at >28 microns. In the HST image the disk appears as a ~30 AU wide ring with a mean radius of ~163 AU and is inclined by 60 degrees from pole-on. At 70 microns it appears partially resolved and is elongated in the same direction and with nearly the same size as seen with HST in scattered light. At 0.6 microns the ring shows no significant brightness asymmetry, implying little or no forward scattering by its constituent dust. With a mean surface brightness of V=23.7 mag per square arcsec, it is the faintest disk imaged to date in scattered light.Comment: 28 pages, 8 figure

NICMOS Imaging of a Damped Lyman-alpha Absorber at z=1.89 toward LBQS 1210+1731 : Constraints on Size and Star Formation Rate

We report results of a high-resolution imaging search (in rest frame H-$\alpha$ and optical continuum) for the galaxy associated with the damped Lyman-$\alpha$ (DLA) absorber at $z=1.892$ toward the $z_{em}=2.543$ quasar LBQS 1210+1731, using HST/NICMOS. After PSF subtraction, a feature is seen in both the broad-band and narrow-band images, at a projected separation of 0.25\arcsec from the quasar. If associated with the DLA, the object would be $\approx 2-3$ $h_{70}^{-1}$ kpc in size with a flux of $9.8 \pm 2.4$ $\mu$Jy in the F160W filter, implying a luminosity at $\lambda_{central}=5500$ {\AA} in the rest frame of $1.5 \times 10^{10}$ $h_{70}^{-2}$ L$_{\odot}$ at $z=1.89$, for $q_{0}=0.5$. However, no significant H-$\alpha$ emission is seen, suggesting a low star formation rate (SFR) (3 $\sigma$ upper limit of 4.0 $h_{70}^{-2}$ M$_{\odot}$ yr$^{-1}$), or very high dust obscuration. Alternatively, the object may be associated with the host galaxy of the quasar. H-band images obtained with the NICMOS camera 2 coronagraph show a much fainter structure $\approx 4-5$ $h_{70}^{-1}$ kpc in size and containing four knots of continuum emission, located 0.7\arcsec away from the quasar. We have probed regions far closer to the quasar sight-line than in most previous studies of high-redshift intervening DLAs. The two objects we report mark the closest detected high-redshift DLA candidates yet to any quasar sight line. If the features in our images are associated with the DLA, they suggest faint, compact, somewhat clumpy objects rather than large, well-formed proto-galactic disks or spheroids.Comment: 52 pages of text, 19 figures, To be published in Astrophysical Journal (accepted Dec. 8, 1999

Recent Enhancements to the Development of CFD-Based Aeroelastic Reduced-Order Models

Recent enhancements to the development of CFD-based unsteady aerodynamic and aeroelastic reduced-order models (ROMs) are presented. These enhancements include the simultaneous application of structural modes as CFD input, static aeroelastic analysis using a ROM, and matched-point solutions using a ROM. The simultaneous application of structural modes as CFD input enables the computation of the unsteady aerodynamic state-space matrices with a single CFD execution, independent of the number of structural modes. The responses obtained from a simultaneous excitation of the CFD-based unsteady aerodynamic system are processed using system identification techniques in order to generate an unsteady aerodynamic state-space ROM. Once the unsteady aerodynamic state-space ROM is generated, a method for computing the static aeroelastic response using this unsteady aerodynamic ROM and a state-space model of the structure, is presented. Finally, a method is presented that enables the computation of matchedpoint solutions using a single ROM that is applicable over a range of dynamic pressures and velocities for a given Mach number. These enhancements represent a significant advancement of unsteady aerodynamic and aeroelastic ROM technology