public static async Task Aggregator(
FPGA.InputSignal <bool> RXD,
FPGA.OutputSignal <bool> TXD
)
{
Sequential handler = () =>
{
FPU.FPUScopeNoSync();
while (true)
{
FPGA.Optimizations.AddOptimizer <DefaultOptimizer>();
const uint baud = 115200;
float[] buff = new float[4];
for (byte i = 0; i < 4; i++)
{
UART.ReadFloat(baud, RXD, out float tmp);
buff[i] = tmp;
}
var deltaV = FPGAOrbitalCalc.DeltaVInclinationOrbit(
buff[0],
buff[1],
buff[2],
buff[3]
);
UART.WriteFloat(baud, deltaV, TXD);
}
};
FPGA.Config.OnStartup(handler);
}
public static async Task Aggregator(
FPGA.InputSignal <bool> RXD,
FPGA.OutputSignal <bool> TXD
)
{
Sequential handler = () =>
{
FPU.FPUScopeNoSync();
while (true)
{
FPGA.Optimizations.AddOptimizer <DefaultOptimizer>();
const uint baud = 115200;
UART.ReadFloat(baud, RXD, out float mass);
UART.ReadFloat(baud, RXD, out float radius);
var vOrbit = FPGAOrbitalCalc.TOrbit(mass, radius);
UART.WriteFloat(baud, vOrbit, TXD);
}
};
FPGA.Config.OnStartup(handler);
}
public void Quokka_VOrbit()
{
using (var port = new QuokkaPort())
{
var r0 = 6371e+3f;
// orbit velocities around the earth
var cases = new[]
{
// low orbit
new
{
R = r0 + 200e+3f,
M = 5.972e+24f,
V = 7800f
},
// ISS orbit
new
{
R = r0 + 409e+3f,
M = 5.972e+24f,
V = 7660f
},
// geostationary orbit
new
{
R = r0 + 35786e+3f,
M = 5.972e+24f,
V = 3100f
},
// moon orbit
new
{
R = r0 + 399000e+3f,
M = 5.972e+24f,
V = 991f
},
};
foreach (var testCase in cases)
{
float calculated = FPGAOrbitalCalc.VOrbit(testCase.M, testCase.R);
var err = Math.Abs((calculated - testCase.V) / testCase.V);
Assert.IsTrue(err < 0.01);
port.WriteFloat(testCase.M);
port.WriteFloat(testCase.R);
var actualBytes = port.Read(4, true, port.DefaultTimeout);
var actual = TestConverters.FloatFromByteArray(actualBytes);
Assert.AreEqual(calculated, actual, $"Failed for {testCase}");
}
}
}
public void VInclination()
{
var R = 6371e+3f;
var M = 5.972e+24f;
var Inc = 0.4974188f;
var nonOptimized = FPGAOrbitalCalc.DeltaVInclinationOrbit(M, R + 100e+3f, R + 1000e3f, Inc);
var optimized = FPGAOrbitalCalc.DeltaVInclinationOrbitOptimized(M, R + 100e+3f, R + 1000e3f, Inc);
Assert.AreEqual(nonOptimized, optimized);
}
public void Quokka_DeltaVOuterOrbit()
{
using (var port = new QuokkaPort())
{
var r0 = 6371e+3f;
// http://www.satsig.net/orbit-research/delta-v-geo-injection-calculator.htm
var cases = new[]
{
// low orbit
new
{
R1 = r0 + 176e+3f,
R2 = r0 + 350e+3f,
M = 5.972e+24f,
DV = 50.58f
},
// to ISS
new
{
R1 = r0 + 100e+3f,
R2 = r0 + 410e+3f,
M = 5.972e+24f,
DV = 90.06f
},
// low orbit to GSO
new
{
R1 = r0 + 200e+3f,
R2 = r0 + 35786e+3f,
M = 5.972e+24f,
DV = 1477.27f
},
};
foreach (var testCase in cases)
{
float calculated = FPGAOrbitalCalc.DeltaVOuterOrbit(testCase.M, testCase.R1, testCase.R2);
var err = Math.Abs((calculated - testCase.DV) / testCase.DV);
Assert.IsTrue(err < 0.01);
port.WriteFloat(testCase.M);
port.WriteFloat(testCase.R1);
port.WriteFloat(testCase.R2);
var actualBytes = port.Read(4, true, port.DefaultTimeout);
var actual = TestConverters.FloatFromByteArray(actualBytes);
Assert.AreEqual(calculated, actual, $"Failed for {testCase}");
}
}
}
public void Quokka_VEsc()
{
var rnd = new Random(Environment.TickCount);
using (var port = new QuokkaPort())
{
// escape velocities from the surface of objects
var cases = new[]
{
// sun
new
{
// radius, in meters
R = 695510e+3f,
// mass, in kgs
M = 1.989e+30f,
// VEsc, is m/s
V = 617500f
},
// earth
new
{
R = 6371e+3f,
M = 5.972e+24f,
V = 11186f
},
// mars
new
{
R = 3389.5e+3f,
M = 6.39e+23f,
V = 5030f
},
};
foreach (var testCase in cases)
{
float calculated = FPGAOrbitalCalc.VEsc(testCase.M, testCase.R);
var err = Math.Abs((calculated - testCase.V) / testCase.V);
Assert.IsTrue(err < 0.01);
port.WriteFloat(testCase.M);
port.WriteFloat(testCase.R);
var actualBytes = port.Read(4, true, port.DefaultTimeout);
var actual = TestConverters.FloatFromByteArray(actualBytes);
Assert.AreEqual(calculated, actual, $"Failed for {testCase}");
}
}
}
public void Quokka_DeltaVInclinationOrbit()
{
using (var port = new QuokkaPort())
{
var r0 = 6371e+3f;
// http://www.satsig.net/orbit-research/delta-v-geo-injection-calculator.htm
var cases = new[]
{
// low orbit to GSO from Kennedy
new
{
R1 = r0 + 200e+3f,
R2 = r0 + 35786e+3f,
M = 5.972e+24f,
Inc = 0.4974188f, // ~28.5*
DV = 1836.49f
},
};
foreach (var testCase in cases)
{
float calculated = FPGAOrbitalCalc.DeltaVInclinationOrbit(testCase.M, testCase.R1, testCase.R2, testCase.Inc);
var err = Math.Abs((calculated - testCase.DV) / testCase.DV);
Assert.IsTrue(err < 0.01);
port.WriteFloat(testCase.M);
port.WriteFloat(testCase.R1);
port.WriteFloat(testCase.R2);
port.WriteFloat(testCase.Inc);
var actualBytes = port.Read(4, true, port.DefaultTimeout);
var actual = TestConverters.FloatFromByteArray(actualBytes);
Assert.AreEqual(calculated, actual, $"Failed for {testCase}");
}
}
}
