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

▶
fluid interface detection with petsc and donlp2
SMART_READER_LITE
LIVE PREVIEW

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

Apr 09, 2023 122 likes •193 views

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

slide-1
SLIDE 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

slide-2
SLIDE 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

slide-3
SLIDE 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

slide-4
SLIDE 4

initial guess inter ace boundary(ri, i=1...n) Computational domain generation (1D + 2D) 2D flow computation min F max F* Optimum interface boundary (ri, i=1...n)

PROBLEM PARAMETERS

Optimum 2D flow solution 2D boundaries (x,r) + interface boundary (ri, i=1...n) + boundary conditions

PROBLEM DATA

Objective function

NO YES

SQP m ethod

DONLP2 TRIANGLE F, F* new interface boundary (ri, i=1...n) (1) (2) (3) (4) (5) (6) (8) (9) (10) interface boundary Mesh generation

(1D + 2D)

PETSc PETSc GNU GSL (7) TECPLOT

SWIRL2D

f

  • Eng. Alin-Adrian ANTON, PhD

PETSc User Meeting Vienna 2016 4/6

slide-5
SLIDE 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

slide-6
SLIDE 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