<p>Rational Krylov methods for model order reduction are studied. A dual rational Arnoldi method for model order reduction and a rational Krylov method for model order reduction and eigenvalue computation have been implemented. It is shown how to deflate redundant or unwanted vectors and how to obtain moment matching. Both methods are designed for generalised state space systems---the former for multiple-input-multiple-output (MIMO) systems from finite element discretisations and the latter for single-input-single-output (SISO) systems---and applied to relevant test problems. The dual rational Arnoldi method is designed for generating real reduced order systems using complex shift points and stabilising a system that happens to be unstable. For the rational Krylov method, a forward error in the recursion and an estimate of the error in the approximation of the transfer function are studie.</p><p>A stability analysis of a heat exchanger model is made. The model is a nonlinear partial differential-algebraic equation (PDAE). Its well-posedness and how to prescribe boundary data is investigated through analysis of a linearised PDAE and numerical experiments on a nonlinear DAE. Four methods for generating reduced order models are applied to the nonlinear DAE and compared: a Krylov based moment matching method, balanced truncation, Galerkin projection onto a proper orthogonal decomposition (POD) basis, and a lumping method.</p>

This thesis is concerned with the finite element analysis and the a posteriori error estimation for eigenvalue problems for general operator pencils on two-dimensional manifolds. A specific application of the presented theory is the computation of corner singularities. Engineers use the knowledge of the so-called singularity exponents to predict the onset and the propagation of cracks. All results of this thesis are explained for two model problems, the Laplace and the linear elasticity problem, and verified by numerous numerical results.

The research work presented in this thesis concern to the stability analysis of linear time-delay systems with low-order controllers. This thesis is divided into three parts.The first part of the thesis focus on the study of linear SISO (single-input/single-output) systems with input/output delays, where the feedback loop is closed with a controller of PID-type. Inspired by the geometrical approach developed by Gu et al. we propose an analytical method to find the stability regions of all stabilizing controllers of PID-type for the time-delay system. Based on this same approach, we propose an algorithm to calculate the degree of fragility of a given controller of PID- type (PI, PD and PID).The second part of the thesis focuses on the stability analysis of linear systems under an NCS (Networked System Control) based approach. More precisely, we first focus in the stabilization problem by taking into account the induced network delays and the effects induced by the sampling period. To carry out such an analysis we have adopted an eigenvalue perturbation-based approach.Finally, in the third part of the thesis we tackle certain problems concerning to the behavior of the zeros of a certain class of sampled-data SISO systems. More precisely, given a continuous-time system, we obtain the sampling intervals guaranteeing the invariance of the number of unstable zeros in each interval. To perform such an analysis, we adopt an eigenvalue perturbation-based approach.

The linear stability of high speed-viscous plane Couette and Couette-Poiseuille flows are investigated numerically. The conservation equations along with Sutherland&amp / #65533 / s viscosity law are studied using a second-order finite difference scheme. The basic velocity and temperature distributions are perturbed by a small-amplitude normalmode disturbance. The small-amplitude disturbance equations are solved numerically using a global method using QZ algorithm to find all the eigenvalues at finite Reynolds numbers, and the incompressible limit of these equations is investigated for Couette-Poiseuille flow. It is found that the instabilities occur, although the corresponding growth rates are often small. Two families of wave modes, Mode I (odd modes) and Mode II (even modes), were found to be unstable at finite Reynolds numbers, where Mode II is the dominant instability among the unstable modes for plane Couette flow. The most unstable mode for plane Couette &amp / #65533 / Poiseuille flow is Mode 0, which is not a member of the even modes. Both even and odd modes are acoustic modes created by acoustic reflections between a will and a relative sonic line. The necessary condition for the existence of such acoustic wave modes is that there is a region of locally supersonic mean flow relative to the phase speed of the instability wave. The effects of viscosity and compressibility are also investigated and shown to have a stabilizing role in all cases studied. Couette-Poiseuille flow stability is investigated in case of a choked channel flow, where the maximum velocity in the channel corresponds to sonic velocity. Neutral stability contours were obtained for this flow as a function if the wave number,Reynolds number and the upper wall Mach number. The critical Reynolds number is found as 5718.338 for an upper wall Mach number of 0.0001, corresponding to the fully Poiseuille case.

This thesis consists of two different parts, related to two different fields in mathematical physics: wave propagation and kinetic theory of gases. Various mathematical and computational problems for equations from these areas are treated.   The first part is devoted to high order finite difference methods for the Helmholtz equation and the wave equation. Compact schemes with high order accuracy are obtained from an investigation of the function derivatives in the truncation error. With the help of the equation itself, it is possible to transfer high order derivatives to lower order or to transfer time derivatives to space derivatives. For the Helmholtz equation, a compact scheme based on this principle is compared to standard schemes and to deferred correction schemes, and the characteristics of the errors for the different methods are demonstrated and discussed. For the wave equation, a finite difference scheme with fourth order accuracy in both space and time is constructed and applied to a problem in discontinuous media.   The second part addresses some problems related to kinetic equations. A direct simulation Monte-Carlo method is constructed for the Landau-Fokker-Planck equation, and numerical tests are performed to verify the accuracy of the algorithm. A formal derivation of the method from the Boltzmann equation with grazing collisions is performed. The linear and linearized Boltzmann collision operators for the hard sphere molecular model are studied using exact reduction of integral equations to ordinary differential equations. It is demonstrated how the eigenvalues of the operators are found from these equations, and numerical values are computed. A proof of existence of non-zero discrete eigenvalues is given. The ordinary diffential equations are also used for investigation of the Chapman-Enskog distribution function with respect to its asymptotic behavior.

This thesis investigates topology optimization techniques for periodic continuum structures at the macroscopic level. Periodic structures are increasingly used in the design of structural systems and sub-systems of buildings, vehicles, aircrafts, etc. The duplication of identical or similar modules significantly reduces the manufacturing cost and greatly simplifies the assembly process. Optimization of periodic structures in the micro level has been extensively researched in the context of material design, while research on topology optimization for macrostructures is very limited and has great potential both economically and intellectually. In the present thesis, numerical algorithms based on the bi-directional evolutionary structural optimization method (BESO) are developed for topology optimization for various objectives and constraints. Soft-kill (replacing void elements with soft elements) formulations of topology optimization problems for solid-void solutions are developed through appropriate material interpolation schemes. Incorporating the optimality criteria and algorithms for mesh-independence and solution-convergence, the present BESO becomes a reliable gradient based technique for topology optimization. Additionally, a new combination of genetic algorithms (GAs) with BESO is developed in order to stochastically search for the global optima. These enhanced BESO algorithms are applied to various optimization problems with the periodicity requirement as an extra constraint aiming at producing periodicity in the layout. For structures under static loading, the present thesis addresses minimization of the mean compliance and explores the applications of conventional stiffness optimization for periodic structures. Furthermore, this thesis develops a volume minimization formulation where the maximum deflection is constrained. For the design of structures subject to dynamic loading, this thesis develops two different approaches (hard-kill and soft-kill) to resolving the problem of localized or artificial modes. In the hard-kill (completely removing void elements) approach, extra control measures are taken in order to eliminate the localized modes in an explicit manner. In the soft-kill approach, a modified power low material model is presented to prevent the occurrence of artificial and localized modes. Periodic stress and strain fields cannot be assumed in structures under arbitrary loadings and boundaries at the macroscopic level. Therefore being different from material design, no natural base cell can be directly extracted from macrostructures. In this thesis, the concept of an imaginary representative unit cell (RUC) is presented. For situations when the structure cannot be discretized into equally-sized elements, the concept of sensitivity density is developed in order for mesh-independent robust solutions to be produced. The RUC and sensitivity density based approach is incorporated into various topology optimization problems to obtain absolute or scaled periodicities in structure layouts. The influence of this extra constraint on the final optima is investigated based on a large number of numerical experiments. The findings shown in this thesis have established appropriate techniques for designing and optimizing periodic structures. The work has provided a solid foundation for creating a practical design tool in the form of a user-friendly computer program suitable for the conceptual design of a wide range of structures.

Made available in DSpace on 2014-06-11T19:26:56Z (GMT). No. of bitstreams: 0 Previous issue date: 2009-02-27Bitstream added on 2014-06-13T20:16:04Z : No. of bitstreams: 1 nunes_jb_me_sjrp.pdf: 1005590 bytes, checksum: 7da54a97a1f2ab452a315062071f2c4e (MD5) / Este trabalho trata do estudo da localização dos zeros dos polinômios gerados por uma determinada relação de recorrência de três termos. O objetivo principal é estudar limitantes, em termos dos coeficientes da relação de recorrência, para as regiões onde os zeros estão localizados. Os zeros são explorados atravé do problema de autovalor associado a uma matriz de Hessenberg. As aplicações são consideradas para polinômios de Szego fSng, alguns polinômios para- ortogonais ½Sn(z) + S¤n (z) 1 + Sn(0) ¾ e ½Sn(z) ¡ S¤n (z) 1 ¡ Sn+1(0) ¾, especialmente quando os coeficientes de reflexão são reais. Um outro caso especial considerado são os zeros do polinômio Pn(z) = n Xm=0 bmzm, onde os coeficientes bm; para m = 0; 1; : : : ; n, são complexos e diferentes de zeros. / In this work we studied the localization the zeros of polynomials generated by a certain three term recurrence relation. The main objective is to study bounds, in terms of the coe±cients of the recurrence relation, for the regions where the zeros are located. The zeros are explored through an eigenvalue representation associated with a Hessenberg matrix. Applications are considered to Szeg}o polynomials fSng, some para-orthogonal polyno- mials ½Sn(z) + S¤n (z) 1 + Sn(0) ¾and ½Sn(z) ¡ S¤n (z) 1 ¡ Sn+1(0) ¾, especially when the re°ection coe±cients are real. As another special case, the zeros of the polynomial Pn(z) = n Xm=0 bmzm, where the non-zero complex coe±cients bm for m = 0; 1; : : : ; n, were considered.

Title: Numerical modeling of free oscillations applied to superconducting-gravimeter data in a low-frequency seismic range Author: Eliška Zábranová Department: Department of Geophysics Supervisor: Doc. RNDr. Ctirad Matyska, DrSc. Abstract: Deformations and changes of the gravitational potential of prestressed selfgravitating elastic bodies caused by free oscillations are described by means of the momentum and Poisson equations and the constitutive relation. For spheri- cally symmetric bodies we transform the equations and boundary conditions into ordinary differential equations of the second order by the spherical harmonic de- composition and further discretize the equations by highly accurate pseudospectral difference schemes on Chebyshev grids. We thus receive a series of matrix eigenvalue problems for eigenfrequencies and eigenfunctions of the free oscillations. Since elas- tic parameters are frequency dependent, we solve the problem for several fiducial frequencies and interpolate the results. Both the mode frequencies and the eigen- functions are benchmarked against the output from the Mineos software package based on Runge-Kutta integration techniques. Subsequently, we use our method to calculate low-frequency synthetic accelerograms of the recent megathrust events and compare them with the observed...

Brake squeal is an old and well-known problem in the vehicle industry and is a frequent source for customer complain. Although, brake squeal is not usually affecting the performance of the brakes, it is still important to address the problem and to predict the brakes tendency to squeal on an early stage in the design process. Brake squeal is usually defined as a sustained, high-frequency vibration of the brake components, due to the braking action. By using simulation in finite element (FE) method it should be possible to predict at what frequencies the brakes tend to emit sound. The method chosen for the analysis was the complex eigenvalues analysis (CEA) method, since it is a well-known tool to predict unstable modes in FE analysis. The results from the CEA were evaluated against measured data from an earlier study. Even though there are four main mechanism formulated in order to explain the up come of squeal, the main focus in this project was modal coupling, since it is the main mechanism in the CEA. A validation of the key components in model was performed before the analysis, in order to achieve better correlation between the FE model and reality. A parametric study was conducted with the CEA, to investigate how material properties and operating parameters effected the brakes tendency to squeal. The following parameters was included in the analysis; coefficient of friction, brake force, damping, rotational velocity, and Young’s modulus for different components. The result from the CEA did not exactly reproduce the noise frequencies captured in experimental tests. The discrepancy is believed to mainly be due to problems in the calibration process of the components in the model. The result did however show that the most effective way to reduce the brakes tendency for squeal was to lower the coefficient of friction. The effect of varying the Young’s modulus different components showed inconsistent results on the tendency to squeal. By adding damping one of the main disadvantages for the CEA, which the over-prediction of the number of unstable modes, where minimized.

Bound state systems consisting of three nonrelativistic particles are numerically studied. Calculations are performed employing two-body and three-body forces as input in the Hamiltonian in order to study the role or contribution of three-body forces to the binding in these systems. The resulting differential Faddeev equations are solved as three-dimensional equations in the two Jacobi coordinates and the angle between them, as opposed to the usual partial wave expansion approach. By expanding the wave function as a sum of the products of spline functions in each of the three coordinates, and using the orthogonal collocation procedure, the equations are transformed into an eigenvalue problem. The matrices in the aforementioned eigenvalue equations are generally of large order. In order to solve these matrix equations with modest and optimal computer memory and storage, we employ the iterative Restarted Arnoldi Algorithm in conjunction with the so-called tensor trick method. Furthermore, we incorporate a polynomial accelerator in the algorithm to obtain rapid convergence. We applied the method to obtain the binding energies of Triton, Carbon-12, and Ozone molecule. / Physics / M.Sc (Physics)