Community Profile

photo

Yanni


Last seen: 5 months ago Active since 2022

Programming Languages:
MATLAB

Statistics

  • Thankful Level 3
  • First Answer

View badges

Content Feed

View by

Question


How to use iteration and error for crank nicolson type to converge spatial and temporal discretization?
I'm using iteration and error for spatial discretization and temporal discretization respectively. I want apply iteration and er...

7 months ago | 1 answer | 0

1

answer

Question


How to obtain convergence of the curve to my boundary condtion?
I want a curve in each 'j'th iteration. But, I didn't get convergence of the curve to plot each 'j'th iterations between limit 1...

7 months ago | 1 answer | 0

1

answer

Question


How to stop time loop when steady state is reached?
I'm using unsteady case. so, it will reach a steady state at a certain time level. I fixed time at 'j'th column wise and it ran ...

7 months ago | 1 answer | 0

1

answer

Question


What is the steady state criteria in crank nicolson scheme, especially impilicit finite difference method? and How to stop time level when steady state is reached?
I'm dealing with unsteady case in my code. so, the CN term updated in each time step and it ran upto final time (tmax=100). but ...

8 months ago | 1 answer | 0

1

answer

Question


How to use indices in A matrix while conditonal and for loop condtion?
How to use indices in LHS of A matrix while using for and if, esleif and else condtions? clc clear all m=5;n=7; A=zeros(m,m)...

8 months ago | 2 answers | 0

2

answers

Question


How to change diagonal, subdiagonal and superdiagonal values with respect time while using loop and conditional statement?
xmax=1; ymax=7; m=20; n=29; dx=xmax/m; dy=ymax/n; dt=0.2; UOLD=zeros(m,n);VOLD=zeros(m,n); A=zeros([1,m]); B=A; C=A; whi...

8 months ago | 2 answers | 0

2

answers

Question


How to describe right hand side of energy equation in tridiagonal system? Didn't get right hand side values of tridiagonal matrix.
xmax=1; ymax=20; dx=0.05; dy=0.25; dt=0.02; m=xmax/dx; n=ymax/dy; U=zeros(m,n); V=zeros(m,n); for j =1:n D(i)=-U(i,j)*...

11 months ago | 1 answer | 0

1

answer

Question


There is a mistake while using condtional and for loop condtion.i got the value zero only. but i want values and here, A is subdiagonal, B is diagonal and C is superdiagonal.
xmax=1; ymax=20; dx=0.05; dy=0.25; dt=0.02;m=xmax/dx; n=ymax/dy; UOLD=zeros(m,n); VOLD=zeros(m,n); for j=1:n if i==2:m...

11 months ago | 1 answer | 0

1

answer

Question


here, i couldn't get the A,B,C values while using if, elseif conditions. how to get values while using conditional statement?
dx=0.05; dy=0.25; dt=0.2; for j=1:n if i==2:m C(i)=V(i,j)*dt/4*dy-dt/2*(dy)^2; elseif i==1:m-1 A(i)...

11 months ago | 1 answer | 0

1

answer