private Planet GetTestPlanetAtmosphere()
{
var planet = new Planet();
planet.Star = GetTestStar();
planet.Star.EcosphereRadiusAU = System.Math.Sqrt(planet.Star.Luminosity);
planet.SemiMajorAxisAU = 0.723332;
planet.Eccentricity = 0.0067;
planet.AxialTiltDegrees = 2.8;
planet.OrbitZone = Environment.OrbitalZone(planet.Star.Luminosity, planet.SemiMajorAxisAU);
planet.DayLengthHours = 2802;
planet.OrbitalPeriodDays = 225;
planet.MassSM = 0.000002447;
planet.GasMassSM = 2.41E-10;
planet.DustMassSM = planet.MassSM - planet.GasMassSM;
planet.RadiusKM = 6051.8;
planet.DensityGCC = Environment.EmpiricalDensity(planet.MassSM, planet.SemiMajorAxisAU, planet.Star.EcosphereRadiusAU, true);
planet.ExosphereTempKelvin = GlobalConstants.EARTH_EXOSPHERE_TEMP / Utilities.Pow2(planet.SemiMajorAxisAU / planet.Star.EcosphereRadiusAU);
planet.SurfaceAccelerationCMSec2 = Environment.Acceleration(planet.MassSM, planet.RadiusKM);
planet.EscapeVelocityCMSec = Environment.EscapeVelocity(planet.MassSM, planet.RadiusKM);
planet.IsGasGiant = false;
planet.Atmosphere.SurfacePressure = 92000;
planet.DaytimeTempKelvin = 737;
planet.NighttimeTempKelvin = 737;
planet.SurfaceTempKelvin = 737;
planet.SurfaceGravityG = 0.9;
planet.MolecularWeightRetained = Environment.MinMolecularWeight(planet);
return(planet);
}
public void GetAngularVelocity()
{
Mass massSM = Mass.FromEarthMasses(1.0);
Length radiusKM = Length.FromKilometers(10000.0);
Density densityGCC = Density.FromGramsPerCubicCentimeter(2.0);
Length semiMajorAxisAU = Length.FromAstronomicalUnits(1.0);
bool isGasGiant = false;
Mass largeMassSM = Mass.FromSolarMasses(1.0);
Duration largeAgeYears = Duration.FromYears365(1e10);
Assert.AreEqual(Env.AngularVelocity(massSM.SolarMasses,
radiusKM.Kilometers,
densityGCC.GramsPerCubicCentimeter,
semiMajorAxisAU.AstronomicalUnits,
isGasGiant,
largeMassSM.SolarMasses,
largeAgeYears.Years365),
_phy.Dynamics.GetAngularVelocity(massSM,
radiusKM,
densityGCC,
semiMajorAxisAU,
isGasGiant,
largeMassSM,
largeAgeYears).RadiansPerSecond,
LowDelta);
}
public void GetOrbitalZone()
{
Luminosity luminosity = Luminosity.FromSolarLuminosities(1.0);
Length orbitalRadius = Length.FromAstronomicalUnits(1.0);
Assert.AreEqual(Env.OrbitalZone(luminosity.SolarLuminosities, orbitalRadius.AstronomicalUnits),
_phy.Astronomy.GetOrbitalZone(luminosity, orbitalRadius)
);
}
public void GetMinimumIllumination()
{
Length a = Length.FromAstronomicalUnits(1.0);
Luminosity l = Luminosity.FromSolarLuminosities(1.0);
Assert.AreEqual(Env.MinimumIllumination(a.AstronomicalUnits, l.SolarLuminosities),
_phy.Astronomy.GetMinimumIllumination(a, l).DecimalFractions,
LowDelta);
}
public void GetPeriod()
{
Length separation = Length.FromAstronomicalUnits(1.0);
Mass smallMass = Mass.FromEarthMasses(1.0);
Mass largeMass = Mass.FromSolarMasses(1.0);
Assert.AreEqual(Env.Period(separation.AstronomicalUnits, smallMass.SolarMasses, largeMass.SolarMasses),
_phy.Astronomy.GetPeriod(separation, smallMass, largeMass).Days,
LowDelta);
}
public void GetHillSphere()
{
Mass sunMass = Mass.FromSolarMasses(1.0);
Mass massSM = Mass.FromEarthMasses(1.0);
Length semiMajorAxisAU = Length.FromAstronomicalUnits(1.0);
Assert.AreEqual(Env.SimplifiedHillSphereAU(sunMass.SolarMasses, massSM.SolarMasses, semiMajorAxisAU.AstronomicalUnits),
_phy.Astronomy.GetHillSphere(sunMass, massSM, semiMajorAxisAU).AstronomicalUnits,
LowDelta);
}
public void GetDayLength()
{
RotationalSpeed angularVelocityRadSec = RotationalSpeed.FromRadiansPerSecond(10.0);
Duration orbitalPeriod = Duration.FromDays(600.0);
Ratio eccentricity = Ratio.FromDecimalFractions(0.5);
Assert.AreEqual(Env.DayLength(angularVelocityRadSec.RadiansPerSecond, orbitalPeriod.Days, eccentricity.Value),
_phy.Astronomy.GetDayLength(angularVelocityRadSec, orbitalPeriod, eccentricity).Hours,
LowDelta);
}
public void GetRMSVelocity()
{
Mass molecularWeight = Mass.FromGrams(14.0);
Temperature exoTemp = Temperature.FromKelvins(5000.0);
Assert.AreEqual(Env.RMSVelocity(molecularWeight.Grams,
exoTemp.Kelvins),
_phy.Physics.GetRMSVelocity(molecularWeight,
exoTemp).CentimetersPerSecond,
LowDelta);
}
public void GetEscapeVelocity()
{
Mass massSM = Mass.FromEarthMasses(2.0);
Length radius = Length.FromKilometers(10000.0);
Assert.AreEqual(Env.EscapeVelocity(massSM.SolarMasses,
radius.Kilometers) / 1e5,
_phy.Dynamics.GetEscapeVelocity(massSM,
radius).KilometersPerSecond,
VeryHighDelta);
}
public void GetInitialMolecularWeightRetained()
{
SatelliteBody p1 = GetTestPlanetAtmosphere();
Data.Planet p2 = GetTestPlanetAtmosphereStarform();
Assert.AreEqual(Env.MoleculeLimit(p2.MassSM, p2.RadiusKM, p2.ExosphereTempKelvin),
_phy.Physics.GetInitialMolecularWeightRetained(p1.Mass,
p1.Radius,
p1.ExosphereTemperature).Grams,
VeryHighDelta);
}
public void GetMolecularWeightRetained()
{
SatelliteBody p1 = GetTestPlanetAtmosphere();
Data.Planet p2 = GetTestPlanetAtmosphereStarform();
Assert.AreEqual(Env.MinMolecularWeight(p2),
_phy.Physics.GetMolecularWeightRetained(p1.SurfaceAcceleration,
p1.Mass,
p1.Radius,
p1.ExosphereTemperature,
p1.Parent.Age).Grams,
UnitDelta);
}
public void GetBaseAngularVelocity()
{
Mass massSM = Mass.FromSolarMasses(1.0);
Length radiusKM = Length.FromKilometers(5000.0);
bool isGasGiant = false;
Assert.AreEqual(Env.BaseAngularVelocity(massSM.SolarMasses,
radiusKM.Kilometers,
isGasGiant),
_phy.Dynamics.GetBaseAngularVelocity(massSM,
radiusKM,
isGasGiant).RadiansPerSecond,
HighDelta);
}
public void GetGasLife()
{
//SatelliteBody p1 = GetTestPlanetAtmosphere();
//Data.Planet p2 = GetTestPlanetAtmosphereStarform();
Mass mwr = Mass.FromGrams(14.5);
Temperature t = Temperature.FromKelvins(6000.0);
Acceleration a = Acceleration.FromStandardGravity(0.9);
Length r = Length.FromKilometers(6500.0);
Assert.AreEqual(Env.GasLife(mwr.Grams,
t.Kelvins,
a.StandardGravity,
r.Kilometers),
_phy.Physics.GetGasLife(mwr,
t,
a,
r).Years365,
LowDelta);
}
public void GetChangeInAngularVelocity()
{
Density densityGCC = Density.FromGramsPerCubicCentimeter(2.5);
Mass massSM = Mass.FromEarthMasses(0.75);
Length radiusKM = Length.FromKilometers(2500.0);
Length semiMajorAxisAU = Length.FromAstronomicalUnits(1.4);
Mass largeMassSM = Mass.FromSolarMasses(1.5);
Assert.AreEqual(Env.ChangeInAngularVelocity(densityGCC.GramsPerCubicCentimeter,
massSM.SolarMasses,
radiusKM.Kilometers,
semiMajorAxisAU.AstronomicalUnits,
largeMassSM.SolarMasses),
_phy.Dynamics.GetChangeInAngularVelocity(densityGCC,
massSM,
radiusKM,
semiMajorAxisAU,
largeMassSM).RadiansPerSecond,
LowDelta);
}