Bursting bubble with symmetry conditions

[Antonella Succar]

Dear basilisk users,

I am trying to simulate a bursting bubble at the interface of a cylindrical water body. The initial geometry consists of a spherical bubble intersecting with the interface of a cylinder. To try and alleviate the ressources requirements, I split the domain in four parts and simulated a quarter of it with the default symmetry conditions on the boundaries. The cavity collapse seems normal, but the evolution of the jet that forms is not: the liquid breakup is really slow and the jet and droplets are stretched a lot before breakup, as you can see in the attached figures (the time is scaled with the capillary time and interface in white is that at a lower level of refinement).  It seems to me that this can stem from the curvature computation at the boundary. Is that the case? Am I missing something?

Regards,

Antonella

[Anil Kumar]

Hi,

I was also facing the same issue. I was able to resolve it by putting more refinement near the axis of symmetry (near jet). But I don't know the reason why it happened. 

Regards,

Anil Kumar,

PhD Scholar,

IIT Bombay, India

[Shyam Sunder Yadav]

Dear Antonella and Anil Kumar

Please check how you are calculating the specific volume alpha and viscosity

The header two-phase-generic.h has the following

    alphav.x[] = fm.x[]/rho(ff);   and   # define rho(f) (clamp(f,0.,1.)*(rho1 - rho2) + rho2)

One can also have the following

alphav.x[] = fm.x[]*(f/rhol + (1-f)/rhog)

Plot the alpha with the above two formulations, you can notice the difference...

Similarly, viscosity can be based on simple averaging or based on harmonic mean...

These small changes will give different results...

Hope this helps

--
You received this message because you are subscribed to the Google Groups "basilisk-fr" group.
To unsubscribe from this group and stop receiving emails from it, send an email to basilisk-fr...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/basilisk-fr/31672757-a036-4a56-8de9-5479a0f5e4b1n%40googlegroups.com.

--

Dr. Shyam Sunder Yadav
Associate Professor
Mechanical Engineering
BITS Pilani
09902346342
http://www.bits-pilani.ac.in/pilani/ssyadav/Profile

The information contained in this electronic communication is intended solely for the individual(s) or entity to which it is addressed. It may contain proprietary, confidential and/or legally privileged information. Any review, retransmission, dissemination, printing, copying or other use of, or taking any action in reliance on the contents of this information by person(s) or entities other than the intended recipient is strictly prohibited and may be unlawful. If you have received this communication in error, please notify us by responding to this email or telephone and immediately and permanently delete all copies of this message and any attachments from your system(s). The contents of this message do not necessarily represent the views or policies of BITS Pilani.

[Antonella Succar]

Thank you both for your answers. For the first solution, refining a lot near the jet sort of removes the advantage of adaptive refinement. As for the second, I am using the same formulation is for the full 3D domain and for the quarter domain with symmetry conditions. It really seems to me as if it is stemming from the curvature calculation on the boundary? 

Thank you again for both for taking the time to reply!

Antonella

[yogesh bhaskar]

hello Sir , this is yogesh , I am new to basilisk and trying to run the burstingBubble.c case but I am getting segmentation fault (core dump) , can u suggest anything for it .

best regard

yogesh

M.tech

IIT kanpur ,india

[Karnia Werry]

Hello Sir，I am trying to simulate a droplet under Couette flow, and encoutered the same question as you. But I used the whole domain and the periodic boundary. I wonder that did you find where the question is ?Or did you solve the problem?