Changeset 25514 for issm/trunk-jpl/src/m/classes/flowequation.m

Timestamp:

09/02/20 12:00:13 (5 years ago)

Author:

Cheng Gong

Message:

NEW: Galerkin Least Squares stabilization for P1P1 elements in FS is added together with a Nitsches method for weakly imposing Dirichlet boundary condition on the bottom surface of the ice. Current issue: the nonlinear iteration may fail to converge if the number of vertical layers is too small(~5).

File:

• issm/trunk-jpl/src/m/classes/flowequation.m (modified) (5 diffs)

issm/trunk-jpl/src/m/classes/flowequation.m

@@ -11 +11 @@
                 isHO                           = 0;
                 isFS                           = 0;
+                isNitscheBC                    = 0;
+      FSNitscheGamma                 = 1e6;
                 fe_SSA                         = '';
                 fe_HO                          = '';
@@ -98 +100 @@
                         md = checkfield(md,'fieldname','flowequation.isHO','numel',[1],'values',[0 1]);
                         md = checkfield(md,'fieldname','flowequation.isFS','numel',[1],'values',[0 1]);
+                        md = checkfield(md,'fieldname','flowequation.isNitscheBC','numel',[1],'values',[0 1]);
+                        md = checkfield(md,'fieldname','flowequation.FSNitscheGamma','numel',[1], '>=', 0.);
                         md = checkfield(md,'fieldname','flowequation.fe_SSA','values',{'P1','P1bubble','P1bubblecondensed','P2','P2bubble'});
                         md = checkfield(md,'fieldname','flowequation.fe_HO' ,'values',{'P1','P1bubble','P1bubblecondensed','P1xP2','P2xP1','P2','P2bubble','P1xP3','P1xP4','P2xP4'});
@@ -144 +148 @@
                         fielddisplay(self,'isHO','is the Higher-Order (HO) approximation used ?');
                         fielddisplay(self,'isFS','are the Full-FS (FS) equations used ?');
+                        fielddisplay(self,'isNitscheBC','is weakly imposed condition used?');
+                        fielddisplay(self,'FSNitscheGamma','Gamma value for the Nitsche term, by default gamma=1e6?');
                         fielddisplay(self,'fe_SSA','Finite Element for SSA  ''P1'', ''P1bubble'' ''P1bubblecondensed'' ''P2''');
                         fielddisplay(self,'fe_HO', 'Finite Element for HO   ''P1'' ''P1bubble'' ''P1bubblecondensed'' ''P1xP2'' ''P2xP1'' ''P2''');
@@ -160 +166 @@
                         WriteData(fid,prefix,'object',self,'fieldname','isHO','format','Boolean');
                         WriteData(fid,prefix,'object',self,'fieldname','isFS','format','Boolean');
+                        WriteData(fid,prefix,'object',self,'fieldname','isNitscheBC','format','Boolean');
+                        WriteData(fid,prefix,'object',self,'fieldname','FSNitscheGamma','format','Double');
                         WriteData(fid,prefix,'object',self,'fieldname','fe_SSA','data',self.fe_SSA,'format','String');
                         WriteData(fid,prefix,'object',self,'fieldname','fe_HO' ,'data',self.fe_HO,'format','String');
@@ -182 +190 @@
                         writejsdouble(fid,[modelname '.flowequation.isHO'],self.isHO);
                         writejsdouble(fid,[modelname '.flowequation.isFS'],self.isFS);
-                        writejsstring(fid,[modelname '.flowequation.fe_SSA'],self.fe_SSA);
+         writejsstring(fid,[modelname '.flowequation.isNitscheBC'],self.isNitscheBC);
+         writejsstring(fid,[modelname '.flowequation.FSNitscheGamma'],self.FSNitscheGamma);
+         writejsstring(fid,[modelname '.flowequation.fe_SSA'],self.fe_SSA);
                         writejsstring(fid,[modelname '.flowequation.fe_HO'],self.fe_HO);
                         writejsstring(fid,[modelname '.flowequation.fe_FS'],self.fe_FS);