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

SlidingInterface problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 5, 2009, 09:48
Default Hello ! I have installed the
  #21
New Member
 
Naish
Join Date: Mar 2009
Location: Germany
Posts: 11
Rep Power: 8
naish is on a distinguished road
Hello !
I have installed the 1.5-dev version on my laptop and the IcoDyMesh works but I Still have the same problem on my 64bit machine where can I get a 1.5-dev version for 64bit ?

thanks N.
naish is offline   Reply With Quote

Old   January 5, 2009, 16:57
Default Hi guys I really need someo
  #22
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
Hi guys

I really need someone that explains me how to create a sliding interface case and the solver/BC settings.

Thanks in advance
antonio_ing is offline   Reply With Quote

Old   January 9, 2009, 17:54
Default Hi Antonio, which version o
  #23
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hi Antonio,

which version of OF do you intend to run a sliding interface on? It would be preferable for you to use 1.5-dev.

regards,
-Louis
louisgag is offline   Reply With Quote

Old   January 9, 2009, 17:56
Default Hi Naish, I think you can s
  #24
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Hi Naish,

I think you can simply compile the 1.5-dev version on your 64bit machine and it will give 64bit binaries.. That's what I did.

cheers,
-Louis
louisgag is offline   Reply With Quote

Old   January 10, 2009, 07:02
Default actually i'm working with open
  #25
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
actually i'm working with openFoam1.5. Should I switch to openfoam1.5-dev? Where can I find all the files that i need. I'm sorry.... i'm new of OpenFoam.

In the meantime, i'm trying to work with the tutorial of icoDyMFoam that doesn't use the sliding interface method but the Automatic motion of the Mesh (movingcone tutorial), but in my case i need a BC that is not constant, but depending on the position Vel_patch=omega*(-y,x,0). I took the parabolic inlet code and slightly modified it in this form:

void parabolicVelocityFvPatchVectorField::updateCoeffs( )
{
// Get range and orientation
boundBox bb(patch().patch().localPoints(), false);

vector ctr = 0.5*(bb.max() + bb.min());

const vectorField& c = patch().Cf();

// Calculate local 1-D coordinate for the parabolic profile
scalarField coord = 2*((c - ctr) & y_)/((bb.max() - bb.min()) & y_);
scalarField cc1 = c & n_;
scalarField cc2 = c & y_;

/*vectorField::operator=(n_*maxValue_*(1.0 - sqr(coord)));*/
vectorField::operator=-maxValue_*cc2*n_+maxValue_*cc1*y_;
}

where n_=(1,0,0) and y_=(1,0,0). But when i compile with wmake libso it gives me an error.
Did i Do something stupid? i'm not really good in c++... sorry :-(. Has anybody something ready for me? thanks very much
antonio_ing is offline   Reply With Quote

Old   January 10, 2009, 13:45
Default http://www.cfd-online.com/cgi-
  #26
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
http://www.cfd-online.com/cgi-bin/Op...0168#POST20168

Take the version from the powerlab website, it's the one that works with GGI.


Cheers,
-Louis
louisgag is offline   Reply With Quote

Old   January 13, 2009, 12:14
Default Finally, i have been able to c
  #27
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
Finally, i have been able to compile the dev version and run the mixed2D case. Once i turned to my turbine again i have defined two patches in gmsh (insideSlider and outsideSlider) that are coincident and with 0 faces (i have seen this in the boundary) but these are circular cylinders in my 2d domain. When i compiled i have obtained the following error:


/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5-dev |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : icoDyMFoam
Date : Jan 13 2009
Time : 17:08:20
Host : antonio-laptop
PID : 8754
Case : /home/cfduser-dev/OpenFOAM/cfduser-dev-1.5-dev/darrieus_icoDyMFoam
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh mixerFvMesh


Not all zones and patches needed in the definition have been found. Please check your mesh definition.

From function void slidingInterface::checkDefinition()
in file polyMeshModifiers/slidingInterface/slidingInterface.C at line 96.

FOAM aborting

Aborted


did I do something wrong?
antonio_ing is offline   Reply With Quote

Old   January 13, 2009, 14:21
Default Antonio, In gmsh I make two
  #28
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Antonio,

In gmsh I make two different meshes and merge them with mergeMeshes utility and it works fine.

cheers,
-Louis
louisgag is offline   Reply With Quote

Old   January 13, 2009, 14:27
Default you are great Louis. I'll try
  #29
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
you are great Louis. I'll try it immediatly and let you know soon :-)
antonio_ing is offline   Reply With Quote

Old   January 14, 2009, 13:14
Default i tried that but i got this er
  #30
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
i tried that but i got this error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5-dev |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : icoDyMFoam
Date : Jan 14 2009
Time : 18:12:31
Host : antonio-laptop
PID : 12195
Case : /home/cfduser-dev/OpenFOAM/cfduser-dev-1.5-dev/test_ext
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh mixerFvMesh
Rotating region marker point: (0 0 0)
Attach-detach action = false
void mixerFvMesh::addZonesAndModifiers() : Zones and modifiers already present. Skipping.


Mesh modifiers not read properly

From function void mixerFvMesh::addZonesAndModifiers()
in file mixerFvMesh/mixerFvMesh.C at line 66.

FOAM aborting

Aborted
antonio_ing is offline   Reply With Quote

Old   January 14, 2009, 13:28
Default Go to constant/polyMesh and de
  #31
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Go to constant/polyMesh and delete all Zones files and meshModifiers. Then, compare your boundary file with the definition of the mixer, especially the sliding patches (in constant/dynamicMeshDict). If all is well, try again.

It looks like you are referring to patches in constant/dynamicMeshDict) that do not exist in the boundary.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   January 14, 2009, 14:52
Default I have deleted the zones files
  #32
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
I have deleted the zones files and the meshmodifiers but i got this error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.5-dev |
| \ / A nd | Web: http://www.OpenFOAM.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : icoDyMFoam
Date : Jan 14 2009
Time : 19:50:44
Host : antonio-laptop
PID : 20238
Case : /home/cfduser-dev/OpenFOAM/cfduser-dev-1.5-dev/test_ext
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh mixerFvMesh
Rotating region marker point: (0 0 0)
Attach-detach action = false
void mixerFvMesh::addZonesAndModifiers() : Zones and modifiers already present. Skipping.
Mixer mesh
origin: (0 0 0)
axis : (0 0 1)
rpm : 10
Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Reading field rAU if present



This mesh contains patches of type empty but is not 1D or 2D
by virtue of the fact that the number of faces of this
empty patch is not divisible by the number of cells.

From function emptyFvPatchField<type>::updateCoeffs()
in file fields/fvPatchFields/constraint/empty/emptyFvPatchField.C at line 148.

FOAM exiting
antonio_ing is offline   Reply With Quote

Old   January 14, 2009, 15:10
Default Antonio, setting your empty
  #33
Senior Member
 
louisgag's Avatar
 
Louis Gagnon
Join Date: Mar 2009
Location: Québec, QC, Canada
Posts: 208
Rep Power: 9
louisgag is on a distinguished road
Send a message via ICQ to louisgag
Antonio,

setting your empty patches to the same name in both gmsh meshes might help.

cheers,
-Louis
louisgag is offline   Reply With Quote

Old   January 14, 2009, 15:51
Default thanks Louis and Hrvoje. my ca
  #34
Member
 
antonio segalini
Join Date: Mar 2009
Posts: 75
Rep Power: 8
antonio_ing is on a distinguished road
thanks Louis and Hrvoje. my case is running now!!!!
antonio_ing is offline   Reply With Quote

Old   August 10, 2009, 06:45
Default
  #35
Senior Member
 
Philippose Rajan
Join Date: Mar 2009
Location: Germany
Posts: 530
Rep Power: 16
philippose will become famous soon enough
Hello there,

A Good Afternoon to everyone!

I have three questions regarding the current implementation of the sliding interface in OpenFOAM 1.5-dev:

1. Do I need to define a separate sliding interface mesh modifier for each pair of physical sliding patches, or can I group all the faces which I consider to be part of the "master" patch and all the faces which I consider to be part of the "slave" patch into just two patch names "Master" and "Slave" even though the faces do not form physically continuous / connected surfaces (for example, the slave patch represents two physically disconnected surfaces but both of which need to slide along a master patch)?

2. Do I need to have a very small distance between the master and the slave patches to ensure that no parts of the slave patch at any point intersect the master patch? (Even if the intersections are caused by the meshing process?)

3. Are there any special requirements which need to be met with regard to the surface meshes on the master and slave patches..... such as.... approximately similar element sizes, etc... and is there any connection between the mesh element sizes, and the velocity with which the sliding interface can be moved?

Have a great day ahead !!

Regards,
Philippose
philippose is offline   Reply With Quote

Old   September 15, 2009, 02:15
Default
  #36
New Member
 
shyam prasad
Join Date: Mar 2009
Posts: 25
Rep Power: 8
shyam is on a distinguished road
Hi Foamers,
I tried to run the above mixer2D tutorial from Jassi using pimpleDyMFoam in OF1.6, I get the following error. Any suggestions/advise will be helpful.

Unknown dynamicFvMesh type mixerFvMesh
Valid dynamicFvMesh types are :
5
(
dynamicMotionSolverFvMesh
solidBodyMotionFvMesh
dynamicRefineFvMesh
dynamicInkJetFvMesh
staticFvMesh
)

From function dynamicFvMesh::New(const IOobject&)
in file dynamicFvMesh/newDynamicFvMesh.C at line 79.
shyam is offline   Reply With Quote

Old   December 31, 2010, 16:37
Question turbDyMFoam for rotation
  #37
New Member
 
Join Date: Dec 2010
Posts: 26
Rep Power: 6
prashant.A is on a distinguished road
Send a message via Skype™ to prashant.A
Hello,

this is my first case of actual mesh rotation using "turbDyMFoam" and I am facing stability issues. The calculation goes fine for first 50 deg of rotation but eventually solution blows up!
I tried limiting Courant No. <=1 but still doesn't help..
Though going below this value may render some stability, but that looks like long long calculation times.

I have switched to PCG from PBiCG for the moment and it is running too slow..

Any suggestions !
prashant.A is offline   Reply With Quote

Old   January 2, 2011, 13:12
Default
  #38
New Member
 
Join Date: Dec 2010
Posts: 26
Rep Power: 6
prashant.A is on a distinguished road
Send a message via Skype™ to prashant.A
I think I found a workaround to below, run it for smaller timesteps, meaning bring down the rotation speed for first few degrees and once it gets going, it can take bigger time steps.
I even made the rpm = 0 in one of the trials, but that may not be required though !
I am also trying transientSimpleDyMFoam to see how that goes.

Thanks,
Prashant

Quote:
Originally Posted by prashant.A View Post
Hello,

this is my first case of actual mesh rotation using "turbDyMFoam" and I am facing stability issues. The calculation goes fine for first 50 deg of rotation but eventually solution blows up!
I tried limiting Courant No. <=1 but still doesn't help..
Though going below this value may render some stability, but that looks like long long calculation times.

I have switched to PCG from PBiCG for the moment and it is running too slow..

Any suggestions !
prashant.A is offline   Reply With Quote

Old   November 7, 2011, 15:01
Default
  #39
Member
 
s.rasoul_varedi
Join Date: Feb 2010
Posts: 82
Rep Power: 6
desert_1250 is an unknown quantity at this point
Send a message via Yahoo to desert_1250
dear foamers
i have one question about sliding interface and moving mesh. i study about VAWT with 3 blade using GGI connectors and i simulated it with consideration of fixed blade. now i wanna simulate it with considering variable pitch for every blade so that extract maximum power from it.
1-how to define multi sliding interface so that each blade rotate from its' axis by a definition function, as rotate simultaneously around shaft

2- how to definition a function for every blades so that at every time step change its the angle of attack

tanx
______
Rasoul
desert_1250 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
VOF diffuse problem on two fluids problem Fang Jin FLUENT 5 February 17, 2011 06:38
Problem in Modelling Heat Transfer Problem Deepak R FLUENT 1 December 6, 2007 10:37
SlidingInterface Tolerance graser OpenFOAM Running, Solving & CFD 2 November 27, 2007 07:53
Problem in Tutorial problem of fluent Phanindra FLUENT 5 April 17, 2007 09:57
problem in solving "wave generation" problem san FLUENT 2 April 3, 2006 23:37


All times are GMT -4. The time now is 14:50.