Added normal BCs in input file

CaNS-World · Nov 7, 2024 · c2f0956 · c2f0956
1 parent 6442d7f
commit c2f0956
Show file tree

Hide file tree

Showing 21 changed files with 388 additions and 210 deletions.
diff --git a/examples/capillary_wave/input.nml b/examples/capillary_wave/input.nml
@@ -31,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'cap-wav-1d'
-cbcpsi(0:1,1:3) = 'P','P',  'P','P',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1., 10.
-mu12(1:2)       = 0.1, 0.01
-sigma           = 1.
-ka12(1:2)       = 0., 0.
-cp12(1:2)       = 0., 0.
-beta12(1:2)     = 0., 0.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 0.51
+inipsi            = 'cap-wav-1d'
+cbcpsi(0:1,1:3)   = 'P','P',  'P','P',  'N','N'
+cbcnor(0:1,1:3,1) = 'P','P',  'P','P',  'N','N'
+cbcnor(0:1,1:3,2) = 'P','P',  'P','P',  'N','N'
+cbcnor(0:1,1:3,3) = 'P','P',  'P','P',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1., 10.
+mu12(1:2)         =  0.1, 0.01
+sigma             =  1.
+ka12(1:2)         =  0., 0.
+cp12(1:2)         =  0., 0.
+beta12(1:2)       =  0., 0.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  0.51
 \
 
 &cudecomp

diff --git a/examples/deforming_circle/input.nml b/examples/deforming_circle/input.nml
@@ -31,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'bub2'
-cbcpsi(0:1,1:3) = 'N','N',  'P','P',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1., 1.
-mu12(1:2)       = 100., 100.
-sigma           = 1.
-ka12(1:2)       = 0., 0.
-cp12(1:2)       = 0., 0.
-beta12(1:2)     = 0., 0.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 1.0
+inipsi            = 'bub2'
+cbcpsi(0:1,1:3)   = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,1) = 'D','D',  'P','P',  'N','N'
+cbcnor(0:1,1:3,2) = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,3) = 'N','N',  'P','P',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1., 1.
+mu12(1:2)         =  100., 100.
+sigma             =  1.
+ka12(1:2)         =  0., 0.
+cp12(1:2)         =  0., 0.
+beta12(1:2)       =  0., 0.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  1.0
 \
 
 &cudecomp

diff --git a/examples/differentially_heated_cavity/input.nml b/examples/differentially_heated_cavity/input.nml
@@ -31,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'uni'
-cbcpsi(0:1,1:3) = 'N','N',  'P','P',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1., 1.
-mu12(1:2)       = 0.0008426149773176359,0.0008426149773176359 ! Pr/(Pr*Ra)**.5
-sigma           = 0.
-ka12(1:2)       = 0.0011867816581938534,0.0011867816581938534 ! 1./(Pr*Ra)**.5
-cp12(1:2)       = 1., 1.
-beta12(1:2)     = 1., 1.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 1.
+inipsi            = 'uni'
+cbcpsi(0:1,1:3)   = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,1) = 'D','D',  'P','P',  'N','N'
+cbcnor(0:1,1:3,2) = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,3) = 'N','N',  'P','P',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1., 1.
+mu12(1:2)         =  0.0008426149773176359,0.0008426149773176359 ! Pr/(Pr*Ra)**.5
+sigma             =  0.
+ka12(1:2)         =  0.0011867816581938534,0.0011867816581938534 ! 1./(Pr*Ra)**.5
+cp12(1:2)         =  1., 1.
+beta12(1:2)       =  1., 1.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  1.
 \
 
 &cudecomp

diff --git a/examples/droplet_in_couette_flow/input.nml b/examples/droplet_in_couette_flow/input.nml
@@ -21,6 +21,7 @@ bcpre(0:1,1:3)   =  0.,0.,   0.,0.,   0. ,0.
 bforce(1:3) = 0., 0., 0.
 gacc(1:3) = 0., 0., 0.
 dims(1:2) = 0, 0
+\
 
 &scalar
 inisca          = 'zer'
@@ -30,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'bub3'
-cbcpsi(0:1,1:3) = 'P','P',  'P','P',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1.,1.
-mu12(1:2)       = 0.003125, 0.003125! 0.25**2/20 (Re = 20)
-sigma           = 0.0078125 ! mu/0.25/0.1 (Ca = 0.1)
-ka12(1:2)       = 0.001,0.001
-cp12(1:2)       = 1., 1.
-beta12(1:2)     = 1., 1.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 1.
+inipsi            = 'bub3'
+cbcpsi(0:1,1:3)   = 'P','P',  'P','P',  'N','N'
+cbcnor(0:1,1:3,1) = 'P','P',  'P','P',  'N','N'
+cbcnor(0:1,1:3,2) = 'P','P',  'P','P',  'N','N'
+cbcnor(0:1,1:3,3) = 'P','P',  'P','P',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1.,1.
+mu12(1:2)         =  0.003125, 0.003125! 0.25**2/20 (Re = 20)
+sigma             =  0.0078125 ! mu/0.25/0.1 (Ca = 0.1)
+ka12(1:2)         =  0.001,0.001
+cp12(1:2)         =  1., 1.
+beta12(1:2)       =  1., 1.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  1.
 \
 
 &cudecomp

diff --git a/examples/forced_hit/input.nml b/examples/forced_hit/input.nml
@@ -21,6 +21,7 @@ bcpre(0:1,1:3)   =  0.,0.,   0.,0.,   0.,0.
 bforce(1:3) = 0., 0., 0.
 gacc(1:3) = 0., 0., 0.
 dims(1:2) = 0, 0
+\
 
 &scalar
 inisca          = 'zer'
@@ -30,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'uni'
-cbcpsi(0:1,1:3) = 'P','P',  'P','P',  'P','P'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1., 1.
-mu12(1:2)       = .002,.002
-sigma           = 0.
-ka12(1:2)       = 0.001,0.001
-cp12(1:2)       = 1., 1.
-beta12(1:2)     = 1., 1.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 0.51
+inipsi            = 'uni'
+cbcpsi(0:1,1:3)   = 'P','P',  'P','P',  'P','P'
+cbcnor(0:1,1:3,1) = 'P','P',  'P','P',  'P','P'
+cbcnor(0:1,1:3,2) = 'P','P',  'P','P',  'P','P'
+cbcnor(0:1,1:3,3) = 'P','P',  'P','P',  'P','P'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1., 1.
+mu12(1:2)         =  .002,.002
+sigma             =  0.
+ka12(1:2)         =  0.001,0.001
+cp12(1:2)         =  1., 1.
+beta12(1:2)       =  1., 1.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  0.51
 \
 
 &cudecomp

diff --git a/examples/lid_driven_cavity/input.nml b/examples/lid_driven_cavity/input.nml
@@ -21,6 +21,7 @@ bcpre(0:1,1:3)   =  0.,0.,   0.,0.,   0.,0.
 bforce(1:3) = 0., 0., 0.
 gacc(1:3) = 0., 0., 0.
 dims(1:2) = 0, 0
+\
 
 &scalar
 inisca          = 'zer'
@@ -30,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'uni'
-cbcpsi(0:1,1:3) = 'P','P',  'N','N',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1., 1.
-mu12(1:2)       = 0.01,0.01
-sigma           = 0.
-ka12(1:2)       = 0.01,0.01
-cp12(1:2)       = 1., 1.
-beta12(1:2)     = 1., 1.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 1.
+inipsi            = 'uni'
+cbcpsi(0:1,1:3)   = 'P','P',  'N','N',  'N','N'
+cbcnor(0:1,1:3,1) = 'P','P',  'N','N',  'N','N'
+cbcnor(0:1,1:3,2) = 'P','P',  'D','D',  'N','N'
+cbcnor(0:1,1:3,3) = 'P','P',  'N','N',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1., 1.
+mu12(1:2)         =  0.01,0.01
+sigma             =  0.
+ka12(1:2)         =  0.01,0.01
+cp12(1:2)         =  1., 1.
+beta12(1:2)       =  1., 1.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  1.
 \
 
 &cudecomp

diff --git a/examples/rising_bubble_2d/input-case-01.nml b/examples/rising_bubble_2d/input-case-01.nml
@@ -21,6 +21,7 @@ bcpre(0:1,1:3)   =  0.,0.,   0.,0.,   0.,0.
 bforce(1:3) = 0., 0., 0.
 gacc(1:3) = 0., 0., -0.98
 dims(1:2) = 0, 0
+\
 
 &scalar
 inisca          = 'zer'
@@ -30,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'bub2'
-cbcpsi(0:1,1:3) = 'N','N',  'P','P',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1000.,100.
-mu12(1:2)       = 10.,1.
-sigma           = 24.5
-ka12(1:2)       = 0.001,0.001
-cp12(1:2)       = 1., 1.
-beta12(1:2)     = 1., 1.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 1.
+inipsi            = 'bub2'
+cbcpsi(0:1,1:3)   = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,1) = 'D','D',  'P','P',  'N','N'
+cbcnor(0:1,1:3,2) = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,3) = 'N','N',  'P','P',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1000.,100.
+mu12(1:2)         =  10.,1.
+sigma             =  24.5
+ka12(1:2)         =  0.001,0.001
+cp12(1:2)         =  1., 1.
+beta12(1:2)       =  1., 1.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  1.
 \
 
 &cudecomp

diff --git a/examples/rising_bubble_2d/input-case-02.nml b/examples/rising_bubble_2d/input-case-02.nml
@@ -21,6 +21,7 @@ bcpre(0:1,1:3)   =  0.,0.,   0.,0.,   0.,0.
 bforce(1:3) = 0., 0., 0.
 gacc(1:3) = 0., 0., -0.98
 dims(1:2) = 0, 0
+\
 
 &scalar
 inisca          = 'zer'
@@ -30,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'bub2'
-cbcpsi(0:1,1:3) = 'N','N',  'P','P',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1000.,1.
-mu12(1:2)       = 10.,0.1
-sigma           = 1.96
-ka12(1:2)       = 0.001,0.001
-cp12(1:2)       = 1., 1.
-beta12(1:2)     = 1., 1.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 1.
+inipsi            = 'bub2'
+cbcpsi(0:1,1:3)   = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,1) = 'D','D',  'P','P',  'N','N'
+cbcnor(0:1,1:3,2) = 'N','N',  'P','P',  'N','N'
+cbcnor(0:1,1:3,3) = 'N','N',  'P','P',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1000.,1.
+mu12(1:2)         =  10.,0.1
+sigma             =  1.96
+ka12(1:2)         =  0.001,0.001
+cp12(1:2)         =  1., 1.
+beta12(1:2)       =  1., 1.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  1.
 \
 
 &cudecomp

diff --git a/examples/rising_bubble_3d/input-case-01.nml b/examples/rising_bubble_3d/input-case-01.nml
@@ -21,6 +21,7 @@ bcpre(0:1,1:3)   =  0.,0.,   0.,0.,   0.,0.
 bforce(1:3) = 0., 0., 0.
 gacc(1:3) = 0., 0., -0.98
 dims(1:2) = 0, 0
+\
 
 &scalar
 inisca          = 'zer'
@@ -30,17 +31,23 @@ ssource         = 0.
 \
 
 &two_fluid
-inipsi          = 'bub3'
-cbcpsi(0:1,1:3) = 'N','N',  'N','N',  'N','N'
-bcpsi(0:1,1:3)  =  0.,0. ,   0.,0. ,   0.,0.
-rho12(1:2)      = 1000.,100.
-mu12(1:2)       = 10.,1.
-sigma           = 24.5
-ka12(1:2)       = 0.001,0.001
-cp12(1:2)       = 1., 1.
-beta12(1:2)     = 1., 1.
-acdi_gam_factor = 1., acdi_gam_min = 0.
-acdi_eps_factor = 1.
+inipsi            = 'bub3'
+cbcpsi(0:1,1:3)   = 'N','N',  'N','N',  'N','N'
+cbcnor(0:1,1:3,1) = 'D','D',  'N','N',  'N','N'
+cbcnor(0:1,1:3,2) = 'N','N',  'D','D',  'N','N'
+cbcnor(0:1,1:3,3) = 'N','N',  'N','N',  'D','D'
+bcpsi(0:1,1:3)    =  0.,0. ,   0.,0. ,   0.,0.
+bcnor(0:1,1:3,1)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,2)  =  0.,0.,    0.,0.,    0.,0.
+bcnor(0:1,1:3,3)  =  0.,0.,    0.,0.,    0.,0.
+rho12(1:2)        =  1000.,100.
+mu12(1:2)         =  10.,1.
+sigma             =  24.5
+ka12(1:2)         =  0.001,0.001
+cp12(1:2)         =  1., 1.
+beta12(1:2)       =  1., 1.
+acdi_gam_factor   =  1., acdi_gam_min = 0.
+acdi_eps_factor   =  1.
 \
 
 &cudecomp