Skip to content

resolves 160: update 2dbuilding ns testcase #166

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 3 commits into
base: develop
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion examples/modules/heatfluid/examples/2Dbuilding_NS/aero.cfg
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@ case.dimension=2
[heat-fluid]
mesh.filename=$cfgdir/aero.geo
gmsh.hsize=0.025#0.01#0.02#0.07#0.1
filename=$cfgdir/aero.json
json.filename=$cfgdir/aero.json

use-natural-convection=1
Boussinesq.ref-temperature=280#293.15
Expand Down
48 changes: 32 additions & 16 deletions examples/modules/heatfluid/examples/2Dbuilding_NS/aero.json
Original file line number Diff line number Diff line change
Expand Up @@ -2,14 +2,31 @@
{
"Name": "Thermo dynamics",
"ShortName":"ThermoDyn",
"Metadata":
{
"version": "1",
"author": ["vincentchabannes","prudhomm","lberti","jobinhio"],
"date": "28/11/2022",
"references":[""]
},
"Models":
{
"use-model-name":1,
"fluid":
{
"equations": "Navier-Stokes"
}
},
"Meshes":
{
"heatfluid":
{
"Import":
{
"mesh.filename":"$cfgdir/aero.geo",
"gmsh.hsize":"0.025" // 0.025#0.01#0.02#0.07#0.1
}
}
},
"Materials":
{
"air":
Expand Down Expand Up @@ -39,19 +56,19 @@
},
"BoundaryConditions":
{
"velocity":
"fluid":
{
"Dirichlet":
"velocity":
{
"exterior-walls": { "expr":"{0,0}" },
"internal-surfaces": { "expr":"{0,0}" },
"internal-surfaces": { "expr":"{0,0}" },
"heater1": { "expr":"{0,0}" },
"heater2": { "expr":"{0,0}" }
}
},
"temperature":
"heat":
{
"Dirichlet":
"temperature":
{
"heater1": { "expr":"310"/*"330"*/ },
"heater2": { "expr":"300"/*"320"*/ }
Expand All @@ -63,7 +80,7 @@
"expr1":"1.0/(0.06+0.01/0.5 + 0.3/0.8 + 0.20/0.032 +0.016/0.313 +0.14)",// h coeff
"expr2":"280"// temperature exterior
},
"exterior-walls-nofluid":
"exterior-walls-nofluid":
{
"expr1":"1.0/(0.06+0.01/0.5 + 0.3/0.8 + 0.20/0.032 +0.016/0.313 +0.14)",// h coeff
"expr2":"280"// temperature exterior
Expand All @@ -74,19 +91,18 @@
},
"PostProcess":
{
"use-model-name":1,
"heat-fluid":
"heatfluid":
{
"Exports":
{
"fields":["fluid.velocity","fluid.pressure","heat.temperature","fluid.pid"]
}
},
"fluid":
{
},
"heat":
{
},
"fluid":
{
},
"heat":
{
}
}
}
}
12 changes: 7 additions & 5 deletions examples/modules/heatfluid/pages/2Dbuilding_NS/README.adoc
Original file line number Diff line number Diff line change
Expand Up @@ -2,6 +2,8 @@
:page-vtkjs: true
:uri-data: https://github.com/feelpp/toolbox/blob/master/examples/modules/heatfluid/examples
:uri-data-edit: https://github.com/feelpp/toolbox/edit/master/examples/modules/heatfluid/examples
:imagesprefix:
ifdef::env-github,env-browser,env-vscode[:imagesprefix: ../../assets/images/]
:page-tags: case
:page-illustration: 2Dbuilding_NS/velocity.png
:description: We simulate the phenomena of natural convection in an apartment of 2 rooms.
Expand Down Expand Up @@ -33,12 +35,12 @@ stem:[\quad]The 2D thermal model of two parts separated by plaster is considered


.Geometry of stem:[\Omega = \cup^3_{i=1} \Omega_i] with stem:[\partial \Omega = \Gamma_{D} \cup \Gamma_R]

image::2Dbuilding_NS/geo.png[]





=== Boundary condition:

==== Condition of velocity
Expand Down Expand Up @@ -131,7 +133,7 @@ T = T_0

.Mesh with GMSH
|====
a| image:2Dbuilding_NS/mesh.png[]
a| image:{imagesprefix}2Dbuilding_NS/mesh.png[]
|====


Expand Down Expand Up @@ -184,20 +186,20 @@ We observe that the flow circulate from the bottom to the top. The cooler fluid

.Heat fluid over times
|====
a| image:2Dbuilding_NS/image1.png[] a| image:2Dbuilding_NS/image2.png[]
a| image:2Dbuilding_NS/image3.png[] a| image:2Dbuilding_NS/image4.png[]
a| image:{imagesprefix}2Dbuilding_NS/image1.png[] a| image:{imagesprefix}2Dbuilding_NS/image2.png[]
a| image:{imagesprefix}2Dbuilding_NS/image3.png[] a| image:{imagesprefix}2Dbuilding_NS/image4.png[]
|====

.Average of temperature over times
image:2Dbuilding_NS/graph1.png[]


The graph show the steady state of temperature at 300K in all domains.


.Fluid velocity and pressure
image:2Dbuilding_NS/velocity.png[] a|image:2Dbuilding_NS/pressure.png[]


The following video shows the transient state of temperature in the domain.

.Video of the temperature field in the cabin
Expand Down
Loading