photo

Yanni


Last seen: 8 Tage vor Aktiv seit 2022

Followers: 0   Following: 0

Programming Languages:
MATLAB

Statistik

  • Thankful Level 3
  • First Answer

Abzeichen anzeigen

Feeds

Anzeigen nach

Frage


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

etwa ein Jahr vor | 1 Antwort | 0

1

Antwort

Frage


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

etwa ein Jahr vor | 1 Antwort | 0

1

Antwort

Frage


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

mehr als ein Jahr vor | 1 Antwort | 0

1

Antwort

Frage


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

mehr als ein Jahr vor | 1 Antwort | 0

1

Antwort

Frage


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

mehr als ein Jahr vor | 2 Antworten | 0

2

Antworten

Frage


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

mehr als ein Jahr vor | 2 Antworten | 0

2

Antworten

Frage


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

mehr als ein Jahr vor | 1 Antwort | 0

1

Antwort

Frage


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

mehr als ein Jahr vor | 1 Antwort | 0

1

Antwort

Frage


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

mehr als ein Jahr vor | 1 Antwort | 0

1

Antwort