Multifluid Hydrodynamic Simulation of Metallic-Plate Collision Using the VOF Method

This study presents a one-dimensional multifluid hydrodynamic simulation of lead and steel plate collisions in explosive welding using a Godunov-type algorithm with pressure relaxation and the Volume-of-Fluid method, successfully capturing tensile stresses and validating the unloading wave arrival time against experimental data.

The Gist

Imagine you are watching a high-speed movie of two metal plates slamming into each other. One is made of lead (soft and heavy), and the other is steel (hard and strong). They are moving so fast—500 meters per second—that when they hit, they don't just bounce; they squish, heat up, and behave like a thick, super-fast fluid for a split second. This is the science of explosive welding, a technique used to join different metals together.

The paper you shared is about building a super-accurate computer simulation to watch this crash happen, frame by frame, without actually blowing anything up in a lab.

Here is the story of how they did it, explained simply:

1. The Problem: Tracking the "Invisible" Line

When the lead and steel hit, they mix at the boundary. In a computer, this is tricky. Imagine a grid of boxes (like a chessboard) covering the scene.

• The Issue: Most computer programs treat a box as being filled with either lead or steel. But right at the collision point, a single box might be half lead and half steel.
• The Old Way: Older methods would blur this line, making the lead look like it's slowly turning into steel, or vice versa. It's like trying to paint a sharp line between red and blue paint, but your brush is too big, so you end up with purple in the middle.
• The New Way (VOF Method): The authors used a technique called VOF (Volume of Fluid). Think of this like a very smart, digital "liquid tracker." Instead of guessing what's in the box, it tracks exactly how much of the box is lead and how much is steel. It keeps the line between them razor-sharp, just like a high-definition camera.

2. The Physics: The "Squeezed Sponge" Analogy

Metals are usually thought of as solid, but when hit this hard, they act like fluids. However, they are "squishy" fluids.

• Compressibility: Imagine a sponge. If you squeeze it, it gets smaller. Lead is like a soft sponge; steel is like a stiff sponge. Air is like a very loose, airy sponge.
• The Challenge: When the plates collide, the "sponges" get squeezed. The air compresses easily, but the steel resists. The computer has to calculate how the pressure changes in each material simultaneously.
• The Solution: The authors created a set of rules (mathematical equations) that treat the whole mix as a single "effective" fluid for movement, but keeps the individual "sponge" properties (density, pressure) separate for each material. This allows them to see exactly how the lead squishes differently than the steel.

3. The "Negative Pressure" Magic

This is the coolest part of the paper.

• The Scenario: When the shockwave from the collision hits the back of the lead plate, it bounces back. Imagine pulling a rubber band too fast; it snaps back. In physics, this creates tension (pulling apart) rather than compression (pushing together).
• The Problem: In many computer models, "negative pressure" (tension) causes the math to crash or break. It's like a calculator trying to divide by zero.
• The Breakthrough: This new simulation handles negative pressure perfectly. It allows the computer to say, "Hey, the metal is being pulled apart right now," without breaking the program. This is crucial because it helps scientists understand why the metals might crack or weld together in specific ways.

4. The Result: A Perfect Match

The team ran their simulation and compared it to:

1. Real-world experiments (actual explosive welding tests).
2. Other computer models (older, less accurate methods).

The Verdict: Their method was the winner.

• It was sharper: It didn't blur the line between lead and steel.
• It was faster: It stabilized the results quickly.
• It was accurate: It predicted exactly when the "unloading wave" (the bounce-back effect) would hit the interface, matching real-life experiments within a tiny fraction of a microsecond.

The Big Picture Analogy

Think of the old simulation methods as a low-resolution pixelated video of a car crash. You can see the cars hit, but the details are blurry, and you can't tell exactly when the bumper bends.

This new paper presents a 4K, slow-motion video of the crash. It shows exactly how the metal flows, how the pressure builds, and even how the metal stretches (negative pressure) before it settles.

Why does this matter?
By understanding exactly how these metals behave at the moment of impact, engineers can design better ways to weld dissimilar metals (like joining aluminum to steel for electric cars or aerospace parts) without wasting expensive materials or risking failure. They can "simulate the crash" on a computer to find the perfect recipe for a strong bond before ever lighting a fuse.

Technical Summary

1. Problem Statement

The study addresses the numerical simulation of explosive welding, specifically the high-speed collision of dissimilar metal plates (lead and steel). The goal is to model the complex dynamic effects occurring during the collision, where metals transition into a "pseudo-fluid" condensed state for several microseconds.

• Physical Setup: A one-dimensional problem involving a lead plate (2 mm thick, 11,300 kg/m³) impacting a steel plate (3 mm thick, 7,900 kg/m³) at a relative velocity of 500 m/s. Both are surrounded by air layers.
• Key Challenge: Accurately tracking the interface between immiscible phases (lead, steel, air) while accounting for compressibility, different equations of state (EOS), and the emergence of tensile stresses (negative pressures) which often cause numerical instability in standard hydrodynamic codes.
• Objective: To develop a robust multifluid algorithm that captures wave propagation, interface dynamics, and negative pressure regions with low numerical diffusion, validated against experimental data and previous simulations.

2. Methodology

The authors propose a modified multifluid Godunov-type (finite-volume) algorithm based on the mechanical-equilibrium Euler equations.

A. Mathematical Model

• Effective Phase Approach: The mixed cell is treated as a single "effective" phase with shared velocity ($u$) and pressure ($p$), but distinct densities and internal energies for each constituent phase.
• Governing Equations:
  • Conservation Laws: Standard conservation of mass, momentum, and total energy are applied to the effective phase.
  • Volume Fraction Advection: A specific advection equation is derived for the volume fraction ($f^{(\alpha)}$) of each phase. Crucially, this equation accounts for compressibility differences between phases. It ensures that more compressible fluids change volume more readily during pressure equilibration, derived using isentropic bulk moduli ($K_S$).
  • Energy Equations: Separate energy equations are solved for each phase to accommodate different Equations of State (EOS). The internal energy equation includes a source term representing the work done by pressure during phase expansion/compression.
• Equations of State: A stiffened-gas EOS is used for all phases (lead, steel, air), allowing for the modeling of condensed media under high pressure.
• Pressure Relaxation: The model assumes fast pressure equilibration (mechanical equilibrium) at the interface. An iterative relaxation step is performed at each time step to redistribute volume fractions and energies until a common pressure is achieved within mixed cells.

B. Numerical Scheme

• Interface Tracking: The Volume-of-Fluid (VOF) method is employed to track the sharp interfaces between immiscible phases. It uses a donor-acceptor scheme for volume advection and material-interface reconstruction.
• Riemann Solver: The HLLC (Harten–Lax–van Leer–Contact) solver is used to resolve local Riemann problems at cell interfaces, providing high resolution for contact discontinuities.
• Reconstruction: The authors compare results using MUSCL-Hancock reconstruction (second-order accuracy) against first-order upwind schemes. MUSCL significantly reduces numerical diffusion.
• Handling Negative Pressure: A key feature of the algorithm is its ability to handle tensile stresses (negative pressure) naturally without special stabilization procedures, which is critical for modeling the unloading waves in the lead plate.

3. Key Contributions

1. Modified CGF Algorithm: The authors adapt the Colella-Glaz-Ferguson (CGF) method (originally for gases) for condensed media by incorporating compressibility effects into the volume fraction advection equation and using stiffened-gas EOS.
2. Robustness with Tensile Stresses: The method successfully simulates the formation of negative pressures and density dips in the lead plate without numerical breakdown, a common failure point in other multiphase solvers.
3. Low Numerical Diffusion: By combining the VOF method with the HLLC solver and MUSCL reconstruction, the scheme achieves high resolution of the contact interface, outperforming previous methods (like the Baer–Nunziato model) on similar grid resolutions.
4. Validation: The model is rigorously validated against experimental data regarding the arrival time of unloading waves and compared against high-resolution simulations from literature.

4. Results

The simulation was run up to $t = 2$ $\mu$s, covering the collision, shock wave propagation, and unloading wave interactions.

• Wave Dynamics:
  • Shock waves formed at the impact interface and propagated to the free surfaces.
  • Unloading waves reflected from the free surfaces (air-metal interfaces) and traveled back toward the collision interface.
  • The unloading wave from the steel reached the interface at $t \approx 1.13$ $\mu$s, consistent with experimental data and previous numerical results ($1.05$ $\mu$s).
• Interface Velocity:
  • The contact velocity increased from 191 m/s to 425.2 m/s as the unloading wave passed through.
  • The numerical solution with MUSCL reconstruction on a coarse grid (3,601 cells) agreed with a fine-grid reference (10,000 cells) within 4.6%, demonstrating high efficiency.
• Negative Pressure:
  • At $t \approx 1.4$ $\mu$s, intersecting unloading waves in the lead plate created regions of negative pressure and reduced density.
  • The proposed method captured these features accurately, whereas previous methods (e.g., Baer–Nunziato) struggled with stability in this regime.
• Conservation: The total error in effective density and integral total energy was negligible (0.06% and 0.02%, respectively) over the simulation duration.

5. Significance

This work demonstrates that a multifluid Eulerian approach with VOF interface tracking is highly effective for simulating explosive welding processes.

• Scientific Impact: It provides a validated numerical tool for studying the Rayleigh–Taylor instability and pseudo-fluid states in metal collisions, which are difficult to observe experimentally.
• Engineering Application: The ability to accurately model negative pressures and interface dynamics without excessive numerical diffusion is crucial for optimizing explosive welding parameters and predicting joint quality.
• Future Work: The authors plan to extend this 1D methodology to 2D formulations to further investigate the instability mechanisms at the contact interface during welding.

In summary, the paper presents a robust, low-diffusion numerical framework that successfully bridges the gap between theoretical multiphase fluid dynamics and the complex realities of high-speed metal collision, particularly in handling the challenging physics of tensile stresses in condensed matter.