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

Pressure Dirichlet value in BiCG STAB

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 24, 2016, 00:56
Post Pressure Dirichlet value in BiCG STAB
  #1
New Member
 
Jawahar
Join Date: Nov 2016
Posts: 4
Rep Power: 10
jawapasu is on a distinguished road
I am solving the pressure poisson equation in BiCG STAB, I am using finite element formulation. the matrix is like this AX=B. Where A is the matrix where each row represents the weight of the surrounding nodes of that node considered(can be taken as the center node which is unknown). X is the unknown node phi and the corresponding B is the known value of the grad times the provisional velocity.
Suppose if there are 9 nodes in 2D, the nodes be numbered from 1 to 9, let 9th node be the center of all other nodes, then the matrix A the row will look like,
[weight of node1 weight of node2 .. .. .. .. .. .. .. weight of node 8]
the X will be
[phi of 9]
and B will be
[grad times velocity]

The problem is when we consider the pressure nodes with known values(i.e., Dirichlet nodes) , then i remove the weight of the nodes from A and add that to B with the given dirichlet values multiplied with the weights of that nodes.
This makes the solution to diverge when i give any values other than zero(dirichlet), Can anyone help me how to give this constant dirichlet value in the matrix??
jawapasu is offline   Reply With Quote

Old   November 24, 2016, 04:41
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Are you using a library or is your own-made solver?
FMDenaro is offline   Reply With Quote

Old   November 24, 2016, 04:46
Default
  #3
New Member
 
Jawahar
Join Date: Nov 2016
Posts: 4
Rep Power: 10
jawapasu is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Are you using a library or is your own-made solver?
its custom made solver, it works fine when pressure dirichlet value is given as zero, but when i give some numerical value it diverges, and the preconditioner is ILU for the solver
jawapasu is offline   Reply With Quote

Old   November 24, 2016, 04:53
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Check if your solver converges with a source term such that the solution is constant, for example 1
FMDenaro is offline   Reply With Quote

Reply

Tags
dirichlet pressure, finite element method, iterative pressure solver, poisson equation

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
Wind tunnel Boundary Conditions in Fluent metmet FLUENT 6 October 30, 2019 13:23
question regarding LES of pipe flow - pimpleFoam Dan1788 OpenFOAM Running, Solving & CFD 37 December 26, 2017 15:42
Outlet pressure aja1345 FLUENT 5 August 18, 2015 12:02
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
Unsteady pressure differential between inlet and outlet of the pipe for single phase joshi20h FLUENT 0 September 26, 2012 13:41


All times are GMT -4. The time now is 20:43.