A combination of implicit and adaptative upwind tools for the numerical solution of incompressible free surface flows

This paper is concerned with the numerical solutions of time dependent two-dimensional incompressible flows. By using the primitive variables of velocity and pressure, the Navier-Stokes and mass conservation equations are solved by a semi-implicit finite difference projection method. A new bounded higher order upwind convection scheme is employed to deal with the non-linear (advective) terms. The procedure is an adaptation of the GENSMAC (J. Comput. Phys. 1994; 110: 171-186) methodology for calculating confined and free surface fluid flows at both low and high Reynolds numbers. The calculations were performed by using the 2D version of the Freeflow simulation system (J. Comp. Visual. Science 2000; 2:199-210). In order to demonstrate the capabilities of the numerical method, various test cases are presented. These are the fully developed flow in a channel, the flow over a backward facing step, the die-swell problem, the broken dam flow, and an impinging jet onto a flat plate. The numerical results compare favourably with the experimental data and the analytical solutions. Copyright (c) 2006 John Wiley & Sons, Ltd.

A simple interferometric method to measure the calibration factor and displacement amplification in piezoelectric flextensional actuators

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

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Differential geometry of intersection curves in R(4) of three implicit surfaces

We present algorithms for computing the differential geometry properties of intersection Curves of three implicit surfaces in R(4), using the implicit function theorem and generalizing the method of X. Ye and T. Maekawa for 4-dimension. We derive t, n, b(1), b(2) vectors and curvatures (k(1), k(2), k(3)) for transversal intersections of the intersection problem. (C) 2008 Elsevier B.V. All rights reserved.

Modeling a spatially restricted distribution in the Neotropics: How the size of calibration area affects the performance of five presence-only methods

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Implicit differential equations with impasse singularities and singular perturbation problems

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

Extension of the Nambu-Jona-Lasinio model predictions at high densities and temperatures using an implicit regularization scheme

Traditional cutoff regularization schemes of the Nambu-Jona-Lasinio model limit the applicability of the model to energy-momentum scales much below the value of the regularizing cutoff. In particular, the model cannot be used to study quark matter with Fermi momenta larger than the cutoff. In the present work, an extension of the model to high temperatures and densities recently proposed by Casalbuoni, Gatto, Nardulli, and Ruggieri is used in connection with an implicit regularization scheme. This is done by making use of scaling relations of the divergent one-loop integrals that relate these integrals at different energy-momentum scales. Fixing the pion decay constant at the chiral symmetry breaking scale in the vacuum, the scaling relations predict a running coupling constant that decreases as the regularization scale increases, implementing in a schematic way the property of asymptotic freedom of quantum chromodynamics. If the regularization scale is allowed to increase with density and temperature, the coupling will decrease with density and temperature, extending in this way the applicability of the model to high densities and temperatures. These results are obtained without specifying an explicit regularization. As an illustration of the formalism, numerical results are obtained for the finite density and finite temperature quark condensate and applied to the problem of color superconductivity at high quark densities and finite temperature.

Determination of jet energy calibration and transverse momentum resolution in CMS

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A generalized alternating-direction implicit scheme for incompressible magnetohydrodynamic viscous flows at low magnetic Reynolds number

This paper presents numerical simulations of incompressible fluid flows in the presence of a magnetic field at low magnetic Reynolds number. The equations governing the flow are the Navier-Stokes equations of fluid motion coupled with Maxwell's equations of electromagnetics. The study of fluid flows under the influence of a magnetic field and with no free electric charges or electric fields is known as magnetohydrodynamics. The magnetohydrodynamics approximation is considered for the formulation of the non-dimensional problem and for the characterization of similarity parameters. A finite-difference technique is used to discretize the equations. In particular, an extension of the generalized Peaceman and Rachford alternating-direction implicit (ADI) scheme for simulating two-dimensional fluid flows is presented. The discretized conservation equations are solved in stream function-vorticity formulation. We compare the ADI and generalized ADI schemes, and show that the latter is more efficient in simulating low Reynolds number and magnetic Reynolds number problems. Numerical results demonstrating the applicability of this technique are also presented. The simulation of incompressible magneto hydrodynamic fluid flows is illustrated by numerical solution for two-dimensional cases. (c) 2007 Elsevier B.V. All rights reserved.

An implicit BEM formulation to model strong discontinuities in solids

A boundary element method (BEM) formulation to predict the behavior of solids exhibiting displacement (strong) discontinuity is presented. In this formulation, the effects of the displacement jump of a discontinuity interface embedded in an internal cell are reproduced by an equivalent strain field over the cell. To compute the stresses, this equivalent strain field is assumed as the inelastic part of the total strain. As a consequence, the non-linear BEM integral equations that result from the proposed approach are similar to those of the implicit BEM based on initial strains. Since discontinuity interfaces can be introduced inside the cell independently on the cell boundaries, the proposed BEM formulation, combined with a tracking scheme to trace the discontinuity path during the analysis, allows for arbitrary discontinuity propagation using a fixed mesh. A simple technique to track the crack path is outlined. This technique is based on the construction of a polygonal line formed by segments inside the cells, in which the assumed failure criterion is reached. Two experimental concrete fracture tests were analyzed to assess the performance of the proposed formulation.

Construction and calibration of permeation tubes.

A dynamic system for stablishing a known standard gas mixture is necessary for many studies such as development and testing of analytical methods. A permeation tube can be used for this purpose. The work described here shows the construction, operation and calibration of a simple permeation tube which can be used to obtain large amounts of a standard gas.

A homological version of the implicit function theorem

Our objective in this paper is to prove an Implicit Function Theorem for general topological spaces. As a consequence, we show that, under certain conditions, the set of the invertible elements of a topological monoid X is an open topological group in X and we use the classical topological group theory to conclude that this set is a Lie group.

Comparison of proton energy loss in thick absorbers in terms of a reduced calibration curve

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Calibration of therapeutic ultrasound equipment

Temporal and spatial acoustic intensity (SATA) of therapeutic ultrasound (US) equipment should be monitored periodically. In order to evaluate the conditions of US equipment in use in the city of Piracicaba-Sao Paulo, Brazil, 31 machines - representing all Brazilian manufacturers - were analysed under continuous and pulsed conditions at a frequency of 1 MHz. Data about temporal and spatial acoustic intensity were collected and the use of equipment was surveyed. Intensities of 0.1, 0.2, 0.5, 0.8, 1.0, 1.5, 2.0, 2.5 and 3.0 Wcm -2, indicated on the equipment panel were analysed using a previously calibrated digital radiation pressure scale, model UPM-DT-1 (Ohmic Instruments Co). The acoustic intensity (I) results were expressed as superior and inferior quartile ranges for transducers with metal surfaces of 9 cm 2 and an effective radiation area (ERA) Of 4 cm 2. The results under continuous conditions were: I 0.1 = -20.0% and -96%. I 0.2 = -3.1% and -83.7%. I 0.5 = -35.0% and -86.5%. I 0.8 = -37.5% and -71.0%. I 2.5 = -49.0% and -69.5%. I 3.0 = -58.1% and -77.6%. For pulsed conditions, intensities were: I 0.1 = -40.0% and -86.2%. I 1.0 = -50.0% and -86.5%. I 1.5 = -62.5% and -82.5%. I 2.0 = -62.5% and -81.6%. I 2.5 = -64.7% and -88.8%. I 3.0 = -87.1% and -94.8%. In reply to the questionnaire drawn up to check the conditions of use of equipment, all users reported the use of hydrosoluble gel as a coupling medium and none had carried out previous calibrations. Most users used intensities in the range of 0.4. to 1.0 Wcm -2 and used machines for 300 to 400 minutes per week. The majority of machines had been bought during the previous seven years and weekly use ranged from less than 100 minutes to 700 minutes (11 hours 40 minutes). Findings confirm previous observations of discrepancy between the intensity indicated on the equipment panel and that emitted by the transducer and highlight the necessity for periodic evaluations of US equipment.

Line based camera calibration in machine vision dynamic applications

The problem of dynamic camera calibration considering moving objects in close range environments using straight lines as references is addressed. A mathematical model for the correspondence of a straight line in the object and image spaces is discussed. This model is based on the equivalence between the vector normal to the interpretation plane in the image space and the vector normal to the rotated interpretation plane in the object space. In order to solve the dynamic camera calibration, Kalman Filtering is applied; an iterative process based on the recursive property of the Kalman Filter is defined, using the sequentially estimated camera orientation parameters to feedback the feature extraction process in the image. For the dynamic case, e.g. an image sequence of a moving object, a state prediction and a covariance matrix for the next instant is obtained using the available estimates and the system model. Filtered state estimates can be computed from these predicted estimates using the Kalman Filtering approach and based on the system model parameters with good quality, for each instant of an image sequence. The proposed approach was tested with simulated and real data. Experiments with real data were carried out in a controlled environment, considering a sequence of images of a moving cube in a linear trajectory over a flat surface.