examples.dg.advection.advection Namespace Reference

Variables
str	torder = 'Dormand45'
str	sorder = 'thirdOrder'
list	velocity = [1.0, 0.0, 0.0]
bool	periodic_boundary_conditions = True
bool	use_limiters = False
warpy	mesh
warpy	fluid = warpy.variable(name='fluid', components=['rho'], basis_array_set=sorder)
warpy	fluid_ic
warpy	ic_va
warpy	ic_ha = warpy.host_actions.va_runner(name='ic_ha', variable_adjusters=[ic_va])
list	apps = [warpy.apps.advection(name='advection', variable=fluid, velocity=velocity)]
list	bc_apps = []
warpy	writer = warpy.host_actions.writer(name='writer', ReadVars=[fluid])
warpy	spatial_solver
list	variable_adjusters = []
warpy	temporal_solver
warpy	integrator
warpy	diag_writer
int	dt = 1e-4
float	dt_final = 2.0
warpy	dt_controller
	else :
int	write_steps = 100
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

apps

list examples.dg.advection.advection.apps = [warpy.apps.advection(name='advection', variable=fluid, velocity=velocity)]

bc_apps

list examples.dg.advection.advection.bc_apps = []

detect_nonscalar

examples.dg.advection.advection.detect_nonscalar

diag_writer

warpy examples.dg.advection.advection.diag_writer

Initial value:

=  warpy.host_actions.diagnostics_writer(
    name='diagnostics',
    write_interval=1e-3,
    diagnostics={integrator: ['integral_rho']},
    output_file='diagnostics.csv')

dt

int examples.dg.advection.advection.dt = 1e-4

dt_controller

warpy examples.dg.advection.advection.dt_controller

Initial value:

=  temporal_solver.create_integrating_controller(
        warpy.integrators.dg_max_integrator,
        warpy.apps.integrate.max_diff,
        init_dt=dt,
        max_dt=1e-2,
        max_growth=10,
        error_tol=1e-10,
        max_error=1e-1,
        stability_limit=False)

dt_final

float examples.dg.advection.advection.dt_final = 2.0

else

examples.dg.advection.advection.else :

fluid

warpy examples.dg.advection.advection.fluid = warpy.variable(name='fluid', components=['rho'], basis_array_set=sorder)

fluid_ic

warpy examples.dg.advection.advection.fluid_ic

Initial value:

=  warpy.apps.functions.arbitrary.fourier(
    name='fluid_ic', sine_amplitudes=[1], sine_frequencies=[2 * np.pi], variable=fluid)

gen_xdmf

examples.dg.advection.advection.gen_xdmf

ic_ha

warpy examples.dg.advection.advection.ic_ha = warpy.host_actions.va_runner(name='ic_ha', variable_adjusters=[ic_va])

ic_va

warpy examples.dg.advection.advection.ic_va

Initial value:

=  warpy.variable_adjusters.function_evaluation.function_evaluator(
    name='ic_va',
    priority=0,
    spatial_order=sorder,
    on_subdomains=['all'],
    applications=[fluid_ic])

integrator

warpy examples.dg.advection.advection.integrator

Initial value:

=  warpy.host_actions.integrator(
    name='mass_integrator',
    time_integrator=temporal_solver,
    patch_process_integrators=[
        warpy.integrators.dg_integrator(
            name='integrator',
            spatial_order=sorder,
            applications=[
                warpy.apps.integrate.int_q_dv(
                    name='int_q_dv',
                    variables=fluid,
                    stage=0,
                    variables_components=['rho'])
            ])
    ])

mesh

warpy examples.dg.advection.advection.mesh

Initial value:

=  warpy.mesh.block(
    Bounds=[-0.5, 0.5],
    NumCells=[128],
    NodeSets=['Left', 'Right'],
    NumLayers=1,
    PeriodicBoundaries=['Left', 'Right'] if periodic_boundary_conditions else None,
    basis_array_set=sorder)

periodic_boundary_conditions

bool examples.dg.advection.advection.periodic_boundary_conditions = True

sim

warpy examples.dg.advection.advection.sim

Initial value:

=  warpy.dg_sim(
    name='advection',
    meshes=[mesh],
    initial_conditions=[ic_ha],
    temporal_solvers=[temporal_solver],
    writers=[writer, diag_writer],
    time=[0, dt_final],
    dt_controller=dt_controller,
    flexible_writeout=False,
    write_steps=write_steps,
    verbosity='info')

sorder

str examples.dg.advection.advection.sorder = 'thirdOrder'

spatial_solver

warpy examples.dg.advection.advection.spatial_solver

Initial value:

=  warpy.spatial_solvers.dg(
    name="dg", spatial_order=sorder, applications=apps, on_subdomains=['all'])

temporal_solver

warpy examples.dg.advection.advection.temporal_solver

Initial value:

=  warpy.host_actions.erk(
    name='rk',
    scheme=torder,
    spatial_solvers=[spatial_solver],
    variable_adjusters=variable_adjusters)

torder

str examples.dg.advection.advection.torder = 'Dormand45'

True

examples.dg.advection.advection.True

use_limiters

bool examples.dg.advection.advection.use_limiters = False

variable_adjusters

list examples.dg.advection.advection.variable_adjusters = []

velocity

list examples.dg.advection.advection.velocity = [1.0, 0.0, 0.0]

write_steps

int examples.dg.advection.advection.write_steps = 100

writer

warpy examples.dg.advection.advection.writer = warpy.host_actions.writer(name='writer', ReadVars=[fluid])