[deepsterblue]

Isn't dieselEngineFoam obsolete? That's what other forum posts seem to hint at. You might be stuck with a buggy version. You may also want to test out Hrv's branch at 1.6-ext, since he's been making a few changes to the source-tree recently.

You might want to stick to a simple setup before delving into sprays with topo-changes.

[Abhishekd18]

Ok. I am using OF 1.6-ext and it still has dieselEngineFoam solver. I am using the kind of mesh depicted in the presentations from Prof. Lucchini (Polimi). The part of mesh where spray is injected is kept stationary and spray hasn't penetrated beyond stationary mesh till the solver crash.

[Abhishekd18]

Hi,

I am trying to use mesquite solver for engine mesh motion. I am getting following error.
I downloaded the latest files from git repository https://github.com/smenon/dynamicTopoFvMesh and followed the instructions given in install.txt.


Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: mesquiteMotionSolver
--> FOAM Warning :
From function dlLibraryTable:pen(const dictionary& dict, const word& libsEntry, const TablePtr tablePtr)
in file /data1/install/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/dlLibraryTableTemplates.C at line 68
library "libmesquiteMotionSolver.so" did not introduce any new entries

Selecting quality metric: InverseMeanRatio
Selecting objective function: LPtoP
Selecting optimization algorithm: FeasibleNewton
Outer termination criterion (tcOuter) was not found. Using default values.
Time = -179


--> FOAM FATAL IO ERROR:
Unknown patchField type timeVaryingDisplacement for patch type patch

Valid patchField types are :

21
(
zeroGradient
processor
generic
global
wedge
surfaceDisplacement
value
fixedValue
slip
empty
oscillatingFixedValue
mixed
calculated
angularOscillatingDisplacement
oscillatingVelocity
symmetryPlane
oscillatingDisplacement
cyclic
angularOscillatingVelocity
surfaceSlipDisplacement
timeVaryingUniformFixedValue
)


file: /data/kur/cfd008/OpenFOAM/RUN/Mesquite/constant/dynamicMeshDict::mesquiteOptions::fixedValuePatche s:iston from line 94 to line 97.

From function PointPatchField<PatchField, Mesh, PointPatch, MatrixType, Type>::New(const PointPatch&, const DimensionedField<Type, Mesh>&, const dictionary&)
in file /data1/install/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/newPointPatchField.C at line 221.

FOAM exiting

[deepsterblue]

Make sure that the compiled library "libmesquiteMotionSolver.so" is the updated one from github (which should overwrite the existing version in OF-1.6-ext), because it seems like there's a conflict. Or maybe you have a duplicate libs entry in dynamicMeshDict / controlDict.

The timeVaryingDisplacement pointPatchField is a custom BC that I wrote, which compiles and writes to $FOAM_USER_LIBBIN. You need to add a line like this to system/controlDict:

libs
(
"libenginePointPatchFields.so"
);

[Abhishekd18]

Yes. It works perfectly in serial. But it takes a lot of time when I turn on dynamicTopoFvMesh. So, I tried to do it in parallel but having problems with parallel. I could not find parallelDynamicTopoFvMesh on the git hub as reported in earlier posts.

[deepsterblue]

The version on github should work in parallel. What's the problem exactly?

It would help if you provided more detail - like a stack-trace or a small test-case that can reproduce the problem.

[Abhishekd18]

Create time

Create dynamic mesh for time = 180

Selecting dynamicFvMesh dynamicTopoFvMesh
Selecting metric Knupp
Selecting motion solver: mesquiteMotionSolver
--> FOAM Warning :
From function dlLibraryTable:pen(const dictionary& dict, const word& libsEntry, const TablePtr tablePtr)
in file /opt/OpenFOAM/OpenFOAM-1.6-ext/src/OpenFOAM/lnInclude/dlLibraryTableTemplates.C at line 68
library "libmesquiteMotionSolver.so" did not introduce any new entries

Selecting quality metric: InverseMeanRatio
Selecting objective function: LPtoP
Selecting optimization algorithm: FeasibleNewton
Outer termination criterion (tcOuter) was not found. Using default values.
Time = 180.25
Solving for point motion: Initial residual: 1 Final residual: 0.0094732 No Iterations: 29
Solving for point motion: Initial residual: 1 Final residual: 0.009793 No Iterations: 55

~~~ Mesh Quality Statistics ~~~
Min: 0.313026
Max: 0.998508
Mean: 0.843209
Cells: 95154
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

[0] finement Progress: 82.3878% : Bisections: 1674, Collapses: 0, Total: 1674
[0]
[0] --> FOAM FATAL ERROR:
[0]
[0]
FOAM parallel run aborting
[0]
[0] Failed to match edge: 24274 :: (4349 627) using comparison edge: (2194 70)--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun has exited due to process rank 0 with PID 28484 on
node cfd011 exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------

[deepsterblue]

Can I get my hands on a test-case? Preferably on a time-step where I can reproduce the problem...

[Abhishekd18]

Hi,

I have attached the simple cylinder case. Geometry is given as .geo file (gmsh format) due to upload size limitation.

[deepsterblue]

I'm not really sure how to get a mesh out of the geometry file you've given me. Care to give a few tips?

[Abhishekd18]

Hi,

Sorry for that. I have sent you the whole setup with mesh in openFoam format on your gmail.

[deepsterblue]

Thanks for that. I was able to reproduce the problem, and I'll get back once I have a fix (and some time on my hands).

[sasanghomi]

Hi abminternet
I am running this case (simpleEngineStem) but I have an error..please help me

ERROR: In /home/sasan/OpenFOAM/OpenFOAM-1.6-ext/ThirdParty/rpmBuild/BUILD/ParaView-3.12.0/VTK/Filtering/vtkExecutive.cxx, line 756
vtkPVCompositeDataPipeline (0x1ffe600): Algorithm vtkPOpenFOAMReader(0x1e37230) returned failure for request: vtkInformation (0x2a41ed0)
Debug: Off
Modified Time: 362516
Reference Count: 1
Registered Events: (none)
Request: REQUEST_DATA
FROM_OUTPUT_PORT: 0
ALGORITHM_AFTER_FORWARD: 1
FORWARD_DIRECTION: 0

[u22]

Hi foamers!
Someone can write the guide step by step for perfoming calculation of diesel. (Read 7 pages of posts. But still didn't get it.)
1) combustion chamber without valves with moving piston. Cyclic BC or full 3D model.
2) diesel spray, evaporation
3) combustion

Which solvers should use?
Which version of OF (1.6 ext, 2.1) ?
Which dynamic mesh library should use?
Which tutorial cases should use for startup? Or some working cases...
Is it works in mpi parallel?

something like that http://files.nequam.se/scania_movie.gif
and http://www.foamcfd.org/Nabla/main/spray.html#spray

[sasanghomi]

I think best solver is sonicTurbDyMEngineFoam . It is available in OP-ext 1.6
I couldn't run it with mpirun.


sasan.

[u22]

So, i have searched forum and now i guess that for OpenFOAM 2.1.1 correct solver for diesel engine is sprayEngineFoam. Correct me please, if i wrong. But for this solver there is no tutorial case. Does anybody have working case for this solver? Share it please..

[u22]

So, a tutorial case for dieselEngineFoam had to do it myself. I took the Niklas Nordin case (http://files.nequam.se/dieselEngineFoam_scania.tgz) and modify it for 1.6-ext. Here is (without mesh and initial U field, use it from original)

[u22]

Originally Posted by sasanghomi

I think best solver is sonicTurbDyMEngineFoam . It is available in OP-ext 1.6
I couldn't run it with mpirun.


sasan.

Is sonicTurbDyMEngineFoam cold flow solver? It does not support combustion, i guess... Correct me please, if i wrong.

[sasanghomi]

yes,It's a solver for cold flow and It does not support combustion...

[sasanghomi]

Originally Posted by Peter_600

Works like a charm Thank you!

Hi peter
I have a question about simpleEngineStem.tar.gz case...this case uses PISO method and have relaxation factors....I think the results are not correct...please correct me...is my idea wrong?