mirror/OpenSpace
1
0
Fork 0
mirror of https://github.com/OpenSpace/OpenSpace.git synced 2026-01-07 04:00:37 -06:00
Code Issues Packages Projects Releases Wiki Activity

Update boundary cases for habitable zone

Browse Source
This commit is contained in:
Emma Broman
2020-12-14 10:10:35 +01:00
parent 92085847f2
commit 9898714478
3 changed files with 19 additions and 13 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
modules/exoplanets/exoplanetshelper.cpp
View File

@@ -136,10 +136,16 @@ glm::dmat3 computeSystemRotation(glm::dvec3 starPosition) {
);
}
std::optional<glm::vec2> computeHabitableZone(float teff, float luminosity) {
// The formula only consider stars with teff in range [2600, 7200] K
if (teff > 7200.f || teff < 2600.f) {
return std::nullopt;
glm::vec2 computeHabitableZone(float teff, float luminosity) {
// Kopparapu's formula only considers stars with teff in range [2600, 7200] K.
// However, we want to use the formula for more stars, so add some flexibility to
// the teff boundaries
if (teff > 8000.f || teff < 2000.f) {
// For the other stars, use a method by Tom E. Morris:
// https://www.planetarybiology.com/calculating_habitable_zone.html
float inner = sqrtf(luminosity / 1.1f);
float outer = sqrtf(luminosity / 0.53f);
return glm::vec2(inner, outer);
}
struct Coefficients {
@@ -151,9 +157,9 @@ std::optional<glm::vec2> computeHabitableZone(float teff, float luminosity) {
// https://depts.washington.edu/naivpl/sites/default/files/HZ_coefficients.dat
constexpr Coefficients coefficients[] = {
// Inner boundary - Runaway greenhouse
{1.10700E+00, 1.33200E-04, 1.58000E-08, -8.30800E-12, -1.93100E-15},
{1.10700E+00f, 1.33200E-04f, 1.58000E-08f, -8.30800E-12f, -1.93100E-15f},
// Outer boundary - Maximum greenhouse
{3.56000E-01, 6.17100E-05, 1.69800E-09, -3.19800E-12, -5.57500E-16}
{3.56000E-01f, 6.17100E-05f, 1.69800E-09f, -3.19800E-12f, -5.57500E-16f}
};
const float tstar = teff - 5780.f;
@@ -168,7 +174,7 @@ std::optional<glm::vec2> computeHabitableZone(float teff, float luminosity) {
distances[i] = std::pow(luminosity / seff, 0.5f);
}
return std::optional<glm::vec2>(std::move(distances));
return distances;
}
std::string createIdentifier(std::string name) {

2
modules/exoplanets/exoplanetshelper.h
View File

@@ -117,7 +117,7 @@ glm::dmat3 computeSystemRotation(glm::dvec3 starPosition);
* \return A vec2 with the lower and upper boundary in atronomical units, if a habitable
zone could be computed. Otherwise an std::nullopt
*/
std::optional<glm::vec2> computeHabitableZone(float teff, float luminosity);
glm::vec2 computeHabitableZone(float teff, float luminosity);
// Create an identifier without whitespaces
std::string createIdentifier(std::string name);

10
modules/exoplanets/exoplanetsmodule_lua.inl
View File

@@ -413,13 +413,13 @@ void createExoplanetSystem(const std::string& starName) {
// Habitable Zone
bool hasTeff = !std::isnan(system.starData.teff);
bool hasLuminosity = !std::isnan(system.starData.luminosity);
std::optional<glm::vec2> zone = std::nullopt;
if (hasTeff && hasLuminosity) {
zone = computeHabitableZone(system.starData.teff, system.starData.luminosity);
}
const glm::vec2 zone = computeHabitableZone(
system.starData.teff,
system.starData.luminosity
);
if (zone.has_value()) {
float meanInclination = 0.f;
for (const ExoplanetDataEntry& p : system.planetsData) {
meanInclination += p.i;
@@ -428,7 +428,7 @@ void createExoplanetSystem(const std::string& starName) {
const glm::dmat4 rotation = computeOrbitPlaneRotationMatrix(meanInclination);
const glm::dmat3 rotationMat3 = static_cast<glm::dmat3>(rotation);
glm::vec2 limitsInMeter = zone.value() * AU;
glm::vec2 limitsInMeter = zone * AU;
float half = 0.5f * (limitsInMeter[1] - limitsInMeter[0]);
float center = limitsInMeter[0] + half;
float relativeOffset = half / center;