physmap 0.2.0__py3-none-any.whl

physmap/corpus/calibration.py

@@ -0,0 +1,543 @@
+"""PhysMAP Calibration-Boundary Corpus — schema, validators, loader, CLI.
+
+The data asset that makes the causal-layer guardrail real: a geometry-agnostic
+catalog of CFD closure validity envelopes, indexed by `closure_id`, with per-bound
+provenance and a `confirmed/claimed/extrapolated` status that is the asset.
+
+Spec: docs/specs/PhysMAP_CalibrationBoundary_Corpus_Spec_v0_1.md
+Storage: a single JSONL file (one closure entry per line, line-diffable).
+This module: schema (dataclasses + controlled vocabs), validators (the 7 QC gates +
+corpus-wide invariants), loader API (load / index_by_id / get_validated_range /
+is_in_calibration), and a CLI (`python -m physmap.corpus.calibration validate`).
+
+The integration seam: `Mechanism.in_calibration()` in stage1_ingest.py uses hardcoded
+`calib_lo/hi`. The corpus-backed replacement is `is_in_calibration(corpus_index,
+closure_id, coord, operating_value)`. The swap is mechanical once the corpus has
+Tier 1 coverage — deliberately not done yet.
+
+Discipline: correct-and-partial beats complete-and-unverified. A `confirmed` bound
+requires a primary-source citation (gate 7); the validator REJECTS any `confirmed`
+without it. Mass-loading training-knowledge entries as `claimed` is allowed and
+expected; passing them off as `confirmed` is forbidden.
+
+Torch-free; depends only on the standard library.
+"""
+
+from __future__ import annotations
+
+import argparse
+import importlib.resources as _ir
+import importlib.util as _iu
+import json
+import os
+import sys
+from dataclasses import asdict, dataclass, field
+from datetime import date
+from pathlib import Path
+
+# ── controlled vocabularies (mirror stage1_ingest.TRUTH_SOURCES pattern) ──────
+
+CLOSURE_FAMILIES = {
+    # two-phase / boiling
+    "flow-boiling", "pool-boiling", "chf-dnb",
+    "rpi-sub-closure", "interfacial-force",
+    # single-phase
+    "single-phase-convection", "rans-turbulence", "wall-function",
+    # medical
+    "hemolysis",
+    # multiphysics
+    "cavitation", "condensation",
+    "combustion-kinetics", "radiation",
+    "non-newtonian-rheology", "porous-media", "species-diffusion",
+    # hypersonic / high-speed
+    "boundary-layer-transition",
+}
+
+BOUND_STATUSES = {"confirmed", "confirmed-contested", "claimed", "extrapolated"}
+
+# The open SEED tier collapses every bound to this single status (the firewall
+# strips the confirmed/claimed narration). It is NOT a premium status — the
+# strict validator rejects it unless seed_tier=True is passed (opt-in), so the
+# premium corpus's confirmed-citation discipline stays byte-identical.
+SEED_BOUND_STATUS = "seed"
+
+SOLVER_NAMES = {"openfoam", "fluent", "starccm+", "cfx", "none"}
+
+# Geometry descriptors that MUST NOT appear in physics_coordinates per gate 6
+# (geometry lives only in regime_context as traceability, never as a coord).
+_FORBIDDEN_COORD_NAMES = {
+    "geometry", "shape", "diameter", "length", "channel_width",
+    "pipe_diameter", "channel_height", "aspect_ratio_geometry",
+}
+
+# Provenance shape requirement: a `confirmed` bound MUST carry both a primary_source
+# flag AND a non-empty citation (gate 7 — handbook restatement ≠ primary confirmation).
+# `provenance` is a dict; the expected keys are documented here, not strict-validated
+# beyond the gate-7 requirement so curators can carry extra context fields.
+_PROVENANCE_PRIMARY_KEY = "primary_source"
+_PROVENANCE_CITATION_KEYS = ("citation", "doi", "page", "table", "figure", "url", "year")
+
+# DEFAULT_PATH is the CANONICAL PREMIUM location and DOES NOT EXIST IN THIS
+# REPOSITORY. The premium corpus is the commercial moat and is not published
+# here, so resolution below always falls through to the bundled 15-closure seed
+# unless $PHYSMAP_CORPUS points somewhere else. That is the intended public
+# behaviour, not a misconfiguration. The path is kept so a premium holder can
+# drop the file in and have it picked up with no code change.
+# DEFAULT_PATH is the CANONICAL PREMIUM location — the full moat corpus in the
+# dev/editable tree. It is NOT shipped in the core wheel (the package-data glob
+# excludes results/calibration_corpus). Curation tooling and dev tests operate
+# on this path directly; runtime loading goes through resolve_corpus_path().
+DEFAULT_PATH = Path(__file__).resolve().parent.parent / "results" / "calibration_corpus" / "corpus.jsonl"
+
+
+# ── open-core corpus resolution (env → premium → bundled seed) ─────────────────
+PHYSMAP_CORPUS_ENV = "PHYSMAP_CORPUS"
+_PREMIUM_PACKAGE = "physmap_corpus_premium"       # optional private package (full corpus)
+_PREMIUM_FILENAME = "corpus_premium.jsonl"
+_SEED_PACKAGE = "physmap.corpus.data"             # bundled in the core wheel
+_SEED_FILENAME = "corpus_seed.jsonl"
+
+
+def _seed_path() -> Path:
+    """The bundled, firewalled SEED corpus (always present in the core wheel)."""
+    return Path(str(_ir.files(_SEED_PACKAGE) / _SEED_FILENAME))
+
+
+def _premium_package_path() -> Path | None:
+    """Full corpus from the optional `physmap_corpus_premium` package, or None."""
+    if _iu.find_spec(_PREMIUM_PACKAGE) is None:
+        return None
+    try:
+        p = Path(str(_ir.files(_PREMIUM_PACKAGE) / _PREMIUM_FILENAME))
+        return p if p.is_file() else None
+    except (FileNotFoundError, ModuleNotFoundError, NotADirectoryError):
+        return None
+
+
+def resolve_corpus_path() -> Path:
+    """Resolve the active calibration corpus, in precedence order:
+
+      1. ``$PHYSMAP_CORPUS`` — explicit path override.
+      2. Premium if present — the canonical dev/editable file (DEFAULT_PATH), or
+         the installed ``physmap_corpus_premium`` package.
+      3. The bundled seed — ``physmap/corpus/data/corpus_seed.jsonl``.
+
+    No network calls, no license checks: tier detection is purely "is the
+    premium file present". In an editable/dev checkout step 2 hits DEFAULT_PATH,
+    so every existing consumer sees the full corpus unchanged; in a core-only
+    wheel install the canonical file is absent and resolution falls to the seed.
+    """
+    env = os.environ.get(PHYSMAP_CORPUS_ENV)
+    if env:
+        return Path(env).expanduser()
+    if DEFAULT_PATH.is_file():
+        return DEFAULT_PATH
+    pkg = _premium_package_path()
+    if pkg is not None:
+        return pkg
+    return _seed_path()
+
+
+def active_tier(path: str | Path | None = None) -> str:
+    """Return "seed" or "premium" for the resolved (or given) corpus.
+
+    Content-based: the seed firewall collapses EVERY bound to status "seed", so
+    a corpus whose first entry is all-"seed" is the seed tier. This is correct
+    even for an explicit ``$PHYSMAP_CORPUS`` override regardless of filename.
+    """
+    p = Path(path) if path is not None else resolve_corpus_path()
+    try:
+        with open(p, encoding="utf-8") as f:
+            for line in f:
+                line = line.strip()
+                if not line:
+                    continue
+                statuses = {b.get("bound_status")
+                            for b in json.loads(line).get("validated_range", [])}
+                return "seed" if statuses == {SEED_BOUND_STATUS} else "premium"
+    except (FileNotFoundError, json.JSONDecodeError):
+        pass
+    return "premium"
+
+
+# ── schema dataclasses ────────────────────────────────────────────────────────
+
+@dataclass
+class CoordinateBound:
+    """One per coordinate of a closure's validity envelope."""
+    coord: str                                  # e.g. "mass_flux", "reynolds_number"
+    min: float | None                           # None = unbounded below
+    max: float | None                           # None = unbounded above
+    unit: str                                   # SI preferred; "dimensionless" for Re/Pr
+    bound_status: str                           # ∈ BOUND_STATUSES
+    provenance: dict                            # see _PROVENANCE_* constants above
+    contested_note: str | None = None           # if sources disagree on this coord
+
+
+@dataclass
+class ClosureEntry:
+    """One named closure: validity envelope + provenance + status + scope filter."""
+    closure_id: str                             # immutable slug: family-name-year
+    closure_family: str                         # ∈ CLOSURE_FAMILIES
+    closure_name: str                           # human name
+    physics_coordinates: list[str]              # ordered list of coord names
+    validated_range: list[CoordinateBound]      # one per coord; order matches
+    fluid_context: list[str]                    # ["water","R12","blood","air",...]
+    regime_context: str | None                  # traceability only; geometry-OK here
+    contested: bool                             # whole-entry disagreement flag
+    solver_availability: list[str]              # ⊂ SOLVER_NAMES
+    weakener_link: list[str]                    # UofA weakener-pattern ids
+    last_reviewed: str                          # ISO date YYYY-MM-DD
+    reviewer: str                               # who curated/verified
+
+
+# ── (de)serialization ─────────────────────────────────────────────────────────
+
+def entry_from_dict(d: dict) -> ClosureEntry:
+    """Build a ClosureEntry from a JSON-loaded dict (handles nested CoordinateBound)."""
+    bounds_raw = d.get("validated_range", [])
+    bounds = [CoordinateBound(**b) for b in bounds_raw]
+    payload = {**d, "validated_range": bounds}
+    return ClosureEntry(**payload)
+
+
+def entry_to_dict(entry: ClosureEntry) -> dict:
+    return asdict(entry)
+
+
+# ── per-entry validation (the 7 QC gates) ─────────────────────────────────────
+
+def validate(entry: ClosureEntry, *, seed_tier: bool = False) -> list[str]:
+    """Return a list of human-readable error strings; empty list = entry is valid.
+
+    Implements the spec's 7 per-entry QC gates plus shape checks.
+
+    seed_tier=True (opt-in) validates a sanitized SEED file: it additionally
+    accepts bound_status "seed". Everything else is unchanged, so premium
+    validation is byte-identical to before (gate 7 only ever fires on
+    confirmed/confirmed-contested, which the seed never carries).
+    """
+    errs: list[str] = []
+    allowed_statuses = BOUND_STATUSES | {SEED_BOUND_STATUS} if seed_tier else BOUND_STATUSES
+    eid = entry.closure_id or "<no closure_id>"
+
+    # Shape: closure_id present, immutable slug
+    if not entry.closure_id or not isinstance(entry.closure_id, str):
+        errs.append(f"{eid}: closure_id missing or not a string")
+    elif " " in entry.closure_id or entry.closure_id != entry.closure_id.lower():
+        errs.append(f"{eid}: closure_id must be a lowercase slug (no spaces)")
+
+    # Controlled vocab: family
+    if entry.closure_family not in CLOSURE_FAMILIES:
+        errs.append(
+            f"{eid}: closure_family '{entry.closure_family}' not in CLOSURE_FAMILIES "
+            f"({sorted(CLOSURE_FAMILIES)})"
+        )
+
+    # Shape: closure_name present
+    if not entry.closure_name:
+        errs.append(f"{eid}: closure_name missing")
+
+    # Shape: physics_coordinates and validated_range length match
+    if len(entry.physics_coordinates) != len(entry.validated_range):
+        errs.append(
+            f"{eid}: physics_coordinates ({len(entry.physics_coordinates)}) and "
+            f"validated_range ({len(entry.validated_range)}) length mismatch"
+        )
+
+    # Shape + Gate 6: coord names in validated_range match physics_coordinates list,
+    # and no geometry descriptors appear (geometry lives only in regime_context).
+    for coord_name in entry.physics_coordinates:
+        if coord_name.lower() in _FORBIDDEN_COORD_NAMES:
+            errs.append(
+                f"{eid}: physics_coordinates contains forbidden geometry descriptor "
+                f"'{coord_name}' (gate 6 — geometry lives only in regime_context)"
+            )
+    coord_names_in_bounds = {b.coord for b in entry.validated_range}
+    coord_names_declared = set(entry.physics_coordinates)
+    if coord_names_in_bounds != coord_names_declared:
+        only_in_bounds = coord_names_in_bounds - coord_names_declared
+        only_in_decl = coord_names_declared - coord_names_in_bounds
+        errs.append(
+            f"{eid}: physics_coordinates vs validated_range.coord mismatch "
+            f"(only in bounds: {sorted(only_in_bounds)}; only in declaration: "
+            f"{sorted(only_in_decl)})"
+        )
+
+    # Per-bound gates 1, 2, 3, 7
+    for b in entry.validated_range:
+        # Gate 1: every bound carries a provenance
+        if not isinstance(b.provenance, dict) or not b.provenance:
+            errs.append(
+                f"{eid}: coord '{b.coord}' missing provenance "
+                f"(gate 1 — every bound carries a citation)"
+            )
+        else:
+            # Gate 7: confirmed (and confirmed-contested) bounds require
+            # primary_source + a non-handbook citation. confirmed-contested
+            # is still a confirmation — source(s) disagree, but the primary-
+            # source provenance discipline applies the same way.
+            if b.bound_status in ("confirmed", "confirmed-contested"):
+                if not b.provenance.get(_PROVENANCE_PRIMARY_KEY):
+                    errs.append(
+                        f"{eid}: coord '{b.coord}' is '{b.bound_status}' but provenance "
+                        f"missing primary_source=true "
+                        f"(gate 7 — handbook restatement ≠ primary confirmation)"
+                    )
+                has_citation = any(
+                    b.provenance.get(k) for k in _PROVENANCE_CITATION_KEYS
+                )
+                if not has_citation:
+                    errs.append(
+                        f"{eid}: coord '{b.coord}' is '{b.bound_status}' but provenance "
+                        f"carries no citation (one of {_PROVENANCE_CITATION_KEYS})"
+                    )
+
+        # Gate 2: bound_status ∈ BOUND_STATUSES (claimed default until upgraded);
+        # seed_tier additionally allows the collapsed "seed" status.
+        if b.bound_status not in allowed_statuses:
+            errs.append(
+                f"{eid}: coord '{b.coord}' bound_status '{b.bound_status}' not in "
+                f"{sorted(allowed_statuses)} (gate 2)"
+            )
+
+        # Gate 3: a contested_note on a bound means contested=True at entry level
+        if b.contested_note and not entry.contested:
+            errs.append(
+                f"{eid}: coord '{b.coord}' has contested_note but entry.contested=False "
+                f"(gate 3 — record disagreement, do NOT silently average)"
+            )
+
+        # Shape: unit present
+        if not b.unit:
+            errs.append(f"{eid}: coord '{b.coord}' missing unit")
+
+        # Shape: at least one of min/max set
+        if b.min is None and b.max is None:
+            errs.append(
+                f"{eid}: coord '{b.coord}' has both min and max unset "
+                f"(at least one bound must be declared)"
+            )
+        if b.min is not None and b.max is not None and b.min > b.max:
+            errs.append(f"{eid}: coord '{b.coord}' min ({b.min}) > max ({b.max})")
+
+    # Gate 4: solver_availability set; if {"none"} alone, require handbook adoption note
+    if not entry.solver_availability:
+        errs.append(f"{eid}: solver_availability empty (gate 4 — scope filter)")
+    else:
+        for s in entry.solver_availability:
+            if s not in SOLVER_NAMES:
+                errs.append(
+                    f"{eid}: solver_availability includes '{s}' not in "
+                    f"{sorted(SOLVER_NAMES)}"
+                )
+        if entry.solver_availability == ["none"]:
+            # Gate 4 requires an explicit handbook-ADOPTION marker, not just any mention
+            # of "handbook" — otherwise "no handbook reference here" trivially passes.
+            adoption_markers = (
+                "handbook-adopted", "handbook adoption",
+                "handbook-canonical", "canonical handbook",
+                "widely re-cited", "widely cited", "widely adopted",
+            )
+            ctx = (entry.regime_context or "").lower()
+            if not any(marker in ctx for marker in adoption_markers):
+                errs.append(
+                    f"{eid}: solver_availability=['none'] requires an explicit handbook-"
+                    f"adoption marker in regime_context (one of: "
+                    f"{', '.join(repr(m) for m in adoption_markers)}); gate 4 — keeps "
+                    f"the count at thousands not hundreds-of-thousands"
+                )
+
+    # Gate 5: weakener_link populated
+    if not entry.weakener_link:
+        errs.append(
+            f"{eid}: weakener_link empty (gate 5 — bridge to UofA weakener vocabulary)"
+        )
+
+    # Shape: last_reviewed parseable ISO date
+    if not entry.last_reviewed:
+        errs.append(f"{eid}: last_reviewed missing")
+    else:
+        try:
+            date.fromisoformat(entry.last_reviewed)
+        except ValueError as e:
+            errs.append(
+                f"{eid}: last_reviewed '{entry.last_reviewed}' not parseable ISO date "
+                f"({e})"
+            )
+
+    # Shape: reviewer present
+    if not entry.reviewer:
+        errs.append(f"{eid}: reviewer missing")
+
+    return errs
+
+
+# ── corpus-wide invariants ────────────────────────────────────────────────────
+
+def validate_corpus(entries: list[ClosureEntry], *, seed_tier: bool = False) -> dict[str, list[str]]:
+    """Per-entry errors keyed by closure_id, plus '__corpus__' for cross-entry issues.
+
+    Cross-entry checks: unique closure_id. seed_tier is forwarded to validate()."""
+    out: dict[str, list[str]] = {}
+    for e in entries:
+        errs = validate(e, seed_tier=seed_tier)
+        if errs:
+            out[e.closure_id or "<no closure_id>"] = errs
+
+    # Unique closure_id across the corpus
+    seen: dict[str, int] = {}
+    for e in entries:
+        seen[e.closure_id] = seen.get(e.closure_id, 0) + 1
+    dups = [cid for cid, n in seen.items() if n > 1]
+    if dups:
+        out.setdefault("__corpus__", []).append(
+            f"duplicate closure_id values: {sorted(dups)}"
+        )
+
+    return out
+
+
+# ── loader API (minimal; future-causal-layer-integration-ready) ───────────────
+
+def load_corpus(path: str | Path | None = None) -> list[ClosureEntry]:
+    """Parse the JSONL file into a list of ClosureEntry. Raises on parse errors;
+    schema validity must be checked separately via validate_corpus().
+
+    path=None resolves the active corpus via resolve_corpus_path() (env → premium
+    → bundled seed). In a dev/editable tree that is the canonical premium file, so
+    behaviour is unchanged; in a core-only install it is the bundled seed."""
+    p = Path(path) if path is not None else resolve_corpus_path()
+    if not p.exists():
+        raise FileNotFoundError(f"corpus file not found: {p}")
+    entries: list[ClosureEntry] = []
+    with open(p, "r", encoding="utf-8") as f:
+        for lineno, line in enumerate(f, start=1):
+            stripped = line.strip()
+            if not stripped:
+                continue
+            try:
+                d = json.loads(stripped)
+            except json.JSONDecodeError as e:
+                raise ValueError(f"corpus JSONL parse error at line {lineno}: {e}") from e
+            try:
+                entries.append(entry_from_dict(d))
+            except TypeError as e:
+                raise ValueError(
+                    f"corpus schema mismatch at line {lineno} "
+                    f"(closure_id={d.get('closure_id', '?')!r}): {e}"
+                ) from e
+    return entries
+
+
+def index_by_id(entries: list[ClosureEntry]) -> dict[str, ClosureEntry]:
+    return {e.closure_id: e for e in entries}
+
+
+def get_validated_range(
+    entries_index: dict[str, ClosureEntry], closure_id: str, coord: str,
+) -> CoordinateBound | None:
+    """Return the CoordinateBound for (closure_id, coord), or None if unknown."""
+    entry = entries_index.get(closure_id)
+    if entry is None:
+        return None
+    for b in entry.validated_range:
+        if b.coord == coord:
+            return b
+    return None
+
+
+def is_in_calibration(
+    entries_index: dict[str, ClosureEntry], closure_id: str, coord: str, value: float,
+) -> bool | None:
+    """Return True/False if the value is inside [min, max] for the (closure_id, coord)
+    bound, or None if the closure_id or coord is not in the corpus. None signals the
+    caller to fall back to its own bounds (Mechanism.calib_lo/hi today)."""
+    b = get_validated_range(entries_index, closure_id, coord)
+    if b is None:
+        return None
+    if b.min is not None and value < b.min:
+        return False
+    if b.max is not None and value > b.max:
+        return False
+    return True
+
+
+# ── CLI ───────────────────────────────────────────────────────────────────────
+
+def _print_errors(errs_by_id: dict[str, list[str]]) -> int:
+    """Print errors grouped by closure_id; return non-zero exit code if any."""
+    if not errs_by_id:
+        return 0
+    total = sum(len(v) for v in errs_by_id.values())
+    print(f"corpus validation FAILED — {total} errors across {len(errs_by_id)} entries:")
+    for cid in sorted(errs_by_id):
+        print(f"  [{cid}]")
+        for err in errs_by_id[cid]:
+            print(f"      • {err}")
+    return 1
+
+
+def main(argv: list[str] | None = None) -> int:
+    parser = argparse.ArgumentParser(
+        description="PhysMAP Calibration-Boundary Corpus — validator + summary CLI",
+    )
+    sub = parser.add_subparsers(dest="cmd", required=True)
+
+    p_val = sub.add_parser("validate", help="run the 7 QC gates + corpus invariants")
+    p_val.add_argument("--path", type=Path, default=None,
+                       help="path to corpus.jsonl (default: the active corpus via "
+                            "resolve_corpus_path — env → premium → bundled seed)")
+    p_val.add_argument("--tier", choices=("auto", "seed", "premium"), default="auto",
+                       help="validation tier; 'auto' detects 'seed' (all-seed bounds) "
+                            "and applies the seed-tier gates")
+
+    p_sum = sub.add_parser("summary", help="print a status summary of the corpus")
+    p_sum.add_argument("--path", type=Path, default=None)
+
+    args = parser.parse_args(argv)
+
+    if args.cmd == "validate":
+        path = args.path or resolve_corpus_path()
+        try:
+            entries = load_corpus(path)
+        except (FileNotFoundError, ValueError) as e:
+            print(f"FAILED to load corpus: {e}", file=sys.stderr)
+            return 2
+        seed_tier = args.tier == "seed" or (args.tier == "auto" and active_tier(path) == "seed")
+        errs_by_id = validate_corpus(entries, seed_tier=seed_tier)
+        if not errs_by_id:
+            tier = "seed" if seed_tier else "premium"
+            print(f"corpus validation PASSED — {len(entries)} entries clean "
+                  f"[{tier}] ({path})")
+            # Status summary by bound_status (useful audit signal)
+            status_counts: dict[str, int] = {}
+            for e in entries:
+                for b in e.validated_range:
+                    status_counts[b.bound_status] = status_counts.get(b.bound_status, 0) + 1
+            total_bounds = sum(status_counts.values())
+            print(f"  total bounds: {total_bounds}")
+            display_statuses = sorted(BOUND_STATUSES | ({SEED_BOUND_STATUS} if seed_tier else set()))
+            for s in display_statuses:
+                print(f"    {s:>12s} : {status_counts.get(s, 0)}")
+            print(f"  contested entries: "
+                  f"{sum(1 for e in entries if e.contested)}/{len(entries)}")
+            return 0
+        return _print_errors(errs_by_id)
+
+    if args.cmd == "summary":
+        path = args.path or resolve_corpus_path()
+        entries = load_corpus(path)
+        family_counts: dict[str, int] = {}
+        for e in entries:
+            family_counts[e.closure_family] = family_counts.get(e.closure_family, 0) + 1
+        print(f"corpus: {len(entries)} entries  [{active_tier(path)}] ({path})")
+        for fam in sorted(family_counts):
+            print(f"  {fam:>26s} : {family_counts[fam]}")
+        return 0
+
+    return 0
+
+
+if __name__ == "__main__":
+    raise SystemExit(main())

physmap/corpus/data/__init__.py

@@ -0,0 +1,12 @@
+"""Bundled open-tier corpus data (package resources).
+
+Holds the sanitized SEED corpus (`corpus_seed.jsonl`) and the premium COVERAGE
+map (`premium_coverage.json`, closure_ids only — no bounds). Both are generated
+by `dev/tools/split_corpus.py` and shipped in the core wheel as package data;
+the loader resolves them via `importlib.resources` when neither an explicit
+`$PHYSMAP_CORPUS` override nor a premium corpus is present.
+
+This package carries NO moat content: the seed is firewalled (every bound
+collapsed to status "seed", all confirmation/disposition stripped) and the
+coverage map is ids-only.
+"""

physmap/corpus/data/corpus_seed.jsonl

@@ -0,0 +1,15 @@
+{"closure_id": "aung-worku-mixed-convection-1986", "closure_family": "single-phase-convection", "closure_name": "Aung-Worku mixed-convection Nusselt correlation (Eq13 family)", "physics_coordinates": ["richardson_number"], "validated_range": [{"coord": "richardson_number", "min": 0.1, "max": 10.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Aung, W. & Worku, G. (1986). 'Theory of fully developed combined convection including flow reversal.' J. Heat Transfer 108:485-488", "doi": "10.1115/1.3246958", "primary_source": false, "note": "mixed-convection regime is typically reported as Ri \u2208 [0.1, 10]; classical 'mixed' regime boundaries"}}], "fluid_context": ["air", "water"], "regime_context": "vertical channel, fully developed mixed forced-natural convection; Eq13 family of correlations used in conv_step3_final validation", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-03", "reviewer": "seed"}
+{"closure_id": "blasius-pohlhausen-flat-plate-forced-1921", "closure_family": "single-phase-convection", "closure_name": "Blasius/Pohlhausen flat-plate laminar forced convection", "physics_coordinates": ["reynolds_number", "prandtl_number", "richardson_number"], "validated_range": [{"coord": "reynolds_number", "min": 0.0, "max": 500000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Pohlhausen, E. (1921). 'Der Waermeaustausch zwischen festen Koerpern und Fluessigkeiten mit kleiner Reibung und kleiner Waermeleitung.' Zeitschrift fuer Angewandte Mathematik und Mechanik 1(2):115-121. Velocity-field foundation: Blasius, H. (1908). 'Grenzschichten in Fluessigkeiten mit kleiner Reibung.' Z. Math. Phys. 56:1-37. Restated in: Bergman, T. L. et al. (2011). Fundamentals of Heat and Mass Transfer (7th ed.), Section 7.2; Schlichting, H. & Gersten, K. (2017). Boundary-Layer Theory (9th ed.), Ch. 9.", "primary_source": true, "year": 1921, "note": "Analytical correlation. Pohlhausen's thermal extension of Blasius's velocity boundary layer. Local Nu_x = 0.332*Re_x^(1/2)*Pr^(1/3); length-averaged Nu_L = 0.664*Re_L^(1/2)*Pr^(1/3). Re upper bound 5e5 is the standard laminar-to-turbulent boundary-layer transition on a smooth isothermal flat plate; above this is transitional/turbulent regime where Blasius/Pohlhausen do not apply."}}, {"coord": "prandtl_number", "min": 0.6, "max": 50.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Bergman, T. L. et al. (2011). Fundamentals of Heat and Mass Transfer (7th ed.), Section 7.2, Eqs. 7.30 and 7.31. Lower bound 0.6 derives from Pohlhausen's asymptotic analysis (assumes Pr = O(1) or larger; for liquid metals with Pr << 1, a separate Pohlhausen low-Pr correlation applies). Upper bound 50 is pragmatic textbook usage.", "primary_source": true, "year": 2011, "page": "Section 7.2", "note": "Lance & Smith vehicle uses air (Pr ~ 0.71), well inside the validated range."}}, {"coord": "richardson_number", "min": null, "max": 0.225, "unit": "dimensionless", "bound_status": "seed", "provenance": {"primary_source": false, "citation": "Sparrow, E.M. & Gregg, J.L. (1959). 'Buoyancy Effects in Forced-Convection Flow and Heat Transfer.' J. Appl. Mech. 26: 133-134", "year": 1959, "discovery_path": {"handbook": "Rohsenow, W.M., Hartnett, J.P. & Cho, Y.I. (eds.) (1998). Handbook of Heat Transfer, 3rd ed. McGraw-Hill. ISBN 0-07-053555-8", "handbook_table": "\u00a74.9 Mixed Convection \u2014 External Flow / Vertical Plates, Eq. 4.158", "handbook_page": "4.75 (PDF p.275)"}, "note": "Buoyancy-validity envelope of the forced-convection flat-plate result: the pure forced-convection Nusselt number stays within 5% of the combined value only while the Richardson number Gr/Re^2 remains at or below 0.225; beyond that, assisting or opposing buoyancy must be accounted for and a mixed-convection correlation used. Per Sparrow & Gregg's perturbation analysis, valid for Prandtl 0.01 to 10. Gr/Re^2 is the Richardson number for mixed convection."}}], "fluid_context": ["air", "water", "fluid-Pr-ge-0.6"], "regime_context": "external flow over an isothermal flat plate, laminar boundary layer (Re_x <= 5e5); handbook-canonical analytical correlation widely cited across all major heat-transfer textbooks (Bergman, Incropera, Schlichting, Kays & Crawford)", "contested": false, "solver_availability": ["none"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-05", "reviewer": "seed"}
+{"closure_id": "dittus-boelter-1930", "closure_family": "single-phase-convection", "closure_name": "Dittus-Boelter forced-convection Nusselt", "physics_coordinates": ["reynolds_number", "prandtl_number"], "validated_range": [{"coord": "reynolds_number", "min": 10000.0, "max": 1200000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Dittus, F.W. & Boelter, L.M.K. (1930). 'Heat transfer in automobile radiators of tubular type.' Univ. California Publ. Engr. 2(13):443-461", "primary_source": false, "note": "Nu = 0.023 Re^0.8 Pr^n; the canonical single-phase forced-convection correlation, ranges per handbook restatement"}}, {"coord": "prandtl_number", "min": 0.7, "max": 160.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Dittus-Boelter 1930", "primary_source": false, "note": "Pr range per handbook restatement; primary-source verification pending"}}], "fluid_context": ["water", "air", "oils"], "regime_context": "fully developed turbulent flow in smooth tubes; small temperature differences (no property-variation correction)", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-03", "reviewer": "seed"}
+{"closure_id": "dittus-boelter-buoyancy-sco2", "closure_family": "single-phase-convection", "closure_name": "Constant-property Dittus-Boelter (supercritical CO2, buoyancy-bounded)", "physics_coordinates": ["reynolds_number", "prandtl_number", "liu_buoyancy_parameter"], "validated_range": [{"coord": "reynolds_number", "min": 10000.0, "max": 1200000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Dittus-Boelter (McAdams form); Re turbulent window per dittus-boelter-1930.", "primary_source": false, "note": "Same Re window as dittus-boelter-1930; broad sCO2 turbulent envelope, NOT the discriminator (Jin Re_b ~3.5e4-9.3e4 is interior)."}}, {"coord": "prandtl_number", "min": 0.5, "max": 2000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Broad sCO2 turbulent envelope (Pr spikes near the pseudo-critical point).", "primary_source": false, "note": "Broad on purpose so the buoyancy coordinate is this regime-specific entry's discriminator (mirrors gnielinski-constprop-sco2); Jin Pr_b ~2.5-16 is interior."}}, {"coord": "liu_buoyancy_parameter", "min": 0.0, "max": 1.3e-05, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Liu et al. 2017 buoyancy parameter Bu; Jin et al. 2023 (Annals of Nuclear Energy 188:109825) Eq. 21 + buoyancy criteria (deterioration onset near Bu = 1.3e-5). Banked evidence Claim jin-2023-correct-dittus-boelter-buoyancy-bound. Set A PRIORI from the paper's stated HTD-onset threshold, NOT tuned to data.", "primary_source": false, "note": "Constant-property Dittus-Boelter validity bound for BUOYANCY. Fires (signal>0) when Bu exceeds 1.3e-5 (buoyancy-significant / heat-transfer-deterioration regime). Bu (Liu, wall-aware) is a DISTINCT physical group from richardson_number (Gr/Re^2) \u2014 the same-name trap is avoided by this dedicated coordinate."}}], "fluid_context": ["supercritical CO2"], "regime_context": "Constant-property turbulent pipe Nu (Dittus-Boelter) applied to supercritical CO2 in a VERTICAL heated tube; valid when buoyancy is negligible (Liu Bu < 1.3e-5), INVALID once buoyancy is significant (upward flow -> heat-transfer deterioration; downward -> enhancement) which the direction-blind bare correlation cannot predict. Jin et al. 2023.", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-10", "reviewer": "seed"}
+{"closure_id": "gnielinski-1976", "closure_family": "single-phase-convection", "closure_name": "Gnielinski forced-convection Nusselt", "physics_coordinates": ["reynolds_number", "prandtl_number", "x_over_D", "prandtl_ratio_bulk_wall", "temperature_ratio_bulk_wall"], "validated_range": [{"coord": "reynolds_number", "min": 3000.0, "max": 5000000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Gnielinski, V. (1976). 'New equations for heat and mass transfer in turbulent pipe and channel flow.' Int. Chem. Engr. 16:359-368", "primary_source": true, "note": "Gnielinski's modified Petukhov correlation; transition + turbulent; ranges per handbook restatement"}}, {"coord": "prandtl_number", "min": 0.5, "max": 2000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Gnielinski 1976 ICE", "primary_source": true, "note": "wider Pr range than Dittus-Boelter; primary-source verification pending"}}, {"coord": "x_over_D", "min": 10.0, "max": 1000000000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Bergman, T. L., Lavine, A. S., Incropera, F. P., DeWitt, D. P. (2011). Fundamentals of Heat and Mass Transfer (7th ed.), Section 8.5, p. 519. Cites the L/D > 10 criterion for fully-developed turbulent flow in circular tubes, applicable to Gnielinski-style correlations. The criterion is textbook-stable: Kays & Crawford 'Convective Heat and Mass Transfer' 4th ed. Ch. 14 gives the same L/D > 10 for thermal entry length in turbulent gas flow.", "primary_source": false, "year": 2011, "note": "L/D > 10 is a handbook usage guideline for the fully-developed turbulent flow assumption underlying Gnielinski 1976. Status 'claimed' (not 'confirmed') because the cutoff is a handbook judgment about when the Hausen entrance-correction factor [1 + (D/L)^(2/3)] becomes negligible (< 5% effect at L/D > 10 for gases). Hausen 1943 gave the correction; the L/D > 10 threshold is the handbook's interpretation of 'negligibly small entrance effect.' For high-Pr liquids the thermal entry length is longer (handbook quote L/D > 60 for water); this bound uses 10 as the lower-end conservative cutoff appropriate for gas-flow vehicles like NACA TN-1451 (air, Pr~0.71)."}}, {"coord": "prandtl_ratio_bulk_wall", "min": 0.1, "max": 10.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"primary_source": false, "citation": "Hufschmidt, W. & Burck, E. (1968). 'Einfluss temperaturabhaengiger Stoffwerte auf den Waermeuebergang bei turbulenter Stroemung von Fluessigkeiten in Rohren bei hohen Waermestromdichten und Prandtlzahlen.' Int. J. Heat Mass Transfer 11:1041-1048 -- origin of the (Pr/Pr_w)^0.11 factor, applied here to the turbulent Gnielinski 1976 correlation; validity window stated in VDI G1 sec.4.4.", "year": 1968, "discovery_path": {"handbook": "VDI Heat Atlas (2nd English ed., Springer-Verlag Berlin Heidelberg, 2010; ed. VDI-GVC)", "handbook_table": "G1 'Heat Transfer in Pipe Flow' (V. Gnielinski), sec.4.4 'Effect of Temperature Dependent Property Variations', Eq.(40)", "handbook_page": "697 (PDF p720)"}}}, {"coord": "temperature_ratio_bulk_wall", "min": 0.5, "max": 1.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"primary_source": false, "citation": "Gnielinski, V. (1976). 'New equations for heat and mass transfer in turbulent pipe and channel flow.' Int. Chem. Engr. 16:359-368 -- fitted the gas correction exponent n=0.45 over the heating range 0.5<T/Tw<1; restated in VDI G1 sec.4.4.", "year": 1976, "discovery_path": {"handbook": "VDI Heat Atlas (2nd English ed., Springer-Verlag Berlin Heidelberg, 2010; ed. VDI-GVC)", "handbook_table": "G1 'Heat Transfer in Pipe Flow' (V. Gnielinski), sec.4.4 'Effect of Temperature Dependent Property Variations', Eq.(41)", "handbook_page": "697 (PDF p720)"}}}], "fluid_context": ["water", "air", "oils", "molten salts"], "regime_context": "transition + fully turbulent pipe/channel flow; smooth walls; fully-developed thermal profile (x/D > 10 lower bound). For x/D < 10 (entrance region) the Hausen correction factor [1 + (D/L)^(2/3)] is required; without it the correlation under-predicts Nu by factor 1.2-2.0 in the developing region.", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-05", "reviewer": "seed"}
+{"closure_id": "gnielinski-constprop-sco2", "closure_family": "single-phase-convection", "closure_name": "Constant-property Gnielinski (supercritical CO2, property-variation-bounded)", "physics_coordinates": ["reynolds_number", "prandtl_number", "viscosity_ratio_wall_bulk"], "validated_range": [{"coord": "reynolds_number", "min": 3000.0, "max": 5000000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Gnielinski 1976 ICE; Re range per Bergman/Incropera sec 8.5", "primary_source": false, "note": "Same Re validity window as gnielinski-1976 (constant-property turbulent pipe); claimed here for this regime-specific entry."}}, {"coord": "prandtl_number", "min": 0.5, "max": 2000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Gnielinski 1976 ICE", "primary_source": false, "note": "Same Pr validity window as gnielinski-1976."}}, {"coord": "viscosity_ratio_wall_bulk", "min": 0.9, "max": 1.1, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Kays, W.M. & Crawford, M.E., 'Convective Heat and Mass Transfer' (property-ratio method). Constant-property correlations (Gnielinski/Dittus-Boelter) omit the (mu_b/mu_w)^n property correction and are stated for SMALL property variation. Constant-property validity is conventionally taken for viscosity-ratio departures within ~+/-10% of unity (mu_w/mu_b in [0.9, 1.1]); beyond that (e.g. near the supercritical pseudo-critical point) the property-ratio/reference-temperature correction is required and the bare correlation is out of validity. Status 'claimed' (handbook judgment, like the gnielinski-1976 x/D>10 bound), set A PRIORI from the property-ratio-method threshold, NOT tuned to the Velazquez data.", "primary_source": false, "note": "Property-variation validity bound for the constant-property assumption. Fires (signal>0) when mu_w/mu_b leaves [0.9,1.1], i.e. strong property variation near pseudo-critical."}}], "fluid_context": ["supercritical CO2"], "regime_context": "Constant-property turbulent pipe Nu applied to supercritical CO2; valid for small property variation (viscosity ratio near unity), INVALID near the pseudo-critical point where rho/mu/cp/lambda vary sharply and the constant-property assumption breaks (heat-transfer deterioration/enhancement the bare correlation cannot predict).", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-07", "reviewer": "seed"}
+{"closure_id": "laminar-fully-developed-tube-isothermal-nu366", "closure_family": "single-phase-convection", "closure_name": "Laminar fully-developed circular-tube Nusselt, uniform surface temperature (Nu=3.66, analytical)", "physics_coordinates": ["reynolds_number"], "validated_range": [{"coord": "reynolds_number", "min": null, "max": 2300.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"primary_source": false, "citation": "Shah, R.K. & London, A.L. (1978). Laminar Flow Forced Convection in Ducts (Supplement 1 to Advances in Heat Transfer). Academic Press, New York", "year": 1978, "page": "circular-tube fully-developed Nusselt numbers (Nu_T=3.657 const-T; Nu_H1=4.364 const-q)", "discovery_path": {"handbook": "Incropera, F.P., DeWitt, D.P., Bergman, T.L. & Lavine, A.S. (2007). Fundamentals of Heat and Mass Transfer, 6th ed. Wiley", "handbook_table": "Eq. 8.55 (Nu_D=3.66) in Table 8.4 Summary of Convection Correlations for Flow in a Circular Tube; laminar limit Re_D,c~2300 from Eq. 8.2", "handbook_page": "p.532 (Table 8.4); p.487 (Eq. 8.2 transition Reynolds)"}, "note": "Laminar-regime validity limit for the analytical fully-developed constant-Nusselt result Nu_D=3.66 (uniform/constant surface temperature) in a circular tube: it applies while the tube Reynolds number stays laminar, up to the transition-onset value near 2300. The 3.66 is the fully-developed analytical asymptote."}}], "fluid_context": ["any single-phase fluid"], "regime_context": "laminar, hydrodynamically and thermally fully-developed flow in a circular tube held at a uniform (constant) surface temperature; the analytical fully-developed Nusselt number is the constant 3.66; valid below the laminar-to-turbulent transition (critical Reynolds ~2300)", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure", "closure-set-mismatch"], "last_reviewed": "2026-06-08", "reviewer": "seed"}
+{"closure_id": "laminar-fully-developed-tube-uniform-flux-nu436", "closure_family": "single-phase-convection", "closure_name": "Laminar fully-developed circular-tube Nusselt, uniform surface heat flux (Nu=4.36, analytical)", "physics_coordinates": ["reynolds_number"], "validated_range": [{"coord": "reynolds_number", "min": null, "max": 2300.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"primary_source": false, "citation": "Shah, R.K. & London, A.L. (1978). Laminar Flow Forced Convection in Ducts (Supplement 1 to Advances in Heat Transfer). Academic Press, New York", "year": 1978, "page": "circular-tube fully-developed Nusselt numbers (Nu_T=3.657 const-T; Nu_H1=4.364 const-q)", "discovery_path": {"handbook": "Incropera, F.P., DeWitt, D.P., Bergman, T.L. & Lavine, A.S. (2007). Fundamentals of Heat and Mass Transfer, 6th ed. Wiley", "handbook_table": "Eq. 8.53 (Nu_D=4.36) in Table 8.4 Summary of Convection Correlations for Flow in a Circular Tube; laminar limit Re_D,c~2300 from Eq. 8.2", "handbook_page": "p.532 (Table 8.4); p.487 (Eq. 8.2 transition Reynolds)"}, "note": "Laminar-regime validity limit for the analytical fully-developed constant-Nusselt result Nu_D=4.36 (uniform/constant surface heat flux) in a circular tube: it applies while the tube Reynolds number stays laminar, up to the transition-onset value near 2300. The 4.36 is the fully-developed analytical asymptote."}}], "fluid_context": ["any single-phase fluid"], "regime_context": "laminar, hydrodynamically and thermally fully-developed flow in a circular tube with a uniform (constant) surface heat flux; the analytical fully-developed Nusselt number is the constant 4.36; valid below the laminar-to-turbulent transition (critical Reynolds ~2300)", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure", "closure-set-mismatch"], "last_reviewed": "2026-06-08", "reviewer": "seed"}
+{"closure_id": "marineau-entropy-layer-shock-interaction-2014", "closure_family": "boundary-layer-transition", "closure_name": "Marineau entropy-layer / shock-interaction hypersonic transition validity boundary", "physics_coordinates": ["entropy_layer_shock_ratio"], "validated_range": [{"coord": "entropy_layer_shock_ratio", "min": 0.1, "max": 1000000000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Marineau, E.C. et al. (2014), 'Investigation of Mach 10 Boundary-Layer Transition ...', SAND2014-4326C / AIAA; Table 3 'Transition Parameters at 0-deg AoA' (S_T/X_SW).", "primary_source": true, "note": "S_T/X_SW >= 0.1 delineates the e^N / 2nd-mode regime (benign); large-bluntness deploy runs fall to S_T/X_SW < 0.1 (entropy-layer dominated, 2nd mode absent) -> OUTSIDE (below). max=1e9 keeps the benign large-S_T/X_SW rows (up to ~113) inside."}}], "fluid_context": ["air"], "regime_context": "hypersonic boundary-layer transition on a blunt cone (Mach~10); entropy-layer/shock-interaction (bluntness)-driven; the e^N/2nd-mode framework is widely cited; transcribed directly from the primary source (Marineau Table 3)", "contested": false, "solver_availability": ["none"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-07", "reviewer": "seed"}
+{"closure_id": "mcadams-vertical-plate-natural-1954", "closure_family": "single-phase-convection", "closure_name": "McAdams vertical-plate laminar natural convection", "physics_coordinates": ["rayleigh_number", "prandtl_number"], "validated_range": [{"coord": "rayleigh_number", "min": 10000.0, "max": 1000000000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "McAdams, W. H. (1954). 'Heat Transmission' (3rd ed.). McGraw-Hill, New York, Ch. 7. Restated in: Bergman, T. L. et al. (2011). Fundamentals of Heat and Mass Transfer (7th ed.), Section 9.6, Eq. 9.27; Churchill, S. W. and Chu, H. H. S. (1975). 'Correlating equations for laminar and turbulent free convection from a vertical plate.' Int. J. Heat Mass Transfer 18(11):1323-1329 (the refined Churchill-Chu form reduces to 0.59*Ra^(1/4) in the laminar limit).", "primary_source": true, "year": 1954, "page": "Ch. 7", "note": "Simple form Nu_L = 0.59 * Ra_L^(1/4) for laminar natural convection from an isothermal vertical heated plate. Ra upper bound 1e9 is the standard laminar-to-turbulent transition for vertical plate natural convection; above this the turbulent form (Nu ~ Ra^(1/3)) applies."}}, {"coord": "prandtl_number", "min": 0.6, "max": 7.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Bergman, T. L. et al. (2011). Fundamentals of Heat and Mass Transfer (7th ed.), Section 9.6. The simple McAdams form 0.59*Ra^(1/4) is calibrated for moderate-Pr fluids (air ~0.71, water ~7); Churchill-Chu 1975 provides an explicit Pr correction (1 + (0.492/Pr)^(9/16))^(-4/9) that extends to liquid metals (Pr<<1) and other fluids. The simple form is conservatively stated for Pr in roughly [0.6, 7] in textbook restatements; for wider Pr, the Churchill-Chu corrected form is required.", "primary_source": true, "year": 2011, "page": "Section 9.6", "note": "Air at Pr=0.71 (Lance & Smith's fluid) sits well inside the validated range. The Pr range stated here is for the simple McAdams form, not the more general Churchill-Chu corrected form."}}], "fluid_context": ["air", "water", "fluid-Pr-ge-0.6-le-7"], "regime_context": "external natural convection from an isothermal vertical heated plate, laminar regime (Ra_L <= 1e9); handbook-canonical correlation widely cited across all major heat-transfer textbooks (McAdams, Bergman, Incropera, Churchill-Chu)", "contested": false, "solver_availability": ["none"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-05", "reviewer": "seed"}
+{"closure_id": "modified-sparrow-cur-asym-narrow-rect-channel-2014", "closure_family": "single-phase-convection", "closure_name": "Modified Sparrow-Cur asymmetrically-heated narrow rectangular channel Nusselt", "physics_coordinates": ["reynolds_number", "prandtl_number"], "validated_range": [{"coord": "reynolds_number", "min": 10000.0, "max": 70000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Original: Sparrow, E. M. & Cur, N. (1982). 'Turbulent Heat Transfer in a Symmetrically or Asymmetrically Heated Flat Rectangular Duct With Flow Separation at Inlet.' Journal of Heat Transfer 104:82-89. Modified form: Forrest, E. C., Hu, L.-W., Buongiorno, J., & McKrell, T. J. (2014). 'Convective Heat Transfer in a High Aspect Ratio Mini-Channel Heated on One Side.' Journal of Heat Transfer 138(2):021704, Eq. (7). DOI: 10.1115/1.4031646. SAND2014-18834J / OSTI 1295764.", "primary_source": true, "year": 2014, "note": "Modified form: Nu = 0.036 Re^0.76 Pr^(1/3). Original Sparrow-Cur (1982): Nu = 0.0464 Re^0.76 at fixed Pr=2.5, isothermal wall, alpha*=0.0556, Re 10000-45000. Forrest 2014 modification: (a) Chilton-Colburn analogy applied to extend Pr to 2.2-5.4 via explicit Pr^(1/3) term; (b) coefficient adjusted from 0.0464 to 0.036 consistent with Pr^(1/3)|_{Pr=2.5}; (c) applied to uniform-heat-flux wall BC (vs original isothermal); (d) Re range extended to 70,000 on Forrest's data. GEOMETRY-MATCHED to one-sided-heated narrow rectangular mini-channels with high aspect ratio (alpha* < 0.1); the only standard corpus closure with explicit one-sided-heating geometry dependence. Forrest 2014 Table 4 reports MAE 6.1% for this Modified form over Re 10000-70000 at Pr 2.2-5.4 -- LOWER than every other circular-pipe closure evaluated on the same one-sided narrow-rect geometry: Dittus-Boelter 6.4%, Gnielinski 8.4%, Petukhov 11.4%, Sieder-Tate 15.2%. The modified form is the most accurate published closure for this geometry."}}, {"coord": "prandtl_number", "min": 2.2, "max": 5.4, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Forrest et al. 2014 J. Heat Transfer 138(2):021704 Table 4 and Eq. (7) (Pr 2.2-5.4 tested range, water at 30-80 degC). Modification base: Sparrow & Cur 1982 J. Heat Transfer 104:82-89 (original Pr=2.5 single value, mass-transfer experiment converted to heat transfer via Chilton-Colburn analogy).", "primary_source": true, "year": 2014, "note": "Pr 2.2-5.4 is Forrest 2014's tested range, corresponding to water at bulk temperatures 30-80 degC. Original Sparrow-Cur was at Pr=2.5 only (single value). The Pr^(1/3) term in the Modified form is a Forrest-2014 Chilton-Colburn-analogy modification, not a directly validated Pr dependence in the original Sparrow-Cur experiments."}}], "fluid_context": ["water", "Pr-2-to-6"], "regime_context": "Fully developed turbulent flow in one-sided heated narrow rectangular mini-channels (high aspect ratio, alpha* < 0.1), uniform heat flux on one wall and insulated/adiabatic other walls; characteristic length = hydraulic diameter Dhyd. Modified form by Forrest 2014 extends original Sparrow & Cur 1982 (fixed Pr=2.5 isothermal mass-transfer with Chilton-Colburn) to variable Pr uniform-heat-flux via explicit Pr^(1/3) factor. Geometry-matched closure for the Forrest mini-channel substrate. Widely re-cited in narrow-channel heat-transfer validation literature for asymmetric one-sided heating (e.g., Forrest et al. 2016 J. Nucl. Eng. & Rad. Sci. on surface-oxidation effects; Tano & Marshall 2022 INL RELAP5-3D MITR/NBSR safety analyses) -- the standard closure choice for one-sided heated narrow rectangular channels in MTR fuel-channel safety analyses.", "contested": false, "solver_availability": ["none"], "weakener_link": ["extrapolated-closure", "geometry-mismatch-risk"], "last_reviewed": "2026-06-05", "reviewer": "seed"}
+{"closure_id": "pate-stainback-freestream-noise-hypersonic-1980", "closure_family": "boundary-layer-transition", "closure_name": "Pate-Stainback freestream-noise hypersonic transition validity boundary", "physics_coordinates": ["freestream_noise_rms_pitot_pct"], "validated_range": [{"coord": "freestream_noise_rms_pitot_pct", "min": 0.5, "max": 10.0, "unit": "percent", "bound_status": "seed", "provenance": {"citation": "Pate, S.R. & Stainback, P.C. \u2014 freestream-disturbance effects on boundary-layer transition in hypersonic wind tunnels (conventional-tunnel transition correlations validated for RMS-Pitot >= ~0.5%; flight/quiet ~0.05-0.1%). Casper et al. quiet-vs-noisy BAM6QT/HWT data.", "primary_source": false, "note": "0.5% RMS-Pitot is the conventionally-cited lower edge of the noisy-tunnel regime the transition correlation was built on; the flight-like quiet regime (~0.05%) is OUTSIDE it (below). max=10% keeps the noisy training (1.0-3.7%) inside -> do-no-harm."}}], "fluid_context": ["air"], "regime_context": "hypersonic boundary-layer transition on a sharp cone; disturbance(freestream-noise)-driven; Pate-Stainback freestream-noise correlation is widely cited in the hypersonic-transition literature", "contested": false, "solver_availability": ["none"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-07", "reviewer": "seed"}
+{"closure_id": "petukhov-1970", "closure_family": "single-phase-convection", "closure_name": "Petukhov forced-convection Nusselt", "physics_coordinates": ["reynolds_number", "prandtl_number"], "validated_range": [{"coord": "reynolds_number", "min": 10000.0, "max": 5000000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Petukhov, B.S. (1970). 'Heat Transfer and Friction in Turbulent Pipe Flow with Variable Physical Properties.' Advances in Heat Transfer 6:503-564. Re range 10000-5e6 per Bergman, T. L., Lavine, A. S., Incropera, F. P., DeWitt, D. P. (2011). Fundamentals of Heat and Mass Transfer (7th ed.), Section 8.5 restatement of Petukhov's stated validity range.", "primary_source": false, "year": 1970, "note": "High-accuracy turbulent-pipe correlation; superseded for transition coverage by Gnielinski 1976. Status 'claimed' (not 'confirmed') because the primary paper (Petukhov 1970 in Advances in Heat Transfer 6:503-564) has not been directly verified; bounds are cited via the Bergman/Incropera 2011 restatement chain."}}, {"coord": "prandtl_number", "min": 0.5, "max": 2000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Petukhov, B.S. (1970). Advances in Heat Transfer 6:503-564. Pr range 0.5-2000 per Bergman/Incropera 2011 sec. 8.5 restatement.", "primary_source": false, "year": 1970, "note": "Same Pr range as Gnielinski 1976 extended-coverage form."}}], "fluid_context": ["water", "air", "oils"], "regime_context": "fully developed turbulent flow in smooth tubes; high-accuracy correlation form using the friction factor as a free parameter; superseded for transition (Re 2300-10000) coverage by Gnielinski 1976.", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-05", "reviewer": "seed"}
+{"closure_id": "sieder-tate-1936", "closure_family": "single-phase-convection", "closure_name": "Sieder-Tate forced-convection with property-variation correction", "physics_coordinates": ["reynolds_number", "prandtl_number"], "validated_range": [{"coord": "reynolds_number", "min": 10000.0, "max": 1000000.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Sieder, E.N. & Tate, G.E. (1936). 'Heat transfer and pressure drop of liquids in tubes.' Ind. Eng. Chem. 28(12):1429-1435", "doi": "10.1021/ie50324a027", "primary_source": false, "note": "Sieder-Tate adds (\u03bc/\u03bc_w)^0.14 viscosity-ratio correction; ranges per handbook restatement"}}, {"coord": "prandtl_number", "min": 0.7, "max": 16700.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"citation": "Sieder-Tate 1936 IEC", "doi": "10.1021/ie50324a027", "primary_source": false, "note": "much wider Pr range than Dittus-Boelter due to viscosity correction; primary-source verification pending"}}], "fluid_context": ["water", "oils", "high-viscosity liquids"], "regime_context": "turbulent flow in tubes with significant property variation (large wall-bulk temperature difference)", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure"], "last_reviewed": "2026-06-03", "reviewer": "seed"}
+{"closure_id": "sieder-tate-laminar-combined-entry-1936", "closure_family": "single-phase-convection", "closure_name": "Sieder-Tate laminar combined-entry-length Nusselt (property-variation correction)", "physics_coordinates": ["prandtl_number", "viscosity_ratio_wall_bulk"], "validated_range": [{"coord": "prandtl_number", "min": 0.6, "max": 5.0, "unit": "dimensionless", "bound_status": "seed", "provenance": {"primary_source": false, "citation": "Sieder, E.N. & Tate, G.E. (1936). 'Heat transfer and pressure drop of liquids in tubes.' Ind. Eng. Chem. 28(12):1429-1435", "doi": "10.1021/ie50324a027", "year": 1936, "discovery_path": {"handbook": "Incropera, F.P., DeWitt, D.P., Bergman, T.L. & Lavine, A.S. (2007). Fundamentals of Heat and Mass Transfer, 6th ed. Wiley", "handbook_table": "Eq. 8.57, \u00a78.4.2 The Entry Region (laminar combined-entry-length correlation)", "handbook_page": 513}, "note": "Prandtl applicability of the laminar combined-entry Sieder-Tate correlation as restated by Incropera."}}, {"coord": "viscosity_ratio_wall_bulk", "min": 0.0044, "max": 9.75, "unit": "dimensionless", "bound_status": "seed", "provenance": {"primary_source": false, "citation": "Sieder, E.N. & Tate, G.E. (1936). 'Heat transfer and pressure drop of liquids in tubes.' Ind. Eng. Chem. 28(12):1429-1435", "doi": "10.1021/ie50324a027", "year": 1936, "discovery_path": {"handbook": "Incropera, F.P., DeWitt, D.P., Bergman, T.L. & Lavine, A.S. (2007). Fundamentals of Heat and Mass Transfer, 6th ed. Wiley", "handbook_table": "Eq. 8.57, \u00a78.4.2 The Entry Region (laminar combined-entry-length Sieder-Tate)", "handbook_page": 513}, "note": "Validity envelope of the Sieder-Tate viscosity-ratio correction for the laminar combined-entry tube correlation: wall-to-bulk viscosity ratio spanning 0.0044 up to 9.75, with Prandtl between 0.60 and 5; wall viscosity is taken at the surface temperature and the remaining properties at the mean."}}], "fluid_context": ["oils", "water", "liquids with large property variation"], "regime_context": "laminar combined thermal+velocity entry length in a circular tube; large wall-to-bulk property variation handled by the Sieder-Tate viscosity-ratio correction", "contested": false, "solver_availability": ["openfoam", "fluent", "starccm+", "cfx"], "weakener_link": ["extrapolated-closure", "closure-set-mismatch"], "last_reviewed": "2026-06-07", "reviewer": "seed"}