Application of the PTT Model to Axisymmetric Free Surface Flows

This work is concerned with numerical simulation of axisymmetric viscoelastic free surface flows using the Phan-Thien-Tanner (PTT) constitutive equation. A finite difference technique for solving the governing equations for unsteady incompressible flows written in Cylindrical coordinates on a staggered grid is described. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are applied. The numerical method is verified by comparing numerical predictions of fully developed flow in a pipe with the corresponding analytic solutions. To demonstrate that the numerical method can simulate axisymmetric free surface flows governed by the PTT model, numerical results of the flow evolution of a drop impacting on a rigid dry plate are presented. In these simulations, the rheological effects of the parameters epsilon and xi are investigated.

Control of instabilities in non-Newtonian free surface fluid flows

Free surface flows in inclined channels can develop periodic instabilities that are propagated downstream as shock waves with well-defined wavelengths and amplitudes. Such disturbances are called roll waves and are common in channels, torrential lava, landslides, and avalanches. The prediction and detection of such waves over certain types of structures and environments are useful for the prevention of natural risks. In this work, a mathematical model is established using a theoretical approach based on Cauchy's equations with the Herschel-Bulkley rheological model inserted into the viscous part of the stress tensor. This arrangement can adequately represent the behavior of muddy fluids, such as water-clay mixture. Then, taking into account the shallow water and the Rankine-Hugoniot's (shock wave) conditions, the equation of the roll wave and its properties, profile, and propagation velocity are determined. A linear stability analysis is performed with an emphasis on determining the condition that allows the generation of such instabilities, which depends on the minimum Froude number. A sensitivity analysis on the numerical parameters is performed, and numerical results including the influence of the Froude number, the index flow and dimensionless yield stress on the amplitude, the wavelength of roll waves and the propagation velocity of roll waves are shown. We show that our numerical results were in agreement with Coussot's experimental results (1994).

Modeling of piezoelectric energy harvesting considering the dependence of the rectifier circuit

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

Modeling 3-D desiccation soil crack networks using a mesh fragmentation technique

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

On the contact stiffness and nonlinear vibration of an elastic body with a rough surface in contact with a rigid flat surface

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

Inverse Modeling of GPS Multipath for Snow Depth Estimation-Part I: Formulation and Simulations

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

Modelling the surface circulation and thermal structure of a tropical reservoir using three-dimensional hydrodynamic lake model and remote-sensing data

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

Modeling of multiple cracks in reinforced concrete members using solid finite elements with high aspect ratio

This paper presents a new technique to model interfaces by means of degenerated solid finite elements, i.e., elements with a very high aspect ratio, with the smallest dimension corresponding to the thickness of the interfaces. It is shown that, as the aspect ratio increases, the element strains also increase, approaching the kinematics of the strong discontinuity. A tensile damage constitutive relation between strains and stresses is proposed to describe the nonlinear behavior of the interfaces associated with crack opening. To represent crack propagation, couples of triangular interface elements are introduced in between all regular (bulk) elements of the original mesh. With this technique the analyses can be performed integrally in the context of the continuum mechanics and complex crack patterns involving multiple cracks can be simulated without the need of tracking algorithms. Numerical tests are performed to show the applicability of the proposed technique, studding also aspects related to mesh objectivity.

Dinâmica de rotação de exoplanetas tipo terrestre com perturbação de terceiro corpo

In this work we study some topics of Celestial Mechanics, namely the problem of rigid body rotation and “spin-orbit” resonances. Emphasis is placed on the problem formulation and applications to some exoplanets with physical parameters (e.g. mass and radius) compatible with a terrestrial type constitution (e.g. rock) belonging to multiple planetary systems. The approach is both analytical and numerical. The analytical part consists of: i) the deduction of the equation of motion for the rotation problem of a spherical body with no symmetry, disturbed by a central body; ii) modeling the same problem by including a third-body in the planet-star system; iii) formulation of the concept of “spin-orbit” resonance in which the orbital period of the planet is an integer multiple of the rotation’s period. Topics of dynamical systems (e.g. equilibrium points, chaos, surface sections, etc.) will be included at this stage. In the numerical part simulations are performed with numerical models developed in the previous analytical section. As a first step we consider the orbit of the planet not perturbed by a third-body in the star-planet system. In this case the eccentricity and orbital semi-major axis of the planet are constants. Here the technique of surface sections, widely used in dynamical systems are applied. Next, we consider the action of a third body, developing a more realistic model for planetary rotation. The results in both cases are compared. Since the technique of disturbed surface sections is no longer applicable, we quantitatively evaluate the evolution of the characteristic angles of rotation (e.g. physical libration) by studying the evolution of individual orbits in the dynamically important regions of phase space, the latter obtained in the undisturbed case

Modeling of As2S3 raman fiber lasers operating in the mid-infrared

A numerical study of As2S3 Raman fiber lasers is carried out to show their potential for the entire coverage of the 3–4-m spectral band. Experimental results are first obtained from such a laser operated under controlled conditions in order to set the fiber parameters (i.e., gain and attenuation coefficients) to be used in the numerical model. An exhaustive numerical analysis is then performed to establish the conditions for optimum lasing performances over the entire 3–4-m spectral band.

Surface coverage in wireless sensor networks based on delaunay tetrahedralization

In this work is presented a new method for sensor deployment on 3D surfaces. The method was structured on different steps. The first one aimed discretizes the relief of interest with Delaunay algorithm. The tetrahedra and relative values (spatial coordinates of each vertex and faces) were input to construction of 3D Voronoi diagram. Each circumcenter was calculated as a candidate position for a sensor node: the corresponding circular coverage area was calculated based on a radius r. The r value can be adjusted to simulate different kinds of sensors. The Dijkstra algorithm and a selection method were applied to eliminate candidate positions with overlapped coverage areas or beyond of surface of interest. Performance evaluations measures were defined using coverage area and communication as criteria. The results were relevant, once the mean coverage rate achieved on three different surfaces were among 91% and 100%.

Coupling remote sensing bio-optical and three-dimensional hydrodynamic modeling to study the phytoplankton dynamics in a tropical hydroelectric reservoir

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

Computer graphics applied to anatomy: a study of two bio-cad modeling methods on finite element analysis of human edentulous hemi-mandible

Modeling is a step to perform a finite element analysis. Different methods of model construction are reported in literature, as the Bio-CAD modeling. The purpose of this study was to perform a model evaluation and application using two methods of Bio-CAD modeling from human edentulous hemi-mandible on the finite element analysis. From CT scans of dried human skull was reconstructed a stereolithographic model. Two methods of modeling were performed: STL conversion approach (Model 1) associated to STL simplification and reverse engineering approach (Model 2). For finite element analysis was used the action of lateral pterygoid muscle as loading condition to assess total displacement (D), equivalent von-Mises stress (VM) and maximum principal stress (MP). Two models presented differences on the geometry regarding surface number (1834 (model 1); 282 (model 2)). Were observed differences in finite element mesh regarding element number (30428 nodes/16683 elements (model 1); 15801 nodes/8410 elements (model 2). D, VM and MP stress areas presented similar distribution in two models. The values were different regarding maximum and minimum values of D (ranging 0-0.511 mm (model 1) and 0-0.544 mm (model 2), VM stress (6.36E-04-11.4 MPa (model 1) and 2.15E-04-14.7 MPa (model 2) and MP stress (-1.43-9.14 MPa (model 1) and -1.2-11.6 MPa (model 2). From two methods of Bio-CAD modeling, the reverse engineering presented better anatomical representation compared to the STL conversion approach. The models presented differences in the finite element mesh, total displacement and stress distribution.

Influence of summertime mesoscale convective systems on the heat balance and surface mixed layer dynamics of a large Amazonian hydroelectric reservoir

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

Modeling the formation and propagation of desiccation cracks in soils

The problem of desiccation cracks in soils has received increasing attention in the last few years, in both experimental investigations and modeling. Experimental research has been mainly focused on the behavior of slurries subjected to drying in plates of different shapes, sizes and thickness. The main objectives of these studies were to learn about the process of crack formation under controlled environmental conditions, and also to better understand the effect of different factors (e.g. soil type, boundary conditions, soil thickness) on the morphology of the crack network. As for the numerical modeling, different approaches have been suggested lately to describe the behavior of drying cracks in soils. One aspect that it is still difficult to describe properly is the crack pattern observed in desiccated soils. This work presents a novel technique to model the behavior of drying soils. The crack patter observed in desiccation tests on circular plates are simulated with the main objective of predicting the effect of soil thickness on crack pattern. Good agreement between experimental results and model prediction are observed.