[Saleh Abuhanieh]

Dear Foamer,


The subject solver has been implemented in OpenFOAM v-1812 (overRhoPimpleDyMFoam).


I wonder if somebody was able to develop the same using foam-extend 4.1.


I tried without success until now, specially that no compressible solver example using overset is available in foam-extend 4.1.


Regards,
Saleh

[Saleh Abuhanieh]

Hi again


Any idea?


Regards
Saleh

[Zane]

Hi Saleh,


My issue is the same as yours. Have you solved it?

By the way, I only have experience implementing rhoPimDyMFoam in foam-extend 4.1.


I'd appreciate it if you could help me.

[Saleh Abuhanieh]

Hi Zhang,


What you are trying to solve?


Regards
Saleh

[Zane]

Dear Saleh,


Sorry for my late reply.

I am trying to simulate transient process of pump-turbine. To simulate water hammer, the change of density need to be considered, so I choose rhoPimpleFoam as the solver. In addition, to realize the movement of guide vanes, I added the code related to dynamicMesh into rhoPimpleFoam and got my own rhoPimpleDyMFoam in foam-extend-4.1. However, morphing mesh like displacementLaplacian(motionSolver) can not meet my requirements, because the movement of guide vanes is too large for morphing mesh. As a result, I choose to use oversetMesh now, but again, I am stuck due to the problem of parallel running of oversetMesh. I guess this problem comes from decompostion.

I am very glad to help you for what I know, and if you know how to solve my problem, please help me. Thanks in advance.


Best regards
Zhang

[Vivio]

Hi Zane,

Did you get an overset version of rhoPimpleFoam working even in serial?

I am trying to develop an overset version of sonicFoam - or any other compressible solver - in foam-extend and testing it with a simple 1D shock-tube case. To do this I looked at how interFoam compares to interOversetFoam, and modifying sonicFoam on the same basis.

Mostly this seems to consist of using oversetInterpolate(field) after solving each field, but I am missing something and this solver does not run.

However, I have a strange experience with foam-extend, which is that for multiphase I get better results using standard interFoam and including "liboversetMesh.so" in the controlDict, rather than using interOversetFoam. It still does the overset assembly automatically, and in particular there is a really bad pressure anomaly at the overset boundary with interOversetFoam, which is not present using interFoam even with an overset mesh. You can see this pressure anomaly in the dam break tutorial too.

However, I don't get how this works from the source code, since there is nothing in the solver that does the interpolation?? Also, if you try this approach with sonicFoam, it immediately throws a segmentation fault, so it's not universal.

[Zane]

Hi Vivio,


I'm sorry, I'm actually not familiar with the overset mesh in foam-extend. However, I think the overset mesh in openfoam.com (e.g. OpenFOAM-v2112) may help you. It has an overset version of rhoPimpleFoam.

As for your strange problem about interOversetFoam, I suggest you check the paths $FOAM_USER_APPBIN and $FOAM_USER_LIBBIN. Maybe that's bacause you accidentally called the wrong solver and libs. In other words, you ran a solver named "interFoam", but actually it was compiled into "interOversetFoam". In general, if a solver doesn't contain codes related to overset mesh, when you run it, it can not cope with overset mesh even including "liboversetMesh.so".


Best regards

[Saleh Abuhanieh]

Hi,


If the geometry is not so complex and the mesh size within 1M, you may try the overRhoPimpleFoam (ESI version). A better alternative can be the HiSA solver (free and openfoam based) which can utilize the ESI overset library.


You may check those:


https://hisa.gitlab.io/doc.html


https://journals.sagepub.com/doi/abs...44100221080771





Regards,
Saleh

[Vivio]

Thanks Saleh, I got the ESI version of overRhoPimpleDyMFoam working on the case and it works a treat, at least its comparable to standard rhoPimpleFoam.



Haven't tried HiSA yet but I will check it out and see if it gives better performance, the pressure-based solvers aren't the best when it comes to shock waves anyway.

[johku]

Originally Posted by Saleh Abuhanieh

- - If the geometry is not so complex and the mesh size within 1M, you may try the overRhoPimpleFoam (ESI version). A better alternative can be the HiSA solver (free and openfoam based) which can utilize the ESI overset library.

You may check those:

https://hisa.gitlab.io/doc.html

https://journals.sagepub.com/doi/abs...44100221080771
- -

Hi Saleh,

Thanks for the link to your publication. How do you enable the ESI overset libraries for the HiSA solver? I tried by adding the following to the system/controlDict and to the constant/dynamicFvMesh files for a static overset case as per the ESI-OpenFOAM tutorials:

controlDict

libs    (overset fvMotionSolvers);

dynamicMeshDict

dynamicFvMesh    dynamicOversetMesh
motionSolverLibs    (fvMotionSolvers);
solver    displacementLaplacian;
displacementLaplacianCoeffs
{
diffusivity    uniform 1;
}

These make the static case solution work "sort of": HiSA starts looping the (pseudo) time iterations as expected. However, it also creates the cell-to-cell addressing on each time iteration which I'm not really expecting. I would expect the cell-to-cell addressing to be done only on the first time step for a static case as in the default static overset mesh cases of the ESI-OpenFOAM. For a transient case with moving mesh, this would be acceptable. Am I missing something?

[Saleh Abuhanieh]

Hi,


Yes, for steady-state cases, the Overset Grid Assembly (OGA) algorithm shall work only at the first iteration.


Without going into the implementation details, for now, HiSA is suitable for the transient cases. For steadystate cases, without any code modifications, the OGA will work at each pseudo iteration.

I hope that was useful.

Regards,
Saleh

[applekiller]

Originally Posted by Saleh Abuhanieh

Hi,


Yes, for steady-state cases, the Overset Grid Assembly (OGA) algorithm shall work only at the first iteration.


Without going into the implementation details, for now, HiSA is suitable for the transient cases. For steadystate cases, without any code modifications, the OGA will work at each pseudo iteration.

I hope that was useful.

Regards,
Saleh

Hello, excuse me. I have been watching the Hisa solver recently， but I have encountered a question and I am looking forward to your answers.

I have a problem about code, I found this code in the file jacobianMatrix.H

typedef jacobianMatrix<2,1> compressibleJacobianMatrix;

What do the two numbers 2 and 1 in this code mean? Can this solver only solve two-dimensional problems?

Best wishes