[luisfeliperojas95]

Good Nigth

My name is Luis Felipe Rojas Avila , Colombian Environmental Engineer. Right now I am finishing my Thesis Project. That's about simulate the DBO concentration across . I have to compile the algorithm is down this message. The issue I am having is ... ¿ Can anybody help me to find the cause for this error?

Following the description of the issue is

*The DBO Equation has been developed by me building the code in C++ Language , using a global balance with VOF method.

Is presented Here



The main factors of the equation are

phi = Flux in the main Finite Element Volume
U = Velocity on the Control Volume
k = Kinetic Constant of DBO
DBO = Biological Oxygen Demand
W = External Inlet of DBO

Then I made the principal code for my problem that is presented Here

*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.

OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.

You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.

Application
icoFoam

Description
Transient solver for incompressible, laminar flow of Newtonian fluids in terms of DBO balance over a geometry that contains water free flow.

\*---------------------------------------------------------------------------*/

#include "fvCFD.H"
#include "pisoControl.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"

pisoControl piso(mesh);

#include "createFields.H"
#include "initContinuityErrs.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nStarting time loop\n" << endl;

while (runTime.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

#include "CourantNo.H"

// Momentum predictor

fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
- fvm::laplacian(nu, U)
);

if (piso.momentumPredictor())
{
solve(UEqn == -fvc::grad(p));
}

// --- PISO loop
while (piso.correct())
{
volScalarField rAU(1.0/UEqn.A());
volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p));
surfaceScalarField phiHbyA
(
"phiHbyA",
fvc::flux(HbyA)
+ fvc::interpolate(rAU)*fvc::ddtCorr(U, phi)
);

adjustPhi(phiHbyA, U, p);

// Update the pressure BCs to ensure flux consistency
constrainPressure(p, U, phiHbyA, rAU);

// Non-orthogonal pressure corrector loop
while (piso.correctNonOrthogonal())
{
// Pressure corrector

fvScalarMatrix pEqn
(
fvm::laplacian(rAU, p) == fvc::div(phiHbyA)
);

pEqn.setReference(pRefCell, pRefValue);

pEqn.solve(mesh.solver(p.select(piso.finalInnerIte r())));

if (piso.finalNonOrthogonalIter())
{
phi = phiHbyA - pEqn.flux();
}

}

#include "continuityErrs.H"

U = HbyA - rAU*fvc::grad(p);
U.correctBoundaryConditions();

}

{

// --- DBO Solver

#include "DBOeq.H"
#include "Eeq.H"

fvScalarMatrix DBOeq

(
(((fvm::((phi)/(U))*DBO())-((fvm::(phi)/(U)))+((((fvm::((phi))/(U))*(DBO-DBO*k)+(fvm::ddt(DBO)))))))-(W)/fvm::(U) == (((((E)/(U))*(DBO-DBO*k)+((E)/(U))*(((DBO-DBO*k)))+(((S))))))

&nbsp;);

runTime.write();

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}

Info<< "End\n" << endl;

return 0;
}

// ************************************************** *********************** //

Then I present the Header files (Auxiliary Files)

//Solve the DBOEquation
//Takes the DBO Value for the General Balance Equation
//Leer datos de entrada

//Calcular la ecuacion general a partir de los datos que se tienen

fvScalarMatrix DBOeq
(
fvScalarMatrix DBOeq
(
// (((phi)/(U))*DBO())
//-(((phi)/(U))*ddt(DBO*k)-
//+((phi)/(U))*(DBO-DBO*k)
//+((Etotal)/(U))*(DBO-DBO*k)
//+((Etotal)/(U))*(DBO-DBO*k)
//-fvm::ddt(DBO)
//== fvOptions(DBO)


(((fvm::((phi)/(U))*DBO())-((fvm::(phi)/(U)))+((((fvm::((phi))/(U))*(DBO-DBO*k)+(fvm::ddt(DBO)))))))+(W)/(U) == (((((E)/(U))*(DBO-DBO*k)+((E)/(U))*(((DBO-DBO*k)))+(((S))))))

&nbsp;);

solve(DBOeq == fvOptions(DBO));

DBOeq.relax();

Then I have the Code for Eeq.H

Ep as double
Ucorr as double
En as double
Etotal as double
Dxm as double
Dx as double
S as double

//Dx == //Es el valor de Delta T para la Celda

{

//#include <checkMesh>

Dxm = min(magSf());

}

// Calcular factor de dispersion longitudinal

{

//#include <checkMesh.H>
//#include "fvCFD.H"

Dxm = min(magSf());
Dx = min(sqrt(magSf()));

Ep = 0.011*(((pow(U,2)*(pow(Dxm,2))/((Dx*U)));

}

//Calculo del factor de correccion para Ep

{

//#include <checkMesh>
//#include "fvCFD.H"

H = sqrt(magSf());
S = 0.50;

Ucorr = (pow(9.81*H*S,0.5));

}

{

//#include <checkMesh>
//#include "fvCFD.H"
//#include <iostream>

Dx = min(sqrt(magSf()));
Ucorr = Ucorr&

E = (Ucorr*Dx)/(2);

}

if (En < Ep);
{

E = Ep-En;

}

else
{

Etotal = 0

}

}

return 0;

Having all these the compiler wmake shows me the next information



The main problem I see is

*¿ What does mean in C++ error
''make: *** [Make/linux64Gcc46DPOpt/DBOQUAL2KFOAM.o] Error 1''?

*¿Is this Related with the fact that when I try to type the name of the Solver OpenFOAM doesn't return any result like this Image.



Because I want to run it , But any reason that doesn't let me work my solver.

I will Appreciate any help with my problem. I attach the Solver folder with the log file on a txt file.

[agustinvo]

Hola Luis Felipe,

it seems your log file is cut. Could you put it completelly? Usualy, in the first line of that bunch of lines you have the explanation of your problem.

[luisfeliperojas95]

Originally Posted by agustinvo

Hola Luis Felipe,

it seems your log file is cut. Could you put it completelly? Usualy, in the first
line of that bunch of lines you have the explanation of your problem.

Thanks for your help . Here is the Completely log file when I execute te compilation process. Thanks for your Help