Numerical study of Mach number and thermal effects on sound radiation by a mixing layer

Mach number and thermal effects on the mechanisms of sound generation and propagation are investigated in spatially evolving two-dimensional isothermal and non-isothermal mixing layers at Mach number ranging from 0.2 to 0.4 and Reynolds number of 400. A characteristic-based formulation is used to solve by direct numerical simulation the compressible Navier-Stokes equations using high-order schemes. The radiated sound is directly computed in a domain that includes both the near-field aerodynamic source region and the far-field sound propagation. In the isothermal mixing layer, Mach number effects may be identified in the acoustic field through an increase of the directivity associated with the non-compactness of the acoustic sources. Baroclinic instability effects may be recognized in the non-isothermal mixing layer, as the presence of counter-rotating vorticity layers, the resulting acoustic sources being found less efficient. An analysis based on the acoustic analogy shows that the directivity increase with the Mach number can be associated with the emergence of density fluctuations of weak amplitude but very efficient in terms of noise generation at shallow angle. This influence, combined with convection and refraction effects, is found to shape the acoustic wavefront pattern depending on the Mach number.

Simulation results and applications of an advection bounded scheme to practical flows

This paper reports experiments on the use of a recently introduced advection bounded upwinding scheme, namely TOPUS (Computers & Fluids 57 (2012) 208-224), for flows of practical interest. The numerical results are compared against analytical, numerical and experimental data and show good agreement with them. It is concluded that the TOPUS scheme is a competent, powerful and generic scheme for complex flow phenomena.

On the radial expansion of tubular structures in a quark-gluon plasma

We study the radial expansion of cylindrical tubes in a hot QGP. These tubes are treated as perturbations in the energy density of the system which is formed in heavy ion collisions at RHIC and LHC. We start from the equations of relativistic hydrodynamics in two spatial dimensions and cylindrical symmetry and perform an expansion of these equations in a small parameter, conserving the nonlinearity of the hydrodynamical formalism. We consider both ideal and viscous fluids and the latter are studied with a relativistic Navier-Stokes equation. We use the equation of state of the MIT bag model. In the case of ideal fluids we obtain a breaking wave equation for the energy density fluctuation, which is then solved numerically. We also show that, under certain assumptions, perturbations in a relativistic viscous fluid are governed by the Burgers equation. We estimate the typical expansion time of the tubes. (C) 2012 Elsevier B.V. All rights reserved.

A new enrichment space for the treatment of discontinuous pressures in multi-fluid flows

In this work, a new enrichment space to accommodate jumps in the pressure field at immersed interfaces in finite element formulations, is proposed. The new enrichment adds two degrees of freedom per element that can be eliminated by means of static condensation. The new space is tested and compared with the classical P1 space and to the space proposed by Ausas et al (Comp. Meth. Appl. Mech. Eng., Vol. 199, 10191031, 2010) in several problems involving jumps in the viscosity and/or the presence of singular forces at interfaces not conforming with the element edges. The combination of this enrichment space with another enrichment that accommodates discontinuities in the pressure gradient has also been explored, exhibiting excellent results in problems involving jumps in the density or the volume forces. Copyright (c) 2011 John Wiley & Sons, Ltd.

URANS calculations for smooth circular cylinder flow in a wide range of Reynolds Numbers: solution verification and validation

The flow around circular smooth fixed cylinder in a large range of Reynolds numbers is considered in this paper. In order to investigate this canonical case, we perform CFD calculations and apply verification & validation (V&V) procedures to draw conclusions regarding numerical error and, afterwards, assess the modeling errors and capabilities of this (U)RANS method to solve the problem. Eight Reynolds numbers between Re = 10 and Re 5 x 10(5) will be presented with, at least, four geometrically similar grids and five discretization in time for each case (when unsteady), together with strict control of iterative and round-off errors, allowing a consistent verification analysis with uncertainty estimation. Two-dimensional RANS, steady or unsteady, laminar or turbulent calculations are performed. The original 1994 k - omega SST turbulence model by Menter is used to model turbulence. The validation procedure is performed by comparing the numerical results with an extensive set of experimental results compiled from the literature. [DOI: 10.1115/1.4007571]

Average control of Markov decision processes with Feller transition probabilities and general action spaces

This paper studies the average control problem of discrete-time Markov Decision Processes (MDPs for short) with general state space, Feller transition probabilities, and possibly non-compact control constraint sets A(x). Two hypotheses are considered: either the cost function c is strictly unbounded or the multifunctions A(r)(x) = {a is an element of A(x) : c(x, a) <= r} are upper-semicontinuous and compact-valued for each real r. For these two cases we provide new results for the existence of a solution to the average-cost optimality equality and inequality using the vanishing discount approach. We also study the convergence of the policy iteration approach under these conditions. It should be pointed out that we do not make any assumptions regarding the convergence and the continuity of the limit function generated by the sequence of relative difference of the alpha-discounted value functions and the Poisson equations as often encountered in the literature. (C) 2012 Elsevier Inc. All rights reserved.

Dynamics of phytoplankton associations in three reservoirs in northeastern Brazil assessed using Reynolds' theory

The aim of the present study was to evaluate the influence of seasonality on the behavior of phytoplankton associations in eutrophic reservoirs with different depths in northeastern Brazil. Five collections were carried out at each of the reservoirs at two depths (0.1 m and near the sediment) at three-month intervals in each season (dry and rainy). The phytoplankton samples were preserved in Lugol's solution and quantified under an inverted microscope for the determination of density values, which were subsequently converted to biomass values based on cellular biovolume and classified in phytoplankton associations. The following abiotic variables were analyzed: water temperature, dissolved oxygen, pH, turbidity, water transparency, total phosphorus, total dissolved phosphorus, orthophosphate and total nitrogen. The data were investigated using canonical correspondence analysis. The influence of seasonality on the dynamics of the phytoplankton community was lesser in the deeper reservoirs. Depth affected the behavior of the algal associations. Variation in light availability was a determinant of changes in the phytoplankton structure. Urosolenia and Anabaena associations were more abundant in shallow ecosystems with a larger eutrophic zone, whereas the Microcystis association was more related to deep ecosystems with adequate availability of nutrients. The distribution of Cyclotella, Geitlerinema, Planktothrix, Pseudanabaena and Cylindrospermopsis associations was different from that seen in subtropical regions and the substitution of these associations was related to a reduction in the eutrophic zone rather than the mixture zone. Published by Elsevier GmbH.

Mapping the Energy Distribution of SERRS Hot Spots from Anti-Stokes to Stokes Intensity Ratios

The anomalies in the anti-Stokes to Stokes intensity ratios in single-molecule surface-enhanced resonance Raman scattering were investigated. Brilliant green and crystal violet dyes were the molecular probes, and the experiments were carried out on an electrochemically activated Ag surface. The results allowed new insights into the origin of these anomalies and led to a new method to confirm the single-molecule regime in surface-enhanced Raman scattering. Moreover, a methodology to estimate the distribution of resonance energies that contributed to the imbalance in the anti-Stokes to Stokes intensity ratios at the electromagnetic hot spots was proposed. This method allowed the local plasmonic resonance energies on the metallic surface to be spatially mapped.

Numerical simulation of an ethanol turbulent spray flame with RANS and diffusion combustion model

A detailed numerical simulation of ethanol turbulent spray combustion on a rounded jet flame is pre- sented in this article. The focus is to propose a robust mathematical model with relatively low complexity sub- models to reproduce the main characteristics of the cou- pling between both phases, such as the turbulence modulation, turbulent droplets dissipation, and evaporative cooling effect. A RANS turbulent model is implemented. Special features of the model include an Eulerian– Lagrangian procedure under a fully two-way coupling and a modified flame sheet model with a joint mixture fraction– enthalpy b -PDF. Reasonable agreement between measured and computed mean profiles of temperature of the gas phase and droplet size distributions is achieved. Deviations found between measured and predicted mean velocity profiles are attributed to the turbulent combustion modeling adopted

Finite element methods for the Stokes system based on a Zienkiewicz type N-simplex

Hermite interpolation is increasingly showing to be a powerful numerical solution tool, as applied to different kinds of second order boundary value problems. In this work we present two Hermite finite element methods to solve viscous incompressible flows problems, in both two- and three-dimension space. In the two-dimensional case we use the Zienkiewicz triangle to represent the velocity field, and in the three-dimensional case an extension of this element to tetrahedra, still called a Zienkiewicz element. Taking as a model the Stokes system, the pressure is approximated with continuous functions, either piecewise linear or piecewise quadratic, according to the version of the Zienkiewicz element in use, that is, with either incomplete or complete cubics. The methods employ both the standard Galerkin or the Petrov–Galerkin formulation first proposed in Hughes et al. (1986) [18], based on the addition of a balance of force term. A priori error analyses point to optimal convergence rates for the PG approach, and for the Galerkin formulation too, at least in some particular cases. From the point of view of both accuracy and the global number of degrees of freedom, the new methods are shown to have a favorable cost-benefit ratio, as compared to velocity Lagrange finite elements of the same order, especially if the Galerkin approach is employed.