Hi, I am currently using a RK solver to solve the ST. Venant equations for a square channel given by \frac{\partial A}{\partial t}+\frac{\partial Q}{\partial s}=0 and \frac{\partial Q}{\partial t} + \frac{\partial}{\partial s}(\frac{Q^2}{A}+g A^2/(2 w)) + gA(S_f -S)=0 Here Q is the discharge, A…

This is a TVD limited code that is relatively quick, but only works for p=1: from mpi4py import MPI from petsc4py import PETSc import numpy as np import ufl from dolfinx import default_real_type, fem, has_adios2, io, mesh from dolfinx.plot import vtk_mesh from dolfinx.fem.petsc import LinearProblem…

One of the few projects I know which has implemented slope limiters for DG methods with FEniCS is Ocellaris. See, e.g., [image] Slope limiting for DG variational formulation Hi, has anyone attempted to use slope limiters for DG schemes in fenics? Right now I’…

I got it working now, but massively decreased the efficiency of the code. I will try to optimize it and once it works decently upload it here so people can use it as a reference in case they ever need it. [image]

Slope limiting for the St.Venant equations in DOLFINx

dolfinx

nate May 18, 2026, 3:25pm 2

One of the few projects I know which has implemented slope limiters for DG methods with FEniCS is Ocellaris. See, e.g.,

Topic		Replies	Views
Slope limiting for DG variational formulation	10	876	November 18, 2024
Does Dolfinx's solver have an issue or require some preconditioning when handling certain meshes like the Icosahedral mesh? General dolfinx , legacy-dolfin	10	178	September 6, 2025
Mixed formulation of stokes equations dolfinx	12	678	June 23, 2024
Solve unsteady NS General dolfinx	16	688	February 20, 2024
Imposing a discontinuity at interface using DG method dolfinx	27	2361	January 19, 2024

Slope limiting for the St.Venant equations in DOLFINx

Related topics