Wrong behaviour for static AD (Manni et.al 2011)

[VikingScientist]

<?xml version="1.0" encoding="UTF-8" standalone="yes"?>                        
                                                                               
<simulation>                                                                   
                                                                               
  <!-- General - geometry definitions !-->                                     
  <geometry dim="2" sets="True">                                               
    <raiseorder patch="1" u="2" v="2"/>                                        
    <refine type="uniform" patch="1" u="9" v="9"/>                             
    <topologysets/>                                                            
  </geometry>                                                                  
                                                                               
  <!-- General - boundary conditions !-->                                      
  <boundaryconditions>                                                         
    <dirichlet set="Boundary" comp="1"/>                                       
  </boundaryconditions>                                                        
                                                                               
  <discretization>                                                             
    <nGauss> 5 </nGauss>                                                       
  </discretization>                                                            
                                                                               
  <!-- Problem-specific block !-->                                             
  <advectiondiffusion>                                                         
    <advectionfield>1|0</advectionfield>                                       
    <!--                                                                       
    <stabilization> supg </stabilization>                                      
    <reactionfield> 1 </reactionfield>                                         
    -->                                                                        
    <source type="expression">2*x*cos(PI*y)</source>                           
    <kappa> 1e-2 </kappa>                                                      
  </advectiondiffusion>                                                        
                                                                               
  <postprocessing>                                                             
    <vtfformat nviz="2">ascii</vtfformat>                                      
  </postprocessing>                                                            
                                                                               
</simulation>                                                                  

Two immedate errors:

• The resulting solution field looks like it is using stabilization (I requested no stabilization)
• The resulting solution field has the wrong magnitude: 0.02, I was expecting 0.6

[VikingScientist]

Three fixes:

1. Syntax error on stabilization term. It should be formatted as

<stabilization type="supg"/>

2. Also <kappa> should be wrapped into a <fluidproperties> tag.

3. The SUPG stabilization is for some reason scaled by Xi in https://github.com/OPM/IFEM-AdvectionDiffusion/blob/master/AdvectionDiffusion.C#L253