atl_stabilize_module

Routines for stabilization of a spectral scheme.

Uses

Help

Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.

Used by

Help

Graph Key

Nodes of different colours represent the following:

Solid arrows point from a submodule to the (sub)module which it is descended from. Dashed arrows point from a module or program unit to modules which it uses.

Derived Types

type, private :: atl_adaptive_orders_type

Components

Type	Visibility	Attributes		Name		Initial
real(kind=rk),	public,	allocatable	::	orders(:)

Functions

private function solve_admissible_state(isen_coef, my_eps, mean, limitedDens, momEnerg) result(t)

Solve for admissible state of the conservative positivity preserving limiter.

Arguments

Type	Intent		Name
real(kind=rk),	intent(in)	::	isen_coef
real(kind=rk),	intent(in)	::	my_eps
real(kind=rk),	intent(in)	::	mean(5)
real(kind=rk),	intent(in)	::	limitedDens
real(kind=rk),	intent(in)	::	momEnerg(2:5)

Return Value real(kind=rk)

private function solve_admissible_state_2d(isen_coef, my_eps, mean, limitedDens, momEnerg) result(t)

Solve for admissible state of the conservative positivity preserving limiter.

Arguments

Type	Intent		Name
real(kind=rk),	intent(in)	::	isen_coef
real(kind=rk),	intent(in)	::	my_eps
real(kind=rk),	intent(in)	::	mean(4)
real(kind=rk),	intent(in)	::	limitedDens
real(kind=rk),	intent(in)	::	momEnerg(2:4)

Return Value real(kind=rk)

Subroutines

public subroutine atl_stabilize(minlevel, maxlevel, statedata_list, statedata_stab_list, mesh_list, scheme_list, poly_proj_pos, poly_proj_list, equation, tree, bc, boundary, general, commStateTimer, material_list)

Subroutine to apply the stabilization procedure to the state vector.

Arguments

Type	Intent		Name
integer,	intent(in)	::	minlevel
integer,	intent(in)	::	maxlevel
type(atl_statedata_type),	intent(inout)	::	statedata_list(minlevel:maxlevel)
type(atl_statedata_type),	intent(inout)	::	statedata_stab_list(minlevel:maxlevel)
type(atl_cube_elem_type),	intent(inout)	::	mesh_list(minlevel:maxlevel)
type(atl_scheme_type),	intent(inout)	::	scheme_list(minlevel:maxlevel)
integer,	intent(inout)	::	poly_proj_pos(minlevel:maxlevel)
type(ply_poly_project_type),	intent(inout)	::	poly_proj_list(:)
type(atl_Equations_type),	intent(in)	::	equation
type(treelmesh_type),	intent(in)	::	tree	treelm mesh
type(atl_boundary_type),	intent(in)	::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)	::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)	::	general	General treelm settings.
integer,	intent(inout)	::	commStateTimer	Timer for measuring the communication time inside this routine.
type(atl_material_type),	intent(in)	::	material_list(minlevel:maxlevel)

public subroutine atl_cons_positivity_preserv(state, mesh, filter, euler, poly_proj)

Apply conservative limitation of denisty and energy to ensure positivity for density and pressure.

Arguments

Type	Intent		Name
type(atl_statedata_type),	intent(inout)	::	state
type(atl_cube_elem_type),	intent(in)	::	mesh
type(atl_cons_positivity_preserv_type),	intent(in)	::	filter
type(atl_euler_type),	intent(in)	::	euler
type(ply_poly_project_type),	intent(inout)	::	poly_proj

public subroutine atl_cons_positivity_preserv_2d(state, mesh, filter, euler, poly_proj)

Apply conservative limitation of denisty and energy to ensure positivity for density and pressure.

Arguments

Type	Intent		Name
type(atl_statedata_type),	intent(inout)	::	state
type(atl_cube_elem_type),	intent(in)	::	mesh
type(atl_cons_positivity_preserv_type),	intent(in)	::	filter
type(atl_euler_type),	intent(in)	::	euler
type(ply_poly_project_type),	intent(inout)	::	poly_proj

public subroutine atl_positivity_preserv(state, mesh, filter, euler, poly_proj)

Apply pointwise limitation of denisty and energy to ensure positivity for density and pressure.

Arguments

Type	Intent		Name
type(atl_statedata_type),	intent(inout)	::	state
type(atl_cube_elem_type),	intent(in)	::	mesh
type(atl_positivity_preserv_type),	intent(in)	::	filter
type(atl_euler_type),	intent(in)	::	euler
type(ply_poly_project_type),	intent(inout)	::	poly_proj

public subroutine atl_positivity_preserv_2d(state, mesh, filter, euler, poly_proj)

Apply pointwise limitation of denisty and energy to ensure positivity for density and pressure.

Arguments

Type	Intent		Name
type(atl_statedata_type),	intent(inout)	::	state
type(atl_cube_elem_type),	intent(in)	::	mesh
type(atl_positivity_preserv_type),	intent(in)	::	filter
type(atl_euler_type),	intent(in)	::	euler
type(ply_poly_project_type),	intent(inout)	::	poly_proj

private subroutine atl_covolume(minlevel, maxlevel, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Covolume filtering for 3D equation.

Arguments

Type	Intent	Optional		Name
integer,	intent(in)		::	minlevel	The minimal refinement level of the mesh.
integer,	intent(in)		::	maxlevel	The maximal refinement level of the mesh.
type(atl_statedata_type),	intent(inout)		::	state(minlevel:maxlevel)	State to be filtered (input and output)
type(atl_statedata_type),	intent(inout)		::	state_stab(minlevel:maxlevel,1:3)
type(atl_cube_elem_type),	intent(inout)		::	mesh(minlevel:maxlevel)	Mesh information for all the levels.
type(atl_covolume_type),	intent(in)		::	filter	The actual co-volume filter to be applied.
type(atl_scheme_type),	intent(inout)		::	scheme(minlevel:maxlevel)	List of numerical schemes for all the levels.
type(atl_Equations_type),	intent(in)		::	equation	Equation kind information
type(treelmesh_type),	intent(in)		::	tree	treelm mesh
type(ply_poly_project_type),	intent(inout)		::	poly_proj(:)	The list of projections.
integer,	intent(inout)		::	poly_proj_pos(minlevel:maxlevel)	The mapping from each level to the projections.
type(atl_boundary_type),	intent(in)		::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)		::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)		::	general	General treelm settings.
integer,	intent(inout)		::	commStateTimer	Timer for measuring the communication time inside this routine.
type(atl_adaptive_orders_type),	intent(in),	optional	::	adaptive_orders(minlevel:maxlevel)	The filters orders, if adaptive filter is applied

private subroutine atl_covolume_tocovolume(minlevel, maxlevel, iDir, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Recursive routine to project the state from primal mesh to covolume mesh.

Arguments

Type	Intent	Optional		Name
integer,	intent(in)		::	minlevel	The minimal refinement level of the mesh.
integer,	intent(in)		::	maxlevel	The maximal refinement level of the mesh.
integer,	intent(in)		::	iDir	The spatial direction for the projection: 1 -> x direction \n 2 -> y direction \n 3 -> z direction \n
type(atl_statedata_type),	intent(inout)		::	state(minlevel:maxlevel)	State to be filtered
type(atl_statedata_type),	intent(inout)		::	state_stab(minlevel:maxlevel,1:3)
type(atl_cube_elem_type),	intent(inout)		::	mesh(minlevel:maxlevel)	Mesh information for all the levels.
type(atl_covolume_type),	intent(in)		::	filter	The actual co-volume filter to be applied.
type(atl_scheme_type),	intent(inout)		::	scheme(minlevel:maxlevel)	List of numerical schemes for all the levels.
type(atl_Equations_type),	intent(in)		::	equation	Equation kind information
type(treelmesh_type),	intent(in)		::	tree	treelm mesh
type(ply_poly_project_type),	intent(inout)		::	poly_proj(:)	The list of projections.
integer,	intent(inout)		::	poly_proj_pos(minlevel:maxlevel)	The mapping from each level to the projections.
type(atl_boundary_type),	intent(in)		::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)		::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)		::	general	General treelm settings.
integer,	intent(inout)		::	commStateTimer	Timer for measuring the communication time inside this routine.
type(atl_adaptive_orders_type),	intent(in),	optional	::	adaptive_orders(minlevel:maxlevel)	The filters orders, if adaptive filter is applied

private recursive subroutine atl_interpolate_elemstate(minlevel, maxlevel, iLevel, iDir, state, state_stab, mesh, scheme, equation, tree, dimen, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer)

Recursive interpolation of element states among the levels of the mesh.

Arguments

Type	Intent		Name
integer,	intent(in)	::	minlevel	The minimal refinement level of the mesh.
integer,	intent(in)	::	maxlevel	The maximal refinement level of the mesh.
integer,	intent(in)	::	iLevel	The current level.
integer,	intent(in)	::	iDir	The direction to interpolate
type(atl_statedata_type),	intent(inout)	::	state(minlevel:maxlevel)	State to be filtered (input and output)
type(atl_statedata_type),	intent(inout)	::	state_stab(minlevel:maxlevel,1:3)
type(atl_cube_elem_type),	intent(inout)	::	mesh(minlevel:maxlevel)	Mesh information for all the levels.
type(atl_scheme_type),	intent(inout)	::	scheme(minlevel:maxlevel)	List of numerical schemes for all the levels.
type(atl_Equations_type),	intent(in)	::	equation	Equation kind information
type(treelmesh_type),	intent(in)	::	tree	treelm mesh
integer,	intent(in)	::	dimen	The spatial dimension under consideration.
type(ply_poly_project_type),	intent(inout)	::	poly_proj(:)	The list of projections.
integer,	intent(inout)	::	poly_proj_pos(minlevel:maxlevel)	The mapping from each level to the projections.
type(atl_boundary_type),	intent(in)	::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)	::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)	::	general	General treelm settings.
integer,	intent(inout)	::	commStateTimer	Timer for measuring the communication time inside this routine.

private subroutine atl_covolume_2d(minlevel, maxlevel, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Covolume filtering for 2D equation.

Arguments

Type	Intent	Optional		Name
integer,	intent(in)		::	minlevel	The minimal refinement level of the mesh.
integer,	intent(in)		::	maxlevel	The maximal refinement level of the mesh.
type(atl_statedata_type),	intent(inout)		::	state(minlevel:maxlevel)	State to be filtered (input and output)
type(atl_statedata_type),	intent(inout)		::	state_stab(minlevel:maxlevel,1:3)
type(atl_cube_elem_type),	intent(inout)		::	mesh(minlevel:maxlevel)	Mesh information for all the levels.
type(atl_covolume_type),	intent(in)		::	filter	The actual co-volume filter to be applied.
type(atl_scheme_type),	intent(inout)		::	scheme(minlevel:maxlevel)	List of numerical schemes for all the levels.
type(atl_Equations_type),	intent(in)		::	equation	Equation kind information
type(treelmesh_type),	intent(in)		::	tree	treelm mesh
type(ply_poly_project_type),	intent(inout)		::	poly_proj(:)	The list of projections.
integer,	intent(inout)		::	poly_proj_pos(minlevel:maxlevel)	The mapping from each level to the projections.
type(atl_boundary_type),	intent(in)		::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)		::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)		::	general	General treelm settings.
integer,	intent(inout)		::	commStateTimer	Timer for measuring the communication time inside this routine.
type(atl_adaptive_orders_type),	intent(in),	optional	::	adaptive_orders(minlevel:maxlevel)	The filters orders, if adaptive filter is applied

private subroutine atl_covolume_tocovolume_2d(minlevel, maxlevel, iDir, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Recursive routine to project the state from primal mesh to covolume mesh.

Arguments

Type	Intent	Optional		Name
integer,	intent(in)		::	minlevel	The minimal refinement level of the mesh.
integer,	intent(in)		::	maxlevel	The maximal refinement level of the mesh.
integer,	intent(in)		::	iDir	The spatial direction for the projection: 1 -> x direction \n 2 -> y direction \n 3 -> z direction \n
type(atl_statedata_type),	intent(inout)		::	state(minlevel:maxlevel)	State to be filtered
type(atl_statedata_type),	intent(inout)		::	state_stab(minlevel:maxlevel,1:3)
type(atl_cube_elem_type),	intent(inout)		::	mesh(minlevel:maxlevel)	Mesh information for all the levels.
type(atl_covolume_type),	intent(in)		::	filter	The actual co-volume filter to be applied.
type(atl_scheme_type),	intent(inout)		::	scheme(minlevel:maxlevel)	List of numerical schemes for all the levels.
type(atl_Equations_type),	intent(in)		::	equation	Equation kind information
type(treelmesh_type),	intent(in)		::	tree	treelm mesh
type(ply_poly_project_type),	intent(inout)		::	poly_proj(:)	The list of projections.
integer,	intent(inout)		::	poly_proj_pos(minlevel:maxlevel)	The mapping from each level to the projections.
type(atl_boundary_type),	intent(in)		::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)		::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)		::	general	General treelm settings.
integer,	intent(inout)		::	commStateTimer	Timer for measuring the communication time inside this routine.
type(atl_adaptive_orders_type),	intent(in),	optional	::	adaptive_orders(minlevel:maxlevel)	The filters orders, if adaptive filter is applied

private subroutine atl_covolume_1d(minlevel, maxlevel, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer)

Covolume filtering for 1D equation.

Arguments

Type	Intent		Name
integer,	intent(in)	::	minlevel	The minimal refinement level of the mesh.
integer,	intent(in)	::	maxlevel	The maximal refinement level of the mesh.
type(atl_statedata_type),	intent(inout)	::	state(minlevel:maxlevel)	State to be filtered (input and output)
type(atl_statedata_type),	intent(inout)	::	state_stab(minlevel:maxlevel,1:3)
type(atl_cube_elem_type),	intent(inout)	::	mesh(minlevel:maxlevel)	Mesh information for all the levels.
type(atl_covolume_type),	intent(in)	::	filter	The actual co-volume filter to be applied.
type(atl_scheme_type),	intent(inout)	::	scheme(minlevel:maxlevel)	List of numerical schemes for all the levels.
type(atl_Equations_type),	intent(in)	::	equation	Equation kind information
type(treelmesh_type),	intent(in)	::	tree	treelm mesh
type(ply_poly_project_type),	intent(inout)	::	poly_proj(:)	The list of projections.
integer,	intent(inout)	::	poly_proj_pos(minlevel:maxlevel)	The mapping from each level to the projections.
type(atl_boundary_type),	intent(in)	::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)	::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)	::	general	General treelm settings.
integer,	intent(inout)	::	commStateTimer	Timer for measuring the communication time inside this routine.

private subroutine atl_covolume_tocovolume_1d(minlevel, maxlevel, iDir, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer)

Recursive routine to project the state from primal mesh to covolume mesh.

Arguments

Type	Intent		Name
integer,	intent(in)	::	minlevel	The minimal refinement level of the mesh.
integer,	intent(in)	::	maxlevel	The maximal refinement level of the mesh.
integer,	intent(in)	::	iDir	The spatial direction for the projection: 1 -> x direction \n 2 -> y direction \n 3 -> z direction \n
type(atl_statedata_type),	intent(inout)	::	state(minlevel:maxlevel)	State to be filtered
type(atl_statedata_type),	intent(inout)	::	state_stab(minlevel:maxlevel,1:3)
type(atl_cube_elem_type),	intent(inout)	::	mesh(minlevel:maxlevel)	Mesh information for all the levels.
type(atl_covolume_type),	intent(in)	::	filter	The actual co-volume filter to be applied.
type(atl_scheme_type),	intent(inout)	::	scheme(minlevel:maxlevel)	List of numerical schemes for all the levels.
type(atl_Equations_type),	intent(in)	::	equation	Equation kind information
type(treelmesh_type),	intent(in)	::	tree	treelm mesh
type(ply_poly_project_type),	intent(inout)	::	poly_proj(:)	The list of projections.
integer,	intent(inout)	::	poly_proj_pos(minlevel:maxlevel)	The mapping from each level to the projections.
type(atl_boundary_type),	intent(in)	::	bc(:)	The global description of the boundaries.
type(atl_level_boundary_type),	intent(in)	::	boundary(minlevel:maxlevel)	The levelwise collection of boundary elements and boundary faces.
type(tem_general_type),	intent(inout)	::	general	General treelm settings.
integer,	intent(inout)	::	commStateTimer	Timer for measuring the communication time inside this routine.

private subroutine atl_spectral_visc_1d(state, mesh, filter, maxPolyDeg)

Damp )the modal coefficients of the state vector by a given spectral viscosity method.

Arguments

Type	Intent		Name
type(atl_statedata_type),	intent(inout)	::	state
type(atl_cube_elem_type),	intent(in)	::	mesh
type(atl_spectral_visc_type),	intent(in)	::	filter
integer,	intent(in)	::	maxPolyDeg

private subroutine atl_spectral_visc_2d(state, mesh, filter, maxPolyDeg, orders)

Damp )the modal coefficients of the state vector by a given spectral viscosity method.

Arguments

Type	Intent	Optional		Name
type(atl_statedata_type),	intent(inout)		::	state
type(atl_cube_elem_type),	intent(in)		::	mesh
type(atl_spectral_visc_type),	intent(in)		::	filter
integer,	intent(in)		::	maxPolyDeg
real(kind=rk),	intent(in),	optional	::	orders(:)

private subroutine atl_spectral_visc_3d(state, mesh, filter, maxPolyDeg, orders)

Damp )the modal coefficients of the state vector by a given spectral viscosity method.

Arguments

Type	Intent	Optional		Name
type(atl_statedata_type),	intent(inout)		::	state
type(atl_cube_elem_type),	intent(in)		::	mesh
type(atl_spectral_visc_type),	intent(in)		::	filter
integer,	intent(in)		::	maxPolyDeg
real(kind=rk),	intent(in),	optional	::	orders(:)

private subroutine atl_cheb_spectral_visc_1d(state, mesh, filter, poly_proj, maxPolyDeg)

Damp the modal coefficients of the state vector by a given spectral viscosity method.

Arguments

Type	Intent		Name
type(atl_statedata_type),	intent(inout)	::	state
type(atl_cube_elem_type),	intent(in)	::	mesh
type(atl_spectral_visc_type),	intent(in)	::	filter
type(ply_poly_project_type),	intent(inout)	::	poly_proj
integer,	intent(in)	::	maxPolyDeg

atl_stabilize_module Module

Contents

Derived Types

Functions

Subroutines

Uses

Used by

Derived Types

type, private :: atl_adaptive_orders_type

Components

Functions

private function solve_admissible_state(isen_coef, my_eps, mean, limitedDens, momEnerg) result(t)

Arguments

Return Value real(kind=rk)

private function solve_admissible_state_2d(isen_coef, my_eps, mean, limitedDens, momEnerg) result(t)

Arguments

Return Value real(kind=rk)

Subroutines

public subroutine atl_stabilize(minlevel, maxlevel, statedata_list, statedata_stab_list, mesh_list, scheme_list, poly_proj_pos, poly_proj_list, equation, tree, bc, boundary, general, commStateTimer, material_list)

Arguments

public subroutine atl_cons_positivity_preserv(state, mesh, filter, euler, poly_proj)

Arguments

public subroutine atl_cons_positivity_preserv_2d(state, mesh, filter, euler, poly_proj)

Arguments

public subroutine atl_positivity_preserv(state, mesh, filter, euler, poly_proj)

Arguments

public subroutine atl_positivity_preserv_2d(state, mesh, filter, euler, poly_proj)

Arguments

private subroutine atl_covolume(minlevel, maxlevel, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Arguments

private subroutine atl_covolume_tocovolume(minlevel, maxlevel, iDir, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Arguments

private recursive subroutine atl_interpolate_elemstate(minlevel, maxlevel, iLevel, iDir, state, state_stab, mesh, scheme, equation, tree, dimen, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer)

Arguments

private subroutine atl_covolume_2d(minlevel, maxlevel, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Arguments

private subroutine atl_covolume_tocovolume_2d(minlevel, maxlevel, iDir, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer, adaptive_orders)

Arguments

private subroutine atl_covolume_1d(minlevel, maxlevel, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer)

Arguments

private subroutine atl_covolume_tocovolume_1d(minlevel, maxlevel, iDir, state, state_stab, mesh, filter, scheme, equation, tree, poly_proj, poly_proj_pos, bc, boundary, general, commStateTimer)

Arguments

private subroutine atl_spectral_visc_1d(state, mesh, filter, maxPolyDeg)

Arguments

private subroutine atl_spectral_visc_2d(state, mesh, filter, maxPolyDeg, orders)

Arguments

private subroutine atl_spectral_visc_3d(state, mesh, filter, maxPolyDeg, orders)

Arguments

private subroutine atl_cheb_spectral_visc_1d(state, mesh, filter, poly_proj, maxPolyDeg)

Arguments