Performance analysis of FIR digital filter design: RNS versus traditional

Residue Number System (RNS) based Finite Impulse Response (FIR) digital filters and traditional FIR filters. This research is motivated by the importance of an efficient filter implementation for digital signal processing. The comparison is done in terms of speed and area requirement for various filter specifications. RNS based FIR filters operate more than three times faster and consumes only about 60% of the area than traditional filter when number of filter taps is more than 32. The area for RNS filter is increasing at a lesser rate than that for traditional resulting in lower power consumption. RNS is a nonweighted number system without carry propogation between different residue digits.This enables simultaneous parallel processing on all the digits resulting in high speed addition and multiplication in the RNS domain

Modelización constitutiva y computacional del daño y la fractura de materiales compuestos

En el trabajo se definen modelos constitutivos que permiten reproducir el proceso de fallo de estructuras de materiales compuestos en distintas escalas bajo cargas estáticas. Se define un modelo constitutivo para determinar la respuesta de estructuras de materiales compuestos mediante la teoría de laminados. El modelo es validado mediante un programa de ensayos experimentales con probetas con un agujero central geométricamente similares. Se muestra la capacidad del modelo de detectar el efecto tamaño. Se define un modelo constitutivo para materiales transversalmente isótropos bajo estados tridimensionales de tensión. El modelo se valida analizando numéricamente el proceso de agrietamiento de la matriz. Finalmente se desarrolla un modelo analítico para determinar el agrietamiento de la matriz y la delaminación entre las capas.

A cell by cell anisotropic adaptive mesh ALE scheme for the numerical solution of the Euler equations

In this paper a cell by cell anisotropic adaptive mesh technique is added to an existing staggered mesh Lagrange plus remap finite element ALE code for the solution of the Euler equations. The quadrilateral finite elements may be subdivided isotropically or anisotropically and a hierarchical data structure is employed. An efficient computational method is proposed, which only solves on the finest level of resolution that exists for each part of the domain with disjoint or hanging nodes being used at resolution transitions. The Lagrangian, equipotential mesh relaxation and advection (solution remapping) steps are generalised so that they may be applied on the dynamic mesh. It is shown that for a radial Sod problem and a two-dimensional Riemann problem the anisotropic adaptive mesh method runs over eight times faster.

A combined BIE-FE method for the Stokes equations

A numerical algorithm for the biharmonic equation in domains with piecewise smooth boundaries is presented. It is intended for problems describing the Stokes flow in the situations where one has corners or cusps formed by parts of the domain boundary and, due to the nature of the boundary conditions on these parts of the boundary, these regions have a global effect on the shape of the whole domain and hence have to be resolved with sufficient accuracy. The algorithm combines the boundary integral equation method for the main part of the flow domain and the finite-element method which is used to resolve the corner/cusp regions. Two parts of the solution are matched along a numerical ‘internal interface’ or, as a variant, two interfaces, and they are determined simultaneously by inverting a combined matrix in the course of iterations. The algorithm is illustrated by considering the flow configuration of ‘curtain coating’, a flow where a sheet of liquid impinges onto a moving solid substrate, which is particularly sensitive to what happens in the corner region formed, physically, by the free surface and the solid boundary. The ‘moving contact line problem’ is addressed in the framework of an earlier developed interface formation model which treats the dynamic contact angle as part of the solution, as opposed to it being a prescribed function of the contact line speed, as in the so-called ‘slip models’. Keywords: Dynamic contact angle; finite elements; free surface flows; hybrid numerical technique; Stokes equations.

Improving flood models using LiDAR: progress and problems

Airborne scanning laser altimetry (LiDAR) is an important new data source for river flood modelling. LiDAR can give dense and accurate DTMs of floodplains for use as model bathymetry. Spatial resolutions of 0.5m or less are possible, with a height accuracy of 0.15m. LiDAR gives a Digital Surface Model (DSM), so vegetation removal software (e.g. TERRASCAN) must be used to obtain a DTM. An example used to illustrate the current state of the art will be the LiDAR data provided by the EA, which has been processed by their in-house software to convert the raw data to a ground DTM and separate vegetation height map. Their method distinguishes trees from buildings on the basis of object size. EA data products include the DTM with or without buildings removed, a vegetation height map, a DTM with bridges removed, etc. Most vegetation removal software ignores short vegetation less than say 1m high. We have attempted to extend vegetation height measurement to short vegetation using local height texture. Typically most of a floodplain may be covered in such vegetation. The idea is to assign friction coefficients depending on local vegetation height, so that friction is spatially varying. This obviates the need to calibrate a global floodplain friction coefficient. It’s not clear at present if the method is useful, but it’s worth testing further. The LiDAR DTM is usually determined by looking for local minima in the raw data, then interpolating between these to form a space-filling height surface. This is a low pass filtering operation, in which objects of high spatial frequency such as buildings, river embankments and walls may be incorrectly classed as vegetation. The problem is particularly acute in urban areas. A solution may be to apply pattern recognition techniques to LiDAR height data fused with other data types such as LiDAR intensity or multispectral CASI data. We are attempting to use digital map data (Mastermap structured topography data) to help to distinguish buildings from trees, and roads from areas of short vegetation. The problems involved in doing this will be discussed. A related problem of how best to merge historic river cross-section data with a LiDAR DTM will also be considered. LiDAR data may also be used to help generate a finite element mesh. In rural area we have decomposed a floodplain mesh according to taller vegetation features such as hedges and trees, so that e.g. hedge elements can be assigned higher friction coefficients than those in adjacent fields. We are attempting to extend this approach to urban area, so that the mesh is decomposed in the vicinity of buildings, roads, etc as well as trees and hedges. A dominant points algorithm is used to identify points of high curvature on a building or road, which act as initial nodes in the meshing process. A difficulty is that the resulting mesh may contain a very large number of nodes. However, the mesh generated may be useful to allow a high resolution FE model to act as a benchmark for a more practical lower resolution model. A further problem discussed will be how best to exploit data redundancy due to the high resolution of the LiDAR compared to that of a typical flood model. Problems occur if features have dimensions smaller than the model cell size e.g. for a 5m-wide embankment within a raster grid model with 15m cell size, the maximum height of the embankment locally could be assigned to each cell covering the embankment. But how could a 5m-wide ditch be represented? Again, this redundancy has been exploited to improve wetting/drying algorithms using the sub-grid-scale LiDAR heights within finite elements at the waterline.

Second order accurate upwind difference schemes for scalar conservation laws with source terms

A second order accurate, characteristic-based, finite difference scheme is developed for scalar conservation laws with source terms. The scheme is an extension of well-known second order scalar schemes for homogeneous conservation laws. Such schemes have proved immensely powerful when applied to homogeneous systems of conservation laws using flux-difference splitting. Many application areas, however, involve inhomogeneous systems of conservation laws with source terms, and the scheme presented here is applied to such systems in a subsequent paper.

Classification of instability modes in a model of aluminium reduction cells with a uniform magnetic field

A unified view on the interfacial instability in a model of aluminium reduction cells in the presence of a uniform, vertical, background magnetic field is presented. The classification of instability modes is based on the asymptotic theory for high values of parameter β, which characterises the ratio of the Lorentz force based on the disturbance current, and gravity. It is shown that the spectrum of the travelling waves consists of two parts independent of the horizontal cross-section of the cell: highly unstable wall modes and stable or weakly unstable centre, or Sele’s modes. The wall modes with the disturbance of the interface being localised at the sidewalls of the cell dominate the dynamics of instability. Sele’s modes are characterised by a distributed disturbance over the whole horizontal extent of the cell. As β increases these modes are stabilized by the field.

Plane wave approximation in linear elasticity

We consider the approximation of solutions of the time-harmonic linear elastic wave equation by linear combinations of plane waves. We prove algebraic orders of convergence both with respect to the dimension of the approximating space and to the diameter of the domain. The error is measured in Sobolev norms and the constants in the estimates explicitly depend on the problem wavenumber. The obtained estimates can be used in the h- and p-convergence analysis of wave-based finite element schemes.

Sparse space-time boundary element methods for the heat equation

The goal of this work is the efficient solution of the heat equation with Dirichlet or Neumann boundary conditions using the Boundary Elements Method (BEM). Efficiently solving the heat equation is useful, as it is a simple model problem for other types of parabolic problems. In complicated spatial domains as often found in engineering, BEM can be beneficial since only the boundary of the domain has to be discretised. This makes BEM easier than domain methods such as finite elements and finite differences, conventionally combined with time-stepping schemes to solve this problem. The contribution of this work is to further decrease the complexity of solving the heat equation, leading both to speed gains (in CPU time) as well as requiring smaller amounts of memory to solve the same problem. To do this we will combine the complexity gains of boundary reduction by integral equation formulations with a discretisation using wavelet bases. This reduces the total work to O(h

Mapping the cytoskeletal prestress

Park CY, Tambe D, Alencar AM, Trepat X, Zhou EH, Millet E, Butler JP, Fredberg JJ. Mapping the cytoskeletal prestress. Am J Physiol Cell Physiol 298: C1245-C1252, 2010. First published February 17, 2010; doi: 10.1152/ajpcell.00417.2009.-Cell mechanical properties on a whole cell basis have been widely studied, whereas local intracellular variations have been less well characterized and are poorly understood. To fill this gap, here we provide detailed intracellular maps of regional cytoskeleton (CSK) stiffness, loss tangent, and rate of structural rearrangements, as well as their relationships to the underlying regional F-actin density and the local cytoskeletal prestress. In the human airway smooth muscle cell, we used micropatterning to minimize geometric variation. We measured the local cell stiffness and loss tangent with optical magnetic twisting cytometry and the local rate of CSK remodeling with spontaneous displacements of a CSK-bound bead. We also measured traction distributions with traction microscopy and cell geometry with atomic force microscopy. On the basis of these experimental observations, we used finite element methods to map for the first time the regional distribution of intracellular prestress. Compared with the cell center or edges, cell corners were systematically stiffer and more fluidlike and supported higher traction forces, and at the same time had slower remodeling dynamics. Local remodeling dynamics had a close inverse relationship with local cell stiffness. The principal finding, however, is that systematic regional variations of CSK stiffness correlated only poorly with regional F-actin density but strongly and linearly with the regional prestress. Taken together, these findings in the intact cell comprise the most comprehensive characterization to date of regional variations of cytoskeletal mechanical properties and their determinants.

Um modelo analítico e computacional para consideração de efeitos de envelhecimento em estruturas de material compósito

Os materiais compósitos reforçados por fibras apresentam vantagens quando comparados aos materiais de construção mais tradicionais como concreto e aço. Por outro lado, devido ao fato destes materiais serem relativamente recentes no mercado, questões a respeito de sua durabilidade são ainda objeto de discussão e faz-se necessária intensa pesquisa sobre o envelhecimento dos compósitos. Como conseqüência, recentemente têm surgido inúmeros trabalhos à respeito da degradação dos compósitos considerando efeitos como temperatura, oxidação, radiação UV, condições de carregamento, etc. A maioria destas pesquisas, no entanto, são realizadas a nível de material e não são diretamente aplicáveis à situações de projeto. Desta forma, existe grande demanda por novos estudos e dados compatíveis com aplicações estruturais. Neste trabalho apresenta-se um modelo analítico-numérico adequado para, interpretação e aplicação destes dados experimentais em análise e projeto de estrutural. A formulação proposta inclui relações constitutivas elásticas anisotrópicas com envelhecimento, relações constitutivas viscoelásticas anisotrópicas com envelhecimento em termos de variáveis de estado, análise de falhas com critério de degradação ajustado à idade do material e considera-se grandes deslocamentos e pequenas deformações. As diferenças essenciais entre os processos de envelhecimento em endurecimento e amolecimento são descritos juntamente com as relações constitutivas para cada caso. Estas equações são deduzidas na forma adequada para análise numérica via método dos elementos finitos usando uma solução incremental-iterativa com consideração de efeitos pos-críticos. Vários exemplos são apresentados, incluindo análises elásticas, viscoelásticas e de falha com envelhecimento.

Ensino de Entropia: um enfoque histórico e epistemológico.

The rational construction necessary to systematize scientific knowledge in physics, introduces difficulties of understanding in some of its concepts. One of these concepts which exemplify properly this difficulty in learning or teaching is entropy. This thesis propose the construction of a didactic route which constitute itself a historical and epistemological course to entropy, intending to contribute for teaching this concept as well as other physics concepts. The basic assumption to build this route is that through the historical review of the development of this concept in the way suggested by Bachelard s (1884-1962) epistemology it is possible to make subjects, to be taught and learned, more meaningful. Initially I composed a brief biographical note to give the reader an idea about the issues, interests and reflections, related to science, and how I dealt with them in my private and professional life, as well as the role they played to lead me to write this thesis. The strategy to construct the route to entropy was to split the usual contents of basic thermodynamics in three moments in a way they can constitute epistemological units , which can be identified by the way of thinking in the corresponding moments of scientific knowledge production: a technical and empiricist moment, a rationalist and positivist moment and a post-positivist rationalist one. The transition between each moment is characterized by a rupture with the former way of thinking; however the progress in the construction of knowledge in the area is evident. As the final part of this work I present an analysis based on elements of Bachelard s epistemology that are present in each moment. This analysis is the basic component of the didactic route that I propose myself to build. The way I made this route guide to entropy could contribute to the construction of other didactic routes in physics and other sciences, in a way to unveil hidden meanings and as a tool to humanize scientific knowledge.

Modelos para análise e dimensionamento de painéis de alvenaria estrutural

The main objective of this thesis was the study of bracing panels of structural masonry, by applying the Finite Element Method and Strut and Tie Method. It was analyzed the following aspects: the effect of orthotropy on the behavior of the panels; distribution of horizontal forces between panels for buildings; comparison between Equivalent Frame and Finite Elements models; panels design with the Strut and Tie Method. The results showed that one should not disregard the orthotropy, otherwise this can lead to models stiffer than the real. Regarding the distribution of horizontal forces, showed that the disregard of lintels and shear deformation leads to significant differences in the simplified model. The results showed also that the models in Finite Element and Equivalent Frame exhibit similar behavior in respect to stiffness of panels and stress distribution over the sessions requested. It was discussing criteria for designing Strut and Tie Method models in one floor panels. Then, the theoretical strength these panels was compared with the rupture strength of panels tested in the literature. The theoretical maximum strength were always less than the rupture strength of the panels obtained in tests, due to the fact that the proposed model cannot represent the behavior of the masonry after the start of the panel cracking due to plasticization of the reinforcement

Análise não-linear de treliças pelo método dos elementos finitos posicional

This work presents the positional nonlinear geometric formulation for trusses using different strain measures. The positional formulation presents an alternative approach for nonlinear problems. This formulation considers nodal positions as variables of the nonlinear system instead of displacements (widely found in literature). The work also describes the arc-length method used for tracing equilibrium paths with snap-through and snap-back. Numerical applications for trusses already established in the literature and comparisons with other studies are provided to prove the accuracy of the proposed formulation