[aussm]

Hello,


i am new to OpenFoam and i want to solve the equation:


fvm::ddt(rho) = fvc::curl(vel x rho)


*rho and vel are volVectorField's
*i want to solve it with euler time steps in the loops
*vel is always normal to rho (magnitude is constant =1)
*initial conditions of rho should look like a circle, so i want to set the initialFieldValues with setFields



i created a own solver by copying "scalarTransportFoam", what do i have to take care of to change my solver the right way?


i would be so happy about some help


LG aussm

[SHUBHAM9595]

Hey Manuel, I assume that rho is not density here...
Well, Just use the createFields.H of scalarTransportFoam and declare both the volVectorFields similar to "U".
And for the latter part where you want to solve for rho using fvm and fvc simultaneously at the same timestep... which is kind of difficult with the current syntax.....U should guide the solver to use the old values of rho in the fvc part.....So replace

fvc::curl(vel ^ rho);

with something like

fvc::curl(vel ^ rho.oldTime() );

and it should work.......