Rayleigh benard cylindrical

[alexandr...@gmail.com]

Hello,

I'm trying to simulate Rayleigh Bénard convection inside a 3D cylinder :

http://www.basilisk.fr/sandbox/Vierron/Natural-convection/cylindrical/rbcyl-emb.c

#include "grid/octree.h"
#include "embed.h"
#include "navier-stokes/centered.h"
#include "tracer.h"
#include "diffusion.h"
#include "radial.h"

Here a snapshot of basilisk view :

I have a problem with my videos (T.mp4 and T2.mp4).

As you can see the area of square or isosurface is too small and I only get half of my domain.

I don't understand why ?

To find the view parameters, I used the basilisk view export and simply added map=radial.

Best regards,

Alexandre

[alexandr...@gmail.com]

Hello again,

I've updated the source code: darcs pull

But now my code doesn't work anymore !! :

/home/alexandre/basilisk/src/qcc -DDISPLAY=-1 -O2 -g -Wall -pipe -D_FORTIFY_SOURCE=2 -MD -o rbcyl2.s.d rbcyl2.c
*updating Makefile.deps
qcc -autolink -DDISPLAY=-1 -O2 -g -Wall -pipe -D_FORTIFY_SOURCE=2 -o rbcyl2/rbcyl2 rbcyl2.c -lm
qcc -autolink -DDISPLAY=-1 -O2 -g -Wall -pipe -D_FORTIFY_SOURCE=2 -o rbcyl2/rbcyl2 rbcyl2.c  -lm
/home/alexandre/basilisk/src/grid/tree-common.h:182: warning: Basilisk C parse error near `v v v v v v v
    scalar * listr = list_concat (p.slist, {cm})'
/home/alexandre/basilisk/src/grid/tree-common.h:184: warning: Basilisk C parse error near `=============
    scalar * listr = list_concat (p.slist, {cm})'
/home/alexandre/basilisk/src/grid/tree-common.h:186: warning: Basilisk C parse error near `*************
    scalar * listr = list_concat (slist, {cm})'
/home/alexandre/basilisk/src/grid/tree-common.h:188: warning: Basilisk C parse error near `^ ^ ^ ^ ^ ^ ^
    restriction (listr)'
/home/alexandre/basilisk/src/grid/tree-common.h: In function ‘adapt_wavelet’:
/home/alexandre/basilisk/src/grid/tree-common.h:181:1: error: unknown type name ‘v’
/home/alexandre/basilisk/src/grid/tree-common.h:181:5: error: expected ‘=’, ‘,’, ‘;’, ‘asm’ or ‘__attribute__’ before ‘v’
/home/alexandre/basilisk/src/grid/tree-common.h:182:48: error: expected statement before ‘)’ token
/home/alexandre/basilisk/src/grid/tree-common.h:183:1: error: expected expression before ‘==’ token
/home/alexandre/basilisk/src/grid/tree-common.h:186:5: error: expected expression before ‘scalar’
/home/alexandre/basilisk/src/grid/tree-common.h:187:1: error: expected expression before ‘^’ token
/home/alexandre/basilisk/src/grid/tree-common.h:189:12: error: ‘listr’ undeclared (first use in this function); did you mean ‘list’?
/home/alexandre/basilisk/src/common.h:365:39: note: in definition of macro ‘pfree’
/home/alexandre/basilisk/src/grid/tree-common.h:189:12: note: each undeclared identifier is reported only once for each function it appears in
/home/alexandre/basilisk/src/common.h:365:39: note: in definition of macro ‘pfree’
make: *** [rbcyl2.s] Error 1
make: *** [/home/alexandre/basilisk/src/Makefile.defs:66 : rbcyl2.s] Erreur 1

Best regards,

Alexandre

[Stephane Popinet]

> /home/alexandre/basilisk/src/grid/tree-common.h:182: warning: Basilisk C
> parse error near `v v v v v v v
>     scalar * listr = list_concat (p.slist, {cm})'
> /home/alexandre/basilisk/src/grid/tree-common.h:184: warning: Basilisk C
> parse error near `=============

You have modified your basilisk code locally and you have a 'darcs pull'
conflict. Darcs has warned you about this, have a look at the message.

The conflicting lines are marked with 'v v v v v v' and '======' etc.

You need to fix the conflict, see:

https://darcsbook.acmelabs.space/chapter06.html#conflict-markers

Stephane

[alexandr...@gmail.com]

Ok thank you Stephane,

But now :

qcc -autolink -DDISPLAY=-1 -O2 -g -Wall -pipe -D_FORTIFY_SOURCE=2 -o rbcyl2/rbcyl2 rbcyl2.c -lm
qcc -autolink -DDISPLAY=-1 -O2 -g -Wall -pipe -D_FORTIFY_SOURCE=2 -o rbcyl2/rbcyl2 rbcyl2.c  -lm

/home/alexandre/basilisk/src/grid/tree-common.h: In function ‘adapt_wavelet’:
/home/alexandre/basilisk/src/grid/tree-common.h:176:27: error: ‘cs’ undeclared (first use in this function); did you mean ‘cm’?
/home/alexandre/basilisk/src/common.h:50:26: note: in definition of macro ‘is_constant’
/home/alexandre/basilisk/src/grid/tree-common.h:176:27: note: each undeclared identifier is reported only once for each function it appears in
/home/alexandre/basilisk/src/common.h:50:26: note: in definition of macro ‘is_constant’
/home/alexandre/basilisk/src/grid/tree-common.h:177:67: error: expected expression before ‘{’ token
/home/alexandre/basilisk/src/grid/tree-mpi.h: In function ‘restore_mpi’:
/home/alexandre/basilisk/src/grid/tree-mpi.h:1305:22: error: ‘cs’ undeclared (first use in this function); did you mean ‘cm’?
/home/alexandre/basilisk/src/grid/tree-mpi.h:1306:66: error: ‘fs’ undeclared (first use in this function); did you mean ‘fp’?
/home/alexandre/basilisk/src/grid/tree-mpi.h:1306:70: warning: excess elements in scalar initializer
/home/alexandre/basilisk/src/grid/tree-mpi.h:1306:70: note: (near initialization for ‘(anonymous)’)

make: *** [rbcyl2.s] Error 1
make: *** [/home/alexandre/basilisk/src/Makefile.defs:66 : rbcyl2.s] Erreur 1

Best regards,

Alexandre

[alexandr...@gmail.com]

Ok,

I had the same problem before : https://groups.google.com/g/basilisk-fr/c/zWm0ygsx6is/m/S5EFvS7ZAAAJ

So : darcs rollback -p "embed default variables before metric ones reorderedin lists for adaption"

And It works again.

But I still have my problem with my videos ...

Thanks,

Alexandre