stardis-spk

Scripts and datasets for running Stardis simulations
git clone git://git.meso-star.fr/stardis-spk.git
Log | Files | Refs | README | LICENSE

commit 0e6ca9dba92f2638b22ce65d12684aeb4bdf3f05
parent f84c3792c58e44379d5e38c632efbd6c651d277e
Author: Christophe Coustet <christophe.coustet@meso-star.com>
Date:   Mon,  4 Jan 2021 17:18:58 +0100

Quelques changements dans le README

Diffstat:
MREADME.md | 42+++++++++++++++++++-----------------------
Mcube/run_green_evaluation.sh | 2+-
2 files changed, 20 insertions(+), 24 deletions(-)

diff --git a/README.md b/README.md @@ -46,8 +46,8 @@ in these boundary STL files which coincide with the triangles in the solid STL files must be *rigorously* the same. This *conformal mesh* constraint is also required for adjacent solids that share a common interface. -Finally you can watch the `model.txt` file which is the `stardis` input data -file. In this file we *connect*: +Finally you can have a look at the `model.txt` file which is the `stardis` input +data file. In this file we *connect*: - physical properties (thermal conductivity, thermal capacity, ...) to the geometrical data (here `solid.stl` only); @@ -79,13 +79,13 @@ You can also simply invoke the `stardis` program in order to compute the temperature at the center of the cube at steady state, by using the following command: - ~/Stardis-StarterPack/cube $ stardis -M model.txt -p 0.5,0.5,0.5,INF -e + ~/Stardis-StarterPack/cube $ stardis -M model.txt -p 0.5,0.5,0.5 -e Please refer to the [stardis](https://www.meso-star.com/projects/stardis/man/man1/stardis.1.html) man page for an explanation about command-line options. The Bash script can be edited and modified. For instance, in section "USER PARAMETERS", the number of -Monte-Carlo samples or the value of the probe time can be updated. +Monte-Carlo samples or the value of the probe time can be changed. #### Dump some "thermal paths" @@ -103,7 +103,7 @@ The script `run_dump_paths.sh` invokes `stardis` twice: yet allows the visualization of thermal paths in a complex geometry where radiative, convective and conductive heat transfers are coupled. -You can update the numbers of paths or the probe position by editing the script. +You can modify the number of paths or the probe position by editing the script. #### Green function evaluation @@ -111,10 +111,8 @@ The last script `run_green_evaluation.sh` shows how to estimate the Green function with `stardis` and how to use it with the `sgreen` program. On first invocation, the script runs `stardis` to evaluate the Green function by generating the required number of thermal paths and store some data (the end -position of each path) in a binary file. This Green function is evaluated only -for a probe position (defined in the `USER PARAMETERS SECTION`: whenever -`X`, `Y`, `Z` or the number of Monte-Carlo samples `NREAL` are -modified, a new Green function will be evaluated. +position of each path) in a binary file. This Green function is evaluated +for the probe position defined in the `USER PARAMETERS SECTION`. This Green function is independent of the value of the sources (values of boundary and initial conditions, as well as the volumetric source term). If you @@ -125,20 +123,18 @@ the current values of sources (no matter whether they have been modified or not); in order to modify the values of the sources, the following line should be modified in the script: - SOURCES_AND_BOUNDARIES="CUBE.VP = 12 LTEMP.T = 290 RTEMP.T = 310 ADIA.F = 5.2" + SOURCES_AND_BOUNDARIES="CUBE.VP=12 LTEMP.T=290 RTEMP.T=310 ADIA.F=5.2" with `CUBE.VP` the value of the volumetric source term (in W/m^3); `LTEMP.T` and `RTEMP.T` the values of the temperature on the left (`LTEMP`) and right -(`RTEMP`) boundaries respectively; and `ADIA.F` the value of the heat flux -density imposed to the boundary `ADIA` (in W/m^2). The names `CUBE`, `RTEMP`, -`LTEMP` and `ADIA` refer to the names given in the file `model.txt` +(`RTEMP`) boundaries respectively (in K); and `ADIA.F` the value of the heat flux +density imposed to the boundary `ADIA` (in W/m^2). The names of the sources (here +`CUBE`, `RTEMP`, `LTEMP` and `ADIA`) are simply the names used in the `model.txt` +file. -We point out that it is not currently possible to define the value of `ADIA.T`: -in the `model.txt` file, the boundary `ADIA` has been assigned a value of flux -density using the `F_BOUNDARY_FOR_SOLID` keyword. If this boundary should be -assigned a given temperature (Dirichlet condition), a different `model.txt` -file should be generated using the `T_BOUNDARY_FOR_SOLID` keyword for this -boundary. +Please note that `ADIA` being a boundary of type flux (defined by the +`F_BOUNDARY_FOR_SOLID` keyword in the `model.txt` file) the corresponding source +value is named `ADIA.F` and not `ADIA.T` as for boundaries of type temperature. ### The heatsink @@ -171,9 +167,9 @@ the resolution of the image and the number of samples per pixel. For each pixel of the image, the luminance is computed by Monte-Carlo and the number of realizations is the specified number of samples per pixel. Computing a high-definition image with little statistical noise can therefore take a long -time (many hours). The values of the parameters that are provided in the script -should result in a computational time of about a dozen minutes on a correct -desktop computer. +time (many hours). The values of the parameters provided in the script +should result in a computation time of about a dozen minutes on a recent +low-end desktop computer. More information about the rendering is provided in the `stardis` man page (such as the parameters associated with the point of view). @@ -181,7 +177,7 @@ as the parameters associated with the point of view). Acknowledgement to Cyril Caliot who designed the porous model for the Optisol project (ANR-11-SEED-0009, PROMES-CNRS, CIRIMAT, SICAT, LTN). This model represents an ideal metallic or SiC foam. This type of foam is used in the -design of heat exchangers in concentrated solar processes, in order to transfer +design of heat exchangers in concentrated solar processes in order to transfer incoming solar radiation energy to a working fluid. ## Copyright notice diff --git a/cube/run_green_evaluation.sh b/cube/run_green_evaluation.sh @@ -20,7 +20,7 @@ NREAL=10000 X=0.5 Y=0.5 Z=0.5 -SOURCES_AND_BOUNDARIES="CUBE.VP = 12 LTEMP.T = 290 RTEMP.T = 310 ADIA.F = 5.2" +SOURCES_AND_BOUNDARIES="CUBE.VP=12 LTEMP.T=290 RTEMP.T=310 ADIA.F=5.2" GREEN_FILE=probe_X${X}-Y${Y}-Z${Z}_N${NREAL}.green SETTINGS_FILE=settings.txt ### END USER PARAMETERS SECTION