July 16, 2019 at 1:09 amWeiqiang LiuSubscriber
I am doing a porous soot combustion project. I reproduced the base case results in literature. The temperature magnitude is just the same with in the paper, though temperature distribution is slightly different. At first, I did not really pay attention to the difference of temperature distribution. I thought it was because we used different fluent solvers.
Then I continued to increase soot consuming rate like the author did. I found temperature would go extremely high at about 70 seconds and divergence happened. I checked my source code especially enthalpy source for several times and found no error. Then I tried to use default value of viscous resistance coefficients for porous media in fluent which is much smaller than the value recommended by the author. Then divergence never happens.
Therefore, I guess divergence in my case is because negative Y velocity namely filtration velocity is too small in my case and further convection is the only heat dissipation mechanism, so heat would accumulates in the corner and yields extremely high temperature.
I put filtration velocity contour in the paper below:
I put filtration velocity contour of my simulation below:
You can see the Y velocity contour is totally different. There is a very obvious negative Y velocity near the end wall of inlet channel which dissipates heat very quickly. However, for my case, the Y velocity is very small and uniform. I think this is the reason divergence happens when I increase soot consuming rate.
I checked almost every details of my model. I still can not get similar Y velocity contour.
Can anybody give me some suggestions?
July 16, 2019 at 4:54 amDrAmineAnsys EmployeeWhy new topic ?
Is the paper 2d or 3d?
What're the differences?
July 16, 2019 at 4:59 amDrAmineAnsys EmployeeMoreover it is working with default resistances and with the ones it shows no filtration and it diverged later. Have you already investigated if there values are correct?
July 16, 2019 at 1:21 pmWeiqiang LiuSubscriber
you're right. I think something wrong with my viscous resistance UDF. In fluent, F_PROFILE is used to store calculated viscous resistance value. I am wondering how fluent deal with the interface between two porous zones like the picture below:
The upper layer is porous soot cake while the lower layer is porous wall. I need to calculate viscous resistance coefficients for these two zones. I am wondering how fluent deal with the interface between them. I mean viscous resistance coefficient is a face type data, right?
July 16, 2019 at 2:11 pmRobAnsys Employee
In a porous media the resistance is on the cell. In a porous jump the resistance is on the face, but won't effect flow parallel to the face.
Between the two zones is an interior assuming you have a conformal mesh: that shouldn't do anything to the flow.
July 16, 2019 at 3:07 pmWeiqiang LiuSubscriber
So should I use C_PROFILE to store viscous value or F_PROFILE. Like you said, resistance is on the cell in porous zones.
July 16, 2019 at 3:31 pmDrAmineAnsys EmployeeYes as a C_profile.
July 16, 2019 at 8:38 pmWeiqiang LiuSubscriber
Thanks very much! I think I have a lot of problems of my case. Now I know that the Y-velocity contour in the paper is definitely right! Because soot was consumed at end of inlet channel at 70 seconds and most gas was forced to flow through this small opening, which yielded very large negative Y-velocity.
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Suppress Fluent to open with GUI while performing in journal file
- Floating point exception in Fluent
- What are the differences between CFX and Fluent?
- Heat transfer coefficient
- Getting graph and tabular data from result in workbench mechanical
- The solver failed with a non-zero exit code of : 2
- Difference between K-epsilon and K-omega Turbulence Model
- Time Step Size and Courant Number
- Mesh Interfaces in ANSYS FLUENT
- error in cfd post
© 2023 Copyright ANSYS, Inc. All rights reserved.