matthias weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/weigelt2012d.pdf · matthias...
TRANSCRIPT
![Page 1: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/1.jpg)
Workshop on Regional Gravity and Geomagnetic Field Modelling Bavarian Academy of Sciences and Humanities, Munich
February 23-24, 2012
Isoparametric Boundary Elements
Matthias Weigelt
Geodätisches Institut Stuttgart [email protected]
![Page 2: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/2.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Primary objective: local gravity field determination from CHAMP, GRACE (and GOCE)
Interest: improved static and time-variable gravity field
(make use of residual signal and data distribution)
Input data: potential values along the orbit
gravity gradient along the orbit
relative gradient along the orbit
Data type: currently only single data type
However: no restrictions on the combination of different data types (combination on normal equation level)
Application
Matthias Weigelt 2
![Page 3: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/3.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Decomposition of the surface into elements with finite extend (boundary elements)
Basic principle
Matthias Weigelt 3
e.g. blocks
e.g. triangles
![Page 4: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/4.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Example: Consider the single layer potential
Separation of the surface into elements:
Assumption:
The boundary elements are called isoparametric if
Basic principle
Matthias Weigelt 4
![Page 5: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/5.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Basic principle
Matthias Weigelt 5
Example: Consider the single layer potential
Including transformation to the normal triangle and spherical integration
Integration by Gaussian quadrature
Integration is exact for a polynomial of order 2L.
![Page 6: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/6.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Base functions are
– strictly space-limited (i.e. band-unlimited)
– compact
– continuous but not differentiable at the edges
– singular if the point of interest lies inside or on the corner of the element
• Weak singularity for potential
• Strong singularity for gradient
Mathematical properties
Matthias Weigelt 6
![Page 7: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/7.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Example: linear triangle
Base function
Matthias Weigelt 7
![Page 8: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/8.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Other base functions
Matthias Weigelt 8
six node triangle
four node quadrilateral
eight node quadrilateral: nine node quadrilateral: infinite quadrilateral:
![Page 9: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/9.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Point grid: any (as long as a proper tessellation is possible)
Currently based on the maxima and minima of an a priori field
Point Grids
Matthias Weigelt 9
![Page 10: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/10.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
• Search for maxima and minima of an a priori field
• Triangulation by Delaunay tessellation
• Least-squares adjustment
– brute-force
– assembly of the normal matrix (singularity!)
– full consideration of the stochastic information
• No regularization
– objective: avoid regularization by proper grid
– iterative search for vertices
Solution approach
Matthias Weigelt 10
![Page 11: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/11.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Closed loop simulation: noisefree and h=0km
Simulation study
Matthias Weigelt 11
![Page 12: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/12.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Closed loop simulation: noise=1% and h=400km
Simulation study
Matthias Weigelt 12
![Page 13: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/13.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Simulation study
Matthias Weigelt 13
![Page 14: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/14.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
• Optimal placement of the surface masses
– considering other quantities, e.g. curvature
– considering quadratic base elements
• Optimal tessellation using also combinations of elements
• Edge effect
– RCR technique
– Infinite elements ?
• Regularization
– avoidable by choosing the proper grid (?)
Open questions
Matthias Weigelt 14
![Page 15: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/15.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
ICCT study group JSG 0.6
Matthias Weigelt 15
Applicability of current GRACE solution strategies to the next generation of inter-satellite range
observations
Chairs: Matthias Weigelt, Co-Chair Adrian Jäggi
In collaboration with JSG 0.3
Comparison of methodologies for regional gravity field modelling
Chairs: M. Schmidt, Co-Chair: Ch. Gerlach
![Page 16: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/16.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Background
Matthias Weigelt 16
![Page 17: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/17.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
Approaches
17
Variational equations In-situ observations
Classical (Reigber 1989, Tapley 2004)
Celestial mechanics approach (Beutler et al. 2010, Jäggi 2007)
Short-arc method (Mayer-Gürr 2006)
…
Energy Integral (Han 2003, Ramillien et al. 2010)
Differential gravimetry (Liu 2010)
LoS Gradiometry (Keller and Sharifi 2005)
…
![Page 18: Matthias Weigelt - orbilu.uni.luorbilu.uni.lu/bitstream/10993/14345/1/Weigelt2012d.pdf · Matthias Weigelt Geodätisches ... The boundary elements are called isoparametric if Basic](https://reader031.vdocument.in/reader031/viewer/2022021900/5b6df0f97f8b9aed178df098/html5/thumbnails/18.jpg)
RegGrav Workshop Feb. 23-24, 2012, Munich
The objectives of the study group are therefore to:
• investigate each solution strategy, identify approximations and linearizations and test them for their permissibility to the next generation of inter-satellite range observations,
• identify limitations or the necessity for additional and/or more accurate measurements,
• quantify the sensitivity to error sources, e.g. in tidal or non-gravitational force modeling,
• investigate the interaction with global and local modeling,
• extend the applicability to planetary satellite mission, e.g. GRAIL
• establish a platform for the discussion and in-depth understanding of each approach and provide documentation.
Objectives
Matthias Weigelt 18