A C2-continuous high-resolution upwind convection scheme

A bounded upwinding scheme for numerical solution of hyperbolic conservation laws and Navier-Stokes equations is presented. The scheme is based on convection boundedness criterion and total variation diminishing stability criteria and developed by employing continuously differentiable functions. The accuracy of the scheme is verified by assessing the error and observed convergence rate on 1-D benchmark test cases. A comparative study between the new scheme and conventional total variation diminishing/convection boundedness criterion-based upwind schemes to solve standard nonlinear hyperbolic conservation laws is also accomplished. The scheme is then examined in the simulation of Newtonian and non-Newtonian fluid flows of increasing complexity; a satisfactory agreement has been observed in terms of the overall behavior. Finally, the scheme is used to study the hydrodynamics of a gas-solid flow in a bubbling fluidized bed. © 2013 John Wiley & Sons, Ltd.

Bicomplexes and conservation laws in non-Abelian Toda models

A bicomplex structure is associated with the Leznov-Saveliev equation of integrable models. The linear problem associated with the zero-curvature condition is derived in terms of the bicomplex linear equation. The explicit example of a non-Abelian conformal affine Toda model is discussed in detail and its conservation laws are derived from the zero-curvature representation of its equation of motion.

Conservation laws for non-potential operators

A study was developed in order to build a function M invariant in time, by means of Hamiltonian's formulation, taking into account the equation associated to the problem, showing that starting from this function the equation of motion of the system with the contour conditions for non-conservative considered problems can be obtained. The Hamiltonian method is extended for these kind of systems in order to validate for non-potential operators through variational approach.

DE SITTER TRANSITIVITY, CONFORMAL TRANSFORMATIONS AND CONSERVATION LAWS

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

Conservation laws for Z(N) symmetric quantum spin models and their exact ground state energies

We derive an infinite set of conserved charges for some Z(N) symmetric quantum spin models by constructing their Lax pairs. These models correspond to the Potts model, Ashkin-Teller model and the particular set of self-dual Z(N) models solved by Fateev and Zamolodchikov [6]. The exact ground state energy for this last family of hamiltonians is also presented. © 1986.

Leis de Conservação Hiperbólicas 2D com Termo Fonte Stiff

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

Bound state solution of the grassmannian nonlinear sigma model with fermions

We construct the S-matrix for bound state (gauge-invariant) scattering for nonlinear sigma models defined on the manifold SU(n) S(U(p)⊗U(n-p)) with fermions. It is not possible to compute gauge non-singlet matrix elements. In the present language, constraints from higher conservation laws determine the bound state solution. An alternative derivation is also presented. © 1988.

Illegal trade of the guitarfish Rhinobatos horkelii on the coasts of central and southern Brazil: genetic identification to aid conservation

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

Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

A continuously differentiable upwinding scheme for the simulation of fluid flow problems

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

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.

Renormalized new solutions for the massless Thirring model

We present a nonperturbative study of the (1 + 1)-dimensional massless Thirring model by using path integral methods. The regularization ambiguities - coming from the computation of the fermionic determinant - allow to find new solution types for the model. At quantum level the Ward identity for the 1PI 2-point function for the fermionic current separates such solutions in two phases or sectors, the first one has a local gauge symmetry that is implemented at quantum level and the other one without this symmetry. The symmetric phase is a new solution which is unrelated to the previous studies of the model and, in the nonsymmetric phase there are solutions that for some values of the ambiguity parameter are related to well-known solutions of the model. We construct the Schwinger-Dyson equations and the Ward identities. We make a detailed analysis of their UV divergence structure and, after, we perform a nonperturbative regularization and renormalization of the model.

Massless DKP field in a Lyra manifold

Massless scalar and vector fields are coupled to the Lyra geometry by means of the Duffin-Kemmer-Petiau (DKP) theory. Using the Schwinger variational principle, the equations of motion, conservation laws and gauge symmetry are implemented. We find that the scalar field couples to the anholonomic part of the torsion tensor, and the gauge symmetry of the electromagnetic field does not break by the coupling with torsion.

Spinorial field and Lyra geometry

The Dirac field is studied in a Lyra space-time background by means of the classical Schwinger Variational Principle. We obtain the equations of motion, establish the conservation laws, and get a scale relation relating the energy-momentum and spin tensors. Such scale relation is an intrinsic property for matter fields in Lyra background.

Metric-affine approach to teleparallel gravity

The teleparallel gravity theory, treated physically as a gauge theory of translations, naturally represents a particular case of the most general gauge-theoretic model based on the general affine group of spacetime. on the other hand, geometrically, the Weitzenbock spacetime of distant parallelism is a particular case of the general metric-affine spacetime manifold. These physical and geometrical facts offer a new approach to teleparallelism. We present a systematic treatment of teleparallel gravity within the framework of the metric-affine theory. The symmetries, conservation laws and the field equations are consistently derived, and the physical consequences are discussed in detail. We demonstrate that the so-called teleparallel GR-equivalent model has a number of attractive features which distinguishes it among the general teleparallel theories, although it has a consistency problem when dealing with spinning matter sources.