Rho = Constant(mesh, PETSc.ScalarType(1.0)) # Densidad (cambiar) Mu = Constant(mesh, PETSc.ScalarType(1.0e-4)) # Viscosidad dinámica (cambiar) Inlet_marker, outlet_marker, wall_marker, obstacle_marker = 2, 3, 4, 5ĭt = Constant(mesh, PETSc.ScalarType(dt)) Mesh, _, ft = gmshio.read_from_msh('cyllinder.msh', MPI.COMM_WORLD, rank=0, gdim=gdim) Sorry for the comments in Spanish, they are not too meaningfullįrom dolfinx.fem import (Constant, Function, FunctionSpace,Īssemble_scalar, dirichletbc, form, locate_dofs_topological, set_bc)įrom import (apply_lifting, assemble_matrix, assemble_vector,įrom dolfinx.io import (VTXWriter, distribute_entity_data, gmshio)įrom ufl import (FacetNormal, FiniteElement, Identity, Measure, TestFunction, TrialFunction, VectorElement,Īs_vector, div, dot, ds, dx, inner, lhs, grad, nabla_grad, rhs, sym)
Code (the mesh generation is the same as in the tutorial with the exception of some modifications to L, H).
