An Eulerian Immersed Boundary Method for flow simulations over stationary and moving rigid bodies

The fluid flow over bodies with complex geometry has been the subject of research of many scientists and widely explored experimentally and numerically. The present study proposes an Eulerian Immersed Boundary Method for flows simulations over stationary or moving rigid bodies. The proposed method allows the use of Cartesians Meshes. Here, two-dimensional simulations of fluid flow over stationary and oscillating circular cylinders were used for verification and validation. Four different cases were explored: the flow over a stationary cylinder, the flow over a cylinder oscillating in the flow direction, the flow over a cylinder oscillating in the normal flow direction, and a cylinder with angular oscillation. The time integration was carried out by a classical 4th order Runge-Kutta scheme, with a time step of the same order of distance between two consecutive points in x direction. High-order compact finite difference schemes were used to calculate spatial derivatives. The drag and lift coefficients, the lock-in phenomenon and vorticity contour plots were used for the verification and validation of the proposed method. The extension of the current method allowing the study of a body with different geometry and three-dimensional simulations is straightforward. The results obtained show a good agreement with both numerical and experimental results, encouraging the use of the proposed method.

Numerical investigation of the three-dimensional secondary instabilities in the time-developing compressible mixing layer

Mixing layers are present in very different types of physical situations such as atmospheric flows, aerodynamics and combustion. It is, therefore, a well researched subject, but there are aspects that require further studies. Here the instability of two-and three-dimensional perturbations in the compressible mixing layer was investigated by numerical simulations. In the numerical code, the derivatives were discretized using high-order compact finite-difference schemes. A stretching in the normal direction was implemented with both the objective of reducing the sound waves generated by the shear region and improving the resolution near the center. The compact schemes were modified to work with non-uniform grids. Numerical tests started with an analysis of the growth rate in the linear regime to verify the code implementation. Tests were also performed in the non-linear regime and it was possible to reproduce the vortex roll-up and pairing, both in two-and three-dimensional situations. Amplification rate analysis was also performed for the secondary instability of this flow. It was found that, for essentially incompressible flow, maximum growth rates occurred for a spanwise wavelength of approximately 2/3 of the streamwise spacing of the vortices. The result demonstrated the applicability of the theory developed by Pierrehumbet and Widnall. Compressibility effects were then considered and the maximum growth rates obtained for relatively high Mach numbers (typically under 0.8) were also presented.

Effects of Variations in Nonlinear Damping Coefficients on the Parametric Vibration of a Cantilever Beam with a Lumped Mass

Uncertainties in damping estimates can significantly affect the dynamic response of a given flexible structure. A common practice in linear structural dynamics is to consider a linear viscous damping model as the major energy dissipation mechanism. However, it is well known that different forms of energy dissipation can affect the structure's dynamic response. The major goal of this paper is to address the effects of the turbulent frictional damping force, also known as drag force on the dynamic behavior of a typical flexible structure composed of a slender cantilever beam carrying a lumped-mass on the tip. First, the system's analytical equation is obtained and solved by employing a perturbation technique. The solution process considers variations of the drag force coefficient and its effects on the system's response. Then, experimental results are presented to demonstrate the effects of the nonlinear quadratic damping due to the turbulent frictional force on the system's dynamic response. In particular, the effects of the quadratic damping on the frequency-response and amplitude-response curves are investigated. Numerically simulated as well as experimental results indicate that variations on the drag force coefficient significantly alter the dynamics of the structure under investigation. Copyright (c) 2008 D. G. Silva and P. S. Varoto.

TURBULENCE AND THE FORMATION OF FILAMENTS, LOOPS, AND SHOCK FRONTS IN NGC 1275

NGC 1275, the central galaxy in the Perseus cluster, is the host of gigantic hot bipolar bubbles inflated by active galactic nucleus (AGN) jets observed in the radio as Perseus A. It presents a spectacular H alpha-emitting nebulosity surrounding NGC 1275, with loops and filaments of gas extending to over 50 kpc. The origin of the filaments is still unknown, but probably correlates with the mechanism responsible for the giant buoyant bubbles. We present 2.5 and three-dimensional magnetohydrodynamical (MHD) simulations of the central region of the cluster in which turbulent energy, possibly triggered by star formation and supernovae (SNe) explosions, is introduced. The simulations reveal that the turbulence injected by massive stars could be responsible for the nearly isotropic distribution of filaments and loops that drag magnetic fields upward as indicated by recent observations. Weak shell-like shock fronts propagating into the intracluster medium (ICM) with velocities of 100-500 km s(-1) are found, also resembling the observations. The isotropic outflow momentum of the turbulence slows the infall of the ICM, thus limiting further starburst activity in NGC 1275. As the turbulence is subsonic over most of the simulated volume, the turbulent kinetic energy is not efficiently converted into heat and additional heating is required to suppress the cooling flow at the core of the cluster. Simulations combining the MHD turbulence with the AGN outflow can reproduce the temperature radial profile observed around NGC 1275. While the AGN mechanism is the main heating source, the SNe are crucial to isotropize the energy distribution.

Spin-charge coupling in quantum wires at zero magnetic field

We discuss an approximation for the dynamic charge response of nonlinear spin-1/2 Luttinger liquids in the limit of small momentum. Besides accounting for the broadening of the charge peak due to two-holon excitations, the nonlinearity of the dispersion gives rise to a two-spinon peak, which at zero temperature has an asymmetric line shape. At finite temperature the spin peak is broadened by diffusion. As an application, we discuss the density and temperature dependence of the Coulomb drag resistivity due to long-wavelength scattering between quantum wires.

Joint forces and torques when walking in shallow water

This study reports for the first time an estimation of the internal net joint forces and torques on adults` lower limbs and pelvis when walking in shallow water, taking into account the drag forces generated by the movement of their bodies in the water and the equivalent data when they walk on land. A force plate and a video camera were used to perform a two-dimensional gait analysis at the sagittal plane of 10 healthy young adults walking at comfortable speeds on land and in water at a chest-high level. We estimated the drag force on each body segment and the joint forces and torques at the ankle, knee, and hip of the right side of their bodies using inverse dynamics. The observed subjects` apparent weight in water was about 35% of their weight on land and they were about 2.7 times slower when walking in water. When the subjects walked in water compared with walking on land, there were no differences in the angular displacements but there was a significant reduction in the joint torques which was related to the water`s depth. The greatest reduction was observed for the ankle and then the knee and no reduction was observed for the hip. All joint powers were significantly reduced in water. The compressive and shear joint forces were on average about three times lower during walking in water than on land. These quantitative results substantiate the use of water as a safe environment for practicing low-impact exercises, particularly walking. (C) 2011 Elsevier Ltd. All rights reserved.

Evaluation of hydrodynamic parameters of a fluidized-bed reactor with immobilized yeast

The fluidized bed reactor has successfully been used to perform biotechnological processes addressed to the production of high added value. The present work evaluates hydrodynamic parameters of a bench-scale fluidized bed reactor with cells of the yeast Candida guilliermondii immobilized either in calcium alginate beads or in polyvinyl alcohol (PVA). The effects of the following variables on cell immobilization were evaluated at 30 degrees C and feeding a synthetic medium containing 50 g L-1 xylose: total particle density (cells plus support), terminal velocity, particle drag force, minimum fluidization velocity and bed porosity. According to the results obtained, the reactor was shown to operate like a fixed-bed bioreactor at xi < 0.5 and a fluidized bed bioreactor at xi > 0.5. The maximum flow rate needed to obtain maximum bed fluidization in the reactor was equal to the terminal velocity of the immobilized cell particles. Particles of cells immobilized within these supports showed values of drag coefficient lower than those reported for other high-density supports. The evaluation of these hydrodynamic characteristics lead to an adequate bed fluidization inside the reactor, thus improving oxygen transference and availability in the fermentation medium, making the process more viable for future scale-up. (c) 2008 Society of Chemical Industry.

Enhanced aeroelastic energy harvesting by exploiting combined nonlinearities: theory and experiment

Converting aeroelastic vibrations into electricity for low power generation has received growing attention over the past few years. In addition to potential applications for aerospace structures, the goal is to develop alternative and scalable configurations for wind energy harvesting to use in wireless electronic systems. This paper presents modeling and experiments of aeroelastic energy harvesting using piezoelectric transduction with a focus on exploiting combined nonlinearities. An airfoil with plunge and pitch degrees of freedom (DOF) is investigated. Piezoelectric coupling is introduced to the plunge DOF while nonlinearities are introduced through the pitch DOF. A state-space model is presented and employed for the simulations of the piezoaeroelastic generator. A two-state approximation to Theodorsen aerodynamics is used in order to determine the unsteady aerodynamic loads. Three case studies are presented. First the interaction between piezoelectric power generation and linear aeroelastic behavior of a typical section is investigated for a set of resistive loads. Model predictions are compared to experimental data obtained from the wind tunnel tests at the flutter boundary. In the second case study, free play nonlinearity is added to the pitch DOF and it is shown that nonlinear limit-cycle oscillations can be obtained not only above but also below the linear flutter speed. The experimental results are successfully predicted by the model simulations. Finally, the combination of cubic hardening stiffness and free play nonlinearities is considered in the pitch DOF. The nonlinear piezoaeroelastic response is investigated for different values of the nonlinear-to-linear stiffness ratio. The free play nonlinearity reduces the cut-in speed while the hardening stiffness helps in obtaining persistent oscillations of acceptable amplitude over a wider range of airflow speeds. Such nonlinearities can be introduced to aeroelastic energy harvesters (exploiting piezoelectric or other transduction mechanisms) for performance enhancement.

Secondary instabilities in the flow around two circular cylinders in tandem

Direct stability analysis and numerical simulations have been employed to identify and characterize secondary instabilities in the wake of the flow around two identical circular cylinders in tandem arrangements. The centre-to-centre separation was varied from 1.2 to 10 cylinder diameters. Four distinct regimes were identified and salient cases chosen to represent the different scenarios observed, and for each configuration detailed results are presented and compared to those obtained for a flow around an isolated cylinder. It was observed that the early stages of the wake transition changes significantly if the separation is smaller than the drag inversion spacing. The onset of the three-dimensional instabilities were calculated and the unstable modes are fully described. In addition, we assessed the nonlinear character of the bifurcations and physical mechanisms are proposed to explain the instabilities. The dependence of the critical Reynolds number on the centre-to-centre separation is also discussed.

Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

The Equilibrium Flux Method for the calculation of flows with non-equilbrium chemical reactions

The Equilibrium Flux Method [1] is a kinetic theory based finite volume method for calculating the flow of a compressible ideal gas. It is shown here that, in effect, the method solves the Euler equations with added pseudo-dissipative terms and that it is a natural upwinding scheme. The method can be easily modified so that the flow of a chemically reacting gas mixture can be calculated. Results from the method for a one-dimensional non-equilibrium reacting flow are shown to agree well with a conventional continuum solution. Results are also presented for the calculation of a plane two-dimensional flow, at hypersonic speed, of a dissociating gas around a blunt-nosed body.

Comparison of experimental and numerical studies of ionizing flow over a cylinder

Comparisons are made between experimental measurements and numerical simulations of ionizing flows generated in a superorbital facility. Nitrogen, with a freestream velocity of around 10 km/s, was passed over a cylindrical model, and images were recorded using two-wavelength holographic interferometry. The resulting density, electron concentration, and temperature maps were compared with numerical simulations from the Langley Research Center aerothermodynamic upwind relaxation algorithm. The results showed generally good agreement in shock location and density distributions. Some discrepancies were observed for the electron concentration, possibly, because simulations were of a two-dimensional flow, whereas the experiments were likely to have small three-dimensional effects.

New calibration technique for multiple-component stress wave force balances

Force measurement in hypervelocity expansion tubes is not possible using conventional techniques. The stress wave force balance technique can be applied in expansion tubes to measure forces despite the short test times involved. This paper presents a new calibration technique for multiple-component stress wave force balances where an impulse response created using a load distribution is required and no orthogonal surfaces on the model exist.. This new technique relies on the tensorial superposition of single-component impulse responses analogous to the vectorial superposition of the calibration loads. The example presented here is that of a scale model of the Mars Pathfinder, but the technique is applicable to any geometry and may be useful for cases where orthogonal loads cannot be applied.

Neighborhood and Efficiency in Manufacturing in Brazilian Regions: A Spatial Markov Chain Analysis

This paper analyzes the geography of regional competitiveness in manufacturing in Brazil. The authors estimate stochastic frontiers to calculate regional efficiency of representative firms in 137 regions in the period 2000-2006, in four sectors defined by technological intensity. The efficiency results are analyzed using Markov Spatial Transition Matrices to provide insights into the transition of regions between efficiency levels, considering their local spatial context. The results indicate that geography plays an important role in manufacturing competitiveness. In particular, regions with more competitive neighbors are more likely to improve their relative efficiency (pull effect) over time, and regions with less competitive neighbors are more likely to lose relative efficiency (drag effect). The authors find that the pull effect is stronger than the drag effect.

The nymph of Tortopus harrisi Traver (Ephemeroptera: Polymitarcyidae)

Polymitarcyidae is a family of burrowing mayflies (Ephemeroptera: Ephemeroidea) distributed throughout the world but with highest diversity in the Neotropics. Tortopus Needham & Murphy, with a Panamerican distribution, is known from twelve species described in the adult stage. Nymphs are only known for three species: T. puella (Pictet), T. obscuripennis Dominguez and T. sarae Dominguez, and present a rather homogeneous morphology (Molineri 2008). They were firstly described for T. puella by Scott et al. (1959) and later Molineri (2008) described the other two. Both studies reported that these species burrow U-shaped tunnels in clay banks of rivers and streams, thus preventing them from being sampled in most limnological studies (that use surbers, drags, or drift nets). The aim of the present contribution is to describe and illustrate the previously unknown nymph of Tortopus harrisi Traver that shows important anatomical differences with the other nymphs known in the genus. This morphological differentiation suggests a different habitat use by these nymphs, sampled with drag and surber samplers in sandy substrate. New locality records are given for T. harrisi in Brazil. The nymphs are preserved in alcohol, mouthparts, legs and genital rudiments were mounted in microscope slides with Canada Balsam. Drawings were made with a camera lucida attached to a stereo microscope. The material is deposited in CUIC (Cornell University Insect Collection, Ithaca, NY), IML (Instituto Miguel Lillo, Tucuman) and in MZSP (Museu de Zoologia da Universidade de Sao Paulo, Sao Paulo). Catalogs and bibliography were consulted at Ephemeroptera Galactica (Hubbard 2009).