The open Delft Advanced Research Terra Simulator (open-DARTS) framework is an open-source reservoir simulation software. The open-DARTS focused on energy transition applications, such as geothermal energy production and carbon sequestration. It enables the modeling of compositional thermal flow, coupled with a geomechanical solver based on the Finite Volume discretization and adjoints method for inverse modeling. The open-DARTS supports different grid types (structured, corner-point geometry, and unstructured), discrete fracture networks, contact mechanics, and various thermal-chemical interactions. The recently proposed generic nonlinear formulation supports the most general nonlinear PDEs designed for various energy transition applications. The open-DARTS has been implemented in C++ and Python to optimize hardware utilization while ensuring flexibility. The most computationally expensive part is written in C++ and compiled into libraries, which are subsequently exposed to Python using pybind11. This allows the extension and overriding of C++ functions by user-defined Python code. For example, using only a Python interface, one can adjust a timestep strategy, nonlinear solver, or properties output. Besides, the Python interface of open-DARTS provides straightforward coupling with other Python-based numerical modeling packages, including the meshing, file storage, caching, and visualization modules. The open-DARTS core uses the advantages of C++ language, such as efficient low-level memory management, object-oriented programming, compile-time polymorphism, and parallelization with OpenMP. One of the advantages of open-DARTS is the Operator-Based Linearization (OBL) technique, which can resolve challenges associated with complex physics and reduce the computation time, especially for ensemble-based simulations. We would also like to share our experience on the project, repository, and the development workflow configuration using gitlab.com, including the build system (cmake), handling merge requests, automated testing in CI/CD pipelines, documentation management (gitlab.io), wiki utilization, and release publishing. Additionally, Python’s integration into open-DARTS offers the advantage of straightforward installation via PyPI and simplifies defining requirements for users who prefer to avoid compiling code from source files.

In this work, we present a mass, energy, enstrophy and vorticity conserving (MEEVC) mixed finite element discretization for two-dimensional incompressible Navier-Stokes equations as an alternative to the original MEEVC scheme proposed in A. Palha and M. Gerritsma (2017) [5]. The present method can incorporate general boundary conditions. Conservation properties are proven. Supportive numerical experiments with both exact and inexact quadrature are provided.

In this work we use algebraic dual spaces with a domain decomposition method to solve the Darcy equations. We define the broken Sobolev spaces and their finite dimensional counterparts. A global trace space is defined that connects the solution between the broken spaces. Use of algebraic dual spaces results in a sparse, metric-free representation of the incompressibility constraint, the pressure gradient term, and on the continuity constraint between the sub domains. To demonstrate this, we solve two test cases: (i) a manufactured solution case, and (ii) an industrial benchmark reservoir modelling problem SPE10. The results demonstrate that the dual spaces can be used for domain decomposition formulation, and despite having more unknowns, requires less simulation time compared to the continuous Galerkin formulation, without compromising on the accuracy of the solution.

We introduce a mimetic dual-field discretization which conserves mass, kinetic energy and helicity for three-dimensional incompressible Navier-Stokes equations. The discretization makes use of a conservative dual-field mixed weak formulation where two evolution equations of velocity are employed and dual representations of the solution are sought for each variable. A temporal discretization, which staggers the evolution equations and handles the nonlinearity such that the resulting discrete algebraic systems are linear and decoupled, is constructed. The spatial discretization is mimetic in the sense that the finite dimensional function spaces form a discrete de Rham complex. Conservation of mass, kinetic energy and helicity in the absence of dissipative terms is proven at the discrete level. Proper dissipation rates of kinetic energy and helicity in the viscous case are also proven. Numerical tests supporting the method are provided.

In ℝⁿ, let Λ^k(Ω) represent the space of smooth differential k-forms in Ω. The de Rham complex consists of a sequence of spaces, Λ^k(Ω), k = 0, 1…, n, connected by the exterior derivative, d: Λ^k(Ω) → Λ^k+1(Ω). Appropriately chosen B-spline spaces together with their associated dual B-spline spaces form a discrete double de Rham complex. In practical applications, this discrete double de Rham complex leads to very sparse systems. In this paper, this construction will be explained and illustrated by means of a non-trivial three-dimensional example.

A new horizontally explicit/vertically implicit (HEVI) time splitting scheme for atmospheric modelling is introduced, for which the horizontal divergence terms are applied within the implicit vertical substep. The new HEVI scheme is implemented in conjunction with a mixed mimetic spectral element spatial discretisation and semi-implicit vertical time stepping scheme that both preserve the skew-symmetric structure of the non-canonical Hamiltonian form of the equations of motion. Within this context the new HEVI scheme allows for the exact balance of all energetic exchanges in space and time. However since the choice of horizontal fluxes for which this balance is satisfied is not consistent with the horizontal velocity at the end of the time level the scheme still admits a temporal energy conservation error. Linearised eigenvalue analysis shows that similar to a fully implicit method, the new HEVI scheme is neutrally stable for all buoyancy modes, and unlike a second order trapezoidal HEVI scheme is stable for all acoustic modes below a certain horizontal CFL number. The scheme is validated against standard test cases for both planetary and non-hydrostatic regimes. For the planetary scale baroclinic instability test case, the new formulation exhibits a secondary oscillation in the potential to kinetic energy power exchanges, with a temporal frequency approximately four times that exhibited by a horizontally third order, vertically second order trapezoidal scheme. For the non-hydrostatic test case, the vertical upwinding of the potential temperature diagnostic equation is shown to reduce spurious oscillations without altering the energetics of the solution, since this upwinding is performed in an energetically consistent manner. For this test case, which is configured on an affine geometry, the exact balance of energy exchanges allows the model to run stably without any form of dissipation.

A model of the three-dimensional rotating compressible Euler equations on the cubed sphere is presented. The model uses a mixed mimetic spectral element discretization which allows for the exact exchanges of kinetic, internal and potential energy via the compatibility properties of the chosen function spaces. A Strang carryover dimensional splitting procedure is used, with the horizontal dynamics solved explicitly and the vertical dynamics solved implicitly so as to avoid the CFL restriction of the vertical sound waves. The function spaces used to represent the horizontal dynamics are discontinuous across vertical element boundaries, such that each horizontal layer is solved independently so as to avoid the need to invert a global 3D mass matrix, while the function spaces used to represent the vertical dynamics are similarly discontinuous across horizontal element boundaries, allowing for the serial solution of the vertical dynamics independently for each horizontal element. The model is validated against standard test cases for baroclinic instability within an otherwise hydrostatically and geostrophically balanced atmosphere, and a non-hydrostatic gravity wave as driven by a temperature perturbation.

In this paper, we present a hybrid mimetic method which solves the mixed formulation of the Poisson problem on curvilinear quadrilateral meshes. The method is hybrid in the sense that the domain is decomposed into multiple disjoint elements and the interelement continuity is enforced using a Lagrange multiplier. The method is mimetic in the sense that the discrete divergence operator is exact. By using the mimetic basis functions and their dual representations, various metric-free discrete terms are obtained. The discrete system can be efficiently solved by first solving a reduced system for the Lagrange multiplier. Numerical experiments which validate the method are presented.

In this paper, we will show that the equivalence of a div-grad Neumann problem and a grad-div Dirichlet problem can be preserved at the discrete level in 3-dimensional curvilinear domains if algebraic dual polynomial representations are employed. These representations will be introduced. Proof of the equivalence at the discrete level follows from the construction of the algebraic dual representations. A 3-dimensional test problem in curvilinear coordinates will illustrate this approach.

In this paper we will consider two curl-curl equations in two dimensions. One curl-curl problem for a scalar quantity F and one problem for a vector field E. For Dirichlet boundary conditions n× E= Ê_⊣ on E and Neumann boundary conditions n×curlF=Ê⊣, we expect the solutions to satisfy E = curl F. When we use algebraic dual polynomial representations, these identities continue to hold at the discrete level. Equivalence will be proved and illustrated with a computational example.

Given a sequence of finite element spaces which form a de Rham sequence, we will construct dual representations of these spaces with associated differential operators which connect these spaces such that they also form a de Rham sequence. The dual representations also need to satisfy the de Rham sequence on the domain boundary. The matrix which converts primal representations to dual representations – the Hodge matrix – is the mass or Gram matrix. It will be shown that a bilinear form of a primal and a dual representation is equal to the vector inner product of the expansion coefficients (degrees of freedom) of both representations. This leads to very sparse system matrices, even for high order methods. The differential operators for dual representations will be defined. Vector operations, grad, curl and div, for primal and dual representations are both topological and do not depend on the metric, i.e. the size and shape of the mesh or the order of the numerical method. Derivatives are evaluated by applying sparse incidence and inclusion matrices to the expansion coefficients of the representations. As illustration of the use of dual representations, the method will be applied to (i) a mixed formulation for the Poisson problem in 3D, (ii) it will be shown that this approach allows one to preserve the equivalence between Dirichlet and Neumann problems in the finite dimensional setting and, (iii) the method will be applied to the approximation of grad–div eigenvalue problem on affine and non-affine meshes.

We present a hybrid mimetic spectral element formulation for Darcy flow. The discrete representations for (1) conservation of mass, and (2) inter-element continuity, are topological relations that lead to sparse matrix systems. These constraints are independent of the element size and shape, and thus invariant under mesh transformations. The resultant algebraic system is extremely sparse even for high degree polynomial basis. Furthermore, the system can be efficiently assembled and solved for each element separately.

The research program of the TCV tokamak ranges from conventional to advanced-tokamak scenarios and alternative divertor configurations, to exploratory plasmas driven by theoretical insight, exploiting the device's unique shaping capabilities. Disruption avoidance by real-time locked mode prevention or unlocking with electron-cyclotron resonance heating (ECRH) was thoroughly documented, using magnetic and radiation triggers. Runaway generation with high-Z noble-gas injection and runaway dissipation by subsequent Ne or Ar injection were studied for model validation. The new 1 MW neutral beam injector has expanded the parameter range, now encompassing ELMy H-modes in an ITER-like shape and nearly non-inductive H-mode discharges sustained by electron cyclotron and neutral beam current drive. In the H-mode, the pedestal pressure increases modestly with nitrogen seeding while fueling moves the density pedestal outwards, but the plasma stored energy is largely uncorrelated to either seeding or fueling. High fueling at high triangularity is key to accessing the attractive small edge-localized mode (type-II) regime. Turbulence is reduced in the core at negative triangularity, consistent with increased confinement and in accord with global gyrokinetic simulations. The geodesic acoustic mode, possibly coupled with avalanche events, has been linked with particle flow to the wall in diverted plasmas. Detachment, scrape-off layer transport, and turbulence were studied in L- and H-modes in both standard and alternative configurations (snowflake, super-X, and beyond). The detachment process is caused by power 'starvation' reducing the ionization source, with volume recombination playing only a minor role. Partial detachment in the H-mode is obtained with impurity seeding and has shown little dependence on flux expansion in standard single-null geometry. In the attached L-mode phase, increasing the outer connection length reduces the in-out heat-flow asymmetry. A doublet plasma, featuring an internal X-point, was achieved successfully, and a transport barrier was observed in the mantle just outside the internal separatrix. In the near future variable-configuration baffles and possibly divertor pumping will be introduced to investigate the effect of divertor closure on exhaust and performance, and 3.5 MW ECRH and 1 MW neutral beam injection heating will be added.

This work presents three methods for enforcing tangential velocity boundary conditions for the MEEVC scheme, which was shown to be mass, enstrophy, energy and vorticity conserving scheme in the case of inviscid flow [1]. While the normal velocity component can be strongly imposed in a div-conforming formulation for the velocity field, inclusion of the tangential velocity needs to be set through an appropriate choice of vorticity boundary conditions. Three methods to impose the tangential velocity boundary condition will be discussed: The kinematic Dirichlet formulation, the kinematic Neumann formulation and the dynamic Neumann formulation. The conservation properties of each of the resulting schemes are analyzed and numerical results are shown for the Taylor–Green vortex and for the dipole collision test cases. These confirm that kinematic Neumann vorticity boundary conditions perform best.

Within the EUROfusion MST1 work package, a series of experiments has been conducted on AUG and TCV devices to disentangle the role of plasma fueling and plasma shape for the onset of small ELM regimes. On both devices, small ELM regimes with high confinement are achieved if and only if two conditions are fulfilled at the same time. Firstly, the plasma density at the separatrix must be large enough (n_e,sep/n_G ∼ 0.3), leading to a pressure profile flattening at the separatrix, which stabilizes type-I ELMs. Secondly, the magnetic configuration has to be close to a double null (DN), leading to a reduction of the magnetic shear in the extreme vicinity of the separatrix. As a consequence, its stabilizing effect on ballooning modes is weakened.

In this paper, we will use algebraic dual polynomials to set up a discrete Steklov-Poincaré operator for the mixed formulation of the Poisson problem. The method will be applied in curvilinear coordinates and to a test problem which contains a singularity. Exponential convergence of the trace variable in H 1 / 2 {H^{1/2}} -norm will be shown.

A mixed mimetic spectral element method is applied to solve the rotating shallow water equations. The mixed method uses the recently developed spectral element histopolation functions, which exactly satisfy the fundamental theorem of calculus with respect to the standard Lagrange basis functions in one dimension. These are used to construct tensor product solution spaces which satisfy the generalized Stokes theorem, as well as the annihilation of the gradient operator by the curl and the curl by the divergence. This allows for the exact conservation of first order moments (mass, vorticity), as well as higher moments (energy, potential enstrophy), subject to the truncation error of the time stepping scheme. The continuity equation is solved in the strong form, such that mass conservation holds point wise, while the momentum equation is solved in the weak form such that vorticity is globally conserved. While mass, vorticity and energy conservation hold for any quadrature rule, potential enstrophy conservation is dependent on exact spatial integration. The method possesses a weak form statement of geostrophic balance due to the compatible nature of the solution spaces and arbitrarily high order spatial error convergence.

This chapter addresses the topological structure of steady, anisotropic, inhomogeneous diffusion problems. Differential operators are represented by sparse incidence matrices, while weighted mass matrices play the role of metric-dependent Hodge matrices. The resulting mixed formulation is point-wise divergence-free if the right hand side function f = 0. The method is inf-sup stable; no stabilization is required and the method displays optimal convergence on orthogonal and deformed grids.

One of the most cited disadvantages of least-squares formulations is its lack of conservation. By a suitable choice of least-squares functional and the use of appropriate conforming finite dimensional function spaces, this drawback can be completely removed. Such a mimetic least-squares method is applied to a curl-curl system. Conservation properties will be proved and demonstrated by test results on two-dimensional curvilinear grids.