Hi everyone,
i think i've n
Hi everyone,
i think i've noticed an error in makeAxialMesh. Comparing a wedge realized with gmsh and another with makeAxialMesh, i've noticed that the wedge angle is 2.5°, not 5°. Checking the code, i can see a scalar variable "angle" set to 2.5° in order to revolve -2.5° and +2.5°, but i also see a 0.5 factor in the "factor" variable. I think that the factor should be set to 1.0 or the angle to 5°, am I wrong ? Sincerely, Etienne. |
Hi!
Etienne convinced me on
Hi!
Etienne convinced me on the IRC that he is right (and I am wrong). The mistake is now corrected in the 1.4.1-version on the SVN. In the 1.5 version it was already corrected Bernhard |
Hi everyone !
I've found an
Hi everyone !
I've found another problem in makeAxialMesh ;) (sorry Bernhard to bother you ;)). This one is present in the 2 versions. But I think we have to discuss if we correct it or not. In fact, makeAxialMesh does not revolve the mesh, it projects the mesh on wedge planes. I agree that the error is really small (depending on the radius of the problem), but it caused me great problems for coincidence with another mesh. I've added 2 lines and modified one to solve the problem in 1.4.1, the same modif can be applied in the 1.5 version. The question is : "should we consider this point as a modelization problem ?" Regards, Etienne. |
Hello,
I've modified the co
Hello,
I've modified the code and added a revolve option. This change should be transparent for users who don't want this feature. The 2 versions 1.4.1 and 1.5 have been updated. To activate it, add '-revolve' to the command line or change it in the dictionary if you use the rotationDict. Check the output when running makeAxialMesh, you'll see 'Revolving nodes' or 'Projecting nodes' Etienne. |
Just tried to compile makeAxia
Just tried to compile makeAxialMesh on 1.5-dev. It is looking for repatchPolyTopoChanger.H, which was removed in 1.5-dev (I can't find it in src/dynamicMesh). Given that makeAxialMesh is distributed through the OpenFOAM-extend repository, I figured it would have been tested and fixed for the -dev versions. Anyone have a solution?
Eric |
Hello Eric,
you're right, m
Hello Eric,
you're right, makeAxialMesh does not compile with 1.5-dev because of slight divergences between the official release and the svn release. I've searched quickly but I've not found a way to solve this problem easily. Sincerely, Etienne. |
Quote:
Bernhard |
makeAxialMake and collapseEdges
Dear Bernhard,
I was able to download and compile the version of makeAxialMesh for OpenFOAM v1.5. It seems to run OK, but I am having trouble running the test case axialCavity. I think the problem is that I am unsure how to use the routines. Firstly, I run blockMesh, then I run the command makeAxialMesh -axis fixedWalls -wedge frontAndBack This seems to give no errors, but writes the mesh to time 0.00125 (why?). In the same directory I run the command collapseEdges 0.0001 175 (small length and large angle as suggested in the wiki) no small edges are collapsed, but this is looking in the time=0 directory. The new mesh is in the time=0.00125 directory. However, when I force collapseEdges to look at the time=0.00125 directory, I get 21 collapsed small edges, the results of which are written in the time=0.0025 directory (again why?). Is the idea to copy the keep copying the files generated by makeAxialMesh and collapseEdges back into constant/polyMesh? Thanks, Peter |
Quote:
If you can do without undo: most mesh utilities have an option -overwrite .... Quote:
Bernhard |
Quote:
Quote:
Quote:
|
I simply used the standard v1.5 and not v1.5-dev. I did not need to do anything special. However, I guess that is not much use to you. Sorry.
Regards, Peter. |
Quote:
Bernhard |
[/QUOTE]
Quote:
All right , i i used the latest from the SVN . It compiled successfullyy . But what is the command . If the use the command: makeAxialMesh . . axis frontAndBackPlanes as in 1.4 , it doesnt not work .. What are the arguments that have to be specified ? |
Mihir,
Take a look in cd $WM_PROJECT_USER_DIR/applications/MakeAxialMesh for the TestCases directory. Look at how the parameters are now passed to makeAxialMesh through the rotation dictionary: system/rotationDict |
Hi,
I hope this is the right post. I'm trying to use makeaxialmesh but I found some difficulties. I have this 2d mesh (1 cell thick) converted from a 2d fluent mesh using fluentofoam, with the axis edge named axis and the opposite edge (the one that has to be splitted in two edges) named wall. In the terminal I went into the folder of the case and I gave this command: makeAxialMesh axis wall and I get this error: Create time Create mesh for time = 0 #0 Foam::error::printStack(Foam::Ostream&) in "/home/lions85/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigSegv::sigSegvHandler(int) in "/home/lions85/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib/libc.so.6" #3 std::basic_string<char, std::char_traits<char>, std::allocator<char> >::basic_string(std::string const&) in "/home/lions85/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux64/lib64/libstdc++.so.6" #4 main in "/home/lions85/OpenFOAM/lions85-1.5/applications/bin/linux64GccDPOpt/makeAxialMesh" #5 __libc_start_main in "/lib/libc.so.6" #6 _start in "/home/lions85/OpenFOAM/lions85-1.5/applications/bin/linux64GccDPOpt/makeAxialMesh" Segmentation fault Probably I'm making a mistake, but I can not find it simply. Is someone able to help me? Thank you very much EManuele |
Quote:
makeAxialMesh axis frontAndBackPlanes |
Sorry,
may you explain it better?? I didn't understand very well. Anyway i tried with the command makeAxialMesh axis frontAndBackPlanes and i got Create time Create mesh for time = 0 #0 Foam::error::printStack(Foam::Ostream&) in "/home/lions85/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigSegv::sigSegvHandler(int) in "/home/lions85/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib/libc.so.6" #3 std::basic_string<char, std::char_traits<char>, std::allocator<char> >::basic_string(std::string const&) in "/home/lions85/OpenFOAM/ThirdParty/gcc-4.3.1/platforms/linux64/lib64/libstdc++.so.6" #4 main in "/home/lions85/OpenFOAM/lions85-1.5/applications/bin/linux64GccDPOpt/makeAxialMesh" #5 __libc_start_main in "/lib/libc.so.6" #6 _start in "/home/lions85/OpenFOAM/lions85-1.5/applications/bin/linux64GccDPOpt/makeAxialMesh" Segmentation fault Thank you Emanuele P.S. consider that i have my boundary edges that I created with gambit and then I have two other boundaries created during the conversion from fluent to foam and one of this is named automatically FrontAndBackPlanes |
Quote:
|
Quote:
The other thing: is there any other application that you successfully compiled before? Is your version of OF a downloaded binary or self-compiled? Because I remotely remember having seen a similar stack-trace in a previous discussion and it was a compilation problem or something like that Bernhard |
Hi, dear all,
I tried to revolve (select withMesh option)my face in Gambit, and I can only get a volumn, but the mesh is not created in the volumn. Does anybody know the answer? sorry, I am away from the topic. |
All times are GMT -4. The time now is 20:41. |