CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   SlidingInterface problem (https://www.cfd-online.com/Forums/openfoam-solving/58107-slidinginterface-problem.html)

naish January 5, 2009 08:48

Hello ! I have installed the
 
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.

antonio_ing January 5, 2009 15:57

Hi guys I really need someo
 
Hi guys

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

Thanks in advance

louisgag January 9, 2009 16:54

Hi Antonio, which version o
 
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 January 9, 2009 16:56

Hi Naish, I think you can s
 
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

antonio_ing January 10, 2009 06:02

actually i'm working with open
 
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

louisgag January 10, 2009 12:45

http://www.cfd-online.com/cgi-
 
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

antonio_ing January 13, 2009 11:14

Finally, i have been able to c
 
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?

louisgag January 13, 2009 13:21

Antonio, In gmsh I make two
 
Antonio,

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

cheers,
-Louis

antonio_ing January 13, 2009 13:27

you are great Louis. I'll try
 
you are great Louis. I'll try it immediatly and let you know soon :-)

antonio_ing January 14, 2009 12:14

i tried that but i got this er
 
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

hjasak January 14, 2009 12:28

Go to constant/polyMesh and de
 
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

antonio_ing January 14, 2009 13:52

I have deleted the zones files
 
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

louisgag January 14, 2009 14:10

Antonio, setting your empty
 
Antonio,

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

cheers,
-Louis

antonio_ing January 14, 2009 14:51

thanks Louis and Hrvoje. my ca
 
thanks Louis and Hrvoje. my case is running now!!!!

philippose August 10, 2009 06:45

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

shyam September 15, 2009 02:15

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.

prashant.A December 31, 2010 15:37

turbDyMFoam for rotation
 
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:confused: !

prashant.A January 2, 2011 12:12

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 (Post 288947)
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:confused: !


desert_1250 November 7, 2011 14:01

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


All times are GMT -4. The time now is 15:57.