From ac71f797bb437eff2c43f911f0c3570a6e6cd285 Mon Sep 17 00:00:00 2001 From: Daniel Douglas Date: Mon, 4 Nov 2024 11:23:53 -0700 Subject: [PATCH] fix the units for the darcy velocity postprocessor --- source/postprocess/visualization/darcy_velocity.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/source/postprocess/visualization/darcy_velocity.cc b/source/postprocess/visualization/darcy_velocity.cc index 20fc20a78d3..c0485bc8831 100644 --- a/source/postprocess/visualization/darcy_velocity.cc +++ b/source/postprocess/visualization/darcy_velocity.cc @@ -87,9 +87,9 @@ namespace aspect const double fluid_density = fluid_out->fluid_densities[q]; const double permeability = fluid_out->permeabilities[q]; const Tensor<1,dim> solid_velocity = in.velocity[q]; - const Tensor<1,dim> darcy_velocity = solid_velocity - + const Tensor<1,dim> darcy_velocity = (solid_velocity - permeability / fluid_viscosity / porosity * gravity * - (solid_density - fluid_density) * velocity_scaling_factor; + (solid_density - fluid_density)) * velocity_scaling_factor; for (unsigned int k=0; k