fluid interface detection with petsc and donlp2
play

Fluid Interface Detection with PETSc and DONLP2 PETSc User Meeting - PowerPoint PPT Presentation

Fluid Interface Detection with PETSc and DONLP2 PETSc User Meeting Vienna 2016 Poster Session Alin-Adrian Anton and *Sebastian Muntean Politehnica University of Timi soara, Romania *Romanian Academy Timi soara Branch


  1. Fluid Interface Detection with PETSc and DONLP2 PETSc User Meeting Vienna 2016 – Poster Session Alin-Adrian Anton and *Sebastian Muntean Politehnica University of Timi¸ soara, Romania *Romanian Academy Timi¸ soara Branch alin.anton@cs.upt.ro June 28, 2016 Eng. Alin-Adrian ANTON, PhD PETSc User Meeting Vienna 2016 1/6

  2. Vortex rope (self-induced instability) taking place in swirling flows the stalled region is filled with stagnant water 3D unsteady flow is modeled considering a 2D axisymmetric steady flow the flowing-stalled fluid interface can be determined using interface capturing techniques (ICaT) and interface tracking techniques (ITrT) Eng. Alin-Adrian ANTON, PhD PETSc User Meeting Vienna 2016 2/6

  3. Two-dimensional axisymmetric flow numerical simulation with stagnant region computed with ICaT Eng. Alin-Adrian ANTON, PhD PETSc User Meeting Vienna 2016 3/6

  4. SWIRL2D 2D boundaries ( x,r ) + PROBLEM (1) interface boundary ( r i , i=1...n ) + DATA boundary conditions initial guess (2) PROBLEM inter ace boundary (r i , i=1...n ) f PARAMETERS Computational domain (3) generation (1D + 2D) PETSc interface boundary (4) Mesh generation TRIANGLE new (1D + 2D) interface boundary (5) 2D flow ( r i , i=1...n ) PETSc computation (8) (6) GNU GSL Objective F, F* function (7) min F NO DONLP2 max F* SQP m ethod YES (10) Optimum Optimum (9) interface boundary 2D flow solution TECPLOT ( r i , i=1...n ) Eng. Alin-Adrian ANTON, PhD PETSc User Meeting Vienna 2016 4/6

  5. the moving knots are used for SQP interface optimization the interface is obtained through interpolation using cubic splines a new mesh is generated using TRIANGLE and the 2D section axisymmetric flow is solved using the finite element method (FEM) Eng. Alin-Adrian ANTON, PhD PETSc User Meeting Vienna 2016 5/6

  6. SWIRL2D solution (upper meridian half-plane) and FLUENT2D axisymmetric inviscid solution (lower meridian half-plane) Eng. Alin-Adrian ANTON, PhD PETSc User Meeting Vienna 2016 6/6

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend


More recommend