CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Neumann boundary condition discretization

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 7, 2019, 19:16
Default Neumann boundary condition discretization
  #1
New Member
 
Karnauhov Valery
Join Date: Dec 2013
Posts: 20
Rep Power: 12
kveki is on a distinguished road
Hello. I apologize in advance for the Google translator.

I need help in understanding how in the FVM the discretization of the Neumann condition occurs at the boundary. I use the book F. Moukalled, L. Mangani, M. Darwish, The Finite Volume Method in Computational Fluid Dynamics. It contains a description of how, for example, to find pressure at the boundary of symmetry (fig.). But I do not quite understand how expressions (15.152) and 15.153) are obtained. Can someone explain to me how they are received? Are there any other books where the discretization of the boundary conditions is also painted?

Last edited by kveki; December 8, 2019 at 16:11.
kveki is offline   Reply With Quote

Old   December 7, 2019, 19:33
Default
  #2
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,192
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
I don't know what was your starting idiom but, chapeau to Google translator, never got such level of result.

Unfortunately, I don't have the book now, so can't be of any help until Monday.
sbaffini is offline   Reply With Quote

Old   December 8, 2019, 05:54
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,892
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by kveki View Post
Hello. I apologize in advance for the Google translator.

I need help in understanding how in the FVM the discretization of the Neumann condition occurs at the boundary. I use the book F. Moukalled, L. Mangani, M. Darwish, The Finite Volume Method in Computational Fluid Dynamics. It contains a description of how, for example, to find pressure at the boundary of symmetry (Fig. 1). But I do not quite understand how expressions (15.152) and 15.153) are obtained. Can someone explain to me how they are received? Are there any other books where the discretization of the boundary conditions is also painted?



I have the book, I just read that page and it seems the Eq.(15.153) is just a linear Taylor expansion along the direction d
FMDenaro is offline   Reply With Quote

Old   December 8, 2019, 16:25
Default
  #4
New Member
 
Karnauhov Valery
Join Date: Dec 2013
Posts: 20
Rep Power: 12
kveki is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
I have the book, I just read that page and it seems the Eq.(15.153) is just a linear Taylor expansion along the direction d
thanks for the answer.
Yes, expression (15.153) is obviously Teylor expansion. Although I don’t understand why in this case the gradient is used on the border, and not in the center of cell C. I still don’t understand why the extrapolation is carried out in the direction d _Cb, and not in the direction n.
My main question was how expression (15.152) was obtained. It is a difference of vectors, but I don’t understand how it was formed. From the gradient vector at point C, the directional derivative vector is subtracted, and as a result we get the gradient vector at the border?
kveki is offline   Reply With Quote

Old   December 8, 2019, 16:35
Default
  #5
New Member
 
Karnauhov Valery
Join Date: Dec 2013
Posts: 20
Rep Power: 12
kveki is on a distinguished road
Quote:
Originally Posted by sbaffini View Post
I don't know what was your starting idiom but, chapeau to Google translator, never got such level of result.

Unfortunately, I don't have the book now, so can't be of any help until Monday.
thanks for the answer.
attached a picture. I need to understand how the expression (15.152) is derived
kveki is offline   Reply With Quote

Old   December 8, 2019, 16:37
Default
  #6
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,892
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by kveki View Post
thanks for the answer.
Yes, expression (15.153) is obviously Teylor expansion. Although I don’t understand why in this case the gradient is used on the border, and not in the center of cell C. I still don’t understand why the extrapolation is carried out in the direction d _Cb, and not in the direction n.
My main question was how expression (15.152) was obtained. It is a difference of vectors, but I don’t understand how it was formed.



I don't have enough time to read the previous pages, in your case the geometry for the unstructured-based grid reconstruction depends on the colocation of the variables as well as the construction of the boundary values.
However, I would simplify the the issue, the symmetry for a scalar function can be prescribed by setting d phi /dn = n.Grad phi = 0 on the symmetry line. The value of of the function phi on the line of symmetry can be extrapolated by setting the problem of computing the unknown value for a known discrete derivatives, computed from the interior using an asymmetric stencil.
FMDenaro is offline   Reply With Quote

Old   December 8, 2019, 17:17
Default
  #7
New Member
 
Karnauhov Valery
Join Date: Dec 2013
Posts: 20
Rep Power: 12
kveki is on a distinguished road
thanks for the answer.
attached a picture. I need to understand how the expression (15.152) is derived
kveki is offline   Reply With Quote

Old   December 8, 2019, 17:31
Default
  #8
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,892
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by kveki View Post
thanks for the answer.
attached a picture. I need to understand how the expression (15.152) is derived



Start from the Eq.(15.151). The Eq.(15.152) is written in such a way that if you project that along the normal direction n, the 151 is fulfilled.
FMDenaro is offline   Reply With Quote

Old   December 9, 2019, 06:00
Default
  #9
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,192
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
Ok, probably the matter in the text is not stated in a sufficiently explicit manner for novices. Yet, the operations are quite clear. The very point here is that you should not really look for derivations, as there is no single correct way of applying bcs, as they all are approximations of the true physics.

For example, a very common approach is to simply set

a) p_b = p_c.

What they do in the text is to actually extrapolate the pressure gradient from c to b, and then constrain it to be parallel to the wall (eq. 15.152).

Then, the pressure at a symmetry boundary (subscript b) is extrapolated from the interior (subscript c) using this corrected gradient (eq. 15.153).

You can probably see that, with respect to my simpler formula in (a), eq. 15.153 takes into account the tangential pressure gradient ONLY when the center C is not along the boundary normal passing trough b. In all the other cases they are equal.

My personal taste is against using cell gradients for specifying boundary values whenever possible and resort instead to simpler formulas. Because cell gradients are computed using those boundary values in the first place. So, in order to correctly restart a simulation, you would need to save on file also gradients, which are going to waste a lot of space.

Of course, there can be tradeoffs. For example, one could use formula (a) for gradients and then use gradients for the true bc (like in eq. 15.153). But I dislike inconsistencies even more.
sbaffini is offline   Reply With Quote

Old   December 9, 2019, 15:22
Default
  #10
New Member
 
Karnauhov Valery
Join Date: Dec 2013
Posts: 20
Rep Power: 12
kveki is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Start from the Eq.(15.151). The Eq.(15.152) is written in such a way that if you project that along the normal direction n, the 151 is fulfilled.
Quote:
Originally Posted by sbaffini View Post
Ok, probably the matter in the text is not stated in a sufficiently explicit manner for novices. Yet, the operations are quite clear. The very point here is that you should not really look for derivations, as there is no single correct way of applying bcs, as they all are approximations of the true physics.

For example, a very common approach is to simply set

a) p_b = p_c.

What they do in the text is to actually extrapolate the pressure gradient from c to b, and then constrain it to be parallel to the wall (eq. 15.152).

Then, the pressure at a symmetry boundary (subscript b) is extrapolated from the interior (subscript c) using this corrected gradient (eq. 15.153).

You can probably see that, with respect to my simpler formula in (a), eq. 15.153 takes into account the tangential pressure gradient ONLY when the center C is not along the boundary normal passing trough b. In all the other cases they are equal.

My personal taste is against using cell gradients for specifying boundary values whenever possible and resort instead to simpler formulas. Because cell gradients are computed using those boundary values in the first place. So, in order to correctly restart a simulation, you would need to save on file also gradients, which are going to waste a lot of space.

Of course, there can be tradeoffs. For example, one could use formula (a) for gradients and then use gradients for the true bc (like in eq. 15.153). But I dislike inconsistencies even more.

I apologize for not immediately noticing and many thanks for the clarification. The expression (15.152) is similar to the expression (15.121) for a similar situation. We have grad pC, which has the normal and tangential components with respect to the symmetry boundary. At the boundary, the normal component of the gradient is zero, then only the tangential component remains, which is calculated based on the elementary difference of the vectors.
kveki is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Constant mass flow rate boundary condition sahm OpenFOAM 0 June 20, 2018 23:45
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 09:44
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 07:28
Low Mixing time Problem Mavier CFX 5 April 29, 2013 01:00
NEUMANN BOUNDARY CONDITION IN MESHLESS METHODS Kémelli Main CFD Forum 2 July 15, 2006 02:41


All times are GMT -4. The time now is 04:56.