add pressure drop for unsteady

2019-10-22 13:03:00 +02:00
parent ecbf2581a0
commit d2f62bd4dc
1 changed files with 7 additions and 0 deletions
--- a/analysis/SubPhase.cpp
+++ b/analysis/SubPhase.cpp
@@ -273,6 +273,13 @@ void SubPhase::Basic(){
 		double dir_x = Fx/force_mag;
 		double dir_y = Fy/force_mag;
 		double dir_z = Fz/force_mag;
+		
+		if (Dm->BoundaryCondition > 0 ){
+			// compute the pressure drop
+			double pressure_drop = (Pressure(Nx*Ny + Nx + 1) - 1.0) / 3.0;
+			force_mag += pressure_drop;
+		}
+		
 		if (force_mag == 0.0){
 			// default to z direction
 			dir_x = 0.0;