vulcan::gravity::j2j4 Namespace Reference

Functions
template<typename Scalar>
Vec3< Scalar >	acceleration (const Vec3< Scalar > &r_ecef, double mu=constants::earth::mu, double J2_coeff=constants::earth::J2, double J3_coeff=constants::earth::J3, double J4_coeff=constants::earth::J4, double R_eq=constants::earth::R_eq)
	J2/J3/J4 gravitational acceleration.
template<typename Scalar>
Scalar	potential (const Vec3< Scalar > &r_ecef, double mu=constants::earth::mu, double J2_coeff=constants::earth::J2, double J3_coeff=constants::earth::J3, double J4_coeff=constants::earth::J4, double R_eq=constants::earth::R_eq)
	J2/J3/J4 gravitational potential.

Function Documentation

acceleration()

template<typename Scalar>

Vec3< Scalar > vulcan::gravity::j2j4::acceleration	(	const Vec3< Scalar > &	r_ecef,
		double	mu = constants::earth::mu,
		double	J2_coeff = constants::earth::J2,
		double	J3_coeff = constants::earth::J3,
		double	J4_coeff = constants::earth::J4,
		double	R_eq = constants::earth::R_eq )

J2/J3/J4 gravitational acceleration.

Includes first three zonal harmonics for high-fidelity gravity modeling.

• J2 captures Earth's equatorial bulge (oblateness)
• J3 captures the slight pear shape of Earth (north/south asymmetry)
• J4 captures additional higher-order oblateness

Template Parameters

Scalar

double or casadi::MX

Parameters

r_ecef	Position in ECEF frame [m]
mu	Gravitational parameter [m³/s²]
J2_coeff	J2 zonal harmonic coefficient
J3_coeff	J3 zonal harmonic coefficient
J4_coeff	J4 zonal harmonic coefficient
R_eq	Equatorial radius [m]

Returns

Gravitational acceleration in ECEF [m/s²]

potential()

template<typename Scalar>

Scalar vulcan::gravity::j2j4::potential	(	const Vec3< Scalar > &	r_ecef,
		double	mu = constants::earth::mu,
		double	J2_coeff = constants::earth::J2,
		double	J3_coeff = constants::earth::J3,
		double	J4_coeff = constants::earth::J4,
		double	R_eq = constants::earth::R_eq )

J2/J3/J4 gravitational potential.

U = -μ/r · [1 - Σ Jn·(R_eq/r)^n·Pn(sin φ)]

Template Parameters

Scalar

double or casadi::MX

Parameters

r_ecef	Position in ECEF frame [m]
mu	Gravitational parameter [m³/s²]
J2_coeff	J2 zonal harmonic coefficient
J3_coeff	J3 zonal harmonic coefficient
J4_coeff	J4 zonal harmonic coefficient
R_eq	Equatorial radius [m]

Returns

Gravitational potential [m²/s²]