CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[Other] simulation of closing the gate using moving mesh

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

Like Tree2Likes
  • 1 Post By fabian_roesler
  • 1 Post By alfa_8C

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 7, 2014, 16:32
Default simulation of closing the gate using moving mesh
  #1
New Member
 
simin
Join Date: Apr 2012
Posts: 13
Rep Power: 14
simin_ds is on a distinguished road
Hi Foamers,

I am going to simulate closing a gate in a channel using moving mesh. I do not know which method I should use to close the gate. Does any one has any suggestion? is there any tutorial for closing a valve or gate in 2d or 3d?

Thank you in advance,
Simin
simin_ds is offline   Reply With Quote

Old   November 10, 2014, 01:28
Default closing gate using oscillatingInletACMI2D
  #2
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 18
fabian_roesler is on a distinguished road
Hi,

it depends on whether you have a rotating or a translatory motion or even a mixture of both in several directions. I'd start with a simple patch moving aside. You could use a sliding interface like in incompressible/pimpleDyMFoam/oscillatingInletACMI2D tutorial. However I rarely used dynamic meshes.

Cheers

Fabian
saatt likes this.
fabian_roesler is offline   Reply With Quote

Old   November 10, 2014, 01:37
Default
  #3
New Member
 
simin
Join Date: Apr 2012
Posts: 13
Rep Power: 14
simin_ds is on a distinguished road
Hi Fabian,
Actually i used this tutorial. But for the gate closure it has no correct mesh motion by closing the gate and it represent the bad mesh skewness near the moving interface. The gate in my case is moving from a side and it is not a moving interface.
Regards,
Simin
simin_ds is offline   Reply With Quote

Old   November 13, 2014, 11:26
Default inapplicability of moving mesh in gate closure
  #4
New Member
 
simin
Join Date: Apr 2012
Posts: 13
Rep Power: 14
simin_ds is on a distinguished road
Hi Foamers,

I sent a post a few days ago for my problem in using moving mesh in OF but I could not solve my problem. Today I explain my problem in detail. Hope the experts in OF may help me.

I am going to simulate the gate closure in a channel. for my case I assume a 2D channel. you can see the geometry in Fig. 1.

I used the structure of movingcone and cyclic-ACMI for my moving mesh simulation. the dynamic mesh dict file is shown here.
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solver velocityComponentLaplacian;

velocityComponentLaplacianCoeffs
{
component y;
diffusivity directional ( 1 200 0 );
}


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


also the pressure file is shown here
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
walls
{
type zeroGradient;
}
valve
{
type zeroGradient;
}
couple1
{
type zeroGradient;
}
couple2
{
type zeroGradient;
}
couple3
{
type zeroGradient;
}
couple4
{
type zeroGradient;
}
defaultFaces
{
type empty;
}
ACMI1_blockage
{
type zeroGradient;
}
ACMI1_couple
{
type cyclicACMI;
value uniform 0;
}
ACMI2_blockage
{
type zeroGradient;
}
ACMI2_couple
{
type cyclicACMI;
value uniform 0;
}
ACMI3_blockage
{
type zeroGradient;
}
ACMI3_couple
{
type cyclicACMI;
value uniform 0;
}
ACMI4_blockage
{
type zeroGradient;
}
ACMI4_couple
{
type cyclicACMI;
value uniform 0;
}
}


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

and velocity file:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (1 0 0);
}
outlet
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
walls
{
// type fixedValue;
type movingWallVelocity;
value uniform (0 0 0);
}
valve
{
// type fixedValue;
type movingWallVelocity;
value uniform (0 0 0);
}
defaultFaces
{
type empty;
}
ACMI1_blockage
{
type fixedValue;
value uniform (0 0 0);
}
ACMI1_couple
{
type cyclicACMI;
value uniform (0 0 0);
}
ACMI2_blockage
{
type fixedValue;
value uniform (0 0 0);
}
ACMI2_couple
{
type cyclicACMI;
value uniform (0 0 0);
}
ACMI3_blockage
{
type fixedValue;
value uniform (0 0 0);
}
ACMI3_couple
{
type cyclicACMI;
value uniform (0 0 0);
}
ACMI4_blockage
{
type fixedValue;
value uniform (0 0 0);
}
ACMI4_couple
{
type cyclicACMI;
value uniform (0 0 0);
}
}


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

and also the pointMotionUy file

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class pointScalarField;
object pointMotionU;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{

inlet
{
type uniformFixedValue;
uniformValue constant 0;
}
outlet
{
type uniformFixedValue;
uniformValue constant 0;
}
walls
{

type uniformFixedValue;
uniformValue constant 0;
}
valve
{
type uniformFixedValue;
uniformValue constant -0.01;
}
defaultFaces
{
type empty;
}
couple1
{
type slip;
}
couple2
{
type slip;
}
couple3
{
type slip;
}
couple4
{
type slip;
}
ACMI1_blockage
{
type fixedValue;
value uniform 0;
}
ACMI1_couple
{
type slip;
//type cyclicACMI;
// value uniform 0;
}
ACMI2_blockage
{
type slip;
//type fixedValue;
//value uniform 0;
}
ACMI2_couple
{
type slip;
//type cyclicACMI;
// value uniform 0;
}
ACMI3_blockage
{
type slip;
//type fixedValue;
//value uniform 0;
}
ACMI3_couple
{
type slip;
//type cyclicACMI;
// value uniform 0;
}
ACMI4_blockage
{
type fixedValue;
value uniform 0;
}
ACMI4_couple
{
type slip;
//type cyclicACMI;
// value uniform 0;
}
}

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

I put also the toposet dict file here:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object topoSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

actions
(
// Get both sides of ami
// ~~~~~~~~~~~~~~~~~~~~~

// Create faceZone for patch couple1
{
name couple1Faces;
type faceSet;
action new;
source patchToFace;
sourceInfo
{
name couple1;
}
}
{
name couple1Faces;
type faceZoneSet;
action new;
source setToFaceZone;
sourceInfo
{
faceSet couple1Faces;
}
}

// Create faceZone for patch couple2
{
name couple2Faces;
type faceSet;
action new;
source patchToFace;
sourceInfo
{
name couple2;
}
}
{
name couple2Faces;
type faceZoneSet;
action new;
source setToFaceZone;
sourceInfo
{
faceSet couple2Faces;
}
}

{
name couple3Faces;
type faceSet;
action new;
source patchToFace;
sourceInfo
{
name couple3;
}
}
{
name couple3Faces;
type faceZoneSet;
action new;
source setToFaceZone;
sourceInfo
{
faceSet couple3Faces;
}
}

// Create faceZone for patch couple2
{
name couple4Faces;
type faceSet;
action new;
source patchToFace;
sourceInfo
{
name couple4;
}
}
{
name couple4Faces;
type faceZoneSet;
action new;
source setToFaceZone;
sourceInfo
{
faceSet couple4Faces;
}
}




// Create cellZone for moving cells in inlet channel
{
name inletChannel;
type cellSet;
action new;
source boxToCell;
sourceInfo
{
box (2.00001 -100 -100) (2.9999 100 100);
}
}
{
name inletChannel;
type cellZoneSet;
action new;
source setToCellZone;
sourceInfo
{
set inletChannel;
}
}
);

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


My problem is that when the gate starts to travel down the nodes near the interfaces do not move with the valve movement. so the negative cell volume occurs in the simulation. you can see the the last step of the simulation in Fig 2.

Can any one help me what boundary condition or dynamic mesh method I should use to solve this problem?

Many thanks for your attention.
Attached Images
File Type: jpg Screenshot from 2014-11-13 18:06:35.jpg (52.4 KB, 127 views)
File Type: jpg Screenshot from 2014-11-13 18:06:58.jpg (64.8 KB, 127 views)
simin_ds is offline   Reply With Quote

Old   November 13, 2014, 14:44
Default
  #5
Senior Member
 
mturcios777's Avatar
 
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 740
Rep Power: 28
mturcios777 will become famous soon enough
If you look at the oscillatingInletACMI2D case in the pimpleDyMFoam tutorials, you might be able to use that approach. I tested changing the amplitude of the motion so that the inlet would become disconnected from the main volume, and for incompressible flow it worked ok.
mturcios777 is offline   Reply With Quote

Old   March 16, 2016, 08:32
Default
  #6
Member
 
Antonio Liggieri
Join Date: Aug 2010
Posts: 76
Rep Power: 15
alfa_8C is an unknown quantity at this point
Hello simin

it's quite a while ago since you posted your issue.

Could you resolve this?

Best regards
Antonio
alfa_8C is offline   Reply With Quote

Old   March 21, 2017, 20:31
Default
  #7
New Member
 
Join Date: Mar 2017
Posts: 28
Rep Power: 9
saatt is on a distinguished road
hey ! i have a same question , could you done well ??
saatt is offline   Reply With Quote

Old   March 22, 2017, 02:34
Default not solved yet
  #8
Member
 
Antonio Liggieri
Join Date: Aug 2010
Posts: 76
Rep Power: 15
alfa_8C is an unknown quantity at this point
Hi,

you may follow the progress on this issue on the following link:

http://www.openfoam.org/mantisbt/view.php?id=2057

alfa_8C
saatt likes this.
alfa_8C is offline   Reply With Quote

Old   April 12, 2019, 05:49
Default
  #9
Member
 
JuanMi
Join Date: Nov 2017
Posts: 41
Rep Power: 8
keitaro7_14 is on a distinguished road
Luckly, it seems the problem has been solved, but only for the OpenFOAM-dev version.
keitaro7_14 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 37 July 20, 2020 05:38
Fatal overflow error, singularity, due to local mesh form? Thomasbr CFX 7 July 8, 2016 09:57
moving mesh simulation: cavity with a moving lid cover badboyz31 CFX 7 December 22, 2015 17:29
[mesh manipulation] Importing Multiple Meshes thomasnwalshiii OpenFOAM Meshing & Mesh Conversion 18 December 19, 2015 18:57
Question on moving mesh, mesh velocity is really small! ripperjack Main CFD Forum 2 April 28, 2014 13:37


All times are GMT -4. The time now is 20:47.