|
[Sponsors] |
[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 1, 2013, 07:13 |
|
#81 | ||
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Hi Danny,
Quote:
The accuracy for these solid solvers is essentially determined by the accuracy of the gradient operator, so use extendedLeastSquares gradScheme. And also make sure you use the included boundary conditions which have boundary non-orthogonal correction, which is imperative for unstructured meshes with these solvers. Quote:
Best regards, Philip |
|||
August 7, 2013, 20:40 |
|
#82 |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
@ bigphil
I am trying to use icoFsiElasticNonLinULSolidFoam in parallel. I am running the tutorial case you have provided. but I get the following error after some time: Code:
Time = 0.065, iteration: 3 Current fsi under-relaxation factor (Aitken): 0.362946 Maximal accumulated displacement of interface points: 0.00876687 [0] [1] [1] [1] --> FOAM FATAL ERROR: [1] face 41 area does not match neighbour by 0.0103705% -- possible face ordering problem. patch: procBoundary1to0 my area:0.000239418 neighbour area: 0.000239443 matching tolerance: 0.0001 Mesh face: 42054 vertices: 4((0.574801 0.214706 0.025334) (0.574801 0.214706 -0.025334) (0.574854 0.209981 -0.025334) (0.574854 0.209981 0.025334)) Rerun with processor debug flag set for more information. [1] [1] From function processorPolyPatch::calcGeometry() [1] in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 216. [1] FOAM parallel run exiting [1] [0] [0] --> FOAM FATAL ERROR: [0] face 41 area does not match neighbour by 0.0103705% -- possible face ordering problem. patch: procBoundary0to1 my area:0.000239443 neighbour area: 0.000239418 matching tolerance: 0.0001 Mesh face: 42167 vertices: 4((0.574801 0.214706 0.025334) (0.574854 0.20998 0.025334) (0.574854 0.20998 -0.025334) (0.574801 0.214706 -0.025334)) Rerun with processor debug flag set for more information. [0] [0] From function processorPolyPatch::calcGeometry() [0] in file meshes/polyMesh/polyPatches/constraint/processor/processorPolyPatch.C at line 216. [0] FOAM parallel run exiting [0] http://www.rodfile.com/tgp64725opo7 |
|
August 12, 2013, 12:18 |
|
#83 | ||
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Hi Daniel,
Quote:
So yes the solver should work with flow inside a deformable pipe. Quote:
http://powerlab.fsb.hr/ped/kturbo/Op...Zagreb2007.pdf http://www.tfd.chalmers.se/~hani/kur.../FSIslides.pdf The partitioned FSI procedure is described in a number of references, for example in Aleksandar Karac's thesis: http://powerlab.fsb.hr/ped/kturbo/Op...arKaracPhD.pdf For your error above, I am not entirely sure of the reason, people have previously mentioned this issue. Is the model converging to the defined tolerance? Also you are correct that the IQN-ILS scheme coupling scheme is not currently implemented, even though it is given as an option. Philip |
|||
August 12, 2013, 14:03 |
|
#84 | |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
Quote:
everything is just fine till t=0.065 sec.. before this time the solution was converging just like the serial run. I was simulating the provided case without any changes. |
||
August 21, 2013, 06:46 |
|
#85 |
New Member
Join Date: Mar 2013
Posts: 17
Rep Power: 13 |
Hi Philip,
Just a small question, could you post the blockMeshDict and an easy way to generate the faceZones file for the FSI test case? This would be of great help. I'm currently developing a FSI solver with space mapping and multi level acceleration, so I need an easy way to generate new grids. Looking forward hearing from you! Best regards, David Blom |
|
August 21, 2013, 07:56 |
|
#86 | ||
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Quote:
Unfortunately the mesh for the FSI test case was created (not by me) using Gambit, so there is no blockMeshDict. If you are familiar with blockMesh, it should be relatively straight-forward to create the blockMeshDict for this case, albeit time consuming. Quote:
I have had a look at this processor faces error: could you try out the attached "moveSolidMeshLeastSquares.H" (place this in the solver directory and execute "wmake"). I have added a few lines that should make sure that the processor boundaries are exactly the same. Let me know if it works for you. Best regards, Philip |
|||
August 21, 2013, 08:39 |
|
#87 | |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
Quote:
I have recompiled the the code with "moveSolidMeshLeastSquares.H" you've provided. I get floating point exception in both parallel and serial runs after only one time step. regards |
||
August 21, 2013, 08:56 |
|
#88 |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
I think I got the error fixed.
I have changed the following in the code you've attached: for newPoints I changed "mesh" to "stressMesh" (solid mesh) I will test and report if the main error is eliminated or not. regards edit: I can confirm the problem has been fixed. I have attached the modification. thank you very much Philip |
|
August 21, 2013, 09:55 |
|
#89 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Quote:
Yep, sorry you are correct 'mesh' is meant to be 'stressMesh' - I tried out this code on an UL solid solver which calls the solid mesh 'mesh'. I will add this change to the git. Philip |
||
August 21, 2013, 10:06 |
|
#90 | |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
Quote:
does it need huge modification on the base code? regards |
||
August 21, 2013, 10:16 |
|
#91 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Quote:
As far as I know, there will be quite a bit of work getting it working, but don't let me deter you Philip |
||
August 21, 2013, 13:20 |
Error
|
#92 |
New Member
Soheil nazmdeh
Join Date: Feb 2013
Posts: 15
Rep Power: 13 |
Hi,
i try to add solid mechanics but this error seen, I've already download openFoam 1.6 ext from openFoam wiki. please help me to add solidMechanics fluidman@fluidman-Aspire-5755G:~/OpenFOAM/OpenFOAM-1.6-ext$ git checkout -b feature/solidMechanics remotes/origin/feature/solidMechanics fatal: Cannot update paths and switch to branch 'feature/solidMechanics' at the same time. Did you intend to checkout 'remotes/origin/feature/solidMechanics' which can not be resolved as commit? |
|
August 21, 2013, 13:45 |
|
#93 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Quote:
I am far from an expert in git, but try the following: Code:
git remote update git fetch git checkout -b feature/solidMechanics origin/feature/solidMechanics |
||
August 21, 2013, 14:27 |
|
#94 | |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
Quote:
The problem is when I change rhoSolid=1000 and rhoFluid=998 after a few iterations I get floating point exception error on solidSolver for DU. thanks for the help man |
||
August 21, 2013, 14:40 |
solid mechanic solver
|
#95 |
New Member
Soheil nazmdeh
Join Date: Feb 2013
Posts: 15
Rep Power: 13 |
thank you philip.
It does not work... Can you send me the Solvers ? Soheil.nazmdeh@gmail.com thank you |
|
August 21, 2013, 17:31 |
|
#96 |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
@ Philip
I am very sorry but I was mistaken, the problem is still there. I was in hurry and entered the parallel command with out noticing that I have drop the main part: I used "mpirun -np 2 icoFsiE...." instead of "mpirun -np 2 icoFsi... -parallel" I have noticed that when I got home. regards |
|
August 23, 2013, 14:06 |
|
#97 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,093
Rep Power: 34 |
Quote:
Philip |
||
August 26, 2013, 15:31 |
|
#98 |
Senior Member
Daniel
Join Date: Mar 2013
Location: Noshahr, Iran
Posts: 348
Rep Power: 21 |
@ bigphil
today I have compiled OpenFOAM-1.6-ext with openmpi-1.5 but the problem is still there. could something be wrong with the tutorial case setup or you don't have this problem on your own system? regards |
|
September 2, 2013, 06:11 |
|
#99 |
New Member
Join Date: Mar 2013
Posts: 17
Rep Power: 13 |
Hi Philip,
I noticed another small problem with the mesh deformation code of the structure solver. The fluid-structure interface isn't matching anymore after a couple of time steps at the attachment location to the cylinder. Is it possible that the plateFix patch is not taken into account by the mesh motion algorithm? Thanks in advance! Best regards, David |
|
September 4, 2013, 02:49 |
FSI in parallel
|
#100 |
New Member
Join Date: Jul 2012
Posts: 4
Rep Power: 14 |
Hi everyone,
I had a look at Philips FSI solver and am especially excited about the possibility to run it in parallel, as icoFsiFoam is apparently not running in parallel yet. However, I managed to run the HronTurekFsi case in parallel but I can't reconstruct the data.. the results of the processors are all stored under the fluid/processor*/"time" folder where also the solid results are located within another folder (fluid/processor*/"time"/solid). When I run reconstructPar that folder is ignored and only the fluid data is reconstructed. Probably I just make a stupid fault.. Any help would be appreciated! Regards, Anna Update: problem solved. not sure whether this is the most elegant way, but I wrote a batch sorting things out! Last edited by sloti; September 4, 2013 at 03:55. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
GPU Linear Solvers for OpenFOAM | gocarts | OpenFOAM Announcements from Other Sources | 37 | August 17, 2022 15:22 |
[Virtualization] OpenFOAM oriented tutorial on using VMware Player - support thread | wyldckat | OpenFOAM Installation | 2 | July 11, 2012 17:01 |
New OpenFOAM Forum Structure | jola | OpenFOAM | 2 | October 19, 2011 07:55 |
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 | wyldckat | OpenFOAM Announcements from Other Sources | 3 | September 8, 2010 07:25 |
OpenFOAM Debian packaging current status problems and TODOs | oseen | OpenFOAM Installation | 9 | August 26, 2007 14:50 |