Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Robust tori in a double-waved Hamiltonian model

A Hamiltonian system perturbed by two waves with particular wave numbers can present robust tori, which are barriers created by the vanishing of the perturbed Hamiltonian at some defined positions. When robust tori exist, any trajectory in phase space passing close to them is blocked by emergent invariant curves that prevent the chaotic transport. Our results indicate that the considered particular solution for the two waves Hamiltonian model shows plenty of robust tori blocking radial transport. (C) 2010 Elsevier B.V. All rights reserved.

Plasma confinement in tokamaks with robust torus

We present a non-linear symplectic map that describes the alterations of the magnetic field lines inside the tokamak plasma due to the presence of a robust torus (RT) at the plasma edge. This RT prevents the magnetic field lines from reaching the tokamak wall and reduces, in its vicinity, the islands and invariant curve destruction due to resonant perturbations. The map describes the equilibrium magnetic field lines perturbed by resonances created by ergodic magnetic limiters (EMLs). We present the results obtained for twist and non-twist mappings derived for monotonic and non-monotonic plasma current density radial profiles, respectively. Our results indicate that the RT implementation would decrease the field line transport at the tokamak plasma edge. (C) 2010 Elsevier B.V. All rights reserved.

Utilização de modelos de programação não-linear para sistematização de terras para irrigação em áreas regulares

Objetivou-se com este trabalho, desenvolver modelos de programação não-linear para sistematização de terras, aplicáveis para áreas com formato regular e que minimizem a movimentação de terra, utilizando o software GAMS para o cálculo. Esses modelos foram comparados com o Método dos Quadrados Mínimos Generalizado, desenvolvido por Scaloppi & Willardson (1986), sendo o parâmetro de avaliação o volume de terra movimentado. Concluiu-se que, ambos os modelos de programação não-linear desenvolvidos nesta pesquisa mostraram-se adequados para aplicação em áreas regulares e forneceram menores valores de movimentação de terra quando comparados com o método dos quadrados mínimos.

Modelagem linear e identificação do modelo dinâmico de um robô móvel com acionamento diferencial

This work presents a modelling and identification method for a wheeled mobile robot, including the actuator dynamics. Instead of the classic modelling approach, where the robot position coordinates (x,y) are utilized as state variables (resulting in a non linear model), the proposed discrete model is based on the travelled distance increment Delta_l. Thus, the resulting model is linear and time invariant and it can be identified through classical methods such as Recursive Least Mean Squares. This approach has a problem: Delta_l can not be directly measured. In this paper, this problem is solved using an estimate of Delta_l based on a second order polynomial approximation. Experimental data were colected and the proposed method was used to identify the model of a real robot

Controlador preditivo não linear aplicado ao controle de golfadas em processos de produção de petróleo

Slugging is a well-known slugging phenomenon in multiphase flow, which may cause problems such as vibration in pipeline and high liquid level in the separator. It can be classified according to the place of its occurrence. The most severe, known as slugging in the riser, occurs in the vertical pipe which feeds the platform. Also known as severe slugging, it is capable of causing severe pressure fluctuations in the flow of the process, excessive vibration, flooding in separator tanks, limited production, nonscheduled stop of production, among other negative aspects that motivated the production of this work . A feasible solution to deal with this problem would be to design an effective method for the removal or reduction of the system, a controller. According to the literature, a conventional PID controller did not produce good results due to the high degree of nonlinearity of the process, fueling the development of advanced control techniques. Among these, the model predictive controller (MPC), where the control action results from the solution of an optimization problem, it is robust, can incorporate physical and /or security constraints. The objective of this work is to apply a non-conventional non-linear model predictive control technique to severe slugging, where the amount of liquid mass in the riser is controlled by the production valve and, indirectly, the oscillation of flow and pressure is suppressed, while looking for environmental and economic benefits. The proposed strategy is based on the use of the model linear approximations and repeatedly solving of a quadratic optimization problem, providing solutions that improve at each iteration. In the event where the convergence of this algorithm is satisfied, the predicted values of the process variables are the same as to those obtained by the original nonlinear model, ensuring that the constraints are satisfied for them along the prediction horizon. A mathematical model recently published in the literature, capable of representing characteristics of severe slugging in a real oil well, is used both for simulation and for the project of the proposed controller, whose performance is compared to a linear MPC

Resposta térmica de um compósito PEEK+PTFE+Fibra de carbono+grafite

Composites based on PEEK + PTFE + CARBON FIBER + Graphite (G_CFRP) has increased application in the top industries, as Aerospace, Aeronautical, Petroleum, Biomedical, Mechanical and Electronics Engineering challenges. A commercially available G_CFRP was warmed up to three different levels of thermal energy to identify the main damage mechanisms and some evidences for their intrinsic transitions. An experimental test rig for systematize a heat flux was developed in this dissertation, based on the Joule Effect. It was built using an isothermal container, an internal heat source and a real-time measurement system for test a sample by time. A standard conical-cylindrical tip was inserted into a soldering iron, commercially available and identified by three different levels of nominal electrical power, 40W (manufacturer A), 40W (manufacturer B), 100W and 150W, selected after screening tests: these power levels for the heat source, after one hour of heating and one hour of cooling in situ, carried out three different zones of degradation in the composite surface. The bench was instrumented with twelve thermocouples, a wattmeter and a video camera. The twelve specimens tested suffered different degradation mechanisms, analyzed by DSC (Differential Scanning Calorimetry) and TG (Thermogravimetry) techniques, Scanning Electron Microscopy (SEM) and Energy-Dispersive X-Rays (EDX) Analysis. Before and after each testing, it was measured the hardness of the sample by HRM (Hardness Rockwell M). Excellent correlations (R2=1) were obtained in the plots of the evaporated area after one hour of heating and one hour of cooling in situ versus (1) the respective power of heat source and (2) the central temperature of the sample. However, as resulting of the differential degradation of G_CFRP and their anisotropy, confirmed by their variable thermal properties, viscoelastic and plastic properties, there were both linear and non-linear behaviour between the temperature field and Rockwell M hardness measured in the radial and circumferential directions of the samples. Some morphological features of the damaged zones are presented and discussed, as, for example, the crazing and skeletonization mechanism of G_CFRP

An Anisotropic Model of Damage and Unilateral Effect for Brittle Materials

This work deals with the initial applications and formulation of an aniscitropic plastic-damage constitutive model proposed for non-linear analysis of reinforced concrete structures submitted to a loading with change of the sign. The original constitutive model is based on the fundamental hypothesis of energy equivalence between real and continuous medium following the concepts of the Continuum Damage Mechanics. The concrete is assumed as an initial elastic isotropic medium presenting anisotropy, permanent strains and bimodularity (distinct elastic responses whether traction or compression stress states prevail) induced by damage evolution. In order to take into account the bimodularity, two damage tensors governing the rigidity in tension or compression regimes are introduced. Then, some conditions are introduced in the original version of the model in order to simulate the damage unilateral effect. The three-dimensional version of the proposed model is analyzed in order to validate its formulation when compared to micromechanical theory. The one-dimensional version of the model is applied in the analyses of a reinforced concrete beam submitted to a loading with change of the sign. Despite the parametric identification problems, the initial applications show the good performance of the model.

Two-dimensional modeling of material failure in reinforced concrete by means of a continuum strong discontinuity approach

The paper presents a new methodology to model material failure, in two-dimensional reinforced concrete members, using the Continuum Strong Discontinuity Approach (CSDA). The mixture theory is used as the methodological approach to model reinforced concrete as a composite material, constituted by a plain concrete matrix reinforced with two embedded orthogonal long fiber bundles (rebars). Matrix failure is modeled on the basis of a continuum damage model, equipped with strain softening, whereas the rebars effects are modeled by means of phenomenological constitutive models devised to reproduce the axial non-linear behavior, as well as the bondslip and dowel effects. The proposed methodology extends the fundamental ingredients of the standard Strong Discontinuity Approach, and the embedded discontinuity finite element formulations, in homogeneous materials, to matrix/fiber composite materials, as reinforced concrete. The specific aspects of the material failure modeling for those composites are also addressed. A number of available experimental tests are reproduced in order to illustrate the feasibility of the proposed methodology. (c) 2007 Elsevier B.V. All rights reserved.

A mixture theory based method for three-dimensional modeling of reinforce concrete members with embedded crack finite elements

The paper presents a methodology to model three-dimensional reinforced concrete members by means of embedded discontinuity elements based on the Continuum Strong Discontinuous Approach (CSDA). Mixture theory concepts are used to model reinforced concrete as a 31) composite material constituted of concrete with long fibers (rebars) bundles oriented in different directions embedded in it. The effects of the rebars are modeled by phenomenological constitutive models devised to reproduce the axial non-linear behavior, as well as the bond-slip and dowel action. The paper presents the constitutive models assumed for the components and the compatibility conditions chosen to constitute the composite. Numerical analyses of existing experimental reinforced concrete members are presented, illustrating the applicability of the proposed methodology.

Attitude stability of artificial satellites subject to gravity gradient torque

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

Secular dynamics and family identification among highly inclined asteroids in the Euphrosyne region

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

The stable archipelago in the region of the Pallas and Hansa dynamical families

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

Dynamical erosion of asteroid groups in the region of the Pallas family

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

Conventional and microwave sintering of CaCu(3)Ti(4)O(12)/CaTiO(3) ceramic composites: non-ohmic and dielectric properties

The non-ohmic and dielectric properties as well as the dependence on the microstructural features of CaCu(3)Ti(4)O(12)/CaTiO(3) ceramic composites obtained by conventional and microwave sintering were investigated. It was demonstrated that the non-ohmic and dielectric properties depend strongly on the sintering conditions. It was found that the non-linear coefficient reaches values of 65 for microwave-sintered samples and 42 for samples sintered in a conventional furnace when a current density interval of 1-10 mA cm(-2) is considered. The non-linear coefficient value of 65 is equivalent to 1500 for samples sintered in the microwave if a current interval of 5-30 mA is considered as is shortly discussed by Chung et al (2004 Nature Mater. 3 774). Due to a high non-linear coefficient and a low leakage current (90 mu A) under both processing conditions, these samples are promising for varistor applications. The conventionally sintered samples exhibit a higher relative dielectric constant at 1 kHz (2960) compared with the samples sintered in the microwave furnace (2100). At high frequencies, the dielectric constant is also larger in the samples sintered in the conventional furnace. Depending on the application, one or another synthesis methodology is recommended, that is, for varistor applications sintered in a microwave furnace and for dielectric application sintered in a conventional furnace.