fenicsx.fluid_dynamics package

Submodules

fenicsx.fluid_dynamics.steady_ns module

class fenicsx.fluid_dynamics.steady_ns.fixedValueVelocity(u_fixed: ndarray, t: float, gdim: int)

Bases: object

fenicsx.fluid_dynamics.steady_ns.lid(x)

fenicsx.fluid_dynamics.steady_ns.lid_velocity_expression(x)

fenicsx.fluid_dynamics.steady_ns.noslip_boundary(x)

class fenicsx.fluid_dynamics.steady_ns.ns_steady_nl(domain: dolfinx.mesh.Mesh, ct: dolfinx.cpp.mesh.MeshTags_int32, ft: dolfinx.cpp.mesh.MeshTags_int32, bound_markers: dict, degree_u: int = 2)

Bases: object

F_fun(snes: petsc4py.PETSc.SNES, x: petsc4py.PETSc.Vec, F_vec: petsc4py.PETSc.Vec) → None

Assemble the residual.

J_fun(snes: petsc4py.PETSc.SNES, x: petsc4py.PETSc.Vec, J_mat: petsc4py.PETSc.Mat, P_mat: petsc4py.PETSc.Mat) → None

Assemble the jacobian.

assemble(maxIter=20, verbose=False)

create_snes_solution() → petsc4py.PETSc.Vec

Create a petsc4py.PETSc.Vec to be passed to petsc4py.PETSc.SNES.solve.

The returned vector will be initialized with the initial guess provided in self._solution.

obj_fun(snes: petsc4py.PETSc.SNES, x: petsc4py.PETSc.Vec) → float64

Compute the norm of the residual.

parameters(nu: float)

set_bc(boundary_type: dict, boundary_value: dict)

set_cavity_bc()

solve()

update_solution(x: petsc4py.PETSc.Vec) → None

Update self._solution with data in x.

fenicsx.fluid_dynamics.transient_ns module

class fenicsx.fluid_dynamics.transient_ns.drag_lift(domain: dolfinx.mesh.Mesh, ft: dolfinx.cpp.mesh.MeshTags_int32, physical_param: dict, obstacle_mark: int, points=[[0.15, 0.2, 0], [0.25, 0.2, 0]])

Bases: object

compute(t: float, dt: float, u_new: dolfinx.fem.Function, p_new: dolfinx.fem.Function)

class fenicsx.fluid_dynamics.transient_ns.pressure_projection(domain: dolfinx.mesh.Mesh, ct: dolfinx.cpp.mesh.MeshTags_int32, ft: dolfinx.cpp.mesh.MeshTags_int32, physical_param: dict, bound_markers: dict, time_adv='EI')

Bases: object

advance(u_tilde: dolfinx.fem.Function)

assembleForm(direct=False)

class fenicsx.fluid_dynamics.transient_ns.tentative_velocity(domain: dolfinx.mesh.Mesh, ct: dolfinx.cpp.mesh.MeshTags_int32, ft: dolfinx.cpp.mesh.MeshTags_int32, inlet, physical_param: dict, bound_markers: dict, time_adv='EI')

Bases: object

advance(t)

assembleForm(direct=False)

class fenicsx.fluid_dynamics.transient_ns.update_velocity(domain: dolfinx.mesh.Mesh, ct: dolfinx.cpp.mesh.MeshTags_int32, ft: dolfinx.cpp.mesh.MeshTags_int32, physical_param: dict, time_adv='EI')

Bases: object

advance(u_tilde: dolfinx.fem.Function, phi)

assembleForm(direct=False)

Module contents

models/fenicsx/fluid_dynamics.

OFELIA: FEniCSx models for Fluid Dynamics.