Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

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.

Criteria for mixed grids in computational fluid dynamics

The finite volume method is used as a numerical method for solving the fluid flow equations. This method is appropriate to employ under structured and unstructured meshes. Mixed grids, combining both types of grids, are investigated. The coupling of different grids is done by overlapping strategy. The computational effort for the mixed grid is evaluated by the CPU-time, with different percentage of covering area of the unstructured mesh. The present scheme is tested for the driven cavity problem, where the incompressible fluid is integrated by calculating the velocity fields and computing the pressure field in each time step. Several schemes for unstructured grid are examined, and the compatibility condition is applied to check their consistency. A scheme to verify the compatibility condition for the unstructured grids is presented. (c) 2006 IMACS. Published by Elsevier B.V. All rights reserved.

Evolution of magnetic field and spin period in accreting neutron stars

Based on the accretion-induced magnetic field decay model, in which a frozen field and an incompressible fluid are assumed, we obtain the following results: (1) an analytic relation between the magnetic field and spin period, if the fastness parameter of the accretion disk is neglected (The evolutionary tracks of accreting neutron stars in the P-B diagram in our model are different from the equilibrium period lines when the influence of the fastness parameter is taken into account.); (2) the theoretical minimum spin period of an accreting neutron star is max(1.1ms (DeltaM/M(circle dot))(-1)R(6)(-5/14) I(45)(M/M(circle dot))(-1/2),1.1ms (M/M(circle dot))(-1/2) R(6)(17/14)), independent of the accretion rate (X-ray luminosity) but dependent on the total accretion mass, DeltaM; however, the minimum magnetic field depends on the accretion rate; (3) the magnetic field strength decreases faster with time than does the period.

Simulação numérica de escoamentos de fluidos com efeitos de rotação

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

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos baseado em volumes de controle em malhas não estruturadas

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

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos de mínimos quadrados

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

Implementação do método das características na modelagem de problemas de convecção natural em cavidades cilíndricas

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

Modelagem do fluxo sanguíneo na aorta abdominal utilizando interação fluido-estrutura

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

Desenvolvimento de uma metodologia numérica para escoamentos viscoelásticos não-isotérmicos

Pós-graduação em Matematica Aplicada e Computacional - FCT

Solução das equações de Navier-Stokes e da equação de transporte por um método de elementos finitos estabilizado pela técnica de separação baseada na característica

Pós-graduação em Engenharia Mecânica - FEIS

Analytical solution for transient one-dimensional couette flow considering constant and time-dependent pressure gradients

This paperaims to determine the velocity profile, in transient state, for a parallel incompressible flow known as Couette flow. The Navier-Stokes equations were applied upon this flow. Analytical solutions, based in Fourier series and integral transforms, were obtained for the one-dimensional transient Couette flow, taking into account constant and time-dependent pressure gradients acting on the fluid since the same instant when the plate starts it´s movement. Taking advantage of the orthogonality and superposition properties solutions were foundfor both considered cases. Considering a time-dependent pressure gradient, it was found a general solution for the Couette flow for a particular time function. It was found that the solution for a time-dependent pressure gradient includes the solutions for a zero pressure gradient and for a constant pressure gradient.

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 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.

LEAST ACTION PRINCIPLE and THE INCOMPRESSIBLE EULER EQUATIONS WITH VARIABLE DENSITY

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