View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0003686 | OpenFOAM | Bug | public | 2021-06-11 15:30 | 2021-06-11 16:11 |
| Reporter | fabio1068 | Assigned To | henry | ||
| Priority | normal | Severity | minor | Reproducibility | have not tried |
| Status | closed | Resolution | suspended | ||
| Platform | Windows 10 WSL | OS | Ubuntu | OS Version | 18.04 LTS |
| Summary | 0003686: SnappyHexMesh castellation leads to poor snapped mesh quality | ||||
| Description | Sometimes the castellation process leaves some "baffle boundary faces" that prevent the following smoothing and snapping process to produce good quality mesh. In attachment you find some Paraview screenshots of the mesh after castellation, smoothing and snapping steps, with the troubling faces highlighted in red. I used the planingHullW3 tutorial with a slightly modified SnappyHexMeshDict (more refined feature edges) to reproduce the issue. I think those faces should not belong to the "hull" boundary, but be normal internal faces instead. | ||||
| Steps To Reproduce | - Copy the OpenFOAM-8/tutorials/multiphase/interFoam/RAS/planingHullW3/ tutorial folder - Substitute the snappyHexMeshDict.2 dictionary with the one in attachment - Run Allmesh.2 - Run ParaFoam to see the mesh (look at the bottom edge of the transom, in the aft part of the boat). If the bug is reproduced, it should look like the ones in attachment | ||||
| Tags | No tags attached. | ||||
|
|
log.snappyHexMesh (115,374 bytes)
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 8-1c9b5879390b
Exec : snappyHexMesh -dict system/snappyHexMeshDict.2 -overwrite
Date : Jun 11 2021
Time : 15:29:10
Host : "DESKTOP-EI0UH3G"
PID : 1981
I/O : uncollated
Case : /home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Read mesh in = 0.03 s
Overall mesh bounding box : (-14.16 -4.8 -6.6) (6 0 3)
Relative tolerance : 1e-06
Absolute matching distance : 2.28391e-05
Reading refinement surfaces.
Read refinement surfaces in = 0.26 s
Reading refinement shells.
Refinement level according to distance to hull
level 3 for all cells within 0.05 metre.
level 2 for all cells within 0.1 metre.
level 2 for all cells within 0.2 metre.
Refinement level 2 for all cells inside wakeRegion
Read refinement shells in = 0 s
Setting refinement level of surface to be consistent with shells.
Checked shell refinement in = 0 s
Reading features.
Read edgeMesh w3.eMesh
points : 1611
edges : 1618
boundingBox : (-0.001479 -0.212 -0.034244) (1.89517 0.212 0.277918)
Refinement level according to distance to "w3.eMesh" (1611 points, 1618 edges).
level 4 for all cells within 0.02 metre.
Read features in = 0 s
Determining initial surface intersections
-----------------------------------------
Edge intersection testing:
Number of edges : 93051
Number of edges to retest : 93051
Number of intersected edges : 372
Calculated surface intersections in = 0.28 s
Initial mesh : cells(local):29397 faces(local):93051 points(local):34430
Initial mesh : cells:29397 faces:93051 points:34430
Cells per refinement level:
0 29397
Adding patches for surface regions
----------------------------------
Patch Type Region
----- ---- ------
hull:
6 wall hull
Added patches in = 0.02 s
Selecting decompositionMethod none
Layer thickness specified as final layer and expansion ratio.
Refinement phase
----------------
Found point (-0.7 0 0) in cell 24184 on processor 0
Feature refinement iteration 0
------------------------------
Marked for refinement due to explicit features : 78 cells.
Determined cells to refine in = 0.15 s
Selected for feature refinement : 78 cells (out of 29397)
getChangedFaces : Detected local:3383 global:3383 changed faces out of 95112
getChangedFaces : Writing 3383 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 95112
Number of edges to retest : 3383
Number of intersected edges : 721
Writing refined but unbalanced "feature refinement iteration 0" mesh to time 0.0001
Dumped debug data in = 0.09 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.02 s
After refinement feature refinement iteration 0 : cells(local):29943 faces(local):95112 points(local):35462
After refinement feature refinement iteration 0 : cells:29943 faces:95112 points:35462
Cells per refinement level:
0 29319
1 624
Feature refinement iteration 1
------------------------------
Marked for refinement due to explicit features : 175 cells.
Determined cells to refine in = 0.15 s
Selected for feature refinement : 311 cells (out of 29943)
getChangedFaces : Detected local:11636 global:11636 changed faces out of 102969
getChangedFaces : Writing 11636 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 102969
Number of edges to retest : 11636
Number of intersected edges : 1218
Writing refined but unbalanced "feature refinement iteration 1" mesh to time 0.0002
Dumped debug data in = 0.08 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.03 s
After refinement feature refinement iteration 1 : cells(local):32120 faces(local):102969 points(local):39109
After refinement feature refinement iteration 1 : cells:32120 faces:102969 points:39109
Cells per refinement level:
0 29183
1 1537
2 1400
Feature refinement iteration 2
------------------------------
Marked for refinement due to explicit features : 351 cells.
Determined cells to refine in = 0.17 s
Selected for feature refinement : 759 cells (out of 32120)
getChangedFaces : Detected local:29265 global:29265 changed faces out of 121902
getChangedFaces : Writing 29265 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 121902
Number of edges to retest : 29265
Number of intersected edges : 3056
Writing refined but unbalanced "feature refinement iteration 2" mesh to time 0.0003
Dumped debug data in = 0.11 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.02 s
After refinement feature refinement iteration 2 : cells(local):37433 faces(local):121902 points(local):47724
After refinement feature refinement iteration 2 : cells:37433 faces:121902 points:47724
Cells per refinement level:
0 29064
1 2200
2 3361
3 2808
Feature refinement iteration 3
------------------------------
Marked for refinement due to explicit features : 699 cells.
Determined cells to refine in = 0.28 s
Selected for feature refinement : 1582 cells (out of 37433)
getChangedFaces : Detected local:62828 global:62828 changed faces out of 161244
getChangedFaces : Writing 62828 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 161244
Number of edges to retest : 62828
Number of intersected edges : 8399
Writing refined but unbalanced "feature refinement iteration 3" mesh to time 0.0004
Dumped debug data in = 0.17 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.06 s
After refinement feature refinement iteration 3 : cells(local):48507 faces(local):161244 points(local):65590
After refinement feature refinement iteration 3 : cells:48507 faces:161244 points:65590
Cells per refinement level:
0 29009
1 2434
2 4387
3 7085
4 5592
Feature refinement iteration 4
------------------------------
Marked for refinement due to explicit features : 0 cells.
Determined cells to refine in = 0.34 s
Selected for feature refinement : 0 cells (out of 48507)
Stopping refining since too few cells selected.
Surface refinement iteration 0
------------------------------
Marked for refinement due to surface intersection : 2297 cells.
Determined cells to refine in = 0.02 s
Selected for refinement : 2410 cells (out of 48507)
getChangedFaces : Detected local:90437 global:90437 changed faces out of 215157
getChangedFaces : Writing 90437 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 215157
Number of edges to retest : 90437
Number of intersected edges : 12227
Writing refined but unbalanced "surface refinement iteration 0" mesh to time 0.0005
Dumped debug data in = 0.21 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.09 s
After refinement surface refinement iteration 0 : cells(local):65377 faces(local):215157 points(local):85655
After refinement surface refinement iteration 0 : cells:65377 faces:215157 points:85655
Cells per refinement level:
0 28723
1 4034
2 8455
3 18573
4 5592
Surface refinement iteration 1
------------------------------
Marked for refinement due to surface intersection : 3429 cells.
Determined cells to refine in = 0.01 s
Selected for refinement : 3656 cells (out of 65377)
getChangedFaces : Detected local:124978 global:124978 changed faces out of 297624
getChangedFaces : Writing 124978 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 297624
Number of edges to retest : 124978
Number of intersected edges : 18474
Writing refined but unbalanced "surface refinement iteration 1" mesh to time 0.0006
Dumped debug data in = 0.31 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.13 s
After refinement surface refinement iteration 1 : cells(local):90969 faces(local):297624 points(local):116847
After refinement surface refinement iteration 1 : cells:90969 faces:297624 points:116847
Cells per refinement level:
0 28654
1 3488
2 14750
3 38485
4 5592
Surface refinement iteration 2
------------------------------
Marked for refinement due to surface intersection : 3726 cells.
Determined cells to refine in = 0.02 s
Selected for refinement : 4066 cells (out of 90969)
getChangedFaces : Detected local:130974 global:130974 changed faces out of 390465
getChangedFaces : Writing 130974 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 390465
Number of edges to retest : 130974
Number of intersected edges : 25745
Writing refined but unbalanced "surface refinement iteration 2" mesh to time 0.0007
Dumped debug data in = 0.38 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.16 s
After refinement surface refinement iteration 2 : cells(local):119431 faces(local):390465 points(local):152744
After refinement surface refinement iteration 2 : cells:119431 faces:390465 points:152744
Cells per refinement level:
0 28616
1 3489
2 13449
3 68285
4 5592
Surface refinement iteration 3
------------------------------
Marked for refinement due to surface intersection : 6 cells.
Determined cells to refine in = 0.01 s
Selected for refinement : 7 cells (out of 119431)
getChangedFaces : Detected local:701 global:701 changed faces out of 390618
getChangedFaces : Writing 701 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 390618
Number of edges to retest : 701
Number of intersected edges : 25745
Writing refined "surface refinement iteration 3" mesh to time 0.0008
Dumped debug data in = 0.3 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.13 s
After refinement surface refinement iteration 3 : cells(local):119480 faces(local):390618 points(local):152799
After refinement surface refinement iteration 3 : cells:119480 faces:390618 points:152799
Cells per refinement level:
0 28616
1 3488
2 13451
3 68333
4 5592
Surface refinement iteration 4
------------------------------
Marked for refinement due to surface intersection : 0 cells.
Determined cells to refine in = 0.03 s
Selected for refinement : 0 cells (out of 119480)
Stopping refining since too few cells selected.
Removing mesh beyond surface intersections
------------------------------------------
Found point (-0.7 0 0) in global region 0 out of 3 regions.
Keeping all cells in region 0 containing point (-0.7 0 0)
Selected for keeping : 101795 cells.
Edge intersection testing:
Number of edges : 339571
Number of edges to retest : 19708
Number of intersected edges : 25745
Writing subsetted mesh to time 0.0009.
Dumped mesh in = 1.8 s
Shell refinement iteration 0
----------------------------
Marked for refinement due to distance to explicit features : 10855 cells.
Marked for refinement due to refinement shells : 12936 cells.
Dumping 23791 cells to cellSet candidateCellsFromShells.
Determined cells to refine in = 0.2 s
Selected for internal refinement : 25229 cells (out of 101795)
getChangedFaces : Detected local:746839 global:746839 changed faces out of 894646
getChangedFaces : Writing 746839 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 894646
Number of edges to retest : 746839
Number of intersected edges : 40550
Writing refined "shell refinement iteration 0" mesh to time 0.001
Dumped debug data in = 1.37 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.35 s
After refinement shell refinement iteration 0 : cells(local):278398 faces(local):894646 points(local):339139
After refinement shell refinement iteration 0 : cells:278398 faces:894646 points:339139
Cells per refinement level:
0 26532
1 16632
2 21643
3 111645
4 101946
Shell refinement iteration 1
----------------------------
Marked for refinement due to distance to explicit features : 3471 cells.
Marked for refinement due to refinement shells : 19443 cells.
Dumping 22914 cells to cellSet candidateCellsFromShells.
Determined cells to refine in = 0.41 s
Selected for internal refinement : 23677 cells (out of 278398)
getChangedFaces : Detected local:800699 global:800699 changed faces out of 1403689
getChangedFaces : Writing 800699 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 1403689
Number of edges to retest : 800699
Number of intersected edges : 40550
Writing refined "shell refinement iteration 1" mesh to time 0.0011
Dumped debug data in = 1.52 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.49 s
After refinement shell refinement iteration 1 : cells(local):444137 faces(local):1403689 points(local):516775
After refinement shell refinement iteration 1 : cells:444137 faces:1403689 points:516775
Cells per refinement level:
0 25736
1 8578
2 132031
3 148078
4 129714
Shell refinement iteration 2
----------------------------
Marked for refinement due to distance to explicit features : 2 cells.
Marked for refinement due to refinement shells : 775 cells.
Dumping 777 cells to cellSet candidateCellsFromShells.
Determined cells to refine in = 0.28 s
Selected for internal refinement : 1283 cells (out of 444137)
getChangedFaces : Detected local:80259 global:80259 changed faces out of 1430911
getChangedFaces : Writing 80259 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 1430911
Number of edges to retest : 80259
Number of intersected edges : 40550
Writing refined "shell refinement iteration 2" mesh to time 0.0012
Dumped debug data in = 1.26 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.48 s
After refinement shell refinement iteration 2 : cells(local):453118 faces(local):1430911 points(local):526036
After refinement shell refinement iteration 2 : cells:453118 faces:1430911 points:526036
Cells per refinement level:
0 25230
1 11865
2 138105
3 148188
4 129730
Shell refinement iteration 3
----------------------------
Marked for refinement due to distance to explicit features : 0 cells.
Marked for refinement due to refinement shells : 0 cells.
Dumping 0 cells to cellSet candidateCellsFromShells.
Determined cells to refine in = 0.28 s
Selected for internal refinement : 86 cells (out of 453118)
getChangedFaces : Detected local:6186 global:6186 changed faces out of 1432774
getChangedFaces : Writing 6186 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 1432774
Number of edges to retest : 6186
Number of intersected edges : 40550
Writing refined "shell refinement iteration 3" mesh to time 0.0013
Dumped debug data in = 1.18 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.47 s
After refinement shell refinement iteration 3 : cells(local):453720 faces(local):1432774 points(local):526694
After refinement shell refinement iteration 3 : cells:453720 faces:1432774 points:526694
Cells per refinement level:
0 25144
1 12553
2 138105
3 148188
4 129730
Shell refinement iteration 4
----------------------------
Marked for refinement due to distance to explicit features : 0 cells.
Marked for refinement due to refinement shells : 0 cells.
Dumping 0 cells to cellSet candidateCellsFromShells.
Determined cells to refine in = 0.28 s
Selected for internal refinement : 0 cells (out of 453720)
Stopping refining since too few cells selected.
Dangling coarse cells refinement iteration 0
--------------------------------------------
Dumping 9 cells to cellSet candidateCellSet.
Determined cells to refine in = 0.02 s
Selected for refinement : 9 cells (out of 453720)
getChangedFaces : Detected local:1176 global:1176 changed faces out of 1432909
getChangedFaces : Writing 1176 changed faces to faceSet changedFaces
Edge intersection testing:
Number of edges : 1432909
Number of edges to retest : 1176
Number of intersected edges : 40550
Writing refined "coarse cell refinement iteration 0" mesh to time 0.0014
Dumped debug data in = 1.2 s
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh in = 0.49 s
After refinement coarse cell refinement iteration 0 : cells(local):453783 faces(local):1432909 points(local):526712
After refinement coarse cell refinement iteration 0 : cells:453783 faces:1432909 points:526712
Cells per refinement level:
0 25137
1 12609
2 138105
3 148186
4 129746
Dangling coarse cells refinement iteration 1
--------------------------------------------
Dumping 0 cells to cellSet candidateCellSet.
Determined cells to refine in = 0.02 s
Selected for refinement : 0 cells (out of 453783)
Stopping refining since too few cells selected.
Dangling coarse cells refinement iteration 0
--------------------------------------------
Dumping 0 cells to cellSet candidateCellSet.
Determined cells to refine in = 0.01 s
Selected for refinement : 0 cells (out of 453783)
Stopping refining since too few cells selected.
Splitting mesh at surface intersections
---------------------------------------
Introducing baffles for 40550 faces that are intersected by the surface.
Edge intersection testing:
Number of edges : 1473346
Number of edges to retest : 261130
Number of intersected edges : 77306
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 40437 duplicate pairs of faces.
Created baffles in = 1.82 s
After introducing baffles : cells(local):453783 faces(local):1473346 points(local):526712
After introducing baffles : cells:453783 faces:1473346 points:526712
Cells per refinement level:
0 25137
1 12609
2 138105
3 148186
4 129746
Writing baffled mesh to time 0.0014
Dumped debug data in = 0.03 s
Introducing baffles to block off problem cells
----------------------------------------------
markFacesOnProblemCells : marked 36876 additional internal faces to be converted into baffles.
Analyzed problem cells in = 0.73 s
Dumping 36876 problem faces to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/0.0014/polyMesh/sets/problemFaces"
Introducing baffles to delete problem cells.
Edge intersection testing:
Number of edges : 1510222
Number of edges to retest : 129471
Number of intersected edges : 77331
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 77313 duplicate pairs of faces.
Created baffles in = 1.81 s
After introducing baffles : cells(local):453783 faces(local):1510222 points(local):526712
After introducing baffles : cells:453783 faces:1510222 points:526712
Cells per refinement level:
0 25137
1 12609
2 138105
3 148186
4 129746
Writing extra baffled mesh to time 0.0015
Dumped debug data in = 0.06 s
Remove unreachable sections of mesh
-----------------------------------
Keeping all cells in region 0 containing point (-0.7 0 0)
Selected for keeping : 397667 cells.
Edge intersection testing:
Number of edges : 1264645
Number of edges to retest : 0
Number of intersected edges : 39819
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 3 duplicate pairs of faces.
Split mesh in = 7.77 s
After subsetting : cells(local):397667 faces(local):1264645 points(local):470739
After subsetting : cells:397667 faces:1264645 points:470739
Cells per refinement level:
0 25137
1 12609
2 138105
3 129066
4 92750
Writing subsetted mesh to time 0.0016
Dumped debug data in = 0.03 s
Handling cells with snap problems
---------------------------------
Introducing baffles for 39819 faces that are intersected by the surface.
Edge intersection testing:
Number of edges : 1264645
Number of edges to retest : 130361
Number of intersected edges : 39819
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 3 duplicate pairs of faces.
Created baffles in = 1.06 s
After introducing baffles : cells(local):397667 faces(local):1264645 points(local):470739
After introducing baffles : cells:397667 faces:1264645 points:470739
Cells per refinement level:
0 25137
1 12609
2 138105
3 129066
4 92750
Writing baffled mesh to time 0.0017
Dumped debug data in = 0.05 s
Introducing baffles to block off problem cells
----------------------------------------------
markFacesOnProblemCells : marked 0 additional internal faces to be converted into baffles.
Analyzed problem cells in = 0.61 s
Dumping 0 problem faces to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/0.0017/polyMesh/sets/problemFaces"
Introducing baffles to delete problem cells.
Edge intersection testing:
Number of edges : 1264645
Number of edges to retest : 0
Number of intersected edges : 39819
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 3 duplicate pairs of faces.
Created baffles in = 0.9 s
After introducing baffles : cells(local):397667 faces(local):1264645 points(local):470739
After introducing baffles : cells:397667 faces:1264645 points:470739
Cells per refinement level:
0 25137
1 12609
2 138105
3 129066
4 92750
Writing extra baffled mesh to time 0.0018
Dumped debug data in = 0.03 s
Remove unreachable sections of mesh
-----------------------------------
Keeping all cells in region 0 containing point (-0.7 0 0)
Selected for keeping : 397667 cells.
Edge intersection testing:
Number of edges : 1264645
Number of edges to retest : 0
Number of intersected edges : 39819
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 3 duplicate pairs of faces.
Split mesh in = 6.74 s
After subsetting : cells(local):397667 faces(local):1264645 points(local):470739
After subsetting : cells:397667 faces:1264645 points:470739
Cells per refinement level:
0 25137
1 12609
2 138105
3 129066
4 92750
Writing subsetted mesh to time 0.0019
Dumped debug data in = 0.03 s
Merge free-standing baffles
---------------------------
freeStandingBaffles : detected 3 free-standing baffles out of 3
freeStandingBaffles : detected 0 planar (within 45 degrees) free-standing baffles out of 3
Detected free-standing baffles : 0
Merged free-standing baffles in = 0.19 s
dupNonManifoldPoints : Found : 5 non-manifold points (out of 470739)
Edge intersection testing:
Number of edges : 1264645
Number of edges to retest : 0
Number of intersected edges : 39819
Detected unsplittable baffles : 0
Writing handleProblemCells mesh to time 0.002.
Merge refined boundary faces
----------------------------
Merging 1252 sets of faces.
Writing all faces to be merged to set "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/0.002/polyMesh/sets/allFaceSets"
Edge intersection testing:
Number of edges : 1262203
Number of edges to retest : 6331
Number of intersected edges : 37381
Checking sync after initial merging 1252 faces.
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Writing initial merged-faces mesh to time 0.0021
Undo iteration 0
----------------
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Merging all points on surface that
- are used by only two boundary faces and
- make an angle with a cosine of more than 0.707107.
Removing 1774 straight edge points ...
Edge intersection testing:
Number of edges : 1262203
Number of edges to retest : 8418
Number of intersected edges : 37381
Checking sync after removing points.
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Undo iteration 0
----------------
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing merged-edges mesh to time 0.0022
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Refined mesh : cells(local):397667 faces(local):1262203 points(local):468375
Refined mesh : cells:397667 faces:1262203 points:468375
Cells per refinement level:
0 25137
1 12609
2 138105
3 129066
4 92750
Writing mesh to time 0.0022
Wrote mesh in = 5.53 s.
Mesh refined in = 42.97 s.
Morphing phase
--------------
Snapping to features in 10 iterations ...
Constructing mesh displacer ...
Using mesh parameters
{
maxNonOrtho 70;
maxBoundarySkewness 20;
maxInternalSkewness 4;
maxConcave 80;
minVol 1e-13;
minTetQuality 1e-30;
minArea -1;
minTwist 0.05;
minDeterminant 0.001;
minFaceWeight 0.05;
minVolRatio 0.01;
minTriangleTwist -1;
nSmoothScale 4;
errorReduction 0.75;
relaxed
{
maxNonOrtho 75;
}
}
Checking initial mesh ...
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Detected 0 illegal faces (concave, zero area or negative cell pyramid volume)
Checked initial mesh in = 1.14 s
Smoothing patch points ...
Smoothing iteration 0
Found 10 non-manifold point(s).
Scaling iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Smoothing iteration 1
Found 0 non-manifold point(s).
Scaling iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Smoothing iteration 2
Found 0 non-manifold point(s).
Scaling iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Scaling iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Scaling iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Writing patch smoothed mesh to time 0.0023.
Dumped mesh in = 5.04 s
Patch points smoothed in = 0 s
Morph iteration 0
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.00109263 min:1.56963e-08 max:0.00394387
Calculated surface displacement in = 0.05 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 102 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 2 baffle edges out of 77416 edges.
Initially selected 2417 points out of 39336 for reverse attraction.
Selected 6460 points out of 39336 for reverse attraction.
Stringing feature edges : changed 262 points
Stringing feature edges : changed 30 points
Stringing feature edges : changed 11 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000239052 -3.07874e-07 -0.00393662) avg:(-1.85565e-05 0.000341564 8.61804e-05)
feature : max:(-0.00286619 0.000362817 0.00259219) avg:(-3.1642e-06 -2.33217e-06 6.62357e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 7
feature edge : 2136
nearest surface : 0
rest : 37193
--> FOAM Warning : Displacement (-2.66255e-05 -4.86169e-10 -4.20969e-06) at mesh point 83059 coord (0.00599616 -0.12 0.0668062) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.48 s
Writing smoothed mesh to time 0.0024
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0024.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 1.82 s
Writing scaled mesh to time 0.0025
Writing displacement field ...
Morph iteration 1
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.00101598 min:1.41602e-07 max:0.00367559
Calculated surface displacement in = 0.11 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 109 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 2 baffle edges out of 77416 edges.
Initially selected 2413 points out of 39336 for reverse attraction.
Selected 6449 points out of 39336 for reverse attraction.
Stringing feature edges : changed 253 points
Stringing feature edges : changed 33 points
Stringing feature edges : changed 10 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000224046 -4.39094e-07 -0.00366875) avg:(-1.69558e-05 0.000312372 7.9133e-05)
feature : max:(-0.00283752 0.000326535 0.00256627) avg:(-2.92714e-06 -3.53888e-06 6.26745e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 7
feature edge : 2133
nearest surface : 0
rest : 37196
--> FOAM Warning : Displacement (-8.62988e-05 -8.55577e-10 -1.43939e-05) at mesh point 83059 coord (0.00596954 -0.12 0.066802) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.49 s
Writing smoothed mesh to time 0.0026
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0026.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0027
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0028
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0029
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 3.29 s
Writing scaled mesh to time 0.003
Writing displacement field ...
Morph iteration 2
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.000872485 min:2.98076e-09 max:0.00320711
Calculated surface displacement in = 0.11 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 129 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 2 baffle edges out of 77416 edges.
Initially selected 2413 points out of 39336 for reverse attraction.
Selected 6440 points out of 39336 for reverse attraction.
Stringing feature edges : changed 239 points
Stringing feature edges : changed 30 points
Stringing feature edges : changed 10 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000195474 -3.97682e-07 -0.00320115) avg:(-1.40545e-05 0.000258126 6.59166e-05)
feature : max:(-0.00272402 0.000261228 0.00246362) avg:(-2.34135e-06 -4.81376e-06 6.20116e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 7
feature edge : 2127
nearest surface : 0
rest : 37202
--> FOAM Warning : Displacement (-0.000173645 1.1936e-09 -2.60342e-05) at mesh point 83059 coord (0.00588324 -0.12 0.0667876) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.45 s
Writing smoothed mesh to time 0.0031
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0031.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0032
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0033
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0034
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0035
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0036
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0037
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 4.89 s
Writing scaled mesh to time 0.0038
Writing displacement field ...
Morph iteration 3
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.000688141 min:9.72468e-09 max:0.00263082
Calculated surface displacement in = 0.14 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 136 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 0 baffle edges out of 77416 edges.
Initially selected 2410 points out of 39336 for reverse attraction.
Selected 6431 points out of 39336 for reverse attraction.
Stringing feature edges : changed 228 points
Stringing feature edges : changed 30 points
Stringing feature edges : changed 9 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000160349 -3.2623e-07 -0.00262593) avg:(-1.03739e-05 0.00019065 4.8942e-05)
feature : max:(-0.00247886 0.00018286 0.00224189) avg:(-1.86573e-06 -6.12886e-06 5.87365e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 6
feature edge : 2126
nearest surface : 0
rest : 37204
--> FOAM Warning : Displacement (-0.000264126 5.09645e-09 -3.27942e-05) at mesh point 83059 coord (0.00570959 -0.12 0.0667616) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.48 s
Writing smoothed mesh to time 0.0039
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0039.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.004
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0041
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0042
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0043
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0044
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0045
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 4.81 s
Writing scaled mesh to time 0.0046
Writing displacement field ...
Morph iteration 4
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.000493952 min:2.81185e-08 max:0.00206054
Calculated surface displacement in = 0.09 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 136 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 0 baffle edges out of 77416 edges.
Initially selected 2407 points out of 39336 for reverse attraction.
Selected 6425 points out of 39336 for reverse attraction.
Stringing feature edges : changed 225 points
Stringing feature edges : changed 26 points
Stringing feature edges : changed 6 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000125145 3.05467e-07 -0.00205673) avg:(-6.76954e-06 0.000123538 3.18964e-05)
feature : max:(-0.00208224 0.000109716 0.00188319) avg:(-1.22483e-06 -6.65149e-06 5.32994e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 6
feature edge : 2118
nearest surface : 0
rest : 37212
--> FOAM Warning : Displacement (1.04754e-05 -5.09215e-06 1.79805e-06) at mesh point 102244 coord (0.603113 -0.0260289 0.0116902) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.61 s
Writing smoothed mesh to time 0.0047
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0047.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 11
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0048
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0049
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.005
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0051
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0052
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0053
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 4.86 s
Writing scaled mesh to time 0.0054
Writing displacement field ...
Morph iteration 5
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.000316551 min:4.75065e-08 max:0.0019965
Calculated surface displacement in = 0.13 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 156 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 0 baffle edges out of 77416 edges.
Initially selected 2405 points out of 39336 for reverse attraction.
Selected 6404 points out of 39336 for reverse attraction.
Stringing feature edges : changed 205 points
Stringing feature edges : changed 14 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000121442 9.11976e-08 -0.00199281) avg:(-3.91726e-06 6.85216e-05 1.77387e-05)
feature : max:(-0.00156168 5.48579e-05 0.00141239) avg:(-7.62892e-07 -5.76265e-06 4.12508e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 6
feature edge : 2100
nearest surface : 0
rest : 37230
--> FOAM Warning : Displacement (1.56616e-05 -3.12643e-05 -6.53766e-05) at mesh point 102244 coord (0.603123 -0.026034 0.011692) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.55 s
Writing smoothed mesh to time 0.0055
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0055.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 3
faces with face-decomposition tet quality < 1e-30 : 52
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 55
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0056
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 12
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 79
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0057
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 6
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 55
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0058
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 3
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 40
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0059
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 35
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.006
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 2
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 31
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0061
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 5.79 s
Writing scaled mesh to time 0.0062
Writing displacement field ...
Morph iteration 6
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:0.000176279 min:2.82842e-10 max:0.0019965
Calculated surface displacement in = 0.11 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 169 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 0 baffle edges out of 77416 edges.
Initially selected 2405 points out of 39336 for reverse attraction.
Selected 6404 points out of 39336 for reverse attraction.
Stringing feature edges : changed 171 points
Stringing feature edges : changed 14 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000121442 9.11976e-08 -0.00199281) avg:(-2.03501e-06 3.21853e-05 7.89214e-06)
feature : max:(0.000155321 -0.0011382 0.000840045) avg:(-4.10847e-07 -4.59459e-06 2.96676e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 6
feature edge : 2100
nearest surface : 0
rest : 37230
--> FOAM Warning : Displacement (1.38599e-05 -4.00253e-05 -9.3821e-05) at mesh point 102244 coord (0.603139 -0.0260652 0.0116266) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.66 s
Writing smoothed mesh to time 0.0063
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0063.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 44
faces with face-decomposition tet quality < 1e-30 : 87
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 47
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0064
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 19
faces with face-decomposition tet quality < 1e-30 : 50
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 78
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0065
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 10
faces with face-decomposition tet quality < 1e-30 : 40
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 72
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0066
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 1
faces with face-decomposition tet quality < 1e-30 : 34
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 65
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0067
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 29
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 66
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0068
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 23
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 63
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0069
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 6.03 s
Writing scaled mesh to time 0.007
Writing displacement field ...
Morph iteration 7
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:8.23038e-05 min:8.88692e-09 max:0.0019965
Calculated surface displacement in = 0.08 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 195 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 0 baffle edges out of 77416 edges.
Initially selected 2403 points out of 39336 for reverse attraction.
Selected 6400 points out of 39336 for reverse attraction.
Stringing feature edges : changed 155 points
Stringing feature edges : changed 14 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000121442 9.11976e-08 -0.00199281) avg:(-9.80132e-07 1.21115e-05 2.75879e-06)
feature : max:(0.000711162 -0.00127841 0.00013439) avg:(-1.59025e-07 -3.34252e-06 1.13769e-06)
Feature analysis : total master points:39336 attraction to :
feature point : 6
feature edge : 2099
nearest surface : 0
rest : 37231
--> FOAM Warning : Displacement (8.16724e-06 -3.12009e-05 -7.75327e-05) at mesh point 102244 coord (0.603153 -0.0261053 0.0115328) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.48 s
Writing smoothed mesh to time 0.0071
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0071.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 83
faces with face-decomposition tet quality < 1e-30 : 192
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 30
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0072
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 49
faces with face-decomposition tet quality < 1e-30 : 124
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 52
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0073
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 19
faces with face-decomposition tet quality < 1e-30 : 75
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 92
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0074
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 8
faces with face-decomposition tet quality < 1e-30 : 58
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 97
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0075
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 1
faces with face-decomposition tet quality < 1e-30 : 54
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 81
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0076
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 49
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 79
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0077
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 6.36 s
Writing scaled mesh to time 0.0078
Writing displacement field ...
Morph iteration 8
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:3.1902e-05 min:5.88172e-10 max:0.0019965
Calculated surface displacement in = 0.11 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 203 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 7 baffle edges out of 77416 edges.
Initially selected 2405 points out of 39336 for reverse attraction.
Selected 6404 points out of 39336 for reverse attraction.
Stringing feature edges : changed 154 points
Stringing feature edges : changed 12 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000121442 9.11976e-08 -0.00199281) avg:(-5.24131e-07 3.46549e-06 1.13773e-06)
feature : max:(0.000168305 0.000236246 -0.00122216) avg:(-4.16546e-09 -2.34688e-06 -1.74544e-07)
Feature analysis : total master points:39336 attraction to :
feature point : 6
feature edge : 2101
nearest surface : 0
rest : 37229
--> FOAM Warning : Displacement (3.38212e-06 -1.59437e-05 -4.09579e-05) at mesh point 102244 coord (0.603161 -0.0261357 0.0114571) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.46 s
Writing smoothed mesh to time 0.0079
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0079.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 107
faces with face-decomposition tet quality < 1e-30 : 299
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 34
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.008
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 86
faces with face-decomposition tet quality < 1e-30 : 214
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 46
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0081
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 44
faces with face-decomposition tet quality < 1e-30 : 181
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 72
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0082
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 19
faces with face-decomposition tet quality < 1e-30 : 129
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 117
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0083
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 6
faces with face-decomposition tet quality < 1e-30 : 107
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 114
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0084
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 89
faces with concavity > 80 degrees : 2
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 98
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0085
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 6.86 s
Writing scaled mesh to time 0.0086
Writing displacement field ...
Morph iteration 9
-----------------
Calculating patchDisplacement as distance to nearest surface point ...
Wanted displacement : average:1.15328e-05 min:1.10413e-11 max:0.0019965
Calculated surface displacement in = 0.1 s
Detecting near surfaces ...
Overriding nearest with intersection of close gaps at 168 out of 39336 points.
Overriding displacement on features :
implicit features : false
explicit features : true
multi-patch features : false
Detected 27 baffle edges out of 77416 edges.
Initially selected 2403 points out of 39336 for reverse attraction.
Selected 6402 points out of 39336 for reverse attraction.
Stringing feature edges : changed 147 points
Stringing feature edges : changed 12 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 5 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 3 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 4 points
Stringing feature edges : changed 1 points
Stringing feature edges : changed 0 points
Attraction:
linear : max:(0.000121442 9.11976e-08 -0.00199281) avg:(-2.86983e-07 1.0217e-06 9.21301e-07)
feature : max:(0.000168305 0.000236246 -0.00122216) avg:(2.85559e-08 -1.89343e-06 -6.47618e-07)
Feature analysis : total master points:39336 attraction to :
feature point : 6
feature edge : 2099
nearest surface : 0
rest : 37231
--> FOAM Warning : Displacement (9.27938e-07 -5.11847e-06 -1.34066e-05) at mesh point 102244 coord (0.603164 -0.0261504 0.0114195) points through the surrounding patch faces
Smoothing displacement ...
Iteration 0
Iteration 10
Iteration 20
Iteration 30
Iteration 40
Iteration 50
Iteration 60
Iteration 70
Iteration 80
Iteration 90
Displacement smoothed in = 2.53 s
Writing smoothed mesh to time 0.0087
Writing displacement field ...
Writing actual patch displacement ...
Dumping move direction to "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/actualPatchDisplacement_0.0087.obj"
Moving mesh ...
Iteration 0
Moving mesh using displacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 165
faces with face-decomposition tet quality < 1e-30 : 315
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 13
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0088
Writing displacement field ...
Iteration 1
Moving mesh using displacement scaling : min:0.75 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 123
faces with face-decomposition tet quality < 1e-30 : 296
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 30
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0089
Writing displacement field ...
Iteration 2
Moving mesh using displacement scaling : min:0.5625 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 88
faces with face-decomposition tet quality < 1e-30 : 238
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 53
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.009
Writing displacement field ...
Iteration 3
Moving mesh using displacement scaling : min:0.421875 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 51
faces with face-decomposition tet quality < 1e-30 : 209
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 73
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0091
Writing displacement field ...
Iteration 4
Moving mesh using displacement scaling : min:0.316406 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 20
faces with face-decomposition tet quality < 1e-30 : 165
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 133
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0092
Writing displacement field ...
Iteration 5
Displacement scaling for error reduction set to 0.
Moving mesh using displacement scaling : min:0.237305 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 1
faces with face pyramid volume < 1e-13 : 7
faces with face-decomposition tet quality < 1e-30 : 135
faces with concavity > 80 degrees : 1
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 133
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0093
Writing displacement field ...
Iteration 6
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 1
faces on cells with determinant < 0.001 : 0
Writing scaled mesh to time 0.0094
Writing displacement field ...
Iteration 7
Moving mesh using displacement scaling : min:0 max:1
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Successfully moved mesh
Moved mesh in = 7.9 s
Writing scaled mesh to time 0.0095
Writing displacement field ...
Repatching faces according to nearest surface ...
Repatched 0 faces in = 0.23 s
Edge intersection testing:
Number of edges : 1262203
Number of edges to retest : 190297
Number of intersected edges : 38081
Merging 6916 sets of faces.
Writing all faces to be merged to set "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/0.002/polyMesh/sets/allFaceSets"
Edge intersection testing:
Number of edges : 1253472
Number of edges to retest : 29375
Number of intersected edges : 29350
Checking sync after initial merging 6916 faces.
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Writing initial merged-faces mesh to time 0.0096
Undo iteration 0
----------------
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 9
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 3
faces on cells with determinant < 0.001 : 0
Writing all faces in error to faceSet "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/0.0096/polyMesh/sets/errorFaces"
Masters that need to be restored:16
Writing all 16 masterfaces to be restored to set "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/0.0096/polyMesh/sets/mastersToRestore"
Edge intersection testing:
Number of edges : 1253493
Number of edges to retest : 106
Number of intersected edges : 29371
Checking sync after restoring 37 faces.
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Writing merged-faces mesh to time 0.0097
Undo iteration 1
----------------
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Merging all points on surface that
- are used by only two boundary faces and
- make an angle with a cosine of more than 0.707107.
Removing 5126 straight edge points ...
Edge intersection testing:
Number of edges : 1253493
Number of edges to retest : 27685
Number of intersected edges : 29371
Checking sync after removing points.
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Undo iteration 0
----------------
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 1
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
**Writing all faces in error to faceSet "/home/fabio/OpenFOAM/fabio-8/run/planingHullW3_original/0.0097/polyMesh/sets/errorFaces"
Detected 0 error faces on boundaries that have been merged. These will be restored to their original faces.
Detected 1 error faces in mesh. Restoring neighbours of faces in error.
Edge intersection testing:
Number of edges : 1253493
Number of edges to retest : 45
Number of intersected edges : 29371
Checking sync after restoring points on 10 faces.
meshRefinement::checkData() : Checking refinement structure.
meshRefinement::checkData() : Checking refinement levels.
meshRefinement::checkData() : Checking synchronization.
meshRefinement::checkData() : Counting duplicate faces.
meshRefinement::checkData() : Found 0 duplicate pairs of faces.
Writing merged-edges mesh to time 0.0098
Writing patchFace merged mesh to time 0.0099
Snapped mesh : cells(local):397667 faces(local):1253493 points(local):461449
Snapped mesh : cells:397667 faces:1253493 points:461449
Cells per refinement level:
0 25137
1 12609
2 138105
3 129066
4 92750
Writing mesh to time 0.0099
Wrote mesh in = 8.25 s.
Mesh snapped in = 93.49 s.
Checking final mesh ...
Checking faces in error :
non-orthogonality > 70 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with face-decomposition tet quality < 1e-30 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.05 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.05 : 0
faces on cells with determinant < 0.001 : 0
Finished meshing without any errors
Finished meshing in = 137.98 s.
End
snappyHexMeshDict.2 (8,637 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object snappyHexMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Which of the steps to run
castellatedMesh true;
snap true;
addLayers false;
// Geometry. Definition of all surfaces. All surfaces are of class
// searchableSurface.
// Surfaces are used
// - to specify refinement for any mesh cell intersecting it
// - to specify refinement for any mesh cell inside/outside/near
// - to 'snap' the mesh boundary to the surface
geometry
{
hull
{
type triSurfaceMesh;
file "w3.stl";
patchInfo
{
type wall;
}
}
wakeRegion
{
type searchableBox;
min (-1.5 -0.4 -0.15);
max (0.1 0.0 0.15);
}
};
// Settings for the castellatedMesh generation.
castellatedMeshControls
{
// Refinement parameters
// ~~~~~~~~~~~~~~~~~~~~~
// If local number of cells is >= maxLocalCells on any processor
// switches from from refinement followed by balancing
// (current method) to (weighted) balancing before refinement.
maxLocalCells 100000;
// Overall cell limit (approximately). Refinement will stop immediately
// upon reaching this number so a refinement level might not complete.
// Note that this is the number of cells before removing the part which
// is not 'visible' from the keepPoint. The final number of cells might
// actually be a lot less.
maxGlobalCells 20000000;
// The surface refinement loop might spend lots of iterations refining just a
// few cells. This setting will cause refinement to stop if <= minimumRefine
// are selected for refinement. Note: it will at least do one iteration
// (unless the number of cells to refine is 0)
minRefinementCells 0;
// Number of buffer layers between different levels.
// 1 means normal 2:1 refinement restriction, larger means slower
// refinement.
nCellsBetweenLevels 3;
// Explicit feature edge refinement
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
// Specifies a level for any cell intersected by its edges.
// This is a featureEdgeMesh, read from constant/triSurface for now.
features
(
{
file "w3.eMesh";
levels ((0.02 4));
}
);
// Surface based refinement
// ~~~~~~~~~~~~~~~~~~~~~~~~
// Specifies two levels for every surface. The first is the minimum level,
// every cell intersecting a surface gets refined up to the minimum level.
// The second level is the maximum level. Cells that 'see' multiple
// intersections where the intersections make an
// angle > resolveFeatureAngle get refined up to the maximum level.
refinementSurfaces
{
hull
{
// Surface-wise min and max refinement level
level (3 3);
}
}
resolveFeatureAngle 45;
// Region-wise refinement
// ~~~~~~~~~~~~~~~~~~~~~~
// Specifies refinement level for cells in relation to a surface. One of
// three modes
// - distance. 'levels' specifies per distance to the surface the
// wanted refinement level. The distances need to be specified in
// descending order.
// - inside. 'levels' is only one entry and only the level is used. All
// cells inside the surface get refined up to the level. The surface
// needs to be closed for this to be possible.
// - outside. Same but cells outside.
refinementRegions
{
hull
{
mode distance;
levels ((0.05 3) (0.1 2) (0.2 2));
}
wakeRegion
{
mode inside;
levels ((1.0 2));
}
}
// Mesh selection
// ~~~~~~~~~~~~~~
// After refinement patches get added for all refinementSurfaces and
// all cells intersecting the surfaces get put into these patches. The
// section reachable from the locationInMesh is kept.
// NOTE: This point should never be on a face, always inside a cell, even
// after refinement.
locationInMesh (-0.7 0 0);
// Whether any faceZones (as specified in the refinementSurfaces)
// are only on the boundary of corresponding cellZones or also allow
// free-standing zone faces. Not used if there are no faceZones.
allowFreeStandingZoneFaces true;
}
// Settings for the snapping.
snapControls
{
//- Number of patch smoothing iterations before finding correspondence
// to surface
nSmoothPatch 3;
//- Relative distance for points to be attracted by surface feature point
// or edge. True distance is this factor times local
// maximum edge length.
// tolerance 4.0;
tolerance 1.0;
//- Number of mesh displacement relaxation iterations.
nSolveIter 100;
//- Maximum number of snapping relaxation iterations. Should stop
// before upon reaching a correct mesh.
nRelaxIter 5;
nFeatureSnapIter 10;
}
// Settings for the layer addition.
addLayersControls
{
// Are the thickness parameters below relative to the undistorted
// size of the refined cell outside layer (true) or absolute sizes (false).
relativeSizes true;
// Per final patch (so not geometry!) the layer information
layers
{
hull
{
nSurfaceLayers 3;
}
}
// Expansion factor for layer mesh
expansionRatio 1.5;
// Wanted thickness of final added cell layer. If multiple layers
// is the thickness of the layer furthest away from the wall.
// Relative to undistorted size of cell outside layer.
// See relativeSizes parameter.
finalLayerThickness 0.7;
// Minimum thickness of cell layer. If for any reason layer
// cannot be above minThickness do not add layer.
// See relativeSizes parameter.
minThickness 0.25;
// If points get not extruded do nGrow layers of connected faces that are
// also not grown. This helps convergence of the layer addition process
// close to features.
// Note: changed(corrected) w.r.t 17x! (didn't do anything in 17x)
nGrow 0;
// Advanced settings
// When not to extrude surface. 0 is flat surface, 90 is when two faces
// are perpendicular
featureAngle 60;
// Maximum number of snapping relaxation iterations. Should stop
// before upon reaching a correct mesh.
nRelaxIter 5;
// Number of smoothing iterations of surface normals
nSmoothSurfaceNormals 1;
// Number of smoothing iterations of interior mesh movement direction
nSmoothNormals 3;
// Smooth layer thickness over surface patches
nSmoothThickness 10;
// Stop layer growth on highly warped cells
maxFaceThicknessRatio 0.5;
// Reduce layer growth where ratio thickness to medial
// distance is large
maxThicknessToMedialRatio 0.3;
// Angle used to pick up medial axis points
// Note: changed(corrected) w.r.t 17x! 90 degrees corresponds to 130 in 17x.
minMedianAxisAngle 90;
// Create buffer region for new layer terminations
nBufferCellsNoExtrude 0;
// Overall max number of layer addition iterations. The mesher will exit
// if it reaches this number of iterations; possibly with an illegal
// mesh.
nLayerIter 50;
// Max number of iterations after which relaxed meshQuality controls
// get used. Up to nRelaxIter it uses the settings in meshQualityControls,
// after nRelaxIter it uses the values in meshQualityControls::relaxed.
nRelaxedIter 20;
}
// Generic mesh quality settings. At any undoable phase these determine
// where to undo.
meshQualityControls
{
#include "meshQualityDict"
}
// Advanced
// Flags for optional output
// 0 : only write final meshes
// 1 : write intermediate meshes
// 2 : write volScalarField with cellLevel for postprocessing
// 4 : write current intersections as .obj files
debug 1;
// Merge tolerance. Is fraction of overall bounding box of initial mesh.
// Note: the write tolerance needs to be higher than this.
mergeTolerance 1E-6;
// ************************************************************************* //
|
|
|
Can you provide a patch to snappyHexMesh to resolve this problem or funding for a core OpenFOAM developer to work on it for you? |
|
|
I think I am not skilled enough in programming to fix this bug and at the moment I am not able to fund a developer to do it (maybe one day I will be). But I thought it was good to report it anyway. |
|
|
I would expect this issue to be very difficult and time-consuming to resolve, someone will have to pay for the work. |
|
|
Pending funding |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2021-06-11 15:30 | fabio1068 | New Issue | |
| 2021-06-11 15:30 | fabio1068 | File Added: 99_snapped_.png | |
| 2021-06-11 15:30 | fabio1068 | File Added: log.snappyHexMesh | |
| 2021-06-11 15:30 | fabio1068 | File Added: snappyHexMeshDict.2 | |
| 2021-06-11 15:30 | fabio1068 | File Added: 22_castellated_.png | |
| 2021-06-11 15:30 | fabio1068 | File Added: 23_smoothed_.png | |
| 2021-06-11 15:46 | henry | Note Added: 0012063 | |
| 2021-06-11 16:00 | fabio1068 | Note Added: 0012064 | |
| 2021-06-11 16:11 | henry | Note Added: 0012065 | |
| 2021-06-11 16:11 | henry | Assigned To | => henry |
| 2021-06-11 16:11 | henry | Status | new => closed |
| 2021-06-11 16:11 | henry | Resolution | open => suspended |
| 2021-06-11 16:11 | henry | Note Added: 0012066 |
