[martn]

Hi guys,

I use OF 1712 and want to simulate a engine motion with a moving piston & valves. The aim is to do it with topological changes of the mesh.
I created a case with the movingConeTopoFvMesh class. The mesh is build up like a moving piston in a cylinder which compresses air. The creation of a meshModifier file within the layerAdditionRemoval type complete my case. It works fine! The air is compressed like it should and the layers are removed (compression) and added (expansion).

Now I want to go further with the layerAdditionRemoval to the valves but without using the movingConeTopoFvMesh class. In the first steps I just replace the movingConeTopoFvMesh with the dynamicMotionSolverTopoFvMesh and the displacementInterpolation solver. So I am able to use a motion function for the moving parts. But despite using the meshModifier with the layerAdditionRemoval it is not working. No layers are chaniging. The cells are compressed by the motion.

Is this combination just not working or what should I change to use layerAR with another motion class except of the movingConeTopoFvMesh. The aim is to use a motion solver which I can predefine a motion function.

Best wishes,
Martn

[simrego]

Hi!


You should check the tutorial:
tutorials/multiphase/compressibleInterDyMFoam/laminar/sphereDrop
I think you will need something like that.


Here you can find an example video about the same motion solver +LAR + AMI (This is my dummy video, just made for some explanation earier...) :
https://www.youtube.com/watch?v=mcWCuTLolFg


But you can also write your own motion solver which will create the topology modifiers for you on the fly.

[martn]

tutorials/multiphase/compressibleInterDyMFoam/laminar/sphereDrop

This looks great!
Thank you!

I am just trying this tutorial but I want to define a motion in the pointDisplacement file with the type of "uniformInterpolatedDisplacement". I think this is the right type to predefine a motion with a table of time & coordinates but I do not know how to create this "table"-file.
I always get the error:

Did not find any times with table

I looked in the "uniformInterpolatedDiscplacementPointPatchVectorF ield.C" file

public fixedValuePointPatchField<vector>
{
    // Private data

        //- Name of displacement field
        const word fieldName_;

        const word interpolationScheme_;

        //- Times with pre-specified displacement
        wordList timeNames_;

        //- Times with pre-specified displacement
        scalarField timeVals_;

        //- User-specified interpolator
        autoPtr<interpolationWeights> interpolatorPtr_;


        //- Cached interpolation times
        labelList currentIndices_;

        //- Cached interpolation weights
        scalarField currentWeights_;

But I do not find out in which form the file with the time table must be. I used something like this:

(
(0.00000000 (0 0 0))
(0.00008333 (5.483e-06 0 0))
(0.00016667 (2.193e-05 0 0))
(0.00025000 (4.9337e-05 0 0))
(0.00033333 (8.7694e-05 0 0))
(0.00041667 (0.000136991 0 0))
(0.00050000 (0.000197212 0 0))
(0.00058333 (0.000268339 0 0))
(0.00066667 (0.00035035 0 0))
(0.00075000 (0.00044322 0 0))

with the time in the first column and the coordinates of the motion in the second column but it didn't work.
Does anybody know more about this?

[goeran]

Hey Martn,
just found this Thread. I actually can't help you and I assume you finished this work


But would it be possible, that you provide the whole case? I am working on a moving piston simulation as well, but I can't get the layerAdditionRemoval modifier to work as supossed.


Thanks!

[2538sukham]

The source code is your friend. The faces which will move are hard coded and the movement is specific for x direction only. It uses motion masking and moving the points nin the motion mask. I did layer removal which works fine but after some removals the face begins collapsing onto itself which should not occur. I can happily share the code for you to check if it works on your machine.
With regards