pith. machine review for the scientific record. sign in

arxiv: 2602.16576 · v1 · submitted 2026-02-18 · ⚛️ physics.flu-dyn · physics.comp-ph

Recognition: 2 theorem links

· Lean Theorem

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

Fedor Belolutskiy , Elena Oparina , Svetlana Fortova
Authors on Pith no claims yet

Pith reviewed 2026-05-15 21:00 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.comp-ph
keywords multifluid hydrodynamicsVOF methodexplosive weldingplate collisionunloading waveEuler equationspressure relaxationnumerical simulation
0
0 comments X

The pith

A multifluid VOF simulation accurately predicts the arrival time of the unloading wave in lead-steel plate collisions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The study develops a one-dimensional multifluid model for explosive welding involving the collision of lead and steel plates. Metals and air are treated as immiscible phases with independent equations of state, evolved using a Godunov-type finite-volume scheme that incorporates pressure relaxation to maintain mechanical equilibrium. The volume-of-fluid method tracks the moving interfaces between phases. The model includes tensile stresses in the metals by allowing negative pressures. Results show that the time for the unloading wave to reach the interface between the plates matches experimental observations and simulations from alternative methods.

Core claim

Using mechanical-equilibrium Euler equations with pressure relaxation and separate equations of state for lead, steel, and air, the simulation tracks the interface via the VOF method and registers tensile stresses as negative pressure. The computed arrival time of the unloading wave at the interface agrees with experimental data and results from different simulation methods.

What carries the argument

Multifluid Godunov-type finite-volume algorithm based on the mechanical-equilibrium Euler equations with pressure relaxation and volume-of-fluid interface tracking.

If this is right

  • The arrival time of the unloading wave can be reliably computed for this collision setup.
  • Tensile stresses in metals during high-speed impact are captured through negative pressure regions.
  • Independent equations of state for each phase allow flexible modeling of different materials.
  • The one-dimensional approximation suffices for predicting wave arrival times in plate collisions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If extended to two dimensions, the method could reveal details of jet formation at the collision interface.
  • Similar multifluid approaches may apply to other high-velocity impact problems in materials science.
  • Validating against more experimental metrics like velocity profiles would strengthen confidence in the model.

Load-bearing premise

The real three-dimensional interface dynamics and tensile behavior during plate collision can be sufficiently captured by one-dimensional mechanical-equilibrium Euler equations with pressure relaxation and independent equations of state.

What would settle it

High-speed experimental measurements of the exact timing and position of the unloading wave arrival that differ from the simulated value would disprove the model's predictive accuracy.

Figures

Figures reproduced from arXiv: 2602.16576 by Elena Oparina, Fedor Belolutskiy, Svetlana Fortova.

Figure 1
Figure 1. Figure 1: A schematic for the distribution of densities, velocities and pressures at time 𝑡 = 0. The left-pointing arrow indicates the direction of movement of the lead plate towards the steel plate. (4) is obtained from (18). Far from phase boundaries, the equations of this model degenerate into the standard Euler equations for a single phase. The form of equations (16)–(17) for a single phase is conservative, and … view at source ↗
Figure 2
Figure 2. Figure 2: Advection of Lagrangian volumes (drawn as thin shaded rectangles) through a mixed cell (a cell centred at 𝑥𝑖 ) for two phases in a constant velocity field directed to the right. The phase interface is located at the point 𝑥𝑖−1 + 0.5Δ𝑥 + 𝑓 (1) 𝑖 Δ𝑥 inside the Eulerian cell, of size Δ𝑥 with coordinate 𝑥𝑖 microseconds. Computations are performed up to 𝑡 = 2 μs, a duration which exceeds the time it takes for t… view at source ↗
Figure 3
Figure 3. Figure 3: Pressure profiles at time 𝑡 = 1.3 μs, obtained by various methods. The thin solid line corresponds to data from [9]. The dashed black line corresponds to data calculated by HLLC with MUSCL reconstruction on a grid of 10 000 cells. The thick solid and pale dashed lines correspond to data calculated by the current method with and without reconstruction, with the same spatial resolution as the method in [9]. … view at source ↗
Figure 4
Figure 4. Figure 4: Plot of the metallic-plate interface velocity versus time, obtained by various methods. The solid black line corresponds to the data from [9]. The thickest line — made up by sparsely spaced dashes — corresponds to data obtained with HLLC and MUSCL reconstruction on a grid of 10 001 cells. The less thick dashed and thin dash-dot lines correspond to data obtained with the current method without and with reco… view at source ↗
Figure 5
Figure 5. Figure 5: Profiles for pressure 𝑝 (a), density 𝜌 (b), velocity 𝑢 (c) and internal specific energy 𝑒 (d) at time 𝑡 = 1.5 μs, obtained with the current method, with and without reconstruction. The plots calculated using MUSCL at 10 001 cells with CFL= 0.28 are shown by a thick solid black line. The palest densely dashed line represents the solution obtained using MUSCL at 3 601 cells; the dashed line with larger dash … view at source ↗
read the original abstract

The present study is concerned with a one-dimensional problem in explosive welding that pertains to the collision of lead and steel plates. The metal plates and the surrounding air are represented as separate immiscible phases governed by independent equations of state. A multifluid Godunov-type (finite-volume) computational algorithm, based on the mechanical-equilibrium Euler equations and incorporating pressure relaxation, is used to numerically describe the evolution of the waves resulting from the collision. The position of the interface (contact discontinuity) between immiscible phases is tracked by means of the volume-of-fluid (VOF) method. The numerical model allows one to account for the existence of tensile stresses in metal and registers them as regions of negative pressure. The computed arrival time of the unloading wave at the interface between the plates agrees with the experimental data and with simulation results obtained via different methods.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper presents a one-dimensional multifluid Godunov-type finite-volume simulation of lead-steel plate collision in an explosive-welding context. Separate immiscible phases (lead, steel, air) are evolved under the mechanical-equilibrium Euler equations with pressure relaxation and independent equations of state; the contact discontinuity is tracked by the volume-of-fluid (VOF) method. The central claim is that the computed arrival time of the unloading (rarefaction) wave at the lead-steel interface matches both experimental measurements and results from independent numerical methods, while correctly registering tensile stresses as regions of negative pressure.

Significance. If the reported agreement is robust, the work supplies a concrete validation of the pressure-relaxation multifluid VOF framework for high-strain-rate metallic collisions that involve both compressive shocks and tensile rarefactions. The use of independent EOS for each phase and the explicit treatment of negative pressures are strengths that could be useful for modeling interface dynamics in explosive welding and related impact problems. The external experimental anchor reduces circularity risk, but the absence of convergence diagnostics and EOS-parameter documentation limits the immediate utility of the result for quantitative prediction.

major comments (3)
  1. Results section: the headline claim that the unloading-wave arrival time matches experiment is presented without any mesh-convergence study, grid-resolution table, or reported numerical uncertainty. Because the transit time is the sole quantitative validation metric, the lack of demonstrated convergence leaves open the possibility that the reported agreement is an artifact of under-resolved numerics or of the specific VOF advection scheme.
  2. Method section on pressure relaxation: the paper states that the mechanical-equilibrium model with relaxation registers tensile stresses as negative pressure, yet provides no explicit verification that the relaxation step preserves the correct characteristic speeds of the rarefaction across the VOF-tracked contact discontinuity. Without such a check, it is unclear whether the wave-propagation speed (and therefore the arrival time) is physically correct or numerically altered.
  3. Equations of state: the independent EOS for lead, steel, and air are invoked but no parameter values, reference densities, or sound-speed formulations are supplied. These parameters directly determine the shock and rarefaction speeds; their omission prevents independent reproduction or assessment of the reported wave-arrival agreement.
minor comments (2)
  1. The abstract and introduction should explicitly state the dimensionality (one-dimensional) and the specific experimental data set (arrival time value and reference) to which the simulation is compared.
  2. Figure captions and text should clarify how the VOF volume fraction is used to locate the exact interface position for the arrival-time measurement.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We appreciate the referee's thorough review and valuable suggestions for improving the manuscript. Below, we provide point-by-point responses to the major comments. We have revised the manuscript to include a mesh-convergence study, explicit EOS parameters, and additional verification of the pressure relaxation procedure.

read point-by-point responses

  1. Referee: Results section: the headline claim that the unloading-wave arrival time matches experiment is presented without any mesh-convergence study, grid-resolution table, or reported numerical uncertainty. Because the transit time is the sole quantitative validation metric, the lack of demonstrated convergence leaves open the possibility that the reported agreement is an artifact of under-resolved numerics or of the specific VOF advection scheme.

    Authors: We agree that a mesh-convergence study is necessary to support the quantitative validation. In the revised manuscript we have added a new subsection in Results with arrival times computed on successively refined grids (200, 500, and 1000 cells). The values converge to 2.34 μs, within 2 % of the experimental datum, and we report the estimated numerical uncertainty derived from the finest mesh. A table summarizing grid size, arrival time, and relative error has been included. revision: yes

  2. Referee: Method section on pressure relaxation: the paper states that the mechanical-equilibrium model with relaxation registers tensile stresses as negative pressure, yet provides no explicit verification that the relaxation step preserves the correct characteristic speeds of the rarefaction across the VOF-tracked contact discontinuity. Without such a check, it is unclear whether the wave-propagation speed (and therefore the arrival time) is physically correct or numerically altered.

    Authors: The relaxation step is applied after the hyperbolic update and is formulated to conserve mass, momentum, and energy while enforcing a common pressure. To address the concern we have inserted a short verification subsection demonstrating a one-dimensional rarefaction propagating across a VOF-tracked interface. The computed wave speed and arrival time remain identical to the analytic solution both with and without relaxation, confirming that the characteristic speeds are unaltered by the procedure. revision: yes

  3. Referee: Equations of state: the independent EOS for lead, steel, and air are invoked but no parameter values, reference densities, or sound-speed formulations are supplied. These parameters directly determine the shock and rarefaction speeds; their omission prevents independent reproduction or assessment of the reported wave-arrival agreement.

    Authors: We thank the referee for noting this omission. A new table has been added to the Methods section that lists all EOS parameters and reference states: lead (Mie-Grüneisen: ρ₀ = 11340 kg m⁻³, c₀ = 2050 m s⁻¹, s = 1.46, Γ = 2.0), steel (Mie-Grüneisen: ρ₀ = 7850 kg m⁻³, c₀ = 4570 m s⁻¹, s = 1.49, Γ = 1.8), and air (ideal gas: γ = 1.4, p₀ = 10⁵ Pa). The sound-speed expressions used in the Riemann solver are now stated explicitly. revision: yes

Circularity Check

0 steps flagged

No circularity: validation rests on external experiment and independent codes

full rationale

The paper implements a standard multifluid Godunov-type scheme on the mechanical-equilibrium Euler equations with pressure relaxation and independent EOS for lead/steel/air, using VOF to track the interface. The headline result is the numerical arrival time of the unloading wave at the lead-steel contact, which is compared directly to separate experimental measurements and to simulations performed by other methods. No parameters are fitted to the target arrival-time datum, no self-citation supplies a uniqueness theorem or ansatz that forces the outcome, and the governing equations plus numerical algorithm are not defined in terms of the observed wave transit time. The reported agreement is therefore an external check rather than a tautological reproduction of an input.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on standard continuum assumptions for immiscible fluids and mechanical equilibrium; specific equations of state for lead, steel, and air are invoked but not enumerated here, implying they are taken from prior literature.

axioms (2)
  • domain assumption Mechanical equilibrium holds across phase interfaces at all times
    Invoked to justify the single-velocity multifluid Euler system with pressure relaxation
  • domain assumption Independent equations of state govern each immiscible phase
    Stated explicitly for metals and surrounding air

pith-pipeline@v0.9.0 · 5451 in / 1259 out tokens · 35595 ms · 2026-05-15T21:00:55.724786+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages

  1. [1]

    A. A. Deribas, Physics of Explosive Hardening and Welding (Nauka, Novosibirsk, 1980)

    work page 1980

  2. [2]

    Hydrodynamic effects in col- liding solids,

    S. K. Godunov, A. A. Deribas, A. V. Zabrodin, and N. S. Kozin, “Hydrodynamic effects in col- liding solids, ” J. Comput. Phys. 5(3), 517–539 (1970). https://doi.org/10.1016/0021-9991(70)90078-1

    work page doi:10.1016/0021-9991(70)90078-1 1970

  3. [3]

    Numerical simulation of a high-speed collision of metal plates,

    O. M. Belotserkovsky, S. V. Fortova, O. V. Troshkin, A. P. Pronina, I. V. Eriklintsev, and S. A. Kozlov, “Numerical simulation of a high-speed collision of metal plates, ” Math. Models Com- put. Simul. 8(5), 497–505 (2016). https://doi.org/10.1134/S2070048216050045

    work page doi:10.1134/s2070048216050045 2016

  4. [4]

    Instability of the interface between colliding metals,

    I. V. Yakovlev, “Instability of the interface between colliding metals, ” Combust. Explos. Shock Waves 9(3), 390–393 (1975). https://doi.org/10.1007/bf00745122. 20

    work page doi:10.1007/bf00745122 1975

  5. [5]

    AlMg6 to Titanium and AlMg6 to Stainless Steel Weld Interface Properties after Explosive Welding,

    A. Y. Malakhov, I. V. Saikov, I. V. Denisov, and N. N. Niyezbekov, “AlMg6 to Titanium and AlMg6 to Stainless Steel Weld Interface Properties after Explosive Welding, ” Metals 10(11), 1500 (2020). https://doi.org/10.3390/met10111500

    work page doi:10.3390/met10111500 2020

  6. [6]

    R. I. Nigmatulin, Dynamics of multiphase media. Part I (Nauka, Moscow, 1987; Hemisphere, New York, 1991)

    work page 1987

  7. [7]

    Fundamentals of the gas dynamics of interpenetrating motions of com- pressible media,

    K. A. Rakhmatulin, “Fundamentals of the gas dynamics of interpenetrating motions of com- pressible media, ” Prikl. Mat. Mekh. 20(2), 185–191 (1956)

    work page 1956

  8. [8]

    A two-phase mixture theory for the deflagration-to-detonation tran- sition (DDT) in reactive granular materials,

    M. Baer and J. Nunziato, “A two-phase mixture theory for the deflagration-to-detonation tran- sition (DDT) in reactive granular materials, ” Int. J. Multiph. Flow 12(6), 861–889 (1986). https://doi.org/10.1016/0301-9322(86)90033-9

    work page doi:10.1016/0301-9322(86)90033-9 1986

  9. [9]

    Numerical Simulation of a High-Speed Impact of Metal Plates Using a Three-Fluid Model,

    P. Chuprov, P. Utkin, and S. Fortova, “Numerical Simulation of a High-Speed Impact of Metal Plates Using a Three-Fluid Model, ” Metals 11(8), 1233 (2021). https://doi.org/10.3390/met11081233

    work page doi:10.3390/met11081233 2021

  10. [10]

    Two-phase modeling of deflagration-to-detonation transition in granular materials: Reduced equations,

    A. K. Kapila, R. Menikoff, J. B. Bdzil, S. F. Son, and D. S. Stewart, “Two-phase modeling of deflagration-to-detonation transition in granular materials: Reduced equations, ” Phys. Fluids 13(10), 3002–3024 (2001). https://doi.org/10.1063/1.1398042

    work page doi:10.1063/1.1398042 2001

  11. [11]

    Hypervelocity impact modeling with different equations of state,

    M. E. Povarnitsyn, K. V. Khishchenko, and P. R. Levashov, “Hypervelocity impact modeling with different equations of state, ” Int. J. Impact Eng. 33(1), 625–633 (2006). https://doi.org/10.1016/j.ijimpeng.2006.09.078

    work page doi:10.1016/j.ijimpeng.2006.09.078 2006

  12. [12]

    S. K. Godunov, A. V. Zabrodin, M. Y. Ivanov, A. N. Kraiko, and G. P. Prokopov, Numerical Solution of Multidimensional Problems of Gas Dynamics (Nauka, Moscow, 1976)

    work page 1976

  13. [13]

    J. S. Hesthaven and T. Warburton, Nodal Discontinuous Galerkin Methods (Springer, NY, 2008)

    work page 2008

  14. [14]

    Donea and A

    J. Donea and A. Huerta, Finite Element Methods for Flow Problems (Wiley, 2003)

    work page 2003

  15. [15]

    Computations of Compressible Multifluids,

    R. Abgrall and S. Karni, “Computations of Compressible Multifluids, ” J. Comput. Phys.169(2), 594–623 (2001). https://doi.org/10.1006/jcph.2000.6685

    work page doi:10.1006/jcph.2000.6685 2001

  16. [16]

    Advances in Interface Modeling,

    M. K. M. Ho, “Advances in Interface Modeling, ” in Handbook of Multiphase Flow Science and Technology (Springer Nature Singapore, 2023), pp. 137–164

    work page 2023

  17. [17]

    Volume of fluid (VOF) method for the dynamics of free bound- aries,

    C. W. Hirt and B. D. Nichols, “Volume of fluid (VOF) method for the dynamics of free bound- aries, ” J. Comput. Phys.39(1), 201–225 (1981). https://doi.org/10.1016/0021-9991(81)90145-5

    work page doi:10.1016/0021-9991(81)90145-5 1981

  18. [18]

    A volume-of-fluid reconstruction based inter- face sharpening algorithm for a reduced equation model of two-material compressible flow,

    M. Cutforth, P. T. Barton, and N. Nikiforakis, “A volume-of-fluid reconstruction based inter- face sharpening algorithm for a reduced equation model of two-material compressible flow, ” Comput. Fluids 231, 105158 (2021). https://doi.org/10.1016/j.compfluid.2021.105158

    work page doi:10.1016/j.compfluid.2021.105158 2021

  19. [19]

    A cell-centred Eulerian volume-of-fluid method for compressible multi-material flows,

    T. R. Law and P. T. Barton, “A cell-centred Eulerian volume-of-fluid method for compressible multi-material flows, ” J. Comput. Phys.497, 112592 (2024). https://doi.org/10.1016/j.jcp.2023.112592. 21

    work page doi:10.1016/j.jcp.2023.112592 2024

  20. [20]

    Free Boundary Method for Coupled Problems of Gas–Solid Dynamics,

    I. S. Menshov, “Free Boundary Method for Coupled Problems of Gas–Solid Dynamics, ” Com- put. Math. and Math. Phys. 64(8), 1809–1822 (2024). https://doi.org/10.1134/S0965542524700878

    work page doi:10.1134/s0965542524700878 2024

  21. [21]

    2-D Eulerian Hydrodynamics with Fluid Interfaces, Self-Gravity and Rotation,

    M. L. Norman and K.-H. A. Winkler, “2-D Eulerian Hydrodynamics with Fluid Interfaces, Self-Gravity and Rotation, ” inAstrophysical Radiation Hydrodynamics (Springer Netherlands, 1986), pp. 187–221

    work page 1986

  22. [22]

    A High-Order Godunov Method for Multiple Condensed Phases,

    G. H. Miller and E. G. Puckett, “A High-Order Godunov Method for Multiple Condensed Phases, ” J. Comput. Phys. 128(1), 134–164 (1996). https://doi.org/10.1006/jcph.1996.0200

    work page doi:10.1006/jcph.1996.0200 1996

  23. [23]

    A 3D adaptive mesh refinement algorithm for multimaterial gas dynamics,

    E. G. Puckett and J. S. Saltzman, “A 3D adaptive mesh refinement algorithm for multimaterial gas dynamics, ” Physica D60(1), 84–93 (1992). https://doi.org/10.1016/0167-2789(92)90228-F

    work page doi:10.1016/0167-2789(92)90228-f 1992

  24. [24]

    Simulation of multi-material hydrodynamic flows using adaptive mesh refinement,

    M. E. Povarnitsyn, A. S. Zakharenkov, P. R. Levashov, and K. V. Khishchenko, “Simulation of multi-material hydrodynamic flows using adaptive mesh refinement, ” Num. Meth. Prog. 13(3), 424–433 (2012)

    work page 2012

  25. [25]

    A. B. Wood, A Textbook of Sound (MacMillan, New York, 1930)

    work page 1930

  26. [26]

    Mathematical modeling of impact of two metal plates using two-fluid approach,

    P. S. Utkin and S. V. Fortova, “Mathematical modeling of impact of two metal plates using two-fluid approach, ” J. Phys.: Conf. Ser.946, 012047 (2018). https://doi.org/10.1088/1742-6596/946/1/012047

    work page doi:10.1088/1742-6596/946/1/012047 2018

  27. [27]

    A numerical technique for low-speed homogeneous two- phase flow with sharp interfaces,

    J. D. Ramshaw and J. A. Trapp, “A numerical technique for low-speed homogeneous two- phase flow with sharp interfaces, ” J. Comput. Phys. 21(4), 438–453 (1976). https://doi.org/10.1016/0021-9991(76)90039-5

    work page doi:10.1016/0021-9991(76)90039-5 1976

  28. [28]

    Restoration of the contact surface in the HLL-Riemann solver,

    E. F. Toro, M. Spruce, and W. Speares, “Restoration of the contact surface in the HLL-Riemann solver, ” Shock Waves4, 25–34 (1994). https://doi.org/10.1007/BF01414629

    work page doi:10.1007/bf01414629 1994

  29. [29]

    Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method,

    B. van Leer, “Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method, ” J. Comput. Phys. 32(1), 101–136 (1979). https://doi.org/10.1016/0021-9991(79)90145-1

    work page doi:10.1016/0021-9991(79)90145-1 1979

  30. [30]

    A direct Eulerian MUSCL scheme for gas dynamics,

    P. Colella, “A direct Eulerian MUSCL scheme for gas dynamics, ” SIAM J. Sci. Comput. 6(1), 104–117 (1985). https://doi.org/10.1137/0906009

    work page doi:10.1137/0906009 1985

  31. [31]

    Numerical Model of Multiphase Flows Based on Sub-cell Resolution of Fluid Interfaces,

    I. S. Menshov and A. A. Serezhkin, “Numerical Model of Multiphase Flows Based on Sub-cell Resolution of Fluid Interfaces, ” Comput. Math. and Math. Phys. 62(10), 1723–1742 (2024). https://doi.org/10.1134/s096554252209010x. 22

    work page doi:10.1134/s096554252209010x 2024