Hello developer, I am a beginner about this software.I encountered difficulties in solving the following problem
The code is as follows
func u0=0;
func v0=0;
real rhol=1.46,rhog=0.58,delta=0.05;
func rho0=0.5*(rhol+rhog)-0.5*(rhol-rhog)*tanh((((x-1)^2+(y-1)^2)^(1/2)-0.4)/delta);
real T=10,dt=0.1,Re=100,We=1000;
real sigma=0.95;
mesh Th = square(30,30,[2x,2y]);
fespace Uh(Th,P2,periodic=[[1,x],[3,x],[2,y],[4,y]]); Uh u=u0,v=v0,rho=rho0,uu,vv,rhot,uold,vold,rhoold;
fespace Ph(Th,P1,periodic=[[1,x],[3,x],[2,y],[4,y]]);Ph p,pp;
problem cons(rho,rhot)=
int2d(Th)(rhorhot/dt)-int2d(Th)((convect([uold,vold],-dt,rhoold)rhot/dt))+int2d(Th)(rhodx(uold)rhot+rhody(vold)rhot);
problem vdw(p,pp)=
int2d(Th)((p-8rhosigma/(3-rho))+3rho^2)pp)
problem NSK([u,v],[uu,vv])=
int2d(Th)(rho(uuu/dt+vvv/dt))-int2d(Th)(rho((convect([uold,vold],-dt,uold)*uu+convect([uold,vold],-dt,vold)vv)/dt))
+int2d(Th)((1/Re)(((4/3)*dx(u)-(2/3)*dy(v))*dx(uu)+(dx(u)+dy(v))*dy(uu)+(dx(v)+dy(u))*dx(vv)+((4/3)*dy(v)-(2/3)*dx(u))*dy(vv))
-int2d(Th)((1/We)rho(dx(dx(dx(rho))+dy(dy(rho)))*uu+dy(dx(dx(rho))+dy(dy(rho)))*vv))
+int2d(Th)(dx(p)*uu+dy(p)*vv);
for(real t=0;t<T;t+=dt)
{
uold=u;
vold=v;
rhoold=rho;
cons;
vdw;
NSK;
plot([u,v,rho],wait=false,value=true);
}
The error is as follows
14 : problem vdw(p,pp)=
15 : int2d(Th)((p-8rhosigma/(3-rho)) error operator - <10LinearCombI7MGauche4C_F0E>,
List of choices
( : , )
( <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> : <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> )
( <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> : <14Matrice_CreuseISt7complexIdEE>, <14Matrice_CreuseISt7complexIdEE> )
( <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> : <14Matrice_CreuseISt7complexIdEE>, <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> )
( <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> : <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE>, <14Matrice_CreuseISt7complexIdEE> )
( <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> : <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE>, <NSt7__cxx114listISt5tupleIJSt7complexIdEP13VirtualMatrixIiS3_EbEESaIS7_EEE> )
( <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> : <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> )
( <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> : <14Matrice_CreuseIdE>, <14Matrice_CreuseIdE> )
( <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> : <14Matrice_CreuseIdE>, <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> )
( <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> : <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE>, <14Matrice_CreuseIdE> )
( <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> : <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE>, <NSt7__cxx114listISt5tupleIJdP13VirtualMatrixIidEbEESaIS5_EEE> )
( <10FormLinear> : <10FormLinear> )
( <12FormBilinear> : <12FormBilinear> )
( <6C_args> : <6C_args>, <10FormLinear> )
( <6C_args> : <6C_args>, <12FormBilinear> )
( <10LinearCombI7MGauche4C_F0E> : <10LinearCombI7MGauche4C_F0E>, <10LinearCombI7MGauche4C_F0E> )
( <10LinearCombI6MDroit4C_F0E> : <10LinearCombI6MDroit4C_F0E>, <10LinearCombI6MDroit4C_F0E> )
( <10LinearCombISt4pairI7MGauche6MDroitE4C_F0E> : <10LinearCombISt4pairI7MGauche6MDroitE4C_F0E>, <10LinearCombISt4pairI7MGauche6MDroitE4C_F0E> )
( <10LinearCombI7MGauche4C_F0E> : <10LinearCombI7MGauche4C_F0E> )
( <10LinearCombI6MDroit4C_F0E> : <10LinearCombI6MDroit4C_F0E> )
( <10LinearCombISt4pairI7MGauche6MDroitE4C_F0E> : <10LinearCombISt4pairI7MGauche6MDroitE4C_F0E> )
( : , )
( <4C_F0> : <12TransE_Array>, <12TransE_Array> )
( <4C_F0> : <7E_Array>, <12TransE_Array> )
( <4C_F0> : <7E_Array>, <7E_Array> )
( <4C_F0> : <12TransE_Array>, <7E_Array> )
( <12Add_Mulc_KN_ISt7complexIdEE> : <8Mulc_KN_ISt7complexIdEE>, <3KN_ISt7complexIdEE> )
( <12Add_Mulc_KN_ISt7complexIdEE> : <3KN_ISt7complexIdEE>, <8Mulc_KN_ISt7complexIdEE> )
( <12Add_Mulc_KN_ISt7complexIdEE> : <8Mulc_KN_ISt7complexIdEE>, <8Mulc_KN_ISt7complexIdEE> )
( <8Mulc_KN_ISt7complexIdEE> : <3KN_ISt7complexIdEE> )
( <7Sub_KN_ISt7complexIdEE> : <3KN_ISt7complexIdEE>, <3KN_ISt7complexIdEE> )
( <12Add_Mulc_KN_IdE> : <8Mulc_KN_IdE>, <3KN_IdE> )
( <12Add_Mulc_KN_IdE> : <3KN_IdE>, <8Mulc_KN_IdE> )
( <12Add_Mulc_KN_IdE> : <8Mulc_KN_IdE>, <8Mulc_KN_IdE> )
( <8Mulc_KN_IdE> : <3KN_IdE> )
( <7Sub_KN_IdE> : <3KN_IdE>, <3KN_IdE> )
( <12Add_Mulc_KN_IlE> : <8Mulc_KN_IlE>, <3KN_IlE> )
( <12Add_Mulc_KN_IlE> : <3KN_IlE>, <8Mulc_KN_IlE> )
( <12Add_Mulc_KN_IlE> : <8Mulc_KN_IlE>, <8Mulc_KN_IlE> )
( <8Mulc_KN_IlE> : <3KN_IlE> )
( <7Sub_KN_IlE> : <3KN_IlE>, <3KN_IlE> )
( : )
( : )
( : )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
( : , )
Error line number 15, in file vdw1.edp, before token )
current line = 15
Compile error :
line number :15, )
error Compile error :
line number :15, )
code = 1 mpirank: 0
How should I modify it