Simulating forced time-periodic flows in industrial applications presents significant computational challenges, partly due to the need to overcome costly transients before achieving time-periodicity. Reduced-order modelling emerges as a promising method to speed-up computations. We extend upon the work of Lotz et al. (2024) where a time-periodic space–time model is introduced. We present a time-periodic reduced-order model that directly finds the time-periodic solution without requiring extensive time integration. The reduced-order model gives a reduction in variables in both space and time. Our approach involves a POD-Galerkin reduced-order model based on a time-periodic full-order model that employs isogeometric analysis, residual-based variational multiscale turbulence modelling and weak boundary conditions. The projection-based reduced-order model inherits these features. We evaluate the reduced-order model with numerical experiments on moving hydrofoils. The motion is known a priori and we restrict ourselves to two spatial dimensions. In these experiments we vary the Strouhal and Reynolds numbers, and the motion profile respectively. Reduced-order model solutions agree well with those of the full-order model. The errors over the entire time period of thrust and lift forces are less than 0.2%. This includes complex scenarios such as the transition from drag to thrust production with increasing Strouhal number. Our time-periodic reduced-order model offers speed-ups ranging from O(10²) to O(10³) compared to the full-order model, depending upon the basis size. This makes it an appealing solution for prescribed time-periodic problems, with potential for additional speedup through nonlinear reduction techniques such as hyper-reduction.

@en

In this paper, we build on the work of Hughes and Sangalli (2007) dealing with the explicit computation of the Fine-Scale Greens’ function. The original approach chooses a set of functionals associated with a projector to compute the Fine-Scale Greens’ function. The construction of these functionals, however, does not generalise to arbitrary projections, higher dimensions, or Spectral Element methods. We propose to generalise the construction of the required functionals by using dual functions. These dual functions can be directly derived from the chosen projector and are explicitly computable. We show how to find the dual functions for both the L² and the H₀¹ projections. We then go on to demonstrate that the Fine-Scale Greens’ functions constructed with the dual basis functions consistently reproduce the unresolved scales removed by the projector. The methodology is tested using one-dimensional Poisson and advection–diffusion problems, as well as a two-dimensional Poisson problem. We present the computed components of the Fine-Scale Greens’ function, and the Fine-Scale Greens’ function itself. These results show that the method works for arbitrary projections, in arbitrary dimensions. Moreover, the methodology can be applied to any Finite/Spectral Element or Isogeometric framework.

@en

The computation of periodic flows is typically conducted over multiple periods. First, a number of periods is used to develop periodic characteristics, and afterwards statistics are collected from averages over multiple periods. As a consequence, it is uncertain whether the numerical results are exactly time-periodic, and additionally, the time domain might be needlessly long. In this article, we circumvent these concerns by using a time-periodic function space. Consequently, the boundary conditions and solutions are exactly periodic. We employ the isogeometric analysis framework to achieve higher-order smoothness in both space and time. The discretization is performed using residual-based variational multiscale modeling and weak boundary conditions are adopted to enhance the accuracy near the moving boundaries of the computational domain. We enforce the time-periodic boundary condition within the isogeometric discretization spaces, which converts the two-dimensional time-dependent problem into a three-dimensional boundary value problem. Furthermore, we determine the boundary velocities of moving hydrofoils directly from the computational mesh and use a conservation methodology for force extraction. Application of the computational setup to heaving and pitching hydrofoils displays very accurate and exactly periodic results for lift and drag.

@en

Turbulent flows past rough surfaces can create substantial energy losses in engineering equipment. During the last decades, developing accurate correlations to predict the thermal and hydrodynamic behavior of rough surfaces has proven to be a difficult challenge. In this work, we investigate the applicability of convolutional neural networks to perform a direct image-to-image translation between the height map of a rough surface and its detailed local skin friction factors and Nusselt numbers. Additionally, we propose the usage of separable convolutional modules to reduce the total number of trainable parameters, and PReLU activation functions to increase the expressivity of the neural networks created. Our final predictions are improved by a new filtering methodology, which is able to combine the results of multiple neural networks while discarding non-physical oscillations likely caused by over-fitting. The main study is based on a new DNS database formed by 80 flow cases at a friction Reynolds number of Re_τ=180 obtained by applying random shifts to the Fourier spectrum of the grit-blasted surface originally scanned by Busse et al. (2015). The results show that machine learning can accurately predict the skin friction values and Nusselt numbers for a rough surface. A detailed comparison with existing correlations in the literature revealed that the maximum errors generated by deep learning were only 8.1% for the global skin friction factors C_f¯ and 2.9% for the Nusselt numbers Nu¯, whereas the best classical correlations identified reached errors of 24.9% and 13.5% for C_f¯ and Nu¯ respectively. The deep learning results also proved stable with respect to rough surfaces with abrupt changes in their roughness elements, and only presented a minor sensitivity with respect to variations in the dataset size.

@en

Over the last decades, many diffuse-interface Navier-Stokes Cahn-Hilliard (NSCH) models with non-matching densities have appeared in the literature. These models claim to describe the same physical phenomena, yet they are distinct from one another. The overarching objective of this work is to bring all of these models together by laying down a unified framework of NSCH models with non-zero mass fluxes. Our development is based on three unifying principles: (1) there is only one system of balance laws based on continuum mixture theory that describes the physical model, (2) there is only one natural energy-dissipation law that leads to quasi-incompressible NSCH models, (3) variations between the models only appear in the constitutive choices. The framework presented in this work now completes the fundamental exploration of alternate non-matching density NSCH models that utilize a single momentum equation for the mixture velocity, but leaves open room for further sophistication in the energy functional and constitutive dependence.

@en

In this work we present a novel monolithic Finite Element method for the hydroelastic analysis of very large floating structures (VLFS) with arbitrary shapes that is stable, energy conserving, and overcomes the need of an iterative algorithm. The new formulation enables a fully monolithic solution of the linear free-surface flow, described by linear potential flow, coupled with floating thin structures, described by the Euler–Bernoulli beam or Poisson–Kirchhoff plate equations. The formulation presented in this work is general in the sense that solutions can be found in the frequency and time domains, it overcomes the need of using elements with (Formula presented.) continuity by employing a continuous/discontinuous Galerkin approach, and it is suitable for finite elements of arbitrary order. We show that the proposed approach can accurately describe the hydroelastic phenomena of VLFS with a variety of tests, including structures with elastic joints, variable bathymetry, and arbitrary structural shapes.

@en

Turbulent flows past rough surfaces can create substantial energy losses in engineering equipment. During the last decades, developing accurate correlations to predict the thermal and hydrodynamic behavior of rough surfaces has proven to be a difficult challenge. In this work, we develop a convolutional neural network architecture to perform a direct image-to-image translation between the height map of a rough surface and its detailed local drag resistance and heat transfer rates. Various techniques are discussed to improve the computational efficiency of the machine learning architecture proposed, and even to reduce its time and space complexity. The main study is based on a new DNS database formed by 24 flow cases at a friction Reynolds number of Re_τ = 180 obtained by applying a random shift to the Fourier spectrum of the grit-blasted surface scanned by Busse et al. (2015,). The results show that machine learning can accurately predict the global values of the drag resistance and heat fluxes across a rough surface. The local predictions for both momentum and heat transfer also show a considerable improvement upon increasing the dataset size. A detailed analysis of the global skin friction values and Stanton numbers predicted by deep learning further reveals that the results surpass the accuracy of traditional correlations by a substantial margin in the dataset analyzed.

@en

MFEM is an open-source, lightweight, flexible and scalable C++ library for modular finite element methods that features arbitrary high-order finite element meshes and spaces, support for a wide variety of discretization approaches and emphasis on usability, portability, and high-performance computing efficiency. MFEM's goal is to provide application scientists with access to cutting-edge algorithms for high-order finite element meshing, discretizations and linear solvers, while enabling researchers to quickly and easily develop and test new algorithms in very general, fully unstructured, high-order, parallel and GPU-accelerated settings. In this paper we describe the underlying algorithms and finite element abstractions provided by MFEM, discuss the software implementation, and illustrate various applications of the library.

@en

Two well-established classes of the interface capturing models are the level-set and phase-field models. Level-set formulations satisfy the maximum principle for the density but are not energy-stable. On the other hand, phase-field models do satisfy the second law of thermodynamics but lack the maximum principle for the density. In this paper we derive a novel model for incompressible immiscible two-phase flow with non-matching densities and surface tension that is both energetically-stable and satisfies the maximum principle for the density. The model finds its place at the intersection of level-set and phase-field models. Its derivation is based on a diffusification of the incompressible two-phase Navier–Stokes equations with non-matching densities and surface tension and involves functional entropy variables. Additionally, we present an associated fully-discrete energy-stable method. Isogeometric analysis is used for the spatial discretization and the temporal-integration is performed with a new time-integration scheme that is a perturbation of the second-order midpoint scheme. The fully-discrete scheme is unconditionally energy-dissipative, pointwise divergence-free and satisfies the maximum principle for the density. Numerical examples in two and three dimensions verify the energetic-stability of the methodology.

@en

We show that in the variational multiscale framework, the weak enforcement of essential boundary conditions via Nitsche's method corresponds directly to a particular choice of projection operator. The consistency, symmetry and penalty terms of Nitsche's method all originate from the fine-scale closure dictated by the corresponding scale decomposition. As a result of this formalism, we are able to determine the exact fine-scale contributions in Nitsche-type formulations. In the context of the advection–diffusion equation, we develop a residual-based model that incorporates the non-vanishing fine scales at the Dirichlet boundaries. This results in an additional boundary term with a new model parameter. We then propose a parameter estimation strategy for all parameters involved that is also consistent for higher-order basis functions. We illustrate with numerical experiments that our new augmented model mitigates the overly diffusive behavior that the classical residual-based fine-scale model exhibits in boundary layers at boundaries with weakly enforced essential conditions.

@en

Wind propulsion for commercial ships has been identified as a key component in the energy transition for the maritime industry. The sailing hybrid ship will operate with leeway (drift) angles to produce a lateral force known as sideforce, for steady operation under sail. In this paper, experimental results for the sailing performance of ships fitted with bilge keel appendages are presented. Systematic variations in appendage height, length, and position were tested, including several special cases (multiple bilge keels). The appendage typology is shown to mitigate the strong ‘destabilizing’ yaw moment that is characteristic of wind-assisted commercial vessels and to promote the non-linear sideforce component. The working principal for bilge keels—promotion of flow separation—can be employed to specify the separation location for components of the vessel vortex wake to improve the sailing performance of the ship.@en

This paper presents a novel variational formulation to simulate linear free-surface flow. The variational formulation is suitable for higher-order finite elements and higher-order and higher-continuity shape functions as employed in Isogeometric Analysis (IGA). The novel formulation combines the interior and free-surface problems in one monolithic formulation. This leads to exact energy conservation and superior performance in terms of accuracy when compared to a traditional segregated formulation. This is confirmed by the numerical computation of traveling waves in a periodic domain and a three-dimensional sloshing problem. The isogeometric approach shows significant improved performance compared to traditional finite elements. Even on very coarse quadratic NURBS meshes the dispersion error is virtually absent.

@en

What it is Broaching-to? When it can occur? What can happen if it happens? What is known about it? How it can be investigated further? Practical ways to prevent it in design and in operation@en

In this paper we show that the variational multiscale method together with the variation entropy concept form the underlying theoretical framework of discontinuity capturing. The variation entropy [M.F.P. ten Eikelder and I. Akkerman, Comput. Methods Appl. Mech. Engrg. 355 (2019) 261-283] is the recently introduced concept that equips total variation diminishing solutions with an entropy foundation. This is the missing ingredient in order to show that the variational multiscale method can capture sharp layers. The novel framework naturally equips the variational multiscale method with a class of discontinuity capturing operators. This class includes the popular YZβ method and methods based on the residual of the variation-entropy. The discontinuity capturing mechanisms do not contain ad hoc devices and appropriate length scales are derived. Numerical results obtained with quadratic NURBS are virtually oscillation-free and show sharp layers, which confirms the viability of the methodology.

@en

A Reynolds-averaged Navier Stokes computational fluid dynamics (RANS-CFD) package will be one of the primary tools used during the development of a performance prediction program for wind-assisted commercial ships. This paper describes the simulation verification exercise, performed in support of the experimental validation presented in Part 1 of this two-part series describing the RANS-CFD method employed in this research. The predominance of large-scale separated flow structures in the wake of the sailing ship, an artefact of sideforce production necessary for sailing, points to a careful verification exercise and estimate for the numerical uncertainty to support the systematic investigation of wind-assisted ship hydromechanics and meshing guidelines within the available computer resources. Methods for CFD uncertainty quantification are defined and implemented for verification cases at leeway angles equal to 0^ᵒ, 6^ᵒ, and 9^ᵒ. Analysis for four sets of grids with different meshing strategies and for varying time steps results in a grid definition and time step for simulation validation. Numerical uncertainty as adopted in Part 1 for validation is defined. Finally, the meshing strategy for full-scale simulation is described, as used for the production runs of the Delft Wind Assist Series.

@en

Computational fluid dynamics (CFD) and fluid–structure interaction (FSI) are growing disciplines in the aeroelastic analysis and design of long-span bridges, which, with their bluff body characteristics, offer major challenges to efficient simulation. In this paper, we employ isogeometric analysis (IGA) based on nonuniform rational B-splines (NURBS) to numerically simulate turbulent flows over moving bridge sections in three dimensions (3D). Stationary and dynamic analyses of two bridge sections, an idealized rectangular shape with aspect ratio 1∶10 and a 1∶50-scale model of the Hardanger bridge, were performed. Wind tunnel experiments and comparative finite-element (FE) analyses of the same sections were also conducted. Studies on the convergence, static dependencies on the angle of attack, and self-excited forces in terms of the aerodynamic derivatives show that IGA successfully captures the bluff-body flow characteristics and exhibits superior per-degree-of-freedom accuracy compared to the more traditional lower-order FE discretizations.@en

This paper presents the notion of a variation entropy. This concept is an entropy framework for the gradient of the solution of a conservation law instead of on the solution itself. It appears that all semi-norms are admissible variation entropies. This provides insight into the total variation diminishing property and justifies it from entropy principles. The framework allows to derive new local variation diminishing properties in the continuous form. This can facilitate the design of new numerical methods for problems that contain discontinuities.@en

This paper presents a new monolithic free-surface formulation that exhibits correct kinetic and potential energy behavior. We focus in particular on the temporal energy behavior of two-fluids flow with varying densities. Correct energy behavior here means that the actual energy evolution of the numerical solution matches the evolution as predicted by the discrete two–fluid equations. We adopt the level-set method to describe the two-fluid surface. To ensure the correct energy behavior we augment the interface convection equation with kinetic and potential energy constraints. We solve the resulting formulation consisting of the fluid and interface equations in a monolithic fashion using a recently proposed level-set method [26]. For the spatial discretization divergence-conforming NURBS are adopted. The resulting discrete equations are solved with a quasi-newton method which partially decouples the constraints from the rest of the problem. As we focus on the energy behavior of time integration in case of varying densities, we restrict ourselves to low-Reynolds-number flow allowing simple Galerkin discretizations. High-Reynolds-number two-fluid flows that require stabilization are beyond the scope of the current paper. The simulation of a dambreak problem numerically supports the correct energy behavior of the proposed methodology. The proposed methodology improves the solution quality significantly upon a more traditional approach. Due to the excellent accuracy per degree of freedom one can suffice with a much lower resolution.@en

This paper presents the construction of a correct-energy stabilized finite element method for the incompressible Navier–Stokes equations. The framework of the methodology and the correct-energy concept have been developed in the convective–diffusive context in the preceding paper [M.F.P. ten Eikelder, I. Akkerman, Correct energy evolution of stabilized formulations: The relation between VMS, SUPG and GLS via dynamic orthogonal small-scales and isogeometric analysis. I: The convective–diffusive context, Comput. Methods Appl. Mech. Engrg. 331 (2018) 259–280]. The current work extends ideas of the preceding paper to build a stabilized method within the variational multiscale (VMS) setting which displays correct-energy behavior. Similar to the convection–diffusion case, a key ingredient is the proper dynamic and orthogonal behavior of the small-scales. This is demanded for correct energy behavior and links the VMS framework to the streamline-upwind Petrov–Galerkin (SUPG) and the Galerkin/least-squares method (GLS). The presented method is a Galerkin/least-squares formulation with dynamic divergence-free small-scales (GLSDD). It is locally mass-conservative for both the large- and small-scales separately. In addition, it locally conserves linear and angular momentum. The computations require and employ NURBS-based isogeometric analysis for the spatial discretization. The resulting formulation numerically shows improved energy behavior for turbulent flows comparing with the original VMS method.

@en