From de092a3a82ea436042321a30e2362d596412c4c8 Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Fri, 12 Jun 2020 19:43:26 +0200
Subject: [PATCH] raise gravity term interpolation degree by one
---
LDDsimulation/domainPatch.py | 8 +++++---
1 file changed, 5 insertions(+), 3 deletions(-)
diff --git a/LDDsimulation/domainPatch.py b/LDDsimulation/domainPatch.py
index 240db36..56ff411 100644
--- a/LDDsimulation/domainPatch.py
+++ b/LDDsimulation/domainPatch.py
@@ -1165,12 +1165,14 @@ class DomainPatch(df.SubDomain):
for phase in has_phases:
rho = df.Constant(self.densities[phase])
+ # V = self.pressure[phase]
+ # element = V.ufl_element()
self.gravity_term.update(
{phase: df.Expression('g*rho*x[1]',
domain=self.mesh, #
- degree=self.interpolation_degree,
- g=g,
- rho=rho
+ degree=self.interpolation_degree+1,
+ g=self.gravity_acceleration,
+ rho=self.densities[phase]
)}
)
--
GitLab