https://github.com/astroshaper/asteroidthermophysicalmodels.jl
Julia-based thermophysical modeling (TPM) tool for asteroids.
https://github.com/astroshaper/asteroidthermophysicalmodels.jl
asteroid asteroids astrodynamics comet comets julia julia-language julialang planet planets simulation solar-system space-science thermophysics
Last synced: 3 months ago
JSON representation
Julia-based thermophysical modeling (TPM) tool for asteroids.
- Host: GitHub
- URL: https://github.com/astroshaper/asteroidthermophysicalmodels.jl
- Owner: Astroshaper
- License: mit
- Created: 2020-10-24T13:01:14.000Z (about 5 years ago)
- Default Branch: main
- Last Pushed: 2025-08-18T15:35:13.000Z (3 months ago)
- Last Synced: 2025-08-19T12:31:25.988Z (3 months ago)
- Topics: asteroid, asteroids, astrodynamics, comet, comets, julia, julia-language, julialang, planet, planets, simulation, solar-system, space-science, thermophysics
- Language: Julia
- Homepage:
- Size: 76.4 MB
- Stars: 49
- Watchers: 4
- Forks: 5
- Open Issues: 8
-
Metadata Files:
- Readme: README.md
- Changelog: CHANGELOG.md
- License: LICENSE
- Roadmap: ROADMAP.md
Awesome Lists containing this project
README
# AsteroidThermoPhysicalModels.jl
[](https://Astroshaper.github.io/AsteroidThermoPhysicalModels.jl/stable)
[](https://Astroshaper.github.io/AsteroidThermoPhysicalModels.jl/dev)
[](https://github.com/Astroshaper/AsteroidThermoPhysicalModels.jl/actions?query=workflow%3ACI+branch%3Amain)
[](https://codecov.io/gh/Astroshaper/AsteroidThermoPhysicalModels.jl)
[](https://github.com/JuliaTesting/Aqua.jl)
`AsteroidThermoPhysicalModels.jl` is a comprehensive Julia-based toolkit for thermophysical modeling (TPM) of asteroids. It allows you to simulate the temperature distribution of asteroids and predict non-gravitational perturbations on their dynamics (Yarkovsky and YORP effects).
## 📚 Documentation
For detailed documentation, please visit:
- [Stable Documentation](https://Astroshaper.github.io/AsteroidThermoPhysicalModels.jl/stable)
- [Development Documentation](https://Astroshaper.github.io/AsteroidThermoPhysicalModels.jl/dev)
Sample notebooks are available in [Astroshaper-examples](https://github.com/Astroshaper/Astroshaper-examples).
## 🚀 Installation
```julia
using Pkg
Pkg.add("AsteroidThermoPhysicalModels")
using AsteroidThermoPhysicalModels
```
Or in the Julia REPL package mode:
```
julia> ] # Press ] to enter package mode
pkg> add AsteroidThermoPhysicalModels
```
## 🔍 Features
### Thermophysical Processes
- **Heat Conduction**: 1-dimensional heat diffusion in depth direction
- Multiple numerical solvers available (explicit Euler, implicit Euler, and Crank-Nicolson methods)
- **Self-Shadowing**: Local shadows cast by topography
- **Self-Heating**: Re-absorption of scattered and radiated photons by surrounding facets
- **Binary Systems**: Support for mutual shadowing (eclipses) and mutual heating between primary and secondary bodies
### Shape Models
- Supports Wavefront OBJ format (*.obj)
### Non-Gravitational Effects
- **Yarkovsky Effect**: Orbital perturbation due to asymmetric thermal emission
- **YORP Effect**: Rotational perturbation due to asymmetric thermal emission
### Coming Soon (v0.2.0)
- Surface roughness modeling integration with `AsteroidShapeModels.jl`
- Enhanced heat conduction solver validation and benchmarks
## 🌟 Example
Temperature distribution of asteroid Didymos and its satellite Dimorphos:
## 📖 Basic Usage
```julia
using AsteroidShapeModels
using AsteroidThermoPhysicalModels
# Load an asteroid shape model
# - `path/to/shape.obj` is the path to your OBJ file (mandatory)
# - `scale` : scale factor for the shape model (e.g., 1000 for km to m conversion)
shape = load_shape_obj("path/to/shape.obj"; scale=1000)
# Set thermal parameters
thermo_params = ThermoParams(
8.0 * 3600, # Rotation period [s]
0.05, # Thermal skin depth [m]
200.0, # Thermal inertia [J m⁻² K⁻¹ s⁻¹/²]
0.1, # Reflectance in visible light [-]
0.0, # Reflectance in thermal infrared [-]
0.9, # Emissivity [-]
0.5, # Depth of lower boundary [m]
0.0125, # Depth step width [m]
41 # Number of depth steps
)
# Create TPM model with solver selection
stpm = SingleAsteroidTPM(shape, thermo_params;
SELF_SHADOWING = true,
SELF_HEATING = true,
SOLVER = CrankNicolsonSolver(thermo_params), # Choose solver
BC_UPPER = RadiationBoundaryCondition(),
BC_LOWER = InsulationBoundaryCondition()
)
# Available solvers:
# - ExplicitEulerSolver(thermo_params) # Fast but requires small time steps
# - ImplicitEulerSolver(thermo_params) # Stable for any time step
# - CrankNicolsonSolver(thermo_params) # Best accuracy
# Run simulation - see documentation for complete examples
```
## 📊 Output
The package produces detailed output files including:
- Surface and subsurface temperature distributions
- Thermal forces and torques
- Energy conservation metrics
## 🤝 Contributing
Contributions are welcome! Please feel free to submit a Pull Request.