Some experiments with high order compact methods using a computer algebra software - Part I

We consider a procedure for obtaining a compact fourth order method to the steady 2D Navier-Stokes equations in the streamfunction formulation using the computer algebra system Maple. The resulting code is short and from it we obtain the Fortran program for the method. To test the procedure we have solved many cavity-type problems which include one with an analytical solution and the results are compared with results obtained by second order central differences to moderate Reynolds numbers. (c) 2005 Elsevier B.V. All rights reserved.

Some experiments with high order compact methods using a computer algebra software - Part II (non-uniform grid)

We generalize a procedure proposed by Mancera and Hunt [P.F.A. Mancera, R. Hunt, Some experiments with high order compact methods using a computer algebra software-Part 1, Appl. Math. Comput., in press, doi: 10.1016/j.amc.2005.05.015] for obtaining a compact fourth-order method to the steady 2D Navier-Stokes equations in the streamfunction formulation-vorticity using the computer algebra system Maple, which includes conformal mappings and non-uniform grids. To analyse the procedure we have solved a constricted stepped channel problem, where a fine grid is placed near the re-entrant corner by transformation of the independent variables. (c) 2006 Elsevier B.V. All rights reserved.

A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.

Effective field theory methods to model compact binaries

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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.

Star formation in the intragroup medium and other diagnostics of the evolutionary stages of compact groups of galaxies

Context. Compact groups of galaxies are entities that have high densities of galaxies and serve as laboratories to study galaxy interactions, intergalactic star formation and galaxy evolution. Aims. The main goal of this study is to search for young objects in the intragroup medium of seven compact groups of galaxies: HCG 2, 7, 22, 23, 92, 100 and NGC 92 as well as to evaluate the stage of interaction of each group. Methods. We used Fabry-Perot velocity fields and rotation curves together with GALEX NUV and FUV images and optical R-band and HI maps. Results. (i) HCG 7 and HCG 23 are in early stages of interaction; (ii) HCG 2 and HCG 22 are mildly interacting; and (iii) HCG 92, HCG 100 and NGC 92 are in late stages of evolution. We find that all three evolved groups contain populations of young blue objects in the intragroup medium, consistent with ages < 100 Myr, of which several are younger than < 10 Myr. We also report the discovery of a tidal dwarf galaxy candidate in the tail of NGC 92. These three groups, besides containing galaxies that have peculiar velocity fields, also show extended HI tails. Conclusions. Our results indicate that the advanced stage of evolution of a group, together with the presence of intragroup HI clouds, may lead to star formation in the intragroup medium. A table containing all intergalactic HII regions and tidal dwarf galaxies confirmed to date is appended.

Vortex configuration flow cell based on low temperature cofired ceramics as a compact chemiluminescence microsystern

The integration of optical detection methods in continuous flow microsystems can highly extend their range of application, as long as some negative effects derived from their scaling down can be minimized. Downsizing affects to a greater extent the sensitivity of systems based on absorbance measurements than the sensitivity of those based on emission ones. However, a careful design of the instrumental setup is needed to maintain the analytical features in both cases. In this work, we present the construction and evaluation of a simple miniaturized optical system, which integrates a novel flow cell configuration to carry out chemiluminescence (CL) measurements using a simple photodiode. It consists of a micro-mixer based on a vortex structure, which has been constructed by means of the low-temperature cofired ceramics (LTCC) technology. This mixer not only efficiently promotes the CL reaction due to the generated high turbulence but also allows the detection to be carried out in the same area, avoiding intensity signal losses. As a demonstration, a flow injection system has been designed and optimized for the detection of cobalt(H) in water samples. It shows a linear response between 2 and 20 mu M with a correlation of r > 0.993, a limit of detection of 1.1 mu M, a repeatability of RSD = 12.4 %, and an analysis time of 17 s. These results demonstrate the suitability of the proposal to the determination of compounds involved in CL reactions by means of an easily constructed versatile device based on low-cost instrumentation.

Compact clinical MRI magnet design using a multilayer current density approach

In this work, a new method of optimization is successfully applied to the theoretical design of compact, actively shielded, clinical MRI magnets. The problem is formulated as a two-step process in which the desired current densities on multiple, cc-axial surface layers are first calculated by solving Fredholm equations of the first kind. Non-linear optimization methods with inequality constraints are then invoked to fit practical magnet coils to the desired current densities. The current density approach allows rapid prototyping of unusual magnet designs. The emphasis of this work is on the optimal design of short, actively-shielded MRI magnets for whole-body imaging. Details of the hybrid numerical model are presented, and the model is used to investigate compact, symmetric, and asymmetric MRI magnets. Magnet designs are presented for actively-shielded, symmetric magnets of coil length 1.0 m, which is considerably shorter than currently available designs of comparable dsv size. Novel, actively-shielded, asymmetric magnet designs are also presented in which the beginning of a 50-cm dsv is positioned just 11 cm from the end of the coil structure, allowing much improved access to the patient and reduced patient claustrophobia. Magn Reson Med 45:331540, 2001. (C) 2001 Wiley-Liss, Inc.

Compact SMB chromatography for binary separation

A thesis submitted for the degree of Doctor of Philosophy

Compact ion-source based on superionic rubidium silver iodide (RbAg4I5) solid electrolyte

Dissertação para obtenção do Grau de Mestre em Engenharia Física

Developing discontinuous Galerkin methods for the resolution of the incompressible Navier-Stokes equations

Informe de investigación elaborado a partir de una estancia en el Laboratorio de Diseño Computacional en Aeroespacial en el Massachusetts Institute of Technology (MIT), Estados Unidos, entre noviembre de 2006 y agosto de 2007. La aerodinámica es una rama de la dinámica de fluidos referida al estudio de los movimientos de los líquidos o gases, cuya meta principal es predecir las fuerzas aerodinámicas en un avión o cualquier tipo de vehículo, incluyendo los automóviles. Las ecuaciones de Navier-Stokes representan un estado dinámico del equilibrio de las fuerzas que actúan en cualquier región dada del fluido. Son uno de los sistemas de ecuaciones más útiles porque describen la física de una gran cantidad de fenómenos como corrientes del océano, flujos alrededor de una superficie de sustentación, etc. En el contexto de una tesis doctoral, se está estudiando un flujo viscoso e incompresible, solucionando las ecuaciones de Navier- Stokes incompresibles de una manera eficiente. Durante la estancia en el MIT, se ha utilizado un método de Galerkin discontinuo para solucionar las ecuaciones de Navier-Stokes incompresibles usando, o bien un parámetro de penalti para asegurar la continuidad de los flujos entre elementos, o bien un método de Galerkin discontinuo compacto. Ambos métodos han dado buenos resultados y varios ejemplos numéricos se han simulado para validar el buen comportamiento de los métodos desarrollados. También se han estudiado elementos particulares, los elementos de Raviart y Thomas, que se podrían utilizar en una formulación mixta para obtener un algoritmo eficiente para solucionar problemas numéricos complejos.

Stabilized Petrov-Galerkin methods for the convection-diffusion-reaction and the Helmholtz equations

We present two new stabilized high-resolution numerical methods for the convection–diffusion–reaction (CDR) and the Helmholtz equations respectively. The work embarks upon a priori analysis of some consistency recovery procedures for some stabilization methods belonging to the Petrov–Galerkin framework. It was found that the use of some standard practices (e.g. M-Matrices theory) for the design of essentially non-oscillatory numerical methods is not feasible when consistency recovery methods are employed. Hence, with respect to convective stabilization, such recovery methods are not preferred. Next, we present the design of a high-resolution Petrov–Galerkin (HRPG) method for the 1D CDR problem. The problem is studied from a fresh point of view, including practical implications on the formulation of the maximum principle, M-Matrices theory, monotonicity and total variation diminishing (TVD) finite volume schemes. The current method is next in line to earlier methods that may be viewed as an upwinding plus a discontinuity-capturing operator. Finally, some remarks are made on the extension of the HRPG method to multidimensions. Next, we present a new numerical scheme for the Helmholtz equation resulting in quasi-exact solutions. The focus is on the approximation of the solution to the Helmholtz equation in the interior of the domain using compact stencils. Piecewise linear/bilinear polynomial interpolation are considered on a structured mesh/grid. The only a priori requirement is to provide a mesh/grid resolution of at least eight elements per wavelength. No stabilization parameters are involved in the definition of the scheme. The scheme consists of taking the average of the equation stencils obtained by the standard Galerkin finite element method and the classical finite difference method. Dispersion analysis in 1D and 2D illustrate the quasi-exact properties of this scheme. Finally, some remarks are made on the extension of the scheme to unstructured meshes by designing a method within the Petrov–Galerkin framework.

Compact matrix expressions for generalized Wald tests of equality of moment vectors

Asymptotic chi-squared test statistics for testing the equality ofmoment vectors are developed. The test statistics proposed aregeneralizedWald test statistics that specialize for different settings by inserting andappropriate asymptotic variance matrix of sample moments. Scaled teststatisticsare also considered for dealing with situations of non-iid sampling. Thespecializationwill be carried out for testing the equality of multinomial populations, andtheequality of variance and correlation matrices for both normal andnon-normaldata. When testing the equality of correlation matrices, a scaled versionofthe normal theory chi-squared statistic is proven to be an asymptoticallyexactchi-squared statistic in the case of elliptical data.