CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

engine simulation with mesh motion and topological changes

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

Like Tree3Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   February 15, 2012, 12:31
Default
  #121
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
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.
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   February 15, 2012, 12:54
Default
  #122
Member
 
Abhishek
Join Date: Dec 2010
Posts: 38
Rep Power: 6
Abhishekd18 is on a distinguished road
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 is offline   Reply With Quote

Old   April 16, 2012, 09:32
Default
  #123
Member
 
Abhishek
Join Date: Dec 2010
Posts: 38
Rep Power: 6
Abhishekd18 is on a distinguished road
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
Abhishekd18 is offline   Reply With Quote

Old   April 16, 2012, 10:49
Default
  #124
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
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"
);
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   April 18, 2012, 11:37
Default
  #125
Member
 
Abhishek
Join Date: Dec 2010
Posts: 38
Rep Power: 6
Abhishekd18 is on a distinguished road
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.
Abhishekd18 is offline   Reply With Quote

Old   April 18, 2012, 11:47
Default
  #126
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
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.
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   April 18, 2012, 22:40
Default
  #127
Member
 
Abhishek
Join Date: Dec 2010
Posts: 38
Rep Power: 6
Abhishekd18 is on a distinguished road
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).
--------------------------------------------------------------------------
Abhishekd18 is offline   Reply With Quote

Old   April 18, 2012, 23:11
Default
  #128
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
Can I get my hands on a test-case? Preferably on a time-step where I can reproduce the problem...
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   April 19, 2012, 08:37
Default
  #129
Member
 
Abhishek
Join Date: Dec 2010
Posts: 38
Rep Power: 6
Abhishekd18 is on a distinguished road
Hi,

I have attached the simple cylinder case. Geometry is given as .geo file (gmsh format) due to upload size limitation.
Attached Files
File Type: gz MesquiteSimpleCylinder.tar.gz (13.9 KB, 70 views)
Abhishekd18 is offline   Reply With Quote

Old   April 19, 2012, 08:48
Default
  #130
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
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?
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   April 19, 2012, 10:20
Default
  #131
Member
 
Abhishek
Join Date: Dec 2010
Posts: 38
Rep Power: 6
Abhishekd18 is on a distinguished road
Hi,

Sorry for that. I have sent you the whole setup with mesh in openFoam format on your gmail.
Abhishekd18 is offline   Reply With Quote

Old   April 19, 2012, 10:36
Default
  #132
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 386
Rep Power: 15
deepsterblue will become famous soon enough
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).
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   July 11, 2012, 08:20
Default
  #133
New Member
 
Join Date: Feb 2012
Posts: 12
Rep Power: 5
daxindamix is on a distinguished road
Hello,

does anyone has an example running with the last version of OpenFOAM? I download files on this forum but I never achieve to run it. There is always an error, probably because the functions change from version to version.

thanks!
daxindamix is offline   Reply With Quote

Old   September 23, 2012, 04:59
Default
  #134
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Iran - Tehran
Posts: 185
Rep Power: 4
sasanghomi is on a distinguished road
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

Last edited by sasanghomi; September 27, 2012 at 15:52.
sasanghomi is offline   Reply With Quote

Old   October 17, 2012, 23:19
Default
  #135
u22
Member
 
Anthony Nitski
Join Date: Aug 2009
Location: Earth
Posts: 34
Rep Power: 8
u22 is on a distinguished road
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
u22 is offline   Reply With Quote

Old   October 18, 2012, 06:58
Default Hi
  #136
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Iran - Tehran
Posts: 185
Rep Power: 4
sasanghomi is on a distinguished road
I think best solver is sonicTurbDyMEngineFoam . It is available in OP-ext 1.6
I couldn't run it with mpirun.


sasan.
sasanghomi is offline   Reply With Quote

Old   November 16, 2012, 01:23
Default
  #137
u22
Member
 
Anthony Nitski
Join Date: Aug 2009
Location: Earth
Posts: 34
Rep Power: 8
u22 is on a distinguished road
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 is offline   Reply With Quote

Old   November 19, 2012, 00:21
Default dieselEngineFoam tutorial for 1.6-ext
  #138
u22
Member
 
Anthony Nitski
Join Date: Aug 2009
Location: Earth
Posts: 34
Rep Power: 8
u22 is on a distinguished road
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)
Attached Files
File Type: zip scania.1.6-ext.zip (95.9 KB, 48 views)
u22 is offline   Reply With Quote

Old   November 19, 2012, 00:31
Default
  #139
u22
Member
 
Anthony Nitski
Join Date: Aug 2009
Location: Earth
Posts: 34
Rep Power: 8
u22 is on a distinguished road
Quote:
Originally Posted by sasanghomi View Post
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.
u22 is offline   Reply With Quote

Old   January 23, 2013, 08:49
Default
  #140
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Iran - Tehran
Posts: 185
Rep Power: 4
sasanghomi is on a distinguished road
yes,It's a solver for cold flow and It does not support combustion...
sasanghomi is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


LinkBacks (?)
LinkBack to this Thread: http://www.cfd-online.com/Forums/openfoam/83177-engine-simulation-mesh-motion-topological-changes.html
Posted By For Type Date
Untitled document This thread Refback February 4, 2014 12:36

Similar Threads
Thread Thread Starter Forum Replies Last Post
Dynamic moving mesh Pei-Ying Hsieh (Hsieh) OpenFOAM Running, Solving & CFD 64 June 7, 2012 10:04
engine simulation with mesh motion and topological changes abminternet OpenFOAM 0 December 16, 2010 12:47
Good mesh for pistoncylinder application Serkan Cetin OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 4 November 3, 2010 08:36
Radiation and miscellaneous enhancements vtk_fan OpenFOAM Running, Solving & CFD 6 February 18, 2008 00:49
Valve action Hrvoje Jasak (Hjasak) OpenFOAM Running, Solving & CFD 0 January 13, 2005 14:23


All times are GMT -4. The time now is 12:22.