just to update the community and people interested in the problem.

I read that if and only if when one apply porous jump as the boundary condition, FLUENT would only consider flows perpendicular to the actual surface. That is definitely not our case. In particular, our system doesn’t possess inlet and/or outlet flows (or it doesn’t have to be perpendicular specifically), since it will be a real earth model simulation which will be representing fluid dynamics and it’s chemical interactions through 3-layered marine sediment. So, we are interested in natural flow of the water liquid in the undermost body, chemical interactions and outcomes of it in the middle body, once again natural flow of the water liquid and products from the reaction in the uppermost body by time, i.e. we solve the system as transient. As you can see, it’s not like a flow through cylindrical objects etc., thus it’s totally different. Today, we conclude that we need to figure out another option for boundary condition, since porous jump seems a bit non-functional.Once again the main fault we encounter when we run the program is most likely arised from enforcing a chemical reaction just in the middle body. Now, I am searching for some answers regarding how to enforce a chemical reaction just in a specific region(in our case in the middle) of a body.

My system is composed of 3 bodies, and the middle body will be the host of a reaction. All those three systems have the porous zone option activated. I want the water liquid which flows through the undermost body and arriving the porous jump boundary, passing through this boundary, reacting with the chemicals I designated in the middle body, creating those other products, and finally once again those products for this time passing through the second porous jump boundary and reaching to the uppermost body.

I hope you could imagine the system. Now I created 2 cell zones just above these porous jump boundaries with a width of 1 grid cell which was 8 m for my case, and also set the Y velocities of those interior cell zones as 0. Logically, since they are just above the porous jump boundaries they are actually in-situ of the reaction area(the middle body). So, I did activate the porous zone and reaction options for my cell zones, and set '0' value for the Y velocities in the fixed values tab (A researcher who encountered the same problem proposed this method, so I wanted to give it a try.).

Now, as for porous jump boundaries, I still set them as porous jump boundaries which I believe where the problem arises. When I didn't create a namedselections for them in mesh section, ANSYS Fluent gives the error for example:

Error: Species are not solved in neighbor cell threads (14 and 16)\n of interior zone 6. \n

Error Object: #f

but when I set them as porous jump boundaries, this time it gives an error of segmentation fault, no flow detected at porous jump 7(exp.).

One more thing, the properties of the porous jump condition:

FACE PERMEABILITY: 1e+10(default) How do I set it, I have a mesh system of 5000x5000 m and 40000 elements(grids) in this system

POROUS MEDIUM THICKNESS: 0 it is.

PRESSURE-JUMP COEFFICIENT: 0(default) I saw some calculations for cylindirical studies but I don't know how to calculate it or what are the values for 2D planar structures(systems).

What is it that I'm doing completely wrong and/or fail to notice?

Any comment would be much appreciated. Thanks in advance.

Cheers,

Selvi