ootk 3.0.0 → 4.0.0-0

package/lib/force/ForceModel.d.ts

@@ -0,0 +1,38 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { J2000, Vector, Vector3D } from 'ootk-core';
+import { Thrust } from './Thrust';
+export declare class ForceModel {
+    private _centralGravity?;
+    private _thirdBodyGravity?;
+    private _solarRadiationPressure?;
+    private _atmosphericDrag?;
+    private _maneuverThrust?;
+    setGravity(mu?: number): this;
+    setEarthGravity(degree: number, order: number): void;
+    setThirdBodyGravity({ moon, sun }: {
+        moon?: boolean;
+        sun?: boolean;
+    }): void;
+    setSolarRadiationPressure(mass: number, area: number, coeff?: number): void;
+    setAtmosphericDrag(mass: number, area: number, coeff?: number, cosine?: number): void;
+    loadManeuver(maneuver: Thrust): void;
+    clearManeuver(): void;
+    acceleration(state: J2000): Vector3D;
+    derivative(state: J2000): Vector;
+}

package/lib/force/ForceModel.js

@@ -0,0 +1,71 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Earth, Vector3D } from 'ootk-core';
+import { AtmosphericDrag } from './AtmosphericDrag';
+import { EarthGravity } from './EarthGravity';
+import { Gravity } from './Gravity';
+import { SolarRadiationPressure } from './SolarRadiationPressure';
+import { ThirdBodyGravity } from './ThirdBodyGravity';
+// / Force model for spacecraft propagation.
+export class ForceModel {
+    setGravity(mu = Earth.mu) {
+        this._centralGravity = new Gravity(mu);
+        return this;
+    }
+    setEarthGravity(degree, order) {
+        this._centralGravity = new EarthGravity(degree, order);
+    }
+    setThirdBodyGravity({ moon = false, sun = false }) {
+        this._thirdBodyGravity = new ThirdBodyGravity(moon, sun);
+    }
+    setSolarRadiationPressure(mass, area, coeff = 1.2) {
+        this._solarRadiationPressure = new SolarRadiationPressure(mass, area, coeff);
+    }
+    setAtmosphericDrag(mass, area, coeff = 2.2, cosine = 4) {
+        this._atmosphericDrag = new AtmosphericDrag(mass, area, coeff, cosine);
+    }
+    loadManeuver(maneuver) {
+        this._maneuverThrust = maneuver;
+    }
+    clearManeuver() {
+        this._maneuverThrust = null;
+    }
+    acceleration(state) {
+        let accVec = Vector3D.origin;
+        if (this._centralGravity) {
+            accVec = accVec.add(this._centralGravity.acceleration(state));
+        }
+        if (this._thirdBodyGravity) {
+            accVec = accVec.add(this._thirdBodyGravity.acceleration(state));
+        }
+        if (this._solarRadiationPressure) {
+            accVec = accVec.add(this._solarRadiationPressure.acceleration(state));
+        }
+        if (this._atmosphericDrag) {
+            accVec = accVec.add(this._atmosphericDrag.acceleration(state));
+        }
+        if (this._maneuverThrust) {
+            accVec = accVec.add(this._maneuverThrust.acceleration(state));
+        }
+        return accVec;
+    }
+    derivative(state) {
+        return state.velocity.join(this.acceleration(state));
+    }
+}
+//# sourceMappingURL=ForceModel.js.map

package/lib/force/ForceModel.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ForceModel.js","sourceRoot":"","sources":["../../src/force/ForceModel.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAiB,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC3D,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AACpD,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,OAAO,EAAE,OAAO,EAAE,MAAM,WAAW,CAAC;AACpC,OAAO,EAAE,sBAAsB,EAAE,MAAM,0BAA0B,CAAC;AAClE,OAAO,EAAE,gBAAgB,EAAE,MAAM,oBAAoB,CAAC;AAGtD,4CAA4C;AAE5C,MAAM,OAAO,UAAU;IAOrB,UAAU,CAAC,KAAa,KAAK,CAAC,EAAE;QAC9B,IAAI,CAAC,eAAe,GAAG,IAAI,OAAO,CAAC,EAAE,CAAC,CAAC;QAEvC,OAAO,IAAI,CAAC;IACd,CAAC;IAED,eAAe,CAAC,MAAc,EAAE,KAAa;QAC3C,IAAI,CAAC,eAAe,GAAG,IAAI,YAAY,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;IACzD,CAAC;IAED,mBAAmB,CAAC,EAAE,IAAI,GAAG,KAAK,EAAE,GAAG,GAAG,KAAK,EAAE;QAC/C,IAAI,CAAC,iBAAiB,GAAG,IAAI,gBAAgB,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;IAC3D,CAAC;IAED,yBAAyB,CAAC,IAAY,EAAE,IAAY,EAAE,KAAK,GAAG,GAAG;QAC/D,IAAI,CAAC,uBAAuB,GAAG,IAAI,sBAAsB,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,CAAC,CAAC;IAC/E,CAAC;IAED,kBAAkB,CAAC,IAAY,EAAE,IAAY,EAAE,KAAK,GAAG,GAAG,EAAE,MAAM,GAAG,CAAC;QACpE,IAAI,CAAC,gBAAgB,GAAG,IAAI,eAAe,CAAC,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,MAAM,CAAC,CAAC;IACzE,CAAC;IAED,YAAY,CAAC,QAAgB;QAC3B,IAAI,CAAC,eAAe,GAAG,QAAQ,CAAC;IAClC,CAAC;IAED,aAAa;QACX,IAAI,CAAC,eAAe,GAAG,IAAI,CAAC;IAC9B,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,IAAI,MAAM,GAAG,QAAQ,CAAC,MAAM,CAAC;QAE7B,IAAI,IAAI,CAAC,eAAe,EAAE;YACxB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,eAAe,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SAC/D;QACD,IAAI,IAAI,CAAC,iBAAiB,EAAE;YAC1B,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,iBAAiB,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SACjE;QACD,IAAI,IAAI,CAAC,uBAAuB,EAAE;YAChC,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,uBAAuB,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SACvE;QACD,IAAI,IAAI,CAAC,gBAAgB,EAAE;YACzB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,gBAAgB,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SAChE;QACD,IAAI,IAAI,CAAC,eAAe,EAAE;YACxB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,eAAe,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SAC/D;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,UAAU,CAAC,KAAY;QACrB,OAAO,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;IACvD,CAAC;CACF"}

package/lib/force/Gravity.d.ts

@@ -0,0 +1,41 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class Gravity implements Force {
+    mu: number;
+    /**
+     * Create a new [Gravity] object with optional gravitational
+     * @param mu Gravitational parameter _(km²/s³)_.
+     */
+    constructor(mu?: number);
+    /**
+     * Calculates the gravitational force in spherical coordinates.
+     *
+     * @param state The J2000 state containing the position and velocity vectors.
+     * @returns The gravitational force vector in spherical coordinates.
+     */
+    private _spherical;
+    /**
+     * Calculates the acceleration due to gravity at a given state.
+     *
+     * @param state The J2000 state at which to calculate the acceleration.
+     * @returns The acceleration vector.
+     */
+    acceleration(state: J2000): Vector3D;
+}

package/lib/force/Gravity.js

@@ -0,0 +1,48 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Earth } from 'ootk-core';
+// / Simple central-body gravity model.
+export class Gravity {
+    /**
+     * Create a new [Gravity] object with optional gravitational
+     * @param mu Gravitational parameter _(km²/s³)_.
+     */
+    constructor(mu = Earth.mu) {
+        this.mu = mu;
+    }
+    /**
+     * Calculates the gravitational force in spherical coordinates.
+     *
+     * @param state The J2000 state containing the position and velocity vectors.
+     * @returns The gravitational force vector in spherical coordinates.
+     */
+    _spherical(state) {
+        const rMag = state.position.magnitude();
+        return state.position.scale(-this.mu / (rMag * rMag * rMag));
+    }
+    /**
+     * Calculates the acceleration due to gravity at a given state.
+     *
+     * @param state The J2000 state at which to calculate the acceleration.
+     * @returns The acceleration vector.
+     */
+    acceleration(state) {
+        return this._spherical(state);
+    }
+}
+//# sourceMappingURL=Gravity.js.map

package/lib/force/Gravity.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Gravity.js","sourceRoot":"","sources":["../../src/force/Gravity.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAmB,MAAM,WAAW,CAAC;AAGnD,uCAAuC;AACvC,MAAM,OAAO,OAAO;IAIlB;;;OAGG;IACH,YAAY,KAAa,KAAK,CAAC,EAAE;QAC/B,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;IACf,CAAC;IAED;;;;;OAKG;IACK,UAAU,CAAC,KAAY;QAC7B,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAExC,OAAO,KAAK,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC;IAC/D,CAAC;IAED;;;;;OAKG;IACH,YAAY,CAAC,KAAY;QACvB,OAAO,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;IAChC,CAAC;CACF"}

package/lib/force/SolarRadiationPressure.d.ts

@@ -0,0 +1,27 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class SolarRadiationPressure extends Force {
+    mass: number;
+    area: number;
+    reflectCoeff: number;
+    constructor(mass: number, area: number, reflectCoeff: number);
+    private static readonly _kRef;
+    acceleration(state: J2000): Vector3D;
+}

package/lib/force/SolarRadiationPressure.js

@@ -0,0 +1,42 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { astronomicalUnit, Sun } from 'ootk-core';
+import { Force } from './Force';
+// / Solar radiation pressure model.
+export class SolarRadiationPressure extends Force {
+    // / Create a new [SolarRadiationPressure] object.
+    constructor(mass, area, reflectCoeff) {
+        super();
+        this.mass = mass;
+        this.area = area;
+        this.reflectCoeff = reflectCoeff;
+    }
+    acceleration(state) {
+        const rSun = Sun.positionApparent(state.epoch);
+        const r = state.position.subtract(rSun);
+        const rMag = r.magnitude();
+        const r2 = rMag * rMag;
+        const ratio = Sun.lightingRatio(state.position, rSun);
+        const p = (ratio * SolarRadiationPressure._kRef) / r2;
+        const flux = r.scale(p / rMag);
+        return flux.scale(((this.area * this.reflectCoeff) / this.mass) * 1e-3);
+    }
+}
+// / Solar pressure _(N/m²)_;
+SolarRadiationPressure._kRef = 4.56e-6 * astronomicalUnit ** 2;
+//# sourceMappingURL=SolarRadiationPressure.js.map

package/lib/force/SolarRadiationPressure.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SolarRadiationPressure.js","sourceRoot":"","sources":["../../src/force/SolarRadiationPressure.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,gBAAgB,EAAS,GAAG,EAAY,MAAM,WAAW,CAAC;AACnE,OAAO,EAAE,KAAK,EAAE,MAAM,SAAS,CAAC;AAEhC,oCAAoC;AACpC,MAAM,OAAO,sBAAuB,SAAQ,KAAK;IAC/C,kDAAkD;IAClD,YAAmB,IAAY,EAAS,IAAY,EAAS,YAAoB;QAC/E,KAAK,EAAE,CAAC;QADS,SAAI,GAAJ,IAAI,CAAQ;QAAS,SAAI,GAAJ,IAAI,CAAQ;QAAS,iBAAY,GAAZ,YAAY,CAAQ;IAEjF,CAAC;IAKD,YAAY,CAAC,KAAY;QACvB,MAAM,IAAI,GAAG,GAAG,CAAC,gBAAgB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QAC/C,MAAM,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;QACxC,MAAM,IAAI,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC;QACvB,MAAM,KAAK,GAAG,GAAG,CAAC,aAAa,CAAC,KAAK,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC;QACtD,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,sBAAsB,CAAC,KAAK,CAAC,GAAG,EAAE,CAAC;QACtD,MAAM,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC;QAE/B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,YAAY,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC;IAC1E,CAAC;;AAbD,6BAA6B;AACL,4BAAK,GAAW,OAAO,GAAG,gBAAgB,IAAI,CAAC,CAAC"}

package/lib/force/ThirdBodyGravity.d.ts

@@ -0,0 +1,27 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class ThirdBodyGravity implements Force {
+    moon: boolean;
+    sun: boolean;
+    constructor(moon?: boolean, sun?: boolean);
+    private static _moonGravity;
+    private static _sunGravity;
+    acceleration(state: J2000): Vector3D;
+}

package/lib/force/ThirdBodyGravity.js

@@ -0,0 +1,56 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Moon, Sun, Vector3D } from 'ootk-core';
+// / Third-body gravity model.
+export class ThirdBodyGravity {
+    // / Create a new [ThirdBodyGravity] object with the selected bodies enabled.
+    constructor(moon = false, sun = false) {
+        this.moon = moon;
+        this.sun = sun;
+        // Nothing to do here.
+    }
+    static _moonGravity(state) {
+        const rMoon = Moon.position(state.epoch);
+        const aNum = rMoon.subtract(state.position);
+        const aDen = aNum.magnitude() ** 3;
+        const bNum = rMoon;
+        const bDen = rMoon.magnitude() ** 3;
+        const gravity = aNum.scale(1 / aDen).add(bNum.scale(-1 / bDen));
+        return gravity.scale(Moon.mu);
+    }
+    static _sunGravity(state) {
+        const rSun = Sun.positionApparent(state.epoch);
+        const aNum = rSun.subtract(state.position);
+        const aDen = aNum.magnitude() ** 3;
+        const bNum = rSun;
+        const bDen = rSun.magnitude() ** 3;
+        const gravity = aNum.scale(1 / aDen).add(bNum.scale(-1 / bDen));
+        return gravity.scale(Sun.mu);
+    }
+    acceleration(state) {
+        let accVec = Vector3D.origin;
+        if (this.moon) {
+            accVec = accVec.add(ThirdBodyGravity._moonGravity(state));
+        }
+        if (this.sun) {
+            accVec = accVec.add(ThirdBodyGravity._sunGravity(state));
+        }
+        return accVec;
+    }
+}
+//# sourceMappingURL=ThirdBodyGravity.js.map

package/lib/force/ThirdBodyGravity.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ThirdBodyGravity.js","sourceRoot":"","sources":["../../src/force/ThirdBodyGravity.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAS,IAAI,EAAE,GAAG,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAGvD,8BAA8B;AAC9B,MAAM,OAAO,gBAAgB;IAC3B,6EAA6E;IAC7E,YAAmB,OAAgB,KAAK,EAAS,MAAe,KAAK;QAAlD,SAAI,GAAJ,IAAI,CAAiB;QAAS,QAAG,GAAH,GAAG,CAAiB;QACnE,sBAAsB;IACxB,CAAC;IAEO,MAAM,CAAC,YAAY,CAAC,KAAY;QACtC,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QACzC,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QAC5C,MAAM,IAAI,GAAG,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,KAAK,CAAC;QACnB,MAAM,IAAI,GAAG,KAAK,CAAC,SAAS,EAAE,IAAI,CAAC,CAAC;QACpC,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAEhE,OAAO,OAAO,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAChC,CAAC;IAEO,MAAM,CAAC,WAAW,CAAC,KAAY;QACrC,MAAM,IAAI,GAAG,GAAG,CAAC,gBAAgB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QAC/C,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QAC3C,MAAM,IAAI,GAAG,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,IAAI,CAAC;QAClB,MAAM,IAAI,GAAG,IAAI,CAAC,SAAS,EAAE,IAAI,CAAC,CAAC;QACnC,MAAM,OAAO,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAEhE,OAAO,OAAO,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IAC/B,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,IAAI,MAAM,GAAG,QAAQ,CAAC,MAAM,CAAC;QAE7B,IAAI,IAAI,CAAC,IAAI,EAAE;YACb,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,gBAAgB,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC;SAC3D;QACD,IAAI,IAAI,CAAC,GAAG,EAAE;YACZ,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,gBAAgB,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC;SAC1D;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF"}

package/lib/force/Thrust.d.ts

@@ -0,0 +1,35 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class Thrust implements Force {
+    center: EpochUTC;
+    radial: number;
+    intrack: number;
+    crosstrack: number;
+    durationRate: number;
+    constructor(center: EpochUTC, radial: number, intrack: number, crosstrack: number, durationRate?: number);
+    deltaV: Vector3D;
+    get magnitude(): number;
+    get duration(): number;
+    get start(): EpochUTC;
+    get stop(): EpochUTC;
+    acceleration(state: J2000): Vector3D;
+    apply(state: J2000): J2000;
+    get isImpulsive(): boolean;
+}

package/lib/force/Thrust.js

@@ -0,0 +1,53 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { RIC, Vector3D } from 'ootk-core';
+// / Thrust force model.
+export class Thrust {
+    constructor(center, radial, intrack, crosstrack, durationRate = 0.0) {
+        this.center = center;
+        this.radial = radial;
+        this.intrack = intrack;
+        this.crosstrack = crosstrack;
+        this.durationRate = durationRate;
+        this.deltaV = new Vector3D(radial * 1e-3, intrack * 1e-3, crosstrack * 1e-3);
+    }
+    get magnitude() {
+        return this.deltaV.magnitude() * 1000.0;
+    }
+    get duration() {
+        return this.magnitude * this.durationRate;
+    }
+    get start() {
+        return this.center.roll(-0.5 * this.duration);
+    }
+    get stop() {
+        return this.center.roll(0.5 * this.duration);
+    }
+    acceleration(state) {
+        const relative = new RIC(Vector3D.origin, this.deltaV.scale(1.0 / this.duration));
+        return relative.toJ2000(state).velocity.subtract(state.velocity);
+    }
+    apply(state) {
+        const relative = new RIC(Vector3D.origin, this.deltaV);
+        return relative.toJ2000(state);
+    }
+    get isImpulsive() {
+        return this.duration <= 0;
+    }
+}
+//# sourceMappingURL=Thrust.js.map

package/lib/force/Thrust.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Thrust.js","sourceRoot":"","sources":["../../src/force/Thrust.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAmB,GAAG,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAG3D,wBAAwB;AACxB,MAAM,OAAO,MAAM;IACjB,YACS,MAAgB,EAChB,MAAc,EACd,OAAe,EACf,UAAkB,EAClB,eAAuB,GAAG;QAJ1B,WAAM,GAAN,MAAM,CAAU;QAChB,WAAM,GAAN,MAAM,CAAQ;QACd,YAAO,GAAP,OAAO,CAAQ;QACf,eAAU,GAAV,UAAU,CAAQ;QAClB,iBAAY,GAAZ,YAAY,CAAc;QAEjC,IAAI,CAAC,MAAM,GAAG,IAAI,QAAQ,CAAC,MAAM,GAAG,IAAI,EAAE,OAAO,GAAG,IAAI,EAAE,UAAU,GAAG,IAAI,CAAC,CAAC;IAC/E,CAAC;IAID,IAAI,SAAS;QACX,OAAO,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,GAAG,MAAM,CAAC;IAC1C,CAAC;IAED,IAAI,QAAQ;QACV,OAAO,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,YAAY,CAAC;IAC5C,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC;IAChD,CAAC;IAED,IAAI,IAAI;QACN,OAAO,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC;IAC/C,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,QAAQ,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;QAElF,OAAO,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;IACnE,CAAC;IAED,KAAK,CAAC,KAAY;QAChB,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,QAAQ,CAAC,MAAM,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAEvD,OAAO,QAAQ,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;IACjC,CAAC;IAED,IAAI,WAAW;QACb,OAAO,IAAI,CAAC,QAAQ,IAAI,CAAC,CAAC;IAC5B,CAAC;CACF"}

package/lib/interpolator/ChebyshevCoefficients.d.ts

@@ -0,0 +1,54 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { PositionVelocity } from 'ootk-core';
+export declare class ChebyshevCoefficients {
+    a: number;
+    b: number;
+    private _cx;
+    private _cy;
+    private _cz;
+    _cxd: Float64Array;
+    _cyd: Float64Array;
+    _czd: Float64Array;
+    constructor(a: number, b: number, _cx: Float64Array, _cy: Float64Array, _cz: Float64Array);
+    /**
+     * Calculates the derivative of a polynomial represented by Chebyshev coefficients.
+     *
+     * @param a - The lower bound of the polynomial's domain.
+     * @param b - The upper bound of the polynomial's domain.
+     * @param c - The Chebyshev coefficients of the polynomial.
+     * @returns The derivative of the polynomial as an array of coefficients.
+     */
+    private static _derivative;
+    get sizeBytes(): number;
+    /**
+     * Evaluates the Chebyshev polynomial represented by the given coefficients at the specified value.
+     *
+     * @param c - The coefficients of the Chebyshev polynomial.
+     * @param t - The value at which to evaluate the polynomial _(POSIX seconds)_.
+     * @returns The result of evaluating the polynomial at the specified value.
+     */
+    evaluate(c: Float64Array, t: number): number;
+    /**
+     * Interpolates the position and velocity at a given time (km, km/s).
+     *
+     * @param t - The time value to interpolate at _(POSIX seconds)_.
+     * @returns An object containing the interpolated position and velocity.
+     */
+    interpolate(t: number): PositionVelocity;
+}

package/lib/interpolator/ChebyshevCoefficients.js

@@ -0,0 +1,90 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Vector3D } from 'ootk-core';
+// / Chebyshev compressed ephemeris coefficients.
+export class ChebyshevCoefficients {
+    // / Create a new [ChebyshevCoefficients] object.
+    constructor(a, b, _cx, _cy, _cz) {
+        this.a = a;
+        this.b = b;
+        this._cx = _cx;
+        this._cy = _cy;
+        this._cz = _cz;
+        this._cxd = ChebyshevCoefficients._derivative(a, b, _cx);
+        this._cyd = ChebyshevCoefficients._derivative(a, b, _cy);
+        this._czd = ChebyshevCoefficients._derivative(a, b, _cz);
+    }
+    /**
+     * Calculates the derivative of a polynomial represented by Chebyshev coefficients.
+     *
+     * @param a - The lower bound of the polynomial's domain.
+     * @param b - The upper bound of the polynomial's domain.
+     * @param c - The Chebyshev coefficients of the polynomial.
+     * @returns The derivative of the polynomial as an array of coefficients.
+     */
+    static _derivative(a, b, c) {
+        const n = c.length;
+        const d = new Float64Array(n);
+        d[n - 1] = 0;
+        d[n - 2] = 2 * (n - 1) * c[n - 1];
+        for (let k = n - 3; k >= 0; k--) {
+            d[k] = d[k + 2] + 2 * (k + 1) * c[k + 1];
+        }
+        for (let k = 0; k < n; k++) {
+            d[k] *= 2 / (b - a);
+        }
+        return d;
+    }
+    // / Return the size _(bytes)_ of this coefficient set's cached data.
+    get sizeBytes() {
+        return (64 * 2 + 64 * 3 * this._cx.length) / 8;
+    }
+    /**
+     * Evaluates the Chebyshev polynomial represented by the given coefficients at the specified value.
+     *
+     * @param c - The coefficients of the Chebyshev polynomial.
+     * @param t - The value at which to evaluate the polynomial _(POSIX seconds)_.
+     * @returns The result of evaluating the polynomial at the specified value.
+     */
+    evaluate(c, t) {
+        const n = c.length;
+        const x = (t - 0.5 * (this.b + this.a)) / (0.5 * (this.b - this.a));
+        const alpha = 2 * x;
+        const beta = -1;
+        let y1 = 0.0;
+        let y2 = 0.0;
+        for (let k = n - 1; k >= 1; k--) {
+            const tmp = y1;
+            y1 = alpha * y1 + beta * y2 + c[k];
+            y2 = tmp;
+        }
+        return x * y1 - y2 + 0.5 * c[0];
+    }
+    /**
+     * Interpolates the position and velocity at a given time (km, km/s).
+     *
+     * @param t - The time value to interpolate at _(POSIX seconds)_.
+     * @returns An object containing the interpolated position and velocity.
+     */
+    interpolate(t) {
+        const pos = new Vector3D(this.evaluate(this._cx, t), this.evaluate(this._cy, t), this.evaluate(this._cz, t));
+        const vel = new Vector3D(this.evaluate(this._cxd, t), this.evaluate(this._cyd, t), this.evaluate(this._czd, t));
+        return { position: pos, velocity: vel };
+    }
+}
+//# sourceMappingURL=ChebyshevCoefficients.js.map

package/lib/interpolator/ChebyshevCoefficients.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ChebyshevCoefficients.js","sourceRoot":"","sources":["../../src/interpolator/ChebyshevCoefficients.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAoB,QAAQ,EAAE,MAAM,WAAW,CAAC;AAEvD,iDAAiD;AACjD,MAAM,OAAO,qBAAqB;IAIhC,iDAAiD;IACjD,YACS,CAAS,EACT,CAAS,EACR,GAAiB,EACjB,GAAiB,EACjB,GAAiB;QAJlB,MAAC,GAAD,CAAC,CAAQ;QACT,MAAC,GAAD,CAAC,CAAQ;QACR,QAAG,GAAH,GAAG,CAAc;QACjB,QAAG,GAAH,GAAG,CAAc;QACjB,QAAG,GAAH,GAAG,CAAc;QAEzB,IAAI,CAAC,IAAI,GAAG,qBAAqB,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QACzD,IAAI,CAAC,IAAI,GAAG,qBAAqB,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QACzD,IAAI,CAAC,IAAI,GAAG,qBAAqB,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;IAC3D,CAAC;IAED;;;;;;;OAOG;IACK,MAAM,CAAC,WAAW,CAAC,CAAS,EAAE,CAAS,EAAE,CAAe;QAC9D,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAE9B,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACb,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;YAC/B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;SAC1C;QACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;SACrB;QAED,OAAO,CAAC,CAAC;IACX,CAAC;IAED,qEAAqE;IACrE,IAAI,SAAS;QACX,OAAO,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACjD,CAAC;IAED;;;;;;OAMG;IACH,QAAQ,CAAC,CAAe,EAAE,CAAS;QACjC,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnB,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACpE,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;QACpB,MAAM,IAAI,GAAG,CAAC,CAAC,CAAC;QAChB,IAAI,EAAE,GAAG,GAAG,CAAC;QACb,IAAI,EAAE,GAAG,GAAG,CAAC;QAEb,KAAK,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;YAC/B,MAAM,GAAG,GAAG,EAAE,CAAC;YAEf,EAAE,GAAG,KAAK,GAAG,EAAE,GAAG,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACnC,EAAE,GAAG,GAAG,CAAC;SACV;QAED,OAAO,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC;IAED;;;;;OAKG;IACH,WAAW,CAAC,CAAS;QACnB,MAAM,GAAG,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;QAC7G,MAAM,GAAG,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,EAAE,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC;QAEhH,OAAO,EAAE,QAAQ,EAAE,GAAG,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC;IAC1C,CAAC;CACF"}