kindred/create
1
1
Fork 0
Code Issues 37 Pull Requests Actions Packages Projects 9 Releases 2 Wiki Activity

[FEM] fix flow equation

Browse Source 
- the density must be output using the correct units. Due to this bug also the viscosity was written with wrong units-
- if there is a dynamic viscosity given for a material, use it directly
This commit is contained in:
Uwe
2023-03-17 02:56:55 +01:00
parent c1751bde4a
commit b1c848b8db
1 changed files with 12 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
src/Mod/Fem/femsolver/elmer/equations/flow_writer.py
View File

@@ -138,33 +138,27 @@ class Flowwriter:
for name in (n for n in refs if n in bodies):
self.write.material(name, "Name", m["Name"])
if "Density" in m:
self.write.material(
name, "Density",
self.write.getDensity(m)
)
density = self.write.convert(m["Density"], "M/L^3")
self.write.material(name, "Density", density)
if "ThermalConductivity" in m:
self.write.material(
name, "Heat Conductivity",
self.write.convert(m["ThermalConductivity"], "M*L/(T^3*O)")
)
if "KinematicViscosity" in m:
density = self.write.getDensity(m)
tConductivity = self.write.convert(m["ThermalConductivity"], "M*L/(T^3*O)")
self.write.material(name, "Heat Conductivity", tConductivity)
if "DynamicViscosity" in m:
dViscosity = self.write.convert(m["DynamicViscosity"], "M/(L*T)")
self.write.material(name, "Viscosity", dViscosity)
elif ("KinematicViscosity" in m) and ("Density" in m):
density = self.write.convert(m["Density"], "M/L^3")
kViscosity = self.write.convert(m["KinematicViscosity"], "L^2/T")
self.write.material(
name, "Viscosity", kViscosity * density)
self.write.material(name, "Viscosity", kViscosity * density)
if "ThermalExpansionCoefficient" in m:
value = self.write.convert(m["ThermalExpansionCoefficient"], "O^-1")
if value > 0:
self.write.material(
name, "Heat expansion Coefficient", value)
self.write.material(name, "Heat expansion Coefficient", value)
if "ReferencePressure" in m:
pressure = self.write.convert(m["ReferencePressure"], "M/(L*T^2)")
self.write.material(name, "Reference Pressure", pressure)
if "SpecificHeatRatio" in m:
self.write.material(
name, "Specific Heat Ratio",
float(m["SpecificHeatRatio"])
)
self.write.material(name, "Specific Heat Ratio", float(m["SpecificHeatRatio"]))
if "CompressibilityModel" in m:
self.write.material(
name, "Compressibility Model",