Implementação paralela de um código de elementos finitos em 2D para as Equações de Navier-Stokes para fluidos incompressíveis com transporte de escalares.

O estudo do fluxo de água e do transporte escalar em reservatórios hidrelétricos é importante para a determinação da qualidade da água durante as fases iniciais do enchimento e durante a vida útil do reservatório. Neste contexto, um código de elementos finitos paralelo 2D foi implementado para resolver as equações de Navier-Stokes para fluido incompressível acopladas a transporte escalar, utilizando o modelo de programação de troca de mensagens, a fim de realizar simulações em um ambiente de cluster de computadores. A discretização espacial é baseada no elemento MINI, que satisfaz as condições de Babuska-Brezzi (BB), que permite uma formulação mista estável. Todas as estruturas de dados distribuídos necessárias nas diferentes fases do código, como pré-processamento, solução e pós-processamento, foram implementadas usando a biblioteca PETSc. Os sistemas lineares resultantes foram resolvidos usando o método da projeção discreto com fatoração LU por blocos. Para aumentar o desempenho paralelo na solução dos sistemas lineares, foi empregado o método de condensação estática para resolver a velocidade intermediária nos vértices e no centróide do elemento MINI separadamente. Os resultados de desempenho do método de condensação estática com a abordagem da solução do sistema completo foram comparados. Os testes mostraram que o método de condensação estática apresenta melhor desempenho para grandes problemas, às custas de maior uso de memória. O desempenho de outras partes do código também são apresentados.

Abordagens do tipo livre de jacobiana na simulação do escoamento de fluidos compressíveis em meios porosos

O desenvolvimento de software livre de Jacobiana para a resolução de problemas formulados por equações diferenciais parciais não-lineares é de interesse crescente para simular processos práticos de engenharia. Este trabalho utiliza o chamado algoritmo espectral livre de derivada para equações não-lineares na simulação de fluxos em meios porosos. O modelo aqui considerado é aquele empregado para descrever o deslocamento do fluido compressível miscível em meios porosos com fontes e sumidouros, onde a densidade da mistura de fluidos varia exponencialmente com a pressão. O algoritmo espectral utilizado é um método moderno para a solução de sistemas não-lineares de grande porte, o que não resolve sistemas lineares, nem usa qualquer informação explícita associados com a matriz Jacobiana, sendo uma abordagem livre de Jacobiana. Problemas bidimensionais são apresentados, juntamente com os resultados numéricos comparando o algoritmo espectral com um método de Newton inexato livre de Jacobiana. Os resultados deste trabalho mostram que este algoritmo espectral moderno é um método confiável e eficiente para a simulação de escoamentos compressíveis em meios porosos.

Precondicionadores baseados na aproximação da inversa.

Neste trabalho de dissertação apresentaremos uma classe de precondicionadores baseados na aproximação esparsa da inversa da matriz de coecientes, para a resolução de sistemas lineares esparsos de grandes portes através de métodos iterativos, mais especificamente métodos de Krylov. Para que um método de Krylov seja eficiente é extremamente necessário o uso de precondicionadores. No contexto atual, onde computadores de arquitetura híbrida são cada vez mais comuns temos uma demanda cada vez maior por precondicionadores paralelizáveis. Os métodos de inversa aproximada que serão descritos possuem aplicação paralela, pois so dependem de uma operação de produto matriz-vetor, que é altamente paralelizável. Além disso, alguns dos métodos também podem ser construídos em paralelo. A ideia principal é apresentar uma alternativa aos tradicionais precondicionadores que utilizam aproximações dos fatores LU, que apesar de robustos são de difícil paralelização.

Balanced parametrizations and canonical forms for system identification

This paper contains a review of recent results concerning the parametrization of asymptotically stable linear systems using balanced realizations. Particular emphasis is given on the application of these results to system identification. This work is part of a continuing programme aimed at elucidating the role of balanced realization in system identification.

The changing face and role of CACSD

Computer Aided Control Engineering involves three parallel streams: Simulation and modelling, Control system design (off-line), and Controller implementation. In industry the bottleneck problem has always been modelling, and this remains the case - that is where control (and other) engineers put most of their technical effort. Although great advances in software tools have been made, the cost of modelling remains very high - too high for some sectors. Object-oriented modelling, enabling truly re-usable models, seems to be the key enabling technology here. Software tools to support control systems design have two aspects to them: aiding and managing the work-flow in particular projects (whether of a single engineer or of a team), and provision of numerical algorithms to support control-theoretic and systems-theoretic analysis and design. The numerical problems associated with linear systems have been largely overcome, so that most problems can be tackled routinely without difficulty - though problems remain with (some) systems of extremely large dimensions. Recent emphasis on control of hybrid and/or constrained systems is leading to the emerging importance of geometric algorithms (ellipsoidal approximation, polytope projection, etc). Constantly increasing computational power is leading to renewed interest in design by optimisation, an example of which is MPC. The explosion of embedded control systems has highlighted the importance of autocode generation, directly from modelling/simulation products to target processors. This is the 'new kid on the block', and again much of the focus of commercial tools is on this part of the control engineer's job. Here the control engineer can no longer ignore computer science (at least, for the time being). © 2006 IEEE.

The changing face and role of CACSD

Computer Aided Control Engineering involves three parallel streams: Simulation and modelling, Control system design (off-line), and Controller implementation. In industry the bottleneck problem has always been modelling, and this remains the case - that is where control (and other) engineers put most of their technical effort. Although great advances in software tools have been made, the cost of modelling remains very high - too high for some sectors. Object-oriented modelling, enabling truly re-usable models, seems to be the key enabling technology here. Software tools to support control systems design have two aspects to them: aiding and managing the work-flow in particular projects (whether of a single engineer or of a team), and provision of numerical algorithms to support control-theoretic and systems-theoretic analysis and design. The numerical problems associated with linear systems have been largely overcome, so that most problems can be tackled routinely without difficulty - though problems remain with (some) systems of extremely large dimensions. Recent emphasis on control of hybrid and/or constrained systems is leading to the emerging importance of geometric algorithms (ellipsoidal approximation, polytope projection, etc). Constantly increasing computational power is leading to renewed interest in design by optimisation, an example of which is MPC. The explosion of embedded control systems has highlighted the importance of autocode generation, directly from modelling/simulation products to target processors. This is the 'new kid on the block', and again much of the focus of commercial tools is on this part of the control engineer's job. Here the control engineer can no longer ignore computer science (at least, for the time being). ©2006 IEEE.

H∞ design procedure using robust stabilization of normalized coprime factors

A two-stage H∞-based design procedure is described which uses a normalized coprime factor approach to robust stabilization of linear systems. A loop-shaping procedure is incroporated to allow the specification of performance characteristics. Theoretical justification of this technique and an outline of the design methodology are given.

Predictive control of a Boeing 747 aircraft using an FPGA

New embedded predictive control applications call for more eficient ways of solving quadratic programs (QPs) in order to meet demanding real-time, power and cost requirements. A single precision QP-on-a-chip controller is proposed, implemented in afield-programmable gate array (FPGA) with an iterative linear solver at its core. A novel offline scaling procedure is introduced to aid the convergence of the reduced precision solver. The feasibility of the proposed approach is demonstrated with a real-time hardware-in-the-loop (HIL) experimental setup where an ML605 FPGA board controls a nonlinear model of a Boeing 747 aircraft running on a desktop PC through an Ethernet link. Simulations show that the quality of the closed-loop control and accuracy of individual solutions is competitive with a conventional double precision controller solving linear systems using a Riccati recursion. © 2012 IFAC.

Power absorption invariance for brownian spring forcing

In [5] it was shown that, for a standard quarter-car vehicle model and a road disturbance whose velocity profile is white noise of intensity A, the mean power dissipated in the suspension is equal to kA/2 where k is the tyre vertical stiffness. It is remarkable that the power dissipation turns out to be independent of all masses and suspension parameters. The proof in [5] makes use of a spectral formulation of white noise and is specific to linear systems. This paper casts the result in a more general form and shows that it follows from a simple application of Ito calculus. © 2012 IEEE.

Stability margins of nonlinear receding-horizon control via inverse optimality

Using the nonlinear analog of the Fake Riccati equation developed for linear systems, we derive an inverse optimality result for several receding-horizon control schemes. This inverse optimality result unifies stability proofs and shows that receding-horizon control possesses the stability margins of optimal control laws. © 1997 Elsevier Science B.V.

What kind of friction model is needed to predict brake squeal?

Predictive models of friction-induced vibration have proved elusive despite decades of research. There are many mechanisms that can cause brake squeal; friction coupled systems can be highly sensitive to small perturbations; and the dynamic properties of friction at the contact zone seem to be poorly understood. This paper describes experimental and theoretical work aimed at identifying the key ingredients of a predictive model. A large-scale experiment was carried out to identify squeal initiations using a pin-on-disc test rig: approximately 30,000 squeal initiations were recorded, covering a very wide range of frequencies. The theoretical model allows for completely general linear systems coupled at a single sliding point by friction: squeal is predicted using a linearised stability analysis. Results will be presented that show that almost all observed squeal events can be predicted within this model framework, but that some subsets require innovative friction modelling: predictions are highly dependent on the particular choice of friction model and its associated parameters. Copyright © 2012 by ASME.

What kind of friction model is needed to predict brake squeal?

Predictive models of friction-induced vibration have proved elusive despite decades of research. There are many mechanisms that can cause brake squeal; friction coupled systems can be highly sensitive to small perturbations; and the dynamic properties of friction at the contact zone seem to be poorly understood. This paper describes experimental and theoretical work aimed at identifying the key ingredients of a predictive model. A large-scale experiment was carried out to identify squeal initiations using a pin-on-disc test rig: approximately 30,000 squeal initiations were recorded, covering a very wide range of frequencies. The theoretical model allows for completely general linear systems coupled at a single sliding point by friction: squeal is predicted using a linearised stability analysis. Results will be presented that show that almost all observed squeal events can be predicted within this model framework, but that some subsets require innovative friction modelling: predictions are highly dependent on the particular choice of friction model and its associated parameters. Copyright © 2012 by ASME.

Modeling of cyanobacterial blooms in hypereutrophic Lake Dianchi, China

Compared with other approaches for modeling and predicting, artificial neural networks are more effective in describing complex and non-linear systems. The occurrence of cyanobacterial blooms has been a continuous and serious problem over the past decades in hypereutrophic Lake Dianchi. Yet, the main factor(s) initiating these blooms remain(s) unclear. During 2001-2002 at 40 sampling sites in Lake Dianchi, physicochemical parameters possibly relating to the blooms were measured. Parameters directly or indirectly relating to the cyanobacterial blooms were used as driving factors in a back-propagation network to model the concentration of chlorophyll a. According to sensitivity analysis, chemical oxygen demand was identified as a very significant environmental factor for algal growth in Lake Dianchi.

轮式移动机器人沙地行驶控制建模与仿真研究

结合车轮沙土相互作用的力学分析,研究解决轮式移动机器人沙地行驶车轮过度滑转下陷问题。考虑纵列式重复通过车轮沙土力学参数的修正,建立六轮式沙地移动机器人的动力学模型,以车轮滑转率为状态变量,设计了移动机器人沙地行驶的滑模驱动控制器跟踪车轮期望滑转率。MATLAB/Simulink仿真结果表明,采用该控制器可以较快地跟踪期望滑转率,有效地限制机器人车轮的滑转,避免车轮的过度下陷。