Unable to compute the curl of a vector field

dokken · June 25, 2024, 8:37pm

You are using a 2D mesh.
The ufl curl operator is only defined in 3D, as written in the documentation:
https://docs.fenicsproject.org/ufl/2024.1.0/manual/form_language.html#curl-and-rot

See for instance

github.com

FEniCS/dolfinx/blob/9f33dabc270b8e3c31e815515ecac5eb2a5414b0/python/demo/demo_pml.py#L378


      
              ky = n_b * k0 * np.sin(theta)
              phi = kx * x[0] + ky * x[1]
              return (-np.sin(theta) * np.exp(1j * phi), np.cos(theta) * np.exp(1j * phi))
          
          
          # -
          
          # For convenience, we define the $\nabla\times$ operator for a 2D vector
          
          
          def curl_2d(a: fem.Function):
              return ufl.as_vector((0, 0, a[1].dx(0) - a[0].dx(1)))
          
          
          # Let's now see how we can implement PMLs for our problem. PMLs are
          # artificial layers surrounding the real domain that gradually absorb
          # waves impinging them. Mathematically, we can use a complex coordinate
          # transformation of this kind to obtain this absorption:
          #
          # $$
          # x^\prime= x\left\{1+j\frac{\alpha}{k_0}\left[\frac{|x|-l_{dom}/2}

or Vector curl in 2D - #2 by nate
for solutions

Topic		Replies	Views
Postprocessing Nedelec element solution General	1	236	September 22, 2023
Cannot solve curl of Lagrange vectorfunctionspace variational formulation	2	855	August 6, 2019
Calculating the curl of a tensor field mathematics	4	722	August 1, 2022
Possible to calculate "H = "curl(E)" without project command or solving SLE? variational formulation	2	339	April 2, 2020
Regarding computation of the error norm General	4	768	December 22, 2022

Unable to compute the curl of a vector field

Related topics