Anisotropic Flow in Event-by-Event Ideal Hydrodynamic Simulations of root s(NN) 200 GeV Au+Au Collisions

We simulate top-energy Au + Au collisions using ideal hydrodynamics in order to make the first comparison to the complete set of midrapidity flow measurements made by the PHENIX Collaboration. A simultaneous calculation of nu(2), nu(3), nu(4), and the first event-by-event calculation of quadrangular flow defined with respect to the nu(2) event plane (nu(4){Psi(2)}) gives good agreement with measured values, including the dependence on both transverse momentum and centrality. This provides confirmation that the collision system is indeed well described as a quark-gluon plasma with an extremely small viscosity and that correlations are dominantly generated from collective effects. In addition, we present a prediction for nu(5).

Charged and strange hadron elliptic flow in Cu plus Cu collisions at root s(NN)=62.4 and 200 GeV

We present the results of an elliptic flow, v(2), analysis of Cu + Cu collisions recorded with the solenoidal tracker detector (STAR) at the BNL Relativistic Heavy Ion Collider at root s(NN) = 62.4 and 200 GeV. Elliptic flow as a function of transverse momentum, v(2)(p(T)), is reported for different collision centralities for charged hadrons h(+/-) and strangeness-ontaining hadrons K-S(0), Lambda, Xi, and phi in the midrapidity region vertical bar eta vertical bar < 1.0. Significant reduction in systematic uncertainty of the measurement due to nonflow effects has been achieved by correlating particles at midrapidity, vertical bar eta vertical bar < 1.0, with those at forward rapidity, 2.5 < vertical bar eta vertical bar < 4.0. We also present azimuthal correlations in p + p collisions at root s = 200 GeV to help in estimating nonflow effects. To study the system-size dependence of elliptic flow, we present a detailed comparison with previously published results from Au + Au collisions at root s(NN) = 200 GeV. We observe that v(2)(p(T)) of strange hadrons has similar scaling properties as were first observed in Au + Au collisions, that is, (i) at low transverse momenta, p(T) < 2 GeV/c, v(2) scales with transverse kinetic energy, m(T) - m, and (ii) at intermediate p(T), 2 < p(T) < 4 GeV/c, it scales with the number of constituent quarks, n(q.) We have found that ideal hydrodynamic calculations fail to reproduce the centrality dependence of v(2)(p(T)) for K-S(0) and Lambda. Eccentricity scaled v(2) values, v(2)/epsilon, are larger in more central collisions, suggesting stronger collective flow develops in more central collisions. The comparison with Au + Au collisions, which go further in density, shows that v(2)/epsilon depends on the system size, that is, the number of participants N-part. This indicates that the ideal hydrodynamic limit is not reached in Cu + Cu collisions, presumably because the assumption of thermalization is not attained.

Detection of Groundwater flow zones in low productivity boreholes (Submitted)

Conventional methods of detecting groundwater flow zones in open boreholes installed in fractured bedrock aquifers rely on either contrasts in water quality or on significant rates of vertical flow. In low productivity boreholes these methods have proced of limited value. Tracer tests completed in six low productivity bedrock boreholes installed into fresh Precambrian crystalline bedrock revealed measurable dilution, despite very low tranmissivities. Fluorescent tracer profiles generated during discharge pumping permitted identification of the principle zones contributing to flow. High resolution acoustic televiewer strike and dip measurements of fractures corresponding to these zones revealed a consistent pattern of regional lineament trends and suggested a strongly anisotropic flow pattern in bedrock.

Fluctuating Initial Conditions and Anisotropic Flows

In this work we study the connection between anisotropic flows and lumpy initial conditions for Au+Au collisions at 200 GeV. We present comparisons between anisotropic flow coefficients and eccentricities up to sixth order, and between initial condition reference angles and azimuthal particle distribution angles. We also present a toy model to justify the lack of connection between flow coefficients and eccentricities for individual events.

Inclusive charged hadron elliptic flow in Au+Au collisions at root s(NN)=7.7-39 GeV

A systematic study is presented for centrality, transverse momentum (p(T)), and pseudorapidity (eta) dependence of the inclusive charged hadron elliptic flow (v(2)) at midrapidity (vertical bar eta vertical bar < 1.0) in Au + Au collisions at root s(NN) = 7.7, 11.5, 19.6, 27, and 39 GeV. The results obtained with different methods, including correlations with the event plane reconstructed in a region separated by a large pseudorapidity gap and four-particle cumulants (v(2){4}), are presented to investigate nonflow correlations and v(2) fluctuations. We observe that the difference between v(2){2} and v(2){4} is smaller at the lower collision energies. Values of v(2), scaled by the initial coordinate space eccentricity, v(2)/epsilon, as a function of p(T) are larger in more central collisions, suggesting stronger collective flow develops in more central collisions, similar to the results at higher collision energies. These results are compared to measurements at higher energies at the Relativistic Heavy Ion Collider (root s(NN) = 62.4 and 200 GeV) and at the Large Hadron Collider (Pb + Pb collisions at root s(NN) = 2.76 TeV). The v(2)(pT) values for fixed pT rise with increasing collision energy within the pT range studied (<2 GeV/c). A comparison to viscous hydrodynamic simulations is made to potentially help understand the energy dependence of v(2)(pT). We also compare the v(2) results to UrQMD and AMPT transport model calculations, and physics implications on the dominance of partonic versus hadronic phases in the system created at beam energy scan energies are discussed.

Directed Flow of Identified Particles in Au plus Au Collisions at root S-NN=200 GeV at RHIC

STAR's measurements of directed flow (v(1)) around midrapidity for pi(+/-), K-+/-, K-S(0), p, and (p) over bar in Au + Au collisions at root s(NN) = 200 GeV are presented. A negative v(1) (y) slope is observed for most of produced particles (pi(+/-), K-+/-, K-S(0), p, and (p) over bar). In 5%-30% central collisions, a sizable difference is present between the v(1)(y) slope of protons and antiprotons, with the former being consistent with zero within errors. The v(1) excitation function is presented. Comparisons to model calculations (RQMD, UrQMD, AMPT, QGSM with parton recombination, and a hydrodynamics model with a tilted source) are made. For those models which have calculations of v(1) for both pions and protons, none of them can describe v(1()y) forpions and protons simultaneously. The hydrodynamics model with a tilted source as currently implemented cannot explain the centrality dependence of the difference between the v(1)(y) slopes of protons and antiprotons.

Drifter observations in the summer time Bay of Biscay slope current

During the summer of 2012, 20 surface drifters drogued at 50 m depth were deployed on the continental slope to the north of the Bay of Biscay. Initially after release the drifters all crossed the slope, with 14 continuing equatorward, parallel to the slope following an absolute dynamic topography feature and 6 returning to the slope, in an eddy, visible in chlorophyll-a maps. Lagrangian statistics show an anisotropic flow field that becomes less tied to the absolute dynamic topography and increasingly dominated by diffusion and eddy processes. A weaker tie to the absolute dynamic topography allowed for total of 8 of the drifters crossed from the deep water onto the shelf, showing pathways for flow across the slope. A combination of drifter trajectories, absolute dynamic topography and chlorophyll-a concentration maps have been used to show that small anticyclonic eddies, tied to the complex slope topography provide a mechanism for on shelf transport. During the summer, the presence of these eddies can be seen in surface chlorophyll-a maps.

Drifter observations in the summer time Bay of Biscay slope current

During the summer of 2012, 20 surface drifters drogued at 50 m depth were deployed on the continental slope to the north of the Bay of Biscay. Initially after release the drifters all crossed the slope, with 14 continuing equatorward, parallel to the slope following an absolute dynamic topography feature and 6 returning to the slope, in an eddy, visible in chlorophyll-a maps. Lagrangian statistics show an anisotropic flow field that becomes less tied to the absolute dynamic topography and increasingly dominated by diffusion and eddy processes. A weaker tie to the absolute dynamic topography allowed for total of 8 of the drifters crossed from the deep water onto the shelf, showing pathways for flow across the slope. A combination of drifter trajectories, absolute dynamic topography and chlorophyll-a concentration maps have been used to show that small anticyclonic eddies, tied to the complex slope topography provide a mechanism for on shelf transport. During the summer, the presence of these eddies can be seen in surface chlorophyll-a maps.

Study of bulk matter properties through strange quark hadrons with the STAR experiment

Relativistic heavy ion collisions are the ideal experimental tool to explore the QCD phase diagram. Several results show that a very hot medium with a high energy density and partonic degrees of freedom is formed in these collisions, creating a new state of matter. Measurements of strange hadrons can bring important information about the bulk properties of such matter. The elliptic flow of strange hadrons such as phi, K(S)(0), Lambda, Xi and Omega shows that collectivity is developed at partonic level and at intermediate p(T) the quark coalescence is the dominant mechanism of hadronization. The nuclear modification factor is an another indicator of the presence of a very dense medium. The comparison between measurements of Au+Au and d+Au collisions, where only cold nuclear matter effects are expected, can shed more light on the bulk properties. In these proceedings, recent results from the STAR experiment on bulk matter properties are presented.

Mapping the hydrodynamic response to the initial geometry in heavy-ion collisions

We investigate how the initial geometry of a heavy-ion collision is transformed into final flow observables by solving event-by-event ideal hydrodynamics with realistic fluctuating initial conditions. We study quantitatively to what extent anisotropic flow (nu(n)) is determined by the initial eccentricity epsilon(n) for a set of realistic simulations, and we discuss which definition of epsilon(n) gives the best estimator of nu(n). We find that the common practice of using an r(2) weight in the definition of epsilon(n) in general results in a poorer predictor of nu(n) than when using r(n) weight, for n > 2. We similarly study the importance of additional properties of the initial state. For example, we show that in order to correctly predict nu(4) and nu(5) for noncentral collisions, one must take into account nonlinear terms proportional to epsilon(2)(2) and epsilon(2)epsilon(3), respectively. We find that it makes no difference whether one calculates the eccentricities over a range of rapidity or in a single slice at z = 0, nor is it important whether one uses an energy or entropy density weight. This knowledge will be important for making a more direct link between experimental observables and hydrodynamic initial conditions, the latter being poorly constrained at present.

Measurement with the ATLAS detector of multi-particle azimuthal correlations in p+Pb collisions at sqrtsNN=5.02 TeV

In order to study further the long-range correlations ("ridge") observed recently in p+Pb collisions at sqrt(s_NN) =5.02 TeV, the second-order azimuthal anisotropy parameter of charged particles, v_2, has been measured with the cumulant method using the ATLAS detector at the LHC. In a data sample corresponding to an integrated luminosity of approximately 1 microb^(-1), the parameter v_2 has been obtained using two- and four-particle cumulants over the pseudorapidity range |eta|<2.5. The results are presented as a function of transverse momentum and the event activity, defined in terms of the transverse energy summed over 3.1anisotropic flow, similar to that observed in Pb+Pb collisions. A comparison is made to results obtained using two-particle correlation methods, and to predictions from hydrodynamic models of p+Pb collisions. Despite the small transverse spatial extent of the p+Pb collision system, the large magnitude of v_2 and its similarity to hydrodynamic predictions provide additional evidence for the importance of final-state effects in p+Pb reactions.

Stability of Poiseuille flow in a fluid overlying an anisotropic and inhomogeneous porous layer

We present the linear stability analysis of horizontal Poiseuille flow in a fluid overlying a porous medium with anisotropic and inhomogeneous permeability. The generalized Darcy model is used to describe the flow in the porous medium with the Beavers-Joseph condition at the interface of the two layers and the eigenvalue problem is solved numerically. The effect of major system parameters on the stability characteristics is addressed in detail. It is shown that the anisotropic and inhomogeneous modulation of the permeability of the underlying porous layer provides an effective means for passive control of the flow stability.

Anisotropic thermopiezoelectric solids with an elliptic inclusion or a hole under uniform heat flow

The two-dimensional problem of a thermopiezoelectric material containing an elliptic inclusion or a hole subjected to a remote uniform heat flow is studied. Based on the extended Lekhnitskii formulation for thermopiezoelectricity, conformal mapping and Laurent series expansion, the explicit and closed-form solutions are obtained both inside and outside the inclusion (or hole). For a hole problem, the exact electric boundary conditions on the hole surface are used. The results show that the electroelastic fields inside the inclusion or the electric field inside the hole are linear functions of the coordinates. When the elliptic hole degenerates into a slit crack, the electroelastic fields and the intensity factors are obtained. The effect of the heat how direction and the dielectric constant of air inside the crack on the thermal electroelastic fields are discussed. Comparison is made with two special cases of which the closed solutions exist and it is shown that our results are valid.