hello everyone please can some body helps me. I’m solving a set of pde and ode in fenics, but when I run a code I get the following error
au1=au0 + DT * (60/70)
^
SyntaxError: invalid syntax
I don’t know what is making this here is the sample of my code
def PressureUpdatePhase3(i,Fevolume):
#dp = 10./100
DT=(t[i]-t[i-1])*1000
iterations=0
Rc=10000
Rp=0.001
Ro=0.06
Rs=1
Rav=Ro/10
Rmv=Rp/10
Ps0=90;
Pc0=8
Qao0=80
Tsys0=0.5
kTsys=0
Tp=60/70
Tsys=Tsys0-kTsys/Tp
err=1
au0=0.000001
Cv=1
Cs=2.75
L0=0.005
Pp=10
while err>0.01:
iteration +=1
Pc1=Pc0+ DT *((Pp-Pv)/(Rv+Rmv)+(Ps0-Pv)/Rc + (Ps0-Pv)/(Rav+Ro))
Ps1= Ps0 + DT * (Qao0+ (Pv-Ps0)/Ro -Ps0/Ro)
Qao1=Qao0 + DT * (Pv-Ps0-Qao0*(Rav+Ro)
au1=au0 + DT * (60/70)
au0 = au0 % 1
if au0<= Tsys/Tp
Phi=(sin(pi * Tp * au0/Tsys))**2
else
Phi=0
pv= Phi * (Fevolume-16.3)+(1-Phi)*1.5*(exp(0.014*Fevolume)-1)
Pc0=Pc1
au0=au
Ps0=Ps1
Qao0=Qao1
print ('\n Error after ' ,iterations, ' iterations = ', err)
print ('\n PRESSURE UPDATE CONVERGED IN ', iterations , ' iterations \n ')
print ('\n Pressure is now', pv, ' \n ')
return Pv