implementation of the height function

IFDL-WSU · Jul 26, 2024 · 915eb86 · 915eb86
1 parent 6b98a14
commit 915eb86
Show file tree

Hide file tree

Showing 5 changed files with 474 additions and 96 deletions.
diff --git a/src/incflo.cpp b/src/incflo.cpp
@@ -156,7 +156,7 @@ void incflo::Evolve()
             }
         }
         get_volume_of_fluid()->output_droplet(m_cur_time,m_nstep);
-        get_volume_of_fluid()->apply_velocity_field(m_cur_time,m_nstep);
+    //    get_volume_of_fluid()->apply_velocity_field(m_cur_time,m_nstep);
         if (writeNow()){
             get_volume_of_fluid()->WriteTecPlotFile (m_cur_time,m_nstep);
             get_volume_of_fluid()->write_tecplot_surface(m_cur_time,m_nstep);

diff --git a/src/prob/prob_init_fluid.cpp b/src/prob/prob_init_fluid.cpp
@@ -1149,10 +1149,10 @@ void incflo::init_droplet (Box const& vbx, Box const& /*gbx*/,
         Real y = Real(j+0.5)*dx[1];
         Real z = Real(k+0.5)*dx[2];
         Real pi = 3.14159265357;
-        vel(i,j,k,0) = 2*sin(2.*pi*y)*sin(pi*x)*sin(pi*x)*sin(2*pi*z)*cos(pi*0./3.);
-        vel(i,j,k,1) = -sin(2.*pi*x)*sin(pi*y)*sin(pi*y)*sin(2*pi*z)*cos(pi*0./3.);
+        vel(i,j,k,0) = 1.;//2*sin(2.*pi*y)*sin(pi*x)*sin(pi*x)*sin(2*pi*z)*cos(pi*0./3.);
+        vel(i,j,k,1) = 0.;//-sin(2.*pi*x)*sin(pi*y)*sin(pi*y)*sin(2*pi*z)*cos(pi*0./3.);
 #if (AMREX_SPACEDIM == 3)
-        vel(i,j,k,2) = -sin(2.*pi*x)*sin(pi*z)*sin(pi*z)*sin(2*pi*y)*cos(pi*0./3.);
+        vel(i,j,k,2) = 0.;//-sin(2.*pi*x)*sin(pi*z)*sin(pi*z)*sin(2*pi*y)*cos(pi*0./3.);
 #endif
     });
 }
diff --git a/src/vof/VolumeOfFluid.H b/src/vof/VolumeOfFluid.H
@@ -18,7 +18,7 @@ public:
                                            amrex::Vector<amrex::MultiFab const*> const& w_mac),
                        amrex::Real dt);
     void tracer_vof_init_fraction(int lev, amrex::MultiFab& a_tracer);
-    void tracer_vof_update(int lev, amrex::MultiFab & a_tracer);
+    void tracer_vof_update(int lev, amrex::MultiFab & a_tracer, amrex::Array<amrex::MultiFab, 2> & a_height);
     void write_tecplot_surface(amrex::Real time, int nstep);
     void WriteTecPlotFile     (amrex::Real time, int nstep);
     void output_droplet (amrex::Real time, int nstep);
@@ -28,6 +28,8 @@ public:
    amrex::Vector<amrex::MultiFab> normal;
 // the plane defined by @m.@x = @alpha
    amrex::Vector<amrex::MultiFab> alpha;
+// the heights used for calculation of curvature
+   amrex::Vector<amrex::Array<amrex::MultiFab, 2>> height;
 // the tag for droplets
    amrex::Vector<amrex::MultiFab> tag;
 private: