Spelling suggestions: "subject:"partielle differentialgleichung."" "subject:"partielle differentialgleichungs.""
71 |
Identification of material parameters in mechanical modelsMeyer, Marcus 04 June 2010 (has links)
Die Dissertation beschäftigt sich mit
Parameteridentifikationsproblemen, wie sie häufig in
Fragestellungen der Festkörpermechanik zu finden sind. Hierbei
betrachten wir die Identifikation von Materialparametern -- die
typischerweise die Eigenschaften der zugrundeliegenden
Materialien repräsentieren -- aus gemessenen Verformungen oder
Belastungen eines Testkörpers. In mathematischem Sinne
entspricht dies der Lösung von Identifikationsproblemen, die
eine spezielle Klasse von inversen Problemen bilden.
Der Inhalt der Dissertation ist folgendermaßen gegliedert. Nach
dem einführenden Abschnitt 1 wird in Abschnitt 2 ein Überblick
von Optimierungs- und Regularisierungsverfahren zur stabilen
Lösung nichtlinearer inverser Probleme diskutiert. In Abschnitt
3 betrachten wir die Identifikation von skalaren und stückweise
konstanten Parametern in linearen elliptischen
Differentialgleichungen. Hierbei werden zwei Testprobleme
erörtert, die Identifikation von Diffusions- und
Reaktionsparameter in einer allgemeinen elliptischen
Differentialgleichung und die Identifikation der
Lame-Konstanten in einem Modell der linearisierten Elastizität.
Die zugrunde liegenden PDE-Modelle und Lösungszugänge werden
erläutert. Insbesondere betrachten wir hier Newton-artige
Algorithmen, Gradientenmethoden, Multi-Parameter
Regularisierung and den evolutionären Algorithmus CMAES.
Abschließend werden Ergebnisse einer numerischen Studie
präsentiert. Im Abschnitt 4 konzentrieren wir uns auf die
Identifikation von verteilten Parametern in hyperelastischen
Materialmodellen. Das nichtlineare Elastizitätsproblem wird
detailiert erläutert und verschiedene Materialmodelle werden
diskutiert (linear elastisches St.-Venant-Kirchhoff Material
und nichtlineare Neo-Hooke, Mooney-Rivlin und Modified-Fung
Materialien. Zur Lösung des resultierenden
Parameteridentifikationsproblems werden Lösungsansätze aus der
optimalen Steuerung in Form eines Newton-Lagrange SQP
Algorithmus verwendet. Die Resultate einer numerischen Studie
werden präsentiert, basierend auf einem zweidimensionales
Testproblem mit einer sogenannten Cook-Mebran. Abschließend
wird im Abschnitt 5 die Verwendung adaptiver FEM für die Lösung
von Parameteridentifikationsproblems kurz erörtert. / The dissertation is focussed on parameter identification
problems arising in the context of structural mechanics. At
this, we consider the identification of material parameters -
which typically represent the properties of an underlying
material - from given measured displacements and forces of a
loaded test body. In mathematical terms such problems denote
identification problems as a special case of general inverse
problems.
The dissertation is organized as follows. After the
introductive section 1, section 2 is devoted to a survey of
optimization and regularization methods for the stable solution
of nonlinear inverse problems. In section 3 we consider the
identification of scalar and piecewise constant parameters in
linear elliptic differential equations and examine two test
problems, namely the identification of diffusion and reaction
parameters in a generalized linear elliptic differential
equation of second order and the identification of the Lame
constants in the linearized elasticity model. The underlying
PDE models are introduced and solution approaches are discussed
in detail. At this, we consider Newton-type algorithms,
gradient methods, multi-parameter regularization, and the
evolutionary algorithm CMAES. Consequently, numerical studies
for a two-dimensional test problem are presented. In section 4
we point out the identification of distributed material
parameters in hyperelastic deformation models. The nonlinear
elasticity boundary value problem for large deformations is
introduced. We discuss several material laws for linear elastic
(St.-Venant-Kirchhoff) materials and nonlinear Neo-Hooke,
Mooney-Rivlin, and Modified-Fung materials. For the solution of
the corresponding parameter identification problem, we focus on
an optimal control solution approach and introduce a
regularized Newton-Lagrange SQP method. The Newton-Lagrange
algorithm is demonstrated within a numerical study. Therefore,
a simplified two-dimensional Cook membrane test problem is
solved. Additionally, in section 5 the application of adaptive
methods for the solution of parameter identification problems
is discussed briefly.
|
72 |
Conditional stability estimates for ill-posed PDE problems by using interpolationTautenhahn, Ulrich, Hämarik, Uno, Hofmann, Bernd, Shao, Yuanyuan January 2011 (has links)
The focus of this paper is on conditional stability estimates for ill-posed inverse problems in partial differential equations. Conditional stability estimates have been obtained in the literature by a couple different methods. In this paper we propose a method called interpolation method, which is based on interpolation in variable Hilbert scales. We are going to work out the theoretical background of this method and show that optimal conditional stability estimates are obtained. The capability of our method is illustrated by a comprehensive collection of different inverse and ill-posed PDE problems containing elliptic and parabolic problems, one source problem and the problem of analytic continuation.
|
73 |
Solving optimal PDE control problems : optimality conditions, algorithms and model reductionPrüfert, Uwe 16 May 2016 (has links)
This thesis deals with the optimal control of PDEs. After a brief introduction in the theory of elliptic and parabolic PDEs, we introduce a software that solves systems of PDEs by the finite elements method. In the second chapter we derive optimality conditions in terms of function spaces, i.e. a systems of PDEs coupled by some pointwise relations. Now we present algorithms to solve the optimality systems numerically and present some numerical test cases. A further chapter deals with the so called lack of adjointness, an issue of gradient methods applied on parabolic optimal control problems. However, since optimal control problems lead to large numerical schemes, model reduction becomes popular. We analyze the proper orthogonal decomposition method and apply it to our model problems. Finally, we apply all considered techniques to a real world problem.:Introduction
The state equation
Optimal control and optimality conditions
Algorithms
The \"lack of adjointness\"
Numerical examples
Efficient solution of PDEs and KKT- systems
A real world application
Functional analytical basics
Codes of the examples
|
74 |
On the quasi-optimal convergence of adaptive nonconforming finite element methods in three examplesRabus, Hella 23 May 2014 (has links)
Eine Vielzahl von Anwendungen in der numerischen Simulation der Strömungsdynamik und der Festkörpermechanik begründen die Entwicklung von zuverlässigen und effizienten Algorithmen für nicht-standard Methoden der Finite-Elemente-Methode (FEM). Um Freiheitsgrade zu sparen, wird in jedem Durchlauf des adaptiven Algorithmus lediglich ein Teil der Gebiete verfeinert. Einige Gebiete bleiben daher möglicherweise verhältnismäßig grob. Die Analyse der Konvergenz und vor allem die der Optimalität benötigt daher über die a priori Fehleranalyse hinausgehende Argumente. Etablierte adaptive Algorithmen beruhen auf collective marking, d.h. die zu verfeinernden Gebiete werden auf Basis eines Gesamtfehlerschätzers markiert. Bei adaptiven Algorithmen mit separate marking wird der Gesamtfehlerschätzer in einen Volumenterm und in einen Fehlerschätzerterm aufgespalten. Da der Volumenterm unabhängig von der diskreten Lösung ist, kann einer schlechten Datenapproximation durch eine lokal tiefe Verfeinerung begegnet werden. Bei hinreichender Datenapproximation wird das Gitter dagegen bezüglich des neuen Fehlerschätzerterms wie üblich level-orientiert verfeinert. Die numerischen Experimente dieser Arbeit liefern deutliche Indizien der quasi-optimalen Konvergenz für den in dieser Arbeit untersuchten adaptiven Algorithmus, der auf separate marking beruht. Der Parameter, der die Verbesserung der Datenapproximation sicherstellt, ist frei wählbar. Dadurch ist es erstmals möglich, eine ausreichende und gleichzeitig optimale Approximation der Daten innerhalb weniger Durchläufe zu erzwingen. Diese Arbeit ermöglicht es, Standardargumente auch für die Konvergenzanalyse von Algorithmen mit separate marking zu verwenden. Dadurch gelingt es Quasi-Optimalität des vorgestellten Algorithmus gemäß einer generellen Vorgehensweise für die drei Beispiele, dem Poisson Modellproblem, dem reinen Verschiebungsproblem der linearen Elastizität und dem Stokes Problem, zu zeigen. / Various applications in computational fluid dynamics and solid mechanics motivate the development of reliable and efficient adaptive algorithms for nonstandard finite element methods (FEMs). To reduce the number of degrees of freedom, in adaptive algorithms only a selection of finite element domains is marked for refinement on each level. Since some element domains may stay relatively coarse, even the analysis of convergence and more importantly the analysis of optimality require new arguments beyond an a priori error analysis. In adaptive algorithms, based on collective marking, a (total) error estimator is used as refinement indicator. For separate marking strategies, the (total) error estimator is split into a volume term and an error estimator term, which estimates the error. Since the volume term is independent of the discrete solution, if there is a poor data approximation the improvement may be realised by a possibly high degree of local mesh refinement. Otherwise, a standard level-oriented mesh refinement based on an error estimator term is performed. This observation results in a natural adaptive algorithm based on separate marking, which is analysed in this thesis. The results of the numerical experiments displayed in this thesis provide strong evidence for the quasi-optimality of the presented adaptive algorithm based on separate marking and for all three model problems. Furthermore its flexibility (in particular the free steering parameter for data approximation) allows a sufficient and optimal data approximation in just a few number of levels of the adaptive scheme. This thesis adapts standard arguments for optimal convergence to adaptive algorithms based on separate marking with a possibly high degree of local mesh refinement, and proves quasi-optimality following a general methodology for three model problems, i.e., the Poisson model problem, the pure displacement problem in linear elasticity and the Stokes equations.
|
75 |
Structural damped sigma-evolution operatorsKainane Mezadek, Mohamed 05 March 2014 (has links)
The subject of the thesis is the investigation of asymptotic properties of solutions of the Cauchy problem for structurally damped sigma-evolution operators with time dependent, monotonous, dissipation term.
An appropriate energy for solutions of the sigma-evolution equations is defined and some estimates for energies of higher order are proved. In the scale invariant case the optimality of these estimates is shown.
Further, the influence of properties of the time dependent dissipation on L^p-L^q estimates for the energy with p and q bigger or equal to 2 and from the conjugate line is clarified.
Also smoothing properties of the operators under consideration are investigated. The connection between the regularity of the data and the regularity of the solution in terms of L^2 based Gevrey spaces is considered.
Finally, L^1-L^1-estimates in the special case delta = sigma/2 and decreasing dissipative coefficient. / Thema der vorliegenden Dissertation ist die Untersuchung asymptotischer Eigenschaften von Lösungen des Cauchy Problems für strukturell gedämpfte sigma-Evolutions-Operatoren mit zeitabhängigem, monotonen Dissipationskoeffizienten.
Es wird eine geeignete Energie definiert und für diese Abschätzungen, auf für entsprechende Energien höherer Ordnung gezeigt. Darüber hinaus wird der Einfluss des Dissipationskoeffizienten auf L^p-L^q Abschätzungen auf und entfernt von der konjugierten Linie untersucht. Im skaleninvarianten Fall wird die Schärfe der Abschätzungen bewiesen.
Weiterhin wird der Zusammenhang zwischen der Regularität der Daten und der der Lösung in Termen von L^2-basierten Gevrey-Räumen untersucht.
Schließlich werden L^1-L^1-Abschätzungen für den Spezialfall delta = sigma/2 und monoton fallenden Dissipationskoeffizienten gezeigt.
|
76 |
Stabilised finite element approximation for degenerate convex minimisation problemsBoiger, Wolfgang Josef 19 August 2013 (has links)
Infimalfolgen nichtkonvexer Variationsprobleme haben aufgrund feiner Oszillationen häufig keinen starken Grenzwert in Sobolevräumen. Diese Oszillationen haben eine physikalische Bedeutung; Finite-Element-Approximationen können sie jedoch im Allgemeinen nicht auflösen. Relaxationsmethoden ersetzen die nichtkonvexe Energie durch ihre (semi)konvexe Hülle. Das entstehende makroskopische Modell ist degeneriert: es ist nicht strikt konvex und hat eventuell mehrere Minimalstellen. Die fehlende Kontrolle der primalen Variablen führt zu Schwierigkeiten bei der a priori und a posteriori Fehlerschätzung, wie der Zuverlässigkeits- Effizienz-Lücke und fehlender starker Konvergenz. Zur Überwindung dieser Schwierigkeiten erweitern Stabilisierungstechniken die relaxierte Energie um einen diskreten, positiv definiten Term. Bartels et al. (IFB, 2004) wenden Stabilisierung auf zweidimensionale Probleme an und beweisen dabei starke Konvergenz der Gradienten. Dieses Ergebnis ist auf glatte Lösungen und quasi-uniforme Netze beschränkt, was adaptive Netzverfeinerungen ausschließt. Die vorliegende Arbeit behandelt einen modifizierten Stabilisierungsterm und beweist auf unstrukturierten Netzen sowohl Konvergenz der Spannungstensoren, als auch starke Konvergenz der Gradienten für glatte Lösungen. Ferner wird der sogenannte Fluss-Fehlerschätzer hergeleitet und dessen Zuverlässigkeit und Effizienz gezeigt. Für Interface-Probleme mit stückweise glatter Lösung wird eine Verfeinerung des Fehlerschätzers entwickelt, die den Fehler der primalen Variablen und ihres Gradienten beschränkt und so starke Konvergenz der Gradienten sichert. Der verfeinerte Fehlerschätzer konvergiert schneller als der Fluss- Fehlerschätzer, und verringert so die Zuverlässigkeits-Effizienz-Lücke. Numerische Experimente mit fünf Benchmark-Tests der Mikrostruktursimulation und Topologieoptimierung ergänzen und bestätigen die theoretischen Ergebnisse. / Infimising sequences of nonconvex variational problems often do not converge strongly in Sobolev spaces due to fine oscillations. These oscillations are physically meaningful; finite element approximations, however, fail to resolve them in general. Relaxation methods replace the nonconvex energy with its (semi)convex hull. This leads to a macroscopic model which is degenerate in the sense that it is not strictly convex and possibly admits multiple minimisers. The lack of control on the primal variable leads to difficulties in the a priori and a posteriori finite element error analysis, such as the reliability-efficiency gap and no strong convergence. To overcome these difficulties, stabilisation techniques add a discrete positive definite term to the relaxed energy. Bartels et al. (IFB, 2004) apply stabilisation to two-dimensional problems and thereby prove strong convergence of gradients. This result is restricted to smooth solutions and quasi-uniform meshes, which prohibit adaptive mesh refinements. This thesis concerns a modified stabilisation term and proves convergence of the stress and, for smooth solutions, strong convergence of gradients, even on unstructured meshes. Furthermore, the thesis derives the so-called flux error estimator and proves its reliability and efficiency. For interface problems with piecewise smooth solutions, a refined version of this error estimator is developed, which provides control of the error of the primal variable and its gradient and thus yields strong convergence of gradients. The refined error estimator converges faster than the flux error estimator and therefore narrows the reliability-efficiency gap. Numerical experiments with five benchmark examples from computational microstructure and topology optimisation complement and confirm the theoretical results.
|
77 |
Wachstumsanalyse amorpher dicker Schichten und Schichtsysteme / Growth analysis of thick amorphous films and multilayersStreng, Christoph 18 May 2004 (has links)
No description available.
|
78 |
Wachstumsanalyse amorpher dicker Schichten und Schichtsysteme / Growth analysis of thick amorphous films and multilayersStreng, Christoph 18 May 2004 (has links)
No description available.
|
79 |
Beiträge zur Regularisierung inverser Probleme und zur bedingten Stabilität bei partiellen DifferentialgleichungenShao, Yuanyuan 17 January 2013 (has links) (PDF)
Wir betrachten die lineare inverse Probleme mit gestörter rechter Seite und gestörtem Operator in Hilberträumen, die inkorrekt sind. Um die Auswirkung der Inkorrektheit zu verringen, müssen spezielle Lösungsmethode angewendet werden, hier nutzen wir die sogenannte Tikhonov Regularisierungsmethode. Die Regularisierungsparameter wählen wir aus das verallgemeinerte Defektprinzip. Eine typische numerische Methode zur Lösen der nichtlinearen äquivalenten Defektgleichung ist Newtonverfahren. Wir schreiben einen Algorithmus, die global und monoton konvergent für beliebige Startwerte garantiert.
Um die Stabilität zu garantieren, benutzen wir die Glattheit der Lösung, dann erhalten wir eine sogenannte bedingte Stabilität. Wir demonstrieren die sogenannte Interpolationsmethode zur Herleitung von bedingten Stabilitätsabschätzungen bei inversen Problemen für partielle Differentialgleichungen.
|
80 |
Beiträge zur Regularisierung inverser Probleme und zur bedingten Stabilität bei partiellen DifferentialgleichungenShao, Yuanyuan 14 January 2013 (has links)
Wir betrachten die lineare inverse Probleme mit gestörter rechter Seite und gestörtem Operator in Hilberträumen, die inkorrekt sind. Um die Auswirkung der Inkorrektheit zu verringen, müssen spezielle Lösungsmethode angewendet werden, hier nutzen wir die sogenannte Tikhonov Regularisierungsmethode. Die Regularisierungsparameter wählen wir aus das verallgemeinerte Defektprinzip. Eine typische numerische Methode zur Lösen der nichtlinearen äquivalenten Defektgleichung ist Newtonverfahren. Wir schreiben einen Algorithmus, die global und monoton konvergent für beliebige Startwerte garantiert.
Um die Stabilität zu garantieren, benutzen wir die Glattheit der Lösung, dann erhalten wir eine sogenannte bedingte Stabilität. Wir demonstrieren die sogenannte Interpolationsmethode zur Herleitung von bedingten Stabilitätsabschätzungen bei inversen Problemen für partielle Differentialgleichungen.
|
Page generated in 0.4953 seconds