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.

Wake instability issues: From circular cylinders to stalled airfoils

Some recent results regarding the global dynamical behaviour of the wake of circular cylinders and airfoils with massive separation are reviewed in this paper. In order to investigate the effect of interference, the three-dimensional instability modes are analysed for the flow around two circular cylinders in tandem. In the same way, the flow around a stalled airfoil is investigated in order to provide a better understanding of the three-dimensional characteristics of wakes forming downstream of a lifting body with massive separation. These results are compared with those found for an isolated cylinder. Some fundamental differences among these flows are discussed. (C) 2011 Elsevier Ltd. All rights reserved.

Floquet stability analysis of the flow around an oscillating cylinder

This investigative work is concerned with the flow around a circular cylinder submitted to forced transverse oscillations. The goal is to investigate how the transition to turbulence is initiated in the wake for cases with different Reynolds numbers (Re) and displacement amplitudes (A). For each Re the motion frequency is kept constant, close to the Strouhal number of the flow around a fixed cylinder at the same Re. Stability analysis of two-dimensional periodic flows around a forced-oscillating cylinder is carried out with respect to three-dimensional infinitesimal perturbations. The procedure consists of performing a Floquet type analysis of time-periodic base flows, computed using the spectral/hp element method. With the results of the Floquet calculations, considerations regarding the stability of the system are drawn, and the form of the instability at its onset is obtained. The critical Reynolds number is observed to change with the amplitude of oscillation. With respect to instabilities, unstable modes with the same symmetry as mode A of a fixed cylinder are observed; however, they present different wavelengths. Also, the instabilities observed for the oscillating cylinder are distinctively stronger in the braid shear layers. Other unstable modes similar to mode B are found. Quasi-periodic modes are observed in the 2S wake, and subharmonic mode occurrences are reported in P + S wakes. (C) 2009 Elsevier Ltd. All rights reserved.

Safety of a railway wheelset - derailment simulation with increasing lateral force

The motivation for this research is to make a comparison between dynamic results of a free railway wheelset derailment and safety limits. For this purpose, a numerical simulation of a wheelset derailment submitted to increasing lateral force is used to compare with the safety limit, using different criteria. A simplified wheelset model is used to simulate derailments with different adhesion conditions. The contact force components, including the longitudinal and spin effects, are identified in a steady-state condition on the verge of a derailment. The contact force ratios are used in a three-dimensional (3D) analytical formula to calculate the safety limits. Simulation results obtained with two contact methods were compared with the published results and the safety limit was identified with the two criteria. Results confirm Nadal`s conservative aspect and show that safety 3D analytical formula presents slightly higher safety limits for lower friction coefficients and smaller limits for high friction, in comparison with the simulation results with Fastsim.

Effect of the domain spanwise periodic length on the flow around a circular cylinder

We assess the effect of the choice of spanwise periodic length on simulations of the flow around a fixed circular cylinder. The Reynolds number is set to 400 because, at this value, both lift coefficient and shedding frequency show significant drop due to three-dimensional flow structures. From the analysis of the three-dimensionalities of the wake and of the integral quantities such as Strouhal number, RMS of lift coefficient and energy contained in the three-dimensional portion of the flow we obtain an estimate of the minimum spanwise length to satisfactorily represent the flow. Furthermore, we observe a distinct wake behavior when the spanwise length is approximately the mode B instability wavelength. (C) 2011 Elsevier Ltd. All rights reserved.

Sphericity of apatite particles determined by gas permeability through packed beds

Because shape is an assessment of the three-dimensional form of a particle, it may be described in terms of sphericity (Psi), which is a measure of how closely a particle approaches a spherical configuration. In this study, Darcy`s law and the Kozeny-Carman model for fluid flow through porous media were applied to packed beds to determine the sphericity (Psi) of apatite particles. The beds were composed of glass spheres or particles of apatite (igneous from Brazil and sedimentary from the United States) of three classes of size (Class 1: -297 +210 mu m; Class 2: -210 +149 mu m; Class 3: -149 +105 mu m). Glass spheres were used to validate the model because of its known sphericity (Psi = 1.00). Apatite particles, either igneous or sedimentary, showed very close values for particle sphericity (Psi approximate to 0.6). Observations on particle images conducted by scanning electron microscopy illustrated that igneous (Psi = 0.623) and sedimentary (Psi = 0.644) particles of apatite of Class 2 predominantly exhibit elongated shape. The close value of particle sphericity (Psi approximate to 0.6) showed by either igneous or sedimentary apatite may be justified by the similarity in particle shape.

Identification of Elastic, Dielectric, and Piezoelectric Constants in Piezoceramic Disks

Three-dimensional modeling of piezoelectric devices requires a precise knowledge of piezoelectric material parameters. The commonly used piezoelectric materials belong to the 6mm symmetry class, which have ten independent constants. In this work, a methodology to obtain precise material constants over a wide frequency band through finite element analysis of a piezoceramic disk is presented. Given an experimental electrical impedance curve and a first estimate for the piezoelectric material properties, the objective is to find the material properties that minimize the difference between the electrical impedance calculated by the finite element method and that obtained experimentally by an electrical impedance analyzer. The methodology consists of four basic steps: experimental measurement, identification of vibration modes and their sensitivity to material constants, a preliminary identification algorithm, and final refinement of the material constants using an optimization algorithm. The application of the methodology is exemplified using a hard lead zirconate titanate piezoceramic. The same methodology is applied to a soft piezoceramic. The errors in the identification of each parameter are statistically estimated in both cases, and are less than 0.6% for elastic constants, and less than 6.3% for dielectric and piezoelectric constants.

Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

Registration of temporal sequences of coronal and sagittal MR images through respiratory patterns

This work discusses the determination of the breathing patterns in time sequence of images obtained from magnetic resonance (MR) and their use in the temporal registration of coronal and sagittal images. The registration is made without the use of any triggering information and any special gas to enhance the contrast. The temporal sequences of images are acquired in free breathing. The real movement of the lung has never been seen directly, as it is totally dependent on its surrounding muscles and collapses without them. The visualization of the lung in motion is an actual topic of research in medicine. The lung movement is not periodic and it is susceptible to variations in the degree of respiration. Compared to computerized tomography (CT), MR imaging involves longer acquisition times and it is preferable because it does not involve radiation. As coronal and sagittal sequences of images are orthogonal to each other, their intersection corresponds to a segment in the three-dimensional space. The registration is based on the analysis of this intersection segment. A time sequence of this intersection segment can be stacked, defining a two-dimension spatio-temporal (2DST) image. The algorithm proposed in this work can detect asynchronous movements of the internal lung structures and lung surrounding organs. It is assumed that the diaphragmatic movement is the principal movement and all the lung structures move almost synchronously. The synchronization is performed through a pattern named respiratory function. This pattern is obtained by processing a 2DST image. An interval Hough transform algorithm searches for synchronized movements with the respiratory function. A greedy active contour algorithm adjusts small discrepancies originated by asynchronous movements in the respiratory patterns. The output is a set of respiratory patterns. Finally, the composition of coronal and sagittal image pairs that are in the same breathing phase is realized by comparing of respiratory patterns originated from diaphragmatic and upper boundary surfaces. When available, the respiratory patterns associated to lung internal structures are also used. The results of the proposed method are compared with the pixel-by-pixel comparison method. The proposed method increases the number of registered pairs representing composed images and allows an easy check of the breathing phase. (C) 2010 Elsevier Ltd. All rights reserved.

Threshold voltages of SOI MuGFETs

The multiple-gate field-effect transistor (MuGFET) is a device with a gate folded on different sides of the channel region. They are one of the most promising technological solutions to create high-performance ultra-scaled SOI CMOS. In this work, the behavior of the threshold voltage in double-gate, triple-gate and quadruple-gate SOI transistors with different channel doping concentrations is studied through three-dimensional numerical simulation. The results indicated that for double-gate transistors, one or two threshold voltages can be observed, depending on the channel doping concentration. However, in triple-gate and quadruple-gate it is possible to observe up to four threshold voltages due to the corner effect and the different doping concentration between the top and bottom of the Fin. (C) 2008 Elsevier Ltd. All rights reserved.

Prospective ultramorphological characterization of human hair by optical coherence tomography

Background/purpose The continuous advancement in cosmetic science has led to an increasing demand for the development of non-invasive, reliable scientific techniques directed toward claim substantiation, which is of utmost relevance, to obtain data regarding the efficacy and safety of cosmetic products. Methods In this work, we used the optical coherence tomography (OCT) technique to produce in vitro transversal section-images of human hair. We also compared the OCT signal before and after chemical treatment with an 18% w/w ammonium thioglycolate solution. Results The mean diameter of the medulla was 29 +/- 7 mu m and the hair diameter was 122 +/- 16 mu m in our samples of standard Afro-ethnic hair. A three-dimensional (3D) image was constructed starting from 601 cross-sectional images (slices). Each slice was taken in steps of 6.0 mu m at eight frames per second, and the entire 3D image was constructed in 60 s. Conclusion It was possible to identify, using the A-scan protocol, the principal structures: the cuticle, cortex and medulla. After chemical treatment, it was not possible to identify the main structures of hair fiber due to index matching promoted by deleterious action of the chemical agent.

Generation and Analysis of Interaction Energy Maps of p-Substituted Benzoic Acid [(5-Nitro-thiophen-2-yl)-methylene]-Hydrazides Active against MRSA

In this preliminary study eighteen p-substituted benzoic acid [(5-nitro-thiophen-2-yl)-methylene]-hydrazides with antimicrobial activity were evaluated against multidrug-resistant Staphylococcus aureus, correlating the three-dimensional characteristics of the ligands with their respective bioactivities. The computer programs Sybyl and CORINA were used, respectively, for the design and three-dimensional conversion of the ligands. Molecular interaction fields were calculated using GRID program. Calculations using Volsurf resulted in a statistically consistent model with 48 structural descriptors showing that hydrophobicity is a fundamental property in the analyzed biological response.

Sketch design drawing, Hill House, Chapel Hill, Brisbane

Three-dimensional projection sketch showing roof forms and wall finishes.

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.

Stimulated Raman adiabatic passage from an atomic to a molecular Bose-Einstein condensate

The process of stimulated Raman adiabatic passage (STIRAP) provides a possible route for the generation of a coherent molecular Bose-Einstein condensate (BEC) from an atomic BEC. We analyze this process in a three-dimensional mean-field theory, including atom-atom interactions and nonresonant intermediate levels. We find that the process is feasible, but at larger Rabi frequencies than anticipated from a crude single-mode lossless analysis, due to two-photon dephasing caused by the atomic interactions. We then identify optimal strategies in STIRAP allowing one to maintain high conversion efficiencies with smaller Rabi frequencies and under experimentally less demanding conditions.