1) adding conj heat transfer to stagflowcatalyst

2) changing blockcat1d refinement criterion 3) changing sundials options in blockcat1d
NREL · Aug 14, 2024 · 0708ddc · 0708ddc
1 parent 76f7a2d
commit 0708ddc
Show file tree

Hide file tree

Showing 9 changed files with 24 additions and 19 deletions.
diff --git a/models/blockCatalyst1d/GNUmakefile b/models/blockCatalyst1d/GNUmakefile
@@ -22,7 +22,7 @@ ifeq ($(USE_SUNDIALS),TRUE)
 # A good check is to see if $(SUNDIALS_ROOT)/lib has a bunch of libsundials_ files
 # To run with sundials, enabled, please compile with USE_SUNDIALS = TRUE
 SUNDIALS_ROOT ?= $(TOP)/../../../../../sundials/instdir
-SUNDIALS_LIB_DIR ?= $(SUNDIALS_ROOT)/lib
+SUNDIALS_LIB_DIR ?= $(SUNDIALS_ROOT)/lib64
 
 USE_CVODE_LIBS ?= TRUE
 USE_ARKODE_LIBS ?= TRUE
@@ -31,7 +31,7 @@ DEFINES += -DAMREX_USE_SUNDIALS
 INCLUDE_LOCATIONS += $(SUNDIALS_ROOT)/include
 LIBRARY_LOCATIONS += $(SUNDIALS_LIB_DIR)
 
-LIBRARIES += -L$(SUNDIALS_LIB_DIR) -lsundials_cvode
-LIBRARIES += -L$(SUNDIALS_LIB_DIR) -lsundials_arkode
+LIBRARIES += -L$(SUNDIALS_LIB_DIR) -lsundials_cvode -lsundials_arkode
+LIBRARIES += -L$(SUNDIALS_LIB_DIR) -lsundials_core -lsundials_nvecserial
 LIBRARIES += -L$(SUNDIALS_LIB_DIR) -lsundials_nvecmanyvector
 endif
diff --git a/models/blockCatalyst1d/inputs_x b/models/blockCatalyst1d/inputs_x
@@ -37,7 +37,7 @@ mflo.lo_bc = 3 6 6
 mflo.hi_bc = 3 6 6
 
 # Tagging
-mflo.tagged_vars = volfrac S1 S2
+mflo.tagged_vars = volfrac
 mflo.S1_refine = 1e20
 mflo.S2_refine = 0.11
 mflo.S1_refinegrad = 0.1

diff --git a/models/blockCatalyst1d/inputs_y b/models/blockCatalyst1d/inputs_y
@@ -37,7 +37,7 @@ mflo.lo_bc = 6 3 6
 mflo.hi_bc = 6 3 6
 
 # Tagging
-mflo.tagged_vars = volfrac S1 S2
+mflo.tagged_vars = volfrac
 mflo.S1_refine = 1e20
 mflo.S2_refine = 0.11
 mflo.S1_refinegrad = 0.1

diff --git a/models/blockCatalyst1d/inputs_z b/models/blockCatalyst1d/inputs_z
@@ -37,7 +37,7 @@ mflo.lo_bc = 6 6 3
 mflo.hi_bc = 6 6 3
 
 # Tagging
-mflo.tagged_vars = volfrac S1 S2
+mflo.tagged_vars = volfrac
 mflo.S1_refine = 1e20
 mflo.S2_refine = 0.11
 mflo.S1_refinegrad = 0.1

diff --git a/models/blockCatalyst1d/verify_spec.py b/models/blockCatalyst1d/verify_spec.py
@@ -43,7 +43,7 @@
 #Plot solutions
 #=======================================
 fig,ax=plt.subplots(2,2,figsize=(8,4))
-ax[0][0].plot(x,exactsoln,'r-',label="Exact solution")
+ax[0][0].plot(x,exactsoln,'r-',label="Exact solution",linewidth=3)
 ax[0][0].plot(x,fld_S1,'k*',label="mflo",markersize=2)
 ax[0][0].legend(loc="best")
 

diff --git a/models/stagFlowCatalyst/Prob.H b/models/stagFlowCatalyst/Prob.H
@@ -89,27 +89,27 @@ void initdata(Box const& bx, Array4<Real> const& phi, GeometryData const& geomda
                 phi(i,j,k,RHOV_INDX) = phi(i,j,k,DENS_INDX)*phi(i,j,k,VELY_INDX);
                 phi(i,j,k,RHOW_INDX) = phi(i,j,k,DENS_INDX)*phi(i,j,k,VELZ_INDX);
 
-                if(vfrac==one)
+                if(vfrac>zeroval)
                 {
                    phi(i,j,k,TEMP_INDX) = mflo_thermo::get_t_from_rpc(
                         phi(i,j,k,DENS_INDX),phi(i,j,k,PRES_INDX),spec);
 
-                        phi(i,j,k,RHOE_INDX)  = phi(i,j,k,DENS_INDX)*
+                    phi(i,j,k,RHOE_INDX)  = phi(i,j,k,DENS_INDX)*
                     mflo_thermo::get_e_from_rpc(phi(i,j,k,DENS_INDX),phi(i,j,k,PRES_INDX),spec);
+
+                    phi(i,j,k,RHOE_INDX) += half*phi(i,j,k,DENS_INDX)*
+                                        ( phi(i,j,k,VELX_INDX)*phi(i,j,k,VELX_INDX)
+                                        + phi(i,j,k,VELY_INDX)*phi(i,j,k,VELY_INDX)
+                                        + phi(i,j,k,VELZ_INDX)*phi(i,j,k,VELZ_INDX) );
                 }
                 else
                 {
                    phi(i,j,k,TEMP_INDX) = mflo_thermo::get_t_from_rpc(
                         phi(i,j,k,DENS_INDX),phi(i,j,k,PRES_INDX),dummyspec);
 
-                   phi(i,j,k,RHOE_INDX)  = phi(i,j,k,DENS_INDX)*
-                    mflo_thermo::get_e_from_rpc(phi(i,j,k,DENS_INDX),phi(i,j,k,PRES_INDX),dummyspec);
+                   phi(i,j,k,RHOE_INDX)  = mflo_thermo::get_solid_rhoe_from_t(phi(i,j,k,TEMP_INDX));
                 }
 
-                phi(i,j,k,RHOE_INDX) += half*phi(i,j,k,DENS_INDX)*
-                                        ( phi(i,j,k,VELX_INDX)*phi(i,j,k,VELX_INDX)
-                                        + phi(i,j,k,VELY_INDX)*phi(i,j,k,VELY_INDX)
-                                        + phi(i,j,k,VELZ_INDX)*phi(i,j,k,VELZ_INDX) );
 
                 if(vfrac>zeroval)
                 {

diff --git a/models/stagFlowCatalyst/inputs b/models/stagFlowCatalyst/inputs
@@ -46,6 +46,8 @@ mflo.volfrac_refinegrad = 0.1
 mflo.vely_refine = 4.0
 mflo.vely_refinegrad = 1e20
 
+mflo.conjugate_heat_transfer=1
+
 # PLOTFILES
 amr.plot_file  = plt    # root name of plot file
 amr.plot_int   = 4000     # number of timesteps between plot files

diff --git a/models/stagFlowCatalyst/thermo.H b/models/stagFlowCatalyst/thermo.H
@@ -20,7 +20,9 @@ namespace mflo_thermo
     {
         //modify if conjugate heat transfer 
         //included
-        amrex::Real temp=300.0;
+        amrex::Real Cv_solid=100.0;
+        amrex::Real rho_solid=2000.0;
+        amrex::Real temp=rhoe/(rho_solid*Cv_solid);
         return(temp);
     }
 
@@ -29,8 +31,8 @@ namespace mflo_thermo
     {
         //modify if conjugate 
         //heat transfer included
-        amrex::Real Cv_solid=1.0;
-        amrex::Real rho_solid=1.0;
+        amrex::Real Cv_solid=100.0;
+        amrex::Real rho_solid=2000.0;
         return(rho_solid*Cv_solid*temp);
     }
     AMREX_GPU_HOST_DEVICE AMREX_INLINE

diff --git a/models/stagFlowCatalyst/userfuncs.H b/models/stagFlowCatalyst/userfuncs.H
@@ -33,8 +33,9 @@ namespace mflo_user_funcs
                 GpuArray<Real, AMREX_SPACEDIM> dx,
                 const Real time)
         {
+            Real volfrac=phi(i,j,k,VFRAC_INDX);
             transpcoeffs(i,j,k,VISC_INDX  ) = 1e-5;       //viscosity
-            transpcoeffs(i,j,k,THCOND_INDX) = 0.01; //thermal conductivity
+            transpcoeffs(i,j,k,THCOND_INDX) = 0.01*volfrac+50.0*(1-volfrac); //thermal conductivity
         }
 
     AMREX_GPU_DEVICE AMREX_INLINE