Número do tipo ponto flutuante com precisão estendida

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

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

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

NUMERICAL SIMULATION of CONVECTION-DIFFUSION PROBLEMS BY THE CONTROL-VOLUME-BASED FINITE-ELEMENT METHOD

This work presents a numerical study of the tri-dimensional convection-diffusion equation by the control-volume-based on finite-element method using quadratic hexahedral elements. Considering that the equation governing this problem in its main variable may represent several properties, including temperature, turbulent kinetic energy, viscous dissipation rate of the turbulent kinetic energy, specific dissipation rate of the turbulent kinetic energy, or even the concentration of a contaminant in a given medium, among others, the wide applicability of this problem is thus evidenced. Three cases of temperature distributions will be studied specifically in this work, in addition to one case of pollutant dispersion upon analysis of the concentration of a contaminant in a fixed flow point. Some comparisons will be carried out against works found in the open literature, while others will be done according to each phenomenon characteristics.

ITERATIVE NUMERICAL-SOLUTION OF SCATTERING PROBLEMS

An iterative Neumann series method, employing a real auxiliary scattering integral equation, is used to calculate scattering lengths and phase shifts for the atomic Yukawa and exponential potentials. For these potentials the original Neumann series diverges. The present iterative method yields results that are far better, in convergence, stability and precision, than other momentum space methods. Accurate result is obtained in both cases with an estimated error of about 1 in 10(10) after some 8-10 iterations.

A bounded upwinding scheme for computing convection-dominated transport problems

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

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

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

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)

Predictor-Corrector Primal-Dual Interior Point Method for Solving Economic Dispatch Problems: A Postoptimization Analysis

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

The strut-and-tie models in reinforced concrete structures analysed by a numerical technique

Os modelos de bielas e tirantes são procedimentos de análise apropriados para projetar elementos de concreto armado em casos de regiões onde há alterações geométricas ou concentrações de tensões, denominadas regiões D. Trata-se de bons modelos de representação da estrutura para avaliar melhor o seu comportamento estrutural e seu mecanismo resistente. O presente artigo aplica a técnica da otimização topológica para identificar o fluxo de tensões nas estruturas, definindo a configuração dos membros de bielas e tirantes, e quantifica seus valores para dimensionamento. Utilizam-se o método ESO, e uma variante desse, o SESO (Smoothing ESO) com o método dos elementos finitos em elasticidade plana. A filosofia do SESO baseia-se na observação de que se o elemento não for necessário à estrutura, sua contribuição de rigidez vai diminuindo progressivamente. Isto é, sua remoção é atenuada nos valores da matriz constitutiva, como se este estivesse em processo de danificação. Para validar a presente formulação, apresentam-se alguns exemplos numéricos onde se comparam suas respostas com as advindas de trabalhos científicos pioneiros sobre o assunto.

A common Tabu search algorithm for the global optimization of engineering problems

A novel common Tabu algorithm for global optimizations of engineering problems is presented. The robustness and efficiency of the presented method are evaluated by using standard mathematical functions and hy solving a practical engineering problem. The numerical results show that the proposed method is (i) superior to the conventional Tabu search algorithm in robustness, and (ii) superior to the simulated annealing algorithm in efficiency. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Lobatto deferred correction for stiff two-point boundary value problems

An iterated deferred correction algorithm based on Lobatto Runge-Kutta formulae is developed for the efficient numerical solution of nonlinear stiff two-point boundary value problems. An analysis of the stability properties of general deferred correction schemes which are based on implicit Runge-Kutta methods is given and results which are analogous to those obtained for initial value problems are derived. A revised definition of symmetry is presented and this ensures that each deferred correction produces an optimal increase in order. Finally, some numerical results are given to demonstrate the superior performance of Lobatto formulae compared with mono-implicit formulae on stiff two-point boundary value problems. (C) 1998 Elsevier B.V. Ltd. All rights reserved.

Bose-Einstein condensation dynamics from the numerical solution of the Gross-Pitaevskii equation

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii (GP) equation with both spherical and axial symmetries. We consider time-evolution problems initiated by suddenly changing the interatomic scattering length or harmonic trapping potential in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent GP equation, where all nonlinear and linear non-derivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.

Wavelet-Galerkin method for one-dimensional elastoplasticity and damage problems: Constitutive modeling and computational aspects

This work presents an analysis of the wavelet-Galerkin method for one-dimensional elastoplastic-damage problems. Time-stepping algorithm for non-linear dynamics is presented. Numerical treatment of the constitutive models is developed by the use of return-mapping algorithm. For spacial discretization we can use wavelet-Galerkin method instead of standard finite element method. This approach allows to locate singularities. The discrete formulation developed can be applied to the simulation of one-dimensional problems for elastic-plastic-damage models. (C) 2007 Elsevier B.V. All rights reserved.

Iterated deferred correction for linear two-point boundary value problems

An analysis of iterated deferred correction based on various classes of implicit Runge-Kutta formulae is given. Out of different possibilities considered, it is shown that those based purely on Lobatto formulae have the best stability. The enhanced stability of Lobatto schemes is very important for the efficient integration of excessively stiff boundary value problems and this is demonstrated by means of some numerical results.

ACCELERATING ITERATIVE ALGORITHMS FOR SYMMETRICAL LINEAR COMPLEMENTARITY-PROBLEMS

We propose a method for accelerating iterative algorithms for solving symmetric linear complementarity problems. The method consists in performing a one-dimensional optimization in the direction generated by a splitting method even for non-descent directions. We give strong convergence proofs and present numerical experiments that justify using this acceleration.