Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. And third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,

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.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.

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.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Comparative analysis of the gas-phase reactions of cylindrospermopsin and the difference in the alkali metal cation mobility

Cylindrospermopsin (CYN) belongs to a group of toxins produced by several strains of freshwater cyanobacteria. It is a compact zwitterionic molecule composed of a uracil section and a tricyclic guanidinium portion with a primarily hepatotoxic effect. Using low multi-stage and high-resolution mass spectrometry, the gas-phase reactions of this toxin have been investigated. Our data show that collision-induced dissociation (CID) spectra of CYN are dominated by neutral losses, and three major initial fragmentation pathways are clearly distinguishable. Interestingly, comparative analysis of protonated and cationizated molecules showed a significant difference in the balance of the SO(3) and terminal ring elimination. These data indicate that the differential ion mobility of H(+), Li(+), Na(+) and K(+) leads to different fragmentation pathways, giving rise to mass spectra with different profiles. Copyright (C) 2008 John Wiley & Sons, Ltd.

A Direct Numerov Sixth-order Numerical Scheme to Accurately Solve the Unidimensional Poisson Equation with Dirichlet Boundary Conditions

In this article, we present an analytical direct method, based on a Numerov three-point scheme, which is sixth order accurate and has a linear execution time on the grid dimension, to solve the discrete one-dimensional Poisson equation with Dirichlet boundary conditions. Our results should improve numerical codes used mainly in self-consistent calculations in solid state physics.

A compact and high-resolution version of a capacitively coupled contactless conductivity detector

An all-in-one version of a capacitively coupled contactless conductivity detector is introduced. The absence of moving parts (potentiometers and connectors) makes it compact (6.5 cm(3)) and robust. A local oscillator, working at 1.1 MHz, was optimized to use capillaries of id from 20 to 100 lam. Low noise circuitry and a high-resolution analog-to-digital converter (ADC) (21 bits effective) grant good sensitivities for capillaries and background electrolytes currently used in capillary electrophoresis. The fixed frequency and amplitude of the signal generator is a drawback that is compensated by the steady calibration curves for conductivity. Another advantage is the possibility of determining the inner diameter of a capillary by reading the ADC when air and subsequently water flow through the capillary. The difference of ADC reading may be converted into the inner diameter by a calibration curve. This feature is granted by the 21-bit ADC, which eliminates the necessity of baseline compensation by hardware. In a typical application, the limits of detection based on the 3 sigma criterion (without baseline filtering) were 0.6, 0.4, 0.3, 0.5, 0.6, and 0.8 mu mol/L for K(+), Ba(2+), Ca(2+), Na(+), Mg(2+), and Li(+), respectively, which is comparable to other high-quality implementations of a capacitively coupled contactless conductivity detector.

Correlation between metal-ceramic bond strength and coefficient of linear thermal expansion difference

The purpose of this study was to evaluate the metal-ceramic bond strength (MCBS) of 6 metal-ceramic pairs (2 Ni-Cr alloys and 1 Pd-Ag alloy with 2 dental ceramics) and correlate the MCBS values with the differences between the coefficients of linear thermal expansion (CTEs) of the metals and ceramics. Verabond (VB) Ni-Cr-Be alloy, Verabond II (VB2), Ni-Cr alloy, Pors-on 4 (P), Pd-Ag alloy, and IPS (I) and Duceram (D) ceramics were used for the MCBS test and dilatometric test. Forty-eight ceramic rings were built around metallic rods (3.0 mm in diameter and 70.0 mm in length) made from the evaluated alloys. The rods were subsequently embedded in gypsum cast in order to perform a tensile load test, which enabled calculating the CMBS. Five specimens (2.0 mm in diameter and 12.0 mm in length) of each material were made for the dilatometric test. The chromel-alumel thermocouple required for the test was welded into the metal test specimens and inserted into the ceramics. ANOVA and Tukey's test revealed significant differences (p=0.01) for the MCBS test results (MPa), with PI showing higher MCBS (67.72) than the other pairs, which did not present any significant differences. The CTE (10-6 oC-1) differences were: VBI (0.54), VBD (1.33), VB2I (-0.14), VB2D (0.63), PI (1.84) and PD (2.62). Pearson's correlation test (r=0.17) was performed to evaluate of correlation between MCBS and CTE differences. Within the limitations of this study and based on the obtained results, there was no correlation between MCBS and CTE differences for the evaluated metal-ceramic pairs.

An Environmental Friendly Procedure for Photometric Determination of Hypochlorite in Tap Water Employing a Miniaturized Multicommuted Flow Analysis Setup

A photometric procedure for the determination of ClO(-) in tap water employing a miniaturized multicommuted flow analysis setup and an LED-based photometer is described. The analytical procedure was implemented using leucocrystal violet (LCV; 4,4', 4 ''-methylidynetris (N, N-dimethylaniline), C(25)H(31)N(3)) as a chromogenic reagent. Solenoid micropumps employed for solutions propelling were assembled together with the photometer in order to compose a compact unit of small dimensions. After control variables optimization, the system was applied for the determination of ClO(-) in samples of tap water, and aiming accuracy assessment samples were also analyzed using an independent method. Applying the paired t-test between results obtained using both methods, no significant difference at the 95% confidence level was observed. Other useful features include low reagent consumption, 2.4 mu g of LCV per determination, a linear response ranging from 0.02 up to 2.0 mg L(-1) ClO(-), a relative standard deviation of 1.0% (n = 11) for samples containing 0.2 mg L(-1) ClO(-), a detection limit of 6.0 mu g L(-1) ClO(-), a sampling throughput of 84 determinations per hour, and a waste generation of 432 mu L per determination.

Desire, Hierarchy, and Agency: Youth, Homosexuality, and Difference Markers in Sao Paulo

In this article, we present the results of ethnographic research undertaken in areas with concentrated recreational venues frequented by diverse youth groups in the city of Sao Paulo, Brazil. This material is part of the larger ""Relationships between race, gender and sexuality in different national and local contexts"" research initiative. Here, we focus upon data collected through ethnographic observation and interviews conducted in various locales where young men meet to engage in homoerotic sociability, demonstrating different degrees of conviviality among groups of different socioeconomic profiles and distinct esthetic preferences, consumer habits, and body types. We explore the production of styles and body presentations that link markers of color/race, gender, and sexuality, as well as the relationships that these maintain with opening or restricting possibilities for the establishment of erotic and affective partnerships involving these boys.

Determination of Three-Dimensional Left Ventricle Motion to Analyze Ventricular Dyssyncrony in SPECT Images

A method to compute three-dimension (3D) left ventricle (LV) motion and its color coded visualization scheme for the qualitative analysis in SPECT images is proposed. It is used to investigate some aspects of Cardiac Resynchronization Therapy (CRT). The method was applied to 3D gated-SPECT images sets from normal subjects and patients with severe Idiopathic Heart Failure, before and after CRT. Color coded visualization maps representing the LV regional motion showed significant difference between patients and normal subjects. Moreover, they indicated a difference between the two groups. Numerical results of regional mean values representing the intensity and direction of movement in radial direction are presented. A difference of one order of magnitude in the intensity of the movement on patients in relation to the normal subjects was observed. Quantitative and qualitative parameters gave good indications of potential application of the technique to diagnosis and follow up of patients submitted to CRT.

Searching for star formation outside galaxies: Multiwavelength analysis of the intragroup medium of Hickson Compact Group 100

We use multiwavelength data (H I, FUV, NUV, R) to search for evidence of star formation in the intragroup medium of the Hickson Compact Group 100. We find that young star-forming regions are located in the intergalactic H I clouds of the compact group which extend to over 130 kpc away from the main galaxies. A tidal dwarf galaxy (TDG) candidate is located in the densest region of the H I tail, 61 kpc from the brightest group member and its age is estimated to be only 3.3 Myr. Fifteen other intragroup H II regions and TDG candidates are detected in the Galaxy Evolution Explorer (GALEX) FUV image and within a field 10' x 10' encompassing the H I tail. They have ages <200 Myr, H I masses of 10(9.2-10.4) M(circle dot), 0.001

Equation of state for the magnetic-color-flavor-locked phase and its implications for compact star models

Using the solutions of the gap equations of the magnetic-color-flavor-locked (MCFL) phase of paired quark matter in a magnetic field, and taking into consideration the separation between the longitudinal and transverse pressures due to the field-induced breaking of the spatial rotational symmetry, the equation of state of the MCFL phase is self-consistently determined. This result is then used to investigate the possibility of absolute stability, which turns out to require a field-dependent ""bag constant"" to hold. That is, only if the bag constant varies with the magnetic field, there exists a window in the magnetic field vs bag constant plane for absolute stability of strange matter. Implications for stellar models of magnetized (self-bound) strange stars and hybrid (MCFL core) stars are calculated and discussed.

Two formation channels of ultra-compact dwarf galaxies in Hickson compact groups

Context. The formation of ultra-compact dwarf galaxies (UCDs) is believed to be driven by interaction, and UCDs are abundant in the cores of galaxy clusters, environments that mark the end-point of galaxy evolution. Nothing is known about the properties of UCDs in compact groups of galaxies, environments where most of galaxy evolution and interaction is believed to occur and where UCDs in an intermediate stage in their evolution may be expected. Aims. The main goal of this study is to detect and characterize, for the first time, the UCD population of compact groups of galaxies. For that, two nearby groups in different evolutionary stages, HCG22 and HCG90, were targeted. Methods. We selected about 40 UCD candidates from pre-existing photometry of both groups, and obtained spectra of these candidates using the VLT FORS2 instrument in MXU mode. Archival HST/ACS imaging was used to measure their structural parameters. Results. We detect 16 and 5 objects belonging to HCG22 and HCG90, respectively, covering the magnitude range -10.0 > M(R) > -11.5 mag. Their integrated colours are consistent with old ages covering a broad range in metallicities (metallicities confirmed by the spectroscopic measurements). Photometric mass estimates put 4 objects in HCG90 and 9 in HCG22 in the mass range of UCDs (> 2 x 10(6) M(circle dot)) for an assumed age of 12Gyr. These UCDs are on average 2-3 times larger than the typical size of Galactic GCs, covering a range of 2 less than or similar to r(h) less than or similar to 21 pc. The UCDs in HCG22 are more concentrated around the central galaxy than in HCG90, at the 99% confidence level. They cover a broad range in [alpha/Fe] abundances from sub-to super-solar. The spectra of 3 UCDs (2 in HCG22, 1 in HCG90) show tentative evidence of intermediate age stellar populations. The clearest example is the largest and most massive UCD (similar to 10(7) M(circle dot)) in our sample, which is detected in HCG22. Its properties are most consistent with a stripped dwarf galaxy nucleus. We calculate the specific frequency (S(N)) of UCDs for both groups, finding that HCG22 has about three times higher S(N) than HCG90. Conclusions. The ensemble properties of the detected UCDs supports two co-existing formation channels: a star cluster origin (low-luminosity, compact sizes, old ages, super-solar alpha/Fe), and an origin as tidally stripped dwarf nuclei (more extended and younger stellar populations). Our results imply that the UCDs detected in both groups do not, in their majority, originate from relatively recent galaxy interactions. Most of the detected UCDs have likely been brought into the group along with their host galaxies.