Search found 86 matches

by Pisolino
Tue Feb 14, 2017 4:53 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Re: Removing mesh volume elements - Evolving (??) mesh

Hi Yvan, thanks ! now it works in single core as it is possible to see in the attached file. Now i'm trying to extend the algorithm to a multi thread process. When i run the same algorithm i found : ==9412== Conditional jump or move depends on uninitialised value(s) ==9412== at 0x401063: bounfunc (c...
by Pisolino
Mon Feb 13, 2017 7:19 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Re: Removing mesh volume elements - Evolving (??) mesh

Hi yvan, thank you, today i have done some experiments with the cs_user_scripts.py. So i wanted to start to define my condition for cell selecting, and again another error. I'm not sure how to use yplus value since it is a boundary value field. Below an extraction of the bounfunc defined in this pos...
by Pisolino
Fri Feb 10, 2017 6:05 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Re: Removing mesh volume elements - Evolving (??) mesh

Hi Yvan, after some debugging the function is working, the boundary condition is transmitted to the i_face. Here the debugged code: static void bounfunc(void *input, cs_lnum_t *n_cells, cs_lnum_t **cell_ids) { cs_lnum_t i, face_id; cs_lnum_t _n_cells = 0; cs_lnum_t *_cell_ids = NULL; cs_lnum_t n_b_f...
by Pisolino
Sun Feb 05, 2017 5:58 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Re: Removing mesh volume elements - Evolving (??) mesh

Hi, thank you again for the big hint. i've started from your function static void bounfunc(void *input, cs_lnum_t *n_cells, cs_lnum_t **cell_ids) { cs_lnum_t i, face_id; cs_lnum_t _n_cells = 0; cs_lnum_t *_cell_ids = NULL; cs_lnum_t n_b_faces = 0; //cs_lnum_t *b_face_ids = NULL; /* import the mesh *...
by Pisolino
Thu Feb 02, 2017 11:45 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Re: Removing mesh volume elements - Evolving (??) mesh

Sorry to bother you again, after a full day in solving this i'm almost losing my mind. I'm continuing to read from fortran and C subroutines and i don't have all the commands to perform something with a language or another one. Despite the logic of what i would do is quite clear to me i don't unders...
by Pisolino
Wed Feb 01, 2017 4:38 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Re: Removing mesh volume elements - Evolving (??) mesh

Hi, again thank you for the support. Just to be sure i understood everything and as a future starting point for other users: in the example cs_user_postprocess.c file : the example function (_t_lt_21_select) select the desired cells with temperature == 21. then a writer is defined, i'll define the m...
by Pisolino
Tue Jan 31, 2017 6:08 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Re: Removing mesh volume elements - Evolving (??) mesh

Thank you for the reply, for my purposes it's ok to restart the simulation from initial values (maybe it's a little slow but it is ok), i'm very happy that it seems possible to perform such simulation. I'm taking a look at cs_user_postprocessing functions. My objective is to automatically: - after N...
by Pisolino
Tue Jan 31, 2017 4:57 pm
Forum: code_saturne usage
Topic: Removing mesh volume elements - Evolving (??) mesh
Replies: 21
Views: 13103

Removing mesh volume elements - Evolving (??) mesh

Hi all, i'm wondering if it is possible, in order to perform shape evolution, in CS to remove some elements (tetra) at the boundary from the mesh and to project the same boundary condition to the new-generated boundary faces ? This would be done every N simulations, then the simulation re-starts aga...
by Pisolino
Sun Nov 20, 2016 11:52 pm
Forum: code_saturne usage
Topic: Topological Optimizazion in Code Saturne
Replies: 12
Views: 13078

Re: Topological Optimizazion in Code Saturne

Hi, every X iterations i increase the head loss in the cells i want to penalise. During this X iterations the flow adapts and then i penalise again. Now i want to add a kind of filter to prevent "chequered" zone.
by Pisolino
Sun Nov 20, 2016 12:25 pm
Forum: code_saturne usage
Topic: Topological Optimizazion in Code Saturne
Replies: 12
Views: 13078

Re: Topological Optimizazion in Code Saturne

Thank you Luciano, for the moment i tried the easy solution : playing with head losses instead of density and increase relaxation. Avoiding storng jumps is quite difficult because of the algorithm nature, i'm using BESO that has only logic conditions ( in solid BESO top-opt it has only density = 1 o...