two-phase jet+embed.h+axi.h

[Xuerao He]

Dear Basilisk users,

I am trying to do this case: an axisymmetric two-phase jet with an embedded needle, as shown in the figure sketch attached. It is modified from the example http://basilisk.fr/src/examples/atomisation.c. The code can run successfully when the two-phase jet flow adds axi.h  or adds a solid needle only, but it comes across errors when it adds both axi.h and the embed.h. I have attached the code file and the error screenshot.  

I want to ask if axi.h embed.h and two-phase.h can be used together? If they can work together, could someone please point out where is wrong in my code?

Best regards,

He Xuerao

[jmlopez...@gmail.com]

Hi,

I have been taking a look to your code. I find out some flaws in the code and a bug on my side (see jet_embed.c below):
1.- I do not know why solid() seems to rotate the solid in axi or get confuse with which is the radial coordinate. It is a little bit strange. It cam be easily solved using vertex field
2.- If you create/modified a solid (cs & fm) you have to update cm & fm

3.-  The upgrade of fm (and cm?) seems to be wrong. The bug appears in the adaption. I think that maybe it is a problem with the list of variables and/or the order they are updated. Comenting "adapt" the code seems to work.

I will try to fix it.

Cheers

[Xuerao He]

Dear Prof. J.M. Lopez-Herrera,

Thank you very much for your response and for offering your assistance. I am a beginner of Basilisk. In fact, I am interested in reproducing your work on electrospray in 2012 which is calculated with GERRIS(Numerical simulation of electrospray in the cone-jet mode). I would like to confirm if it is possible to achieve the same results using the current version of Basilisk.

Regarding the problem of code error, I saw an email you sent in the forum before https://groups.google.com/g/basilisk-fr/c/tBY2Q3cUN-Q/m/EqUaeJWWCwAJ, suggesting using mask matrix to realize solid. I tried to do this, and it seems to be feasible. Please take a look at the code in my attachment. 

Thanks again.

Best regards，

Xuerao

jmlopez...@gmail.com <jmlopez...@gmail.com> 于2023年6月7日周三 22:39写道：

--
You received this message because you are subscribed to a topic in the Google Groups "basilisk-fr" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/basilisk-fr/HuBi2hRkMwE/unsubscribe.
To unsubscribe from this group and all its topics, send an email to basilisk-fr...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/basilisk-fr/d679a0af-920e-40be-a632-bd389d394b34n%40googlegroups.com.