{"id":"b72fc648-934b-4919-bfdc-eeaaf2fa9f48","arxiv_id":"1908.04529","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"THINC-scaling uses the tanh function to convert between level set and volume-of-fluid fields, uniting the two in a single framework that maintains volume conservation and geometric accuracy.","lead":"This paper introduces THINC-scaling, a numerical scheme that converts between volume-of-fluid and level set descriptions of moving interfaces using a hyperbolic tangent function. It aims to combine the volume conservation of VOF with the accurate geometry of level set methods.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Geometric faithfulness rests on an unvalidated polynomial fit to non-signed-distance level set values in interface cells; a quantitative geometry check is needed.","rationale":"The paper's central claim has two parts: VOF conservation and level-set geometric faithfulness. The VOF update (Eq. 19) is a finite-volume scheme, so volume conservation is enforced by construction. The non-trivial part is geometric faithfulness, and that comes entirely from the polynomial PSI P_i+phi_delta. The polynomial P_i is computed in Step 1 from the discrete level set field. If that field is not an accurate signed distance function, or if the polynomial cannot represent the local interface topology (e.g., two branches of a thin filament in one cell), the PSI will be geometrically wrong. The reader's weakest_assumption identified the polynomial reconstruction as critical. I agree, but I sharpen it: the level set values in interface cells are deliberately not reinitialized (Step 5), so they are algebraic distances from the previous PSI, not signed distances. Fitting a polynomial to a stencil that mixes signed-distance and algebraic-distance values can bias the fit, and this is a distinct mechanism not spelled out in the reader's rationale. Since the paper only reports VOF L1 errors and visual contours, the claim of 'geometrical faithfulness' is unverified. The proposed test is a direct, quantitative check of the geometry in the thin-tail regime. Because the concern is specific and addressable by additional validation, the conditional verdict stands; the paper should not be rejected outright, but acceptance should require this geometric evidence or a softened claim.","tokens_in":8541,"tokens_out":12421,"duration_ms":127477,"concrete_test":"For the Rider-Kothe vortex test on a 128x128 grid, at t=T/2 compute the exact interface by backward integration of the velocity field. In the thin-tail cells, measure the error of the PSI zero set: mean normal angle error, mean curvature error, and whether two distinct interface branches are present in cells where the filament is sub-cell. Also compare the level set values in interface cells against the exact signed distance. If the normal error exceeds 5 degrees in the tail, or the PSI returns a single curve where the exact interface has two branches, the claim of geometric faithfulness for under-resolved interfaces is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Step 4.2 updates the level set at cell centers using the inverse THINC scaling of the departure cell's reconstruction H (Eq. 23). In interface cells this makes the new level set value equal to the PSI value P_i(x)+phi_delta, an algebraic distance, not a signed distance. Step 5 reinitializes only cells away from the interface, so these non-signed-distance values persist in interface cells. In the next Step 1, the polynomial P_i is fitted to a stencil that mixes reinitialized signed-distance cells with non-reinitialized algebraic-distance cells. This mismatch can bias the polynomial fit, and Step 2's VOF constraint only shifts the polynomial, not its shape. For under-resolved features such as the thin filament in the vortex test at t=T/2, a single polynomial zero set cannot represent two close interface branches inside one cell; the PSI may place the interface on an averaged curve. The paper's central claim of preserving the level set's geometrical faithfulness is therefore not established in the regime where the method is said to excel, and the paper reports only VOF L1 errors, not normal, curvature, or level set errors.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a novel interface-capturing scheme, THINC-scaling, that unifies VOF (volume of fluid) and level set methods. The key idea is an analytic conversion between a level set function and a continuous Heaviside-like VOF function via the THINC function: the forward THINC scaling maps a level set polynomial to a VOF field, and the inverse THINC scaling maps the VOF field back to a level set field. The paper gives a solution procedure (Section 3) that reconstructs a cell-wise polynomial from the discrete level set, shifts it to satisfy the VOF conservation constraint, transports the VOF function by a finite volume scheme, updates the level set by a semi-Lagrangian step using the inverse THINC scaling, and reinitializes the level set away from the interface. Numerical tests are presented for Zalesak's solid body rotation and the Rider-Kothe single vortex on both Cartesian and triangular unstructured grids, with VOF L1 errors compared against several geometric VOF methods. The central claim is that the scheme simultaneously preserves the mass/volume conservation of VOF and the geometrical faithfulness of level set methods.","tokens_in":8801,"tokens_out":5358,"duration_ms":56480,"significance":"If the central claim is substantiated, the paper offers a practically useful unification: a conservative interface-capturing method with sub-cell geometric information that is simpler than PLIC-based VOF and avoids the conservation defect of pure level set methods. The analytic conversion formulas (6) and (9) are clean and exact, the VOF constraint (16) enforces conservation by construction, and the scheme extends naturally to unstructured grids. The benchmarks show VOF L1 errors competitive with or better than published geometric VOF methods. However, the geometric-fidelity claim is currently supported only by qualitative contour plots; no quantitative level set or geometric error is reported, and the treatment of the level set in interface cells raises a correctness question that the paper does not address. For these reasons the result is promising but not yet fully established.","major_comments":[{"comment":"The quantitative assessment is limited to the VOF L1 error defined in Eq. (26). The abstract and Section 5 claim that THINC-scaling preserves the geometrical faithfulness of the level set method, but no error measure for the level set field, interface position, normal, or curvature is reported. Figures 4-6 show qualitative PSI plots only. Since geometric fidelity is the distinguishing advantage over algebraic VOF methods, please add quantitative geometry diagnostics, such as error in the zero-contour position relative to the exact solution, normal and curvature errors, or a measure of signed-distance accuracy in interface cells, for the Zalesak and vortex tests.","section":"Section 4, Eq. (26)"},{"comment":"In an interface cell, the inverse THINC scaling (23) yields phi(x) = psi_id(x) = P_id(x) + phi_delta, which is an algebraic distance to the PSI, not a signed distance. Step 5 reinitializes only cells away from the interface, so these non-signed-distance values persist in interface cells. The next Step 1 then fits P_i to a stencil that mixes reinitialized signed-distance values with non-reinitialized algebraic-distance values; this mismatch can bias the polynomial fit, and Step 2 only shifts the polynomial rather than correcting its shape. The paper does not quantify this effect. Please verify that the level set remains a consistent signed-distance-like function in interface cells, for example by reporting |grad phi| in interface cells or by testing sensitivity to the reinitialization band.","section":"Section 3, Steps 4.2 and 5"},{"comment":"The sub-cell geometric fidelity claim relies on a single polynomial P_i per cell. For a thin filament whose two interface branches pass through the same cell, as in the vortex test at t=T/2 (Fig. 5), a single polynomial zero set can represent the two branches only in special cases; in general the PSI will place the interface on an averaged curve. The manuscript provides no geometry error data for the thin tail region, so this limitation is not assessed. Please add a quantitative study of the filament region or explicitly discuss the topological limitation of a single polynomial representation.","section":"Section 3, Step 1 and Remark 6"},{"comment":"No comparison is made with the coupled THINC/level set method of Ref. [17], although Remark 5 identifies the semi-Lagrangian inverse-THINC level set update as the essential difference from [17]. A direct comparison with [17] on the same benchmarks is needed to substantiate that the new update rule is at least as accurate as the HJ-WENO update it replaces.","section":"Section 4, Tables 1 and 2"}],"minor_comments":[{"comment":"There are typographical errors: 'samectime' in the abstract and 'unstructured and unstructured grids' in Section 5 should read 'structured and unstructured grids'.","section":"Abstract and Section 5"},{"comment":"The caption of Fig. 1 contains garbled text ('0) x (/g32/g73'), and the phrase 'VOF velue' in Section 2 should be 'VOF value'.","section":"Section 2 and Fig. 1"},{"comment":"The notation in Eq. (8), using dx, dy, dz and D, is unclear; a multi-index notation would make the constraint condition easier to read.","section":"Section 2, Eq. (8)"},{"comment":"Existence and uniqueness of the solution phi_delta to Eq. (16) are not discussed. Since the Newton iteration is referenced to [16], a brief remark on the monotonicity of the left-hand side in phi_delta would help.","section":"Section 3, Step 2"},{"comment":"The header of Table 1 is corrupted ('Methods 642 Order 128 2 Order 256 2'); the formatting should be cleaned up.","section":"Section 4, Table 1"},{"comment":"The polynomial order p is fixed at 2 in all numerical tests, and no study of the dependence on p is reported. Since the paper emphasizes high-order interface representation, a p-convergence test or a discussion of the practical range of p would strengthen the claim.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The central idea is sound and the VOF results are competitive, but the main advertised advantage, geometric faithfulness of the level set, is not quantitatively demonstrated. The issue raised in Steps 4.2 and 5 about algebraic-distance values persisting in interface cells is the most important technical point to resolve; it is fixable with additional diagnostics or a modified reinitialization strategy, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, the core idea is real: they use the THINC function as an exact analytic bridge between a VOF field and a level set field, and they exploit that bridge to update the level set semi-Lagrangianly without any WENO reconstruction. That is a concrete algorithmic simplification over their own coupled THINC/level set method (ref. 17), and it is the main reason to read the paper. Second, the paper does not actually demonstrate the thing its abstract promises most loudly—geometrical faithfulness. They report VOF L1 errors only, not normal, curvature, or distance-to-interface errors, and they never compare against the coupled THINC/level set method from the same group, which is the most relevant baseline. That is the soft spot, and it is real, but it is addressable rather than fatal.\n\nWhat is genuinely new: the inverse THINC scaling formula (Eq. 9) is elementary, but using it to convert the cell-wise THINC reconstruction into a pointwise level set value at the departure point is a neat trick. It avoids the cost and complexity of high-order interpolation for the level set advection. The per-cell polynomial surface, PSI, is also a nice way to preserve sub-cell interface geometry while keeping the VOF conservation constraint exact. On the standard Zalesak and vortex benchmarks, the L1 errors are competitive with geometric VOF methods on both Cartesian and triangular grids. The derivation of the conversion formulas is clean, and the VOF constraint in Eq. (16) is a conservation constraint, not a fitted parameter. I have no circularity concern.\n\nWhere I would push back: the term \"unify\" is doing a lot of work. The scheme is a hybrid that alternates between VOF and level set updates, not a single evolution equation, but that is semantics and not a flaw. The more substantive worry is the one your stress-test flags. In interface cells, Step 4.2 produces level set values that are algebraic distances from the PSI, not signed distances, and Step 5 only reinitializes away from the interface. So the next polynomial fit in Step 1 mixes signed-distance and non-signed-distance values, which can bias the shape of the polynomial, not just its shift. For thin filaments below grid resolution, a single quadratic PSI per cell cannot represent two close interface branches, and the level set will likely smooth them. The paper's own Fig. 5 shows the tail is captured, but they never measure how faithfully. That is a genuine gap, and it weakens the central claim. I would also note the paper would be stronger with code or a reproducibility statement, and with a comparison to ref. 17.\n\nAll that said, the algorithmic idea is sound and the errors on benchmark advection tests are good. The paper deserves a serious peer review—a referee should push for quantitative geometry metrics (normal error, curvature error, reinitialization bias) and a head-to-head against the coupled THINC/level set method. If those come out well, this becomes a citeable contribution. So: send it to review, but insist on the geometry check.","headline":"A clean, genuinely useful unification of VOF and level set via THINC scaling, but the paper undersells its own novelty and oversells geometric faithfulness without measuring it.","tokens_in":9263,"tokens_out":937,"would_cite":true,"duration_ms":11185,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M08","76M12","76T99"],"pacs":[],"model":"deepseek-v4-flash","headline":"The volume-of-fluid and level set interface methods are two views of one object: THINC-scaling converts each into the other, so one procedure keeps volume conservation and geometrical faithfulness at the same time.","keywords":["THINC scaling","volume of fluid","level set method","interface capturing","semi-Lagrangian update","volume conservation","high-order polynomial","unstructured grid"],"falsifier":"Run the Rider\\u2013Kothe single-vortex test on a 128 by 128 grid for the full forward\\u2013reverse cycle, measuring total VOF and the L1 error of the returned circle against the exact initial shape: if volume drifts beyond the reported L1 level or the tail region fails to return to the circle, the claimed simultaneous conservation and geometry are not achieved. A sharper check is to compare the stretched tail at $t=T/2$ on the same grid against a high-resolution reference; if the under-resolved filament is not reproduced, the polynomial reconstruction assumption fails.","tokens_in":8372,"feed_emoji":"🌊","tokens_out":9274,"duration_ms":82016,"temperature":0.7,"pith_summary":"The paper sets out to show that volume-of-fluid (VOF) and level set methods, usually built as completely separate solvers, are really two readings of the same field under a hyperbolic-tangent scaling. In the THINC-scaling scheme, a level set polynomial is passed through a tanh function whose cell average is forced to equal the given VOF value, and the level set value is recovered from that same tanh function by the inverse transformation. The VOF field is advanced by a conservative finite volume update, while the level set field is advanced by a semi-Lagrangian update from the reconstructed tanh function without extra spatial interpolation. On solid-body rotation and vortex-deformation benchmarks, on both Cartesian and triangular grids, the scheme reconstructs interfaces including filaments thinner than a cell, with L1 errors comparable to geometrical VOF methods. If the scheme works as claimed, multiphase flow codes would no longer need to couple two different interface representations and choose between conserving volume and resolving geometry.","feed_headline":"One scaling formula unifies VOF and level-set interface tracking","feed_subtitle":"A single THINC scaling converts level set to a conservative VOF field and back, catching sub-cell filaments.","key_machinery":"The machinery is the THINC scaling identity and its inverse: $H(x)=\\frac{1}{2}(1+\\tanh(\\beta P(x)))$ with $P$ a polynomial matching the level set derivatives, and $\\varphi(x)=\\frac{1}{\\beta}\\tanh^{-1}(2H-1)$. The scheme's working object is the cell-wise THINC function constructed under the exact VOF constraint, whose unknown shift $\\varphi_\\Delta$ is fixed by solving for the cell average, making the piecewise zero contour of $P+\\varphi_\\Delta$ the sub-cell interface. This one function does double duty: the finite volume update of the cell-averaged VOF uses it directly, and the semi-Lagrangian level set update evaluates it at the departure point through the inverse scaling, so no separate reconstruction operator is introduced. The steepness parameter $\\beta$ controls interface thickness, and the paper uses $\\beta\\approx 6$ for a three-cell-thick transition layer.","core_discovery":"The central claim is that the THINC function $H(x)=\\frac{1}{2}(1+\\tanh(\\beta P(x)))$, built from a level set polynomial $P$ constrained by the derivatives of the level set field, is exactly the object that connects the two methods. Given a cell-wise VOF value, the unknown shift $\\varphi_\\Delta$ in the polynomial argument is solved so that the cell average of $H$ matches the VOF value, which enforces volume conservation by construction. The inverse formula $\\varphi(x)=\\frac{1}{\\beta}\\tanh^{-1}(2H-1)$ then supplies the level set value at any point in the cell, including the departure point of a semi-Lagrangian trajectory, so the level set update needs no polynomial interpolation of the level set field itself. The zero contour of the shifted polynomial, called the Polynomial Surface of the Interface, gives the sub-cell interface geometry including position, normal, and curvature inside each interface cell. The paper argues that this two-way scaling unifies VOF and level set into a single conservative and geometrically faithful framework, and supports the argument with advection tests on structured and unstructured grids.","pith_inferences":["A testable extension the paper does not run is a three-dimensional case with topological change; the tanh scaling should in principle handle breakup and merger because the VOF cell average is conserved, but the reinitialization step would need to re-seed the level set after topology changes.","The conservation claim is exact only in the limit of large $\\beta$; at finite $\\beta$ the tanh function is a smeared Heaviside, so an editor would expect small volume drift on very long advection runs unless $\\beta$ is increased with resolution.","The same two-way scaling could be used as a post-processing tool: given a VOF field from any method, invert the tanh to obtain a smooth level set field and extract curvature from the polynomial surface, an application the paper only hints at.","A stress test worth running is vortex deformation with the reverse flow turned off and the filament stretched beyond one cell width; the authors show such a case at $T/2$, but quantitative conservation of the under-resolved tail over long time is not measured."],"forward_implications":["A single implementation can carry both the VOF field and the level set field, eliminating the need for two separate advection solvers and specialized coupling logic.","The sub-cell interface from the polynomial zero contour gives direct access to normals and curvature, so surface-tension and sharp-interface terms can be computed from the same reconstruction that conserves volume.","On unstructured triangular grids, the scheme reaches L1 errors lower than the compared geometrical VOF methods in the single-vortex test, suggesting that the algebraic tanh construction is not a compromise on unstructured meshes.","The semi-Lagrangian level set update removes the high-order WENO reconstruction used in earlier coupled THINC/level set schemes, which simplifies the algorithm while keeping sub-cell resolution.","The framework extends naturally to higher order interface polynomials, since the THINC scaling does not depend on a particular polynomial degree."],"supporting_citations":[{"why":"Provides the quadratic surface reconstruction and Gaussian quadrature formulas used to solve the VOF constraint for the polynomial shift.","marker":"[16]"},{"why":"The prior coupled THINC/level set scheme from which this paper departs by replacing WENO level set advection with the inverse-scaling semi-Lagrangian update.","marker":"[17]"},{"why":"Supplies the steepness parameter estimate that fixes beta from the desired interface thickness.","marker":"[14]"},{"why":"Defines the single-vortex deformation benchmark and the PLIC comparison tradition used for the main accuracy tests.","marker":"[4]"},{"why":"Defines the Zalesak solid-body rotation benchmark used for the first convergence checks.","marker":"[22]"},{"why":"The leading comparison method for the unstructured-grid vortex test and a source of the reported error table.","marker":"[25]"},{"why":"Supplies the fast sweeping method used to reinitialize the level set field in step 5.","marker":"[19]"}],"fun_headline_variants":["THINC scaling: one formula unifies VOF and level set","Single scaling formula merges VOF and level set methods","THINC-scaling: conservative VOF meets faithful level set","One polynomial weds VOF conservation to level set geometry","THINC mapping: convert level set to VOF and back"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the cell-wise polynomial fit to the discrete level set values remains accurate inside the cell even when an interface feature is thinner than the grid; if the fit is poor, the inverse THINC step hands an incorrect level set value to the next cycle, and the reinitialization may not fully repair it.","fun_headline_variants_meta":{"raw":{"variants":["THINC scaling: one formula unifies VOF and level set","Single scaling formula merges VOF and level set methods","THINC-scaling: conservative VOF meets faithful level set","One polynomial weds VOF conservation to level set geometry","THINC mapping: convert level set to VOF and back"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000259,"raw_usage":{"total_tokens":1550,"prompt_tokens":872,"completion_tokens":678,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":488,"completion_tokens_details":{"reasoning_tokens":593}},"tokens_in":488,"tokens_out":678,"duration_ms":6161,"temperature":1.0,"reasoning_tokens":593,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:39:08.101173+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the Rider\\u2013Kothe single-vortex test on a 128 by 128 grid for the full forward\\u2013reverse cycle, measuring total VOF and the L1 error of the returned circle against the exact initial shape: if volume drifts beyond the reported L1 level or the tail region fails to return to the circle, the claimed simultaneous conservation and geometry are not achieved. A sharper check is to compare the stretched tail at $t=T/2$ on the same grid against a high-resolution reference; if the under-resolved filament is not reproduced, the polynomial reconstruction assumption fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the quadratic surface reconstruction and Gaussian quadrature formulas used to solve the VOF constraint for the polynomial shift."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The prior coupled THINC/level set scheme from which this paper departs by replacing WENO level set advection with the inverse-scaling semi-Lagrangian update."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the steepness parameter estimate that fixes beta from the desired interface thickness."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the single-vortex deformation benchmark and the PLIC comparison tradition used for the main accuracy tests."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Zalesak solid-body rotation benchmark used for the first convergence checks."},{"cited_title":"Scheufler, J","cited_arxiv_id":null,"evidence_quote":"The leading comparison method for the unstructured-grid vortex test and a source of the reported error table."},{"cited_title":"Zhao, A fast sweeping method for eikonal equations, M athematics of computation 74 (250) (2005) 603–627","cited_arxiv_id":null,"evidence_quote":"Supplies the fast sweeping method used to reinitialize the level set field in step 5."}],"review_version":1}