CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Journal

Barestrand et al. - 2023 - Modeling Convective Heat Transfer of Air in a Data...

Register Blogs Community New Posts Updated Threads Search

Like Tree2Likes
  • 2 Post By milocheung

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 7, 2023, 03:18
Default Barestrand et al. - 2023 - Modeling Convective Heat Transfer of Air in a Data...
  #1
New Member
 
Miguel Nóbrega
Join Date: Jul 2012
Posts: 21
Rep Power: 14
mnobrega is on a distinguished road
Comments and questions are welcomed on the article Barestrand et al., Modeling Convective Heat Transfer of Air in a Data Center Using OpenFOAM: Evaluation of the Boussinesq Buoyancy, OpenFOAM Journal, https://doi.org/10.51560/ofj.v3.59.
mnobrega is offline   Reply With Quote

Old   August 7, 2023, 15:30
Default
  #2
Senior Member
 
Charles
Join Date: Aug 2016
Location: Vancouver, Canada
Posts: 148
Rep Power: 9
Marpole is on a distinguished road
Thanks for the paper. Can you elaborate "An opening is therefore made in the geometry to act as a pressure reference to alleviate this, enabling the use of a specified flow rate for the CRAH-return areas."? Is this opening the supply area (grey)?
__________________
Charles L.
Marpole is offline   Reply With Quote

Old   August 15, 2023, 03:17
Default Pressure Opening
  #3
New Member
 
Henrik
Join Date: Jan 2016
Posts: 5
Rep Power: 10
broccolibadger is on a distinguished road
Hi Marpole,

The pressure opening acts as a sort of reference, indirectly modelling that the data center is not airtight towards the atmosphere. Here it is done by making a small 10x10 cm hole with a simple pressure boundary condition. The effect on the stiffness of the CFD solution is substantial in how it improves convergence due to the other boundaries using enforced mass flows to preserve continuity across racks.

/Henrik
broccolibadger is offline   Reply With Quote

Old   January 8, 2024, 05:01
Default Floating point exception (core dumped) error in steady_incompressible case
  #4
New Member
 
milo cheung
Join Date: Oct 2023
Posts: 5
Rep Power: 2
milocheung is on a distinguished road
Hi everyone, I am really appreciate this paper, and also provided the case files. Unfortunately, I've encountered some issues while attempting to run the "steady_incompressible" case, which has brought my progress to a halt. The detail error file also attached. It will be nice if anyone can help. Thank you very much.

Here's my setting and changes.
use finest grid no changes,"40mm, UUF" that one
system/meshParDict file. Updated "numberOfSubdomains" to 16
create logs folder
Copy and run the commands from "Allrun" file.
nmesh=16
rm -r constant/polyMesh/
blockMesh | tee logs/log1_blockMesh
surfaceFeatureExtract | tee logs/log1_features
rm -r processor*
decomposePar -decomposeParDict system/meshParDict | tee logs/log1_decompose
mpirun --oversubscribe -np $nmesh renumberMesh -decomposeParDict system/meshParDict -overwrite -parallel | tee logs/log1_rnM
mpirun --oversubscribe -np $nmesh snappyHexMesh -decomposeParDict system/meshParDict -overwrite -parallel | tee logs/log1_snappy
reconstructParMesh -constant -fullMatch | tee logs/log1_reconstruct
rm -r processor*
checkMesh | tee logs/log1_check
rm -r 0
cp -r 0.org 0
rm -r postprocessing
renumberMesh -overwrite | tee logs/log2e_rnM
buoyantBoussinesqSimpleFoam

Error when I run "renumberMesh".No problem when I pass "--oversubscribe" option in every mpirun.
command:
mpirun --oversubscribe -np $nmesh renumberMesh -decomposeParDict system/meshParDict -overwrite -parallel | tee logs/log1_rnM

error message:
There are not enough slots available in the system to satisfy the 16
slots that were requested by the application:

renumberMesh
Error when I run "buoyantBoussinesqSimpleFoam"
command:
buoyantBoussinesqSimpleFoam

error:
Time = 9

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PBiCGStab::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 ? at ??:?
#7 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<do uble> >&, Foam::dictionary const&) const at ??:?
#8 ? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/buoyantBoussinesqSimpleFoam
#9 ? in /lib/x86_64-linux-gnu/libc.so.6
#10 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#11 ? in /usr/lib/openfoam/openfoam2012/platforms/linux64GccDPInt32Opt/bin/buoyantBoussinesqSimpleFoam
Floating point exception (core dumped)

__________________
Milo
Attached Files
File Type: txt errorDetail.txt (23.7 KB, 3 views)

Last edited by milocheung; January 8, 2024 at 20:48.
milocheung is offline   Reply With Quote

Old   January 8, 2024, 14:41
Default
  #5
Senior Member
 
Charles
Join Date: Aug 2016
Location: Vancouver, Canada
Posts: 148
Rep Power: 9
Marpole is on a distinguished road
As you can read from the error message "error message: There are not enough slots available in the system to satisfy the 16", the program asks for using 16 CPU and your computer might have less than 16 CPU.
__________________
Charles L.
Marpole is offline   Reply With Quote

Old   January 8, 2024, 21:20
Default
  #6
New Member
 
milo cheung
Join Date: Oct 2023
Posts: 5
Rep Power: 2
milocheung is on a distinguished road
Hi Marpole,

Thank you for your reply. Apologies for the confusion due to the unclear typesetting; I have two distinct issues to address.

For the first issue regarding MPI, I believe I've found a potential solution. I'm considering using the "oversubscribe" option with the mpirun command, which seems appropriate given my computer's specs: an AMD Ryzen™ 9 7940HS Processor featuring 8 cores and 16 threads. This should allow me to utilize all available threads even if the number of processes exceeds the number of physical cores.

Regarding the second problem, I encountered a "Floating point exception" error while running buoyantBoussinesqSimpleFoam. Any idea for it?


__________________
Milo
milocheung is offline   Reply With Quote

Old   January 8, 2024, 22:27
Default
  #7
Senior Member
 
Charles
Join Date: Aug 2016
Location: Vancouver, Canada
Posts: 148
Rep Power: 9
Marpole is on a distinguished road
Hi Milo,

The second error is obviously a divergence error as you can see the residual of epsilon, k and the error of continuity. When reading the max/min U from the log, they are way large (e.g. 295,721m/s at time step 2).

This can be caused by a number of reasons. But I would review mesh and mesh quality first (file log1_check).
__________________
Charles L.
Marpole is offline   Reply With Quote

Old   January 8, 2024, 22:43
Default
  #8
New Member
 
milo cheung
Join Date: Oct 2023
Posts: 5
Rep Power: 2
milocheung is on a distinguished road
Hi Marpole,

File uploaded. Thank you.
Attached Files
File Type: txt log1_check.txt (6.2 KB, 3 views)
milocheung is offline   Reply With Quote

Old   January 8, 2024, 23:15
Default
  #9
Senior Member
 
Charles
Join Date: Aug 2016
Location: Vancouver, Canada
Posts: 148
Rep Power: 9
Marpole is on a distinguished road
The mesh quality is Okay. Can you display your decomposeParDict?
__________________
Charles L.
Marpole is offline   Reply With Quote

Old   January 8, 2024, 23:35
Default
  #10
New Member
 
milo cheung
Join Date: Oct 2023
Posts: 5
Rep Power: 2
milocheung is on a distinguished road
Hi Marpole,

I pass "system/meshParDict" in "-decomposeParDict" option in following commands:
decomposePar -decomposeParDict system/meshParDict | tee logs/log1_decompose
mpirun -np $nmesh renumberMesh -decomposeParDict system/meshParDict -overwrite -parallel | tee logs/log1_rnM
mpirun -np $nmesh snappyHexMesh -decomposeParDict system/meshParDict -overwrite -parallel | tee logs/log1_snappy
system/meshParDict content

Code:
FoamFile
{
 class dictionary;
 format ascii;
 object decomposeParDict;
 version 2.0;
}

numberOfSubdomains  16;

method          scotch;
system/decomposeParDict

Code:
FoamFile
{
 class dictionary;
 format ascii;
 object decomposeParDict;
 version 2.0;
}

numberOfSubdomains  88;

method          metis;

coeffs
{
    n           (4 4 6);
}


The log1_decompose, log1_rnM, log1_snappy file also uploaded. Thank you.


__________________
Milo
Attached Files
File Type: txt log1_snappy.txt (49.0 KB, 1 views)
File Type: txt log1_rnM.txt (2.5 KB, 1 views)
File Type: txt log1_decompose.txt (18.2 KB, 2 views)
milocheung is offline   Reply With Quote

Old   January 14, 2024, 23:08
Default
  #11
New Member
 
milo cheung
Join Date: Oct 2023
Posts: 5
Rep Power: 2
milocheung is on a distinguished road
Hi,

I can execute the case files now. I downgrade the swak4Foam to v2021.05, and ubuntu from 22.04 to 20.04. Then I am able to run the cases. I think the ubuntu version is the main reason.

Thank you.
__________________
Milo
broccolibadger and Marpole like this.
milocheung is offline   Reply With Quote

Reply

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[General] Extracting ParaView Data into Python Arrays Jeffzda ParaView 30 November 6, 2023 21:00
Moist air condensation modeling for heat exchanger rdp5252 Fluent Multiphase 0 November 18, 2021 14:43
Heat transfer between wall and air in Fluent jimmy871013 FLUENT 6 April 27, 2018 02:00
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
No heat transfer in air thomasyangfly CFX 1 April 18, 2013 06:35


All times are GMT -4. The time now is 02:01.