Forum: Main CFD Forum
September 8, 2011, 13:19
|
|
Replies: 72
Views: 44,962
Hi,
I have employed the Gauss-Seidel methd...
Hi,
I have employed the Gauss-Seidel methd for solving the pressure Poisson equation. For the boundary nodes, you have to take care of the coefficients of different terms according to the...
|
|
Forum: Main CFD Forum
September 7, 2011, 13:06
|
|
Replies: 72
Views: 44,962
Hello,
The modifications I made are the...
Hello,
The modifications I made are the definition of boundary conditions and the schemes for solving the viscosity. The boundary conditions are described previously, and the scheme has been...
|
|
Forum: Main CFD Forum
September 6, 2011, 14:44
|
|
Replies: 72
Views: 44,962
Well, I can't say it is solved. Since I think...
Well, I can't say it is solved. Since I think both explicit and implicit viscosity will work for this problem. However, it turns out both converge now, but the solutions from explicit method are way...
|
|
Forum: Main CFD Forum
September 6, 2011, 13:19
|
|
Replies: 72
Views: 44,962
|
|
Forum: Main CFD Forum
September 5, 2011, 13:12
|
|
Replies: 72
Views: 44,962
Hi Vincent,
I think I can answer your...
Hi Vincent,
I think I can answer your questions about U+Ubc. If you work through the code, you will find that Ubc is not actually for applying the boundary conditions in the normal sense, it is...
|
|
Forum: Main CFD Forum
August 31, 2011, 13:24
|
|
Replies: 72
Views: 44,962
|
|
Forum: Main CFD Forum
August 29, 2011, 12:45
|
|
Replies: 72
Views: 44,962
Hi Prasanth,
I was using the predictor and...
Hi Prasanth,
I was using the predictor and correct method. I have checked all the things I can think of that may result in the failure of convergence, and it turns out the boundary conditions are...
|
|
Forum: Main CFD Forum
August 23, 2011, 17:24
|
|
Replies: 72
Views: 44,962
|
|
Forum: Main CFD Forum
August 21, 2011, 15:42
|
|
Replies: 72
Views: 44,962
Thanks :)
I have noticed the issure about your...
Thanks :)
I have noticed the issure about your second tip before. When solving for the duct flow problem, the pressure goes up to a unreasonably huge value and the solver fails to converge. So...
|
|
Forum: Main CFD Forum
August 8, 2011, 10:24
|
|
Replies: 72
Views: 44,962
|
|
Forum: Main CFD Forum
August 5, 2011, 10:25
|
|
Replies: 72
Views: 44,962
2. plotm.m
%%% Run pcolm ...
2. plotm.m
%%% Run pcolm
[X,Y,XC,YC,FX,FY,U,V,P,F1,F2]=pcolm;
%%% Get size
NI=length(X); NJ=length(Y); sz=[NJ,NI];
%%% Quiver velocity field
figure;...
|
|
Forum: Main CFD Forum
August 5, 2011, 10:24
|
|
Replies: 72
Views: 44,962
function CALCP;
global LTIME LTEST
global IPR...
function CALCP;
global LTIME LTEST
global IPR JPR IU IV IP
global NI NJ NIM NJM NIJ LI
global X Y XC YC FX FY
global DENSIT VISC GDS DT DTR SMALL ALFA ULID
global U V P PP F1 F2 DPX DPY U0 V0...
|
|
Forum: Main CFD Forum
August 5, 2011, 10:23
|
|
Replies: 72
Views: 44,962
Here it is. there are 2 m-files.
1. pcolm.m...
Here it is. there are 2 m-files.
1. pcolm.m (it is splitted in to 2 parts, since too long)
%%% pcolm - a simplified version of pcol.f code (Peric, 1997)
%%%
%%% This version is writen for...
|
|
Forum: Main CFD Forum
August 4, 2011, 15:28
|
|
Replies: 72
Views: 44,962
Hi, I have found a MATLAB code online for the...
Hi, I have found a MATLAB code online for the lid-driven cavity problem, and trying to modify it to see if this code will work for a channel flow or Couette flow. Do you mind having a look at this...
|
|
Forum: Main CFD Forum
August 3, 2011, 23:34
|
|
Replies: 72
Views: 44,962
|
|
Forum: Main CFD Forum
August 2, 2011, 11:30
|
|
Replies: 72
Views: 44,962
Thanks, Vincent.
I am using the...
Thanks, Vincent.
I am using the fractional step method (predictor and corrector) on a staggered grid.
The boundary conditions I defined for the pipe flow are as follows:
Inlet: uniform u...
|
|
Forum: Main CFD Forum
July 6, 2011, 15:00
|
|
Replies: 72
Views: 44,962
Matlab code for pipe flow
Hi All,
I have written a MATLAB code for a 2-D lid driven cavity problem, and it works fast and well, the results are consistent with experimental data. But this code does not work for the duct...
|
