examples.fvm.sodshock Namespace Reference

Variables
int	torder = 4
int	sorder = 4
int	jump = 2
int	gamma = 5/3
warpy	mesh
warpy	fluid
warpy	dens_ic
warpy	energy_ic
warpy	null_func = warpy.functions.set_to(name='null', value=0)
list	ic_applicators
list	ss1
list	ss_sets = [ss1]
warpy	writer
warpy	sim
	None
	gen_xdmf

Variable Documentation

dens_ic

warpy examples.fvm.sodshock.dens_ic

Initial value:

=  warpy.functions.heaviside(name='dens_ic',
                                    direction=[1,0,0],
                                    center=[0,0,0],
                                    neg=1,
                                    pos=jump)

energy_ic

warpy examples.fvm.sodshock.energy_ic

Initial value:

=  warpy.functions.heaviside(name='energy_ic',
                                               direction=[1,0,0],
                                               center=[0,0,0],
                                               neg=1/(gamma*(gamma-1)),
                                               pos=jump/(gamma*(gamma-1)))

fluid

warpy examples.fvm.sodshock.fluid

Initial value:

=  warpy.variable(name='fluid',
                       components=['rho', 'px', 'py', 'pz', 'e'])

gamma

int examples.fvm.sodshock.gamma = 5/3

gen_xdmf

examples.fvm.sodshock.gen_xdmf

ic_applicators

list examples.fvm.sodshock.ic_applicators

Initial value:

=  [
    warpy.applicator(spatial_order=sorder*2,
                            fun=dens_ic,
                            var=fluid,
                            components=['rho'],
                            spatial_scheme='FiniteVolume'),
    warpy.applicator(spatial_order=sorder,
                            fun=null_func,
                            var=fluid,
                            components=['px','py','pz'],
                            spatial_scheme='FiniteVolume'),
    warpy.applicator(spatial_order=sorder*2,
                            fun=energy_ic,
                            var=fluid,
                            components=['e'],
                            spatial_scheme='FiniteVolume')
    ]

jump

int examples.fvm.sodshock.jump = 2

mesh

warpy examples.fvm.sodshock.mesh

Initial value:

=  warpy.mesh.block(Bounds=[-1, 1],
                        NumCells=[100],
                        NodeSets=['Left', 'Right'],
                        NumLayers=sorder,
                        PeriodicBoundaries=['Left','Right'])

None

examples.fvm.sodshock.None

null_func

warpy examples.fvm.sodshock.null_func = warpy.functions.set_to(name='null', value=0)

sim

warpy examples.fvm.sodshock.sim

Initial value:

=  warpy.fvm_sim(name='sodshock',
                    torder=torder,
                    sorder=sorder,
                    meshes=[mesh],
                    initial_conditions=ic_applicators,
                    boundary_conditions=None,
                    writers=[writer],
                    subsolver_sets=ss_sets,
                    time=[0,1],
                    dt=1,
                    write_steps=100,
                    cfl_max=0.5,
                    verbosity='info')

sorder

int examples.fvm.sodshock.sorder = 4

ss1

list examples.fvm.sodshock.ss1

Initial value:

=  [warpy.hyperapps.euler_roe(name='euler_flux',
                                 variable=fluid,
                                 components=None,
                                 gamma=gamma)]

ss_sets

list examples.fvm.sodshock.ss_sets = [ss1]

torder

int examples.fvm.sodshock.torder = 4

writer

warpy examples.fvm.sodshock.writer

Initial value:

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