Assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems

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

LEAST ACTION PRINCIPLE and THE INCOMPRESSIBLE EULER EQUATIONS WITH VARIABLE DENSITY

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

On the regularity of the pressure field of relaxed solutions to Euler equations with variable density

In this paper we improve the regularity in time of the gradient of the pressure field in the solution of relaxed version of variational formulation proposed by V. I. Arnold and by Y. Brenier, for the incompressible Euler equations with variable density. We obtain that the pressure field is not only a measure, but a function in Lloc2((0,T);BVloc(D)) as an extension of the work of Ambrosio and Figalli (2008) in [1] to the variable density case. © 2013 Elsevier Ltd.

A Family of Stationary Solutions to the Euler Equations and Generalized Solutions

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

Dynamics at infinity and other global dynamical aspects of Shimizu-Morioka equations

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

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.

Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Linearizability of the perturbed Burgers equation

We show in this report that the perturbed Burgers equation ut = 2uux + uxx + ε(3 α1u2ux + 3 α2uuxx + 3 α3u2 x + α4uxxx) is equivalent, through a near-identity transformation and up to O(ε), to a linearizable equation if the condition 3 α1 - 3 α3 - 3/2α2 + 3/2α4 = 0 is satisfied. In the case this condition is not fulfilled, a normal form for the equation under consideration is given. We show, furthermore, that nonlinearizable cases lead to perturbative expansions with secular-type behavior. Then, to illustrate our results, we make a linearizability analysis of the equations governing the dynamics of a one-dimensional gas.

Object oriented modeling of piston engines

This work reports a conception phase of a piston engine global model. The model objective is forecast the motor performance (power, torque and specific consumption as a function of rotation and environmental conditions). Global model or Zero-dimensional is based on flux balance through each engine component. The resulting differential equations represents a compressive unsteady flow, in which, all dimensional variables are areas or volumes. A review is presented first. The ordinary differential equation system is presented and a Runge-Kutta method is proposed to solve it numerically. The model includes the momentum conservation equation to link the gas dynamics with the engine moving parts rigid body mechanics. As an oriented to objects model the documentation follows the UML standard. A discussion about the class diagrams is presented, relating the classes with physical model related. The OOP approach allows evolution from simple models to most complex ones without total code rewrite. Copyright © 2001 Society of Automotive Engineers, Inc.

Array of Bose-Einstein condensates under time-periodic Feshbach-resonance management

The investigation of the dynamics of a discrete soliton in an array of Bose-Einstein condensates under the action of a periodically time-modulated atomic scattering length [Feshbach-resonance management (FRM)] was discussed. The slow and rapid modulations, in comparison with the tunneling frequency were considered. An averaged equation, which was a generalized discrete nonlinear Schrödinger equation, including higher-order effective nonlinearities and intersite nonlinear interactions was derived in the case of the rapid modulation. It was demonstrated that the modulations of sufficient strength results in splitting of the soliton by direct simulations.

Dynamics of the deconfinement transition

We estimate the dissipation coefficient Γ that appears in Ginzburg-Landau-Langevin equations that describe phenomenologically the deconfinement transition in QCD. This is done through the implementation of Glauber dynamics of pure SU(3) lattice gauge theory. The coefficient Γ is extracted from the short-time exponential growth of the equal time correlation function of the order parameter. Although the absolute determination of Γ is ambiguous due to the difficulties in relating real time and Monte Carlo time, its relative temperature dependence can be obtained with much less arbitrariness. © 2007 American Institute of Physics.

Equações de estado exóticas na equação de Einstein

Pós-graduação em Física - IFT

Semi-analytical study of the rotational motion stability of artificial satellites using quaternions

This study at aims performing the stability analysis of the rotational motion to artificial satellites using quaternions to describe the satellite attitude (orientation on the space). In the system of rotational motion equations, which is composed by four kinematic equations of the quaternions and by the three Euler equations in terms of the rotational spin components. The influence of the gravity gradient and the direct solar radiation pressure torques have been considered. Equilibrium points were obtained through numerical simulations using the softwares Matlab and Octave, which are then analyzed by the Routh-Hurwitz Stability Criterion.

An integrable evolution equation for surface waves in deep water

In order to describe the dynamics of monochromatic surface waves in deep water, we derive a nonlinear and dispersive system of equations for the free surface elevation and the free surface velocity from the Euler equations in infinite depth. From it, and using a multiscale perturbative method, an asymptotic model for small wave steepness ratio is derived. The model is shown to be completely integrable. The Lax pair, the first conserved quantities as well as the symmetries are exhibited. Theoretical and numerical studies reveal that it supports periodic progressive Stokes waves which peak and break in finite time. Comparison between the limiting wave solution of the asymptotic model and classical results is performed.

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

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