CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Solver for FSI strongly coupled

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree3Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   August 2, 2011, 06:55
Default
  #41
Member
 
Santiago
Join Date: Dec 2009
Posts: 85
Rep Power: 7
gascortado is on a distinguished road
You have to compile the solver first. See this one http://www.cfd-online.com/Forums/ope...tml#post317171
Good luck!
gascortado is offline   Reply With Quote

Old   August 3, 2011, 00:11
Default
  #42
Member
 
Santiago
Join Date: Dec 2009
Posts: 85
Rep Power: 7
gascortado is on a distinguished road
Hi again. I'm trying to set up a FSI between a horizontal plate inside a fluid with the leading edge of the plate fixed. The basic geometry is similar to this one Meshing multiple 3D regions

I want to test Charbel's solver for a 3D case using LES. I successfully modified the original consoleFluid so LES is used instead of laminar flow.

I create 2 meshes in gambit, one for the surrounding fluid (with the plate as a void) and another one with the plate only as my solidMesh. I convert them using fluent3DmeshToFoam. When I run the case I get the following error:

--> FOAM FATAL ERROR:
Problem with patch interpolation definition

From function testVorflexFoam
in file readCouplingProperties.H at line 30.

FOAM aborting

Aborted


the readCouplingProperties.H looks like this:

// Read solid patch data
word solidPatchName(couplingProperties.lookup("solidPat ch"));

label solidPatchID =
stressMesh.boundaryMesh().findPatchID(solidPatchNa me);


// Read fluid patch data
word fluidPatchName(couplingProperties.lookup("fluidPat ch"));

label fluidPatchID =
mesh.boundaryMesh().findPatchID(fluidPatchName);


if (solidPatchID < 0 || fluidPatchID < 0)
{
FatalErrorIn(args.executable())
<< "Problem with patch interpolation definition"
<< abort(FatalError);
}


Do you guys know what is the meaning of this error? Is my procedure right?
I'm attaching my boundary files for both solid and fluid meshes. I will appreciate any comments or information. Thanks

BTW: What does the MotionU file do?
Attached Files
File Type: txt boundaryFluid.txt (1.8 KB, 24 views)
File Type: txt boundary Solid.txt (1.1 KB, 12 views)
File Type: txt Usolid.txt (1.3 KB, 18 views)
File Type: txt U.txt (5.9 KB, 10 views)
gascortado is offline   Reply With Quote

Old   August 3, 2011, 05:08
Default
  #43
New Member
 
charbel's Avatar
 
Charbel
Join Date: Jan 2011
Location: Lebanon
Posts: 27
Rep Power: 6
charbel is on a distinguished road
Hi Santiago,
Go into ~/fluid/constant/couplingProperties and modify the patch names as follows:
Code:
solidPatch      plate;

fluidPatch      plate;

movingRegion    region0;
You can use pimpleFsiFoam where some errors are corrected.
arashfluid likes this.
__________________
--
Charbel
charbel is offline   Reply With Quote

Old   September 9, 2011, 10:43
Default
  #44
Member
 
Santiago
Join Date: Dec 2009
Posts: 85
Rep Power: 7
gascortado is on a distinguished road
Hi,

I'm trying to extract the pressure that the fluid transfers to the plate. I see that at "setPressure.H" there is a scalarField created called "solidPatchPressure". Any idea of what I need to do in order to get these values and print them to a file? I want to take a look at those values and then apply the pressure to a finite element code to compare deformations. Any help?

Thanks
gascortado is offline   Reply With Quote

Old   September 16, 2011, 13:52
Default
  #45
Member
 
Santiago
Join Date: Dec 2009
Posts: 85
Rep Power: 7
gascortado is on a distinguished road
While trying to verify the pressure force exerted to the plate I come across the following line of code in setPressure.h

vector totalPressureForce =
sum
(
p.boundaryField()[fluidPatchID]*
mesh.Sf().boundaryField()[fluidPatchID]
);


While reading the pressure generated by p.boundaryField()[fluidPatchID] i see that the pressure is either negative or positive as expected.

Then while reading
mesh.Sf().boundaryField()[fluidPatchID] I see that we get the area vectors for each boundary cell and each of the component of this vector can be positive and negative.

Now, when the product of the area vector and the pressure is done it seems to me that the final sign of the pressure Force will be screwed up since the scalar pressure (+ or -) is multiplied by the area vector components (also + or -). Shouldn't if be done so the pressure magnitud at the cell center is multiplied by the absolute value of the the area of the cell??

Does anyone know if this is a bug or if I am mistaken? Thanks
gascortado is offline   Reply With Quote

Old   September 19, 2011, 05:12
Default Has anyone faced the same probolem?
  #46
ICL
New Member
 
Ali Sh
Join Date: Jun 2009
Location: London
Posts: 28
Rep Power: 8
ICL is on a distinguished road
I tried to modify fluid solver of icoFSIFoam to ultimately model EHL in ball bearing. The computation stops after a few iterations. In order to re-check the fluid solver with extended version, I tried to solve the fluid region separately. The computation diverged after a number of iteration while the same fluid solver on the same domain converged using OF 2.0.0 and OF 1.5!


ICL is offline   Reply With Quote

Old   October 28, 2011, 06:16
Default
  #47
New Member
 
Michael Bruckner
Join Date: Apr 2009
Location: France
Posts: 27
Rep Power: 8
michaelb is on a distinguished road
Hi Charbel,

first of all thanks for sharing this FSI solver. I successfully compiled and ran it.
Could you confirm that the computation is very slow on a single core : about 25 hours for 1sec of runTime ?
Did you managed to set up a parallel version of it with e.g. in the Allrun file ?:
decomposePar
runParallel $application 4
reconstructPar

I tried including a controlDecomposeParDict in my fluid/system folder but whithout any success : it is not working since the coupling with the solid solver is not done properly.
I'd be very interested if you would came up with a solution like specifying X procs for the fluid and Y procs for the solid !
thanks

Michael

michaelb is offline   Reply With Quote

Old   December 27, 2011, 04:24
Default
  #48
Member
 
Join Date: Nov 2009
Posts: 61
Rep Power: 7
waku2005 is on a distinguished road
Hi Charbel,

Thanks a lot for sharing your FSI solver, pimpleFsiFoam.
I also successfully compiled and ran the test case attached with the previous version in the thread #14.

I'd like to ask you whether this solver can be applied to lower rhof/rhos cases.
My application is deformation of beam like structure (a silicon rubber block) in a water with retatively slow flow velocity of 0.1 m/s.
rho of fluid (water) is around 1000 kg/m3 and silicon rubber has almost the same, thus rhof/rhos is nearly unity.

Even I changed fluid density from 1 to 1000 kg/m3, and adjusted time step, solver diverged in several time steps.

Any comment will be highly appreciated.
Best regards,
waku
waku2005 is offline   Reply With Quote

Old   January 24, 2012, 09:29
Default
  #49
Senior Member
 
Join Date: Jun 2011
Posts: 141
Rep Power: 6
mechy is on a distinguished road
Quote:
Originally Posted by gascortado View Post
Thanks for your reply Matz,

When I use the browser I can see the folder lnInclude inside $WM_PROJECT_DIR/src/tetDecompositionFiniteElement directory. Inside lnInclude I see a whole bunch of symbolic files among all tetFemMatrices.H.

However, when I try accessing
$WM_PROJECT_DIR/src/tetDecompositionFiniteElement/lnInclude by using the terminal the folder is not recognized. Furthermore, by using the "ls" command, the folder lnInclude is not listed:

santiago@santiago-laptop:/usr/lib/OpenFOAM-1.6-ext/src/tetDecompositionFiniteElement$ ls

Allwmake tetFemCalculus tetPolyMesh
elementMesh tetFemMatrix tetPolyMeshCellDecomp
fields tetFemMethod tetPolyMeshFaceDecomp
Make.cellDecomp tetFemSolution tetPolyPatchInterpolation
Make.faceDecomp tetPointMesh


When I tried doing:
cd $WM_PROJECT_DIR/src/tetDecompositionFiniteElement/lnInclude

bash: cd: /usr/lib/OpenFOAM-1.6-ext/src/tetDecompositionFiniteElement/lnInclude: No such file or directory

I then switched to this directory santiago@santiago-laptop:/usr/lib/OpenFOAM-1.6-ext/src/tetDecompositionFiniteElement$ ln -s ../tetFemMatrix/tetFemMatrices.H tetFemMatrices.H

ln: creating symbolic link `tetFemMatrices.H': Permission denied

I wonder why is the folder not showing up in the terminal. Does this mean that there is something wrong with the binary files?

This is very frustrating. Thanks again


dear gascortado
I have the same problem how I can remove this error
I hvae ubuntu 10.04
please tell me how I can solve the Permission denied error ?

best regards
mechy is offline   Reply With Quote

Old   January 25, 2012, 19:07
Default
  #50
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,258
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings mechy,

Quote:
Originally Posted by mechy View Post
please tell me how I can solve the Permission denied error ?
Try using sudo:
Code:
sudo ln -s ../tetFemMatrix/tetFemMatrices.H tetFemMatrices.H
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   May 19, 2012, 01:44
Default
  #51
New Member
 
Masashi Arima
Join Date: Apr 2012
Posts: 2
Rep Power: 0
M.Arima is on a distinguished road
Hi,

I tried to set up pimpleFsiFoam under OF 1.5-dev, however, I failed!
Do I need to install OF1.6 to use this solver?

my procdure is to go into the directrory of pimpleFsiFoam and then typed in wclean and wmake.
the error is following!

:~/Desktop/pimpleFsiFoam$ wclean
:~/Desktop/pimpleFsiFoam$ wmake
Making dependency list for source file pimpleFsiFoam.C
Making dependency list for source file tractionDisplacement/tractionDisplacementFvPatchVectorField.C
could not open file turbulenceModel.H for source file pimpleFsiFoam.C
could not open file dynamicFvMesh.H for source file pimpleFsiFoam.C
could not open file createDynamicFvMesh.H for source file pimpleFsiFoam.C
SOURCE=tractionDisplacement/tractionDisplacementFvPatchVectorField.C ; g++ -m32 -Dlinux -DDP -DFOAM_DEV_REVISION_NUMBER=exported -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-40 -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/dynamicMesh/dynamicFvMesh/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/dynamicMesh/dynamicMesh/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/meshTools/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/turbulenceModels/incompressible/turbulenceModel -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/transportModels -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/transportModels/incompressible/singlePhaseTransportModel -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude -DFACE_DECOMP -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/tetDecompositionFiniteElement/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/dynamicMesh/meshMotion/tetDecompositionMotionSolver/lnInclude -IlnInclude -I. -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/OSspecific/Unix/lnInclude -fPIC -pthread -c $SOURCE -o Make/linuxGccDPOpt/tractionDisplacementFvPatchVectorField.o
SOURCE=pimpleFsiFoam.C ; g++ -m32 -Dlinux -DDP -DFOAM_DEV_REVISION_NUMBER=exported -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-40 -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/dynamicMesh/dynamicFvMesh/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/dynamicMesh/dynamicMesh/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/meshTools/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/turbulenceModels/incompressible/turbulenceModel -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/transportModels -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/transportModels/incompressible/singlePhaseTransportModel -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude -DFACE_DECOMP -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/tetDecompositionFiniteElement/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/dynamicMesh/meshMotion/tetDecompositionMotionSolver/lnInclude -IlnInclude -I. -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/OpenFOAM/lnInclude -I/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/OSspecific/Unix/lnInclude -fPIC -pthread -c $SOURCE -o Make/linuxGccDPOpt/pimpleFsiFoam.o
pimpleFsiFoam.C:37:29: error: turbulenceModel.H: No such file or directory
pimpleFsiFoam.C:38:27: error: dynamicFvMesh.H: No such file or directory
pimpleFsiFoam.C:58:36: error: createDynamicFvMesh.H: No such file or directory
tractionDisplacement/tractionDisplacementFvPatchVectorField.C: In member function ‘virtual void Foam::tractionDisplacementFvPatchVectorField::upda teCoeffs()’:
tractionDisplacement/tractionDisplacementFvPatchVectorField.C:190: warning: suggest parentheses around ‘+’ in operand of ‘&’
tractionDisplacement/tractionDisplacementFvPatchVectorField.C:190: warning: suggest parentheses around ‘+’ in operand of ‘&’
tractionDisplacement/tractionDisplacementFvPatchVectorField.C:191: warning: suggest parentheses around ‘+’ in operand of ‘&’
In file included from pimpleFsiFoam.C:59:
/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude/readPIMPLEControls.H: In function ‘int main(int, char**)’:
/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude/readPIMPLEControls.H:1: error: ‘mesh’ was not declared in this scope
In file included from pimpleFsiFoam.C:66:
createFields.H:128: error: ‘incompressible’ was not declared in this scope
createFields.H:128: error: template argument 1 is invalid
createFields.H:129: error: invalid type in declaration before ‘(’ token
createFields.H:130: error: ‘incompressible’ is not a class or namespace
In file included from pimpleFsiFoam.C:67:
createStressMesh.H:12: error: ‘New’ is not a member of ‘Foam::volPointInterpolation’
In file included from solveFluid.H:29,
from weakCoupling.H:3,
from pimpleFsiFoam.C:103:
UEqn.H:5: error: base operand of ‘->’ is not a pointer
In file included from weakCoupling.H:3,
from pimpleFsiFoam.C:103:
solveFluid.H:116: error: base operand of ‘->’ is not a pointer
In file included from weakCoupling.H:7,
from pimpleFsiFoam.C:103:
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:25: warning: suggest parentheses around ‘+’ in operand of ‘&’
In file included from solveFluid.H:29,
from onlyCFD.H:1,
from pimpleFsiFoam.C:109:
UEqn.H:5: error: base operand of ‘->’ is not a pointer
In file included from onlyCFD.H:1,
from pimpleFsiFoam.C:109:
solveFluid.H:116: error: base operand of ‘->’ is not a pointer
In file included from solveFluid.H:29,
from strongCouplingConstant.H:5,
from pimpleFsiFoam.C:115:
UEqn.H:5: error: base operand of ‘->’ is not a pointer
In file included from strongCouplingConstant.H:5,
from pimpleFsiFoam.C:115:
solveFluid.H:116: error: base operand of ‘->’ is not a pointer
In file included from strongCouplingConstant.H:8,
from pimpleFsiFoam.C:115:
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:25: warning: suggest parentheses around ‘+’ in operand of ‘&’
In file included from solveFluid.H:29,
from strongCouplingAitken.H:9,
from pimpleFsiFoam.C:121:
UEqn.H:5: error: base operand of ‘->’ is not a pointer
In file included from strongCouplingAitken.H:9,
from pimpleFsiFoam.C:121:
solveFluid.H:116: error: base operand of ‘->’ is not a pointer
In file included from strongCouplingAitken.H:12,
from pimpleFsiFoam.C:121:
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:25: warning: suggest parentheses around ‘+’ in operand of ‘&’
In file included from solveFluid.H:29,
from strongCouplingAitken.H:36,
from pimpleFsiFoam.C:121:
UEqn.H:5: error: base operand of ‘->’ is not a pointer
In file included from strongCouplingAitken.H:36,
from pimpleFsiFoam.C:121:
solveFluid.H:116: error: base operand of ‘->’ is not a pointer
In file included from strongCouplingAitken.H:39,
from pimpleFsiFoam.C:121:
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:24: warning: suggest parentheses around ‘+’ in operand of ‘&’
solveSolid.H:25: warning: suggest parentheses around ‘+’ in operand of ‘&’
/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude/readPIMPLEControls.H:3: warning: unused variable ‘nOuterCorr’
/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude/readPIMPLEControls.H:4: warning: unused variable ‘nCorr’
/home/c05arima/OpenFOAM/OpenFOAM-1.5-dev/src/finiteVolume/lnInclude/readTimeControls.H:40: warning: unused variable ‘maxCo’
pimpleFsiFoam.C:76: warning: unused variable ‘minGammaAitken’
pimpleFsiFoam.C:77: warning: unused variable ‘maxGammaAitken’

Last edited by M.Arima; May 20, 2012 at 00:46.
M.Arima is offline   Reply With Quote

Old   May 19, 2012, 17:58
Default
  #52
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,258
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings Masashi,

It's probably best that you try to use 1.6-ext instead, since 1.5-dev is AFAIK no longer under development.

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   July 31, 2013, 17:08
Default
  #53
New Member
 
Join Date: Oct 2011
Posts: 25
Rep Power: 5
danny261083 is on a distinguished road
Hi Dr. Habchi,

I read through one of your recent journal publications on the ability of the pimple fsi solver to deal with large structural displacements and rotations using the StVk constitutive model & ALE formulation of the N-S equations for the 2D case. However, I was curious to know whether the solver could be readily extended to a 3D case involving complicated flow and solid domains with unstructured tri + tetrahedral meshes, (e.g. modeling bio-fluid structure interactions).

Regards
danny261083 is offline   Reply With Quote

Old   August 5, 2013, 04:10
Default
  #54
New Member
 
charbel's Avatar
 
Charbel
Join Date: Jan 2011
Location: Lebanon
Posts: 27
Rep Power: 6
charbel is on a distinguished road
Quote:
Originally Posted by danny261083 View Post
Hi Dr. Habchi,

I read through one of your recent journal publications on the ability of the pimple fsi solver to deal with large structural displacements and rotations using the StVk constitutive model & ALE formulation of the N-S equations for the 2D case. However, I was curious to know whether the solver could be readily extended to a 3D case involving complicated flow and solid domains with unstructured tri + tetrahedral meshes, (e.g. modeling bio-fluid structure interactions).

Regards
Hi
In fact mainly the solver can solve 3D case but the problem is that we cannot run parallel simulation.
Meanwhile you can find in the thread below new solvers including FSI and which may run in parallel...
Solid Mechanics Solvers added to OpenFOAM Extend
Regards
arashfluid likes this.
__________________
--
Charbel
charbel is offline   Reply With Quote

Old   October 7, 2013, 03:30
Default
  #55
Member
 
arash mahboubi
Join Date: Jun 2013
Location: Iran
Posts: 37
Rep Power: 4
arashfluid is on a distinguished road
Send a message via Yahoo to arashfluid
Hi Charbel
I set up a FSI problem for compressible flow over naca0012 airfoil with pitching motion. So i created a compressible solver by combine pimplefsifoam and sonicfoam. Solver has been properly compiled.

I create 2 meshes in gambit, one for the fluid and another one with the airfoil only as my solidMesh. I convert them using fluentmeshToFoam. When I run the case I get the following error:


--> FOAM FATAL ERROR:
Different dimensions for =
dimensions : [0 0 1 0 0 0 0] = [-1 3 1 0 0 0 0]


From function dimensionSet:perator=(const dimensionSet& ds) const
in file dimensionSet/dimensionSet.C at line 143.

FOAM aborting

what's problem?

Last edited by arashfluid; November 4, 2013 at 04:11.
arashfluid is offline   Reply With Quote

Old   October 9, 2013, 02:28
Default
  #56
New Member
 
Vincent
Join Date: Aug 2013
Location: Singapore
Posts: 7
Rep Power: 3
VincentChai is on a distinguished road
Hi Charbel,

I am also a new user to OpenFOAM. May i know have you validate the icoFsiFoam for your application? Thanks.

Regards,
Vincent
VincentChai is offline   Reply With Quote

Old   October 14, 2014, 10:51
Default
  #57
New Member
 
Join Date: Oct 2011
Posts: 25
Rep Power: 5
danny261083 is on a distinguished road
Quote:
Originally Posted by waku2005 View Post
Hi Charbel,

Thanks a lot for sharing your FSI solver, pimpleFsiFoam.
I also successfully compiled and ran the test case attached with the previous version in the thread #14.

I'd like to ask you whether this solver can be applied to lower rhof/rhos cases.
My application is deformation of beam like structure (a silicon rubber block) in a water with retatively slow flow velocity of 0.1 m/s.
rho of fluid (water) is around 1000 kg/m3 and silicon rubber has almost the same, thus rhof/rhos is nearly unity.

Even I changed fluid density from 1 to 1000 kg/m3, and adjusted time step, solver diverged in several time steps.

Any comment will be highly appreciated.
Best regards,
waku
Hi waku,

I would like to know the correct implementation of this solver in OpenFOAM 1.6ext. Also, I am not sure that I have complied the updated solver (and test case). I would like to know whether you could post a link of the same.

Thanks in advance.
danny261083 is offline   Reply With Quote

Old   February 17, 2015, 13:05
Default
  #58
Senior Member
 
ok
Join Date: Oct 2013
Posts: 256
Rep Power: 4
Maimouna is on a distinguished road
Dear All who posted this thread,

many thanks for this really useful thread. I'll try to post my planning in the next period of time to get help and advice from your side.

Now, I'm working in vortex-induced vibration of an elastic circular cylinder using OpenFOAM. This is FSI problem including fluid and solid (elastic cylinder). I installed FSI solver in foam-extend-3.1. But I'm still struggling in the steps how to start working with this case in FSI solver. I already finished the VIV case for rigid cylinder. Now, I'm looking to get the simulation of an elastic one. For the elastic cylinder, the cylinder will change its shape because of some fluid forces affect on it. The equation that I have to solve here is the dynamic elasticity equation. Regarding my problem, I have some questions

1. Which solver is the best for my case, icoFsiFoam or icoFsiElasticNonLinULSolidFoam solver in solidMechanics tutorial?

Do you have any idea which case should follow to help me in my case problem? I'm thinking about flappingConsoleSmall case, would that help?


I'm looking forward to get your advice. I'll be really happy to work with you.

What about mesh generation? Shall I generate mesh for the plate and another one for cylinder separately? Note that, I generated mesh for the VIV for the rigid cylinder case by Gmsh? Any idea please.

Any answer would be really appreciated and many thanks in advanced.

Sorry, I posted my problem in Elastic or Plastic Deformations in OpenFOAM
because I don't know the correct place to post my problem.
Maimouna is offline   Reply With Quote

Reply

Tags
elastic, icofsifoam, strongly weakly coupling

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
NASA Rotor 37 & Pressure-Based Coupled Solver Razvan FLUENT 6 June 11, 2015 03:58
coupled solver wont work in star ccm+ richie CD-adapco 5 November 4, 2008 05:51
coupled solver in FLUENT ztdep FLUENT 1 September 14, 2006 22:34
The seperated solver and coupled one summer FLUENT 0 July 5, 2006 15:06
can use coupled solver to calculate combustion cfdfans FLUENT 0 November 25, 2003 05:48


All times are GMT -4. The time now is 00:21.