Heat equation with Delta Dirac function as boundary condition

[dokken]

Are you now trying to have a single dirac delta function? If you have a dirac delta function at a single point, you could have a look at:

Point Sources (Redux)

Following my post above, with a minor adaptations, a point source can be implemented as follows: # Create a point source for Poisson problem # Author: Jørgen S. Dokken # SPDX-License-Identifier: MIT from mpi4py import MPI from petsc4py import PETSc import dolfinx import dolfinx.fem.petsc import numpy as np import ufl def compute_cell_contributions(V, points): # Determine what process owns a point and what cells it lies within mesh = V.mesh _, _, owning_points, cells = dolfinx.cp…

that uses a single point source in either dolfin or dolfinx.