Element statistics

General discussion about Elmer
Post Reply
Rich_B
Posts: 421
Joined: 24 Aug 2009, 20:18

Element statistics

Post by Rich_B »

Hello,

As a follow up to the EM_HT_LE thread, it seemed that the Element Statistics solver listed in ElmerModelsManual might be useful to check the quality of the mesh.

The attached test case uses a simple angle3d.grd file to create a mesh, and then running elmersolver leads to a seg fault.

Here is the end of the solver log:
SetActiveElementsTable: Creating active element table for: elementstats
SetActiveElementsTable: Number of active elements found : 300
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: elementstats
ElementStats: ------------------------------
ElementStats: Studying elements of the mesh
ElementStats: Going through bulk elements
Rich.
Attachments
quality.zip
(1.45 KiB) Downloaded 47 times
Rich_B
Posts: 421
Joined: 24 Aug 2009, 20:18

Re: Element statistics

Post by Rich_B »

Hello,

After a small edit of ElementStats.F90, see attached, here are the statistics for the EM_HT_LE spiral model, also attached:
ElementStats: ------------------------------
ElementStats: Studying elements of the mesh
ElementStats: Going through bulk elements
Statistics for mesh operator: Element size on bulk
Minimum values Element size: 0.55999881409878982
Maximum values Element size: 4269.8599285751388
Average values Element size: 209.91834385011725
Ratio of element sizes: 7624.7660192756939

Histogram:
1 : 24399 ( 0.55999881409878982 to 427.48999179020279 )
2 : 1923 ( 427.48999179020279 to 854.41998476630681 )
3 : 1428 ( 854.41998476630681 to 1281.3499777424108 )
4 : 530 ( 1281.3499777424108 to 1708.2799707185150 )
5 : 206 ( 1708.2799707185150 to 2135.2099636946191 )
6 : 32 ( 2135.2099636946191 to 2562.1399566707228 )
7 : 6 ( 2562.1399566707228 to 2989.0699496468269 )
8 : 5 ( 2989.0699496468269 to 3415.9999426229310 )
9 : 10 ( 3415.9999426229310 to 3842.9299355990352 )
10 : 37 ( 3842.9299355990352 to 4269.8599285751388 )
Statistics for mesh operator: Element skew (degs) on bulk
Minimum values Element skew (degs): 0.0000000000000000
Maximum values Element skew (degs): 0.0000000000000000
Average values Element skew (degs): 0.0000000000000000

Histogram:
1 : 28576 ( 0.0000000000000000 to 0.0000000000000000 )
Statistics for mesh operator: Element ratio on bulk
Minimum values Element ratio: 1.0219296326372667
Maximum values Element ratio: 10.187464658321593
Average values Element ratio: 1.4517417017050855

Histogram:
1 : 27767 ( 1.0219296326372667 to 1.9384831352056993 )
2 : 455 ( 1.9384831352056993 to 2.8550366377741319 )
3 : 5 ( 2.8550366377741319 to 3.7715901403425645 )
4 : 8 ( 3.7715901403425645 to 4.6881436429109975 )
5 : 16 ( 4.6881436429109975 to 5.6046971454794292 )
6 : 17 ( 5.6046971454794292 to 6.5212506480478627 )
7 : 155 ( 6.5212506480478627 to 7.4378041506162944 )
8 : 11 ( 7.4378041506162944 to 8.3543576531847279 )
9 : 132 ( 8.3543576531847279 to 9.2709111557531614 )
10 : 10 ( 9.2709111557531614 to 10.187464658321593 )
ElementStats: Going through boundary elements
Statistics for mesh operator: Element size on boundary
Minimum values Element size: 1.5375627800450000
Maximum values Element size: 544.44444444755550
Average values Element size: 60.315310373087030
Ratio of element sizes: 354.09574913853027

Histogram:
1 : 30670 ( 1.5375627800450000 to 55.828250946796047 )
2 : 2095 ( 55.828250946796047 to 110.11893911354709 )
3 : 1482 ( 110.11893911354709 to 164.40962728029814 )
4 : 749 ( 164.40962728029814 to 218.70031544704918 )
5 : 775 ( 218.70031544704918 to 272.99100361380022 )
6 : 120 ( 272.99100361380022 to 327.28169178055128 )
7 : 7 ( 327.28169178055128 to 381.57237994730230 )
8 : 5 ( 381.57237994730230 to 435.86306811405336 )
9 : 10 ( 435.86306811405336 to 490.15375628080443 )
10 : 280 ( 490.15375628080443 to 544.44444444755550 )
Statistics for mesh operator: Element skew (degs) on boundary
Minimum values Element skew (degs): 0.0000000000000000
Maximum values Element skew (degs): 0.0000000000000000
Average values Element skew (degs): 0.0000000000000000

Histogram:
1 : 36193 ( 0.0000000000000000 to 0.0000000000000000 )
Statistics for mesh operator: Element ratio on boundary
Minimum values Element ratio: 1.0000000000000000
Maximum values Element ratio: 8.8191710368945575
Average values Element ratio: 1.2124762585015294

Histogram:
1 : 35310 ( 1.0000000000000000 to 1.7819171036894557 )
2 : 457 ( 1.7819171036894557 to 2.5638342073789113 )
3 : 5 ( 2.5638342073789113 to 3.3457513110683674 )
4 : 8 ( 3.3457513110683674 to 4.1276684147578226 )
5 : 16 ( 4.1276684147578226 to 4.9095855184472788 )
6 : 17 ( 4.9095855184472788 to 5.6915026221367349 )
7 : 155 ( 5.6915026221367349 to 6.4734197258261901 )
8 : 47 ( 6.4734197258261901 to 7.2553368295156462 )
9 : 132 ( 7.2553368295156462 to 8.0372539332051023 )
10 : 46 ( 8.0372539332051023 to 8.8191710368945575 )
ElementStats: ------------------------------
Attachments
ElementStats.F90
(13.85 KiB) Downloaded 45 times
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Element statistics

Post by raback »

Hi Rich,

Thank you for the fix! I just applied it to the devel branch.

-Peter
Rich_B
Posts: 421
Joined: 24 Aug 2009, 20:18

Re: Element statistics

Post by Rich_B »

Hello Peter,

Glad to help!

Thanks, Rich.
Post Reply