|
[Sponsors] |
[dRInterfaceLib] dynamicInterfaceRefineFvMesh libary tobias holzmann |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 19, 2017, 13:07 |
dynamicInterfaceRefineFvMesh libary tobias holzmann
|
#1 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10 |
https://www.cfd-online.com/Forums/op...o-regions.html
here our super moderator Tobias developed a fantastic extension for mesh refinement at the interface. i am trying to compile for openfoam-v1612+ i got following error Code:
wmake libso . wmakeLnInclude: linking include files to ./lnInclude Making dependency list for source file dynamicInterfaceRefineFvMesh.C g++ -std=c++11 -m64 -DOPENFOAM_PLUS=1612 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3 -DNoRepository -ftemplate-depth-100 -I/home/hasan/OpenFOAM/OpenFOAM-v1612+/src/triSurface/lnInclude -I/home/hasan/OpenFOAM/OpenFOAM-v1612+/src/meshTools/lnInclude -I/home/hasan/OpenFOAM/OpenFOAM-v1612+/src/dynamicMesh/lnInclude -I/home/hasan/OpenFOAM/OpenFOAM-v1612+/src/finiteVolume/lnInclude -I/home/hasan/OpenFOAM/OpenFOAM-v1612+/src/dynamicFvMesh/lnInclude -IlnInclude -I. -I/home/hasan/OpenFOAM/OpenFOAM-v1612+/src/OpenFOAM/lnInclude -I/home/hasan/OpenFOAM/OpenFOAM-v1612+/src/OSspecific/POSIX/lnInclude -fPIC -c dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C -o Make/linux64GccDPInt32Opt/dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.o dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C: In member function ‘Foam::autoPtr<Foam::mapPolyMesh> Foam::dynamicInterfaceRefineFvMesh::refine(const labelList&)’: dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C:349:17: error: ‘fillSignallingNan’ is not a member of ‘Foam::sigFpe’ sigFpe::fillSignallingNan(phi.internalField()); ^ dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C:389:13: error: ‘GeometricBoundaryField’ is not a member of ‘Foam::surfaceScalarField {aka Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>}’ surfaceScalarField::GeometricBoundaryField& bphi = ^ dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C:389:57: error: ‘bphi’ was not declared in this scope surfaceScalarField::GeometricBoundaryField& bphi = ^ dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C: In member function ‘Foam::autoPtr<Foam::mapPolyMesh> Foam::dynamicInterfaceRefineFvMesh::unrefine(const labelList&)’: dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C:579:13: error: ‘GeometricBoundaryField’ is not a member of ‘Foam::surfaceScalarField {aka Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>}’ surfaceScalarField::GeometricBoundaryField& bphi = ^ dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C:579:57: error: ‘bphi’ was not declared in this scope surfaceScalarField::GeometricBoundaryField& bphi = ^ dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C: In member function ‘virtual bool Foam::dynamicInterfaceRefineFvMesh::writeObject(Foam::IOstream::streamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const’: dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.C:1662:9: error: ‘writeObjects’ is not a member of ‘Foam::dynamicFvMesh’ dynamicFvMesh::writeObjects(fmt, ver, cmp) ^ /home/hasan/OpenFOAM/OpenFOAM-v1612+/wmake/rules/General/transform:25: recipe for target 'Make/linux64GccDPInt32Opt/dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.o' failed make: *** [Make/linux64GccDPInt32Opt/dynamicInterfaceRefineFvMesh/dynamicInterfaceRefineFvMesh.o] Error 1 hasan@hasan-Lenovo-G50-70:~/OpenFOAM/hasan-v1612+/OpenFOAM_extensions/dynamicinterfacerefinefvmesh/src/dynamicFvMesh$ please take a look Last edited by Tobi; August 22, 2017 at 10:58. |
|
August 22, 2017, 11:01 |
|
#2 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52 |
Hi Saddy,
first of all thanks for the first sentence To your question. Unfortunately the library is based on OpenFOAM-2.3.x and I should update it to 5.0 (maybe today in the evening). However, if you just need an interface refinement, you can use the standard library. Just set the refinement between e.g. alpha 0.4 and 0.6. In addition to your question, the field names changes during different releases and even more if different openfoam versions are used. Mine is based on the Foundation one but you are using the one from ESI. I am not sure if it is compatible. You can try to compile the library if I updated it to version 5.x but I don't think that it will work out of the box. But again, an interface refinement can be done with the normal library too.
__________________
Keep foaming, Tobias Holzmann |
|
August 22, 2017, 15:30 |
|
#3 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52 |
Done,
I got the same errors in the 5.x version. So now you can use the library for the Foundation Version and I guess also for the ESI one (because I had the same errors). Cheers.
__________________
Keep foaming, Tobias Holzmann |
|
August 23, 2017, 05:30 |
|
#4 | |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10 |
Quote:
i have also found that openfoam lacks in refinement department. if there was anisotropic/unstructured/tetrahedral refinement available. openfoam would be very robust. thanks so much for updating the code. it means a lot. actually my problem is based on refining the mesh inside a cylinder and i have tried with normal library, trust me. but it doesn't work half as good as yours. i installed your extension in OF2.3.1 and the RESULTS WERE COMPARABLE TO GERRIS OCTREE. i ran the dam break problem and i used refinement level 8 with just a single mesh consisting of two blocks only. your library i feel makes it faster too. i'll try installing for my openfoam-v1612+ and let you know if it works |
||
August 24, 2017, 05:54 |
|
#6 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10 |
i am sure you aware of gerris flow solver. actually gerris uses octree meshing. the simulation starts with a single block and refinement takes place only where the interface is. the adaptive meshing of gerris is by far the best. your code follows the same refinement followed by gerris
there is no need for giving a base mesh like in openfoam. simple geometries can be given by means of equation: cylinder by x2+y2=r2 and extruding in z..simple a good visulaization of gerris octree meshing is here; https://www.youtube.com/watch?v=yJhYENolQeU https://www.youtube.com/watch?v=3Jf_-XMv1wk the discussion on it's tree strucure/adaptivity is given here. in short gerris:solves for each finer subgrid to capture flow physics,to reach kolmogrov length scales https://www.youtube.com/watch?v=DKn9iuD7Ihk i believe openfoam lags in this department. GERRIS adaptivemeshing comes so close to DNS!. |
|
August 24, 2017, 06:30 |
|
#7 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52 |
I am sorry but I never heard about Gerris. Seems that it is not longer in development. Thanks for the information but I will not check it out. If you are so happy with that software, use it.
__________________
Keep foaming, Tobias Holzmann |
|
August 24, 2017, 07:35 |
|
#8 |
Member
Ali Noaman Ibrahim
Join Date: Sep 2015
Location: US_Chicago
Posts: 97
Rep Power: 11 |
That is really beneficial discussion. However, I would like to add that Gerris is not that effective for handling the contact angle phenomena particularly for 3D problems. That why I switched towards OpenFoam after working with gerris for a while.
Now, I have solved very complicated cases of interface tracking. However, I never used the adaptive mesh refinement. Any idea about links, examples , ... etc preferably if there is something works with Salome? Thank you all Sent from my iPhone using CFD Online Forum mobile app |
|
August 24, 2017, 09:05 |
|
#9 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 |
FYI: Gerris has long been superseeded by PARIS: http://www.ida.upmc.fr/~zaleski/paris/index.html
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
August 24, 2017, 09:09 |
|
#10 |
Member
Ali Noaman Ibrahim
Join Date: Sep 2015
Location: US_Chicago
Posts: 97
Rep Power: 11 |
That is true ( I am following Gerris news!) but not for tracking the contact line in its 3d fashion.
Regards |
|
August 24, 2017, 11:27 |
|
#11 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10 |
Yes , Gerris is no longer in development.
it's been followed with basilisk which is still in process of testing and development, paris is also a good. but openFoam needs unstructured adaptive refinement. because adaptive meshing saves time when it comes to simulating big & complex geometric dimensions. it's not possible to figure out snappyhexmesh problems and keep trying until you get it right. and unless you have hexahedral mesh openfoam doesn't refine.it's just not possible to make hexahedral mesh for some geometries at all take a look at these posts, Adaptive mesh and movement |
|
August 24, 2017, 11:41 |
|
#12 |
Member
Ali Noaman Ibrahim
Join Date: Sep 2015
Location: US_Chicago
Posts: 97
Rep Power: 11 |
Thanks Saddy. I will look at take a look at that post.
Sent from my iPhone using CFD Online Forum mobile app |
|
August 24, 2017, 15:04 |
|
#13 | |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10 |
Quote:
and yeah i installed your OpenFOAM5.0 it compiled after 1 error: line1662 i changed writeobjects to writeobject. unfortunately it compiles but doesn't work on the way it should. i'll try to sit down with the codes and try to compile it again.!!! thanks |
||
August 26, 2017, 03:54 |
|
#14 |
Senior Member
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10 |
PROFESSOR Zaleski ia coming to our college in september to give guest lectures
|
|
Tags |
interface refinement, refinefvmesh |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
question about CFD-DTF libary | jackson | Main CFD Forum | 0 | September 12, 2009 01:06 |