Good morning,
I would like to ask around if anyone has a working tutorial case for the ABLTerrainSolver avaliable and is willing to share it?
I wonder in genereral what has to be done to make if work? I think one has to create a mesh above a terrain and let the
ABLTerrainSolver do the trick? What about the boundary patches? They cann’t still be cyclic can they? How do I change
them and to what?
Thanks a lot!
Best,
Thomas
Hello everyone,
I’m still trying to create a ABLTerrainSolver case and it seems there are some bugs in the implementation
or am I doing something wrong? Somebody here who could help me with this case?
dropbox.com/s/d9zzp19qfasym … e.tgz?dl=0
It uses moveDynamicMesh to attach to the terrain and calls setFieldsABL afterwards. The ABLTerrainSolver
crashes with this error message:
PIMPLE: Operating solver in PISO mode
Starting time loop
<U_1> = (-1.25337e-05 15.0001 0) <U_s> = (0.0136353 14.5618 0.680718) <dU/dn> = (-0.000705281 0.0216922 -0.0335939)
DILUPBiCG: Solving for flm, Initial residual = 1, Final residual = 0.0359155, No Iterations 1
bounding flm, min: -0.0230707 max: 0.0223845 average: -4.22079e-05
DILUPBiCG: Solving for fmm, Initial residual = 1, Final residual = 0.0359094, No Iterations 1
bounding fmm, min: -0.00913156 max: 0.133026 average: 0.00141043
TRef = 300
uStarMean = 1.09808 LMean = 1e+30 phiMMean = 1
UParallelMeanMag = 14.5776 UParallelPMeanMag = 14.7777
RwMagMean = 1.20486 RwMeanMag = 1.20447 sqrt(RwMagMean) = 1.09766 sqrt(RwMeanMag) = 1.09748 uStarMean = 1.09808
Time = 21 Time Step = 21
Courant Number mean: 1.55351 max: 2.0861
deltaT = 0.358696
<U_1> = (-1.25337e-05 15.0001 0) <U_s> = (0.0136353 14.5618 0.680718) <dU/dn> = (-0.000705281 0.0216922 -0.0335939)
DILUPBiCG: Solving for Ux, Initial residual = 0.350015, Final residual = 1.96455e-16, No Iterations 6
DILUPBiCG: Solving for Uy, Initial residual = 0.970456, Final residual = 1.77876e-16, No Iterations 6
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.45511e-17, No Iterations 6
DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 5.24417e-16, No Iterations 6
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00935106, No Iterations 11
<U_1> = (-0.0509501 15.9364 0.476419) <U_s> = (-0.0439888 15.6557 0.746193) <dU/dn> = (-0.000636055 0.0133576 -0.0127496)
DILUPBiCG: Solving for T, Initial residual = 0.242254, Final residual = 2.41253e-17, No Iterations 6
GAMG: Solving for p_rgh, Initial residual = 0.119069, Final residual = 0.00113359, No Iterations 8
<U_1> = (-0.0329056 15.9365 0.484592) <U_s> = (-0.0282349 15.7298 0.752945) <dU/dn> = (-0.000600525 0.0116376 -0.0125505)
DILUPBiCG: Solving for T, Initial residual = 0.237378, Final residual = 2.67629e-17, No Iterations 6
GAMG: Solving for p_rgh, Initial residual = 0.0170006, Final residual = 8.93181e-09, No Iterations 35
<U_1> = (-0.0303681 15.9439 0.485109) <U_s> = (-0.0256299 15.7398 0.752529) <dU/dn> = (-0.000603657 0.0115944 -0.0125146)
DILUPBiCG: Solving for T, Initial residual = 0.234728, Final residual = 9.29193e-17, No Iterations 6
Local Flux Continuity Error: Min 6.62286e-16 Max 1.27786 Weighted Mean 1.99251e-05
Total Boundary Flux: -5361.84
Total Boundary Area: 2.64029e+06
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3
at ??:?
#4 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#5
at ??:?I would loooove to find out why this is happening! I am using the current github release
of SOWFA and OpenFOAM 2.3.0 on a 14.04 XUbuntu 64bit system (to setup the case).
The case does only run up to this point due to the fact that I disabled the “fixedHeatingRate”
in the qwall file (0.orig folder) and set it to “zeroGradient”. There seems to be an error in
the patch field source too. Please help!
Best,
Thomas
The problem has been solved … have a look into the “SOWFA support” thread…