Update test scripts

Author: MasanoriKanamaru

test/TPM_Didymos/TPM_Didymos.jl

@@ -54,7 +54,7 @@
     P₂ = SPICE.convrt(11.93 , "hours", "seconds")  # Rotation period of Dimorphos
 
     n_cycle = 2  # Number of cycles to perform TPM
-    n_step_in_cycle = 72  # Number of time steps in one rotation period
+    n_step_in_cycle = 120  # Number of time steps in one rotation period
 
     et_begin = SPICE.utc2et("2027-02-18T00:00:00")  # Start time of TPM
     et_end   = et_begin + P₂ * n_cycle  # End time of TPM
@@ -94,10 +94,6 @@
     ρ  = 2170.0  # Density [kg/m³]
     Cₚ = 600.0   # Heat capacity [J/kg/K]
 
-    l₁ = AsteroidThermoPhysicalModels.thermal_skin_depth(P₁, k, ρ, Cₚ)  # Thermal skin depth for Didymos
-    l₂ = AsteroidThermoPhysicalModels.thermal_skin_depth(P₂, k, ρ, Cₚ)  # Thermal skin depth for Dimorphos
-    Γ = AsteroidThermoPhysicalModels.thermal_inertia(k, ρ, Cₚ)          # Thermal inertia for Didymos and Dimorphos [tiu]
-
     R_vis = 0.059  # Reflectance in visible light [-]
     R_ir  = 0.0    # Reflectance in thermal infrared [-]
     ε     = 0.9    # Emissivity [-]
@@ -106,8 +102,8 @@
     n_depth = 41  # Number of depth steps
     Δz = z_max / (n_depth - 1)  # Depth step width [m]
 
-    thermo_params1 = AsteroidThermoPhysicalModels.ThermoParams(P₁, l₁, Γ, R_vis, R_ir, ε, z_max, Δz, n_depth)
-    thermo_params2 = AsteroidThermoPhysicalModels.ThermoParams(P₂, l₂, Γ, R_vis, R_ir, ε, z_max, Δz, n_depth)
+    thermo_params1 = AsteroidThermoPhysicalModels.ThermoParams(k, ρ, Cₚ, R_vis, R_ir, ε, z_max, Δz, n_depth)
+    thermo_params2 = AsteroidThermoPhysicalModels.ThermoParams(k, ρ, Cₚ, R_vis, R_ir, ε, z_max, Δz, n_depth)
 
     ##= Setting of TPM =##
     stpm1 = AsteroidThermoPhysicalModels.SingleTPM(shape1, thermo_params1;

test/TPM_Ryugu/TPM_Ryugu.jl

@@ -46,7 +46,7 @@
     P = SPICE.convrt(7.63262, "hours", "seconds")  # Rotation period of Ryugu
 
     n_cycle = 2  # Number of cycles to perform TPM
-    n_step_in_cycle = 72  # Number of time steps in one rotation period
+    n_step_in_cycle = 120  # Number of time steps in one rotation period
 
     et_begin = SPICE.utc2et("2018-07-01T00:00:00")  # Start time of TPM
     et_end   = et_begin + P * n_cycle  # End time of TPM
@@ -74,7 +74,7 @@
     k  = 0.1     # Thermal conductivity [W/m/K]
     ρ  = 1270.0  # Density [kg/m³]
     Cₚ = 600.0   # Heat capacity [J/kg/K]
-    
+
     l = AsteroidThermoPhysicalModels.thermal_skin_depth(P, k, ρ, Cₚ)  # Thermal skin depth [m]
     Γ = AsteroidThermoPhysicalModels.thermal_inertia(k, ρ, Cₚ)        # Thermal inertia [tiu]
 
@@ -83,10 +83,10 @@
     ε     = 1.0   # Emissivity [-]
 
     z_max = 0.6   # Depth of the lower boundary of a heat conduction equation [m]
-    n_depth = 41  # Number of depth steps
+    n_depth = 61  # Number of depth steps
     Δz = z_max / (n_depth - 1)  # Depth step width [m]
 
-    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(P, l, Γ, R_vis, R_ir, ε, z_max, Δz, n_depth)
+    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(k, ρ, Cₚ, R_vis, R_ir, ε, z_max, Δz, n_depth)
 
     ##= Setting of TPM =##
     stpm = AsteroidThermoPhysicalModels.SingleTPM(shape, thermo_params;

test/TPM_non-uniform_thermoparams/TPM_non-uniform_thermoparams.jl

@@ -87,21 +87,18 @@
     """
 
     k  = [r[3] > 0 ? 0.1 : 0.3  for r in shape.face_centers]  # Thermal conductivity [W/m/K]
-    ρ  = 1270.0  # Density [kg/m³]
-    Cₚ = 600.0   # Heat capacity [J/kg/K]
-    
-    l = AsteroidThermoPhysicalModels.thermal_skin_depth(P, k, ρ, Cₚ)  # Thermal skin depth [m]
-    Γ = AsteroidThermoPhysicalModels.thermal_inertia(k, ρ, Cₚ)        # Thermal inertia [tiu]
+    ρ  = fill(1270.0, n_face)  # Density [kg/m³]
+    Cₚ = fill(600.0, n_face)   # Heat capacity [J/kg/K]
 
     R_vis = [r[3] > 0 ? 0.04 : 0.1 for r in shape.face_centers]  # Reflectance in visible light [-]
-    R_ir  = [r[3] > 0 ? 1.0 : 0.9  for r in shape.face_centers]  # Reflectance in thermal infrared [-]
+    R_ir  = [r[3] > 0 ? 0.0 : 0.0  for r in shape.face_centers]  # Reflectance in thermal infrared [-]
     ε     = fill(1.0, n_face)   # Emissivity [-]
 
     z_max = 0.6   # Depth of the lower boundary of a heat conduction equation [m]
     n_depth = 41  # Number of depth steps
     Δz = z_max / (n_depth - 1)  # Depth step width [m]
 
-    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(P, l, Γ, R_vis, R_ir, ε, z_max, Δz, n_depth)
+    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(k, ρ, Cₚ, R_vis, R_ir, ε, z_max, Δz, n_depth)
 
     ##= Setting of TPM =##
     stpm = AsteroidThermoPhysicalModels.SingleTPM(shape, thermo_params;

test/TPM_zero-conductivity/TPM_zero-conductivity.jl

@@ -35,8 +35,9 @@
     n_face = length(shape.faces)  # Number of faces
 
     ##= Thermal properties: zero-conductivity case =##
-    l = 0.0  # Thermal skin depth [m]
-    Γ = 0.0  # Thermal inertia [J m⁻² K⁻¹ s⁻¹/2]
+    k  = 0.0     # Thermal conductivity [W/m/K]
+    ρ  = 1270.0  # Density [kg/m³]
+    Cₚ = 600.0   # Heat capacity [J/kg/K]
 
     R_vis = 0.1  # Reflectance in visible light [-]
     R_ir  = 0.0  # Reflectance in thermal infrared [-]
@@ -46,7 +47,7 @@
     n_depth = 41  # Number of depth steps
     Δz = z_max / (n_depth - 1)  # Depth step width [m]
 
-    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(P, l, Γ, R_vis, R_ir, ε, z_max, Δz, n_depth)
+    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(k, ρ, Cₚ, R_vis, R_ir, ε, z_max, Δz, n_depth)
 
     ##= Setting of TPM =##
     stpm = AsteroidThermoPhysicalModels.SingleTPM(shape, thermo_params;

test/heat_conduction_1D.jl

@@ -20,13 +20,9 @@
     )
 
     ##= Thermal properties =##
-    P  = 1.0
-    k  = 1.0
-    ρ  = 1.0
-    Cₚ = 1.0
-    
-    l = AsteroidThermoPhysicalModels.thermal_skin_depth(P, k, ρ, Cₚ)
-    Γ = AsteroidThermoPhysicalModels.thermal_inertia(k, ρ, Cₚ)
+    k  = 0.1     # Thermal conductivity [W/m/K]
+    ρ  = 1270.0  # Density [kg/m³]
+    Cₚ = 600.0   # Heat capacity [J/kg/K]
 
     R_vis = 0.0  # Reflectance in visible light [-]
     R_ir  = 0.0  # Reflectance in thermal infrared [-]
@@ -36,7 +32,7 @@
     n_depth = 101  # Number of depth steps
     Δz = z_max / (n_depth - 1)  # Depth step width [m]
 
-    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(P, l, Γ, R_vis, R_ir, ε, z_max, Δz, n_depth)
+    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(k, ρ, Cₚ, R_vis, R_ir, ε, z_max, Δz, n_depth)
 
     ##= TPMs with different solvers =##
     SELF_SHADOWING = false

test/thermal_radiation.jl

@@ -68,7 +68,7 @@
     n_depth = 41  # Number of depth steps
     Δz = z_max / (n_depth - 1)  # Depth step width [m]
 
-    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(P, l, Γ, R_vis, R_ir, ε, z_max, Δz, n_depth)
+    thermo_params = AsteroidThermoPhysicalModels.ThermoParams(k, ρ, Cₚ, R_vis, R_ir, ε, z_max, Δz, n_depth)
 
     ##= Setting of TPM =##
     stpm = AsteroidThermoPhysicalModels.SingleTPM(shape, thermo_params;