Add tanh initialisation

shkodm · shkodm · commit 36056b3abac1 · 2025-12-15T16:52:35.000+10:00
diff --git a/models/multiphase/d3q27_pf_velocity/Dynamics.c.Rt b/models/multiphase/d3q27_pf_velocity/Dynamics.c.Rt
@@ -400,6 +400,10 @@ CudaDeviceFunction void specialCases_Init()
 		PhaseF = 1 - 0.5 *  ( tanh( 2.0 * ( X - Washburn_start ) / IntWidth ) -
 						   tanh( 2.0 * ( X - Washburn_end  )  / IntWidth ));
 	}
+    if (Tanh_init > 0) {
+        PhaseF = 0.5 + 0.5*( tanh( 2.0 * ( X - Tanh_init) / IntWidth ));
+    }
+
 }
 
 CudaDeviceFunction void Init_distributions()

Original file line number	Diff line number	Diff line change
`@@ -400,6 +400,10 @@ CudaDeviceFunction void specialCases_Init()`
`400`	`400`	`PhaseF = 1 - 0.5 * ( tanh( 2.0 * ( X - Washburn_start ) / IntWidth ) -`
`401`	`401`	`tanh( 2.0 * ( X - Washburn_end ) / IntWidth ));`
`402`	`402`	`}`
	`403`	`+ if (Tanh_init > 0) {`
	`404`	`+ PhaseF = 0.5 + 0.5( tanh( 2.0 ( X - Tanh_init) / IntWidth ));`
	`405`	`+ }`
	`406`	`+`
`403`	`407`	`}`
`404`	`408`
`405`	`409`	`CudaDeviceFunction void Init_distributions()`