examples.dg.advection.multi_subdomain Namespace Reference

Variables
str	torder = 'RK1'
str	sorder = 'secondOrder'
list	velocity = [1., 0., 0.]
warpy	mesh
warpy	fluida
warpy	fluidb
warpy	a_ic
warpy	b_ic
list	applicators
list	apps_a = [warpy.apps.advection(name='advection', variable=fluida,velocity=velocity)]
list	apps_b = [warpy.apps.advection(name='advection', variable=fluidb,velocity=velocity)]
list	bc_apps = []
warpy	writer = warpy.host_actions.writer(name='writer',ReadVars=[fluida, fluidb])
warpy	ss_a
warpy	ss_b
warpy	bc_a_app
warpy	bc_b_app
list	variable_adjusters
warpy	temporal_solver
int	dt = 1e-3
warpy	dt_controller = warpy.dt_calc.fixed_dt(dt)
int	write_steps = 100
warpy	sim
	gen_xdmf
	True
	detect_nonscalar

Variable Documentation

a_ic

warpy examples.dg.advection.multi_subdomain.a_ic

Initial value:

=  warpy.functions.fourier(name='a_ic', sine_amplitudes=[1],
                               cosine_amplitudes=[0],
                               cosine_frequencies=[0],
                               sine_frequencies=[2*np.pi])

applicators

list examples.dg.advection.multi_subdomain.applicators

Initial value:

=  [warpy.applicator(spatial_order=sorder,fun=a_ic,var=fluida,spatial_scheme='Nodal'),
               warpy.applicator(spatial_order=sorder,fun=b_ic,var=fluidb,spatial_scheme='Nodal')]

apps_a

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

apps_b

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

b_ic

warpy examples.dg.advection.multi_subdomain.b_ic

Initial value:

=  warpy.functions.fourier(name='b_ic', sine_amplitudes=[0.5],
                               cosine_amplitudes=[0],
                               cosine_frequencies=[0],
                               sine_frequencies=[2*np.pi])

bc_a_app

warpy examples.dg.advection.multi_subdomain.bc_a_app

Initial value:

=  warpy.apps.multiplier_bc(name='mult_a',
                                    fluid_a=fluida, fluid_b=fluidb,
                                    multiplier=2,
                                    on_boundaries=['lc','cr'])

bc_apps

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

bc_b_app

warpy examples.dg.advection.multi_subdomain.bc_b_app

Initial value:

=  warpy.apps.multiplier_bc(name='mult_b',
                                    fluid_a=fluidb, fluid_b=fluida,
                                    multiplier=0.5,
                                    on_boundaries=['lc','cr'])

detect_nonscalar

examples.dg.advection.multi_subdomain.detect_nonscalar

dt

int examples.dg.advection.multi_subdomain.dt = 1e-3

dt_controller

warpy examples.dg.advection.multi_subdomain.dt_controller = warpy.dt_calc.fixed_dt(dt)

fluida

warpy examples.dg.advection.multi_subdomain.fluida

Initial value:

=  warpy.variable(name='a',
                        components=['a'],
                        basis_array_set=sorder,
                        subdomains=['left','right'],
                        virtual_subdomains=['mid'])

fluidb

warpy examples.dg.advection.multi_subdomain.fluidb

Initial value:

=  warpy.variable(name='b',
                        components=['b'],
                        basis_array_set=sorder,
                        subdomains=['mid'],
                        virtual_subdomains=['left','right'])

gen_xdmf

examples.dg.advection.multi_subdomain.gen_xdmf

mesh

warpy examples.dg.advection.multi_subdomain.mesh

Initial value:

=  warpy.mesh.arbitrary('mixedline_100.inp',
                            NumLayers=1,
                            PeriodicBoundaries=['leftwall','rightwall'],
                            basis_array_set=sorder)

sim

warpy examples.dg.advection.multi_subdomain.sim

Initial value:

=  warpy.dg_sim(name='multi_subdomain',
                   meshes=[mesh],
                   initial_conditions=applicators,
                   temporal_solvers=[temporal_solver],
                   writers=[writer],
                   time=[0,dt*write_steps],
                   dt_controller=dt_controller,
                   flexible_writeout=True,
                   write_steps=write_steps,
                   verbosity='info')

sorder

str examples.dg.advection.multi_subdomain.sorder = 'secondOrder'

ss_a

warpy examples.dg.advection.multi_subdomain.ss_a

Initial value:

=  warpy.spatial_solvers.dg(name="dg_a", spatial_order=sorder, applications=apps_a,
                                on_subdomains=['left', 'right'])

ss_b

warpy examples.dg.advection.multi_subdomain.ss_b

Initial value:

=  warpy.spatial_solvers.dg(name="dg_b", spatial_order=sorder, applications=apps_b,
                                on_subdomains=['mid'])

temporal_solver

warpy examples.dg.advection.multi_subdomain.temporal_solver

Initial value:

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

torder

str examples.dg.advection.multi_subdomain.torder = 'RK1'

True

examples.dg.advection.multi_subdomain.True

variable_adjusters

list examples.dg.advection.multi_subdomain.variable_adjusters

Initial value:

=  [
    warpy.variable_adjusters.boundary_conditions('bc_a', priority=1,
                                                 spatial_order=sorder,
                                                 applications=[bc_a_app],
                                                 copy_all=False,
                                                 is_virtual_bc=True,
                                                 on_subdomains=['mid']),
    warpy.variable_adjusters.boundary_conditions('bc_b', priority=1,
                                                 spatial_order=sorder,
                                                 applications=[bc_b_app],
                                                 copy_all=False,
                                                 is_virtual_bc=True,
                                                 on_subdomains=['left','right'])]

velocity

list examples.dg.advection.multi_subdomain.velocity = [1., 0., 0.]

write_steps

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

writer

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