KLUEDO RSS FeedKLUEDO Dokumente/documents
https://kluedo.ub.uni-kl.de/index/index/
Wed, 21 Jul 2010 14:10:41 +0200Wed, 21 Jul 2010 14:10:41 +0200Identification of Temperature Dependent Parameters in Radiative Heat Transfer
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2216
Laser-induced thermotherapy (LITT) is an established minimally invasive percutaneous technique of tumor ablation. Nevertheless, there is a need to predict the effect of laser applications and optimizing irradiation planning in LITT. Optical attributes (absorption, scattering) change due to thermal denaturation. The work presents the possibility to identify these temperature dependent parameters from given temperature measurements via an optimal control problem. The solvability of the optimal control problem is analyzed and results of successful implementations are shown.Oliver Tse; René Pinnau; Norbert Siedowpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/2216Wed, 21 Jul 2010 14:10:41 +0200Optimal Control of Melt Spinning Processes
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1928
An optimal control problem for a mathematical model of a melt spinning process is considered. Newtonian and non--Newtonian models are used to describe the rheology of the polymeric material, the fiber is made of. The extrusion velocity of the polymer at the spinneret as well as the velocity and temperature of the quench air serve as control variables. A constrained optimization problem is derived and the first--order optimality system is set up to obtain the adjoint equations. Numerical solutions are carried out using a steepest descent algorithm.Thomas Götz; Shyam Pererapreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1928Sun, 20 Jan 2008 12:37:46 +0100Optimal Control of Film Casting Processes
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1927
We present an optimal control approach for the isothermal film casting process with free surfaces described by averaged Navier-Stokes equations. We control the thickness of the film at the take-up point using the shape of the nozzle. The control goal consists in finding an even thickness profile. To achieve this goal, we minimize an appropriate cost functional. The resulting minimization problem is solved numerically by a steepest descent method. The gradient of the cost functional is approximated using the adjoint variables of the problem with fixed film width. Numerical simulations show the applicability of the proposed method.Thomas Götz; K. Selvanayagampreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1927Sun, 20 Jan 2008 12:37:37 +0100Convergent Finite Element Discretizations of the Density Gradient Equation for Quantum Semiconductors
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1864
We study nonlinear finite element discretizations for the density gradient equation in the quantum drift diffusion model. Especially, we give a finite element description of the so--called nonlinear scheme introduced by {it Ancona}. We prove the existence of discrete solutions and provide a consistency and convergence analysis, which yields the optimal order of convergence for both discretizations. The performance of both schemes is compared numerically, especially with respect to the influence of approximate vacuum boundary conditions.Rene Pinnau; Jorge Mauricio Ruizpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1864Fri, 11 May 2007 20:46:05 +0200The Semiconductor Model Hierarchy in Optimal Dopant Profiling
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1781
We consider optimal design problems for semiconductor devices which are simulated using the energy transport model. We develop a descent algorithm based on the adjoint calculus and present numerical results for a ballistic diode. Further, we compare the optimal doping profile with results computed on basis of the drift diffusion model. Finally, we exploit the model hierarchy and test the space mapping approach, especially the aggressive space mapping algorithm, for the design problem. This yields a significant reduction of numerical costs and programming effort.Concetta Drago; Rene Pinnaupreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1781Tue, 26 Sep 2006 19:06:17 +0200Model Reduction Techniques for Frequency Averaging in Radiative Heat
Transfer
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1738
We study model reduction techniques for frequency averaging in radiative heat transfer. Especially, we employ proper orthogonal decomposition in combination with the method of snapshots to devise an automated a posteriori algorithm, which helps to reduce significantly the dimensionality for further simulations. The reliability of the surrogate models is tested and we compare the results with two other reduced models, which are given by the approximation using the weighted sum of gray gases and by an frequency averaged version of the so-called \(\mathrm{SP}_n\) model. We present several numerical results underlining the feasibility of our approach.Rene Pinnau; Alexander Schulzepreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1738Tue, 16 May 2006 14:12:20 +0200Regularized Fixed-Point Iterations for Nonlinear Inverse Problems
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1653
In this paper we introduce a derivative-free, iterative method for solving nonlinear ill-posed problems \(Fx=y\), where instead of \(y\) noisy data \(y_\delta\) with \(|| y-y_\delta ||\leq \delta\) are given and \(F:D(F)\subseteq X \rightarrow Y\) is a nonlinear operator between Hilbert spaces \(X\) and \(Y\). This method is defined by splitting the operator \(F\) into a linear part \(A\) and a nonlinear part \(G\), such that \(F=A+G\). Then iterations are organized as \(A u_{k+1}=y_\delta-Gu_k\). In the context of ill-posed problems we consider the situation when \(A\) does not have a bounded inverse, thus each iteration needs to be regularized. Under some conditions on the operators \(A\) and \(G\) we study the behavior of the iteration error. We obtain its stability with respect to the iteration number \(k\) as well as the optimal convergence rate with respect to the noise level \(\delta\), provided that the solution satisfies a generalized source condition. As an example, we consider an inverse problem of initial temperature reconstruction for a nonlinear heat equation, where the nonlinearity appears due to radiation effects. The obtained iteration error in the numerical results has the theoretically expected behavior. The theoretical assumptions are illustrated by a computational experiment.S.S. Pereverzyev; R. Pinnau; N. Siedowpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1653Wed, 10 Aug 2005 17:10:20 +0200Initial Temperature Reconstruction for a Nonlinear Heat Equation: Application to Radiative Heat Transfer
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1607
Consider a cooling process described by a nonlinear heat equation. We are interested to recover the initial temperature from temperature measurements which are available on a part of the boundary for some time. Up to now even for the linear heat equation such a problem has been usually studied as a nonlinear ill-posed operator equation, and regularization methods involving Frechet derivatives have been applied. We propose a fast derivative-free iterative method. Numerical results are presented for the glass cooling process, where nonlinearity appears due to radiation.Sergiy Pereverzyev; Rene Pinnau; Norbert Siedowpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1607Fri, 04 Feb 2005 10:33:11 +0100On an asymptotic expansion for porous media flow of Carreau fluids
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1568
Porous media flow of polymers with Carreau law viscosities and their application to enhanced oil recovery (EOR) is considered. Applying the homogenization method leads to a nonlinear two-scale problem. In case of a small difference between the Carreau and the Newtonian case an asymptotic expansion based on the small deviation of the viscosity from the Newtonian case is introduced. For uni-directional pressure gradients, which is a reasonable assumption in applications like EOR, auxiliary problems to decouple the micro- from the macrovariables are derived. The microscopic flow field obtained by the proposed approach is compared to the solution of the two-scale problem. Finite element calculations for an isotropic and an anisotropic pore cell geometries are used to validate the accuracy and speed-up of the proposed approach. The order of accuracy has been studied by performing the simulations up to the third order expansion for the isotropic geometry.Thomas Götz; Hanna Parhusippreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1568Thu, 09 Sep 2004 19:06:00 +0200Wave Reflection and Refraction in Triclinic Crystalline media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1414
In this paper, the reflection and refraction of a plane wave at an interface between .two half-spaces composed of triclinic crystalline material is considered. It is shown that due to incidence of a plane wave three types of waves namely quasi-P (qP), quasi-SV (qSV) and quasi-SH (qSH) will be generated governed by the propagation condition involving the acoustic tensor. A simple procedure has been presented for the calculation of all the three phase velocities of the quasi waves. It has been considered that the direction of particle motion is neither parallel nor perpendicular to the direction of propagation. Relations are established between directions of motion and propagation, respectively. The expressions for reflection and refraction coefficients of qP, qSV and qSH waves are obtained. Numerical results of reflection and refraction coefficients are presented for different types of anisotropic media and for different types of incident waves. Graphical representation have been made for incident qP waves and for incident qSV and qSH waves numerical data are presented in two tables.Amares Chattopadhyaypreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1414Tue, 15 Jul 2003 13:10:48 +0200Regularized Multiresolution Recovery of the Mass Density Distribution from Satellite Data of the Earth's Gravitational Field
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1413
The inverse problem of recovering the Earth's density distribution from satellite data of the first or second derivative of the gravitational potential at orbit height is discussed. This problem is exponentially ill-posed. In this paper a multiscale regularization technique using scaling functions and wavelets constructed for the corresponding integro-differential equations is introduced and its numerical applications are discussed. In the numerical part the second radial derivative of the gravitational potential at 200 km orbit height is calculated on a point grid out of the NASA/GSFC/NIMA Earth Geopotential Model (EGM96). Those simulated derived data out of SGG satellite measurements are taken for convolutions with the introduced scaling functions yielding a multiresolution analysis of harmonic density variations in the Earth's crust.Volker Michelpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1413Tue, 15 Jul 2003 13:04:31 +0200A Model for Spherical SH-Wave Propagation in Self-reinforced Linearly Elastic Media
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1412
The original publication is available at www.springerlink.com. This original publication also contains further results. We study a spherical wave propagating in radius- and latitude-direction and oscillating in latitude-direction in case of fibre-reinforced linearly elastic material. A function system solving Euler's equation of motion in this case and depending on certain Bessel and associated Legendre functions is derived.Amares Chattopadhyay; Volker Michelpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1412Tue, 15 Jul 2003 10:45:21 +0200Smoothing Splines in Multiscale Geopotential Determination from Satellite Data
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1396
SST (satellite-to-satellite tracking) and SGG (satellite gravity gradiometry) provide data that allows the determination of the first and second order radial derivative of the earth's gravitational potential on the satellite orbit, respectively. The modeling of the gravitational potential from such data is an exponentially ill-posed problem that demands regularization. In this paper, we present the numerical studies of an approach, investigated in [24] and [25], that reconstructs the potential with spline smoothing. In this case, spline smoothing is not just an approximation procedure but it solves the underlying compact operator equation of the SST-problem and the SGG-problem. The numerical studies in this paper are performed for a simplified geometrical scenario with simulated data, but the approach is designed to handle first or second order radial derivative data on a real satellite orbit.Kerstin Hesse; Martin Guttingpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1396Mon, 26 May 2003 10:49:20 +0200SST-Regularisierung durch Multiresolutionstechniken auf der Basis von CHAMP-Daten
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1377
Die Bestimmung des Erdgravitationspotentials aus den Meßdaten des Forschungssatelliten CHAMP lässt sich als Operatorgleichung formulieren (SST-Problem). Dieser Ansatz geht davon aus, dass ein geometrischer Orbit des Satelliten CHAMP vorliegt. Mittels numerischer Differentiation unter Einsatz eines geeigneten Denoising Verfahrens kann dann aus dem geometrischen Orbit der Gradient des Potentials längs der Bahn bestimmt werden. Damit sind insbesondere die Radialableitung (und der Flächengradient) auf einem Punktgitter auf der Bahn bekannt. In einem erdfesten System stellt sich dies als eine nahezu vollständige Überdeckung der Erde (bis auf Polar Gaps) mit einem ziemlich dichten Datengitter auf Flughöhe des Satelliten dar. Die Lösung der SST-Operatorgleichung (Bestimmung des Potentials auf der Erdoberfläche aus Kenntnis der Radialableitung auf einem Datengitter auf Flughöhe) ist ein schlecht gestelltes inverses Problem, das mit einer geeigneten Regularisierungstechnik gelöst werden muß. Im vorliegenden Fall wurde eine solche Regularisierung mit Hilfe von nicht-bandlimitierten Regularisierungsskalierungsfunktionen und Regularisierungswavelets umgesetzt. Diese sind stark ortslokalisierend und führen daher auf ein Potentialmodell, welches eine Linearkombination stark ortslokalisierender Funktionen ist. Ein solches Modell kann als Lokalmodell auch aus nur lokalen Daten berechnet werden und bietet daher gegenüber Kugelfunktionsmodellen wie EGM96 erhebliche Vorteile für die moderne Geopotentialbestimmung. Die Diskretisierung und numerische Umsetzung der Berechnung eines solchen Modells erfolgt mit Splines, die hier ebenfalls Linearkombinationen stark ortslokalisierender Funktionen sind. Die großen linearen Gleichungssysteme, die zur Berechnung der glättenden oder interpolierenden Splines gelöst werden müssen, können auf schnelle und effiziente Weise mit dem Schwarzschen alternierenden Algorithmus in Verbindung mit schnellen Summationsverfahren (Fast Multipole Methods) gelöst werden. Eine Kombination des Schwarzschen alternierenden Algorithmus mit solchen schnellen Summationsverfahren ermöglicht eine weitere erhebliche Beschleunigung beim Lösen dieser Gleichungssysteme. Zur Bestimmung von Glättungsparametern (Spline-Smoothing) und Regularisierungsparametern kann die L-Curve Method zum Einsatz kommen.Willi Freeden; Thorsten Maierpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1377Mon, 24 Feb 2003 14:15:34 +0100Multiresolution Data Analysis - Numerical Realization by use of Domain Decomposition Methods and Fast Multipole Techniques / Multiscale Solutions of Oblique Boundary-Value Problems by Layer Potentials
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1344
This survey paper deals with multiresolution analysis from geodetically relevant data and its numerical realization for functions harmonic outside a (Bjerhammar) sphere inside the Earth. Harmonic wavelets are introduced within a suit- able framework of a Sobolev-like Hilbert space. Scaling functions and wavelets are defined by means of convolutions. A pyramid scheme provides efficient implementation und economical computation. Essential tools are the multiplicative Schwarz alternating algorithm (providing domain decomposition procedures) and fast multipole techniques (accelerating iterative solvers of linear systems).W. Freeden; C. Mayerpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1344Thu, 12 Sep 2002 10:37:46 +0200New Numerical Schemes based on Relaxation Systems for Conservation Laws
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1337
S. V. Rao Raghuramapreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1337Fri, 19 Jul 2002 00:00:00 +0200Spline Modelling of Geostrophic Flow: Theoretical and Algorithmic Aspects
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1338
Spline functions that approximate (geostrophic) wind field or ocean circulation data are developed in a weighted Sobolev space setting on the (unit) sphere. Two problems are discussed in more detail: the modelling of the (geostrophic) wind field from (i)discrete scalar air pressure data and (ii) discrete vectorial velocity data. Domain decomposition methods based on the Schwarz alternating algorithm for positive definite symmetric matrices are described for solving large linear systems occuring in vectorial spline interpolation or smoothing of geostrophic flow.W. Freeden; K. Hessepreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1338Fri, 19 Jul 2002 00:00:00 +0200Introduction to Systems and Control Theory
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1336
These lecture notes give a completely self-contained introduction to the control theory of linear time-invariant systems. No prior knowledge is requried apart from linear algebra and some basic familiarity with ordinary differential equations. Thus, the course is suited for students of mathematics in their second or third year, and for theoretically inclined engineering students. Because of its appealing simplicity and elegance, the behavioral approch has been adopted to a large extend. A short list of recommended text books on the subject has been added, as a suggestion for further reading.Eva Zerzpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1336Mon, 01 Jul 2002 00:00:00 +0200Multiscale Deformation Analysis by Cauchy-Navier Wavelets
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1281
A geoscientifically relevant wavelet approach is established for the classical (inner) displacement problem corresponding to a regular surface (such as sphere, ellipsoid, actual earth's surface). Basic tools are the limit and jump relations of (linear) elastostatics. Scaling functions and wavelets are formulated within the framework of the vectorial Cauchy-Navier equation. Based on appropriate numerical integration rules a pyramid scheme is developed providing fast wavelet transform (FWT). Finally multiscale deformation analysis is investigated numerically for the case of a spherical boundary.M. K. Abeyratne; W. Freeden; C. Mayerpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1281Thu, 21 Mar 2002 00:00:00 +0100Maximum Entropy Moment Systems and Galilean Invariance
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1276
In this article, we investigate the maximum entropy moment closure in gas dynamics. We show that the usual choice of polynomial weight functions may lead to hyperbolic systems with an unpleasant state space: equilibrium states are boundary points with possibly singular fluxes. In order to avoid singularities, the necessary arises to find weight functions which growing sub-quadratically at infinity. Unfortunately, this requirement leads to a conflict with Galilean invariance of the moment systems because we can show that rotational and translational invariant, finite dimensional function spaces necessarily consist of polynomials.Michael Junk; Andreas Unterreiterpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1276Tue, 11 Dec 2001 00:00:00 +0100Hyperbolic Conservation Laws and Industrial Applications
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1270
Michael Junkpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1270Mon, 10 Dec 2001 00:00:00 +0100A limiter based on kinetic theory
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1271
Mapundi K. Banda; Michael Junk; Axel Klarpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1271Mon, 10 Dec 2001 00:00:00 +0100Do Finite Volume Methods Need a Mesh?
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1272
In this article, finite volume discretizations of hyperbolic conservation laws are considered, where the usual triangulation is replaced of unity on the computational domain.Michael Junkpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1272Mon, 10 Dec 2001 00:00:00 +0100Rigorous Navier-Stokes Limit of the Lattice Boltzmann Equation
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1273
Here we riqorously investigate the diffusive limit of a velocity-discrete Boltzmann equation which is used in the lattice Boltzmann method to construct approximate solutions of the incompressible Navier-Stokes equation.Michael Junk; Wen-An Yongpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1273Mon, 10 Dec 2001 00:00:00 +0100A Hybrid Simulation Method for Radivative Transfer Equations
https://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1274
We consider heat transfer processes where radiation in a large number of frequency bands plays a dominant role.Michael Junk; A. Unterreiter; F. Zingsheimpreprinthttps://kluedo.ub.uni-kl.de/frontdoor/index/index/docId/1274Mon, 10 Dec 2001 00:00:00 +0100