 |
 |
 |
|
|
|
[Sponsors] | ||||
simpleFoam Validation in Urban Environment using AIJ guidelines (openCAE) |
 |
|
|
LinkBack | Thread Tools | Display Modes |
March 20, 2013, 21:53
|
simpleFoam Validation in Urban Environment using AIJ guidelines (openCAE)
|
#1 |
|
Member
Jose Rey
Join Date: Oct 2012
Posts: 84
Rep Power: 6  |
Hi,
I've been trying to run Case-C of the AIJ-PWEAB test cases. It consists of a model of flow of air around a simple building. The model utilizes simpleFoam to solve it. I am using openFoam-2.1.1. I got over a few hurdles to get it to work, but now that it runs, the sequence breaks at iteration #2. The model can be found at: Case C = Simple Building Block. This is what the terminal shows when running the Allrun script : Code:
Running ./makeMesh on /home/admin1/AIJ-PWEAB/CaseC Running setDiscreteFields on /home/admin1/AIJ-PWEAB/CaseC Running simpleFoam on /home/admin1/AIJ-PWEAB/CaseC /opt/openfoam211/bin/tools/RunFunctions: line 37: 4194 Floating point exception(core dumped) $APP_RUN $* > log.$APP_NAME 2>&1 Running foamCalc on /home/admin1/AIJ-PWEAB/CaseC Running foamLog on /home/admin1/AIJ-PWEAB/CaseC Running sample on /home/admin1/AIJ-PWEAB/CaseC Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : simpleFoam
Date : Mar 20 2013
Time : 21:20:32
Host : "admin1-VirtualBox"
PID : 4194
Case : /home/admin1/AIJ-PWEAB/CaseC
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading field p
Reading field U
Reading/calculating face flux field phi
Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
sigmaEps 1.3;
}
No field sources present
SIMPLE: convergence criteria
field p tolerance 0.001
field "(U|k|epsilon|omega)" tolerance 0.0001
Starting time loop
Time = 1
DILUPBiCG: Solving for Ux, Initial residual = 0.375188, Final residual = 0.0147972, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0487866, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.0542991, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.00946159, No Iterations 2
time step continuity errors : sum local = 0.429902, global = -5.28454e-17, cumulative = -5.28454e-17
DILUPBiCG: Solving for epsilon, Initial residual = 0.0471403, Final residual = 0.000449379, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.2091, Final residual = 0.00420194, No Iterations 1
ExecutionTime = 0.22 s ClockTime = 1 s
Time = 2
DILUPBiCG: Solving for Ux, Initial residual = 0.837893, Final residual = 0.035823, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.732454, Final residual = 0.0340259, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.726034, Final residual = 0.034602, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::DICSmoother::DICSmoother(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::lduMatrix::smoother::addsymMatrixConstructorToTable<Foam::DICSmoother>::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::lduMatrix::smoother::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::GAMGSolver::initVcycle(Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::lduMatrix::smoother>&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#8 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/simpleFoam"
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/simpleFoam"
|
|
|
|
|
|
March 21, 2013, 03:18
|
|
#2 | |
|
New Member
Join Date: Nov 2012
Posts: 2
Rep Power: 0 |
You should check carefully the log files produced, especially the mesh part. There is a file setSet.batch missing in this directory. I copied the file from Case-A directory then everything is fine. I think the maintainer should fix this problem since I was puzzled by this problem for a long time.
Quote:
|
||
|
|
|
||
|
March 21, 2013, 07:11
|
|
#3 |
|
Member
Jose Rey
Join Date: Oct 2012
Posts: 84
Rep Power: 6 |
It Worked !! Thanks !!
You were exactly right. The file that is missing in the "setSet.batch". The log.MakeMesh complaints about this 1/3 of the way down: Code:
set points
end points
set face - owner - neigubour
end face - owner - neigubour
set boundary
make Mesh end
end
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : insideCells constant/triSurface/building.stl building
Date : Mar 20 2013
Time : 21:20:24
Host : "admin1-VirtualBox"
PID : 4187
Case : /home/admin1/AIJ-PWEAB/CaseC
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = 0
Reading surface from "constant/triSurface/building.stl"
Selected 7056 of 168912 cells
Writing selected cells to cellSet building
Use this cellSet e.g. with subsetMesh :
subsetMesh building
End
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : setSet -constant -batch setSet.batch
Date : Mar 20 2013
Time : 21:20:27
Host : "admin1-VirtualBox"
PID : 4188
Case : /home/admin1/AIJ-PWEAB/CaseC
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = constant
Time:constant cells:168912 faces:516498 points:178850 patches:6 bb:(-1 -1.5 0) (1.5 1.5 1.8)
cellSets:
building size:7056
Time = constant
mesh not changed.
Reading commands from file "setSet.batch"
--> FOAM FATAL ERROR:
Cannot open file "setSet.batch"
From function setSet
in file setSet.C at line 890.
FOAM exiting
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : subsetMesh building -overwrite
Date : Mar 20 2013
Time : 21:20:29
Host : "admin1-VirtualBox"
PID : 4189
Case : /home/admin1/AIJ-PWEAB/CaseC
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Reading cell set from building
Adding exposed internal faces to a patch called "oldInternalFaces" (created if necessary)
--> FOAM Serious Error :
From function IOobject::readHeader(Istream&)
in file db/IOobject/IOobjectReadHeader.C at line 89
Reading "/home/admin1/AIJ-PWEAB/CaseC/constant/index.html" at line 1
First token could not be read or is not the keyword 'FoamFile'
Check header is of the form:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class IOobject;
location "constant";
object index.html;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Writing subsetted mesh and fields to time 0
End
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : changeDictionary
Date : Mar 20 2013
Time : 21:20:32
Host : "admin1-VirtualBox"
PID : 4190
Case : /home/admin1/AIJ-PWEAB/CaseC
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
Read dictionary changeDictionaryDict with replacements for dictionaries
1
(
boundary
)
Replacing entries in dictionary boundary
Special handling of boundary as polyMesh/boundary file.
Loaded dictionary boundary with entries
7
(
x_
_x
y_
_y
z_
_z
oldInternalFaces
)
Merging entries from
2
(
z_
oldInternalFaces
)
fieldDict:
{
x_
{
type patch;
nFaces 0;
startFace 18872;
}
_x
{
type patch;
nFaces 0;
startFace 18872;
}
y_
{
type patch;
nFaces 0;
startFace 18872;
}
_y
{
type patch;
nFaces 0;
startFace 18872;
}
z_
{
type wall;
nFaces 784;
startFace 18872;
}
_z
{
type patch;
nFaces 0;
startFace 19656;
}
oldInternalFaces
{
type wall;
nFaces 3808;
startFace 19656;
}
}
Writing modified fieldDict boundary
End
|
|
|
|
|
|
|
March 21, 2013, 16:49
|
Linux Question: Link Command
|
#4 |
|
Member
Jose Rey
Join Date: Oct 2012
Posts: 84
Rep Power: 6 |
I am now trying to run the Case B which consists of a validation of several turbulent models for wind around a single building. When you call the Allrun script, it goes through an array of values that direct it to each directory. Each Allrun script for each model contains the following code:
Code:
link ../share/Allrun Code:
link: missing operand after `../share/Allrun' Try `link --help' for more information. The complete Case B is at the SVN repository here: http://www.opencae.jp/svn/OpenFOAM-V...J-PWEAB/trunk/ This is the directory structure of Case B with the files highlighted: . ├── Allclean ├── Allrun ├── index.html ├── kEpsilon │** ├── 0 │** ├── Allclean │** ├── Allrun │** ├── constant │** ├── index.html │** └── system ├── LICENSE.GPL2 ├── NonlinearKEShih │** ├── 0 │** ├── Allclean │** ├── Allrun │** ├── constant │** ├── index.html │** └── system ├── README ├── realizableKE │** ├── 0 │** ├── Allclean │** ├── Allrun │** ├── constant │** ├── index.html │** └── system ├── RNGkEpsilon │** ├── 0 │** ├── Allclean │** ├── Allrun │** ├── constant │** ├── index.html │** └── system └── share ├── 0 ├── Allclean ├── Allrun <-- this one contains the juice ├── box.dat ├── caseBMesh.foam ├── constant ├── index.html ├── log.cellSet ├── log.makeMesh ├── log.setDiscreteFields ├── makeMesh ├── makeStructuredGridMesh.py ├── measured ├── org0 ├── postprocess.py ├── removeCellBoxes.py ├── system ├── x.dat ├── y.dat └── z.dat |
|
|
|
|
|
|
| Tags |
| buildings, simplefoam, urban wind, validation |
| Thread Tools | |
| Display Modes | |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| CFX problem in ubuntu (linux) | Vigneshramaero | CFX | 0 | July 13, 2012 10:22 |
| CFX-Pre problem, pls help!!! | cth_yao | CFX | 0 | February 17, 2012 00:52 |
| Porous jump in urban environment | samygero | FLUENT | 0 | February 14, 2011 18:26 |
1Likes
Linear Mode
