NEED HELP TO DEBUG THE FOLLOWING CODE: LID DRIVEN CAVITY USING SIMPLE APPROACH. REFER CFD BOOK BY JOHN D ANDERSON FOR THE ALGORITHM & THE STAGGERED GRID

$\%$ Two$-$dimensional Flow Solver with Pressure Correction Method

clc;

clear all;

close all;

$\%$ GEOMETRY CONSTANTS

L $=1\mathrm{.}0$; $\%$ length of domain, m

H $=1\mathrm{.}0$; $\%$ height of domain, m

Nxp $=129$; $\%$ # of points in x direction at which pressure is resolved $=$ No$.$ of columns

Nyp $=129$; $\%$ # of points in y direction at which pressure is resolved $=$ No$.$ of rows

Nxu $=$ Nxp $+1$; $\%$ # of points in x direction at which u velocity is resolved

Nyu $=$ Nyp; $\%$ # of points in y direction at which u velocity is resolved

Nxv $=$ Nxp $+2$; $\%$ # of points in x direction at which v velocity is resolved

Nyv $=$ Nyp $+1$; $\%$ # of points in y direction at which v velocity is resolved

dx $=$ L$/($Nxp$-1)$; $\%$ in x direction

dy $=$ H$/($Nyp$-1)$; $\%$ in y direction

dt $=1$e$-3$; $\%$ time resolution, seconds

alpha$\_$p $=0\mathrm{.}1$; $\%$ under$-$relaxation factor for pressure

alpha$\_$u $=0\mathrm{.}1$; $\%$ under$-$relaxation factor for u velocity

alpha$\_$v $=0\mathrm{.}1$; $\%$ under$-$relaxation factor for v velocity

$\%$ FLUID CONSTANTS

rho $=1\mathrm{.}0$; $\%$ density, kg$/$m$^3$

mu $=0\mathrm{.}1$; $\%$ dynamic viscosity, kg$/$m s

$\%$ Reynolds Number, Re $=100$

$\%$ INITIAL CONDITIONS

pstar $=$ ones$($Nyp$,$Nxp$)$;

pstar$1=$ ones$($Nyp$,$Nxp$)$;

ustar $=$ zeros$($Nyu$,$Nxu$)$;

vstar $=$ zeros$($Nyv$,$Nxv$)$;

rho$\_$ustar $=$ rho.$*$ustar;

rho$\_$vstar $=$ rho.$*$vstar;

pprime $=$ zeros$($Nyp$,$Nxp$)$;

pprime$1=$ zeros$($Nyp$,$Nxp$)$;

uprime $=$ zeros$($Nyu$,$Nxu$)$;

vprime $=$ zeros$($Nyv$,$Nxv$)$;

A $=$ zeros$($Nyu$,$Nxu$)$;

B $=$ zeros$($Nyv$,$Nxv$)$;

$\%$ BOUNDARY CONDITIONS

ustar$(1,$:$)=1$;

v $=$ vstar;

u $=$ ustar; $\%$ horizontal velocity, ft$/$s

$\%$ Set boundary conditions for corners

u$($Nyu$,1)=0$; $\%$ Bottom$-$left corner

u$($Nyu$,$Nxu$)=0$; $\%$ Bottom$-$right corner

v$(1,1)=0$; $\%$ Top$-$left corner

v$($Nyv$,1)=0$; $\%$ Bottom$-$left corner

v$(1,$Nxv$)=0$; $\%$ Top$-$right corner

v$($Nyv$,$Nxv$)=0$; $\%$ Bottom$-$right corner

$\%$ Convergence criterion

tolerance $=1$e$-16$;

$\%$ PRESSURE CORRECTION METHOD

b $=-$dt$/($dx$)^2$;

c $=-$dt$/($dy$)^2$;

a $=-2*($b $+$ c$)$;

d $=$ zeros$($Nyp$,$Nxp$)$; $\%$ mass flow term in continuity equation

for iter $=1$:$10000$

for i $=1$:Nxp$-1\%$ column, $1$ to $128$

for j $=2$:Nyu$-1\%$ row, $2$ to $128$

k $=$ Nyp$-(($Nyp$-$j$)+1)+1$; $\%2$ to $128$

kv $=$ Nyv$-(($Nyv$-$j$))+1$; $\%3$ to $129$

ku $=$ Nyu$-(($Nyu$-$j$)+1)+1$; $\%2$ to $128$

vbar$1=0\mathrm{.}5*($v$($kv$-1,$i$+1)+$v$($kv$-1,$i$+2))$;

vbar$2=0\mathrm{.}5*($v$($kv$,$i$+1)+$v$($kv$,$i$+2))$;

A$($ku$,$i$+1)=-(($rho$*($u$($ku$,$i$+2)^2)-$rho$*($u$($ku$,$i$)^2))/(2*$dx$)...$

$+(($rho$*$u$($ku$-1,$i$+1)*$vbar$1)-($rho$*$u$($ku$+1,$i$+1)*$vbar$2))/(2*$dy$))...$

$+$ mu$*(($u$($ku$,$i$+2)-2*$u$($ku$,$i$+1)+$u$($ku$,$i$))/($dx$^2)...$

$+($u$($ku$-1,$i$+1)-2*$u$($ku$,$i$+1)+$u$($ku$+1,$i$+1))/($dy$^2))$;

rho$\_$ustar$($ku$,$i$+1)=$ rho$\_$ustar$($ku$,$i$+1)+$ A$($ku$,$i$+1)*$dt $...$

$-($dt$/$dx$)*($pstar$($k$,$i$+1)-$ pstar$($k$,$i$))$;

if i $==1$

ustar$($ku$,$i$)=-$ustar$($ku$,$i$+1)$;

end

if i $==$ Nxp$-1$

ustar$($ku$,$i$+2)=-$ustar$($ku$,$i$+1)$;

end

u$($ku$,$i$+1)=$ ustar$($ku$,$i$+1)$;

end

end

u $=$ ustar;

for i $=2$:Nxu$-2\%$ column, $1$ to $128$

for j $=1$:Nyu$-1\%$ row, $1$ to $128$

k $=$ Nyp$-(($Nyp$-$j$))+1$; $\%2$ to $129$

kv $=$ Nyv$-(($Nyv$-$j$))+2$; $\%3$ to $130$

ku $=$ Nyu$-(($Nyu$-$j$))+1$; $\%2$ to $129$

ubar$1=0\mathrm{.}5*($u$($ku$,$i$+1)+$u$($ku$-1,$i$+1))$;

ubar$2=0\mathrm{.}5*($u$($ku$,$i$)+$u$($ku$-1,$i$))$;

B$($kv$-1,$i$+1)=-((($rho$*$v$($kv$-1,$i$+2)*$ubar$1)-($rho$*$v$($kv$-1,$i$)*$ubar$2))/(2*$dx$)...$

$+($rho$*($v$($kv$-2,$i$+1)^2)-$rho$*($v$($kv$,$i$+1)^2))/(2*$dy$))...$

$+$ mu$*(($v$($kv$-1,$i$+2)-2*$v$($kv$-1,$i$+1)+$v$($kv$-1,$i$))/($dx$^2)...$

$+($v$($kv$-2,$i$+1)-2*$v$($kv$-1,$i$+1)+$v$($kv$,$i$+1))/($dy$^2))$;

rho$\_$vstar$($kv$-1,$i$+1)=$ rho$\_$vstar$($kv$-1,$i$+1)+$ B$($kv$-1,$i$+1)*$dt $...$

$-($dt$/$dy$)*($pstar$($k$-1,$i$)-$ pstar$($k$,$i$))$;

vstar$($kv$-1,$i$+1)=$ rho$\_$vstar$($kv$-1,$i$+1)/$rho;

vstar$(1,$:$)=-$vstar$(2,$:$)$;

vstar$($Nyv$,$:$)=-$vstar$($Nyv$-1,$:$)$;

v$($kv$-1,$i$+1)=$ vstar$($kv$-1,$i$+1)$;

end

end

v $=$ vstar;

max$\_$divergence $=0$; $\%$ Initia

Step by Step Solution

There are 3 Steps involved in it

1 Expert Approved Answer

Step: 1 Unlock

Question Has Been Solved by an Expert!

Get step-by-step solutions from verified subject matter experts

Step: 2 Unlock

Step: 3 Unlock

Students Have Also Explored These Related Databases Questions!

Q:

answer the question clearly You are building a flight-control system for which a convincing safety case must be made. Would you assign the tasks of safety requirements engineering, test case...

Q:

BCO 313: Negotiation at Work Anne Dwyer: a.dwyer@euruni.edu Negotiation at Work \"Managers and executives negotiate constantly over issues as varied as hiring decisions and purchases, corporate...

Q:

Questions : Why did Palm computing's first product fail? Why should someone finance Palm Computing? Which financing option should Palm Computing choose? What is their cash requirement? Palm...

Q:

Why an international code of business ethics would be good for business Journal of Business Ethics; Dordrecht; Jan 1998; Larry R Smeltzer Marianne M Jennings Source (subtitle): Volume: Issue: Start...

Q:

Chapter One: Valuing Diversity R The wise are as rare as eagles that fly high in the sky. Bantu proverb Managing Workplace Diversity I Chapter One: Valuing Diversity VALUING DIVERSITY Chapter...

Q:

Hi, I need someone to do summary for the article I upload AUDITING: A JOURNAL OF PRACTICE & THEORY Vol. 28, No. 2 November 2009 pp. 1-34 American Accounting Association DOI: 10.2308 / aud.2009.28.2.1...

Q:

I am struggling to do the proper computations for some simple accounting. After i get this down, i will produce a paper on the importance of each number i solved for. Can someone please help me...

Q:

Answer question "Google in China" using image provided G C HB HB G X + X canvas.biola.edu/groups/39504/discussion_topics/545276 Update : Apps YouTube Netflix M Inbox (651) - devon... My Tasks - My...

Q:

GUESTIONS-For purposes of this questions, read the article,Inside Amazon.Address the following: Can Covey's teachings be reconciled with management practices at Amazon? If you were a consultant work...

Q:

Poll indicate that 60% of small businesses surveyed plan to spend as much or more in 2012 than they did in 2011 on online marketing (bit.ly/Oq5n26). Three in four (75%) small businesses have...

Q:

Find all the zeros of the function. (Enter your answers as a comma-separated list.) g(x)=x+ 10x + 22 X= x Write the polynomial as a product of linear factors. Use a graphing utility to verify your...

Q:

Why is using a predetermined overhead rate a preferred method of accounting for overhead costs? Group of answer choices Actual overhead costs can fluctuate from month to month Actual overhead cost...

Q:

15-4B securi 2 Prepare a taible that summarols a table that sammarizes (a) the realized gains and losses and (b) the unrealized gains on for long-term tions and events in its long-termm investments...

Previous Question Next Question