in silico
Phase 1: In Silico Validation
- Objective
- To determine if the EnKF algorithm, when coupled to a simplified but representative TRC model, can achieve state estimation errors below 15% under simulated sensor noise and model errors, and to assess sensitivity to key unknown parameters.
- Estimated cost
- €0-2,000 (primarily software licenses if not available; OpenFOAM is open-source)
- Estimated duration
- 4-8 weeks
- Success criteria
- RMSE of estimated mushy zone thickness · ≤ 15% (relative to simulated “truth”) for ≥80% of the tested operational points (casting speed, roll gap). · (Calculated from the EnKF state vector (solid fraction field) vs. the “truth” simulation output at the same timestep. Thickness defined as distance between 0.2 and 0.8 solid fraction isotherms at the strip centerline.)
- Sensitivity to IHTC error · A ±30% perturbation in the model’s IHTC causes a < ±25% change in the EnKF-estimated average cooling rate. · (Compare the distribution of estimated cooling rates from Monte Carlo runs with nominal vs. perturbed IHTC, using the paired difference analysis from the methodology.)
- Filter stability · EnKF does not diverge (error remains bounded) over a simulated casting duration of 60 seconds for all tested tuning parameters within a reasonable range. · (Visual inspection and analysis of estimation error time series.)
- Go if
- Success Criteria 1 AND 2 are met, demonstrating the core estimation concept is viable in silico and robust to a key model error.
- No-go if
- Success Criterion 1 is not met (RMSE >15%) even under ideal simulated conditions (low noise, perfect model), indicating a fundamental flaw in the EnKF application or model observability.
- Pivot if
- Success Criterion 1 is met but Criterion 2 is not (high sensitivity to IHTC). Pivot to Phase 2 focusing on experimental IHTC characterization or to an EnKF variant more robust to parametric errors (e.g., dual estimation).
- Risks
- Simplified 2D model fails to capture essential 3D effects (edge effects, roll deflection) leading to misleadingly good/bad in silico performance.Probability: mediumMitigation: Include a 3D effect proxy in the 2D model (e.g., adjusted boundary conditions based on literature) and plan for a 3D model check in Phase 1 if results are borderline.
- EnKF tuning (ensemble size, inflation) is overly sensitive and difficult to stabilize, consuming excessive time.Probability: mediumMitigation: Start with established tuning heuristics from literature (e.g., ensemble size = 50-100, adaptive inflation). Use open-source DA libraries with built-in utilities.
- Lack of access to commercial simulation software (COMSOL).Probability: lowMitigation: Use open-source OpenFOAM. The “chtMultiRegionFoam” solver with melting/solidification modifications is a viable, though more complex, starting point.