|
[Sponsors] |
Simple explain for satisfying boundary conditions in Stam and Bridson papers |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
June 1, 2017, 17:10 |
Simple explain for satisfying boundary conditions in Stam and Bridson papers
|
#1 |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
I read Stable Fluids by Jos Stam and FLUID SIMULATION SIGGRAPH 2007 Course Notes by Robert Bridson. The boundary conditions in the two papers are different in my sight. I don't understand them. How are these two ways related to solving a linear system of PDEs with boundary conditions?
In Stam's paper, he uses a function called set_bnd to set boundaries, so is this function just setting the boundaries that I need in my matrix of linear systems of Pressure Poisson and Diffuse Equation? In Numerical Analysis, we don't modify the matrix with the boundary values, but we put the boundary condition either it's a Neumann or Dirichlet. I don't really understand what kind of boundaries is this set_bnd function. My thinking now is that it sets the boundary values which will be used to solve linear systems but what kind of boundaries? In Bridson's paper, I don't find any implementation about boundaries except for speaking about the normal components of velocity. Actually, I can't find any relation between this and Stam's way. In the section on making fluid incompressible, he says something about boundaries but I don't understand. He just satisfies boundary conditions on the boundary cells then he solves the system or he put the boundaries in the matrix? Actually I don't know how to satisfy u*n = u_solid * n which I think is the no stick boundaries for interacting between a fluid and solid boundaries. |
|
June 1, 2017, 17:17 |
|
#2 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Quote:
That's not correct...when you prescribe Neumann BC.s, some entries of the matrix are changed. |
||
June 1, 2017, 17:54 |
|
#3 |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
I mean that if we just have a dirichlet or Neumann condition like u'(0)=1,we can just evaluate u (0) and put it in the first row of the matrix
Another thing, when we just solve the pressure poisson equation, if we just solve at the slip boundary conditions u.n = u.solid , do we use the pressure normal derivative boundary condition (partial P / Partial n =0)? Or we just use this pressure condition in the no slip condition? |
|
June 1, 2017, 18:05 |
|
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
We recently discussed about BC.s for the pressure equation in this thread:
pressure boundary conditions for collocated grid for Navier-Stokes |
|
June 1, 2017, 22:06 |
|
#5 |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
I am reading but my question is a lot simpler. I just want to solve the pressure poisson equation with the slip boundary conditions. The condition for the velocity is u*n = usolid*n which I don't know how to satisfy in the simulation. Another condition for the pressure, will it be the normal derivative of pressure =0? Is it different from the collocation grids and mac grids(staggered grid)? Please I am a beginner and seeking your help please be patient with me and help although my questions are stupid.
|
|
June 1, 2017, 23:27 |
|
#6 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Quote:
For the pressure equation, to be well posed the problem, you prescribe: dp/dn=n.(v*-v) |
||
June 2, 2017, 01:04 |
|
#7 |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
Can you take a look at the paper of Bridson? He didn't apply this boundary you said. I don't know actually how to apply the boundary of the pressure in the equation pf pressure poisson. His way is totally confusing and I can't understand the usage of velocity solid boundaries. He doesn't talk about the pressure boundaries. Please if you have a source or an explanation to clarify all of this confusion about boundaries please send me it. Thanks and cheers. What do you mean that it's prescribed in the momentum?
|
|
June 2, 2017, 03:50 |
|
#8 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Quote:
He describes the bc.s at the wall in terms of Neumann condition, see the section at page 27. |
||
June 2, 2017, 13:22 |
|
#9 |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
||
June 2, 2017, 13:29 |
|
#10 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Quote:
No, the Neumann condition implies only the normal component of the velocity, that means physically permeable/not-permeable wall. The tangential condition is prescribed ONLY in the momentum equations. |
||
June 2, 2017, 15:25 |
|
#11 |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
I read in Stam's ways of implementing a Fluid solvers that we use this normal derivative boundary condition to solve the pressure poisson problem. And he said that we use this zero normal derivative of pressure. Actually I don't understand what you mean with prescribed only in the momentum. Thanks for your time and help cheers.
|
|
June 2, 2017, 15:31 |
|
#12 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Quote:
The pressure equation requires only one boundary condition that is prescribed in term of the normal component of the velocity (by meanso of the Neumann condition) to fulfill the divergence-free constraint. That does not imply that you fix also the tangential component of the velocity, this latter being prescribed when solving the momentum equation. Setting dp/dn=0 implies that n.v*=n.v_n+1 |
||
June 2, 2017, 17:30 |
|
#13 | |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
Quote:
in the code I am totally confused |
||
June 2, 2017, 17:37 |
|
#14 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
The Neumann BC is derived from the momentum equation written in the form of the Hodge decomposition:
v_n+1 + Grad phi = v* you have just to project along the vector unit n. But if you do not have a theoretical background about the numerical formulation for the incompressible form, I suggest to start with a good texbook of CFD. |
|
June 2, 2017, 17:53 |
|
#15 | ||
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
In this source here in section 38.2.4
http://developer.download.nvidia.com...gems_ch38.html He said in the boundary conditions section Quote:
Quote:
|
|||
June 2, 2017, 18:20 |
|
#16 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
I don't see mathematical difference in that... Sec.4 in the Bridson notes explains nothing else that a projection method based on the Hodge decomposition.
|
|
June 2, 2017, 18:57 |
|
#17 | |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
Quote:
Bridson uses a different notation of the pressure at boundaries. What's the difference ? |
||
June 2, 2017, 19:16 |
|
#18 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Prescribing homogeneous or non homogeneous Neumann bc for the pressure equation is largely discussed in literature. It has difference in terms of accuracy, you can see that dp/dn=0 is a buondary layer approximation.
However, both bc.s are mathematically correct for the Poisson problem to be well posed. Again, you cannot set the tangential velocity component in the pressure equation. Do not forget that the pressure equation is nothing else that the continuity equation. Eq.17 is prescribed for the momentum equation. |
|
June 2, 2017, 21:26 |
|
#19 | |
New Member
Join Date: Mar 2017
Posts: 16
Rep Power: 9 |
Quote:
|
||
June 3, 2017, 05:11 |
|
#20 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,849
Rep Power: 73 |
Ok, however the issue is not simple and requires some care. I suggest to read other references about the topic of (exact/approximate) projection methods.
|
|
|
|