|
[Sponsors] |
December 18, 2016, 20:40 |
Help with the compiling process
|
#1 |
New Member
Luis Felipe Rojas
Join Date: Aug 2015
Posts: 13
Rep Power: 11 |
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)))))) ); 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)))))) ); 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. |
|
December 19, 2016, 07:34 |
|
#2 |
Senior Member
Agustín Villa
Join Date: Apr 2013
Location: Alcorcón
Posts: 314
Rep Power: 15 |
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. |
|
December 19, 2016, 19:04 |
The Log File posted
|
#3 |
New Member
Luis Felipe Rojas
Join Date: Aug 2015
Posts: 13
Rep Power: 11 |
Thanks for your help . Here is the Completely log file when I execute te compilation process. Thanks for your Help
|
|
Tags |
bod, concentration, development, pisofoam, solver |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
The fl process could not be started because of UDF | majid_kamyab | Fluent UDF and Scheme Programming | 6 | December 15, 2015 09:42 |
[OpenFOAM.org] MPI compiling and version mismatch | pki | OpenFOAM Installation | 7 | June 15, 2015 17:21 |