code_saturne User's Forum

Dear Code_Saturne members,

We are currently working on a centrifugal blood pump. And we need the data where we should have the velocities (x,y,z) and pressures at each node of the mesh as well as the node reference for each of the elements of the mesh. Basically, we have to give our results in the form which is given here (mxyz and mien format): http://www.cats.rwth-aachen.de/about/so ... rmats/mixd
Do you know if this is possible to obtain in Code Saturne?

Best regards

Hello,

I looked at the format specification, and it seems it only handles one value for "nen" (number of element nodes), meaning it can only handle mesh with a single element type. If you're using the rotor-stator functionality of Code_Saturne, at least with the mesh joining option, you can expect to have more complex polyhedra at the joining interfaces (even if your initial mesh contains a single element type), so export to that format seems difficult.

To add output to a given format, you would need to add a "writer" (see src/fvm/fvm_writer* and src/fvm/fvm_to_ensight.c, which you clone and adapt).

A solution for multiple element types might be to use Code_Saturne's standard output, read it and tetrahedralize it under ParaView, and output it to a format supported by ParaView, which you could then convert to the mien format...

Regards,

Yvan