Skip to content

feat: forced transition (xstrip) support #16

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
aeronauty-flexcompute merged 3 commits into from
Apr 8, 2026
Merged

feat: forced transition (xstrip) support #16

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
316c187
feat(flexfoil): add forced transition (xstrip) support
aeronauty Apr 8, 2026
9f893d2
chore: bump versions to 1.1.4/1.1.6 and add What's New for forced tra…
aeronauty Apr 8, 2026
ea7e1b3
merge: resolve conflict in python-api.mdx (keep xstrip params)
aeronauty Apr 8, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
11 changes: 9 additions & 2 deletions crates/rustfoil-coupling/src/global_newton.rs
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -138,6 +138,11 @@ pub struct GlobalNewtonSystem {
pub current_iteration: usize,
/// ANTE: base thickness contribution at blunt trailing edge (XFOIL WGAP(1))
pub ante: f64,
/// Forced transition arc-length on upper surface (from XSTRIP).
/// `None` falls back to TE arc length (default XFOIL behaviour).
pub xiforc_upper: Option<f64>,
/// Forced transition arc-length on lower surface (from XSTRIP).
pub xiforc_lower: Option<f64>,
}

impl GlobalNewtonSystem {
Expand Down Expand Up @@ -189,6 +194,8 @@ impl GlobalNewtonSystem {
// Current iteration - set in build_global_system
current_iteration: 0,
ante: 0.0,
xiforc_upper: None,
xiforc_lower: None,
}
}

Expand Down Expand Up @@ -790,9 +797,9 @@ impl GlobalNewtonSystem {
let is_transition_interval =
!s1.is_turbulent && !s1.is_wake && s2.is_turbulent && !s2.is_wake;
let x_forced = if surface == 0 {
stations.get(self.iblte_upper).map(|s| s.x)
self.xiforc_upper.or_else(|| stations.get(self.iblte_upper).map(|s| s.x))
} else {
stations.get(self.iblte_lower).map(|s| s.x)
self.xiforc_lower.or_else(|| stations.get(self.iblte_lower).map(|s| s.x))
};
let live_transition = if is_transition_interval {
Some(trchek2_full(
Expand Down
Loading
Loading