The paper presents a monolithic finite element model for the hydro-visco-elastic analysis of floating membranes interacting with ocean waves. The formulation couples linearised potential flow and viscoelastic membrane equations, offering a versatile tool for modelling arbitrarily shaped floating membranes in varying sea-bed topography. The paper also presents a wet modal analysis for the coupled problem, accounting for the added mass and stiffness of the surrounding fluid. This model is used to study the dependence of the wet natural frequencies of floating membranes on the material properties. It is also used to analyse the reflection, transmission, scattering and absorption of ocean wave energy by 1D and 2D floating membranes. Notably, the paper underscores the impact of proportional material damping on these observed phenomena. The results highlight local peaks in the viscoelastic behaviour at the calculated wet natural frequencies, and demonstrate the outward dispersion of incoming wave around finite 2D membranes. Furthermore, the model is employed to examine the interaction of viscoelastic membranes with other structures, such as a monopile, under the influence of ocean waves. This comprehensive investigation contributes to a deeper understanding of the fluid–structure interaction inherent to certain floating solar, wave-energy converter and floating breakwater technologies.@en

The paper presents a monolithic finite element model for the hydro-visco-elastic analysis of floating membranes interacting with ocean waves. The formulation couples linearised potential flow and viscoelastic membrane equations, offering a versatile tool for modelling arbitrarily shaped floating membranes in varying sea-bed topography. The paper also presents a wet modal analysis for the coupled problem, accounting for the added mass and stiffness of the surrounding fluid. This model is used to study the dependence of the wet natural frequencies of floating membranes on the material properties. It is also used to analyse the reflection, transmission, scattering and absorption of ocean wave energy by 1D and 2D floating membranes. Notably, the paper underscores the impact of proportional material damping on these observed phenomena. The results highlight local peaks in the viscoelastic behaviour at the calculated wet natural frequencies, and demonstrate the outward dispersion of incoming wave around finite 2D membranes. Furthermore, the model is employed to examine the interaction of viscoelastic membranes with other structures, such as a monopile, under the influence of ocean waves. This comprehensive investigation contributes to a deeper understanding of the fluid–structure interaction inherent to certain floating solar, wave-energy converter and floating breakwater technologies.@en

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

@en

Offshore wind farms, critical for sustainable energy production, face the challenge of optimization among many parameters influencing key performance indicators in competitive ways. This research introduces the novel Integrative Maximized Aggregated Preference Wind Farm Optimization (IMAP-WFO) framework – a comprehensive tool designed to enhance flexibility, accuracy, and uncertainty quantification in offshore wind farm design and operation. Existing methods often fall short due to limitations in adaptability and precision, especially when modeling complex multi-physical behaviors under uncertain conditions. IMAP-WFO overcomes these limitations by combining advanced statistical techniques and simulation methods. At its core are parametric design performance functions, capturing critical aspects of wind farm behavior, including energy production, material usage, and structural fatigue. These functions rely on Kriging meta-models. To address inherent uncertainty, Monte Carlo simulations provide a probabilistic assessment of outcomes. IMAP-WFO's true innovation lies in translating technical functions into socio-economic objectives, including sustainability metrics, annual energy production, capital expenditure, operational expenditure, model uncertainty, and lifetime fatigue. Stakeholders can dynamically weigh these objectives based on their preferences. A validation process ensures the accuracy of design performance functions, comparing simulated results with real-world data. IMAP-WFO's application is demonstrated through case studies: optimizing the levelized cost of energy and exploring wind farm control strategies.@en

This work aims to develop a low-fidelity model for a lattice support structure for offshore wind applications. The proposed low-fidelity model consists of a sequence of regular Timoshenko beams, each of them characterized by homogenized mechanical and mass properties representative of the single bays of the reference space-frame structure. The homogenized elastic coefficients of the sequence of beams are then computed by means of two alternative procedures: case (a), via analytical expressions available in the literature and accounting for a partially isotropic behaviour; case (b) by means of an optimization procedure, with ad hoc calibration factors. The suggested methods to derive the homogenized elastic coefficients are then tested for both straight and tapered lattice structures. The prediction performance is evaluated in terms of estimation of the first five natural frequencies and mode shapes, response to dynamic loads, and ability to predict rotor-structure interaction phenomena. A parametric study is then performed to evaluate the potential and limitations of the proposed models. To bypass the optimization procedure (b), a data-driven approach is also proposed for the case of straight lattice structures. Overall, the developed low-fidelity model leads to a computational speed-up factor of at least 60. The prediction reliability of the low-fidelity model is discussed for a tapered and regular straight lattice structure. However, for the latter one, a more detailed comparative study between the various modelling assumptions is performed and discussed. With reference to the straight lattice tower, whenever an optimization procedure is used (case (b)), and with reference to a typical subset of the investigated geometrical parameter space, the mean prediction error of the first five natural frequencies is lower than 1%. On the other hand, for case (a) and for the same investigated subset, the mean prediction errors for the first two bending modes and the torsional mode are, 5.2%, 13.3% and 18.8%, respectively. These results are improved in case a data-driven regression model is used to predict the calibration factors, leading to mean prediction errors below 5% for the entire investigated parameter space.

@en

Floating wind farms are a promising solution for offshore wind energy production in deep waters. However, the design optimisation process of these farms is difficult due to their complex and multidisciplinary nature. Furthermore, current optimisation methods: 1) ignore and/or provide no insight into the dynamic interplay between the preference-dominated management domain and the object-performance-dominated engineering domain; 2) are limited to evaluating potentially sub-optimal design alternatives; 3) contain fundamental aggregation modelling errors; 4) do not return a single optimal design point. This paper presents an optimisation framework that overcomes these shortcomings and enables truly integrative multi-objective design optimisation. It includes a surrogate model that interacts with the wind turbine simulation tool OpenFAST to enable preliminary design of the structure’s mooring system. Applied to a demonstration project and validated against real projects in a maritime contractor environment, the workflow shows improvements in tender performance and added value over single-sided cost optimisations.@en

The hydrofoil harnesses wave energy and converts it into thrust. In this paper, we present the results of the first experimental study investigating the dynamic behavior of a fully passive foil with spring-loaded pitch and heave in regular waves. Our study shows that the real-time load signal is multi-harmonic with strong superposition, directly proving the robust energy harvesting performance due to the restoring springs. By interpreting the hydrofoil's pose and path from an image sequence captured underwater, we conclude the dynamic evolution of the fully passive hydrofoil interacting with regular waves. The hydrofoil's dynamics exhibit asymmetric surge, pitch, and heave in a motion cycle. Furthermore, we employ a pixel capturing algorithm with self-correction utility to quantify the hydrofoil's forward displacement from the image sequence of the moving carriage. These findings provide valuable insight into the performance and potential of hydrofoils for marine propulsion.

@en

Flow simulations on porous media, reconstructed from Micro-Computerised Tomography (μCT) scans, is becoming a common tool to compute the permeability of rocks. Still, some conditions need to be met to obtain accurate results. Only if the sample size is equal or larger than the Representative Elementary Volume will the computed effective permeability be representative of the rock at a continuum scale. Moreover, the numerical discretisation of the digital rock needs to be fine enough to reach numerical convergence. In the particular case of using Finite Elements (FE) and cartesian meshes, studies have shown that the meshes should be at least two times finer than the original image resolution in order to reach the simulation's mesh convergence. These two conditions and the increased resolution of μCT-scans to observe finer details of the microstructure, can lead to extremely computationally expensive numerical simulations. In order to reduce this cost, we couple a FE numerical model for Stokes flow in porous media with an unfitted boundary method for cartesian meshes, which allows to improve results precision for coarse meshes. Indeed, this method enables to obtain a definition of the pore–grain interface as precise as for a conformal mesh, without a computationally expensive and complex mesh generation for μCT-scans of rocks. From the benchmark of three different rock samples, we observe a clear improvement of the mesh convergence for the permeability value using the unfitted boundary method on cartesian meshes. An accurate permeability value is obtained for a mesh coarser than the initial image resolution. The method is then applied to a large sample of a high-resolution μCT-scan to showcase its advantage.

@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

We propose a Lagrangian solid mechanics framework for the simulation of salt tectonics and other large-deformation geomechanics problems at the basin scale. Our approach relies on general elastic-viscoplastic constitutive models to characterize the deformation of geologic strata, in contrast with the majority of published works on the subject, which utilize nonlinear Stokes flow models. By means of multiscale asymptotics, we also show that the inertia term in the momentum balance equation can be safely neglected, if the goal is to track the Earth's crust deformation over long periods of time. Our time integration strategy is a blended transient/quasistatic approach, in that it consists of a constitutive stress update, subject to the constraint that the stresses must satisfy static equilibrium. In addition, we use stabilized finite element methods specifically built for triangular and tetrahedral grids, which can also perform well under incompressibility constraints. Our approach offers computational geologists the following advantages: (1) improved flexibility in the choice of subsurface constitutive models with respect to the nonlinear Stokes flow; (2) improved efficiency over transient dynamics algorithms used in this context in the past, which are forced to resolve seismic events over geologic time scales; and (3) improved robustness in large strain computations over quadrilateral/hexahedral finite elements. We demonstrate the performance of the proposed approach with simulations of passive diapirism.

@en

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. We present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

@en

The discontinuous Galerkin (DG) method has found widespread application in elliptic problems with rough coefficients, of which the Darcy flow equations are a prototypical example. One of the long-standing issues of DG approximations is the overall computational cost, and many different strategies have been proposed, such as the variational multiscale DG method, the hybridizable DG method, the multiscale DG method, the embedded DG method, and the Enriched Galerkin method. In this work, we propose a mixed dual-scale Galerkin method, in which the degrees-of-freedom of a less computationally expensive coarse-scale approximation are linked to the degrees-of-freedom of a base DG approximation. We show that the proposed approach has always similar or improved accuracy with respect to the base DG method, with a considerable reduction in computational cost. For the specific definition of the coarse-scale space, we consider Raviart–Thomas finite elements for the mass flux and piecewise-linear continuous finite elements for the pressure. We provide a complete analysis of stability and convergence of the proposed method, in addition to a study on its conservation and consistency properties. We also present a battery of numerical tests to verify the results of the analysis, and evaluate a number of possible variations, such as using piecewise-linear continuous finite elements for the coarse-scale mass fluxes.@en

The numerical simulation of physical phenomena and engineering problems can be affected by numerical errors and various types of uncertainties. Characterizing the former in computational frameworks involving system parameter uncertainties becomes a key issue. In this work, we study the behavior of new variational multiscale (VMS) error estimators for the propagation of parametric uncertainties in a Convection–Diffusion–Reaction (CDR) problem. A sensitivity analysis is performed to assess the performance of the error estimator with respect to the mesh discretization and physical parameters (here, the viscosity value and advection velocity). Three different manufactured analytical solutions are considered as benchmarking tests. Next, the performance of the VMS error estimators is evaluated for the CDR problem with uncertain input parameters. For this purpose, two probabilistic models are constructed for the viscosity and advection direction, and the uncertainties are propagated using a polynomial chaos expansion approach. A convergence analysis is specifically carried out for different configurations, corresponding to regimes where the CDR operator is either smooth or non-smooth. An assessment of the proposed error estimator is finally conducted for the three tests, considering both the viscous- and convection-dominated regimes.

@en

In this work, we define a family of explicit a posteriori error estimators for Finite Volume methods in computational fluid dynamics. The proposed error estimators are inspired by the Variational Multiscale method, originally defined in a Finite Element context. The proposed error estimators are tested in simulations of the incompressible Navier-Stokes equations, the thermally-coupled Navier-Stokes equations, and the fully-coupled compressible large eddy simulation of the HIFiRE Direct Connect Rig Scramjet combustor.

@en

In this article, we develop a dynamic version of the variational multiscale (D-VMS) stabilization for nearly/fully incompressible solid dynamics simulations of viscoelastic materials. The constitutive models considered here are based on Prony series expansions, which are rather common in the practice of finite element simulations, especially in industrial/commercial applications. Our method is based on a mixed formulation, in which the momentum equation is complemented by a pressure equation in rate form. The unknown pressure, displacement, and velocity are approximated with piecewise linear, continuous finite element functions. To prevent spurious oscillations, the pressure equation is augmented with a stabilization operator specifically designed for viscoelastic problems, in that it depends on the viscoelastic dissipation. We demonstrate the robustness, stability, and accuracy properties of the proposed method with extensive numerical tests in the case of linear and finite deformations.

@en

In this work, we develop a high-performance numerical framework for the large eddy simulation (LES) of incompressible flows. The spatial discretization of the nonlinear system is carried out using mixed finite element (FE) schemes supplemented with symmetric projection stabilization of the convective term and a penalty term for the divergence constraint. These additional terms introduced at the discrete level have been proved to act as implicit LES models. In order to perform meaningful wall-unresolved simulations, we consider a weak imposition of the boundary conditions using a Nitsche's-type scheme, where the tangential component penalty term is designed to act as a wall law. Next, segregated Runge–Kutta (SRK) schemes (recently proposed by the authors for laminar flow problems) are applied to the LES simulation of turbulent flows. By the introduction of a penalty term on the trace of the acceleration, these methods exhibit excellent stability properties for both implicit and explicit treatment of the convective terms. SRK schemes are excellent for large-scale simulations, since they reduce the computational cost of the linear system solves by splitting velocity and pressure computations at the time integration level, leading to two uncoupled systems. The pressure system is a Darcy-type problem that can easily be preconditioned using a traditional block-preconditioning scheme that only requires a Poisson solver. At the end, only coercive systems have to be solved, which can be effectively preconditioned by multilevel domain decomposition schemes, which are both optimal and scalable. The framework is applied to the Taylor–Green and turbulent channel flow benchmarks in order to prove the accuracy of the convection-stabilized mixed FEs as LES models and SRK time integrators. The scalability of the preconditioning techniques (in space only) has also been proven for one step of the SRK scheme for the Taylor–Green flow using uniform meshes. Moreover, a turbulent flow around a NACA profile is solved to show the applicability of the proposed algorithms for a realistic problem.

@en

The variational multiscale method thought as an implicit large eddy simulation model for turbulent flows has been shown to be an alternative to the widely used physical-based models. This method is traditionally combined with equal-order velocity-pressure pairs, since it provides pressure stabilization. In this work, we consider a different approach, based on inf-sup stable elements and convection-only stabilization. In order to do so, we consider a symmetric projection stabilization of the convective term using an orthogonal subscale decomposition. The accuracy and efficiency of this method compared with residual-based algebraic subgrid scales and orthogonal subscales methods for equal-order interpolation is assessed in this paper. Moreover, when inf-sup stable elements are used, the grad-div stabilization term has been shown to be essential to guarantee accurate solutions. Hence, a study of the influence of such term in the large eddy simulation of turbulent incompressible flows is also performed. Furthermore, a recursive block preconditioning strategy has been considered for the resolution of the problem with an implicit treatment of the projection terms. Two different benchmark tests have been solved: the Taylor-Green Vortex flow with Re=1600, and the Turbulent Channel Flow at Re_τ=395 and Re_τ=590.

@en

In this work we study the performance of some variational multiscale models (VMS) in the large eddy simulation (LES) of turbulent flows. We consider VMS models obtained by different subgrid scale approximations which include either static or dynamic subscales, linear or nonlinear multiscale splitting, and different choices of the subscale space. After a brief review of these models, we discuss some implementation aspects particularly relevant to the simulation of turbulent flows, namely the use of a skew symmetric form of the convective term and the computation of projections when orthogonal subscales are used. We analyze the energy conservation (and numerical dissipation) of the alternative VMS formulations, which is numerically evaluated. In the numerical study, we have considered three well known problems: the decay of homogeneous isotropic turbulence, the Taylor-Green vortex problem and the turbulent flow in a channel. We compare the results obtained using different VMS models, paying special attention to the effect of using orthogonal subscale spaces. The VMS results are also compared against classical LES scheme based on filtering and the dynamic Smagorinsky closure. Altogether, our results show the tremendous potential of VMS for the numerical simulation of turbulence. Further, we study the sensitivity of VMS to the algorithmic constants and analyze the behavior in the small time step limit. We have also carried out a computational cost comparison of the different formulations. Out of these experiments, we can state that the numerical results obtained with the different VMS formulations (as far as they converge) are quite similar. However, some choices are prone to instabilities and the results obtained in terms of computational cost are certainly different. The dynamic orthogonal subscales model turns out to be best in terms of efficiency and robustness.

@en

In this work, we propose Runge-Kutta time integration schemes for the incompressible Navier-Stokes equations with two salient properties. First, velocity and pressure computations are segregated at the time integration level, without the need to perform additional fractional step techniques that spoil high orders of accuracy. Second, the proposed methods keep the same order of accuracy for both velocities and pressures. The segregated Runge-Kutta methods are motivated as an implicit-explicit Runge-Kutta time integration of the projected Navier-Stokes system onto the discrete divergence-free space, and its re-statement in a velocity-pressure setting using a discrete pressure Poisson equation. We have analysed the preservation of the discrete divergence constraint for segregated Runge-Kutta methods and their relation (in their fully explicit version) with existing half-explicit methods. We have performed a detailed numerical experimentation for a wide set of schemes (from first to third order), including implicit and IMEX integration of viscous and convective terms, for incompressible laminar and turbulent flows. Further, segregated Runge-Kutta schemes with adaptive time stepping are proposed.

@en