Erosion of an active fault scarp leads to drainage capture in the Amazon region, Brazil

Far from the continental margin, drainage basins in Central Amazonia should be in topographic steady state; but they are not. Abandoned remnant fluvial valleys up to hundreds of square kilometers in size are observed throughout Amazonia, and are evidence of significant landscape reorganization. While major Late Miocene drainage shifts occurred due to initiation of the transcontinental Amazon River, local landscape change has remained active until today. Driven either by dynamic topography, tectonism, and/or climatic fluctuations, drainage captures in Amazonia provide a natural experiment for assessing the geomorphic response of low-slope basins to sudden, capture related base-level falls. This paper evaluates the timing of geomorphic change by examining a drainage capture event across the Baependi fault scarp involving the Cuieiras and TarumA-Mirim River basins northwest of the city of Manaus in Brazil. A system of capture-related knickpoints was generated by base-level fall following drainage capture; through numerical modeling of their initiation and propagation, the capture event is inferred to have occurred between the middle and late Pleistocene, consistent with other studies of landscape change in surrounding areas. In low-slope settings like the Amazon River basin, base-level fall can increase erosion rates by more than an order of magnitude, and moderate to large river basins can respond to episodes of base-level fall over timescales of tens to hundreds of thousands of years. Copyright (c) 2013 John Wiley & Sons, Ltd.

Control of Limit Cycle Oscillation in a Three Degrees of Freedom Airfoil Section Using Fuzzy Takagi-Sugeno Modeling

This work presents a strategy to control nonlinear responses of aeroelastic systems with control surface freeplay. The proposed methodology is developed for the three degrees of freedom typical section airfoil considering aerodynamic forces from Theodorsen's theory. The mathematical model is written in the state space representation using rational function approximation to write the aerodynamic forces in time domain. The control system is designed using the fuzzy Takagi-Sugeno modeling to compute a feedback control gain. It useds Lyapunov's stability function and linear matrix inequalities (LMIs) to solve a convex optimization problem. Time simulations with different initial conditions are performed using a modified Runge-Kutta algorithm to compare the system with and without control forces. It is shown that this approach can compute linear control gain able to stabilize aeroelastic systems with discontinuous nonlinearities.

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.

Numerical and experimental investigations of three-dimensional container filling with Newtonian viscous fluids

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

Numerical solution of the FENE-CR model in complex flows

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

A Numerical Study of the FENE-CR Model Applied to a Jet Flow Problem

The FENE-CR model is investigated through a numerical algorithm to simulate the time-dependent moving free surface flow produced by a jet impinging on a flat surface. The objective is to demonstrate that by increasing the extensibility parameter L, the numerical solutions converge to the solutions obtained with the Oldroyd-B model. The governing equations are solved by an established free surface flow solver based on the finite difference and marker-and-cell methods. Numerical predictions of the extensional viscosity obtained with several values of the parameter L are presented. The results show that if the extensibility parameter L is sufficiently large then the extensional viscosities obtained with the FENE-CR model approximate the corresponding Oldroyd-B viscosity. Moreover, the flow from a jet impinging on a flat surface is simulated with various values of the extensibility parameter L and the fluid flow visualizations display convergence to the Oldroyd-B jet flow results.

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