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,

Determination of the numerical parameters of a continuous damage model for the structural analysis of clay brick masonry

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

Finite-Element Analysis of Stress on Dental Implant Prosthesis

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

A planar flexible robotic manipulator

Presents the dynamic modelling of a flexible robotic manipulator with two flexible links and two revolute joints, which rotates in the horizontal plane. The dynamic equations are derived using the Newton-Euler formulation and the finite element method, based on elementary beam theory, which is used to discretize the displacements such that the small motion is represented in terms of nodal displacements. Computer simulation results are presented to illustrate this study. The dynamic model becomes necessary for use in future design and control applications.

Vibration Attenuation in Rotating Machines Using Smart Spring Mechanism

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

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

Influence of Cusp Inclination on Stress Distribution in Implant-Supported Prostheses. A Three-Dimensional Finite Element Analysis

Purpose: The aim of this study was to assess the influence of cusp inclination on stress distribution in implant-supported prostheses by 3D finite element method.Materials and Methods: Three-dimensional models were created to simulate a mandibular bone section with an implant (3.75 mm diameter x 10 mm length) and crown by means of a 3D scanner and 3D CAD software. A screw-retained single crown was simulated using three cusp inclinations (10 degrees, 20 degrees, 30 degrees). The 3D models (model 10d, model 20d, and model 30d) were transferred to the finite element program NeiNastran 9.0 to generate a mesh and perform the stress analysis. An oblique load of 200 N was applied on the internal vestibular face of the metal ceramic crown.Results: The results were visualized by means of von Mises stress maps. Maximum stress concentration was located at the point of application. The implant showed higher stress values in model 30d (160.68 MPa). Cortical bone showed higher stress values in model 10d (28.23 MPa).Conclusion: Stresses on the implant and implant/abutment interface increased with increasing cusp inclination, and stresses on the cortical bone decreased with increasing cusp inclination.

A marker-and-cell approach to free surface 2-D multiphase flows

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

MATHEMATICA notebook for computing tetrahedral finite element shape functions and matrices for the Helmholtz equation

A MATHEMATICA notebook to compute the elements of the matrices which arise in the solution of the Helmholtz equation by the finite element method (nodal approximation) for tetrahedral elements of any approximation order is presented. The results of the notebook enable a fast computational implementation of finite element codes for high order simplex 3D elements reducing the overheads due to implementation and test of the complex mathematical expressions obtained from the analytical integrations. These matrices can be used in a large number of applications related to physical phenomena described by the Poisson, Laplace and Schrodinger equations with anisotropic physical properties.

Osmotic dehydration of apples in sugar/salt solutions: Concentration profiles and effective diffusion coefficients

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

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

Short and long period optimization of drug doses in the treatment of AIDS

Técnicas de otimização numérica são úteis na solução de problemas de determinação da melhor entrada para sistemas descritos por modelos matemáticos e cujos objetivos podem ser expressos de uma maneira quantitativa. Este trabalho aborda o problema de otimizar as dosagens dos medicamentos no tratamento da AIDS em termos de um balanço entre a resposta terapêutica e os efeitos colaterais. Um modelo matemático para descrever a dinâmica do vírus HIV e células CD4 é utilizado para calcular a dosagem ótima do medicamento no tratamento a curto prazo de pacientes com AIDS por um método de otimização direta utilizando uma função custo do tipo Bolza. Os parâmetros do modelo foram ajustados com dados reais obtidos da literatura. Com o objetivo de simplificar os procedimentos numéricos, a lei de controle foi expressa em termos de uma expansão em séries que, após truncamento, permite obter controles sub-ótimos. Quando os pacientes atingem um estado clínico satisfatório, a técnica do Regulador Linear Quadrático (RLQ) é utilizada para determinar a dosagem permanente de longo período para os medicamentos. As dosagens calculadas utilizando a técnica RLQ , tendem a ser menores do que a equivalente terapia de dose constante em termos do expressivo aumento na contagem das células T+ CD4 e da redução da densidade de vírus livre durante um intervalo fixo de tempo.

Interpolating EFGM for computing continuous and discontinuous electromagnetic fields

Purpose - This paper proposes an interpolating approach of the element-free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value problems in domains involving material non-homogeneities. The suitability and efficiency of the proposed implementation are evaluated for a given set of test cases of electrostatic field in domains involving different material interfaces.Design/methodology/approach - the authors combined an interpolating approximation with a modified domain truncation scheme, which avoids additional techniques for enforcing the Dirichlet boundary conditions and for dealing with material interfaces usually employed in meshfree formulations.Findings - the local electric potential and field distributions were correctly described as well as the global quantities like the total potency and resistance. Since, the treatment of the material interfaces becomes practically the same for both the finite element method (FEM) and the proposed EFGM, FEM-oriented programs can, thus, be easily extended to provide EFGM approximations.Research limitations/implications - the robustness of the proposed formulation became evident from the error analyses of the local and global variables, including in the case of high-material discontinuity.Practical implications - the proposed approach has shown to be as robust as linear FEM. Thus, it becomes an attractive alternative, also because it avoids the use of additional techniques to deal with boundary/interface conditions commonly employed in meshfree formulations.Originality/value - This paper reintroduces the domain truncation in the EFGM context, but by using a set of interpolating shape functions the authors avoided the use of Lagrange multipliers as well Mathematics in Engineering high-material discontinuity.

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.

Concentration profiles and effective diffusion coefficients of sucrose and water in osmo-dehydrated apples

The spatial distribution of water and sugars in half-fresh apples dehydrated in sucrose solutions (30% and 50% w/w, 27 degrees C) for 2, 4 and 8 h, was determined. Each half was sliced as from the exposed surface. The density, water and sugar contents were determined for each piece. A mathematical model was fitted to the experimental data of the water and sucrose contents considering the overall flux and tissue shrinkage. A numerical method of finite differences permitted the calculation of the effective diffusion coefficients as a function of concentration, using material coordinates and integrating the two differential equations (for water and sucrose) simultaneously. The coefficients obtained were one or even two orders of magnitude lower than those for pure solutions and presented unusual concentration dependence. The behaviour of the apple tissue was also studied using light microscopy techniques to obtain images of the osmotically treated pieces (20%, 30% and 50% w/w sucrose solutions for 2, 4 and 8 h). (c) 2006 Elsevier Ltd. All rights reserved.