|
[Sponsors] |
[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 1, 2012, 19:38 |
|
#21 |
Member
Join Date: Sep 2012
Posts: 51
Rep Power: 14 |
Hi,
thank you for the development of a large deformation solid solver. I'm planning to use it for fsi in arteries and I was wondering : do you have any idea of the amount of work required to use it on OF 2.1.1 ? (and which functionalities would be lost). Thanks Florian |
|
October 2, 2012, 02:04 |
|
#22 |
New Member
kob
Join Date: Nov 2011
Posts: 28
Rep Power: 15 |
Hi Philip
I always get error when I use the following command in the floder OpenFOAM-1.6-ext: Code:
git checkout -b feature/solidMechanics remotes/origin/feature/solidMechanics Thank you! bryant |
|
October 2, 2012, 06:01 |
|
#23 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
Quote:
I don't think it should be too difficult to get it to work in standard OpenFOAM, but I haven't tried. As regards losing functionality, I don't think anything should be lost with the large strain solver. I recommend you take the solver of interest and take the required bits out of the solidModels library (rheologyModel and boundary conditions - fvPatchFields) and try compile it. It should be possible to get it working, with a little help for the smart forum people. Best regards, Philip |
||
October 2, 2012, 06:03 |
|
#24 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
Quote:
You can browse and download the code manually from the source-forge here. What error do you get when you run the git command? Philip |
||
October 2, 2012, 15:28 |
|
#25 | |
Member
Join Date: Sep 2012
Posts: 51
Rep Power: 14 |
Quote:
Regards, Florian |
||
October 3, 2012, 00:50 |
|
#26 |
New Member
kob
Join Date: Nov 2011
Posts: 28
Rep Power: 15 |
Hi Philip
The error I met just like someone have mentioned before.When I run the following command: Code:
git checkout -b feature/solidMechanics remotes/origin/feature/solidMechanics Code:
fatal: Not a git repository (or any of the parent directories): .git regards! bryant |
|
October 3, 2012, 04:48 |
|
#27 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
Quote:
OK I believe the git command does not work for you because you did not download OpenFOAM-1.6-ext via git, is this the case? You can download OpenFOAM-1.6-ext via git using the command: Code:
git clone git://openfoam-extend.git.sourceforge.net/gitroot/openfoam-extend/openfoam-extend Code:
git checkout -b feature/solidMechanics remotes/origin/feature/solidMechanics Philip |
||
October 27, 2012, 18:45 |
Brilliant tutorials! However...
|
#28 |
New Member
Join Date: Jun 2012
Posts: 7
Rep Power: 14 |
Hi there,
Thanks for the upload, a very good contribution, but I'm having trouble compiling the icoFsiElasticNonLinULSolidFoam code. Please find attached my log file for the compiling process. Any ideas how to fix this? Has anyone got the icoFsiElasticNonLinULSolidFoam code working?? Best Regards, Steven solid_mechanics.txt |
|
October 30, 2012, 05:50 |
|
#29 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
Quote:
The reason for this error is due to a missing lnInclude links directory in $FOAM_SRC/tetDecompositionFiniteElement, so to solve the problem: Code:
cd $FOAM_SRC/tetDecompositionFiniteElement wmakeLnInclude . The solver should then compile. Best regards, Philip |
||
October 30, 2012, 10:06 |
Thank you
|
#30 |
New Member
Join Date: Jun 2012
Posts: 7
Rep Power: 14 |
Dear Philip,
Many thanks for your help! The code is working now. Best Regards, Steven |
|
November 3, 2012, 14:56 |
|
#31 |
Member
Join Date: Sep 2012
Posts: 51
Rep Power: 14 |
[FIXED] it finally seems to bee working after changing the BC
Hi, after trying some of the tutorials I wanted to test your solver on my cases, especially the large deformation FSI solver. For testing purposes I set up a very simple pipe with a fluid inside (the fluid mesh being generated with snappyHexMesh and the solid via extrudeToRegionMesh). The BC are : fluid : inlet : fixed velocity, zero gradient pressure outlet : zero gradient velocity, fixed value pressure wall : movingWallVelocity, zero pressure and fixedValue (0 0 0) for motionU solid : outlet/inlet : fixedDisplacement (0 0 0) inter/outer wall : type solidTraction; nonLinear updatedLagrangian; traction uniform ( 0 0 0 ); pressure uniform 0; value uniform (0 0 0); I also changed the gradient scheme for DU as the code politely asked me from Gauss linear to extendedLeastSquares 0. The problem is, it can't pass the first iteration and crash in tears with a floating point exception : Code:
~$ icoFsiElasticNonLinULSolidFoam /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM Extend Project: Open source CFD | | \\ / O peration | Version: 1.6-ext | | \\ / A nd | Web: www.extend-project.de | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 1.6-ext-959ec266ba5c Exec : icoFsiElasticNonLinULSolidFoam Date : Nov 03 2012 Time : 14:42:05 Host : barbaportable PID : 21978 Case : /home/florian/OpenFOAM/florian-1.6-ext/run/testnoumerooneFSI/HronTurekFsi/fluid nProcs : 1 SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create dynamic mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: laplaceFaceDecomposition Selecting motion diffusivity: quadratic Reading transportProperties Reading field p Reading field U Reading/calculating face flux field phi Reading incremental displacement field DU Patch region0_to_wall4_wall3 Traction boundary field: DU nonLinear set to updated Lagrangian max(detdd) = 17.9352 min(detdd) = 5.69387 average(detdd) = 9.90012 max(detdd) = 17.9352 min(detdd) = 5.69387 average(detdd) = 9.90012 nAddCells/nCells = 0 Reading incremental displacement field DV Reading accumulated velocity field V Reading accumulated stress field sigma Reading incremental stress field DSigma Selecting rheology model linearElastic Reading coupling properties Create fluid-to-solid and solid-to-fluid interpolators Check fluid-to-solid and solid-to-fluid interpolators Fluid-to-solid face interpolation error: 0.00131215 Solid-to-fluid face interpolation error: 0.00131215 Starting time loop Time = 1e-05 Selecting coupling scheme Aitken Time = 1e-05, iteration: 1 Current fsi under-relaxation factor: 0.01 Maximal accumulated displacement of interface points: 0 Courant Number mean: 1.8894e-06 max: 1.64473e-05 velocity magnitude: 0.2 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 7.13165e-12, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 6.85799e-12, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.94807e-13, No Iterations 1 GAMG: Solving for p, Initial residual = 1, Final residual = 5.14394e-07, No Iterations 20 GAMG: Solving for p, Initial residual = 0.0179676, Final residual = 7.3011e-07, No Iterations 8 time step continuity errors : sum local = 5.29294e-12, global = -2.04832e-12, cumulative = -2.04832e-12 GAMG: Solving for p, Initial residual = 0.00179327, Final residual = 6.69856e-07, No Iterations 6 GAMG: Solving for p, Initial residual = 0.000394718, Final residual = 7.74832e-07, No Iterations 3 time step continuity errors : sum local = 5.71226e-12, global = -1.40003e-12, cumulative = -3.44834e-12 Setting traction on solid patch Total traction force = (-0.071344 0.125002 1.21712) Exception en point flottant Both solid and fluid domains can be seen in the attached screenshot. By the way, what is the acceptable interpolation error ? Thank you in advanced for you help, regards, Florian Last edited by Hiroshiman; November 5, 2012 at 13:57. |
|
November 7, 2012, 06:13 |
|
#32 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
Quote:
For future reference, what BC did you have to change to solver your problem? I am not sure what interpolation error is acceptable, it depends on your application and what you are interested in. More information can be found on this solver in the conference paper by Tukovic and Jasak here. Philip |
||
November 7, 2012, 10:30 |
|
#33 | |
Member
Join Date: Sep 2012
Posts: 51
Rep Power: 14 |
Hi, thank you for the answer.
Quote:
motionU Code:
boundaryField { region0_to_wall4_wall3 { type fixedValue; value uniform (0 0 0); } minZ { type fixedValue; value uniform (0 0 0); } maxZ { type fixedValue; value uniform (0 0 0); } } Code:
boundaryField { region0_to_wall4_wall3 { type fixedValue; value uniform 0; } minZ { type fixedValue; value uniform 0; } maxZ { type zeroGradient; } } Code:
boundaryField { region0_to_wall4_wall3 { type movingWallVelocity; value uniform (0 0 0); } minZ { type zeroGradient; } maxZ { type fixedValue; value uniform (0 0 0.2); } } Do have any any if anyone has already implemented 0D BC on the solid part ? Regards, Florian Last edited by Hiroshiman; November 7, 2012 at 13:24. |
||
November 7, 2012, 19:12 |
|
#34 |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
||
November 9, 2012, 21:02 |
|
#35 |
Member
Join Date: Sep 2012
Posts: 51
Rep Power: 14 |
||
November 10, 2012, 13:09 |
|
#36 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
Quote:
I am not aware of such boundary conditions, but I think they should be straight-forward to implement. The BC would be based on solidTraction and then the traction_ could be calculated each increment based on the displacement using what ever spring model you like. Philip |
||
November 11, 2012, 03:50 |
|
#37 | |
Member
Join Date: Sep 2012
Posts: 51
Rep Power: 14 |
Quote:
Thank you for your answer, I'll look at it more closely. Florian |
||
December 16, 2012, 16:24 |
|
#38 |
Member
Johannes Baumann
Join Date: Mar 2009
Location: Baden-Wuerttemberg, Germany
Posts: 43
Rep Power: 17 |
Hi Philip,
first of all, many thanks to you and your group for sharing the Solid Mechanics Solvers. Must have been quite a bit of work... I'd love to have a look at the new solvers, but unfortunately, I get build errors when using gcc-4.7.1. Do you have a gcc-47 update available (maybe in form of a seperate branch) or do I have to switch back to gcc-46, at least for now? Thanks & best regards, Johannes |
|
December 16, 2012, 16:46 |
|
#39 | |
Super Moderator
Philip Cardiff
Join Date: Mar 2009
Location: Dublin, Ireland
Posts: 1,095
Rep Power: 34 |
Quote:
Hopefully you will find the solvers useful. I have not tried gcc-47 yet so I have no experience with it. Changing back to gcc-46 is probably the quickest way to get things working. If you post the errors here then we can try to fix them. Best regards, Philip Last edited by bigphil; December 16, 2012 at 16:46. Reason: typo |
||
December 18, 2012, 05:13 |
|
#40 |
Member
Johannes Baumann
Join Date: Mar 2009
Location: Baden-Wuerttemberg, Germany
Posts: 43
Rep Power: 17 |
Hi Philip,
thanks for your fast response. Luckily, my problem has been a false alarm. I don't know exactly how this happened but somehow during git checking out/compiling/cleaning I managed to kill my OF-extend binaries, so all the errors were just dependency issues. Yesterday I recompiled everything including the solidMechanics branch without a single error so everything is working now. Best regards, Johannes |
|
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 |