atest-ufem-navier-stokes-periodic-channel-2d.py File Reference

Go to the source code of this file.

Namespaces
	atest-ufem-navier-stokes-periodic-channel-2d

Variables
float	h = 0.5
int	re_tau = 50
float	nu = 0.0001
int	a_tau = re_tau**2
tuple	Uc = a_tau/nu*(h**2/2.)
tuple	root = cf.Core.root()
tuple	env = cf.Core.environment()
tuple	model = root.create_component('NavierStokes', 'cf3.solver.ModelUnsteady')
tuple	domain = model.create_domain()
tuple	physics = model.create_physics('cf3.UFEM.NavierStokesPhysics')
tuple	solver = model.create_solver('cf3.UFEM.Solver')
tuple	ns_solver = solver.add_unsteady_solver('cf3.UFEM.NavierStokes')
tuple	blocks = domain.create_component('blocks', 'cf3.mesh.BlockMesh.BlockArrays')
tuple	points = blocks.create_points(dimensions = 2, nb_points = 6)
tuple	block_nodes = blocks.create_blocks(2)
tuple	block_subdivs = blocks.create_block_subdivisions()
float	grading = 1.2
tuple	gradings = blocks.create_block_gradings()
tuple	left_patch = blocks.create_patch_nb_faces(name = 'left', nb_faces = 2)
tuple	bottom_patch = blocks.create_patch_nb_faces(name = 'bottom', nb_faces = 1)
tuple	top_patch = blocks.create_patch_nb_faces(name = 'top', nb_faces = 1)
tuple	right_patch = blocks.create_patch_nb_faces(name = 'right', nb_faces = 2)
tuple	mesh = domain.create_component('Mesh', 'cf3.mesh.Mesh')
tuple	partitioner = domain.create_component('Partitioner', 'cf3.mesh.actions.PeriodicMeshPartitioner')
tuple	link_horizontal = partitioner.create_link_periodic_nodes()
list	u_in = [0., 0.]
tuple	ic_g = solver.InitialConditions.create_initial_condition(builder_name = 'cf3.UFEM.InitialConditionFunction', field_tag = 'body_force')
	bc = ns_solver.BoundaryConditions
tuple	pressure_integral = solver.add_unsteady_solver('cf3.UFEM.SurfaceIntegral')
tuple	bulk_velocity = solver.add_unsteady_solver('cf3.UFEM.BulkVelocity')
	lss = ns_solver.LSS
tuple	write_manager = solver.add_unsteady_solver('cf3.solver.actions.TimeSeriesWriter')
tuple	writer = write_manager.create_component('PVTUWriter', 'cf3.mesh.VTKXML.Writer')
tuple	time = model.create_time()
int	u_max = 0