atest-ufem-les-wale.py File Reference

Go to the source code of this file.

Namespaces
	atest-ufem-les-wale

Variables
int	h = 1
float	nu = 0.0001
int	re_tau = 180
	u_tau = re_tau*nu/h
int	a_tau = re_tau**2
tuple	Uc = a_tau/nu*(h**2/2.)
float	u_ref = 0.5
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	wale = solver.add_unsteady_solver('cf3.UFEM.les.WALE')
tuple	ns_solver = solver.add_unsteady_solver('cf3.UFEM.NavierStokesSemiImplicit')
float	tstep = 0.025
int	y_segs = 32
int	x_size = 4
int	z_size = 4
int	x_segs = 32
int	z_segs = 32
	ungraded_h = h
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()
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	nb_procs = cf.Core.nb_procs()
tuple	mesh = domain.create_component('Mesh', 'cf3.mesh.Mesh')
tuple	create_point_region = domain.create_component('CreatePointRegion', 'cf3.mesh.actions.AddPointRegion')
tuple	partitioner = domain.create_component('Partitioner', 'cf3.mesh.actions.PeriodicMeshPartitioner')
tuple	link_horizontal = partitioner.create_link_periodic_nodes()
tuple	link_spanwise = partitioner.create_link_periodic_nodes()
	lss = ns_solver.VelocityLSS.LSS
tuple	ic_u = solver.InitialConditions.NavierStokes.create_initial_condition(builder_name = 'cf3.UFEM.InitialConditionFunction', field_tag = 'navier_stokes_u_solution')
tuple	ic_g = solver.InitialConditions.NavierStokes.create_initial_condition(builder_name = 'cf3.UFEM.InitialConditionFunction', field_tag = 'body_force')
	bc_u = ns_solver.VelocityLSS.BC
tuple	stats = solver.add_unsteady_solver('cf3.solver.actions.TurbulenceStatistics')
tuple	restart_writer = solver.add_restart_writer()
tuple	dir_avg = solver.TimeLoop.children.WriteRestartManager.create_component('DirectionalAverage', 'cf3.solver.actions.DirectionalAverage')
tuple	skip_director = solver.add_unsteady_solver('cf3.solver.ActionDirectorWithSkip')
tuple	print_timings = skip_director.create_component('PrintTimingTree', 'cf3.common.PrintTimingTree')
tuple	time = model.create_time()