From Coils to Surface Recession: Fully Coupled Simulation of Ablation in ICP Wind Tunnels

This paper presents a fully coupled, multiphysics computational framework that successfully simulates the end-to-end ablation process in the UIUC Plasmatron X ICP wind tunnel, demonstrating high predictive accuracy for heat flux, temperature, and recession rates when validated against experimental data.

The Gist

The Big Picture: Simulating a "Super-Heated Hair Dryer"

Imagine you are designing a spaceship that needs to fly so fast it breaks the sound barrier. When it does, the air in front of it gets crushed and heated to temperatures hotter than the surface of the sun. To survive, the spaceship needs a special "thermal suit" (called a Thermal Protection System, or TPS).

Testing these suits in real flight is incredibly expensive and dangerous. So, scientists use giant machines on the ground called ICP Wind Tunnels. Think of these tunnels as super-powered, electrode-free hair dryers that shoot a stream of super-hot, electrically charged gas (plasma) at a sample of material to see how long it lasts before it melts or burns away.

This paper is about a new, incredibly smart computer program that can simulate exactly what happens inside this "hair dryer" and how the material reacts, without needing to run a single physical experiment first.

---

The Problem: The Old Way Was Like Guessing

Previously, scientists tried to simulate these tests in two separate steps:

1. Step A: They would guess how hot the air was at the surface of the material.
2. Step B: They would feed that guess into a second program to see how the material melted.

The Analogy: Imagine trying to predict how a piece of butter melts on a hot pan.

• The Old Way: You guess the pan is 200°F. You calculate the melting. Then you realize the butter actually cooled the pan down, so you guess again. It's a slow, back-and-forth guessing game that often gets the details wrong.
• The Problem: As the material melts (ablates), it releases gas that pushes back against the hot air, changing the temperature and pressure. The old methods couldn't handle this "feedback loop" well.

The Solution: The "Fully Coupled" Team-Up

The authors built a fully coupled simulation. This means they created a single digital ecosystem where three different physics problems talk to each other instantly and continuously.

Think of it like a three-person relay team passing a baton at the speed of light:

1. The Flow Solver (HEGEL): This is the "Weatherman." It simulates the swirling, super-hot plasma gas, the wind speed, and the pressure.
2. The Magnet Solver (FLUX): This is the "Electrician." It simulates the giant copper coils that create the magnetic field to heat the gas. It tells the Weatherman where the heat is coming from.
3. The Material Solver (CHyPS): This is the "Survivor." It simulates the graphite sample. It feels the heat, gets hot, starts to melt (ablate), and sends a signal back saying, "Hey, I'm melting and releasing gas, which is changing the wind!"

The Magic: In this new framework, when the "Survivor" melts, the "Weatherman" immediately knows the wind has changed, and the "Electrician" adjusts the power. They all update each other thousands of times a second. This is called ab initio, meaning they start from the basic laws of physics and build up, rather than relying on guesses or experimental data to fill in the blanks.

What Did They Do?

They took their new super-simulation and applied it to a real machine at the University of Illinois called Plasmatron X.

• The Setup: They simulated a 350 kW machine (that's a lot of power, like 350 hair dryers running at once) shooting plasma at a 30mm puck of graphite.
• The Test: They ran the simulation for 200 seconds (a long time for a computer simulation) to see how the graphite puck would shrink (recede) and how hot it would get.

The Results: "Spot On"

When they compared their computer predictions to real-world experiments done in the actual lab:

• Temperature: The computer predicted the temperature of the melting puck with less than 12% error. That is incredibly accurate for something this complex.
• Melting Rate: The computer predicted how fast the puck would shrink with less than 10% error.
• The "Hair Dryer" Transition: They also showed their model could predict when the plasma jet would switch from a gentle breeze (subsonic) to a supersonic shockwave, just like in real life.

Why Does This Matter?

This is a game-changer for aerospace engineering.

1. No More Guessing: Before, if a simulation didn't match the experiment, engineers would tweak the numbers until it did. This new method doesn't need that. It predicts the outcome based purely on physics.
2. Save Money and Time: Instead of building expensive physical samples and burning them up in a wind tunnel, engineers can now run these "digital twins" to test new materials. They can design better heat shields faster and cheaper.
3. Safety: By understanding exactly how materials behave in these extreme conditions, we can design spacecraft that are safer for astronauts.

The Bottom Line

The authors have built a digital time machine for hypersonic testing. They successfully connected the dots between the electricity, the hot gas, and the melting material into one seamless story. While there are still tiny imperfections (mostly because real-world materials are messy and hard to measure perfectly), this framework proves that we can now simulate these extreme environments with high confidence, paving the way for the next generation of spacecraft.

Technical Summary

1. Problem Statement

Hypersonic vehicles face extreme aerodynamic heating and chemical reactions during flight, necessitating robust Thermal Protection Systems (TPS). Ground-based testing in Inductively Coupled Plasma (ICP) wind tunnels is critical for evaluating TPS materials because these facilities generate high-enthalpy, electrode-free plasma flows that mimic flight conditions without contaminating the test sample.

However, accurately predicting material response in ICP tunnels is challenging due to:

• Complex Physics: The interaction between electromagnetic fields, plasma flow (Navier-Stokes), and thermo-chemical material ablation.
• Limitations of Current Methods: Traditional "film coefficient" approaches decouple the flow and material, failing to capture non-equilibrium effects, surface blowing, and time-varying boundary conditions caused by geometric changes (recession).
• Lack of Predictive Capability: Existing simulations often rely on experimentally measured boundary conditions (e.g., nozzle exit profiles) or empirical tuning of power coupling efficiency, rather than simulating the facility from the power source to the material surface.

The objective of this work is to develop and validate a fully coupled, end-to-end numerical framework that simulates an entire ICP facility—from electromagnetic induction to plasma generation and finally to time-accurate material ablation—without relying on empirical tuning or experimental input data for boundary conditions.

2. Methodology

The authors developed a partitioned, multi-physics computational framework integrating three distinct solvers:

• Plasma Flow Solver (HEGEL): A high-fidelity finite-volume Navier-Stokes solver. It handles the plasma flow using a two-temperature (2T) Non-Local Thermodynamic Equilibrium (NLTE) model for low-pressure cases and Local Thermodynamic Equilibrium (LTE) for high-pressure cases.
• Electromagnetic Solver (FLUX): A mixed finite-element solver that computes the electromagnetic field generated by the induction coils. It solves Maxwell's equations to determine the induced currents and fields.
• Material Response Solver (CHyPS): A high-order discontinuous Galerkin solver based on the mfem library. It models the thermo-chemical response of the TPS, including heat conduction, pyrolysis, and surface recession (ablation) using a $B'$-table approach for equilibrium chemistry.

Coupling Strategy:
The solvers are coupled via preCICE, an open-source library for partitioned multi-physics simulations.

1. Volume Coupling: HEGEL and FLUX exchange data (Joule heating, Lorentz force, electrical conductivity) throughout the domain.
2. Surface Coupling: A loosely coupled, pseudo-steady-state approach is used between the plasma and material solvers.
   • HEGEL provides heat flux ($q''$) and pressure ($p$) to CHyPS.
   • CHyPS calculates surface temperature ($T$), recession (displacement $\Delta$), and blowing rate ($\dot{m}''$).
   • HEGEL updates the mesh geometry based on $\Delta$ (using Linear Transfinite Interpolation) and applies a blowing wall boundary condition based on surface normal momentum balance.
   • This cycle repeats at dynamic time intervals based on the temperature gradient.

Simulation Setup:
The framework was applied to the 350 kW Plasmatron X facility at the University of Illinois Urbana-Champaign (UIUC). The simulation covers the induction coils, torch, nozzle, chamber, and a 30 mm isoQ graphite TPS sample. The gas mixture is an 11-species air model.

3. Key Contributions

• First Ab Initio, End-to-End Simulation: This is the first study to simulate an entire ICP facility from the electromagnetic power deposition in the coils to the material recession of the test sample, driven solely by independently measured facility inputs (power, pressure, efficiency) without empirical calibration.
• Fully Coupled Multiphysics Framework: It successfully integrates electromagnetic induction, magnetohydrodynamics (MHD), and time-accurate material ablation in a single predictive loop.
• Validation of Low-Pressure Transitions: The model captures the transition from subsonic to supersonic plasma jet behavior at low pressures, a phenomenon driven by Rayleigh flow choking, validating the framework's ability to handle complex flow regimes.
• Elimination of Empirical Tuning: Unlike previous studies that adjust coupling efficiency or effective power to match experiments, this framework uses physics-based models exclusively, establishing a truly predictive capability.

4. Results

The framework was validated against cold-wall calorimetry and graphite ablation experiments in the Plasmatron X facility under three different operating conditions (varying pressure and efficiency).

• Plasma Physics Validation:

  • The model correctly reproduced vortex-mode recirculation, Joule heating, and Lorentz-force-induced flow confinement.
  • It accurately predicted the transition from subsonic to supersonic flow at low pressures (600 Pa), matching experimental shock-diamond structures and spacing within ~10–20% error.

• Cold-Wall Heat Flux:

  • Predicted stagnation-point heat fluxes fell within the 20% experimental uncertainty bounds for all tested cases. The authors note that experimental heat flux measurements have inherent uncertainties (e.g., slug diameter, thermal losses), making this level of agreement significant.

• Material Response (Ablation):

  • Temperature History: The fully coupled simulations reproduced measured stagnation-point temperature histories with high accuracy.
    • Case 1 (20 kPa): 3.3% error.
    • Case 2 (10 kPa): 5.9% error.
    • Case 3 (5.5 kPa): 11.2% error.
  • Recession Rates: The predicted steady-state recession rates matched experimental data with errors below 10% for all cases (4.6% to 9.8%).
  • Transient Behavior: Minor discrepancies were observed during the early transient heating phase (first few seconds), attributed to uncertainties in power-coupling efficiency, simplified equilibrium ablation models, and the use of material property data for a slightly different graphite grade (Mersen 2160 vs. 2340).

5. Significance

This work represents a paradigm shift in the simulation of hypersonic ground testing facilities. By demonstrating that a fully coupled, physics-based model can predict complex ablation phenomena with high fidelity without empirical tuning, the authors have established a robust foundation for predictive design.

The implications include:

• Improved Mission Planning: More reliable interpretation of experimental data and better planning for hypersonic material testing campaigns.
• Design Optimization: The ability to simulate "what-if" scenarios for new TPS materials and facility configurations without costly physical prototyping.
• Scientific Advancement: It bridges the gap between electromagnetic facility modeling and material science, providing a comprehensive tool to understand the intricate gas-surface interactions in high-enthalpy environments.

The authors conclude that while minor uncertainties remain (particularly in transient phases and material property datasets), the framework's predictive capability is sufficient to serve as a primary tool for the design and analysis of future hypersonic missions.