opm/ifem-darcy

Darcy application built using IFEM

ifem-darcy's Issues

support 'dim' keyword in geometry tag

Like we do in the Poisson applilcation

Memory issue in Darcy-object

It seems like the Darcy-object remembers something from the last run, which affects the outcome of the next run.
If I run the following code

I get the following plot

and if I change p1 from p1=np.linspace(0, 4, 10) to p1=np.linspace(-4, 4, 10) the plot becomes

i.e. the QoI values do not change, even though the pressure values are changed.
However if I run the following code, where I make a new Darcy-object for each loop

I get the following plot

If I again change p1 from p1=np.linspace(0, 4, 10) to p1=np.linspace(-4, 4, 10) the plot becomes

Here we can see that the two plots agree in the area where the pressure is equal (i.e. p $\in$ [0,4]).

The XML-file used:

Code used to make geometry:

Secondary solution (grad u) is not computed correctly

<?xml version="1.0" encoding="UTF-8" standalone="yes"?>                        
<simulation>                                                                   
  <!-- General - geometry definitions !-->                                     
  <geometry dim="3">                                                           
    <patchfile>PipeFlow3D.g2</patchfile>                                       
    <raiseorder patch="1" u="0" v="0" w="1"/>                                  
    <refine patch="1" u="3" v="3" w="7"/>                                      
    <topologysets>                                                             
      <set name="Inflow" type="face">                                          
        <item patch="1">5</item>                                               
      </set>                                                                   
      <set name="Outflow" type="face">                                         
        <item patch="1">6</item>                                               
      </set>                                                                   
    </topologysets>                                                            
  </geometry>                                                                  
                                                                               
  <!-- General - boundary conditions !-->                                      
  <boundaryconditions>                                                         
    <dirichlet set="Inflow"  comp="1" type="constant"> 3 </dirichlet>          
    <dirichlet set="Outflow" comp="1" type="constant"> 0 </dirichlet>          
  </boundaryconditions>                                                        
                                                                               
  <!-- Problem-specific block !-->                                             
  <darcy>                                                                      
    <permvalues>1|1|1</permvalues>                                             
    <bodyforce>0|0|0</bodyforce>                                               
    <source>0</source>                                                         
  </darcy>                                                                     
                                                                               
  <postprocessing>                                                             
    <vtfformat nviz="2">ascii</vtfformat>                                      
  </postprocessing>                                                            
                                                                               
</simulation>                                                                  

Where PipeFlow3D.g2 is a straight cylinder:

700 1 0 0                                                                      
3 1                                                                            
3 3                                                                            
0.000000 0.000000 0.000000 1.000000 1.000000 1.000000                          
3 3                                                                            
0.000000 0.000000 0.000000 1.000000 1.000000 1.000000                          
2 2                                                                            
0.000000 0.000000 1.000000 1.000000                                            
-0.000000 -1.414214 1.414214 1.000000                                          
0.000000 -2.000000 0.000000 0.707107                                           
0.000000 -1.414214 -1.414214 1.000000                                          
-0.000000 0.000000 2.000000 0.707107                                           
0.000000 0.000000 0.000000 1.000000                                            
0.000000 0.000000 -2.000000 0.707107                                           
-0.000000 1.414214 1.414214 1.000000                                           
0.000000 2.000000 0.000000 0.707107                                            
0.000000 1.414214 -1.414214 1.000000                                           
18.000000 -1.414214 1.414214 1.000000                                          
12.727922 -2.000000 0.000000 0.707107                                          
18.000000 -1.414214 -1.414214 1.000000                                         
12.727922 0.000000 2.000000 0.707107                                           
18.000000 0.000000 0.000000 1.000000                                           
12.727922 0.000000 -2.000000 0.707107                                          
18.000000 1.414214 1.414214 1.000000                                           
12.727922 2.000000 0.000000 0.707107                                           
18.000000 1.414214 -1.414214 1.000000                                          

The pressure should form a linear drop from left to right with a gradient that is constant [a,0,0] for some value a. Instead, this is what I get: