CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Issue on zeroGradient BC for system rotation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 21, 2017, 19:26
Default Issue on zeroGradient BC for system rotation
  #1
New Member
 
eco
Join Date: Jan 2017
Posts: 5
Rep Power: 9
magnushaese is on a distinguished road
Dear formers,

i have an issue with the zeroGradient BC at flow outlet combined with system rotation for a 2-D (xy) incompressible steady state problem. The whole domain is about to be rotated among the z-axis.

1. Whats working so far:
Plane-turbulent-channel flow setup, good results without system rotation using either SRFSimpleFoam or a modified simpleFoam with Coriolis force in the momentum equation. Both solvers work fine with arbitrary two-equation-turbulence-models (e.g. kEpsilon, kOmegaSST) and give the same results.

BCs velocity/relative velocity:
inlet = fixed value
outlet = zeroGradient
BCs pressure:
inlet = zeroGradient
outlet = fixed value 0

2. The issue:
Starting from a fully converged non rotation state, as soon as the system is rotated among the the z-axis in positive or negative direction the outflow shows unphysical behavior, not holding the condition of zeroGradient for the relative velocity (and absolute velocity plus turbulence quantities, whereas the pressure is fine). The result is the same for the modified simpleFoam and SRFSimpleFoam.

(Picture see attachment)
Top case, positive rotation among z-axis.
Bot case, negative though.

3. My comment:
The issue seems to be independent of the solver and turbulence model (since the used turb-models dont even represent rotational dependency anyways). Furthermore, the picture illustrates the impact of the rotational direction which shows that there might be a problem with the understanding of absolute and relative velocity at the outlet BC.
As far as i understand the "type" of a BC is independent of the basic flow quantities such as U,p,k,omega what ever, it just "zeros" the normal gradient of the quantity it is applied to, right? In my understanding this should not be the problem in the case of "Urel" though.

I am tankful for every advice!

Best regards,

Magnus
Attached Images
File Type: png zoomed-view-outlet.png (42.8 KB, 42 views)

Last edited by magnushaese; February 21, 2017 at 19:34. Reason: highlighting the problem with color...
magnushaese is offline   Reply With Quote

Old   February 28, 2017, 07:54
Default
  #2
New Member
 
eco
Join Date: Jan 2017
Posts: 5
Rep Power: 9
magnushaese is on a distinguished road
solved: new boundary condition with the pressure gradient balances the Coriolis force.

grad(P) = -2 Omega ^ U

Last edited by magnushaese; March 8, 2017 at 11:47.
magnushaese is offline   Reply With Quote

Old   January 28, 2022, 14:43
Default
  #3
Member
 
George
Join Date: Dec 2020
Posts: 31
Rep Power: 6
damon707 is on a distinguished road
Hello magnushaese,

Did you use SRFSimpleFoam for this simulation?

Best,
George
damon707 is offline   Reply With Quote

Reply

Tags
boundary condition u, simplefoam coriolis, srfsimplefoam, zerogradient u


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
FOAM FATAL ERROR Maximum number of iterations exceeded liugx212 OpenFOAM Running, Solving & CFD 48 February 13, 2024 22:42
FATAL ERROR: Max number of iterations exceeded hbrist7 OpenFOAM Running, Solving & CFD 3 May 13, 2022 15:55
interFoam two-phase pipe flow air phase behaviour katete OpenFOAM Running, Solving & CFD 11 February 3, 2021 04:14
time step continuity problem in VAWT simulation lpz_michele OpenFOAM Running, Solving & CFD 5 February 22, 2018 20:50
Declaration of the function TH() Tobi OpenFOAM Programming & Development 10 July 5, 2016 16:17


All times are GMT -4. The time now is 13:09.