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.

Simulation results and applications of an advection bounded scheme to practical flows

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

Charges and fluxes in Maxwell theory on compact manifolds with boundary

We investigate the charges and fluxes that can occur in higher-order Abelian gauge theories defined on compact space-time manifolds with boundary. The boundary is necessary to supply a destination to the electric lines of force emanating from brane sources, thus allowing non-zero net electric charges, but it also introduces new types of electric and magnetic flux. The resulting structure of currents, charges, and fluxes is studied and expressed in the language of relative homology and de Rham cohomology and the corresponding abelian groups. These can be organised in terms of a pair of exact sequences related by the Poincare-Lefschetz isomorphism and by a weaker flip symmetry exchanging the ends of the sequences. It is shown how all this structure is brought into play by the imposition of the appropriately generalised Maxwell's equations. The requirement that these equations be integrable restricts the world-volume of a permitted brane (assumed closed) to be homologous to a cycle on the boundary of space-time. All electric charges and magnetic fluxes are quantised and satisfy the Dirac quantisation condition. But through some boundary cycles there may be unquantised electric fluxes associated with quantised magnetic fluxes and so dyonic in nature.

On the Baker-Akhiezer function in the AKNS scheme

Expressions for the Baker-Akhiezer function and their logarithmic space and time derivatives are derived in terms of the matrix elements of U - V matrices and 'squared basis functions'. These expressions generalize the well known formulas for the KdV equation case and establish links between different forms of the Whitham averaging procedure.

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.

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.

A label-free immunosensor based on recordable compact disk chip for early diagnostic of the dengue virus infection

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

Compact disks, a new source for gold electrodes. Application to the quantification of copper by PSA

In this work we describe a versatile and very sensitive way for copper quantification by potentiometric stripping analysis using gold electrodes obtained from recordable compact disks (CDs). This new source of electrodes (CDtrodes) shown similar performance to the commercial gold electrodes with superior versatility and lower cost. Recordable CDs contains a highly pure gold film with thickness between 50 and 100 nm and superficial area of ca. 100 cm(2). The working electrode developed was used successfully in stationary cell and many experimental parameters have been optimized. For copper, the detection limit attained was 30 ng L-1 (600 s deposition time) with remarkable precision (standard deviation of 1.8 % for 20 repetitive measurements using 25 mu gL(-1) of copper with 60 s of deposition time). The gold electrode developed was used for analysis of copper in sugar cane spirits and tap water samples. The results were compared with those obtained by atomic absorption spectroscopy.

A critical investigation of the Tanford-Kirkwood scheme by means of Monte Carlo simulations

Monte Carlo simulations are used to assess the adequacy of the Tanford-Kirkwood prescription for electrostatic interactions in macromolecules. Within a continuum dielectric framework, the approach accurately describes salt screening of electrostatic interactions for moderately charged systems consistent with common proteins at physiological conditions. The limitations of the Debye-Huckel theory, which forms the statistical mechanical basis for the Tanford-Kirkwood result, become apparent for highly charged systems. It is shown, both by an analysis of the Debye-Huckel theory and by numerical simulations, that the difference in dielectric permittivity between macromolecule and surrounding solvent does not play a significant role for salt effects if the macromolecule is highly charged. By comparison to experimental data, the continuum dielectric model (combined with either an approximate effective Hamiltonian as in the Tanford-Kirkwood treatment or with exact Monte Carlo simulations) satisfactorily predicts the effects of charge mutation on metal ion binding constants, but only if the macromolecule and solvent are assigned the same or similar permittivities.

A case study of compact cogeneration using various fuels

A method for studying the technical and economic feasibility of absorption refrigeration systems in compact cogenerators is presented. The system studied consists of an internal combustion engine, an electric generator and a heat exchanger to recover residual heat from the refrigeration water and exhaust gases. As an application, a computer program simulates the cogeneration system in a building which already has 75 kW of installed electric power. The maximum electric and refrigeration demands are 45 kW and 76 kW respectively. This study simulates the system performance, utilizing diesel oil, sugar cane alcohol and natural gas as possible fuels. (C) 1997 Elsevier B.V. Ltd.

Comment on Casimir force in compact non-commutative extra dimensions and radius stabilization

We call attention to a series of mistakes in a paper by S. Nam recently published in this journal (J. High Energy Phys. 10 (2000) 044).

Exact dynamics of single-qubit-gate fidelities under the measurement-based quantum computation scheme

Measurement-based quantum computation is an efficient model to perform universal computation. Nevertheless, theoretical questions have been raised, mainly with respect to realistic noise conditions. In order to shed some light on this issue, we evaluate the exact dynamics of some single-qubit-gate fidelities using the measurement-based quantum computation scheme when the qubits which are used as a resource interact with a common dephasing environment. We report a necessary condition for the fidelity dynamics of a general pure N-qubit state, interacting with this type of error channel, to present an oscillatory behavior, and we show that for the initial canonical cluster state, the fidelity oscillates as a function of time. This state fidelity oscillatory behavior brings significant variations to the values of the computational results of a generic gate acting on that state depending on the instants we choose to apply our set of projective measurements. As we shall see, considering some specific gates that are frequently found in the literature, the fast application of the set of projective measurements does not necessarily imply high gate fidelity, and likewise the slow application thereof does not necessarily imply low gate fidelity. Our condition for the occurrence of the fidelity oscillatory behavior shows that the oscillation presented by the cluster state is due exclusively to its initial geometry. Other states that can be used as resources for measurement-based quantum computation can present the same initial geometrical condition. Therefore, it is very important for the present scheme to know when the fidelity of a particular resource state will oscillate in time and, if this is the case, what are the best times to perform the measurements.

Feynman integrals with tensorial structure in the negative dimensional integration scheme

The negative-dimensional integration method (NDIM) is revealing itself as a very useful technique for computing massless and/or massive Feynman integrals, covariant and noncovanant alike. Up until now however, the illustrative calculations done using such method have been mostly covariant scalar integrals/without numerator factors. We show here how those integrals with tensorial structures also can be handled straightforwardly and easily. However, contrary to the absence of significant features in the usual approach, here the NDIM also allows us to come across surprising unsuspected bonuses. Toward this end, we present two alternative ways of working out the integrals and illustrate them by taking the easiest Feynman integrals in this category that emerge in the computation of a standard one-loop self-energy diagram. One of the novel and heretofore unsuspected bonuses is that there are degeneracies in the way one can express the final result for the referred Feynman integral.