In this thesis, we present algorithms for local and global minimization of some Procrustes type problems. Typically, these problems are about rotating and scaling a known set of data to fit another set with applications related to determination of rigid body movements, factor analysis and multidimensional scaling. The known sets of data are usually represented as matrices, and the rotation to be determined is commonly a matrix Q with orthonormal columns. The algorithms presented use Newton and Gauss-Newton search directions with optimal step lengths, which in most cases result in a fast computation of a solution. Some of these problems are known to have several minima, e.g., the weighted orthogonal Procrustes problem (WOPP). A study on the maximal amount of minima has been done for this problem. Theoretical results and empirical observations gives strong indications that there are not more than 2n minimizers, where n is the number of columns in Q. A global optimization method to compute all 2n minima is presented. Also considered in this thesis is a cubically convergent iteration method for solving nonlinear equations. The iteration method presented uses second order information (derivatives) when computing a search direction. Normally this is a computational heavy task, but if the second order derivatives are constant, which is the case for quadratic equations, a performance gain can be obtained. This is confirmed by a small numerical study. Finally, regularization of ill-posed nonlinear least squares problems is considered. The quite well known L-curve for linear least squares problems is put in context for nonlinear problems.

The presence of wave motion is the defining feature in many fields of application,such as electro-magnetics, seismics, acoustics, aerodynamics,oceanography and optics. In these fields, accurate numerical simulation of wave phenomena is important for the enhanced understanding of basic phenomenon, but also in design and development of various engineering applications. In general, numerical simulations must be confined to truncated domains, much smaller than the physical space were the wave phenomena takes place. To truncate the physical space, artificial boundaries, and corresponding boundary conditions, are introduced. There are four main classes of methods that can be used to truncate problems on unbounded or large domains: boundary integral methods, infinite element methods, non-reflecting boundary condition methods and absorbing layer methods. In this thesis, we consider different aspects of non-reflecting boundary conditions and absorbing layers. In paper I, we construct discretely non-reflecting boundary conditions for a high order centered finite difference scheme. This is done by separating the numerical solution into spurious and physical waves, using the discrete dispersion relation. In paper II-IV, we focus on the perfectly matched layer method, which is a particular absorbing layer method. An open issue is whether stable perfectly matched layers can be constructed for a general hyperbolic system. In paper II, we present a stable perfectly matched layer formulation for 2 x 2 symmetric hyperbolic systems in (2 + 1) dimensions. We also show how to choose the layer parameters as functions of the coefficient matrices to guarantee stability. In paper III, we construct a new perfectly matched layer for the simulation of elastic waves in an anisotropic media. We present theoretical and numerical results, showing that the stability properties of the present layer are better than previously suggested layers. In paper IV, we develop general tools for constructing PMLs for first order hyperbolic systems. We present a model with many parameters which is applicable to all hyperbolic systems, and which we prove is well-posed and perfectly matched. We also use an automatic method, derived in paper V, for analyzing the stability of the model and establishing energy inequalities. We illustrate our techniques with applications to Maxwell s equations, the linearized Euler equations, as well as arbitrary 2 x 2 systems in (2 + 1) dimensions. In paper V, we use the method of Sturm sequences for bounding the real parts of roots of polynomials, to construct an automatic method for checking Petrowsky well-posedness of a general Cauchy problem. We prove that this method can be adapted to automatically symmetrize any well-posed problem, producing an energy estimate involving only local quantities. / QC 20100830

In this work CFD was used to simulate the flow through a novel valve design for internal combustion engines. CFD is numerical method for simulating the behaviour of systems involving flow processes. A FEM was used for solving the equations. Literature on the topic was studied to gain an understanding of the performance limiters on the Internal combustion engine. This understanding was used to set up models that better would mimic physical phenomena compared to previous studies. The models gave plausible results as to fluid velocities and in-cylinder flow patterns. Comsol Multiphysics 4.1 was used for the computations.

No description available.

<p>In this thesis the two-fluid (Eulerian/Eulerian) formulation for dispersed two-phase flow is considered. Closure laws are needed for this type of models. We investigate both empirically based relations, which we refer to as a nongranular model, and relations obtained from kinetic theory of dense gases, which we refer to as a granular model. For the granular model, a granular temperature is introduced, similar to thermodynamic temperature. It is often assumed that the granular energy is in a steady state, such that an algebraic granular model is obtained. </p><p>The inviscid non-granular model in one space dimension is known to be conditionally well-posed. On the other hand, the viscous formulation is locally in time well-posed for smooth initial data, but with a medium to high wave number instability. Linearizing the algebraic granular model around constant data gives similar results. In this study we consider a couple of issues. </p><p>First, we study the long time behavior of the viscous model in one space dimension, where we rely on numerical experiments, both for the non-granular and the algebraic granular model. We try to regularize the problem by adding second order artificial dissipation to the problem. The simulations suggest that it is not possible to obtain point-wise convergence using this regularization. Introducing a new measure, a concept of 1-D bubbles, gives hope for other convergence than point-wise. </p><p>Secondly, we analyse the non-granular formulation in two space dimensions. Similar results concerning well-posedness and instability is obtained as for the non-granular formulation in one space dimension. Investigation of the time scales of the formulation in two space dimension suggests a sever restriction on the time step, such that explicit schemes are impractical. </p><p>Finally, our simulation in one space dimension show that peaks or spikes form in finite time and that the solution is highly oscillatory. We introduce a model problem to study the formation and smoothness of these peaks.</p>

This thesis consists of an introduction and five papers concerning different numerical and mathematical aspects of viscous shocks. Hyperbolic conservation laws are used to model wave motion and advect- ive transport in a variety of physical applications. Solutions of hyperbolic conservation laws may become discontinuous, even in cases where initial and boundary data are smooth. Shock waves is one important type of discontinu- ity. It is also interesting to study the corresponding slightly viscous system, i.e., the system obtained when a small viscous term is added to the hyper- bolic system of equations. By a viscous shock we denote a thin transition layer which appears in the solution of the slightly viscous system instead of a shock in the corresponding purely hyperbolic problem. A slightly viscous system, a so called modified equation, is often used to model numerical solutions of hyperbolic conservation laws and their beha- vior in the vicinity of shocks. Computations presented elsewhere show that numerical solutions of hyperbolic conservation laws obtained by higher order accurate shock capturing methods in many cases are only first order accurate downstream of shocks. We use a modified equation to model numerical solu- tions obtained by a generic second order shock capturing scheme for a time dependent system in one space dimension. We present analysis that show how the first order error term is related to the viscous terms and show that it is possible to eliminate the first order downstream error by choosing a special viscosity term. This is verified in computations. We also extend the analysis to a stationary problem in two space dimensions. Though the technique of modified equation is widely used, rather little is known about when (for what methods etc.) it is applicable. The use of a modified equation as a model for a numerical solution is only relevant if the numerical solution behaves as a continuous function. We have experimentally investigated a range of high resolution shock capturing methods. Our experiments indicate that for many of the methods there is a continuous shock profile. For some of the methods, however, this not the case. In general the behavior in the shock region is very complicated. Systems of hyperbolic conservation laws with solutions containing shock waves, and corresponding slightly viscous equations, are examples where the available theoretical results on existence and uniqueness of solutions are very limited, though it is often straightforward to find approximate numerical solu- tions. We present a computer-assisted technique to prove existence of solu- tions of non-linear boundary value ODEs, which is based on using an approx- imate, numerical solution. The technique is applied to stationary solutions of the viscous Burgers' equation.We also study a corresponding method suggested by Yamamoto in SIAM J. Numer. Anal. 35(5)1998, and apply also this method to the viscous Burgers' equation.

In this thesis the two-fluid (Eulerian/Eulerian) formulation for dispersed two-phase flow is considered. Closure laws are needed for this type of models. We investigate both empirically based relations, which we refer to as a nongranular model, and relations obtained from kinetic theory of dense gases, which we refer to as a granular model. For the granular model, a granular temperature is introduced, similar to thermodynamic temperature. It is often assumed that the granular energy is in a steady state, such that an algebraic granular model is obtained. The inviscid non-granular model in one space dimension is known to be conditionally well-posed. On the other hand, the viscous formulation is locally in time well-posed for smooth initial data, but with a medium to high wave number instability. Linearizing the algebraic granular model around constant data gives similar results. In this study we consider a couple of issues. First, we study the long time behavior of the viscous model in one space dimension, where we rely on numerical experiments, both for the non-granular and the algebraic granular model. We try to regularize the problem by adding second order artificial dissipation to the problem. The simulations suggest that it is not possible to obtain point-wise convergence using this regularization. Introducing a new measure, a concept of 1-D bubbles, gives hope for other convergence than point-wise. Secondly, we analyse the non-granular formulation in two space dimensions. Similar results concerning well-posedness and instability is obtained as for the non-granular formulation in one space dimension. Investigation of the time scales of the formulation in two space dimension suggests a sever restriction on the time step, such that explicit schemes are impractical. Finally, our simulation in one space dimension show that peaks or spikes form in finite time and that the solution is highly oscillatory. We introduce a model problem to study the formation and smoothness of these peaks. / QC 20101018

<p>This thesis consists of an introduction and five papers concerning different numerical and mathematical aspects of viscous shocks. </p><p>Hyperbolic conservation laws are used to model wave motion and advect- ive transport in a variety of physical applications. Solutions of hyperbolic conservation laws may become discontinuous, even in cases where initial and boundary data are smooth. Shock waves is one important type of discontinu- ity. It is also interesting to study the corresponding slightly viscous system, i.e., the system obtained when a small viscous term is added to the hyper- bolic system of equations. By a viscous shock we denote a thin transition layer which appears in the solution of the slightly viscous system instead of a shock in the corresponding purely hyperbolic problem. </p><p>A slightly viscous system, a so called modified equation, is often used to model numerical solutions of hyperbolic conservation laws and their beha- vior in the vicinity of shocks. Computations presented elsewhere show that numerical solutions of hyperbolic conservation laws obtained by higher order accurate shock capturing methods in many cases are only first order accurate downstream of shocks. We use a modified equation to model numerical solu- tions obtained by a generic second order shock capturing scheme for a time dependent system in one space dimension. We present analysis that show how the first order error term is related to the viscous terms and show that it is possible to eliminate the first order downstream error by choosing a special viscosity term. This is verified in computations. We also extend the analysis to a stationary problem in two space dimensions. </p><p>Though the technique of modified equation is widely used, rather little is known about when (for what methods etc.) it is applicable. The use of a modified equation as a model for a numerical solution is only relevant if the numerical solution behaves as a continuous function. We have experimentally investigated a range of high resolution shock capturing methods. Our experiments indicate that for many of the methods there is a continuous shock profile. For some of the methods, however, this not the case. In general the behavior in the shock region is very complicated.</p><p>Systems of hyperbolic conservation laws with solutions containing shock waves, and corresponding slightly viscous equations, are examples where the available theoretical results on existence and uniqueness of solutions are very limited, though it is often straightforward to find approximate numerical solu- tions. We present a computer-assisted technique to prove existence of solu- tions of non-linear boundary value ODEs, which is based on using an approx- imate, numerical solution. The technique is applied to stationary solutions of the viscous Burgers' equation.We also study a corresponding method suggested by Yamamoto in SIAM J. Numer. Anal. 35(5)1998, and apply also this method to the viscous Burgers' equation.</p>

Gas sensors are today used in many different application areas, and one growing future market is battery operated sensors. As many gas sensor components are heated, one major limit of the operation time is caused by the power dissipated as heat. AppliedSensor is a company that develops and produces gas sensor components, modules and solutions, among which battery operated gas sensors are one targeted market. The aim of the diploma work has been to simulate the heat transfer on a hydrogen gas sensor component and its closest surroundings consisting of a carrier mounted on a printed circuit board. The component is heated in order to improve the performance of the gas sensing element. Power dissipation occurs by all three modes of heat transfer; conduction from the component through bond wires and carrier to the printed circuit board as well as convection and radiation from all the surfaces. It is of interest to AppliedSensor to understand which factors influence the heat transfer. This knowledge will be used to improve different aspects of the gas sensor, such as the power consumption. Modeling and simulation have been performed in FEMLAB, a tool for solving partial differential equations by the finite element method. The sensor system has been defined by the geometry and the material properties of the objects. The system of partial differential equations, consisting of the heat equation describing conduction and boundary conditions specifying convection and radiation, was solved and the solution was validated against experimental data. The convection increases with the increase of hydrogen concentration. A great effort was made to finding a model for the convection. Two different approaches were taken, the first based on known theory from the area and the second on experimental data. When the first method was compared to experiments, it turned out that the theory was insufficient to describe this small system involving hydrogen, which was an unexpected but interesting result. The second method matched the experiments well. For the continuation of the project at the company, a better model of the convection would be a great improvement.

<p>Gas sensors are today used in many different application areas, and one growing future market is battery operated sensors. As many gas sensor components are heated, one major limit of the operation time is caused by the power dissipated as heat. AppliedSensor is a company that develops and produces gas sensor components, modules and solutions, among which battery operated gas sensors are one targeted market.</p><p>The aim of the diploma work has been to simulate the heat transfer on a hydrogen gas sensor component and its closest surroundings consisting of a carrier mounted on a printed circuit board. The component is heated in order to improve the performance of the gas sensing element.</p><p>Power dissipation occurs by all three modes of heat transfer; conduction from the component through bond wires and carrier to the printed circuit board as well as convection and radiation from all the surfaces. It is of interest to AppliedSensor to understand which factors influence the heat transfer. This knowledge will be used to improve different aspects of the gas sensor, such as the power consumption.</p><p>Modeling and simulation have been performed in FEMLAB, a tool for solving partial differential equations by the finite element method. The sensor system has been defined by the geometry and the material properties of the objects. The system of partial differential equations, consisting of the heat equation describing conduction and boundary conditions specifying convection and radiation, was solved and the solution was validated against experimental data.</p><p>The convection increases with the increase of hydrogen concentration. A great effort was made to finding a model for the convection. Two different approaches were taken, the first based on known theory from the area and the second on experimental data. When the first method was compared to experiments, it turned out that the theory was insufficient to describe this small system involving hydrogen, which was an unexpected but interesting result. The second method matched the experiments well. For the continuation of the project at the company, a better model of the convection would be a great improvement.</p>