asd

kalman/input_files/aorta.yaml

@@ -150,4 +150,4 @@ estimation:
         reparameterize: True
         ODV_functional:
             enable: False
-            VENC: 44    # cm/s  ref: 61 for 70% and 44 for 50%
+            VENC: 61    # cm/s  ref: 61 for 70% and 44 for 50%

kalman/input_files/channel3d.yaml

@@ -5,12 +5,12 @@ fluid:
     dynamic_viscosity: 0.035
 
 io:
-    write_path: 'results/channel3d/60/'
+    write_path: 'results/channel3d/'
     restart:
         path: ''  # './projects/nse_coa3d/results/test_restart2/' 
         time: 0
     write_xdmf: True
-    write_checkpoints: False
+    write_checkpoints: True
     write_hdf5_timeseries: False
     write_velocity: 'update'     # update or tentative
 
@@ -26,7 +26,7 @@ boundary_conditions:
         type: 'dirichlet'
         value: ['0','0','U*(1-x[0]*x[0] - x[1]*x[1])*sin(DOLFIN_PI*t/T)']
         parameters:
-            U: 60
+            U: 30
             T: 0.9
             t: 0
     -    
@@ -106,14 +106,17 @@ estimation:
     measurements:
         -
             mesh: './meshes/channel3d.h5'
-            fe_degree: 0
-            xdmf_file: 'results/channel3d/measurements/u_all.xdmf'
-            file_root: 'results/channel3d/measurements/u{i}.h5'
+            fe_degree: 1
+            xdmf_file: 'results/channel3d/measurements/Perturbation/SNRinfV50/u_all.xdmf'
+            file_root: 'results/channel3d/measurements/Perturbation/SNRinfV50/u{i}.h5'
             indices: 0 # indices of checkpoints to be processed. 0 == all
             velocity_direction: ~
-            noise_stddev: 10.0     # standard deviation of Gaussian noise
+            noise_stddev: 2.0     # standard deviation of Gaussian noise
     
     roukf:
         particles: 'simplex' # unique or simplex
         observation_operator: 'postprocessing'         #state or postprocessing
         reparameterize: True
+        ODV_functional:
+            enable: True
+            VENC: 18    # cm/s  ref: 25 for 70% and 18 for 50%