HronTurekFsi3 foam-extend3.1

wyldckat · April 17, 2016, 20:42

Quick answer: It's now 01h40 (AM). I'm nearly exhausted and I need to get up in less than 7h from now.

@Maimouna: Attached is the file "beamInCross_sculptored_v1.tar.gz", which is roughly based on the original tutorial case "HronTurekFsi3".

The new case uses a different meshing strategy, based on a single block in "blockMeshDict" in the fluid region and then it refines multiple times over the region of interest.
You will now have to study the by yourself or with someone else's help, because I probably won't have any time left this week.

Details to keep in mind:

This case is not yet defined to use the properties of the plate from your article. I didn't have the time/patience to do that.
The plate is thicker than you wanted. It's roughly "0.000312 m" thick, because the meshing strategy I used is on one hand simpler to create, but on the other hand it can be bothersome to do when I'm in a hurry.
The folder "fluid/0.org" is now the real initial time step folder. Configure things there if you need to. The folder "fluid/0" is only temporary and it will be automatically deleted.
To run the case, simply do these steps:
Code:
```
cd fluid
./Allrun
```
To restart the case, run from within the folder "fluid" this command:
Code:
```
./Allclean
```
The mesh has over 1 million cells. The first time step takes 173s to run in my machine and roughly 60s on the next time steps.
- This means that for running 6 seconds of simulation time, at 0.001s per time step, equates to: 6/0.001*60 = 360000 seconds to run, which is roughly... 100 hours... 4.17 days to run.

I do not have enough experience with fsiFoam to know how this case can be executed in parallel.

Good luck!

Woj3x · April 18, 2016, 18:23

Dear Bruno
Is fsiFoam prepared to use variable time step or is it required to modify its source code?
I tried to run case like it's showed in this link (I didn't modify fsiFoam source code):
http://www.idurun.com/2008/06/13/run...p-in-openfoam/
but I don't see any effect

My controlDict:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     3.0                                |
|   \\  /    A nd           | Web:         http://www.extend-project.de       |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application fsiFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         6e-2;

deltaT          1e-7;//2e-7;//1e-7;
/*new_change*/
maxCo 0.2; // Or other Courant number you wish

adjustTimeStep yes; // Or no 

maxDeltaT 1e-5; // Maximum deltaT in seconds

writeControl adjustableRunTime;

/*new_change*/

writeControl    runTime;

writeInterval   2e-7;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression uncompressed;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

adjustTimeStep  no;

maxCo           0.2;

functions
(
   beamReport
   {
       type pointHistory;
       functionObjectLibs
       (
         "libpointHistory.so"
       );

       refHistoryPoint (0.0001 0.012 0.012);

       region solid;
   }

   setInletVelocity
   {
       type setInletVelocity;
       functionObjectLibs
       (
         "libsetInletVelocity.so"
       );
   }

   forces
   {
       type                forces;
       functionObjectLibs  ( "libforces.so" );
       outputControl       timeStep;
       outputInterval      1;
       patches             (interface);
       pName               p;
       UName               U;
       rhoName             rhoInf;
       log                 true;
       rhoInf              1000;
       CofR                (0.5 0.1 0);
   }

);

// ************************************************************************* //

Thank you in advance
Wojciech

wyldckat · April 19, 2016, 18:45

@Maimouna: Is it working?

@Woj3x: If it works, it will only work with this option:

Code:

adjustTimeStep yes;

Let me check the code... nope, it's fixed time step. There is no time step update in the code "fsiFoam.C".

Maimouna · April 25, 2016, 07:19

Hi Bruno,

Code:

@Maimouna: Is it working?

it is not yet, I got problem in my old Ubuntu machine. Now, I got a new one. I installed foam-extend-3.1 in the new machine and I'm on the stage to install and compile FluidStructureInteraction package. Not finished yet.

I faced some difficulties on that. It was installed and compiled in the old machine. Could you help me with the steps how to install and compile FSI package for foam-extend-3.1?

Once finished this, I will test the case.

I have some questions please,
1. What did you call your meshing way?
2. What are the differences between blockMeshDict and blockMeshDict.first for the fluid part?
3. What is Reynold's number in your case?

Lots of thanks with kind regards

Maimouna

Maimouna · April 26, 2016, 12:15

Dear Bruno,

following the steps given to install FSI package for foam-extend-3.1, I got the following error

Code:

make: *** [/home/owner/foam/owner-3.1/lib/linux64GccDPOpt/libhronTurekReport.so] Error 1

as it is shown in the screen results below

Code:

owner@owner-SATELLITE-L50-C:~$ fe31
owner@owner-SATELLITE-L50-C:~$ cd foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3$ chmod 755 *
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3$ sed -i s/tcsh/sh/g *Links
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3$ ./removeSerialLinks fluid solid
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3$ ./makeSerialLinks fluid solid
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3$ cd fluid/
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid$ ./Allclean 
Cleaning /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid case
wclean ../hronTurekReport
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid$ ./Allrun 
Running blockMesh on /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid
Running setSet on /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid
Running setsToZones on /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid
Running blockMesh on /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid
Running setSet on /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid
Running setsToZones on /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid
Making dependency list for source file hronTurekReport.C
SOURCE=hronTurekReport.C ;  g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3  -DNoRepository -ftemplate-depth-200 -I../../../../src/fluidStructureInteraction/lnInclude     -I/home/owner/foam/foam-extend-3.1/src/finiteVolume/lnInclude     -I/home/owner/foam/foam-extend-3.1/src/finiteArea/lnInclude     -I/home/owner/foam/foam-extend-3.1/src/meshTools/lnInclude     -I/home/owner/foam/foam-extend-3.1/src/dynamicMesh/dynamicFvMesh/lnInclude     -I/home/owner/foam/foam-extend-3.1/src/triSurface/lnInclude     -I/home/owner/foam/foam-extend-3.1/src/lagrangian/basic/lnInclude     -I/home/owner/foam/foam-extend-3.1/src/solidModels/lnInclude -IlnInclude -I. -I/home/owner/foam/foam-extend-3.1/src/foam/lnInclude -I/home/owner/foam/foam-extend-3.1/src/OSspecific/POSIX/lnInclude   -fPIC -c $SOURCE -o Make/linux64GccDPOpt/hronTurekReport.o
/usr/bin/ld: cannot find -lfluidStructureInteraction
collect2: error: ld returned 1 exit status
make: *** [/home/owner/foam/owner-3.1/lib/linux64GccDPOpt/libhronTurekReport.so] Error 1
Running fsiFoam on /home/owner/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid
owner@owner-SATELLITE-L50-C:~/foam/owner-3.1/FluidStructureInteraction/run/fsiFoam/HronTurekFsi3/fluid

How could solve that? What's missing?

I was working well in my old machine, but unfortunately I didn't save the steps and I'm not sure which website I was using to install FSI package. Could you please help me how to reinstall that package for foam-extend-3.1. I'm not sure, if I sould install it within extend-bazaar directory or within a new directory called owner-3.1 (for example).

I tried the steps given in

HTML Code:

http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2014/ThomasVyzikas/interFlowReport.pdf

and

HTML Code:

http://openfoamwiki.net/index.php/Extend-bazaar/Toolkits/Fluid-structure_interaction

Waiting to get your help. Many thanks in advanced.

Maimouna

Maimouna · May 1, 2016, 05:34

Dear Bruno,

I knew that it is not the best place to post my question anout installation of FSI package for extend-foam-3.1. Could you please help me how to install that package. And please answer my questions in the last post.

The installation takes for more than 5 days and it's still not finish, but wasn't that long time in my old machine. May be I followed wrong steps.

Waiting to get your answer please.

Warm regards

Maimouna

wyldckat · May 1, 2016, 19:53

Quick answers:

Quote:

Originally Posted by Maimouna

Hi Bruno,

Code:

@Maimouna: Is it working?

it is not yet, I got problem in my old Ubuntu machine. Now, I got a new one. I installed foam-extend-3.1 in the new machine and I'm on the stage to install and compile FluidStructureInteraction package. Not finished yet.

Unfortunately, this happens to perhaps 1-15% of students who need to finish assignments... it's when the machine is needed the most, that it breaks down...

Quote:

Originally Posted by Maimouna

I have some questions please,
1. What did you call your meshing way?
2. What are the differences between blockMeshDict and blockMeshDict.first for the fluid part?
3. What is Reynold's number in your case?

There is no specific name for it. It's just a matter of relying on the tools we have in OpenFOAM/foam-extend to create the mesh we need. In this case, create a base mesh, refine it to the point we need to and then remove the cells we don't want in our mesh.
"blockMeshDict.first" isn't needed. It's something from some other case, I don't remember which one.
I have absolutely no idea. I didn't do any calculation myself, I only set-up the case so that it would mesh and run. You will need to check all details and fluid/solid properties.

As for building the toolkit, follow the instructions given here: http://openfoamwiki.net/index.php/Ex...oam-extend_3.1

In addition, I need to know which installation instructions you followed for installing foam-extend 3.1 and which Linux Distribution and version you are using.

Maimouna · May 2, 2016, 17:41

Quote:

Originally Posted by wyldckat

Quick answers:

Unfortunately, this happens to perhaps 1-15% of students who need to finish assignments... it's when the machine is needed the most, that it breaks down...

There is no specific name for it. It's just a matter of relying on the tools we have in OpenFOAM/foam-extend to create the mesh we need. In this case, create a base mesh, refine it to the point we need to and then remove the cells we don't want in our mesh.
"blockMeshDict.first" isn't needed. It's something from some other case, I don't remember which one.
I have absolutely no idea. I didn't do any calculation myself, I only set-up the case so that it would mesh and run. You will need to check all details and fluid/solid properties.

As for building the toolkit, follow the instructions given here: http://openfoamwiki.net/index.php/Ex...oam-extend_3.1

In addition, I need to know which installation instructions you followed for installing foam-extend 3.1 and which Linux Distribution and version you are using.

Hi Bruno,

to inform you FSI package was installed, thank you. But the case doesn't working yet. What do you recommended to change in order to make it run? I tried to attache you the case but it gives me this message

Quote:

Your submission could not be processed because the token has expired.

Please push the back button and reload the previous window.

Which token? How could solve that?

Many thaks.

Maimouna

wyldckat · May 3, 2016, 16:31

Quote:

Originally Posted by Maimouna

Which token? How could solve that?

Quick answer: It says right after that line what you should do:

Quote:

Please push the back button and reload the previous window.

In other words, refresh the page, or close this window and open a new window.

Woj3x · May 8, 2016, 11:11

Dear Bruno
I have a strange problem with my case.
I run it with time step 1e-7 s and Currant number became too high and at Time = 2.8e-05, iteration: 14 I've get:

Code:

 Floating point exception (core dumped)

So I decided to run this case again from Time 2.7e-05 (my last saved step) with new time step: 1e-8 s but Currant number is even bigger right now. How is it possible and how can I fix this? Is it connected with fvSolution's tolerances (I didn't change it) or my time step is still to big, or my mesh isn't good eanough?

Bellow I attachat few selected results from my log.fsiFoam for time step 1e-7 s and 1e-8 s and link for my case and full log files. I will be really gratefull if you could look into this case.
Thank you in advance

Time step 1e-7:

Code:

Time = 2.71e-05, iteration: 20
Courant Number mean: 0.000764766 max: 0.109055 velocity magnitude: 397.31

Time = 2.72e-05, iteration: 20
Courant Number mean: 0.000871482 max: 0.131472 velocity magnitude: 471.698

Time = 2.73e-05, iteration: 20
Courant Number mean: 0.00101104 max: 0.161952 velocity magnitude: 570.313

Time = 2.74e-05, iteration: 20
Courant Number mean: 0.00119925 max: 0.205171 velocity magnitude: 705.737

Time = 2.75e-05, iteration: 20
Courant Number mean: 0.00146355 max: 0.270027 velocity magnitude: 900.267

Time = 2.76e-05, iteration: 20
Courant Number mean: 0.00185561 max: 0.375597 velocity magnitude: 1196.43

Time = 2.77e-05, iteration: 20
Courant Number mean: 0.00248462 max: 0.568582 velocity magnitude: 1681.5

Time = 2.78e-05, iteration: 20
Courant Number mean: 0.00362456 max: 0.999457 velocity magnitude: 2526.32

Time = 2.79e-05, iteration: 20
Courant Number mean: 0.0062092 max: 2.87982 velocity magnitude: 3859.73

Time = 2.8e-05, iteration: 14
Courant Number mean: 0.00902291 max: 36.0155 velocity magnitude: 6439.49
Floating point exception (core dumped)

Time step 1e-8:

Code:

Time = 2.711e-05, iteration: 20
Courant Number mean: 8.68819e-05 max: 0.0127651 velocity magnitude: 464.691

Time = 2.72e-05, iteration: 20
Courant Number mean: 0.000345062 max: 0.0626587 velocity magnitude: 2186.54

Time = 2.722e-05, iteration: 20
Courant Number mean: 0.000505583 max: 0.0960744 velocity magnitude: 3283.61

Time = 2.725e-05, iteration: 20
Courant Number mean: 0.00100025 max: 0.211915 velocity magnitude: 6809.41

Time = 2.727e-05, iteration: 20
Courant Number mean: 0.00181839 max: 0.445003 velocity magnitude: 12879.9

Time = 2.73e-05, iteration: 20
Courant Number mean: 0.00966776 max: 21.7874 velocity magnitude: 68842.3

My case:

https://www.dropbox.com/s/ujbyns2f6s...1a.tar.gz?dl=0

Full logs:
https://www.dropbox.com/sh/il4cgm0vo...98-z65CQa?dl=0

Could you also tell me how to interpret velocity magnitude from log file? Is it total velocity of all fluid? I can't get why is it so big and what it truly refer to.

Maimouna · May 13, 2016, 07:16

Quote:

Originally Posted by Woj3x

Dear Bruno
I have a strange problem with my case.
I run it with time step 1e-7 s and Currant number became too high and at Time = 2.8e-05, iteration: 14 I've get:

Code:

 Floating point exception (core dumped)

So I decided to run this case again from Time 2.7e-05 (my last saved step) with new time step: 1e-8 s but Currant number is even bigger right now. How is it possible and how can I fix this? Is it connected with fvSolution's tolerances (I didn't change it) or my time step is still to big, or my mesh isn't good eanough?

Bellow I attachat few selected results from my log.fsiFoam for time step 1e-7 s and 1e-8 s and link for my case and full log files. I will be really gratefull if you could look into this case.
Thank you in advance

Time step 1e-7:

Code:

Time = 2.71e-05, iteration: 20
Courant Number mean: 0.000764766 max: 0.109055 velocity magnitude: 397.31

Time = 2.72e-05, iteration: 20
Courant Number mean: 0.000871482 max: 0.131472 velocity magnitude: 471.698

Time = 2.73e-05, iteration: 20
Courant Number mean: 0.00101104 max: 0.161952 velocity magnitude: 570.313

Time = 2.74e-05, iteration: 20
Courant Number mean: 0.00119925 max: 0.205171 velocity magnitude: 705.737

Time = 2.75e-05, iteration: 20
Courant Number mean: 0.00146355 max: 0.270027 velocity magnitude: 900.267

Time = 2.76e-05, iteration: 20
Courant Number mean: 0.00185561 max: 0.375597 velocity magnitude: 1196.43

Time = 2.77e-05, iteration: 20
Courant Number mean: 0.00248462 max: 0.568582 velocity magnitude: 1681.5

Time = 2.78e-05, iteration: 20
Courant Number mean: 0.00362456 max: 0.999457 velocity magnitude: 2526.32

Time = 2.79e-05, iteration: 20
Courant Number mean: 0.0062092 max: 2.87982 velocity magnitude: 3859.73

Time = 2.8e-05, iteration: 14
Courant Number mean: 0.00902291 max: 36.0155 velocity magnitude: 6439.49
Floating point exception (core dumped)

Time step 1e-8:

Code:

Time = 2.711e-05, iteration: 20
Courant Number mean: 8.68819e-05 max: 0.0127651 velocity magnitude: 464.691

Time = 2.72e-05, iteration: 20
Courant Number mean: 0.000345062 max: 0.0626587 velocity magnitude: 2186.54

Time = 2.722e-05, iteration: 20
Courant Number mean: 0.000505583 max: 0.0960744 velocity magnitude: 3283.61

Time = 2.725e-05, iteration: 20
Courant Number mean: 0.00100025 max: 0.211915 velocity magnitude: 6809.41

Time = 2.727e-05, iteration: 20
Courant Number mean: 0.00181839 max: 0.445003 velocity magnitude: 12879.9

Time = 2.73e-05, iteration: 20
Courant Number mean: 0.00966776 max: 21.7874 velocity magnitude: 68842.3

My case:

https://www.dropbox.com/s/ujbyns2f6s...1a.tar.gz?dl=0

Full logs:
https://www.dropbox.com/sh/il4cgm0vo...98-z65CQa?dl=0

Could you also tell me how to interpret velocity magnitude from log file? Is it total velocity of all fluid? I can't get why is it so big and what it truly refer to.

Hi Wojciech,

could you please explain more about your case might I could help.

Kind regards

Maimouna

KarthickRajkumar · November 13, 2016, 16:08

Dear all,

I am working on a project involving fsifoam solver in foam-extend 3.2. I copied the BeamInCrossFlow tutorial and tried it for my application. I get the following error everytime. I was able to rectify all my previous errors using this thread, but this error I couldn't.

Could anyone please help me find to prevent this error?
Should you need more details, pls ask me!

Thanks!!

Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Reading stress mesh
Creating traction displacement boundary conditions
Creating fixed displacement boundary condition
Creating fixed displacement boundary condition
Selecting rheology model linearElastic
Creating constitutive model
Selecting coupling scheme Aitken

Starting time loop

Creating beamReport function object.
History point ID: 139
History point coordinates: (0.15 0.045 0)
Reference point coordinates: (0.45 0.15 -0.15)
Seting inlet velocity
Umax = 0.00489435
Time = 0.1

Create extended GGI zone-to-zone interpolator
Checking fluid-to-solid face interpolator
Fluid-to-solid face interpolation error: 0.0126103
Checking solid-to-fluid point interpolator (GGI)
Extended GGI, master point distance, max: 1e+15, avg: 2.30143e+14, min: -0.002
Extended GGI, master point orientation (<0), max: 1, min: -1, nIncorrectPoints: 638/982

--> FOAM FATAL ERROR:
Master point addressing is not correct

From function ExtendedGGIInterpolation::masterToSlavePointInterp olate(const Field<Type> pf)
in file ./numerics/ggi/ExtendedGGIInterpolation/ExtendedGGIInterpolation.C at line 703.

FOAM aborting
Aborted (core dumped)

Woj3x · November 13, 2016, 18:17

Dear KarthickRajkumar

I'm not 100% sure about it but:

Quote:

--> FOAM FATAL ERROR:
Master point addressing is not correct

It means (if I good remember) that solid and fluid geometry isn't compatible enough (I mean there is a gap between solid and fluid interface). So IMO you need to correct your geometry and mesh.

Kind regards

KarthickRajkumar · November 16, 2016, 16:02

Dear Wojciech,

Thank you for the prompt reply. I somehow managed to fix the error by correcting the mesh and the geometry as you suggested.

Thank you very much. But I face following problems after Iteration starts.

1. When I define startTime 0 with DeltaT 0.00001, it says floating point exception (core dumped)

Time = 1e-05, iteration: 1
Current fsi under-relaxation factor: 0.4
Creating accumulated fluid interface displacement
Maximal accumulated displacement of interface points: 0
GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1
Evolving flow model: consistentIcoFlow
Courant Number mean: 2.14175e-07 max: 6e-06 velocity magnitude: 0.0001
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.31844e-07, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 0.999999, Final residual = 3.66704e-07, No Iterations 4
GAMG: Solving for p, Initial residual = 1, Final residual = 5.30983e-07, No Iterations 7
GAMG: Solving for p, Initial residual = 0.0114117, Final residual = 9.28912e-07, No Iterations 67
time step continuity errors : sum local = 3.64821e-14, global = -3.21772e-24
GAMG: Solving for p, Initial residual = 0.266161, Final residual = 8.08369e-07, No Iterations 67
GAMG: Solving for p, Initial residual = 8.50573e-07, Final residual = 5.47487e-07, No Iterations 2
time step continuity errors : sum local = 2.46898e-14, global = 9.27448e-27
GAMG: Solving for p, Initial residual = 0.288257, Final residual = 8.85958e-07, No Iterations 51
GAMG: Solving for p, Initial residual = 9.1568e-07, Final residual = 6.94922e-07, No Iterations 2
time step continuity errors : sum local = 3.34729e-14, global = -3.38972e-24
Setting traction on solid patch
Total force (fluid) = (0.00870746 2.02047e-05 1.30104e-13)
Total force (solid) = (-0.00870746 -2.02047e-05 -1.30104e-13)
Evolving stress model
-1.65018, 1.0115
DICPCG: Solving for D, Initial residula = 1, Final residual = 0.000587252, No outer iterations = 108
Max relative residual = 24.1206, Relative residual = 0.000514779, enforceLinear = 1
Current fsi relative residual norm: 1

Time = 1e-05, iteration: 2
Current fsi under-relaxation factor: 0.4
Maximal accumulated displacement of interface points: 2.95406e+19
GAMG: Solving for cellMotionUx, Initial residual = 1, Final residual = 0.000189476, No Iterations 2
GAMG: Solving for cellMotionUy, Initial residual = 1, Final residual = 0.000189394, No Iterations 2
GAMG: Solving for cellMotionUx, Initial residual = 9.24019e-05, Final residual = 4.92158e-07, No Iterations 4
GAMG: Solving for cellMotionUy, Initial residual = 9.23671e-05, Final residual = 4.91437e-07, No Iterations 4
Evolving flow model: consistentIcoFlow
Floating point exception (core dumped)

2. When I define startTime 0 with DeltaT 0.0000001, it says continuity errors cannot be removed by adjusting the outflow.

Time = 1e-07, iteration: 1
Current fsi under-relaxation factor: 0.4
Creating accumulated fluid interface displacement
Maximal accumulated displacement of interface points: 0
GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1
GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1
Evolving flow model: consistentIcoFlow
Courant Number mean: 2.14175e-09 max: 6e-08 velocity magnitude: 0.0001
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 8.64249e-10, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.999999, Final residual = 7.81295e-07, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 9.35145e-07, No Iterations 20
GAMG: Solving for p, Initial residual = 0.000371717, Final residual = 9.18606e-07, No Iterations 20
time step continuity errors : sum local = 3.90917e-16, global = -3.72126e-26
GAMG: Solving for p, Initial residual = 0.184356, Final residual = 8.65415e-07, No Iterations 32
GAMG: Solving for p, Initial residual = 6.96729e-07, Final residual = 3.55687e-07, No Iterations 2
time step continuity errors : sum local = 2.19023e-16, global = 1.67978e-30

--> FOAM FATAL ERROR:
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Specified mass inflow : 1.29924e-07
Specified mass outflow : 9.35926e-20
Difference : 1.29924e-07
Adjustable mass outflow : 4.03464e-16

From function adjustPhi
(
surfaceScalarField& phi,
const volVectorField& U,
const volScalarField& p
)
in file cfdTools/general/adjustPhi/adjustPhi.C at line 150.

FOAM exiting

Could you please explain whats happening here!!

Thanks!!

Woj3x · November 16, 2016, 23:45

Dear KarthickRajkumar

First of all it would be nice if you could tell me if is it still a BeamInCrossFlow case?
I see that it is 2D case so I am not sure about this.
How much did you change this case?
Without this informations it's really hard to tell whats wrong with it.
IMO it can be boundary conditions, material parameters or even script problem (ex. velocityInlet from Allrun).

Wojciech

KarthickRajkumar · November 17, 2016, 03:02

Dear Wojciech,

I have attached my case folder here. In fluid/blockMeshdict you will be able to see the 2D design I am simulating (Roughly drawn in the comments section).

About the changes I made in the case files:
BlockMeshDict in fluid/constant/polymesh - defined fluid geometry here
BlockMeshDict in solid/constant/polymesh - defined solid geometry here
And controlDict files, I didn't change any of the other files

You can just run Allrun script from fluid folder and see the Mesh. You will be also able to see the errors I described in the previous post.

I just wanted to simulate the deflection of the lip part in the design (beam) due to the flow.

Since I am new to fsi and also to OF, it would also be great if you could describe in short how the work between solid and fluid folder happens!! Or could you please share if there's a link where I can to get learn about fsi simulations in OF from which I can avoid obvious mistakes!!

Thanks!!

KarthickRajkumar · November 24, 2016, 04:27

Quote:

Originally Posted by Maimouna

Dear Bruno,

many thanks for that well explanation in the last post.

As a good news after some changes in the last case ''myCylinderCase-Fsi'' the beam is moving as shown in the attached screenshot. (I tried to attach animation but I don't know which extension is accepted .avi not accepted). The new case with new changes is attached too to have a look.

Dear Maimouna,

I am doing the similar simulation using fsifoam solver. Similar setup as HronTurekfsi tutorial, but with a different geometry. I followed the steps suggested by Bruno until here. I have the same problem as yours. The simulation works fine but there is no change in the movement of the plate like you had.

You had said you did some changes in the set up and beam was moving. Could you please write me about what caused that problem and how exactly did yu fix it??

I see you didn't change the inletVelocity as Parabolic or gradual increase at the inlet as Bruno suggested. I would like to know how did you fix that error so that I may give a thought about it in my case.

Thanks,
Karthick

Maimouna · November 24, 2016, 04:52

Quote:

Originally Posted by KarthickRajkumar

Dear Maimouna,

I am doing the similar simulation using fsifoam solver. Similar setup as HronTurekfsi tutorial, but with a different geometry. I followed the steps suggested by Bruno until here. I have the same problem as yours. The simulation works fine but there is no change in the movement of the plate like you had.

You had said you did some changes in the set up and beam was moving. Could you please write me about what caused that problem and how exactly did yu fix it??

I see you didn't change the inletVelocity as Parabolic or gradual increase at the inlet as Bruno suggested. I would like to know how did you fix that error so that I may give a thought about it in my case.

Thanks,
Karthick

Hi Karthick,

sure I will help, could you please provide me by your case to have a look.

Which extend version yoy use?
What is your geometry and BCs?

Kind regards
Maimouna

Woj3x · November 24, 2016, 05:39

Hi Karthick,
Sorry I did not answer but I had a lots to do in last week.
Could you look into your log.fsifoam file and tell us how big is Courant number in your results?

KarthickRajkumar · November 24, 2016, 07:17

Hai Maimouna and Wojciech,

Thank you for the prompt reply.

I use foam-extend 3.2 for my simulation.
My geometry is of a musical instrument. You can see the geometry I have roughly drawn under comments in the fluid/constant/polymesh/blockMeshdict for better understanding.

I have also specified in the diagram what is cylinder and plate in my case.
My aim is to see the deflection of the plate due to the flow.

Initially I tried with BeamInCross flow and then I realised it has got some other situation for the inlet velcoity.
So now I am using HronTurek case, which has similar setup as mine and I think is the better option for my simulation.

So about the BC as you asked:

U: Only change I made is the inlet velocity set to (-10 0 0) - negative since the flow is in the negative x direction.
P: No changes
pointMotionU: No changes

Additional Info

I ran the HronTurek case with the inlet as outlet and outlet as inlet, so that the flow comes from negative x direction and to ensure that the simulation runs in such a case. It does run and I see the deflection.

Things which I doubt:

1. I ran the simulation only for 2 seconds since it takes more than 10 hours due to fine mesh. So should I run the simulation for more seconds to check for the deflection?

2. The difference between the tutorial and my set up is that, cylinder is not attached to any wall in tutorial but in my case the cylinder is attached to a wall. I believe it should not affect the simulation set up. And may be because it is fixed the deflection can be zero until 2 seconds. What do you think in this case?

3. Deflection is 0 for every time step until 2 seconds and just in case you need the force in force.dat file, here it is:

Force at the end of 2 seconds is
17.0856 527.416 1

@ Wojciech: About the courant number: I didn't save the log file. Didn't know I have to. But my simulation blows up at 2 seconds due to high courant number. (I ran it with 0.00001 time step). I thought of running the simulation in the cluster wth 0.0000001 time step if it shows some deflection.
These are all my thoughts up till now. Could you please help me at it. Any comments or suggestion would really help me.

You can find the my latest case file here. Please avoid any other files I uploaded before.

Thanks,
Karthick

April 19, 2016, 18:45		#83
wyldckat Retired Super Moderator Bruno Santos Join Date: Mar 2009 Location: Lisbon, Portugal Posts: 10,975 Blog Entries: 45 Rep Power: 128	@Maimouna: Is it working? @Woj3x: If it works, it will only work with this option: Code: adjustTimeStep yes; Let me check the code... nope, it's fixed time step. There is no time step update in the code "fsiFoam.C". Woj3x likes this.

April 25, 2016, 07:19		#84
Maimouna Senior Member ok Join Date: Oct 2013 Posts: 346 Rep Power: 13	Hi Bruno, Code: @Maimouna: Is it working? it is not yet, I got problem in my old Ubuntu machine. Now, I got a new one. I installed foam-extend-3.1 in the new machine and I'm on the stage to install and compile FluidStructureInteraction package. Not finished yet. I faced some difficulties on that. It was installed and compiled in the old machine. Could you help me with the steps how to install and compile FSI package for foam-extend-3.1? Once finished this, I will test the case. I have some questions please, 1. What did you call your meshing way? 2. What are the differences between blockMeshDict and blockMeshDict.first for the fluid part? 3. What is Reynold's number in your case? Lots of thanks with kind regards Maimouna

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Other] mesh airfoil NACA0012	anand_30	OpenFOAM Meshing & Mesh Conversion	13	March 7, 2022 17:22
[blockMesh] error message with modeling a cube with a hold at the center	hsingtzu	OpenFOAM Meshing & Mesh Conversion	2	March 14, 2012 09:56
[blockMesh] BlockMesh FOAM warning	gaottino	OpenFOAM Meshing & Mesh Conversion	7	July 19, 2010 14:11
[blockMesh] Axisymmetrical mesh	Rasmus Gjesing (Gjesing)	OpenFOAM Meshing & Mesh Conversion	10	April 2, 2007 14:00
[Gmsh] Import gmsh msh to Foam	adorean	OpenFOAM Meshing & Mesh Conversion	24	April 27, 2005 08:19

May 1, 2016, 05:34		#86
Maimouna Senior Member ok Join Date: Oct 2013 Posts: 346 Rep Power: 13	Dear Bruno, I knew that it is not the best place to post my question anout installation of FSI package for extend-foam-3.1. Could you please help me how to install that package. And please answer my questions in the last post. The installation takes for more than 5 days and it's still not finish, but wasn't that long time in my old machine. May be I followed wrong steps. Waiting to get your answer please. Warm regards Maimouna

November 13, 2016, 16:08		#92
KarthickRajkumar New Member Karthick Join Date: Oct 2016 Location: Munich Posts: 18 Rep Power: 9	Dear all, I am working on a project involving fsifoam solver in foam-extend 3.2. I copied the BeamInCrossFlow tutorial and tried it for my application. I get the following error everytime. I was able to rectify all my previous errors using this thread, but this error I couldn't. Could anyone please help me find to prevent this error? Should you need more details, pls ask me! Thanks!! Create time Create dynamic mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: velocityLaplacian Selecting motion diffusion: quadratic Selecting motion diffusion: inverseDistance Reading stress mesh Creating traction displacement boundary conditions Creating fixed displacement boundary condition Creating fixed displacement boundary condition Selecting rheology model linearElastic Creating constitutive model Selecting coupling scheme Aitken Starting time loop Creating beamReport function object. History point ID: 139 History point coordinates: (0.15 0.045 0) Reference point coordinates: (0.45 0.15 -0.15) Seting inlet velocity Umax = 0.00489435 Time = 0.1 Create extended GGI zone-to-zone interpolator Checking fluid-to-solid face interpolator Fluid-to-solid face interpolation error: 0.0126103 Checking solid-to-fluid point interpolator (GGI) Extended GGI, master point distance, max: 1e+15, avg: 2.30143e+14, min: -0.002 Extended GGI, master point orientation (<0), max: 1, min: -1, nIncorrectPoints: 638/982 --> FOAM FATAL ERROR: Master point addressing is not correct From function ExtendedGGIInterpolation::masterToSlavePointInterp olate(const Field<Type> pf) in file ./numerics/ggi/ExtendedGGIInterpolation/ExtendedGGIInterpolation.C at line 703. FOAM aborting Aborted (core dumped)

November 16, 2016, 16:02		#94
KarthickRajkumar New Member Karthick Join Date: Oct 2016 Location: Munich Posts: 18 Rep Power: 9	Dear Wojciech, Thank you for the prompt reply. I somehow managed to fix the error by correcting the mesh and the geometry as you suggested. Thank you very much. But I face following problems after Iteration starts. 1. When I define startTime 0 with DeltaT 0.00001, it says floating point exception (core dumped) Time = 1e-05, iteration: 1 Current fsi under-relaxation factor: 0.4 Creating accumulated fluid interface displacement Maximal accumulated displacement of interface points: 0 GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1 GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1 GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1 GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1 Evolving flow model: consistentIcoFlow Courant Number mean: 2.14175e-07 max: 6e-06 velocity magnitude: 0.0001 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.31844e-07, No Iterations 5 DILUPBiCG: Solving for Uy, Initial residual = 0.999999, Final residual = 3.66704e-07, No Iterations 4 GAMG: Solving for p, Initial residual = 1, Final residual = 5.30983e-07, No Iterations 7 GAMG: Solving for p, Initial residual = 0.0114117, Final residual = 9.28912e-07, No Iterations 67 time step continuity errors : sum local = 3.64821e-14, global = -3.21772e-24 GAMG: Solving for p, Initial residual = 0.266161, Final residual = 8.08369e-07, No Iterations 67 GAMG: Solving for p, Initial residual = 8.50573e-07, Final residual = 5.47487e-07, No Iterations 2 time step continuity errors : sum local = 2.46898e-14, global = 9.27448e-27 GAMG: Solving for p, Initial residual = 0.288257, Final residual = 8.85958e-07, No Iterations 51 GAMG: Solving for p, Initial residual = 9.1568e-07, Final residual = 6.94922e-07, No Iterations 2 time step continuity errors : sum local = 3.34729e-14, global = -3.38972e-24 Setting traction on solid patch Total force (fluid) = (0.00870746 2.02047e-05 1.30104e-13) Total force (solid) = (-0.00870746 -2.02047e-05 -1.30104e-13) Evolving stress model -1.65018, 1.0115 DICPCG: Solving for D, Initial residula = 1, Final residual = 0.000587252, No outer iterations = 108 Max relative residual = 24.1206, Relative residual = 0.000514779, enforceLinear = 1 Current fsi relative residual norm: 1 Time = 1e-05, iteration: 2 Current fsi under-relaxation factor: 0.4 Maximal accumulated displacement of interface points: 2.95406e+19 GAMG: Solving for cellMotionUx, Initial residual = 1, Final residual = 0.000189476, No Iterations 2 GAMG: Solving for cellMotionUy, Initial residual = 1, Final residual = 0.000189394, No Iterations 2 GAMG: Solving for cellMotionUx, Initial residual = 9.24019e-05, Final residual = 4.92158e-07, No Iterations 4 GAMG: Solving for cellMotionUy, Initial residual = 9.23671e-05, Final residual = 4.91437e-07, No Iterations 4 Evolving flow model: consistentIcoFlow Floating point exception (core dumped) 2. When I define startTime 0 with DeltaT 0.0000001, it says continuity errors cannot be removed by adjusting the outflow. Time = 1e-07, iteration: 1 Current fsi under-relaxation factor: 0.4 Creating accumulated fluid interface displacement Maximal accumulated displacement of interface points: 0 GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1 GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1 GAMG: Solving for cellMotionUx, Initial residual = 0, Final residual = 0, No Iterations 1 GAMG: Solving for cellMotionUy, Initial residual = 0, Final residual = 0, No Iterations 1 Evolving flow model: consistentIcoFlow Courant Number mean: 2.14175e-09 max: 6e-08 velocity magnitude: 0.0001 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 8.64249e-10, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.999999, Final residual = 7.81295e-07, No Iterations 1 GAMG: Solving for p, Initial residual = 1, Final residual = 9.35145e-07, No Iterations 20 GAMG: Solving for p, Initial residual = 0.000371717, Final residual = 9.18606e-07, No Iterations 20 time step continuity errors : sum local = 3.90917e-16, global = -3.72126e-26 GAMG: Solving for p, Initial residual = 0.184356, Final residual = 8.65415e-07, No Iterations 32 GAMG: Solving for p, Initial residual = 6.96729e-07, Final residual = 3.55687e-07, No Iterations 2 time step continuity errors : sum local = 2.19023e-16, global = 1.67978e-30 --> FOAM FATAL ERROR: Continuity error cannot be removed by adjusting the outflow. Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow. Specified mass inflow : 1.29924e-07 Specified mass outflow : 9.35926e-20 Difference : 1.29924e-07 Adjustable mass outflow : 4.03464e-16 From function adjustPhi ( surfaceScalarField& phi, const volVectorField& U, const volScalarField& p ) in file cfdTools/general/adjustPhi/adjustPhi.C at line 150. FOAM exiting Could you please explain whats happening here!! Thanks!!

November 16, 2016, 23:45		#95
Woj3x New Member Wojciech Gołąbek Join Date: Dec 2013 Posts: 29 Rep Power: 12	Dear KarthickRajkumar First of all it would be nice if you could tell me if is it still a BeamInCrossFlow case? I see that it is 2D case so I am not sure about this. How much did you change this case? Without this informations it's really hard to tell whats wrong with it. IMO it can be boundary conditions, material parameters or even script problem (ex. velocityInlet from Allrun). Wojciech

November 24, 2016, 05:39		#99
Woj3x New Member Wojciech Gołąbek Join Date: Dec 2013 Posts: 29 Rep Power: 12	Hi Karthick, Sorry I did not answer but I had a lots to do in last week. Could you look into your log.fsifoam file and tell us how big is Courant number in your results?