[13msmemusman]

Hi Foamers i was trying to implement following code in openFoam
Info<< nl << "Starting time loop" << endl;
while (runTime.loop())
{
#include "readPISOControls.H"
#include "readBPISOControls.H"
Info<< "Time = " << runTime.timeName() << nl << endl;
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"<< " ClockTime = " <<
runTime.elapsedClockTime() << " s"<< nl << endl;
#include "CourantNo.H"
{
volVectorField UcrosB= U ^ B;
volVectorField j=-sigma*fvc::grad(epo) + sigma*UcrosB;
volVectorField jcrosB = j ^ B;
fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
- (jcrosB/rho)
- fvm::laplacian(nu, U)
);
UEqn.relax();
solve(UEqn == -fvc::grad(p));
// --- PISO loop
for (int corr=0; corr<nCorr; corr++)
{
p.storePrevIter();
volScalarField rAU(1.0/UEqn.A());
U = rAU*UEqn.H();
phi = (fvc::interpolate(U) & mesh.Sf())
+ fvc::ddtPhiCorr(rAU, U, phi);
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
fvScalarMatrix pEqn
(
fvm::laplacian(rAU, p) == fvc::div(phi)
);
pEqn.setReference(pRefCell, pRefValue);
pEqn.solve();
if (nonOrth == nNonOrthCorr)
{
phi -= pEqn.flux();
}
p.relax();
}
#include "continuityErrs.H"
U -= rAU*fvc::grad(p);
U.correctBoundaryConditions();
}
}
// --- B-PISO loop
for (int Bcorr=0; Bcorr<nBcorr; Bcorr++)
{
volVectorField UcrosB= U ^ B;
fvScalarMatrix BEqn
(
fvm::laplacian(epo) - fvc::div(UcrosB)
);
BEqn.solve();
// #include "magneticFieldErr.H"
}
runTime.write();
}
Info<< "End\n" << endl;
return 0;
}



The Solver gives me following errors

mhdFoamModified.C:78:35: error: ‘epo’ was not declared in this scope
volVectorField j=-sigma*fvc::grad(epo) + sigma*UcrosB;
^
mhdFoamModified.C:96:3: error: ‘ddtPhiCorr’ is not a member of ‘Foam::fvc’
+ fvc::ddtPhiCorr(rAU, U, phi);
^
mhdFoamModified.C:122:16: error: ‘epo’ was not declared in this scope
fvm::laplacian(epo) - fvc::div(UcrosB)
^
In file included from mhdFoamModified.C:70:0:
/opt/openfoam231/src/finiteVolume/lnInclude/readPISOControls.H:3:15: warning: unused variable ‘nOuterCorr’ [-Wunused-variable]
const int nOuterCorr =
^
/opt/openfoam231/src/finiteVolume/lnInclude/readPISOControls.H:12:16: warning: unused variable ‘momentumPredictor’ [-Wunused-variable]
const bool momentumPredictor =
^
/opt/openfoam231/src/finiteVolume/lnInclude/readPISOControls.H:15:16: warning: unused variable ‘transonic’ [-Wunused-variable]
const bool transonic =
^
In file included from createFields.H:96:0,
from mhdFoamModified.C:62:
createPhiB.H:47:25: warning: unused variable ‘phiB’ [-Wunused-variable]
surfaceScalarField& phiB = *phiBPtr;
^
make: *** [Make/linuxGccDPOpt/mhdFoamModified.o] Error 1



Please help me i got this code from following link

http://www.woodruffscientific.com/files/hcll2Report.pdf