Pith. sign in

REVIEW 4 major objections 6 minor 31 references

Predicting The Evolution of Interfaces with Fourier Neural Operators

T0 review · 4 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash

Pith's one-line read A Fourier Neural Operator trained on volume-of-fluid simulations predicts liquid-vapour interface evolution with $R^2=0.95$ and 1–8 ms inference, making real-time multiphase control plausible.

desk verdict Useful first FNO-for-interface-evolution demo, but the headline R2=0.95 is on a smoothed field, not the interface, so the accuracy claim outruns the evidence. read the letter →

arxiv 2505.13463 v1 pith:63MKC4TL submitted 2025-05-05 cs.LG physics.comp-phphysics.flu-dyn

classification cs.LGphysics.comp-phphysics.flu-dyn
keywords FourierNeuralOperatorsoperatorlearningmultiphaseflowvolumeoffluidinterfacetrackingsigneddistancefunctionsurrogatemodelingreal-timecontrol
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper sets out to show that a Fourier Neural Operator (FNO) can act as a fast, accurate surrogate for volume-of-fluid simulations of liquid-vapour flows, predicting how the phase interface evolves in time. The authors train the operator on simulation data for a dam-break flow and for freely falling liquid blobs, representing the interface by a smoothed signed-distance field. They report validation $R^2 = 0.95$ and inference times of 1–8 ms, a speed that would make the surrogate fast enough to sit inside a real-time control loop, where conventional CFD is too slow. If the claim holds, neural operators become a practical forecasting layer for multiphase industrial processes and digital twins.

What carries the argument

The load-bearing mechanism is the Fourier Neural Operator's spectral convolution layer, which learns a kernel directly in Fourier space, truncates to a fixed number of modes, and thereby captures global spatial dependencies at a cost of $O(N\log N + k_{\max}N)$ per layer. To make the discontinuous volume-fraction field learnable, the paper maps $\alpha$ to a reconstructed distance function through $\zeta = \varepsilon\tanh^{-1}(1-2\alpha)$, a regularized inverse Heaviside transform that converts the sharp interface into a smooth field. The network takes the initial RDF plus a scalar time channel as input and outputs the predicted future RDF. The smoothing transform is what lets a continuous operator-learning architecture represent a discontinuous two-phase interface, while the spectral kernel is what lets the model propagate interface information across the whole domain in each layer.

What would settle it

Reconstruct the volume fraction from a predicted SDF via the inverse of Eq. (31), threshold it at $\alpha = 0.5$, and compare the resulting interface contours to the ground-truth VOF interface at $t = 0.25$ s and $t = 0.50$ s; if contour-position errors are large while SDF $R^2$ is 0.95, the central claim of accurate interface prediction fails.

Watch

Extended reading notes

Core claim

The central claim is that an FNO trained on compressed volume-of-fluid data learns the operator that maps an initial interface configuration to its future evolution, and that the learned operator generalizes both in time and across initial conditions. In the dam-break case, training on the first 60% of the simulated evolution lets the model forecast the remaining 40%, including the collapse and impact dynamics. In the second case, the model trained on two snapshots from 900 random initial conditions predicts the interface at later times for unseen configurations. Accuracy is measured on reconstructed distance fields (RDFs), where the model reaches MSE 9.72, MAE 1.69, and $R^2 = 0.95$, with 1–8 ms inference per prediction. The paper concludes that this makes neural-operator-based prediction fast enough to control multiphase applications.

Load-bearing premise

The accuracy score is computed on the smoothed reconstructed distance field and never on the physical volume fraction or on the interface location, and the smoothing parameter $\varepsilon$ is unspecified.

Editorial extensions

If this is right

  • Predictions at 1–8 ms make the surrogate a candidate for real-time monitoring and feedback control in multiphase systems, a regime where full CFD is too slow.
  • A model trained only on the early dam-break evolution can forecast subsequent impact and nonlinear deformation, suggesting the learned operator captures physics rather than memorizing the training window.
  • Two-snapshot training on random initial conditions generalizes to unseen configurations, which would let users build surrogates from sparse simulation or experimental data.
  • Because the method learns from volume-fraction fields produced by a standard VOF solver, the same pipeline can be re-trained for other two-phase geometries and fluids without changing the architecture.
  • The reported time scale matches the response time of many multiphase processes, so process controllers could embed the FNO as a predictive model.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper validates on smooth RDFs; a natural next check is to threshold the predicted RDF at zero and compare the extracted interface contour against the ground-truth VOF interface, since control actions depend on where the interface is, not just on field correlation.
  • The transform $\alpha \leftrightarrow \zeta$ is invertible for $\alpha \in (0,1)$, so every SDF prediction can be mapped back to a volume-fraction field; the paper does not report errors in that physical space, and doing so on the published dataset would settle how much of the $R^2$ score is inherited from smooth regions away from the interface.
  • The unspecified smoothing parameter $\varepsilon$ controls the sharpness of the reconstructed distance field; reporting its value and testing sensitivity to it would indicate whether the reported accuracy is robust or tuned to a favourable smoothness.
  • If the method extends to 3D and to phase change, as the paper names as next steps, the same operator-learning pipeline would apply to spray atomization and evaporating films, where interface evolution currently limits control bandwidth.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The manuscript proposes using Fourier Neural Operators (FNOs) as a surrogate for two-phase volume-of-fluid simulations. The authors train an FNO on signed-distance fields derived from the volume fraction α via a regularized inverse Heaviside transformation (Eq. 31), taking as input the initial SDF plus a scalar time channel and outputting the predicted SDF at a future time. Two benchmark cases are presented: a dam-break forecasting scenario and a dataset of random liquid-blob initial conditions. The authors report validation R²=0.95 on the SDF, inference times of 1–8 ms, and compare FNO with U-Net and GCN, concluding that FNO-based surrogates are fast enough for real-time control of multiphase processes.

Significance. If the accuracy claim held at the level of the physical interface, the work would be a useful demonstration of neural-operator surrogates for multiphase-flow control. The manuscript has clear strengths: it uses a high-fidelity VOF solver, provides a detailed architectural description, includes a baseline comparison with U-Net and GCN, and reports inference times. However, the reported accuracy is measured on a smooth transformed field rather than on the volume fraction or the extracted interface geometry, and the transformation contains an unreported free parameter. As it stands, the central claim of highly accurate interface prediction is not established.

major comments (4)
  1. [§3, Eq. (31) and Table 2] The central accuracy claim rests on metrics computed on the reconstructed distance function ζ, not on the physical volume fraction α or on the extracted interface. Because Eq. (31), ζ = ε·tanh⁻¹(1−2α), is singular at α=0 and α=1, the implementation must involve clipping or regularization, and the smoothing parameter ε is never reported. The reported MSE=9.72, MAE=1.69, and R²=0.95 are therefore not reproducible and are not interpretable as physical errors. Please report ε and the regularization strategy, and re-evaluate the predictions using α-based metrics (for example, errors restricted to interface cells with 0<α<1) or geometric interface measures such as the distance between predicted and true zero level sets (Hausdorff or contour distance).
  2. [§3, Case 1] The temporal-extrapolation claim for Case 1 is supported only by the qualitative comparison in Figure 1 and the loss curves in Figure 2; no quantitative error metric is reported for the held-out final 40% of the simulation. Since the case is explicitly presented as a test of extrapolation beyond the training horizon, the same numeric metrics used in Table 2 (MSE, MAE, R², and ideally interface-location error) should be reported for Case 1 separately.
  3. [§2 vs §3, Table 2] The validation metrics defined in Eqs. (27)–(30) are written in terms of α, yet Table 2 reports them on SDF fields. This inconsistency matters because a global R² on a smooth field is dominated by bulk regions, so R²=0.95 can coexist with systematically misplaced interfaces. Please state explicitly which field each metric is evaluated on, and add interface-focused error statistics (e.g., errors in the band |ζ|<ε or in cells with 0<α<1) so that the reader can assess the actual interface accuracy.
  4. [§3, Figure 3] The bottom row of Figure 3 shows 'binary volume fractions reconstructed from each RDF', but the reconstruction procedure is not described (for example, is it a threshold at ζ=0? Is there a re-initialization step?). Without this information, the visual comparison cannot be reproduced, and no quantitative error on α is given. Please describe the reconstruction and report error statistics for α as well.
minor comments (6)
  1. [§2, Figures 2 and 4] The 'Lp loss' label in Figures 2 and 4 is not defined; please specify the norm (presumably p=2) and use consistent terminology with Eq. (25).
  2. [§3, Case 1] The sentence 'The dataset comprises T simulations performed using compressibleInterIsoFoam' leaves T undefined; please provide the actual number of training and validation simulations used in Case 1.
  3. [§2, Hardware paragraph] The statement 'Training on a N=900 base simulation dataset (yielding 1800 SDF samples)' appears in the Hardware section but actually describes the Case 2 dataset; this placement is confusing and should be moved to the relevant experimental section.
  4. [Figure 2, caption (b)] The caption says 'Predicted vs. ground-truth volume fraction' while the main text refers to the signed-distance field; please clarify which field is plotted.
  5. [§4, Conclusions] The claim that the model 'successfully captures the underlying physics and extrapolates the evolution beyond the training window' is not substantiated by quantitative results for Case 1; please either soften this statement or add the missing metrics.
  6. [§3, Inference time] The 1–8 ms inference time is reported without specifying the hardware, the batch size, or the number of samples averaged, and there is no comparison with the per-step CFD solve time; please provide these details to support the real-time control claim.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the FNO is trained on held-out VOF simulation data and the RDF transform is a bijective preprocessing step, not a forced equivalence.

full rationale

The paper's derivation chain is a standard supervised operator-learning pipeline: VOF simulations generate volume-fraction fields, these are transformed into smooth reconstructed distance functions via Eq. (31), and an FNO is trained to map an initial RDF field to future RDF fields. Validation is performed on held-out simulation samples (40% temporal extrapolation in Case 1 and a 90/10 split in Case 2), so the reported R2 = 0.95 is an independent holdout metric, not a quantity fitted during training. The RDF transform is invertible in the interior of the domain (alpha = (1 - tanh(zeta/epsilon))/2), so predicting zeta is mathematically equivalent to predicting alpha up to the smoothing parameter; this is preprocessing, not a circular definition of the target. The only self-citation, reference [29], appears in the description of geometric curvature reconstruction in the OpenFOAM solver and is not load-bearing for the neural operator claim. A genuine validation concern exists: the global R2 is computed on the smooth SDF field, where bulk-cell contributions can dominate and the zero level set is not separately scored, so the paper has not directly quantified interface-position error. However, that is a correctness/reproducibility weakness, not circularity, because the predictions are not constructed from the ground-truth outputs. Overall, no step in the claimed derivation reduces to its own inputs by construction.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

Central claim rests on the fidelity of simulation data, the SDF preprocessing, and the generalization of a high-capacity fitted model; no new physical entities or forces are introduced.

free parameters (5)
  • Interface smoothing parameter epsilon
    In Eq. (31), epsilon controls the sharpness of the RDF transform; its value is not reported, and it directly affects the smoothness of the target field and hence the reported SDF errors.
  • Number of Fourier modes kmax = 20 per dimension
    Architecture hyperparameter chosen by hand; part of the model configuration, affecting truncation and accuracy.
  • Hidden channels dv = 96
    Architecture width chosen by hand.
  • Number of Fourier layers = 5
    Depth chosen by hand.
  • Learning rate, weight decay, batch size, epochs = 5e-4, 1e-4, 32, 100
    Optimization choices; not fitted to the data but determine the trained model.
assumptions (4)
  • domain assumption The compressibleInterIsoFoam simulation data provides accurate ground truth for the interface evolution.
    The surrogate's accuracy is measured entirely against these simulations; no experimental or independent validation is performed.
  • domain assumption The reconstructed distance function transform (Eq. 31) is smooth and carries the same interface information as alpha.
    The model is trained and evaluated on this transformed field; if the transform or its inverse loses fidelity at the interface, the reported metrics misrepresent interface accuracy.
  • domain assumption The learned operator G_theta approximates the true operator G in the L2 norm (Eq. 25) with sufficient accuracy.
    Standard neural operator approximation assumption; no theoretical guarantee is invoked for this specific flow.
  • domain assumption The training and test initial conditions are independent draws from the same distribution.
    Case 2 uses a random 90%/10% split; generalization beyond this distribution is not demonstrated.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Predicting The Evolution of Interfaces with Fourier Neural Operators." pith.science (2026). https://pith.science/paper/63MKC4TL

@misc{pith2026250513463,
  author       = {Pith},
  title        = {Pith review of: Predicting The Evolution of Interfaces with Fourier Neural Operators},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/63MKC4TL}},
  note         = {Machine review of arXiv:2505.13463}
}
read the original abstract

Recent progress in AI has established neural operators as powerful tools that can predict the evolution of partial differential equations, such as the Navier-Stokes equations. Some complex problems rely on sophisticated algorithms to deal with strong discontinuities in the computational domain. For example, liquid-vapour multiphase flows are a challenging problem in many configurations, particularly those involving large density gradients or phase change. The complexity mentioned above has not allowed for fine control of fast industrial processes or applications because computational fluid dynamics (CFD) models do not have a quick enough forecasting ability. This work demonstrates that the time scale of neural operators-based predictions is comparable to the time scale of multi-phase applications, thus proving they can be used to control processes that require fast response. Neural Operators can be trained using experimental data, simulations or a combination. In the following, neural operators were trained in volume of fluid simulations, and the resulting predictions showed very high accuracy, particularly in predicting the evolution of the liquid-vapour interface, one of the most critical tasks in a multi-phase process controller.

Figures

Figures reproduced from arXiv: 2505.13463 by the authors.

Figure 1
Figure 1. Time evolution of the liquid interface in the dam-break scenario. The figure shows three predicted time-steps [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Training diagnostics and spatial accuracy of the Fourier Neural Operator on the dam-break bench￾mark. (a) The evolution of the L p loss over 100 training epochs indicates consistent convergence of the surrogate model. (b) Compares the predicted and ground-truth spatial means of the signed-distance field (SDF) across training epochs, showing strong agreement and low variance. (c) Visualizes the spatial distribution o… view at source ↗
Figure 3
Figure 3. Comparison of predicted and ground-truth results at intermediate and final time steps. Each subfigure presents a visual comparison between the initial input, the ground-truth at a given time, and the predicted field from the FNO model. The top row of each panel shows the signed-distance fields (RDF): (left) the initial RDF at t = 0, (centre) the ground truth at the evaluated time, and (right) the FNO prediction. The… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Comparison of L p training and validation loss across architectures. The FNO exhibits superior con￾vergence behavior and consistently achieves the lowest validation loss, underscoring its ability to learn continuous operator mappings in complex flow regimes. While the …

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 27 canonical work pages

  1. [29]

    A computational study of thermally induced secondary atomization in multicomponent droplets

    Paolo Guida, Alberto Ceschin, Saumitra Saxena, Hong G Im, and William L Roberts. A computational study of thermally induced secondary atomization in multicomponent droplets. Journal of Fluid Mechanics, 935, 2022

  2. [1]

    Thermo-fluid dynamics of two-phase flow

    Mamoru Ishii and Takashi Hibiki. Thermo-fluid dynamics of two-phase flow. Springer, 2010

  3. [2]

    Direct numerical simulations of gas–liquid multiphase flows

    Gretar Tryggvason, Ruben Scardovelli, and Stephane Zaleski. Direct numerical simulations of gas–liquid multiphase flows. Cambridge University Press, 2011

  4. [3]

    Computational methods for multiphase flow

    Andrea Prosperetti and Gretar Tryggvason. Computational methods for multiphase flow. Cambridge University Press, 2009

  5. [4]

    A continuum method for modeling surface tension

    JU Brackbill, DB Kothe, and C Zemach. A continuum method for modeling surface tension. J. Comput. Phys., 100(2):335–354, 1992

  6. [5]

    V olume of fluid (vof) method for the dynamics of free boundaries.J

    CW Hirt and BD Nichols. V olume of fluid (vof) method for the dynamics of free boundaries.J. Comput. Phys., 39(1):201–225, 1981

  7. [6]

    Direct numerical simulation of gas–liquid multiphase flows using the volume-of-fluid method

    Ruben Scardovelli and Stephane Zaleski. Direct numerical simulation of gas–liquid multiphase flows using the volume-of-fluid method. Annu. Rev. Fluid Mech., 31(1):567–603, 1999

  8. [7]

    A front-tracking method for the computations of interfacial flows with sharp interfaces

    Gretar Tryggvason, Bernard Bunner, Asghar Esmaeeli, and Damir Juric. A front-tracking method for the computations of interfacial flows with sharp interfaces. J. Comput. Phys., 169(2):708–759, 2001

Show all 31 references
  1. [8]

    A volume-tracking method for incompressible multifluid flows with large density variations

    Murray Rudman. A volume-tracking method for incompressible multifluid flows with large density variations. Int. J. Numer. Meth. Fluids, 28(2):357–378, 1998

  2. [9]

    Gerris: a tree-based adaptive solver for the incompressible euler equations

    Stéphane Popinet. Gerris: a tree-based adaptive solver for the incompressible euler equations. J. Comput. Phys., 190(2):572–600, 2003

  3. [10]

    A new volume-of-fluid method in openfoam

    Johan Roenby, Bjarke Eltard Larsen, Henrik Bredmose, and Hrvoje Jasak. A new volume-of-fluid method in openfoam. In VII International Conference on Computational Methods in Marine Engineering. Nantes: International Center for Numerical Methods in Engineering, 2017

  4. [11]

    Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear pdes

    Maziar Raissi, Paris Perdikaris, and George E Karniadakis. Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear pdes. J. Comput. Phys., 378:686–707, 2019

  5. [12]

    Learning nonlinear operators via deeponet based on the universal approximation theorem

    Lu Lu, Pengzhan Jin, and George E Karniadakis. Learning nonlinear operators via deeponet based on the universal approximation theorem. Nat. Mach. Intell., 3:218–229, 2021. 10 Predicting interface evolution with FNOs

  6. [13]

    Neural operators for accelerating scientific simulations and design

    Kamyar Azizzadenesheli, Nikola Kovachki, Zongyi Li, Miguel Liu-Schiaffini, Jean Kossaifi, and Anima Anand- kumar. Neural operators for accelerating scientific simulations and design. Nature Reviews Physics, 6(5):320–328, 2024

  7. [14]

    Fourier neural operator for parametric partial differential equations

    Zongyi Li, Nikola Kovachki, Kamyar Azizzadenesheli, Burigede Liu, Kaushik Bhattacharya, Andrew M Stuart, and Anima Anandkumar. Fourier neural operator for parametric partial differential equations. arXiv preprint arXiv:2010.08895, 2020

  8. [15]

    Neural operator: Graph kernel network for pdes

    Nikola B Kovachki, Zongyi Li, Burigede Liu, Kamyar Azizzadenesheli, Kaushik Bhattacharya, Andrew M Stuart, and Anima Anandkumar. Neural operator: Graph kernel network for pdes. arXiv preprint arXiv:2103.12806, 2021

  9. [16]

    Physics-based deep learning for fluid flow modeling

    Nils Thuerey and et al. Physics-based deep learning for fluid flow modeling. Comput. Graph. Forum, 40(2):479– 496, 2021

  10. [17]

    Fourcastnet: A global data-driven high-resolution weather model using adaptive fnos

    Jaideep Pathak and et al. Fourcastnet: A global data-driven high-resolution weather model using adaptive fnos. arXiv preprint arXiv:2202.11214, 2022

  11. [18]

    Message passing neural pde solvers

    Johannes Brandstetter, Daniel Worrall, and Max Welling. Message passing neural pde solvers. arXiv preprint arXiv:2202.03376, 2022

  12. [19]

    Turbulence modeling in the age of data

    Karthik Duraisamy, Gianluca Iaccarino, and Heng Xiao. Turbulence modeling in the age of data. Annu. Rev. Fluid Mech., 51:357–377, 2019

  13. [20]

    Numerical models of surface tension

    Stéphane Popinet. Numerical models of surface tension. Annual Review of Fluid Mechanics, 50:49–75, 2018

  14. [21]

    Enhancing computational fluid dynamics with machine learning

    Ricardo Vinuesa and Steven L Brunton. Enhancing computational fluid dynamics with machine learning. Nat. Comput. Sci., 2:358–366, 2022

  15. [22]

    Ai4s: Artificial intelligence for science

    Stefano Markidis. Ai4s: Artificial intelligence for science. Computing in Science & Engineering, 23(5):5–8, 2021

  16. [23]

    Towards physics-informed foundation models for scientific machine learning

    Vishal Gupta and et al. Towards physics-informed foundation models for scientific machine learning. arXiv preprint arXiv:2205.12109, 2022

  17. [24]

    Validation of volume- of-fluid openfoam® isoadvector solvers using single bubble benchmarks

    Lionel Gamet, Marco Scala, Johan Roenby, Henning Scheufler, and Jean-Lou Pierson. Validation of volume- of-fluid openfoam® isoadvector solvers using single bubble benchmarks. Computers & Fluids, 213:104722, 2020

  18. [25]

    Openfoam: A c++ library for complex physics simulations

    Hrvoje Jasak, Aleksandar Jemcov, Zeljko Tukovic, et al. Openfoam: A c++ library for complex physics simulations. In International workshop on coupled methods in numerical dynamics, volume 1000, pages 1–20. IUC Dubrovnik Croatia, 2007

  19. [26]

    Accurate and efficient surface reconstruction from volume fraction data on general meshes

    Henning Scheufler and Johan Roenby. Accurate and efficient surface reconstruction from volume fraction data on general meshes. Journal of Computational Physics, 383:1–23, 2019

  20. [27]

    A continuum method for modeling surface tension

    J.U Brackbill, D.B Kothe, and C Zemach. A continuum method for modeling surface tension. Journal of Computational Physics, 100(2):335 – 354, 1992

  21. [28]

    Cummins, Marianne M

    Sharen J. Cummins, Marianne M. Francois, and Douglas B. Kothe. Estimating curvature from volume fractions. Computers & Structures, 83(6):425 – 434, 2005. Frontier of Multi-Phase Flow Analysis and Fluid-Structure

  22. [30]

    Adam: A method for stochastic optimization

    Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980, 2014

  23. [31]

    Improved adam optimizer for deep neural networks

    Zijun Zhang. Improved adam optimizer for deep neural networks. In 2018 IEEE/ACM 26th international symposium on quality of service (IWQoS), pages 1–2. Ieee, 2018. 11

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.