An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 1: Rods

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

Stochastic gauges in quantum dynamics for many-body simulations

Quantum dynamics simulations can be improved using novel quasiprobability distributions based on non-orthogonal Hermitian kernel operators. This introduces arbitrary functions (gauges) into the stochastic equations. which can be used to tailor them for improved calculations. A possible application to full quantum dynamic simulations of BEC's is presented. (C) 2001 Elsevier Science B.V. All rights reserved.

Quantum interference in optical fields and atomic radiation

We discuss the connection between quantum interference effects in optical beams and radiation fields emitted from atomic systems. We illustrate this connection by a study of the first- and second-order correlation functions of optical fields and atomic dipole moments. We explore the role of correlations between the emitting systems and present examples of practical methods to implement two systems with non-orthogonal dipole moments. We also derive general conditions for quantum interference in a two-atom system and for a control of spontaneous emission. The relation between population trapping and dark states is also discussed. Moreover, we present quantum dressed-atom models of cancellation of spontaneous emission, amplification on dark transitions, fluorescence quenching and coherent population trapping.

Local Binary Patterns in Focal-Plane Processing. Analysis and Applications

Feature extraction is the part of pattern recognition, where the sensor data is transformed into a more suitable form for the machine to interpret. The purpose of this step is also to reduce the amount of information passed to the next stages of the system, and to preserve the essential information in the view of discriminating the data into different classes. For instance, in the case of image analysis the actual image intensities are vulnerable to various environmental effects, such as lighting changes and the feature extraction can be used as means for detecting features, which are invariant to certain types of illumination changes. Finally, classification tries to make decisions based on the previously transformed data. The main focus of this thesis is on developing new methods for the embedded feature extraction based on local non-parametric image descriptors. Also, feature analysis is carried out for the selected image features. Low-level Local Binary Pattern (LBP) based features are in a main role in the analysis. In the embedded domain, the pattern recognition system must usually meet strict performance constraints, such as high speed, compact size and low power consumption. The characteristics of the final system can be seen as a trade-off between these metrics, which is largely affected by the decisions made during the implementation phase. The implementation alternatives of the LBP based feature extraction are explored in the embedded domain in the context of focal-plane vision processors. In particular, the thesis demonstrates the LBP extraction with MIPA4k massively parallel focal-plane processor IC. Also higher level processing is incorporated to this framework, by means of a framework for implementing a single chip face recognition system. Furthermore, a new method for determining optical flow based on LBPs, designed in particular to the embedded domain is presented. Inspired by some of the principles observed through the feature analysis of the Local Binary Patterns, an extension to the well known non-parametric rank transform is proposed, and its performance is evaluated in face recognition experiments with a standard dataset. Finally, an a priori model where the LBPs are seen as combinations of n-tuples is also presented

Theory of laser-induced ultrafast structural changes in solids

The present thesis is a contribution to the study of laser-solid interaction. Despite the numerous applications resulting from the recent use of laser technology, there is still a lack of satisfactory answers to theoretical questions regarding the mechanism leading to the structural changes induced by femtosecond lasers in materials. We provide here theoretical approaches for the description of the structural response of different solids (cerium, samarium sulfide, bismuth and germanium) to femtosecond laser excitation. Particular interest is given to the description of the effects of the laser pulse on the electronic systems and changes of the potential energy surface for the ions. Although the general approach of laser-excited solids remains the same, the potential energy surface which drives the structural changes is calculated with different theoretical models for each material. This is due to the difference of the electronic properties of the studied systems. We use the Falicov model combined with an hydrodynamic method to study photoinduced phase changes in cerium. The local density approximation (LDA) together with the Hubbard-type Hamiltonian (LDA+U) in the framework of density functional theory (DFT) is used to describe the structural properties of samarium sulfide. We parametrize the time-dependent potential energy surface (calculated using DFT+ LDA) of bismuth on which we perform quantum dynamical simulations to study the experimentally observed amplitude collapse and revival of coherent $A_{1g}$ phonons. On the basis of a time-dependent potential energy surface calculated from a non-orthogonal tight binding Hamiltonian, we perform molecular dynamics simulation to analyze the time evolution (coherent phonons, ultrafast nonthermal melting) of germanium under laser excitation. The thermodynamic equilibrium properties of germanium are also reported. With the obtained results we are able to give many clarifications and interpretations of experimental results and also make predictions.

Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tesselations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique.

Ritmo Circadiano de Triiodotironina (T3) e Tiroxina (T4) Plasmática em Ovelhas Ideal durante o Anestro Estacional

The objective was to determine the plasma levels of the triiodothyronine (T-3) and thyroxine (T-4), in six Ideal ewe sheeps during seasonal anoestrus, by radioimmunoassay. The blood samples were collected from the jugular vein catheter over a period of 24 hours, at 2-h intervals (12, 14, 16, IX, 20, 22, 24, OZ, 04, 06, OX, and 10 hours). II was used analysis of variance for data adjustments, and the decomposition of the collected periods in non-orthogonal polynomial. It was observed a crescent linear effect of the T-4 plasma levels, presenting lower level of 49.60+/-11.41 and higher level of 76.99+/-10.0 ng/ml, at fourth and last blood samples collection, respectively. T-3 plasma levels showed a significant cubic response during 24 hours, with the highest value of 0.94+/-0.20 ng/ml at the end of the the study. In conclusion, thyroidal hormones levels were lower in the first blood samples collection, increasing progressively during the night, reaching the highest levels at the morning period.

Variance dispersion graphs for comparing blocked response surface designs

In this paper, we extend the use of the variance dispersion graph (VDG) to experiments in which the response surface (RS) design must be blocked. Through several examples we evaluate the prediction performances of RS designs in non-orthogonal block designs compared with the equivalent unblocked designs and orthogonally blocked designs. These examples illustrate that good prediction performance of designs in small blocks can be expected in practice. Most importantly, we show that the allocation of the treatment set to blocks can seriously affect the prediction properties of designs; thus, much care is needed in performing this allocation.

Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Reconhecimento de padrões via Wavelets: aplicação no monitoramento de estruturas a partir do efeito do multicaminho proveniente de sinais de satélites GPS

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

Desenvolvimento de Interface gráfica como suporte para soluções numéricas das equações de Maxwell em coordenadas gerais – 3D

Neste trabalho, é implementada uma interface gráfica de usuários (GUI) usando a ferramenta Qt da Nokia (versão 3.0). A interface visa simplificar a criação de cenários para a realização de simulações paralelas usando a técnica numérica Local Nonorthogonal Finite Difference Time-Domain (LN-FDTD), aplicada para solucionar as equações de Maxwell. O simulador foi desenvolvido usando a linguagem de programação C e paralelizado utilizando threads. Para isto, a biblioteca pthread foi empregada. A visualização 3D do cenário a ser simulado (e da malha) é realizada por um programa especialmente desenvolvido que utiliza a biblioteca OpenGL. Para melhorar o desenvolvimento e alcançar os objetivos do projeto computacional, foram utilizados conceitos da Engenharia de Software, tais como o modelo de processo de software por prototipagem. Ao privar o usuário de interagir diretamente com o código-fonte da simulação, a probabilidade de ocorrência de erros humanos durante o processo de construção de cenários é minimizada. Para demonstrar o funcionamento da ferramenta desenvolvida, foi realizado um estudo relativo ao efeito de flechas em linhas de baixa tensão nas tensões transitórias induzidas nas mesmas por descargas atmosféricas. As tensões induzidas nas tomadas da edificação também são estudadas.

Uso de coproduto da produção de biodiesel como alimento para bovinos confinados

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

Desempenho de novilhos Nelore confinados alimentados com óleo de linhaça in natura ou protegido

Pós-graduação em Zootecnia - FCAV

The General Curvilinear Environmental Model (GCEM)

[EN] The General Curvilinear Environmental Model is a high-resolution system composed of the General Curvilinear Coastal Ocean Model (GCCOM) and the General Curvilinear Atmospheric Model (GCAM). Both modules are capable of reading a general curvilinear grid, orthogonal as well as non-orthogonal in all three directions. These two modules are weakly coupled using the distributed coupling toolkit (DCT). The model can also be nested within larger models and users are able to interact with the model and run it using a web based computational environment.

Petrography and geochemical composition of the Atlantis II Fracture Zone

SeaBeam echo sounding, seismic reflection, magnetics, and gravity profiles were run along closely spaced tracks (5 km) parallel to the Atlantis II Fracture Zone on the Southwest Indian Ridge, giving 80% bathymetric coverage of a 30- * 170-nmi strip centered over the fracture zone. The southern and northern rift valleys of the ridge were clearly defined and offset north-south by 199 km. The rift valleys are typical of those found elsewhere on the Southwest Indian Ridge, with relief of more than 2200 m and widths from 22 to 38 km. The ridge-transform intersections are marked by deep nodal basins lying on the transform side of the neovolcanic zone that defines the present-day spreading axis. The walls of the transform generally are steep (25°-40°), although locally, they can be more subdued. The deepest point in the transform is 6480 m in the southern nodal basin, and the shallowest is an uplifted wave-cut terrace that exposes plutonic rocks from the deepest layer of the ocean crust at 700 m. The transform valley is bisected by a 1.5-km-high median tectonic ridge that extends from the northern ridge-transform intersection to the midpoint of the active transform. The seismic survey showed that the floor of the transform contains up to 0.5 km of sediment. Piston-coring at two locations on the transform floor recovered more than 1 m of sand and gravel, which appears to be turbidites shed from the walls of the fracture zone. Extensive dredging showed that more than two-thirds of the crust exposed in the transform valley and its walls were plutonic rocks, principally gabbros and residual mantle peridotites. In contrast, based on dredging and seafloor morphology, only relatively undisrupted pillow basalt flows have been exposed on crust of the same age spreading away from the transform. Magnetic anomalies are well defined out to 11 m.y. over the flanking transverse ridges and transform valley, even where layer 2 appears to be absent. The total opening rate is 1.6 cm/yr, but the arrangement of the anomalies indicates that the spreading for each ridge is asymmetric, with the ridge flanks facing the transform spreading at a rate of 1.0 cm/yr. Such an asymmetric spreading pattern requires that both the northern and southern ridges migrate away from each other at 0.2 cm/yr, thus lengthening the transform at 0.4 cm/yr for the last 11 m.y. To the north, the fracture zone valley is oriented differently from the present-day transform, indicating a paleospreading direction change at 17 m.y. from N10°E to due north-south. This change placed the transform into extension for the 11-m.y. period required for simple orthogonal ridge-transform geometry to be reestablished and produced a large transtensional basin within the transform valley. This basin was split by continued transform slip after 11 m.y., with the larger half moving to the north with the African Plate.