Issues Solving Pure Neumann Problems in dolfinx

dokken · February 17, 2023, 11:09pm

You can remove the constant nullspace with PETSc. See for instance

ComputationalPhysiology/oasisx/blob/main/src/oasisx/fracstep.py#L269


      
          self._M.assemble()
          _fem.petsc.assemble_matrix(self._K, self._stiffness_Vi)
          self._K.assemble()
          
          
# Assemble stiffness matrix with boundary conditions
          _fem.petsc.assemble_matrix(self._Ap, self._stiffness_Q, bcs=[
              bc._bc for bc in self._bcs_p])
          self._Ap.assemble()
          
          
if len(self._bcs_p) == 0:
              nullspace = _PETSc.NullSpace().create(constant=True, comm=self._mesh.comm)
              self._Ap.setNullSpace(nullspace)
          
          
if not self._low_memory:
              for i in range(self._mesh.geometry.dim):
                  # Preassemble pressure matrix for tentantive velocity
                  _fem.petsc.assemble_matrix(self._grad_p_M[i], self._grad_p[i])
                  self._grad_p_M[i].assemble()
          
          
        # Assemble pressure RHS matrices
                  _fem.petsc.assemble_matrix(self._p_M[i], self._p_rhs[i])

github.com

ComputationalPhysiology/oasisx/blob/main/src/oasisx/fracstep.py#L449


      
          def pressure_solve(self) -> np.int32:
              """
              Solve pressure correction problem
              """
          
          
    # Set difference vector to previous time step
              self._dp.x.array[:] = -self._ps.x.array[:]
          
          
    if len(self._bcs_p) == 0:
                  nullspace = self._Ap.getNullSpace()
                  nullspace.remove(self._b_c.vector)
              converged = self._solver_p.solve(self._b_c.vector, self._p)
          
          
    # Compute phi  = p^(n-1/2) - p*
              self._dp.vector.axpy(1, self._p.vector)
              return converged
          
          
def velocity_update(self, dt) -> npt.NDArray[np.int32]:
              """
              Compute Velocity update
              """

Topic		Replies	Views
Poisson problem with Neumann boundary condition variational formulation	12	727	September 2, 2024
How can I use Block Forms, Multiphenicsx, and Lagrange Multipliers to Solve Poisson's Equation with Pure Neumann Boundary Conditions? dolfinx	4	446	December 4, 2023
Inconsistent results with different dolfinx/petsc versions	0	432	July 28, 2022
Singular poisson equation in DOLFINx dolfinx	11	538	November 10, 2023
Complex valued Neumann b.c dolfinx	3	522	August 5, 2021

Issues Solving Pure Neumann Problems in dolfinx

Related Topics