Evaluation of a polynya flux model by means of thermal infrared satellite estimates

We test the ability of a two-dimensional flux model to simulate polynya events with narrow open-water zones by comparing model results to ice-thickness and ice-production estimates derived from thermal infrared Moderate Resolution Imaging Spectroradiometer (MODIS) observations in conjunction with an atmospheric dataset. Given a polynya boundary and an atmospheric dataset, the model correctly reproduces the shape of an 11 day long event, using only a few simple conservation laws. Ice production is slightly overestimated by the model, owing to an underestimated ice thickness. We achieved best model results with the consolidation thickness parameterization developed by Biggs and others (2000). Observed regional discrepancies between model and satellite estimates might be a consequence of the missing representation of the dynamic of the thin-ice thickening (e.g. rafting). We conclude that this simplified polynya model is a valuable tool for studying polynya dynamics and estimating associated fluxes of single polynya events.

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,

The concept of quasi-integrability: a concrete example

We use the deformed sine-Gordon models recently presented by Bazeia et al [1] to take the first steps towards defining the concept of quasi-integrability. We consider one such definition and use it to calculate an infinite number of quasi-conserved quantities through a modification of the usual techniques of integrable field theories. Performing an expansion around the sine-Gordon theory we are able to evaluate the charges and the anomalies of their conservation laws in a perturbative power series in a small parameter which describes the ""closeness"" to the integrable sine-Gordon model. We show that in the case of the two-soliton scattering the charges, up to first order of perturbation, are conserved asymptotically, i.e. their values are the same in the distant past and future, when the solitons are well separated. We indicate that this property may hold or not to higher orders depending on the behavior of the two-soliton solution under a special parity transformation. For closely bound systems, such as breather-like field configurations, the situation however is more complex and perhaps the anomalies have a different structure implying that the concept of quasi-integrability does not apply in the same way as in the scattering of solitons. We back up our results with the data of many numerical simulations which also demonstrate the existence of long lived breather-like and wobble-like states in these models.

Exact vortex solutions in a CP(N) Skyrme-Faddeev type model

We consider a four dimensional field theory with target space being CP(N) which constitutes a generalization of the usual Skyrme-Faddeev model defined on CP(1). We show that it possesses an integrable sector presenting an infinite number of local conservation laws, which are associated to the hidden symmetries of the zero curvature representation of the theory in loop space. We construct an infinite class of exact solutions for that integrable submodel where the fields are meromorphic functions of the combinations (x(1) + i x(2)) and (x(3) + x(0)) of the Cartesian coordinates of four dimensional Minkowski space-time. Among those solutions we have static vortices and also vortices with waves traveling along them with the speed of light. The energy per unity of length of the vortices show an interesting and intricate interaction among the vortices and waves.

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.

Conserved currents in gravitational models with quasi-invariant Lagrangians: Application to teleparallel gravity

Conservation laws in gravitational theories with diffeomorphism and local Lorentz symmetry are studied. Main attention is paid to the construction of conserved currents and charges associated with an arbitrary vector field that generates a diffeomorphism on the spacetime. We further generalize previous results for the case of gravitational models described by quasi-invariant Lagrangians, that is, Lagrangians that change by a total derivative under the action of the local Lorentz group. The general formalism is then applied to the teleparallel models, for which the energy and the angular momentum of a Kerr black hole are calculated. The subsequent analysis of the results obtained demonstrates the importance of the choice of the frame.

Comments on regularization ambiguities and local gauge symmetries

We study the regularization ambiguities in an exact renormalized (1 + 1)-dimensional field theory. We show a relation between the regularization ambiguities and the coupling parameters of the theory as well as their role in the implementation of a local gauge symmetry at quantum level.

Loop scattering in two-dimensional QCD

Using the integrability conditions that we recently obtained in two-dimensional QCD with massless fermions we arrive at a sufficient number of conservation laws to fix the scattering amplitudes involving a local version of the Wilson loop operator.

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.

Extended gauge theories

A scheme inspired in Lie algebra extensions is introduced that enlarges gauge models to allow some coupling between space-time and gauge space. Everything may be written in terms of a generalized covariant derivative including usual differential plus purely algebraic terms. A noncovariant vacuum appears, introducing a natural symmetry breaking, but currents satisfy conservation laws alike those found in gauge theories. © 1991 American Institute of Physics.