Broke out heat simulation prep calc

StarCoreSE · Jul 27, 2024 · 39b5cac · 39b5cac
1 parent 383e41f
commit 39b5cac
Show file tree

Hide file tree

Showing 3 changed files with 50 additions and 35 deletions.
diff --git a/ThermalDynamics/Data/Scripts/Thermodynamics/Settings.cs b/ThermalDynamics/Data/Scripts/Thermodynamics/Settings.cs
@@ -16,7 +16,7 @@ public class Settings
 		public const string Filename = "ThermodynamicsConfig.cfg";
 		public const string Name = "Thermodynamics";
 		public const bool Debug = true;
-		public const bool DebugBlockColors = true;
+		public const bool DebugBlockColors = false;
 
 		public static Settings Instance;
 
@@ -77,7 +77,7 @@ public static Settings GetDefaults()
 				EnableDamage = true,
 				Frequency = 60,
 				SimulationSpeed = 1,
-                SolarEnergy = 40000f, 
+                SolarEnergy = 1000f, 
 				EnvironmentalRaycastDistance = 5000f,
 				VaccumeRadiationStrength = 0.05f,
 				PresurizedAtmoConductivity = 0.026f,

diff --git a/ThermalDynamics/Data/Scripts/Thermodynamics/ThermalCell.cs b/ThermalDynamics/Data/Scripts/Thermodynamics/ThermalCell.cs
@@ -395,6 +395,12 @@ private void UpdateHeat()
         /// </summary>
         internal void Update()
         {
+            // TODO: wind/liquid convective cooling
+            //float viscosity = 
+            //float windIntensity = DirectionalRadiationIntensity(ref Grid.FrameWindDirection, ref Grid.WindNode);
+            //float v = Grid.FrameWindDirection.Length();
+            //float hc = 10.45f - v + 10f * (float)Math.Sqrt(v);
+
             UpdateFrameAndLastTemperature();
             float totalRadiation = CalculateTotalRadiation();
             float deltaTemperature = CalculateDeltaTemperature();
@@ -463,7 +469,7 @@ private void HandleCriticalTemperature()
 
         private void UpdateDebugVisuals()
         {
-            if (Settings.Debug && MyAPIGateway.Session.IsServer)
+            if (Settings.DebugBlockColors && MyAPIGateway.Session.IsServer)
             {
                 Vector3 color = Tools.GetTemperatureColor(Temperature);
                 if (Block.ColorMaskHSV != color)

diff --git a/ThermalDynamics/Data/Scripts/Thermodynamics/ThermalGrid.cs b/ThermalDynamics/Data/Scripts/Thermodynamics/ThermalGrid.cs
@@ -390,48 +390,58 @@ public override void UpdateBeforeSimulation()
 
         private void PrepareNextSimulationStep()
         {
-            SolarRadiationNode.Update();
-
-            FrameSolarOccluded = false;
-            FrameSolarDirection = MyVisualScriptLogicProvider.GetSunDirection();
-            FrameMatrix = Grid.WorldMatrix;
 
             Vector3D position = Grid.PositionComp.WorldAABB.Center;
-            PlanetManager.Planet p = PlanetManager.GetClosestPlanet(position);
+            PrepareEnvironmentTemprature(ref position);
+            PrepareSolarEnvironment(ref position);
+        }
+
+        private void PrepareEnvironmentTemprature(ref Vector3D position)
+        {
+            PlanetManager.Planet planet = PlanetManager.GetClosestPlanet(position);
+            if (planet == null) return;
+
 
             bool isUnderground = false;
-            if (p != null)
+            PlanetDefinition def = planet.Definition();
+            Vector3 local = position - planet.Position;
+            Vector3D surfacePointLocal = planet.Entity.GetClosestSurfacePointLocal(ref local);
+            isUnderground = local.LengthSquared() < surfacePointLocal.LengthSquared();
+            float airDensity = planet.Entity.GetAirDensity(position);
+            float windSpeed = planet.Entity.GetWindSpeed(position);
+
+
+            float ambient = def.UndergroundTemperature;
+            if (!isUnderground)
+            {
+                float dot = (float)Vector3D.Dot(Vector3D.Normalize(local), FrameSolarDirection);
+                ambient = def.NightTemperature + ((dot + 1f) * 0.5f * (def.DayTemperature - def.NightTemperature));
+            }
+            else
             {
-                PlanetDefinition def = p.Definition();
-                Vector3 local = position - p.Position;
-                Vector3D surfacePointLocal = p.Entity.GetClosestSurfacePointLocal(ref local);
-                isUnderground = local.LengthSquared() < surfacePointLocal.LengthSquared();
-                float airDensity = p.Entity.GetAirDensity(position);
-                float windSpeed = p.Entity.GetWindSpeed(position);
+                FrameSolarOccluded = true;
+            }
 
+            FrameAmbientTemprature = Math.Max(2.7f, ambient * airDensity);
+            float frameAmbiSquared = FrameAmbientTemprature * FrameAmbientTemprature;
+            FrameAmbientTempratureP4 = frameAmbiSquared * frameAmbiSquared;
+            FrameSolarDecay = 1 - def.SolarDecay * airDensity;
 
-                float ambient = def.UndergroundTemperature;
-                if (!isUnderground)
-                {
-                    float dot = (float)Vector3D.Dot(Vector3D.Normalize(local), FrameSolarDirection);
-                    ambient = def.NightTemperature + ((dot + 1f) * 0.5f * (def.DayTemperature - def.NightTemperature));
-                }
-                else
-                {
-                    FrameSolarOccluded = true;
-                }
+            //FrameWindDirection = Vector3.Cross(planet.GravityComponent.GetWorldGravityNormalized(position), planet.Entity.WorldMatrix.Forward).Normalized() * windSpeed;
+            //MySimpleObjectDraw.DrawLine(position, position + FrameWindDirection, MyStringId.GetOrCompute("Square"), ref color2, 0.1f);
 
-                FrameAmbientTemprature = Math.Max(2.7f, ambient * airDensity);
-                float frameAmbiSquared = FrameAmbientTemprature * FrameAmbientTemprature;
-                FrameAmbientTempratureP4 = frameAmbiSquared * frameAmbiSquared;
-                FrameSolarDecay = 1 - def.SolarDecay * airDensity;
+            //TODO: implement underground core temparatures
+        }
 
-                FrameWindDirection = Vector3.Cross(p.GravityComponent.GetWorldGravityNormalized(position), p.Entity.WorldMatrix.Forward).Normalized() * windSpeed;
+        private void PrepareSolarEnvironment(ref Vector3D position)
+        {
+            if (!Settings.Instance.EnableSolarHeat) return;
 
+            SolarRadiationNode.Update();
 
-                //TODO: implement underground core temparatures
-            }
-            if (FrameSolarOccluded) return;
+            FrameSolarOccluded = false;
+            FrameSolarDirection = MyVisualScriptLogicProvider.GetSunDirection();
+            FrameMatrix = Grid.WorldMatrix;
 
             LineD line = new LineD(position, position + (FrameSolarDirection * 15000000));
             List<MyLineSegmentOverlapResult<MyEntity>> results = new List<MyLineSegmentOverlapResult<MyEntity>>();
@@ -520,7 +530,6 @@ private void PrepareNextSimulationStep()
                 var color = (FrameSolarOccluded) ? Color.Red.ToVector4() : Color.White.ToVector4();
                 var color2 = Color.LightGoldenrodYellow.ToVector4();
                 MySimpleObjectDraw.DrawLine(position, position + (FrameSolarDirection * 15000000), MyStringId.GetOrCompute("Square"), ref color, 0.1f);
-                MySimpleObjectDraw.DrawLine(position, position + FrameWindDirection, MyStringId.GetOrCompute("Square"), ref color2, 0.1f);
             }
         }