# ESA PB05-Style Brouwer–Lyddane Averaged-Element Orbit Propagator (J2–J6 + Luni-solar)

Propagate near-Earth and GNSS-region orbits with a semi-analytical averaged-element model in the Brouwer–Lyddane / PB05 tradition: orbit-averaged secular drift of the mean elements under the EGM96 zonal field J2–J6 (numerical Gauss-equation averaging — the Delaunay-normalisation analogue with the closed-form J2 rates as a built-in cross-check), first-order J2 short-period mean↔osculating reconstruction, long-period J3 frozen-orbit analysis (e = −J3·sin i/(2·J2·p), ω = 90°), double-averaged luni-solar quadrupole secular effects on an elliptical inclined perturber orbit (the Celletti/Kozai–Lidov term with the 39.23° critical mutual inclination), tesseral ground-track resonance checks (1:1, 2:1, 3:2, 5:2 and small-rational commensurabilities with Earth rotation), sun-synchronous inclination solving, and an in-report RK4 numerical-propagation cross-validation of the averaged model.

> Canonical page: https://elysiatools.com/en/tools/esa-pb05-pce-pck-reddy-porter-keplerian-equation-averaged-element-propagator

- **Category:** Astronomy

- **Keywords:** averaged element propagator, brouwer lyddane mean elements, j2 secular rates, raan drift, frozen orbit, kozai lidov, lunisolar perturbation, sun synchronous orbit, tesseral resonance, mean to osculating

## Overview

The propagator integrates the MEAN elements (a, e, i, Ω, ω, M) with an RK4 step whose right-hand side is the orbit-average of the Gauss variational rates over the zonal field to the selected depth — the numerical equivalent of the Brouwer–Lyddane/PB05 Delaunay normalisation, which the tool cross-checks against the closed-form first-order J2 rates (Ω̇ = −3/2·J2·n·(Re/p)²·cos i, ω̇ = 3/4·J2·n·(Re/p)²·(5cos²i−1)). The luni-solar option adds the quadrupole-order double-averaged third-body term (averaged over both the satellite and the perturber orbit), which reproduces the classic GEO inclination growth of roughly 0.75–0.95°/yr and carries the Kozai–Lidov critical mutual inclination arccos(√(3/5)) = 39.23°. Short-period J2 deviations are integrated from the instantaneous rates and can be added back to sample the osculating elements (mean→osculating back-transformation). The report also runs an independent RK4 Cartesian numerical propagation with the same force model and shows the drift between the two, the tesseral ground-track resonance classification, frozen-orbit (J3 long-period) parameters, and the sun-synchronous inclination solution for the given geometry.

## Inputs

- **Semi-major axis a (km)** (number): 8000
- **Eccentricity e** (number): 0.05
- **Inclination i (deg)** (number): 52.5
- **RAAN Ω (deg)** (number): 30
- **Argument of perigee ω (deg)** (number): 70
- **Mean anomaly M (deg)** (number): 0
- **Propagation duration (days)** (number): 30
- **Zonal gravity depth** (select)
- **Third-body perturbations** (select)
- **Add mean→osculating reconstruction panel** (checkbox)

## When to use

- Analyzing long-term orbital drift in RAAN and argument of perigee caused by Earth oblateness and zonal harmonics (J2–J6).
- Predicting multi-month or multi-year inclination growth in geostationary orbits driven by lunar and solar gravitational perturbations.
- Designing frozen orbits or sun-synchronous trajectories and screening for repeating ground-track tesseral resonances.

## How it works

- Accepts classical Keplerian orbital elements (a, e, i, Ω, ω, M) along with propagation duration, zonal gravity depth (J2, J2–J4, or J2–J6), and third-body configuration.
- Performs RK4 integration on the orbit-averaged Gauss variational equations to propagate mean orbital elements over the chosen timespan.
- Calculates short-period J2 variations to reconstruct osculating elements and checks for Kozai–Lidov critical angles, frozen orbit parameters, and tesseral commensurabilities.
- Generates an interactive HTML report containing mean-rate diagnostic cards, element-history trajectory charts, closed-form J2 rate comparisons, and an independent RK4 Cartesian cross-validation table.

## Use cases

- Mission planning for LEO constellations requiring precise nodal precession rates and ground-track resonance identification.
- Station-keeping budget estimation for geostationary satellites assessing annual inclination drift caused by the Sun and Moon.
- Designing Earth observation frozen orbits by matching eccentricity to the J3 long-period equilibrium condition.

## Frequently asked questions

### What is the primary difference between mean and osculating orbital elements?

Mean elements average out short-period periodic fluctuations to capture secular and long-period drift, whereas osculating elements describe the instantaneous Keplerian state.

### How does the tool handle zonal gravity depth?

You can select J2 only, J2–J4, or J2–J6 based on the EGM96 gravity model to evaluate secular nodal and perigee precession.

### What does the double-averaged luni-solar perturbation model compute?

It averages third-body gravitational forces over both the satellite orbit and the perturber orbit to model secular inclination evolution and Kozai–Lidov dynamics.

### How is the averaged propagator cross-validated?

The generated report includes an independent Cartesian RK4 numerical propagation using the same force model to quantify drift between semi-analytical and numerical methods.

### Can this tool calculate sun-synchronous inclinations?

Yes, the diagnostics panel evaluates the geometry and determines the theoretical sun-synchronous inclination matching the orbit's semi-major axis and eccentricity.

## Related tools

- [Correlation Analyzer](https://elysiatools.com/en/tools/correlation-analyzer): Advanced correlation analysis tool that calculates correlation coefficients between variables to measure the strength and direction of their linear relationships. Perfect for statistical analysis, financial modeling, scientific research, and data exploration. Features: - Multiple correlation methods (Pearson, Spearman, Kendall) - Correlation matrix generation - Statistical significance testing (p-values) - Confidence intervals calculation - Heatmap visualization - Scatter plot matrix generation - Missing value handling strategies - Outlier detection and handling - Group analysis capabilities - Detailed statistical reports Common Use Cases: - Financial market analysis and risk assessment - Scientific research and hypothesis testing - Customer behavior and marketing analysis - Healthcare and medical data analysis - Quality control and process optimization - Educational performance evaluation
- [Image Palette to Design Tokens](https://elysiatools.com/en/tools/image-to-design-tokens): Extract a dominant color palette from an image (k-means clustering), then export it as CSS variables, SCSS variables, a Tailwind config, or JSON design tokens — with named colors and an auto-generated shade scale for each
- [Rent vs Buy Total Cost of Ownership Tutor — Mortgage, Tax, Insurance, Maintenance + NPV Break-Even and IRR](https://elysiatools.com/en/tools/rent-vs-buy-total-cost-of-ownership-mortgage-tax-insurance-maintenance-irr-tutor): Side-by-side monthly rent-vs-buy simulation: mortgage amortization, property tax, insurance, HOA and maintenance against rent plus renter insurance; the renter invests the down payment and monthly savings at an opportunity-cost rate. Outputs break-even year on net sale proceeds, per-year differential NPV, SALT-cap-aware after-tax cash flows, a differential IRR, and a ±20% tornado sensitivity chart.
- [Minor (Local) Head Loss Calculator (h = K·v²/2g)](https://elysiatools.com/en/tools/minor-loss-calculator): Compute the minor (local) head loss and pressure drop for pipe fittings, valves and bends: h_m = K·v²/(2g), ΔP_m = K·ρ·v²/2. Ships typical K values for 14 common components (gate/globe/check/angle/ball valves, 90°/45° elbows, tees, sudden contraction/expansion, sharp entrance/exit); choose a preset or supply a custom K. Supports n identical fittings in series. Returns head loss (m, ft) and pressure drop (Pa, kPa, bar). Density in kg/m³, g/cm³ or lb/ft³.
- [Pipe Pressure Drop Calculator (Darcy-Weisbach)](https://elysiatools.com/en/tools/pipe-pressure-drop-darcy): Compute the Darcy-Weisbach major (friction) pressure drop in a straight pipe: ΔP = f·(L/D)·(ρ·v²/2) Pa and head loss h_f = f·(L/D)·v²/(2g) m. The user supplies the Darcy friction factor f (not the Fanning factor). Length in m/km/ft, diameter in m/cm/mm/inch, density in kg/m³/g/cm³/lb/ft³ — all normalised to SI internally. Returns ΔP in Pa, kPa and bar, and head loss in m and ft. Gravity defaults to 9.81 m/s² and may be overridden.
- [Data Boundary Processor](https://elysiatools.com/en/tools/data-boundary-processor): Advanced boundary value processing tool that identifies and handles minimum/maximum values in numerical data. Perfect for data validation, range checking, statistical analysis, and data preprocessing. Features: - Multiple boundary detection methods (absolute, percentile, standard deviation) - Flexible handling strategies (clip, remove, replace, transform) - Custom range validation - Asymmetric boundary handling - Batch processing capabilities - Comprehensive boundary statistics - Data quality assessment - Visual boundary reports Common Use Cases: - Data validation and quality control - Sensor data range checking - Financial data limit enforcement - Statistical data preprocessing - Machine learning feature engineering - Database constraint validation
- [Utility Bill Line-Item Decoder & Hidden Fee Detector (CSR/CFL/MHC/HFA Riders, Peak Demand, kWh↔Btu↔therm↔ton-hours)](https://elysiatools.com/en/tools/csr-cfl-cfl-mhc-hfa-btu-ton-hours-utility-bill-line-item-validator-and-rgp-peak-explainer): Decode home and commercial electric/gas/water bills: classify tariff riders and charges (CSR/CFL/MHC/HFA-style codes, TOU on/off-peak, demand, capacity tag, decoupling), convert kWh↔Btu↔therms↔Ccf↔m³↔ton-hours with EIA factors, expose hidden fixed fees, and project annual cost with load-factor and peak-shaving savings.
- [Density Calculator (ρ = m/V)](https://elysiatools.com/en/tools/density-calculator): Compute density, mass, or volume from the other two, with specific gravity and a float/sink check

## Samples

- [Web Image Processing Python Samples](https://elysiatools.com/en/samples/web-image-processing-python): Web Python image processing examples using PIL/Pillow including reading, saving, resizing, and format conversion
- [Android Image Processing Java Samples](https://elysiatools.com/en/samples/android-image-processing-java): Android Java image processing examples including reading/saving images, scaling, and format conversion
- [Android Image Processing Kotlin Samples](https://elysiatools.com/en/samples/android-image-processing-kotlin): Android Kotlin image processing examples including reading/saving images, scaling, and format conversion
- [macOS Image Processing Objective-C Samples](https://elysiatools.com/en/samples/macos-image-processing-objectivec): macOS Objective-C image processing examples including image reading/saving, image scaling, and format conversion