Values = np.zeros((gdim, x.shape),dtype=PETSc.ScalarType) Scalar_element = FiniteElement("CG", mesh.ufl_cell(), 1) Vector_element = VectorElement("CG", mesh.ufl_cell(), 2) # ESPACIOS DE FUNCIONES Y FUNCIONES: ELEMENTOS TAYLOR-HOOD Nu = Constant(mesh, PETSc.ScalarType(mu/rho)) # Viscosidad cinemáticaį = Constant(mesh, PETSc.ScalarType((0, 0))) # Body forcesĬs = PETSc.ScalarType(0.035) # Constante de Smagorinsky 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)
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |