Community Profile

photo

Yanni


Last seen: circa 2 mesi fa Attivo dal 2022

Programming Languages:
MATLAB

Statistiche

  • Thankful Level 3
  • First Answer

Visualizza badge

Content Feed

Visto da

Domanda


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...

11 mesi fa | 1 risposta | 0

1

risposta

Domanda


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...

11 mesi fa | 1 risposta | 0

1

risposta

Domanda


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 ...

11 mesi fa | 1 risposta | 0

1

risposta

Domanda


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 ...

circa un anno fa | 1 risposta | 0

1

risposta

Domanda


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)...

circa un anno fa | 2 risposte | 0

2

risposte

Domanda


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...

circa un anno fa | 2 risposte | 0

2

risposte

Domanda


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)*...

circa un anno fa | 1 risposta | 0

1

risposta

Domanda


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...

oltre un anno fa | 1 risposta | 0

1

risposta

Domanda


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)...

oltre un anno fa | 1 risposta | 0

1

risposta