Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

A p-Multigrid Accelerated Nodal Spectral Element Method for Free-Surface Incompressible Navier-Stokes Model of Nonlinear Water Waves

T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A σ-transformed spectral element discretization of free-surface Navier-Stokes achieves spectral convergence for steep waves, reproduces the submerged-bar experiment, and solves the pressure Poisson problem at O(n) cost with a geometric…

desk verdict Useful incremental SEM-INSE solver, but the O(n) scaling claim is only tested at small amplitude and the steep-wave 'spectral convergence' wording oversells the results. read the letter →

arxiv 2411.14977 v1 pith:Q7DLKRBY submitted 2024-11-22 math.NA cs.NAphysics.flu-dyn

Anders Melander , Wojciech Laskowski , Spencer J. Sherwin , Allan P. Engsig-Karup This is my paper · ORCID
classification math.NAcs.NAphysics.flu-dyn MSC 65M7065M6065N5576B1576D0565F10
keywords nonlinearwaterwavesfree-surfaceflowincompressibleNavier-Stokesequationsnodalspectralelementmethodsigma-transformgeometricp-multigridpressurePoissonproblemtime-domainsimulation
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

The paper builds a numerical wave tank from the two-dimensional incompressible Navier-Stokes equations with a tracked free surface, and claims that a nodal spectral element discretization on a σ-transformed, time-frozen mesh achieves spectral (exponential) convergence for highly nonlinear waves. It further claims that the model reproduces the classical submerged-bar experiment on harmonic wave generation, and that the pressure Poisson equation—normally the computational bottleneck—can be solved at each Runge-Kutta stage by a geometric p-multigrid preconditioner with iteration counts that stay bounded as resolution grows, giving O(n) total cost. The broader aim is a solver that combines finite-element geometric flexibility with spectral accuracy and near-linear scaling as a stepping stone to three-dimensional wave-structure simulations.

What carries the argument

The central mechanism is the σ-transform, $\sigma=(z+h(x))/(\eta(x,t)+h(x))$, mapping the time-dependent physical fluid domain to a fixed reference strip so that the mesh and operators do not need rebuilding at each time step. On the reference domain, a nodal spectral element method with tensor-product Legendre bases and Legendre-Gauss-Lobatto nodes provides the spatial discretization, with gradient recovery through global $L^2$ projections and a modal spectral filter to stabilize the nonlinear terms. The accuracy claim is carried by the mixed-stage pressure Poisson problem, whose operator blends metric coefficients from two consecutive Runge-Kutta stages; solving it at every stage enforces $\nabla_\sigma\cdot u=0$ and defines the pressure-velocity coupling. The scalability claim is carried by the geometric p-multigrid V-cycle—coarsening by halving polynomial order, additive Schwarz smoothing, and modal transfer operators with restriction as the transpose of prolongation—used as a preconditioner for PDC and GMRES, with its operators assembled once from a linearized small-amplitude wave state rather than recomputed per stage.

What would settle it

Track the pressure-Poisson iteration count at every Runge-Kutta stage while simulating waves near 90% of the maximal breaking steepness for many wave periods, and repeat for the submerged-bar case. If the count climbs as the wave train evolves or steepens, the frozen linearized preconditioner of Section 5.4 is not doing the claimed work, and the $O(n)$ scaling measured at tolerance $10^{-6}$ in Section 6.3 would not persist.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that a nodal spectral element discretization of the σ-transformed incompressible Navier-Stokes equations, coupled with a mixed-stage pressure-velocity projection, achieves spectral error convergence for highly nonlinear waves while preserving mass conservation. The mixed-stage Poisson problem for the dynamic pressure, $\nabla_\sigma^{(k)}\cdot\nabla_\sigma^{(k-1)} p_D^{(k-1)} = \cdots$, mixes metric coefficients from two consecutive Runge-Kutta stages and must be solved at every stage; the paper shows this bottleneck can be handled by a geometric p-multigrid V-cycle used as a preconditioner for PDC and GMRES, with the multigrid operators frozen from a linearized small-amplitude wave state at the start of the simulation. Numerical evidence is offered in three parts: velocity error after one time step falls exponentially with polynomial order $P=2,\ldots,8$ at depths $kh=0.5,2,2\pi$ and steepness up to 90% of the breaking limit; the free-surface signal matches the Beji-Battjes measurements at six gauge positions for harmonic generation over a submerged bar; and pressure-solve iteration counts stay flat between roughly $10^3$ and $10^5$ degrees of freedom, whether resolution is added through element count or polynomial order, which the authors read as $O(n)$ computational scalability.

Load-bearing premise

The $O(n)$ scalability claim depends on the assumption that a preconditioner assembled once from a linearized small-amplitude wave state stays effective for the exact nonlinear, moving-domain pressure Poisson problem at every Runge-Kutta stage; the paper does not test whether iteration counts grow with wave steepness or with simulation time.

Editorial extensions

If this is right

  • Spectral convergence means users can raise polynomial order instead of refining the mesh to reach engineering accuracy, and the p-multigrid preconditioner keeps the added cost near-linear.
  • The frozen preconditioner built from a linearized small-amplitude wave state removes the per-stage operator re-assembly that moving-domain formulations would otherwise require.
  • Close agreement with the submerged-bar data indicates the model captures shoaling and harmonic decomposition of nonlinear waves over uneven bathymetry.
  • The solver strategy is claimed to extend naturally to three space dimensions and to be highly suitable for massively parallel implementations, where the pressure Poisson problem dominates cost.
  • The p-multigrid approach applies equally when resolution grows through element count or through polynomial order, with iteration counts staying flat in both directions.

Reading between the lines

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

  • The paper's own comparison table shows a direct solver beating the multigrid preconditioners at tight tolerances on the small test problem; a crossover study on larger and three-dimensional problems would show where the O(n) advantage takes over, a comparison the paper leaves implicit.
  • The exponential-convergence evidence is a single-step study; whether the spectral rate survives long-time integration, where the modal filter and temporal errors accumulate, is not quantified in the paper.
  • The scaling data come from a thin domain with only two vertical elements; three-dimensional grids or boundary-layer-refined meshes change the sparsity pattern, so flat iteration counts would need re-demonstration there.
  • The natural stress test of the frozen preconditioner is a long steep-wave run near 90% of the breaking limit: per-stage iteration counts would show whether the linearization degrades over time.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper presents a nodal spectral element method for the sigma-transformed free-surface incompressible Navier-Stokes equations, using explicit fourth-order Runge-Kutta time integration with a mixed-stage pressure-velocity coupling and a Poisson problem solved by a geometric p-multigrid preconditioner combined with PDC or GMRES. The model is verified against Rienecker-Fenton streamfunction wave solutions for a range of depths and steepnesses, validated against the Beji-Battjes submerged-bar experiment, and assessed for solver scalability. The central claims are spectral (exponential) convergence for highly nonlinear waves, excellent agreement with the bar experiment, and O(n) computational scalability of the p-multigrid-accelerated Poisson solver.

Significance. If fully substantiated, this is a useful contribution: it extends the single-domain pseudospectral model of Melander et al. to a geometrically flexible multi-domain spectral element setting, with a careful derivation of the mixed-stage Poisson problem and its weak form. The use of external analytical and experimental benchmarks is a strength, as is the explicit attention to the pressure Poisson problem as the main computational bottleneck. The p-multigrid transfer and smoothing framework follows established constructions and is presented in sufficient detail for implementation, apart from the missing parameter values noted below. However, the two headline claims are not fully supported by the evidence as presented: the high-steepness convergence plots do not clearly demonstrate exponential convergence, and the O(n) scalability claim rests on a single small-amplitude test case. The paper would be publishable after these load-bearing points are addressed and the missing numerical parameters are reported.

major comments (3)
  1. [Section 6.1, Figure 3(g)-(i)] The abstract and conclusion claim spectral convergence in the velocity fields for highly nonlinear waves, but the 90% steepness cases in Figure 3(g)-(i) do not show an exponential error decay. The errors decrease from about 1e-2 to 1e-4 over P=2 to P=8, with the curves appearing to flatten at high P, which is more suggestive of algebraic or saturation behavior than exponential convergence. To support the claim, the authors should report the measured convergence rates (e.g., fitted slopes or exponential decay constants) for the 90% cases and should state the error tolerance, time-step size, and spectral filter parameters used in these runs, since any of these could limit the observed convergence.
  2. [Section 5.4 and Section 6.3] The O(n) scalability claim is load-bearing and is not yet established for the nonlinear moving-domain problem. The preconditioner M is built once from a linearized small-amplitude-wave formulation, while the actual mixed-stage Poisson operator changes at every Runge-Kutta stage through the sigma-transform. The only scaling experiments, Table 1 and Figure 6, use kh=1 and (H/L)=0.0301, i.e., 30% of the maximum steepness, which is precisely the small-amplitude regime for which M was designed. No experiment varies the wave steepness toward the 90% cases used in Section 6.1, nor integrates over many periods, to show that iteration counts and O(n) time scaling persist when the domain deformation is large. Without such a robustness study, the advertised O(n) scalability for nonlinear waves remains a conjecture.
  3. [Section 4.4 and Section 5.3] Several numerical parameters that directly control the reported results are not specified, which prevents reproduction and weakens the quantitative claims. Section 4.4 defines the spectral filter S(i) with parameters alpha, beta, and Pc but never gives their values for the convergence or bar tests. Section 5.3 defines the additive Schwarz smoother with a user-determined overlap but does not report the overlap width or the number of smoothing iterations per multigrid level. These values affect both stability and the iteration counts in Table 1 and Figure 6, so they should be reported for every numerical experiment.
minor comments (4)
  1. [Section 6.2] The wave period is stated as T=0.02 s, which is inconsistent with H=0.02 m and L=3.74 m in a wave tank; this is almost certainly a typo for a period near 2 s. Please correct it so the experimental setup is reproducible.
  2. [Figure 3 caption] Subfigures (c), (f), and (i) are labeled 'kh = 2' but the text in Section 6.1 states that the deep-water case is kh = 2π; the caption labels appear to be a typo and should be fixed.
  3. [Section 6.1] The statement that 'the time step size is chosen such that spatial errors dominate' would be more useful if the actual time step or CFL criterion were reported, especially because the convergence study is performed after one time step and temporal and spatial errors interact.
  4. [Throughout] Please correct minor typos and grammar issues, including 'Kroenecker' (should be 'Kronecker'), 'The slope of that bar' (should be 'The slope of the bar'), 'O(n) computational scalability simulations' in the abstract, and 'orthonomal' in Section 4.3.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central accuracy claims are verified against external analytical and experimental data, and the O(n) solver scaling is a measured property, not an input.

full rationale

The paper's central claims are (i) spectral convergence of the nodal SEM discretization for nonlinear free-surface waves, (ii) agreement with the Beji-Battjes submerged-bar experiment, and (iii) O(n) scalability of the geometric p-multigrid accelerated Poisson solver. None of these reduce to the model's own assumptions or to a fitted parameter. The convergence study compares against the external Rienecker-Fenton streamfunction solution, and the bar test compares against external experimental data. The O(n) scalability is reported from timing and iteration-count measurements (Figure 6, Table 1) rather than derived from an assumption that already contains the result. Self-citations such as [1], [19], and [34] are used as building blocks for the discretization and solver strategy, but they do not constitute the verification of the central numerical results; the p-multigrid operators, transfer operators, and smoother are described and tested in the present paper. The only notable concern is that the frozen linearized preconditioner of Section 5.4 is tested only at kh=1 and 30% of maximum steepness, so the robustness of O(n) scaling for steeper or longer simulations is not demonstrated. That is an efficiency-claim and robustness gap, not a circularity: the O(n) claim is not true by construction and would fail if iteration counts grew, but no equation in the paper defines the scaling into existence. The paper is therefore self-contained with respect to its primary validation claims, and no circular step can be exhibited with the required quotation and reduction.

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

The model introduces no new physical entities. Its numerical parameters (filter, multigrid settings, time step) are not fully reported, which limits reproducibility, but they are not fitted to external data for the convergence benchmarks. The linearized preconditioner is an ad hoc modeling choice tied to the O(n) claim.

free parameters (4)
  • Spectral filter parameters (alpha, beta, Pc) = not reported
    Equation (62) defines an exponential cutoff filter applied to eta, u, w. The filter strength is chosen to remove only a few percent of highest-mode energy, but exact values are not given; results may depend on them.
  • Multigrid smoother overlap and iteration counts = not reported
    Additive Schwarz smoother in (77)-(78) requires overlap and iteration counts; these affect the reported iteration counts and timings in Table 1 and Figure 6.
  • Time step size = not reported
    Convergence study says Delta t is chosen so spatial errors dominate; the value is not given.
  • Relaxation zone lengths and generation/absorption setup = not reported
    Appendix A defines piecewise relaxation functions with domain bounds x1, x2; exact values for the bar test are not fully specified.
assumptions (6)
  • domain assumption The sigma transform (7) is a smooth, invertible mapping from the time-dependent physical domain to a fixed reference domain for non-overturning waves.
    Used throughout Sections 2.2 and 3 to derive transformed equations; fails for breaking or overturning waves, which the model excludes.
  • domain assumption Slip (impermeability) boundary conditions at walls and bathymetry (2) and zero dynamic pressure at the free surface (33) define the numerical wave tank.
    Adopted from wave tank practice; the bar benchmark relies on these conditions.
  • ad hoc to paper The linearized small-amplitude wave formulation yields a time-constant preconditioner that remains effective for the nonlinear moving-domain Poisson problem.
    Section 5.4 uses this to avoid recomputing multigrid operators each stage; no proof or steepness sensitivity study is given.
  • ad hoc to paper Spectral filtering with an exponential cutoff does not alter the convergence order of the scheme.
    Section 4.4 applies the filter to stabilize the scheme; no analysis shows that the chosen filter preserves spectral convergence.
  • domain assumption Rienecker-Fenton streamfunction solutions are valid exact references for inviscid, irrotational waves even for steep cases up to 90% of the breaking limit.
    Used as the exact solution in Section 6.1 convergence study; assumes non-viscous, irrotational flow with nu=0.
  • standard math Integration by parts and the chain rule for the sigma transform are applied correctly in deriving weak forms (41)-(43) and the mixed-stage Poisson operator (32).
    The derivation is presented but not formally verified; errors here would affect the whole scheme.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A p-Multigrid Accelerated Nodal Spectral Element Method for Free-Surface Incompressible Navier-Stokes Model of Nonlinear Water Waves." pith.science (2026). https://pith.science/paper/Q7DLKRBY

@misc{pith2026241114977,
  author       = {Pith},
  title        = {Pith review of: A p-Multigrid Accelerated Nodal Spectral Element Method for Free-Surface Incompressible Navier-Stokes Model of Nonlinear Water Waves},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Q7DLKRBY}},
  note         = {Machine review of arXiv:2411.14977}
}
abstract

We present a spectral element model for general-purpose simulation of non-overturning nonlinear water waves using the incompressible Navier-Stokes equations (INSE) with a free surface. The numerical implementation of the spectral element method is inspired by the related work by Engsig-Karup et al. (2016) and is based on nodal Lagrange basis functions, mass matrix-based integration and gradient recovery using global $L^2$ projections. The resulting model leverages the high-order accurate -- possibly exponential -- error convergence and has support for geometric flexibility allowing for computationally efficient simulations of nonlinear wave propagation. An explicit fourth-order accurate Runge-Kutta scheme is employed for the temporal integration, and a mixed-stage numerical discretization is the basis for a pressure-velocity coupling that makes it possible to maintain high-order accuracy in both the temporal and spatial discretizations while preserving mass conservation. Furthermore, the numerical scheme is accelerated by solving the discrete Poisson problem using an iterative solver strategy based on a geometric $p$-multigrid method. This problem constitutes the main computational bottleneck in INSE models. It is shown through numerical experiments, that the model achieves spectral convergence in the velocity fields for highly nonlinear waves, and there is excellent agreement with experimental data for the simulation of the classical benchmark of harmonic wave generation over a submerged bar. The geometric $p$-multigrid solver demonstrates $O(n)$ computational scalability simulations, making it a suitable efficient solver strategy as a candidate for extensions to more complex, real-world scenarios.

Figures

Figures reproduced from arXiv: 2411.14977 by the authors.

Figure 1
Figure 1. Illustration of the physical and σ-transformed domains. The description of the evolution of water waves in the time domain can be described by the incompressible Navier-Stokes equations (INSE) along with a free surface kinematic boundary condition [17, 35]. By assuming that the fluid density is constant, the fluid flow becomes divergence free. In the following, we introduce the fluid domain Ω ⊂ R d (d = 2) be a boun… view at source ↗
Figure 2
Figure 2. Example of V -cycle for geometric p-multigrid on a 2 × 2 mesh. 5.1. Grid strategy An integral part of the geometric p-multigrid method is establishing the grid hierarchy. As mentioned, this is done by reducing the polynomial order p of the SEM discretization. We let GN be the original fine grid of order PN and G1 be the coarsest grid of order P1. The subsequent coarser grids are then generated by reducing the order … view at source ↗
Figure 3
Figure 3. Error convergence for the velocities u and w across varying depth and wave steepness. 6.2. Harmonic generation over a submerged bar To show that the solver is capable of handling both wave-tank setups, as well as non-flat bathymetry, we present results of harmonic generation over a submerged bar, which was performed both experimentally and numerically by Beji and Battjes [44]. The setup of the experiment can be seen… view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Setup of the bar test. 21 [PITH_FULL_IMAGE:figures/full_fig_p021_4.png]
Figure 5
Figure 5. Figure 5: Comparison of numerical and experimental results [PITH_FULL_IMAGE:figures/full_fig_p022_5.png]
Figure 6
Figure 6. Figure 6: shows the iteration count and computational scaling when either the horizontal order or element count is increased. Both methods shows excellent scalability, with the iteration count showing very little variation, with no tendency to increase at higher DoFs. Moreover, …

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. FNPF-SEM: A parallel spectral element model in Firedrake for fully nonlinear water wave simulations

    math.NA 2025-06 conditional novelty 5.0 of 10

    A Firedrake-based spectral element model for fully nonlinear potential flow water waves is verified through convergence and scaling tests and validated against experimental and analytical benchmarks.

Reference graph

Works this paper leans on

48 extracted references · 37 canonical work pages · cited by 1 Pith paper

  1. [19]

    A High-Order Hybrid-Spectral Incompressible Navier-Stokes Model For Nonlinear Water Waves

    A. Melander, M. E. Bitsch, D. Chen, A. P. Engsig-Karup, A high-order hybrid-spectral incom- pressible navier-stokes model for nonlinear water waves, 2 024. doi: 10.48550/arXiv.2406.00991. arXiv:2406.00991

    work page Pith review arXiv doi:10.48550/arxiv.2406.00991

  2. [1]

    Engsig-Karup, C

    A. Engsig-Karup, C. Eskilsson, D. Bigoni, A stabilised n odal spectral element method for fully nonlinear water waves, J. Comput. Phys. 318 (2016) 1–21. doi :10.1016/j.jcp.2016.04.060

    doi:10.1016/j.jcp.2016.04.060 2016

  3. [2]

    Ransley, S

    E. Ransley, S. Yan, S. A. Brown, T. Mai, D. Graham, Q. Ma, P. -H. Musiedlak, A. P. Engsig-Karup, C. Eskilsson, Q. L. J. W ang, Z. Xie, S. Venkatachalam, T. Stoe sser, Y. Zhuang, Q. Li, D. W an, G. Chen, H. Chen, L. Qian, Z. Ma, C. Mingham, D. Causon, I. Gati n, H. Jasak, V. Vukcevic, S. Downie, P. Higuera, E. Buldakov, D. Stagonas, Q. Chen, J. Z ang, D. G...

    work page 2019

  4. [3]

    Sriram, S

    V. Sriram, S. Agarwal, S. Yan, Z. Xie, S. Saincher, T. Schl urmann, Q. Ma, T. Stoesser, Y. Zhuang, B. Han, W. Zhao, X. Yang, Z. Li, D. W an, Y. Zhang, B. Teng, D. Nin g, N. Zhang, X. Zheng, G. Xu, J. Gong, Y. Li, K. Liao, W. Duan, R. Han, W. Asnim, Z. Sulaiman, Z. Zhou, J. Qin, Y. Li, Z. Song, X. Lou, L. Lu, C. Yuan, Y. Ma, C. Ai, G. Dong, H. Sun, Q. W an...

    work page 2021

  5. [4]

    Osher, J

    S. Osher, J. A. Sethian, Fronts propagating with curvatu re-dependent speed: Al- gorithms based on Hamilton-Jacobi formulations, J. Comp. P hys 79 (1988) 12–49. doi:10.1016/0021-9991(88)90002-2

    doi:10.1016/0021-9991(88)90002-2 1988

  6. [5]

    Grooss, J

    J. Grooss, J. Hesthaven, A level set discontinuous galer kin method for free surface flows, Comput. Methods Appl. Mech. Eng. 195 (2006) 3406–3429. doi: 10.1016/j.cma.2005.06.020

    doi:10.1016/j.cma.2005.06.020 2006

  7. [6]

    H. Bihs, A. Kamath, M. Alagan Chella, A. Aggarwal, Øivind A. Arntsen, A new level set numerical wave tank with improved density interpolation for complex w ave hydrodynamics, Comput. Fluids 140 (2016) 191–208. doi: 10.1016/j.compfluid.2016.09.012

    work page doi:10.1016/j.compfluid.2016.09.012 2016

  8. [7]

    C. Hirt, B. Nichols, Volume of fluid (VOF) method for the dy namics of free boundaries, J. Comp. Phys. 39 (1981) 201–225. doi: 10.1016/0021-9991(81)90145-5

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

Show all 48 references
  1. [8]

    W eller, G

    H. W eller, G. Tabor, H. Jasak, C. Fureby, A tensorial appr oach to computational con- tinuum mechanics using object orientated techniques, Comp ut. Phys. 12 (1998) 620–631. doi:10.1063/1.168744. 24

    1998 doi

  2. [9]

    P. Lin, C. Li, A σ -coordinate three-dimensional numerical model for surfac e wave propagation, Int. J. Num. Meth. Fluids. 38 (2002) 1045 – 1068. doi: 10.1002/fld.258

    2002 doi

  3. [10]

    Donea, A

    J. Donea, A. Huerta, J.-P. Ponthot, A. Rodr ´ ıguez-Ferran, Arbitrary Lagrangian–Eulerian Methods, Encyclopedia of computational mechanics (2004). doi: 10.1002/0470091355.ecm009

    2004 doi

  4. [11]

    M. S. Longuet-Higgins, E. D. Cokelet, The deformation o f steep surface waves on water - I. A numer- ical method of computation, Proc. R. Soc. Lond. A 350 (1976) 1 –26. doi: 10.1098/rspa.1976.0092

    1976

  5. [12]

    A. P. Engsig-Karup, C. Monteserin, C. Eskilsson, A mixe d eulerian–lagrangian spectral element method for nonlinear wave interaction with fixed structures , W ater W aves 1 (2019) 315–342. doi:10.1007/s42286-019-00018-5

    2019 doi

  6. [13]

    B. Li, C. A. Fleming, Three-dimensional model of navier -stokes equations for water waves, J. W aterw. Port Coast. Ocean Eng. 127 (2001) 16–25. doi: 10.1061/(ASCE)0733-950X(2001)127:1(16)

    2001 doi

  7. [14]

    Decoene, J

    A. Decoene, J. Gerbeau, Sigma transformation and ale fo rmulation for three-dimensional free surface flows, Int. J. Numer. Methods Fluids 59 (2009) 357 – 38 6. doi: 10.1002/fld.1816

    2009 doi

  8. [15]

    Chorin, Numerical solution of the navier–stokes equ ations, Math

    A. Chorin, Numerical solution of the navier–stokes equ ations, Math. Comput. 22 (1968). doi:10.2307/2004575

    1968 doi

  9. [16]

    F. H. Harlow, J. E. W elch, Numerical calculation of time -dependent viscous incompressible flow of fluid with free surface, Phys. Fluids 8 (1965) 2182–2189. doi :10.1063/1.1761178

    1965 doi

  10. [17]

    Engsig-Karup, A

    A. Engsig-Karup, A. Melander, A high-order finite differ ence incompressible navier-stokes model for water waves and wave-structure interaction, 2024. Proc . 39th Int. W orkshop W ater W aves Float. Bodies (IWWWFB)

    2024

  11. [18]

    W. Pan, S. C. Kramer, M. D. Piggott, A σ -coordinate non-hydrostatic discontinuous finite element coastal ocean model, Ocean Model. 157 (2021) 101732. doi: 10.1016/j.ocemod.2020.101732

    2021

  12. [20]

    P. A. Madsen, H. B. Bingham, H. Liu, A new boussinesq meth od for fully nonlinear waves from shallow to deep water, J. Fluid Mech. 462 (2002) 1–30. doi: 10.1017/S0022112002008467

    2002 doi

  13. [21]

    A. P. Engsig-Karup, J. S. Hesthaven, H. Bingham, P. Mads en, Nodal DG-FEM so- lution of high-order boussinesq-type equations, J. Eng. Ma th. 56 (2006) 351–370. doi:10.1007/s10665-006-9064-z

    2006 doi

  14. [22]

    T. B. Christiansen, H. B. Bingham, A. P. Engsig-Karup, G . Ducrozet, P. Ferrant, Efficient hybrid- spectral model for fully nonlinear numerical wave tank, 201 3. Proc. ASME 2013 32nd Int. Conf. Ocean, Offshore Arctic Eng., 2013, Nantes, France

    2013

  15. [23]

    Klahn, P

    M. Klahn, P. A. Madsen, D. R. Fuhrman, A new sigma-transf orm based fourier-legendre- galerkin model for nonlinear water waves, Int. J. Numer. Met hods Fluids 93 (2020) 220–248. doi:10.1002/fld.4881

    2020 doi

  16. [24]

    Bonnefoy, G

    F. Bonnefoy, G. Ducrozet, D. Le Touz´ e, P. Ferrant, Time domain simulation of nonlinear water waves using spectral methods, W orld Scient ific, 2010, pp. 129–164. doi:10.1142/9789812836502\_0004

    2010 doi

  17. [25]

    Melander, A

    A. Melander, A. P. Engsig-Karup, A p-multigrid hybrid- spectral model for nonlinear water waves, W ater W aves (2024). doi:10.1007/s42286-024-00098-y

    2024 doi

  18. [26]

    A. T. Patera, A spectral element method for fluid dynamic s: Laminar flow in a channel expansion, J. Comput. Phys. 54 (1984) 468–488. doi: 10.1016/0021-9991(84)90128-1

    1984 doi

  19. [27]

    Maday, A

    Y. Maday, A. T. Patera, Spectral element methods for the incompressible navier-stokes equations, State-of-the-art Surveys on Computational Mechanics (198 9) 71–143

  20. [28]

    Robertson, S

    I. Robertson, S. Sherwin, Free-surface flow simulation using hp/spectral elements, J. Comput. Phys. 155 (1999) 26–53. doi: 10.1006/jcph.1999.6328

    1999

  21. [29]

    Kreiss, J

    H.-O. Kreiss, J. Oliger, Comparison of accurate method s for the integration of hyperbolic equations, Tellus 24 (1972) 199–215. doi: 10.1111/j.2153-3490.1972.tb01547.x

    1972

  22. [30]

    Trottenberg, C

    U. Trottenberg, C. W. Oosterlee, A. Schuller, Multigri d Methods: Basics, Parallelism and Adap- tivity, Academic Press, Cambridge, Massachusetts, 2001

    2001

  23. [31]

    Brandt, Barriers to Achieving Textbook Multigrid Effi ciency (TME) in CFD, Technical Report NASA/CR-1998-207647, ICASE Interim Report No

    A. Brandt, Barriers to Achieving Textbook Multigrid Effi ciency (TME) in CFD, Technical Report NASA/CR-1998-207647, ICASE Interim Report No. 32, NASA Lan gley Research Center, Institute for Computer Applications in Science and Engineering (ICAS E), Hampton, V A, 1998. Prepared f...

    1998

  24. [32]

    Thomas, B

    J. Thomas, B. Diskin, A. Brandt, Textbook multigrid effic iency for fluid simulations, Annu. Rev. Fluid Mech 17 (2002) 317–40. doi: 10.1146/annurev.fluid.35.101101.161209

    2002

  25. [33]

    Engsig-Karup, H

    A. Engsig-Karup, H. Bingham, O. Lindberg, An efficient fle xible-order model for 3D nonlinear water waves, J. Comput. Phys. 228 (2009) 2100–2118. doi: 10.1016/j.jcp.2008.11.028. 25

    2009 doi

  26. [34]

    A. P. Engsig-Karup, W. Laskowski, Efficient p-multigrid spectral element model for water waves and marine offshore structures, Int. J. Numer. Methods Fluid s 93 (2021). doi: 10.1002/fld.5011

    2021 doi

  27. [35]

    A. P. Engsig-Karup, L. S. Glimberg, A. S. Nielsen, O. Lin dberg, Fast hydrodynamics on heteroge- nous many-core hardware, 2013

    2013

  28. [36]

    Duarte, R

    F. Duarte, R. Gormaz, S. Natesan, Arbitrary lagrangian –eulerian method for navier–stokes equations with moving boundaries, Comput. Methods Appl. Me ch. Eng. 193 (2004) 4819–4836. doi:10.1016/j.cma.2004.05.003

    2004 doi

  29. [37]

    Carpenter, C

    M. Carpenter, C. Kennedy, Fourth-order 2N-storage Run ge-Kutta schemes, NASA reports TM, 109112 (1994)

    1994

  30. [38]

    E. M. Rønquist, A. T. Patera, Spectral element multigri d. i. formulation and numerical results, J. Sci. Comput. 2 (1987) 389–406. doi: 10.1007/BF01061297

    1987 doi

  31. [39]

    Lottes, P

    J. Lottes, P. Fischer, Hybrid multigrid/schwarz algor ithms for the spectral element method, J. Sci. Comput. 24 (2005) 45–78. doi: 10.1007/s10915-004-4787-3

    2005 doi

  32. [40]

    A. P. Engsig-Karup, Analysis of efficient preconditione d defect correction methods for nonlinear water waves, Int. J. Numer. Methods Fluids 74 (2014) 749–773 . doi: 10.1002/fld.3873

    2014 doi

  33. [41]

    Y. Saad, M. H. Schultz, Gmres: a generalized minimal res idual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Statist. Comput. 7 (1986) 856–8 69. doi: 10.1137/0907058

    1986 doi

  34. [42]

    M. M. Rienecker, J. D. Fenton, A Fourier approximation m ethod for steady water waves, Journal of Fluid Mechanics 104 (1981) 119–137. doi: 10.1017/S0022112081002851

    1981 doi

  35. [43]

    J. A. Battjes, Surf similarity, 1974. 14th Coastal Eng. Conf., 1974, Copenhagen, Denmark

    1974

  36. [44]

    S. Beji, J. A. Battjes, Numerical simulation of nonline ar wave propagation over a bar, Coast. Eng. 23 (1994) 1–16. doi: 10.1016/0378-3839(94)90012-4

    1994 doi

  37. [45]

    Larsen, H

    J. Larsen, H. Dancy, Open boundaries in short wave simul ations — a new approach, Coast. Eng. 7 (1983) 285–297. doi: 10.1016/0378-3839(83)90022-4

    1983 doi

  38. [46]

    A. P. Engsig-Karup, M. G. Madsen, S. L. Glimberg, A massi vely parallel gpu-accelerated model for analysis of fully nonlinear free surface waves, Int. J. N umer. Methods Fluids 70 (2012) 20–36. doi:10.1002/fld.2675

    2012 doi

  39. [47]

    S. L. Glimberg, A. P. Engsig-Karup, L. N. Olson, A massiv ely scalable distributed multigrid framework for nonlinear marine hydrodynamics, Int. J. High Perform. Comput. Appl. 33 (2019) 855–868. doi: 10.1177/1094342019826662

    2019 doi

  40. [48]

    A. P. Engsig-Karup, Unstructured Nodal DG-FEM solutio n of high-order Boussinesq-type equa- tions, Ph.D. thesis, Department of Mechanical Engineering , Technical University of Denmark, 2006. Appendix A. Generation and absorption zones To generate waves through the use generati...

    2006

Pith tools

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