CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

HronTurekFsi3 foam-extend3.1

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

Like Tree15Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   April 17, 2016, 20:42
Default
  #81
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
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:
  1. 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.
  2. 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.
  3. 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.
  4. To run the case, simply do these steps:
    Code:
    cd fluid
    ./Allrun
  5. To restart the case, run from within the folder "fluid" this command:
    Code:
    ./Allclean
  6. 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!
Attached Files
File Type: gz beamInCross_sculptored_v1.tar.gz (13.5 KB, 13 views)
Woj3x likes this.
wyldckat is offline   Reply With Quote

Old   April 18, 2016, 18:23
Default
  #82
New Member
 
Wojciech Gołąbek
Join Date: Dec 2013
Posts: 29
Rep Power: 6
Woj3x is on a distinguished road
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

Last edited by Woj3x; April 18, 2016 at 19:50.
Woj3x is offline   Reply With Quote

Old   April 19, 2016, 18:45
Default
  #83
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
@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.
wyldckat is offline   Reply With Quote

Old   April 25, 2016, 07:19
Default
  #84
Senior Member
 
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 7
Maimouna is on a distinguished road
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 is offline   Reply With Quote

Old   April 26, 2016, 12:15
Default
  #85
Senior Member
 
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 7
Maimouna is on a distinguished road
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 is offline   Reply With Quote

Old   May 1, 2016, 05:34
Default
  #86
Senior Member
 
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 7
Maimouna is on a distinguished road
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
Maimouna is offline   Reply With Quote

Old   May 1, 2016, 19:53
Default
  #87
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quick answers:
Quote:
Originally Posted by Maimouna View Post
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 View Post
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?
  1. 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.
  2. "blockMeshDict.first" isn't needed. It's something from some other case, I don't remember which one.
  3. 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 and Woj3x like this.
wyldckat is offline   Reply With Quote

Old   May 2, 2016, 17:41
Default
  #88
Senior Member
 
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 7
Maimouna is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
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...


  1. 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.
  2. "blockMeshDict.first" isn't needed. It's something from some other case, I don't remember which one.
  3. 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
Maimouna is offline   Reply With Quote

Old   May 3, 2016, 16:31
Default
  #89
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,125
Blog Entries: 39
Rep Power: 110
wyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of lightwyldckat is a glorious beacon of light
Quote:
Originally Posted by Maimouna View Post
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.
wyldckat is offline   Reply With Quote

Old   May 8, 2016, 11:11
Default
  #90
New Member
 
Wojciech Gołąbek
Join Date: Dec 2013
Posts: 29
Rep Power: 6
Woj3x is on a distinguished road
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.

Last edited by Woj3x; May 8, 2016 at 19:35.
Woj3x is offline   Reply With Quote

Old   May 13, 2016, 07:16
Default
  #91
Senior Member
 
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 7
Maimouna is on a distinguished road
Quote:
Originally Posted by Woj3x View Post
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
Maimouna is offline   Reply With Quote

Old   November 13, 2016, 17:08
Default
  #92
New Member
 
Karthick
Join Date: Oct 2016
Posts: 15
Rep Power: 3
KarthickRajkumar is on a distinguished road
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)
KarthickRajkumar is offline   Reply With Quote

Old   November 13, 2016, 19:17
Default
  #93
New Member
 
Wojciech Gołąbek
Join Date: Dec 2013
Posts: 29
Rep Power: 6
Woj3x is on a distinguished road
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 likes this.
Woj3x is offline   Reply With Quote

Old   November 16, 2016, 17:02
Default
  #94
New Member
 
Karthick
Join Date: Oct 2016
Posts: 15
Rep Power: 3
KarthickRajkumar is on a distinguished road
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!!
KarthickRajkumar is offline   Reply With Quote

Old   November 17, 2016, 00:45
Default
  #95
New Member
 
Wojciech Gołąbek
Join Date: Dec 2013
Posts: 29
Rep Power: 6
Woj3x is on a distinguished road
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
Woj3x is offline   Reply With Quote

Old   November 17, 2016, 04:02
Default
  #96
New Member
 
Karthick
Join Date: Oct 2016
Posts: 15
Rep Power: 3
KarthickRajkumar is on a distinguished road
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!!
Attached Files
File Type: zip 1.FSI Project(I).zip (26.7 KB, 2 views)
KarthickRajkumar is offline   Reply With Quote

Old   November 24, 2016, 05:27
Default
  #97
New Member
 
Karthick
Join Date: Oct 2016
Posts: 15
Rep Power: 3
KarthickRajkumar is on a distinguished road
Quote:
Originally Posted by Maimouna View Post
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
KarthickRajkumar is offline   Reply With Quote

Old   November 24, 2016, 05:52
Default
  #98
Senior Member
 
ok
Join Date: Oct 2013
Posts: 346
Rep Power: 7
Maimouna is on a distinguished road
Quote:
Originally Posted by KarthickRajkumar View Post
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
Maimouna is offline   Reply With Quote

Old   November 24, 2016, 06:39
Default
  #99
New Member
 
Wojciech Gołąbek
Join Date: Dec 2013
Posts: 29
Rep Power: 6
Woj3x is on a distinguished road
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?
Woj3x is offline   Reply With Quote

Old   November 24, 2016, 08:17
Default
  #100
New Member
 
Karthick
Join Date: Oct 2016
Posts: 15
Rep Power: 3
KarthickRajkumar is on a distinguished road
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
Attached Files
File Type: gz HronTurekFsi3 (check-5).tar.gz (12.5 KB, 7 views)
KarthickRajkumar is offline   Reply With Quote

Reply

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
error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Native Meshers: blockMesh 2 March 14, 2012 10:56
mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 05:16
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00
Import gmsh msh to Foam adorean Open Source Meshers: Gmsh, Netgen, CGNS, ... 24 April 27, 2005 08:19


All times are GMT -4. The time now is 09:49.