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

engineTopochangerMesh simpleTwoStroke layerAR does not work

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 20, 2017, 09:10
Default engineTopochangerMesh simpleTwoStroke layerAR does not work
  #1
New Member
 
Rainer Kotmann
Join Date: May 2012
Location: Aachen, Germany
Posts: 1
Rep Power: 0
KotRainer11 is on a distinguished road
Dear Foamers,

I want to simulate a simple compressor where a piston is compressing a fluid and transporting it into a channel. Since this is nearly the same as in a two-stroke engine, I used the engineTopoChangerMesh simpleTwoStroke. At first, I want to simulate a simple 2D case. Please see the initial Setup at top bottom center in the figure -180.png attached.
My engineGeometry file looks as follows:

Code:
 
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.0                                |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      engineGeometry;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//engineMesh     layered;
engineTopoChangerMesh      simpleTwoStroke;
//engineTopoChangerMesh      layerAR;

conRodLength    conRodLength [ 0 1 0 0 0 0 0 ] 0.35;
bore            bore [ 0 1 0 0 0 0 0 ] 0.5;
stroke          stroke [ 0 1 0 0 0 0 0 ] 0.15;
clearance       clearance [ 0 1 0 0 0 0 0 ] 0.002;
rpm             rpm [ 0 0 -1 0 0 0 0 ] 735;
cylinderHead cylinderHead;
piston
{
 patch piston;
 coordinateSystem
 {
  origin (0 0 0);
  coordinateRotation
  {
   type axes;
   e1 (1 0 0);
   e2 (0 1 0);
   e3 (0 0 1);
  }
  
 }
 minLayer 0.0005;
 maxLayer 0.002;
}
deformAngle 38;
delta  0.001;
offSet  0.002;
checkMesh on;
scavInCylPatch liner;
scavInPortPatch scavenging;
scavPortsTol 1e-3;
coordinateSystem
    {
        type            cartesian;
        origin          (0 0 0);
        axis            (0 0 1);
        direction       (0 1 0);
    }
}
// ************************************************************************* //
When I start moving the mesh with moveDyMEngineMesh, I get the follwing Output for the first two time steps:

Code:
 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.0                                |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build    : 4.0-64ca3820ca97
Exec     : moveDyMEngineMesh
Date     : Apr 20 2017
Time     : 10:14:19
Host     : up720
PID      : 2610
CtrlDict : "/home/nordent/foam/nordent-4.0/run/ameria_2d_valve/system/controlDict"
Case     : /home/nordent/foam/nordent-4.0/run/ameria_2d_valve
nProcs   : 1
SigFpe   : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create engine time
Create mesh
Selecting engineTopoChangerMesh simpleTwoStroke
checkAndCalculate()
piston
cylinderHead
liner
valve
scavenging
outlet
front
back
deckHeight: 0.152
piston position: 0
Time = -180
Adding zones to the engine mesh
Adding face zone for piston layer addition/removal
Number of head points = 502
Number of cells in the moving region: 38000
Adding 2 point, 4 face zones and 1 cell zones
Adding 2 topology modifiers
piston face index = 0
virtualPistonPosition = 0.002
piston position = 0
Time = -179.5
Sliding interfaces decoupled
Executing layer action
virtualPistonPosition = 0.00200224, deckHeight = 0.152, pistonPosition = 0
Piston layering mode
Coupling sliding interfaces
Sliding interfaces coupled: 1
volume continuity errors : volume = 0.00162, old volume = 0.00162, max error = 2.77933e-14, sum local = 2.50461e-19, global = 6.25311e-21
Mesh Courant Number mean: 0.00137621 max: 0.00448766
    Point usage OK.
    Upper triangular ordering OK.
    Topological cell zip-up check OK.
    Face vertices OK.
    Face-face connectivity OK.
    Mesh topology OK.
    Boundary openness (-4.24201e-18 1.64967e-15 -3.01046e-18) Threshold = 1e-06 OK.
    Max cell openness = 1.05879e-16 OK.
    Max aspect ratio = 20.0449 OK.
    Minumum face area = 1e-06. Maximum face area = 4e-05.  Face area magnitudes OK.
    Min volume = 2e-08. Max volume = 4e-08.  Total volume = 0.00162.  Cell volumes OK.
    Mesh non-orthogonality Max: 0.128562 average: 0.00225346 Threshold = 70
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.00224382 OK.
    Mesh geometry OK.
Mesh OK.
ExecutionTime = 1.24 s  ClockTime = 1 s
Time = -179
Decoupling sliding interfaces
sliding interfaces successfully decoupled!!!
Executing layer action
virtualPistonPosition = 0.00200898, deckHeight = 0.152, pistonPosition = 2.24383e-06
Piston layering mode
Coupling sliding interfaces
Sliding interfaces coupled: 1
volume continuity errors : volume = 0.00162, old volume = 0.00162, max error = 1.86943e-14, sum local = 3.05258e-19, global = 2.81601e-23
Mesh Courant Number mean: 0.00412862 max: 0.0134629
    Point usage OK.
    Upper triangular ordering OK.
    Topological cell zip-up check OK.
    Face vertices OK.
    Face-face connectivity OK.
    Mesh topology OK.
    Boundary openness (-4.24201e-18 1.64967e-15 -3.01046e-18) Threshold = 1e-06 OK.
    Max cell openness = 1.05879e-16 OK.
    Max aspect ratio = 20.1795 OK.
    Minumum face area = 1e-06. Maximum face area = 4e-05.  Face area magnitudes OK.
    Min volume = 2e-08. Max volume = 4e-08.  Total volume = 0.00162.  Cell volumes OK.
    Mesh non-orthogonality Max: 0.514232 average: 0.00901361 Threshold = 70
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.00897456 OK.
    Mesh geometry OK.
Mesh OK.
ExecutionTime = 2.64 s  ClockTime = 3 s
Unfortunately, no layer addition/removal is happening and the whole compressor cylinder domain is moving.The attached Picture 0.png Shows the mesh at the end of the simulation (TDC). I do not know why this is happening. Could you please support me with my case setup? Please let me know if further Information is required.

Thanks in advance.
Attached Images
File Type: png -180.png (159.1 KB, 34 views)
File Type: jpg 0.jpg (57.5 KB, 21 views)
KotRainer11 is offline   Reply With Quote

Reply

Tags
dynamic mesh, topological changes


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
RP_Set_Integer does not work in parallel 86lolo Fluent UDF and Scheme Programming 2 July 3, 2014 11:37
Does CX_Interpret_String work in parallel? 86lolo Fluent UDF and Scheme Programming 2 June 30, 2014 04:36
Companies that lease software & hardware for cloud-based work? Catthan ANSYS 0 June 18, 2014 10:53
Why do the Plant library cases don't work? Alumna Phoenics 6 June 22, 2004 12:08
why my In-Form doesn't work? green Phoenics 2 May 27, 2004 21:03


All times are GMT -4. The time now is 11:46.