[Peter_600]

Hi

I try to set up an engine mesh with canted valves and piston layering. Therefore, I chose the thobois class in the engineTopoChangerMesh class.
I already succeed in running a case with the accordionEngineMesh class but with thobois class, i got the following error. I am running it on 1.6-ext.

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM Extend Project: Open source CFD |
| \\ / O peration | Version: 1.6-ext |
| \\ / A nd | Web: www.extend-project.de |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-ext-f3a01becc315
Exec : moveEngineTopoMeshPeter
Date : Jan 26 2011
Time : 19:53:24
Host : linux-ey53.site
PID : 8079
Case : /home/peter/OpenFOAM/peter-1.6-ext/run/meshMotion/Engine_one_Valve
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh

Selecting engineTopoChangerMesh thobois
// using new solver syntax:
rho
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
ft
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
fu
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
b
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
Xi
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
Su
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
h
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
hu
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

// using new solver syntax:
epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-06;
relTol 0;
}

Selecting motion solver: laplaceCellDecomposition
// using new solver syntax:
motionU
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
relTol 0.0001;
}

Selecting motion diffusivity: quadratic
Time = -540
void Foam::thoboisMesh::addZonesAndModifiers() : Zones and modifiers already present. Skipping.
deckHeight: 0.09
piston position: 0
Point zones found = 4
Face zones found = 2
Cell zones found = 2
Time = -539
bool Foam::layerSmooth::update()
mSolver.updateMesh(topoChangeMap())
deltaZ = 0 Piston at:0
pistonLayerID: 1
**Piston layering mode
setting valve motion b.c.
Valve n. 1 velocity = -0
Valve 0 lift: 1e-05 velocity: -0
valve motion boundary conditions set


--> FOAM FATAL ERROR:
Attempt to cast type fixedValue to type componentMixed

From function refCast<To>(From&)
in file /home/peter/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/typeInfo.H at line 115.

FOAM aborting

Aborted

Any ideas?
Remark: I had to change the addAttachDetachFacesThoboisMesh.H file with addAttachDetachFacesAccordionEngineMesh.H. Otherwise I couldn t create the meshModifiers. The thoboisMesh class looks very similar to the accordionEngineMesh class, but with piston layer.
Peter

[Peter_600]

I solved the problem.
I have changed in the motionU file the boundary conditions for the piston to componentMixed.
Does anyone knows what is meant by this bc§? I don t understand this boundary conditon. How to set it up?
Anyway the case isn t still working properly.

Peter

[abminternet]

Hi Peter,

I am trying now to do the same, tilted valves and piston layering. I am working on the thobois class but started also with the accordionEngine class following your posts, but I can't get it to work. What changes did u need to do? All the point, face and cell Zones for moving and static parts are created, but I am quite confused cause they dont get written as mesh modifiers in any way, so don't know if the code is complete or not :S any ideas? would appreciate it

[dreamilam]

Hi Peter,
I want to do a 3D engine simulation (cold flow) with canted valves. can u please share a tutorial?

[amrinder703]

Hello Peter

I am doing PhD in CFD in Mechanical Eng. Department, I am trying to simulate 4 stroke diesel engine with flat head and canted valves around 4 degree. I am trying to use thobois class in OpenFoam Ext 3.1
Am I using right class?
I think you have successfully implemented this class and i think no would better know about this class than you.

I have some queries about some keywords definitions?

detachFaces
movingPoints
movingInternalPoints
staticCells
movingCells
staticPointsName

any demo case or any kind of information would be really helpful.

Thanks
Amrinder Singh