Vof

src/Make.incflo

@@ -14,6 +14,7 @@ Bdirs += src/projection
 Bdirs += src/rheology
 Bdirs += src/setup
 Bdirs += src/utilities
+Bdirs += src/vof
 
 ifeq ($(USE_PARTICLES), TRUE)
   DEFINES += -DINCFLO_USE_PARTICLES

src/incflo.H

@@ -20,6 +20,7 @@
 
 #include <DiffusionTensorOp.H>
 #include <DiffusionScalarOp.H>
+#include <VolumeOfFluid.H>
 
 enum struct StepType {
     Predictor, Corrector
@@ -31,6 +32,7 @@ public:
 
     friend DiffusionTensorOp;
     friend DiffusionScalarOp;
+    friend VolumeOfFluid;
 
     enum struct FluidModel {
         Newtonian, powerlaw, Bingham, HerschelBulkley, deSouzaMendesDutra
@@ -211,6 +213,18 @@ public:
                            amrex::Vector<amrex::MultiFab const*> const& eta,
                            amrex::Real dt_diff);
 
+    ///////////////////////////////////////////////////////////////////////////
+    //
+    //  tacer advection by VOF method
+    //
+    ////////////////////////////////////////////////////////////////////////////
+
+    void tracer_vof_advection (amrex::Vector<amrex::MultiFab*> const& tracer,
+                           AMREX_D_DECL(amrex::Vector<amrex::MultiFab const*> const& u_mac,
+                                        amrex::Vector<amrex::MultiFab const*> const& v_mac,
+                                        amrex::Vector<amrex::MultiFab const*> const& w_mac));
+
+
     [[nodiscard]] amrex::Array<amrex::MultiFab,AMREX_SPACEDIM>
     average_scalar_eta_to_faces (int lev, int comp, amrex::MultiFab const& cc_eta) const;
 
@@ -675,6 +689,8 @@ private:
     std::unique_ptr<DiffusionTensorOp> m_diffusion_tensor_op;
     std::unique_ptr<DiffusionScalarOp> m_diffusion_scalar_op;
 
+    //vof class pointer
+    std::unique_ptr<VolumeOfFluid> p_volume_of_fluid;
     //
     // end of member variables
     //
@@ -755,6 +771,8 @@ private:
     DiffusionTensorOp* get_diffusion_tensor_op ();
     DiffusionScalarOp* get_diffusion_scalar_op ();
 
+    VolumeOfFluid* get_volume_of_fluid ();
+
     amrex::Vector<amrex::MultiFab*> get_velocity_old () noexcept;
     amrex::Vector<amrex::MultiFab*> get_velocity_new () noexcept;
     amrex::Vector<amrex::MultiFab*> get_velocity_eb () noexcept;

src/incflo.cpp

@@ -164,6 +164,8 @@ void incflo::Evolve()
         if (writeNow())
         {
             WritePlotFile();
+            get_volume_of_fluid()->WriteTecPlotFile (m_cur_time,m_nstep);
+        get_volume_of_fluid()->write_tecplot_surface(m_cur_time,m_nstep);
             m_last_plt = m_nstep;
         }
 

src/incflo_apply_predictor.cpp

@@ -215,4 +215,10 @@ void incflo::ApplyPredictor (bool incremental_projection)
                                    AMREX_D_DECL(GetVecOfConstPtrs(u_mac), GetVecOfConstPtrs(v_mac),
                                    GetVecOfConstPtrs(w_mac)));
 #endif
+
+// use vof to advect tracer
+
+    tracer_vof_advection(get_tracer_new (), AMREX_D_DECL(GetVecOfConstPtrs(u_mac), GetVecOfConstPtrs(v_mac),
+                                         GetVecOfConstPtrs(w_mac)));
+
 }

src/incflo_regrid.cpp

@@ -42,6 +42,8 @@ void incflo::MakeNewLevelFromCoarse (int lev,
     m_diffusion_tensor_op.reset();
     m_diffusion_scalar_op.reset();
 
+    p_volume_of_fluid.reset();
+
     // Note: finest_level has not yet been updated and so we use lev
 #ifdef AMREX_USE_EB
     macproj = std::make_unique<Hydro::MacProjector>(Geom(0,lev),
@@ -92,6 +94,8 @@ void incflo::RemakeLevel (int lev, Real time, const BoxArray& ba,
     m_diffusion_tensor_op.reset();
     m_diffusion_scalar_op.reset();
 
+    p_volume_of_fluid.reset();
+
 #ifdef AMREX_USE_EB
     macproj = std::make_unique<Hydro::MacProjector>(Geom(0,finest_level),
                       MLMG::Location::FaceCentroid,  // Location of mac_vec
@@ -115,5 +119,6 @@ void incflo::ClearLevel (int lev)
     m_factory[lev].reset();
     m_diffusion_tensor_op.reset();
     m_diffusion_scalar_op.reset();
+    p_volume_of_fluid.reset();
     macproj.reset();
 }

src/prob/incflo_prob_I.H

@@ -164,4 +164,13 @@
                     amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> const& dx,
                     amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> const& problo,
                     amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> const& probhi) const;
+
+    void droplet (amrex::Box const& vbx, amrex::Box const& gbx,
+                               amrex::Array4<amrex::Real> const& vel,
+                               amrex::Array4<amrex::Real> const& density,
+                               amrex::Array4<amrex::Real> const& tracer,
+                               amrex::Box const& domain,
+                               amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> const& dx,
+                               amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> const& problo,
+                               amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> const& probhi);
 #endif

src/prob/prob_init_fluid.cpp

@@ -184,12 +184,19 @@ void incflo::prob_init_fluid (int lev)
                                ld.tracer.array(mfi),
                                domain, dx, problo, probhi);
 
+        }
+        else if (1109 == m_probtype)
+        {
         }
         else
         {
             amrex::Abort("prob_init_fluid: unknown m_probtype");
         };
     }
+
+    if (1109 == m_probtype) {
+        tracer_vof_init_fraction(lev, ld.tracer, this);
+    }
 }
 
 void incflo::init_rotating_flow (Box const& vbx, Box const& /*gbx*/,

src/utilities/io.cpp

@@ -343,8 +343,6 @@ void incflo::WritePlotFile()
 
     const std::string& plotfilename = amrex::Concatenate(m_plot_file, m_nstep);
 
-    amrex::Print() << "  Writing plotfile " << plotfilename << " at time " << m_cur_time << std::endl;
-
     int ncomp = 0;
 
     // Velocity components
@@ -693,3 +691,5 @@ void incflo::WritePlotFile()
     particleData.Checkpoint(plotfilename);
 #endif
 }
+
+

src/vof/CMakeLists.txt

@@ -0,0 +1,8 @@
+target_include_directories(incflo PRIVATE ${CMAKE_CURRENT_LIST_DIR})
+
+target_sources(incflo
+   PRIVATE
+   VolumeOfFluid.H
+   VolumeOfFluid.cpp
+   incflo_vof.cpp
+   )

src/vof/Make.package

@@ -0,0 +1,4 @@
+CEXE_sources += VolumeOfFluid.cpp
+CEXE_sources += incflo_vof.cpp
+CEXE_headers += VolumeOfFluid.H
+

src/vof/VolumeOfFluid.H

@@ -0,0 +1,36 @@
+//struct TracerVOF {
+//  GfsVariableTracer parent;
+//
+//  GfsVariable * m[FTT_DIMENSION], * alpha;
+//};
+// Header file (VolumeOfFluid class)
+#ifndef INCFLO_VOF_
+#define INCFLO_VOF_
+class incflo;
+
+class VolumeOfFluid
+{
+public:
+    VolumeOfFluid (incflo* a_incflo);
+    void tracer_vof_advection(amrex::Vector<amrex::MultiFab*> const& tracer,
+                              AMREX_D_DECL(amrex::Vector<amrex::MultiFab const*> const& u_mac,
+                                           amrex::Vector<amrex::MultiFab const*> const& v_mac,
+                                           amrex::Vector<amrex::MultiFab const*> const& w_mac),
+                       amrex::Real dt);
+    void write_tecplot_surface(amrex::Real time, int nstep);
+    void WriteTecPlotFile     (amrex::Real time, int nstep);
+
+// normal vector of interface
+   amrex::Vector<amrex::MultiFab> normal;
+// the plane defined by @m.@x = @alpha
+   amrex::Vector<amrex::MultiFab> alpha;
+private:
+incflo* v_incflo; //incflo object
+int finest_level;
+void tracer_vof_update(amrex::Vector<amrex::MultiFab*> const& tracer);
+
+};
+
+void tracer_vof_init_fraction(int lev, amrex::MultiFab& a_tracer, incflo const* a_incflo);
+
+#endif