|
[Sponsors] |
Flow Patterns Randomly Appearing In Transient Sim |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 15, 2017, 03:42 |
Flow Patterns Randomly Appearing In Transient Sim
|
#1 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
Good evening CFD-Online,
I am attempting to do some analysis on turbulent flow through a tube with baffling inside. Currently, I am meshing in Salome, and setting my case up in HELYX-OS. My issue is that beyond my inlet, at the first time step, every single baffle produces some sort of velocity. I have attached two pictures. One is at time 0, and one is at the immediate following time step. Regardless of inlet velocity, or time step, this issue occurs. I have attempted using DES-SA, in pimpleFoam, pisoFoam, and icoFoam, all producing the same issue. Please let me know what extra information I can provide. Thank you! -Dan |
|
June 15, 2017, 04:03 |
|
#2 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13 |
Well, from the looks of it this is probably a flow feature, so no suprise here. BUT you have to refine your mesh, this is way to coarse for a CFD sim. How many cells do you have?
|
|
June 15, 2017, 04:24 |
|
#3 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
I made the mesh coarser so I could run variations faster and see if I can solve the issue. Are you saying the coarse mesh is causing this problem? I am currently around 200,000 cells. I intend to be much higher once this is sorted.
Edit: I will go ahead and try a finer mesh before I call it a night |
|
June 15, 2017, 04:31 |
|
#4 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13 |
I would not call it an issue, because it's the flow you want to study. Around the baffles velocity must increase because of continuity, therefore you see transient velocity gradients around the edges = vorticies.
In regions with high velocity gradients e.g. your inlet, around the baffles etc. it's good practice to refine the mesh. In other regions it can be coarser. |
|
June 15, 2017, 04:35 |
|
#5 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
The continuity is where I have the problem. I have tried making my velocity incredibly small, and my time step incredibly narrow. However, the first instant of the simulation is always flow through the entire tube. Rather than starting at the inlet and marching forward, it immediately goes across the entire tube. It is as if each individual baffle behaves as an inlet.
Also let me thank you for the effort in replying! I hope I can understand this problem soon |
|
June 15, 2017, 04:46 |
|
#6 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13 |
Well, this is expected behaviour. I assume you have an incompressible fluid. So the very instant you calculate the flow, meaning the velocity, pressure field etc. so in the very instant of calculation, no matter how small or big the time step is, continuity must be obeyed in every area of the domain. There is nothing which needs to propagate first into the domain. Continuity means velocity increase in area's with small diameter. It's an instantenous effect.
Good luck! |
|
June 15, 2017, 04:52 |
|
#7 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
I expect the flow to obey continuity when it reaches the first baffle, however it appears to reach every single baffle and escape the tube at the first time step. I am working on a video that will hopefully highlight the issue better.
My expectations for this simulation were for the flow to begin at the inlet and hit one baffle at a time, and take a given amount of time to exit the tube. Hopefully I can have the video ready soon. The best example I can come up with is, if you are at a sink: The instant you turn on the sink, the water will work its way down and eventually hit the drain. In my simulation, it is as if the water immediately hits the drain |
|
June 15, 2017, 04:55 |
|
#8 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13 |
Ok, try to post the video. Maybe i'm missing something here..
|
|
June 15, 2017, 04:58 |
|
#9 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
Here is a link to a video.
You'll see a a large amount of velocity initially, then it will diminish at the first time step, and eventually drop to 0 as my simulation ends. I would expect the flow to begin at the inlet and slowly move over. No matter what I do, I cannot see the flow propagate through the tube. Edit: I recommend viewing the video at 0.25x speed |
|
June 15, 2017, 05:19 |
|
#10 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
Quick post: it is quite late where I am so I am gonna turn in for the night but I'll be back here tomorrow
I hope I can get this solved soon Sent from my iPhone using CFD Online Forum mobile app |
|
June 15, 2017, 05:51 |
|
#11 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13 |
This yt video is unavaible for me. Black Screen.
|
|
June 15, 2017, 05:53 |
|
#12 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
I apparently marked it as private on accident. It should be working now. If not, I will reupload tomorrow
Thank you again for your time tonight Sent from my iPhone using CFD Online Forum mobile app |
|
June 15, 2017, 05:56 |
|
#13 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13 |
Ok, can see it now and i see the problem. Was not clear to me from the previous posts. Strange and a bit difficult to tell what the issue here is exactly. It looks like your velocity inlet is temporary or sth.. Can you post BC's?
I'm off now too for a bit.. |
|
June 15, 2017, 15:58 |
Flow Patterns Randomly Appearing In Transient Sim
|
#14 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
For that video, my velocity inlet is time-varying. My intention is to turn the inlet on periodically to see the effect of multiple pulses on a system.
The problem still appears with a contant fixed value velocity. If I make the time step or velocity very small, the animation looks much like the first two image i posted. Here is a second video with a constant fixed value inlet For the rest of my boundary conditions, I have a fixed 0 pressure outlet, and everything else is either a symmetry or part of my geometry. Last edited by durndurn; June 16, 2017 at 14:52. |
|
June 19, 2017, 09:55 |
|
#15 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 13 |
Hi,
have you managed to solve the issue? I just had a look at your second video and i think i know what you mean. I looks like your velocity gradient at the inlet is moving into the domain. Are u sure you have a converged solution for every time step? On the other hand i'm still not sure this is can not be a physical transient effect... Maybe somebody else with experience can have a look and try to explain. |
|
June 19, 2017, 12:02 |
|
#16 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
I haven't solved it just yet.
The residuals look relatively clean. Pressure is a little oscillatory, but I would still call it convergence. I have a few more things to try and see if I can solve the problem. Thank you for your effort though! Sent from my iPhone using CFD Online Forum mobile app |
|
June 22, 2017, 03:34 |
|
#17 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
I experimented with boundary conditions and getting rid of some symmetry walls. Still no fix. Hopefully someone can offer some insight
|
|
June 22, 2017, 07:43 |
|
#18 |
Senior Member
|
Hi,
@durndurn Pattern looks like "non-orthogonality induced flow", so, since you have proposed it in the first message, please post the following: 1. checkMesh output. 2. fvSchemes & fvSolution files. 3. your boundary conditions. checkMesh could be posted inside CODE tags, for other files it would be more convenient if you create archive of fvSchemes, fvSolution files and 0 folder. |
|
June 22, 2017, 14:28 |
|
#19 |
New Member
Dan
Join Date: Jun 2017
Posts: 13
Rep Power: 8 |
Alexey,
The archive is attached and the checkMesh is below. Code:
Time = 0 Mesh stats points: 13945 faces: 135584 internal faces: 125064 cells: 65162 faces per cell: 4 boundary patches: 7 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 0 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 65162 polyhedra: 0 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology inlet 5 7 ok (non-closed singly connected) sym1 2708 1562 ok (non-closed singly connected) sym2 2499 1454 ok (non-closed singly connected) outlet 50 36 ok (non-closed singly connected) sym3 97 63 ok (non-closed singly connected) sym4 95 62 ok (non-closed singly connected) smallbarnobaf 5066 2851 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (-0.5 -9.714451465e-17 -0.00254) (0 0.5 0.75) Mesh has 3 geometric (non-empty/wedge) directions (1 1 1) Mesh has 3 solution (non-empty) directions (1 1 1) Boundary openness (-1.194181634e-16 1.069998123e-16 -1.230706011e-16) OK. Max cell openness = 2.202585184e-16 OK. Max aspect ratio = 8.897957786 OK. Minimum face area = 5.997818444e-07. Maximum face area = 0.01165816905. Face area magnitudes OK. Min volume = 2.336042024e-10. Max volume = 0.0003647460945. Total volume = 0.1874906724. Cell volumes OK. Mesh non-orthogonality Max: 57.78184399 average: 15.54904036 Non-orthogonality check OK. Face pyramids OK. Max skewness = 0.6432186553 OK. Coupled point location match (average 0) OK. Mesh OK. |
|
June 22, 2017, 16:29 |
|
#20 |
Senior Member
|
With 58 degrees of non-orthogonality you can:
- improve your mesh - use "leastSquares" instead of "Gauss linear" as gradScheme. I.e. fvSchemes file, sections gradSchemes. Do not know why HELYX-OS decided to put "Gauss linear", though it seems Engys guys are not happy with all grad schemes in OpenFOAM. - use 2 non-orthogonal correctors. fvSolutions, PISO dictionary, nNonOrthogonalCorrectors key. Start with two last modifications. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Issues on the simulation of high-speed compressible flow within turbomachinery | dowlee | OpenFOAM Running, Solving & CFD | 11 | August 6, 2021 06:40 |
initialization of flow field in a transient flow calculation | zhengjg | FLUENT | 1 | January 15, 2014 06:01 |
2D Transient flow using Fidap | Young | FLUENT | 1 | July 2, 2004 12:02 |
computation about flow around a yawed cone | Tylor Xie | Main CFD Forum | 0 | June 9, 1999 07:33 |
Flow patterns in a stirred tank | Glenn Price | Main CFD Forum | 13 | March 20, 1999 10:16 |