Model Predictive Path Planning in Navier-Stokes Flow with POD-Based Reduced-Order Models

Basic Information

• Paper ID: 2511.22123
• Title: Model Predictive Path Planning in Navier-Stokes Flow with POD-Based Reduced-Order Models
• Authors: Adam Waterman, Martin Guay (Queen's University, Canada)
• Classification: math.OC (Optimization and Control), cs.SY (Systems and Control), eess.SY (Systems and Control)
• Publication Date: November 27, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2511.22123

Abstract

This paper proposes a framework for optimal trajectory generation in flow-driven systems governed by the Navier-Stokes equations, combining Proper Orthogonal Decomposition (POD)-based reduced-order models (ROM) with Model Predictive Control (MPC). The approach comprises: (i) approximating velocity fields from data through snapshot POD and orthogonal projection; (ii) deriving Galerkin-projected dynamics models in reduced-order coordinates; (iii) employing MPC to plan control inputs that guide an agent through the predicted flow field while satisfying state and actuator constraints. By leveraging reduced-order modeling, the method achieves real-time control in high-dimensional flow environments. Simulations demonstrate accurate flow field reconstruction and efficient trajectory generation in realistic wind environments.

Research Background and Motivation

Problem Definition

This research addresses the real-time trajectory planning problem in flow-driven systems, particularly for high-dimensional flow environments governed by the Navier-Stokes equations. Typical application scenarios include:

• Station-keeping and navigation of stratospheric balloons
• Path planning for underwater gliders
• Other mobile platforms dominated by environmental flow fields with limited actuator capabilities

Significance

1. Computational Challenges: Direct application of high-fidelity Navier-Stokes models for MPC optimization is computationally prohibitive and cannot meet real-time control requirements
2. Practical Demands: Many real systems (e.g., high-altitude balloons) have horizontal motion primarily driven by environmental wind fields, with only limited control available in the vertical direction
3. Predictive Value: Accurate prediction of flow field evolution and exploitation of flow structures (e.g., vertical wind shear) can significantly improve navigation efficiency

Limitations of Existing Approaches

1. Traditional MPC: Direct application to PDE-constrained systems incurs excessive computational costs
2. Reinforcement Learning Methods: Approaches such as Balloon Learning Environment (BLE) lack physical consistency and interpretability
3. Existing ROM-MPC: Primarily focus on closed-loop stabilization rather than path generation in flow-driven scenarios

Research Motivation

This paper proposes a physics-informed, model-based predictive path planning framework that combines POD-based reduced-order modeling with receding horizon optimization to achieve computationally efficient real-time trajectory generation while maintaining physical consistency.

Core Contributions

1. Integrated Framework: First integration of POD-Galerkin reduced-order models with MPC for path planning in Navier-Stokes flow fields
2. Time-Varying Observer Design: Proposes an Extended Kalman Filter (EKF) combining fixed sensor networks and mobile agent measurements, enabling flow state estimation through time-varying measurement matrix C(t)
3. Flow-Aware Trajectory Generation: Develops finite-horizon optimization that leverages ROM-predicted velocity fields, exploiting flow stratification through vertical control inputs
4. Practical Validation: Verifies the method using ERA5 reanalysis wind field data, demonstrating effectiveness and computational efficiency in realistic meteorological environments

Methodology Details

Task Definition

Inputs:

• Historical flow field snapshot data {u_k(x)}^m_
• Fixed sensor network measurements y_fixed(t)
• Agent current position and velocity measurements y_mobile(t)
• Target location x_ref

Outputs:

• Vertical velocity control sequence {u_{z,k+i|k}}^{N-1}_
• Predicted trajectory {x_{k+i|k}}^N_

Constraints:

• Spatial boundaries x ∈ X ⊂ ℝ³
• Actuator limits u_z ∈ U (e.g., maximum climb rate |u_z| ≤ 1 m/s)
• Terminal constraints x_{k+N|k} ∈ X_f

Model Architecture

1. POD Reduced-Order Modeling

Velocity Field Approximation:

v(x,t) = Σ(i=1 to n) a_i(t)φ_i(x)

where φ_i(x) are POD spatial modes and a_i(t) are time-varying modal coefficients.

POD Mode Computation: Obtained by solving the eigenvalue problem:

Uc = λc, U_ij = (1/m)⟨u_i, u_j⟩

Modes are ordered by energy λ_ℓ in descending order, with the first n modes capturing dominant kinetic energy.

Galerkin-Projected Dynamics: Substituting the velocity field expansion into the Navier-Stokes equation and projecting onto the POD mode space:

ȧ_k(t) = L_k a(t) + a(t)^T Q_k a(t)

where:

• L_k: Linear term with elements L_ = ⟨∇²φ_i(x), φ_k(x)⟩
• Q_k: Quadratic nonlinear term with elements Q_ = ⟨(φ_j(x)·∇)φ_i(x), φ_k(x)⟩

2. Extended Kalman Filter Observer

State-Space Model:

ȧ_k(t) = L_k a(t) + a(t)^T Q_k a(t)  (State equation)
y(t) = C(t)a(t)                        (Observation equation)

Time-Varying Measurement Matrix:

C(t) = [C_fixed; C_mobile(t)]

• C_fixed: POD modal evaluations at fixed sensor locations
• C_mobile(t): POD modal evaluations at agent's current location

Discrete-Time EKF:

Prediction step:

â_{k|k-1} = f_d(â_{k-1})
P_{k|k-1} = A_{k-1}P_{k-1|k-1}A^T_{k-1} + Q

Update step:

K_k = P_{k|k-1}C^T_k(C_k P_{k|k-1}C^T_k + R)^{-1}
â_{k|k} = â_{k|k-1} + K_k(y_k - C_k â_{k|k-1})
P_{k|k} = (I - K_k C_k)P_{k|k-1}

3. Model Predictive Path Planning

Agent Kinematics:

ẋ(t) = v_pred(x(t), t) + u_z(t)e_z

Horizontal motion is entirely flow-driven; vertical direction is controllable.

Finite-Horizon Optimization Problem:

min_{u_z} V_f(x_{k+N|k}) + Σ(i=0 to N-1) ℓ(x_{k+i|k}, u_{z,k+i|k})

s.t. x_{k+i+1|k} = x_{k+i|k} + Δt[v_pred(x_{k+i|k}, t_{k+i}) + u_{z,k+i|k}e_z]
     x_{k+i|k} ∈ X, u_{z,k+i|k} ∈ U, x_{k+N|k} ∈ X_f

Cost Function:

ℓ(x, u_z) = w_p ||x - x_ref||²_2 + w_u u²_z

Balances target tracking and control effort.

Technical Innovations

1. Physical Consistency: Galerkin projection preserves the physical structure of the Navier-Stokes equations rather than employing black-box data-driven models
2. Mobile Sensing: The agent itself serves as a mobile sensor, with time-varying measurement matrix C(t) enabling continuous improvement of flow field estimation
3. Flow Field Exploitation: The optimizer actively exploits wind shear by adjusting altitude rather than passively responding to the environment
4. Receding Horizon Strategy: Only the first control action is executed at each step, with subsequent replanning based on updated flow predictions, providing robustness to prediction uncertainty
5. Computational Efficiency: The reduced-order model reduces the high-dimensional PDE to an n-dimensional ODE system (where n is much smaller than the original mesh dimension), enabling real-time optimization

Experimental Setup

Dataset

ERA5 Reanalysis Data:

• Spatial Domain: Centered at 49°N, 81°W, spanning 2 degrees in both longitude and latitude (~146 km × 220 km)
• Vertical Range: Pressure levels from 125 hPa to 10 hPa (~15-30 km altitude)
• Temporal Sampling: 36 randomly selected days over 6 months, uniformly sampled hourly
• Snapshot Count: m velocity field snapshots used for POD decomposition
• Preprocessing: Mean deviation removed to allow ROM to capture wind field dynamics

Evaluation Metrics

1. Final Target Distance:

   d_f = ||x(T_f) - x_ref||_2
   

2. Average Flow Alignment:

   γ̄ = (1/T_f)∫[ẋ(t)·v_pred(x,t)/(||ẋ(t)|| ||v_pred(x,t)||)]dt
   

   
   Measures planner's effectiveness in exploiting favorable flow directions
3. Control Cost:

   u_rms = √[(1/T_f)∫u²_z(t)dt]
   

4. Runtime (RT): Total computation time measured using MATLAB tic-toc

Comparison Methods

Primarily through varying prediction horizon lengths:

• T = 1.5h (short horizon)
• T = 3h (medium-short horizon)
• T = 6h (medium-long horizon)
• T = 12h (long horizon)

Implementation Details

Computing Platform:

• MATLAB R2024a
• AMD Ryzen 2600X (6-core, 3.6 GHz)
• 32 GB RAM

Optimizer:

• Sequential Quadratic Programming (SQP)
• Warm-start strategy (initialization with previous step's input sequence)

Simulation Parameters:

• Total simulation time: 24 hours
• Sampling interval: Δt = 10 minutes
• Target region: Cylindrical region with 50 km radius
• Maximum climb rate: u_max = 1 m/s
• Number of sensors: 8 (including agent)

Cost Weights:

• w_p: Position deviation weight
• w_u: Control cost weight (Empirically tuned to balance performance)

Experimental Results

Main Results

Quantitative Performance Comparison (Table I):

Key Findings:

1. Significant Horizon Effect: Increasing from 1.5h to 6h reduces final distance from 54.6 km to 13.5 km (75% improvement)
2. Diminishing Returns: Improvement from 6h to 12h is only 0.1 km, while computation time increases nearly threefold
3. High Flow Alignment: All configurations achieve γ̄ ≥ 0.92, indicating effective flow exploitation
4. Improved Control Efficiency: Longer horizons enable smoother altitude adjustments (u_rms decreases from 0.60 to 0.36 m/s)

Flow Field Reconstruction Accuracy

Figure 3 shows:

• ROM prediction RMSE relative to actual velocity field over time
• Prediction error remains acceptable at the 15-hour mark
• Heatmaps in Figures 1 and 2 demonstrate high spatial structure consistency between predicted and actual fields

Trajectory Analysis

3-Hour Prediction Horizon (Figures 4 and 6):

• Agent starts in west-east convective region
• Frequent altitude adjustments to find favorable wind layers
• Final distance to target: 32.4 km
• Control actions are aggressive (u_rms = 0.56 m/s)

12-Hour Prediction Horizon (Figures 5 and 7):

• Smoother and more strategic trajectory
• Avoids short-term gains with long-term disadvantages
• Altitude changes are gradual and planned
• Final distance: 13.4 km, control cost reduced by 40%

Experimental Findings

1. Prediction Horizon Trade-offs:
   • Short horizon: Locally feasible but globally inefficient, producing "myopic" behavior
   • Long horizon: Globally optimal but computationally expensive, sensitive to long-term prediction uncertainty
   • Optimal horizon: 6-12 hours balances performance and computation
2. Flow Exploitation Strategy:
   • Planner learns to "wait": moves slowly through unfavorable wind layers, awaiting favorable conditions
   • Active wind shear exploitation: adjusts altitude to enter wind layers favorable to target direction
3. Computational Feasibility:
   • Even with 12-hour horizon, total computation time of 1983 seconds remains acceptable for 24-hour simulation
   • Real applications could achieve further acceleration through parallelization and specialized solvers
4. Robustness:
   • Receding horizon replanning provides natural robustness to flow field prediction errors
   • Continuous EKF updates reduce prediction drift

Related Work

Model Predictive Control

• Process Control: Building climate control 1,2
• Flow Control: Cylinder wake control 3,4
• PDE-Constrained Optimization: Surveys 5-7 demonstrate ROM-MPC applications in large-scale systems

Reduced-Order Modeling

• POD Theory: First introduced by Lumley 11, formalized by Berkooz et al. 12
• Fluid Mechanics Applications: Rowley's 9 balanced POD, Noack et al.'s 10 flow control ROM
• Comprehensive Reviews: Yano 8 and Taira et al. 14 provide comprehensive surveys of aerodynamic and fluid modal analysis

Data-Driven Methods

• DMD and Koopman Operators: Proctor et al. 16, Korda and Mezić 15
• Reinforcement Learning: Balloon Learning Environment 17,18 for stratospheric balloon navigation

Advantages of This Work

1. Physics-Informed: Maintains NS equation physical consistency compared to pure learning methods
2. Real-Time Feasible: Achieves orders of magnitude computational efficiency improvement over high-fidelity CFD-MPC
3. Path Generation Focus: Distinct from traditional ROM-MPC stabilization objectives, emphasizing navigation performance

Conclusions and Discussion

Main Conclusions

1. Framework Effectiveness: The combination of POD-Galerkin ROM and MPC successfully enables efficient path planning in Navier-Stokes flow fields
2. Real-Time Capability: Reduced-order modeling reduces computational complexity from high-dimensional PDE to low-dimensional ODE, making receding horizon optimization feasible on practical timescales
3. Physical Consistency: Galerkin projection preserves fundamental fluid dynamics structures, producing physically reasonable trajectories
4. Practical Validation: Simulations using real ERA5 wind field data demonstrate applicability in complex meteorological environments

Limitations

1. Model Assumptions:
   • Assumes agent as point mass, neglecting actual aircraft attitude dynamics and aerodynamic properties
   • Considers only vertical control, not modeling horizontal propulsion capabilities
2. ROM Limitations:
   • POD modes based on historical data may poorly predict extreme flow fields outside training set
   • Linear subspace assumption may fail to capture highly nonlinear or transient phenomena
3. Computational Cost:
   • Long horizon (12h) computation time remains high (~2000 seconds)
   • Lacks direct computational comparison with high-fidelity CFD or alternative methods
4. Feasibility Guarantees:
   • Cannot guarantee strict recursive feasibility due to flow-driven horizontal motion
   • Lacks theoretical stability or convergence proofs
5. Experimental Scope:
   • Testing limited to single geographic region and season
   • Does not account for measurement noise, sensor failures, and other practical factors

Future Directions

1. Dynamics Integration:
   • Incorporate detailed aircraft dynamics models (attitude, aerodynamic drag)
   • Consider multi-input control (horizontal propulsion + vertical control)
2. Hybrid Learning Framework:
   • Integration with Balloon Learning Environment 17,18
   • Combine physics ROM advantages with reinforcement learning benefits
3. Uncertainty Quantification:
   • Integrate uncertainty estimation for ROM predictions
   • Robust or stochastic MPC formulations
4. Adaptive Reduced-Order Modeling:
   • Online POD mode updates to adapt to new flow patterns
   • Dynamic adjustment of retained mode count
5. Multi-Agent Coordination:
   • Extension to coordinated path planning for multiple agents
   • Exploit distributed sensing from agent swarms

In-Depth Evaluation

Strengths

1. Methodological Innovation:

• First systematic combination of POD-Galerkin ROM with MPC for NS flow path planning
• Clever time-varying measurement matrix design incorporating mobile agent as sensor
• Organic fusion of physics-informed and data-driven approaches

2. Theoretical Rigor:

• Clear complete derivation chain from NS equations
• Rigorous mathematical formulation of Galerkin projection
• Solid theoretical foundation for EKF observer design

3. Experimental Sufficiency:

• Real ERA5 data enhances credibility
• Systematic comparison across multiple prediction horizons
• Comprehensive quantitative metrics (distance, alignment, control cost, runtime)

4. Practical Value:

• Addresses computational bottleneck in real applications
• Provides framework directly applicable to stratospheric balloons and similar platforms
• Computational efficiency enables real-time control

5. Presentation Clarity:

• Clear structure with logical flow
• Well-defined mathematical notation
• Effective supporting figures and tables

Weaknesses

1. Method Limitations:

• Overly Strong Simplifications: Point mass model ignores complex aircraft dynamics
• Unidirectional Control: Vertical-only control assumption restricts applicability
• Lack of Theoretical Guarantees: No recursive feasibility or stability proofs

2. Experimental Design Flaws:

• Missing Baselines: No direct comparison with other methods (RL, traditional MPC, heuristics)
• Limited Scenarios: Single geographic region and season; generalization unknown
• Idealized Conditions: Ignores measurement noise, communication delays, actuator dynamics

3. Insufficient Analysis Depth:

• ROM Accuracy Impact: Lacks systematic analysis of mode count n on planning performance
• Missing Sensitivity Analysis: No principled guidance for cost weight selection (w_p, w_u)
• No Failure Cases: Doesn't demonstrate method failure modes or boundary conditions

4. Computational Efficiency Assessment:

• Limited Absolute Meaning: No comparison with high-fidelity methods
• Unknown Scalability: Performance with larger domains or higher resolution untested
• Real-Time Questionable: 2000-second computation for 12-hour horizon may be too long for some applications

5. Reproducibility Issues:

• Code not open-sourced
• Some hyperparameters (POD mode count n, EKF covariances Q and R) not explicitly provided
• Data preprocessing details (mean removal specifics) insufficiently detailed

Impact

Contribution to Field:

• Theoretical: Provides new theoretical framework for path planning in flow-driven systems
• Methodological: Demonstrates practical ROM-MPC potential for complex PDE-constrained problems
• Applied: Offers viable control solution for stratospheric balloons and similar systems

Practical Value:

• Medium-High: Method demonstrates good performance under ideal conditions but requires further engineering
• Technology Readiness Level: Approximately TRL 3-4 (concept validation), significant work needed for deployment

Reproducibility:

• Medium: Clear theoretical description but lacks code and detailed parameters; complete reproduction moderately challenging

Expected Impact:

• May inspire further ROM-MPC combination research
• Provides physics-informed methodology insights for autonomous navigation
• Potential adoption and improvement by stratospheric balloon and ocean glider communities

Applicable Scenarios

Ideal Application Scenarios:

1. Stratospheric Balloon Navigation: Original motivation; horizontal motion wind-dominated
2. Ocean Gliders: Similar vertical control + flow-driven horizontal motion characteristics
3. Atmospheric Sampling UAVs: Require energy-efficient path planning exploiting wind fields

Scenarios Requiring Improvement:

1. High-Maneuverability Aircraft: Need complete dynamics and multi-input control
2. Highly Transient Flow Fields: Turbulence and convection may violate POD linear subspace assumptions
3. Multi-Constraint Tasks: Obstacle avoidance, multiple objectives require extended optimization formulation

Inapplicable Scenarios:

1. Ground Mobile Robots: Flow effects negligible; method advantages unclear
2. High-Precision Point Tasks: Point mass model and unidirectional control insufficient
3. Extreme Real-Time Requirements: Current computation time too long for some applications

References

Key References:

1. POD Theory Foundations:
   • Lumley (1967): First POD application in turbulence analysis
   • Berkooz, Holmes, Lumley (1993): POD formalization in fluid mechanics
2. Reduced-Order Modeling:
   • Rowley (2005): Balanced POD for fluid modeling
   • Noack et al. (2011): Reduced-order modeling for flow control monograph
3. ROM-MPC Applications:
   • Benner, Gugercin, Willcox (2015): Parametric dynamical systems reduction survey
   • Yano (2020): Model reduction in computational aerodynamics
4. Practical Applications:
   • Greaves et al. (2021): Balloon Learning Environment
   • Oldewurtel et al. (2012): MPC in building energy control

Overall Assessment: This is a solid methodological paper that successfully combines classical POD reduced-order modeling with modern MPC technology to solve real-time path planning in flow-driven systems. The theoretical derivation is rigorous, experimental validation is thorough, and practical value is clear. Main limitations include strong simplifying assumptions, lack of direct method comparisons, and insufficient engineering details for deployment. The paper makes positive contributions to the control theory and fluid mechanics intersection, providing a good foundation for subsequent research. Recommended future work should prioritize dynamics integration, uncertainty quantification, and practical system validation.