self_LinearEuler2D_t Module

  1. SELF_LinearEuler2D_t.f90
  2. self_LinearEuler2D_t

This module defines a class that can be used to solve the Linear Euler equations in 2-D. The Linear Euler Equations, here, are the Euler equations linearized about a motionless background state.

The solution variables are

The sound speed is carried as a solution variable so that it can vary in space. Its flux and source are identically zero, so the sound speed is held fixed in time at each spatial location.

The conservative flux is

and the source terms are null.


Uses

  • module~~self_lineareuler2d_t~~UsesGraph module~self_lineareuler2d_t self_LinearEuler2D_t module~self_model SELF_Model module~self_lineareuler2d_t->module~self_model module~self_dgmodel2d SELF_DGModel2D module~self_lineareuler2d_t->module~self_dgmodel2d module~self_mesh SELF_Mesh module~self_lineareuler2d_t->module~self_mesh module~self_boundaryconditions SELF_BoundaryConditions module~self_lineareuler2d_t->module~self_boundaryconditions module~self_supportroutines SELF_SupportRoutines module~self_model->module~self_supportroutines HDF5 HDF5 module~self_model->HDF5 FEQParse FEQParse module~self_model->FEQParse module~self_hdf5 SELF_HDF5 module~self_model->module~self_hdf5 module~self_metadata SELF_Metadata module~self_model->module~self_metadata module~self_dgmodel2d->module~self_boundaryconditions module~self_geometry_2d SELF_Geometry_2D module~self_dgmodel2d->module~self_geometry_2d module~self_gpu SELF_GPU module~self_dgmodel2d->module~self_gpu module~self_dgmodel2d_t SELF_DGModel2D_t module~self_dgmodel2d->module~self_dgmodel2d_t module~self_mesh_2d~2 SELF_Mesh_2D module~self_dgmodel2d->module~self_mesh_2d~2 module~self_gpuinterfaces SELF_GPUInterfaces module~self_dgmodel2d->module~self_gpuinterfaces module~self_constants SELF_Constants module~self_mesh->module~self_constants iso_c_binding iso_c_binding module~self_mesh->iso_c_binding module~self_domaindecomposition~2 SELF_DomainDecomposition module~self_mesh->module~self_domaindecomposition~2 module~self_boundaryconditions->module~self_supportroutines module~self_boundaryconditions->module~self_metadata module~self_boundaryconditions->iso_c_binding module~self_supportroutines->module~self_constants iso_fortran_env iso_fortran_env module~self_supportroutines->iso_fortran_env module~self_geometry_2d->module~self_supportroutines module~self_geometry_2d->module~self_constants module~self_geometry_2d->module~self_mesh_2d~2 module~self_tensor_2d~2 SELF_Tensor_2D module~self_geometry_2d->module~self_tensor_2d~2 module~self_vector_2d SELF_Vector_2D module~self_geometry_2d->module~self_vector_2d module~self_lagrange SELF_Lagrange module~self_geometry_2d->module~self_lagrange module~self_data SELF_Data module~self_geometry_2d->module~self_data module~self_scalar_2d~2 SELF_Scalar_2D module~self_geometry_2d->module~self_scalar_2d~2 module~self_gpu->iso_c_binding module~self_gpu_enums SELF_GPU_enums module~self_gpu->module~self_gpu_enums module~self_constants->iso_c_binding module~self_constants->iso_fortran_env module~self_hdf5->HDF5 module~self_hdf5->module~self_constants module~self_hdf5->iso_fortran_env mpi mpi module~self_hdf5->mpi module~self_metadata->HDF5 module~self_metadata->module~self_hdf5 module~self_dgmodel2d_t->module~self_model module~self_dgmodel2d_t->module~self_boundaryconditions module~self_dgmodel2d_t->module~self_supportroutines module~self_dgmodel2d_t->HDF5 module~self_dgmodel2d_t->FEQParse module~self_dgmodel2d_t->module~self_geometry_2d module~self_dgmodel2d_t->module~self_hdf5 module~self_dgmodel2d_t->module~self_metadata module~self_dgmodel2d_t->module~self_mesh_2d~2 module~self_mappedscalar_2d SELF_MappedScalar_2D module~self_dgmodel2d_t->module~self_mappedscalar_2d module~self_mappedvector_2d~2 SELF_MappedVector_2D module~self_dgmodel2d_t->module~self_mappedvector_2d~2 module~self_mesh_2d~2->module~self_gpu module~self_mesh_2d~2->iso_c_binding module~self_mesh_2d_t SELF_Mesh_2D_t module~self_mesh_2d~2->module~self_mesh_2d_t module~self_gpuinterfaces->module~self_gpu module~self_gpuinterfaces->iso_c_binding module~self_domaindecomposition~2->iso_c_binding module~self_domaindecomposition_t SELF_DomainDecomposition_t module~self_domaindecomposition~2->module~self_domaindecomposition_t module~self_domaindecomposition~2->mpi module~self_tensor_2d~2->module~self_gpu module~self_tensor_2d~2->module~self_constants module~self_tensor_2d~2->iso_c_binding module~self_tensor_2d~2->iso_fortran_env module~self_tensor_2d_t SELF_Tensor_2D_t module~self_tensor_2d~2->module~self_tensor_2d_t module~self_gpublas SELF_GPUBLAS module~self_tensor_2d~2->module~self_gpublas module~self_gpu_enums->iso_c_binding module~self_mappedscalar_2d->module~self_gpu module~self_mappedscalar_2d->module~self_gpuinterfaces module~self_mappedscalar_2d->iso_c_binding module~self_mappedscalar_2d_t SELF_MappedScalar_2D_t module~self_mappedscalar_2d->module~self_mappedscalar_2d_t module~self_mesh_2d_t->module~self_mesh module~self_mesh_2d_t->module~self_supportroutines module~self_mesh_2d_t->HDF5 module~self_mesh_2d_t->module~self_constants module~self_mesh_2d_t->module~self_hdf5 module~self_mesh_2d_t->iso_c_binding module~self_mesh_2d_t->module~self_domaindecomposition~2 module~self_mesh_2d_t->module~self_lagrange module~self_quadrature SELF_Quadrature module~self_mesh_2d_t->module~self_quadrature module~self_domaindecomposition_t->module~self_supportroutines module~self_domaindecomposition_t->module~self_constants module~self_domaindecomposition_t->iso_c_binding module~self_domaindecomposition_t->module~self_lagrange module~self_domaindecomposition_t->mpi module~self_vector_2d->module~self_gpu module~self_vector_2d->module~self_constants module~self_vector_2d->module~self_gpuinterfaces module~self_vector_2d->iso_c_binding module~self_vector_2d_t SELF_Vector_2D_t module~self_vector_2d->module~self_vector_2d_t module~self_vector_2d->module~self_gpublas module~self_lagrange->module~self_gpu module~self_lagrange->module~self_constants module~self_lagrange->iso_c_binding module~self_lagrange->iso_fortran_env module~self_lagrange_t SELF_Lagrange_t module~self_lagrange->module~self_lagrange_t module~self_data->HDF5 module~self_data->FEQParse module~self_data->module~self_constants module~self_data->module~self_hdf5 module~self_data->module~self_metadata module~self_data->iso_c_binding module~self_data->module~self_lagrange module~self_scalar_2d~2->module~self_gpu module~self_scalar_2d~2->module~self_constants module~self_scalar_2d~2->module~self_gpuinterfaces module~self_scalar_2d~2->iso_c_binding module~self_scalar_2d_t SELF_Scalar_2D_t module~self_scalar_2d~2->module~self_scalar_2d_t module~self_scalar_2d~2->module~self_gpublas module~self_mappedvector_2d~2->module~self_gpu module~self_mappedvector_2d~2->module~self_gpuinterfaces module~self_mappedvector_2d~2->iso_c_binding module~self_mappedvector_2d_t SELF_MappedVector_2D_t module~self_mappedvector_2d~2->module~self_mappedvector_2d_t module~self_mappedvector_2d~2->module~self_gpublas module~self_vector_2d_t->HDF5 module~self_vector_2d_t->FEQParse module~self_vector_2d_t->module~self_constants module~self_vector_2d_t->module~self_hdf5 module~self_vector_2d_t->module~self_metadata module~self_vector_2d_t->iso_c_binding module~self_vector_2d_t->module~self_lagrange module~self_vector_2d_t->module~self_data module~self_lagrange_t->module~self_supportroutines module~self_lagrange_t->HDF5 module~self_lagrange_t->module~self_constants module~self_lagrange_t->module~self_hdf5 module~self_lagrange_t->iso_c_binding module~self_lagrange_t->iso_fortran_env module~self_lagrange_t->module~self_quadrature module~self_scalar_2d_t->HDF5 module~self_scalar_2d_t->FEQParse module~self_scalar_2d_t->module~self_constants module~self_scalar_2d_t->module~self_hdf5 module~self_scalar_2d_t->module~self_metadata module~self_scalar_2d_t->iso_c_binding module~self_scalar_2d_t->module~self_lagrange module~self_scalar_2d_t->module~self_data module~self_mappedvector_2d_t->FEQParse module~self_mappedvector_2d_t->module~self_geometry_2d module~self_mappedvector_2d_t->module~self_constants module~self_mappedvector_2d_t->module~self_mesh_2d~2 module~self_mappedvector_2d_t->iso_c_binding module~self_mappedvector_2d_t->module~self_domaindecomposition~2 module~self_mappedvector_2d_t->module~self_tensor_2d~2 module~self_mappedvector_2d_t->module~self_vector_2d module~self_mappedvector_2d_t->module~self_lagrange module~self_tensor_2d_t->HDF5 module~self_tensor_2d_t->FEQParse module~self_tensor_2d_t->module~self_constants module~self_tensor_2d_t->module~self_hdf5 module~self_tensor_2d_t->module~self_metadata module~self_tensor_2d_t->iso_c_binding module~self_tensor_2d_t->module~self_lagrange module~self_tensor_2d_t->module~self_data module~self_quadrature->module~self_constants module~self_quadrature->iso_fortran_env module~self_gpublas->module~self_constants module~self_gpublas->iso_c_binding module~self_gpublas->module~self_gpu_enums module~self_gpublas->iso_fortran_env module~self_mappedscalar_2d_t->FEQParse module~self_mappedscalar_2d_t->module~self_geometry_2d module~self_mappedscalar_2d_t->module~self_constants module~self_mappedscalar_2d_t->module~self_mesh_2d~2 module~self_mappedscalar_2d_t->iso_c_binding module~self_mappedscalar_2d_t->module~self_domaindecomposition~2 module~self_mappedscalar_2d_t->module~self_tensor_2d~2 module~self_mappedscalar_2d_t->module~self_lagrange module~self_mappedscalar_2d_t->module~self_scalar_2d~2

Used by

  • module~~self_lineareuler2d_t~~UsedByGraph module~self_lineareuler2d_t self_LinearEuler2D_t module~self_lineareuler2d_pml_t self_LinearEuler2D_PML_t module~self_lineareuler2d_pml_t->module~self_lineareuler2d_t module~self_lineareuler2d self_LinearEuler2D module~self_lineareuler2d->module~self_lineareuler2d_t module~self_lineareuler2d~2 self_LinearEuler2D module~self_lineareuler2d~2->module~self_lineareuler2d_t module~self_lineareuler2d_pml self_LinearEuler2D_PML module~self_lineareuler2d_pml->module~self_lineareuler2d_pml_t module~self_lineareuler2d_pml~2 self_LinearEuler2D_PML module~self_lineareuler2d_pml~2->module~self_lineareuler2d_pml_t

Contents

Derived Types

Functions

Subroutines


Derived Types

type, public, extends(DGModel2D) :: LinearEuler2D_t

Components

TypeVisibilityAttributesNameInitial
type(MappedScalar2D), public :: dSdt
real(kind=prec), public :: dt
real(kind=prec), public :: entropy
type(MappedVector2D), public :: flux
type(MappedScalar2D), public :: fluxDivergence
real(kind=prec), public :: g =0.0_prec
type(SEMQuad), public, pointer:: geometry
logical, public :: gradient_enabled =.false.
type(BoundaryConditionList), public :: hyperbolicBCs
integer, public :: ioIterate =0
type(Mesh2D), public, pointer:: mesh
integer, public :: nvar
type(BoundaryConditionList), public :: parabolicBCs
logical, public :: prescribed_bcs_enabled =.true.
real(kind=prec), public :: rho0 =1.0_prec
type(MappedScalar2D), public :: solution
type(MappedVector2D), public :: solutionGradient
type(MappedScalar2D), public :: source
real(kind=prec), public :: t
logical, public :: tecplot_enabled =.true.
procedure(SELF_timeIntegrator), public, pointer:: timeIntegrator=> Euler_timeIntegrator
type(MappedScalar2D), public :: workSol

Type-Bound Procedures

procedure, public :: AdditionalFree => AdditionalFree_Model
procedure, public :: AdditionalInit => AdditionalInit_LinearEuler2D_t
procedure, public :: AdditionalOutput => AdditionalOutput_Model
procedure, public :: BoundaryFlux => BoundaryFlux_DGModel2D
procedure, public :: CalculateEntropy => CalculateEntropy_DGModel2D
procedure, public :: CalculateSolutionGradient => CalculateSolutionGradient_DGModel2D
procedure, public :: CalculateTendency => CalculateTendency_DGModel2D
procedure, public :: Euler_timeIntegrator
procedure, public :: FluxMethod => fluxmethod_DGModel2D
procedure, public :: ForwardStep => ForwardStep_Model
procedure, public :: Free => Free_DGModel2D
procedure, public :: GetSimulationTime
procedure, public :: IncrementIOCounter
procedure, public :: Init => Init_DGModel2D
procedure, public :: LowStorageRK2_timeIntegrator
procedure, public :: LowStorageRK3_timeIntegrator
procedure, public :: LowStorageRK4_timeIntegrator
procedure, public :: MapBoundaryConditions => MapBoundaryConditions_DGModel2D_t
procedure, public :: PostStepHook => PostStepHook_Model
procedure, public :: PreStepHook => PreStepHook_Model
procedure, public :: PreTendencyHook => PreTendencyHook_Model
procedure, public :: PrintType => PrintType_Model
procedure, public :: ReadModel => Read_DGModel2D_t
procedure, public :: ReportEntropy => ReportEntropy_Model
procedure, public :: ReportMetrics => ReportMetrics_DGModel2D_t
procedure, public :: ReportUserMetrics => ReportUserMetrics_Model
procedure, public :: SetBoundaryCondition => setboundarycondition_DGModel2D_t
procedure, public :: SetGradientBoundaryCondition => setgradientboundarycondition_DGModel2D_t
procedure, public :: SetMetadata => SetMetadata_LinearEuler2D_t
procedure, public :: SetNumberOfVariables => SetNumberOfVariables_LinearEuler2D_t
procedure, public :: SetSimulationTime
generic, public :: SetSolution => SetSolutionFromChar_DGModel2D_t, SetSolutionFromEqn_DGModel2D_t
generic, public :: SetTimeIntegrator => SetTimeIntegrator_withChar
procedure, public :: SourceMethod => sourcemethod_DGModel2D
procedure, public :: SphericalSoundWave => SphericalSoundWave_LinearEuler2D_t
procedure, public :: UpdateGRK2 => UpdateGRK2_DGModel2D
procedure, public :: UpdateGRK3 => UpdateGRK3_DGModel2D
procedure, public :: UpdateGRK4 => UpdateGRK4_DGModel2D
procedure, public :: UpdateSolution => UpdateSolution_DGModel2D
procedure, public :: WriteModel => Write_DGModel2D_t
procedure, public :: WriteTecplot => WriteTecplot_DGModel2D_t
procedure, public :: entropy_func => entropy_func_LinearEuler2D_t
procedure, public :: flux1D => flux1d_Model
procedure, public :: flux2d => flux2d_LinearEuler2D_t
procedure, public :: flux3D => flux3d_Model
procedure, public :: riemannflux1d => riemannflux1d_Model
procedure, public :: riemannflux2d => riemannflux2d_LinearEuler2D_t
procedure, public :: riemannflux3d => riemannflux3d_Model
procedure, public :: source1d => source1d_Model
procedure, public :: source2d => source2d_Model
procedure, public :: source3d => source3d_Model

Functions

public pure function entropy_func_LinearEuler2D_t(this, s) result(e)

The entropy function is the sum of kinetic and internal energy For the linear model, this is

Read more…

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D_t), intent(in) :: this
real(kind=prec), intent(in) :: s(1:this%nvar)

Return Value real(kind=prec)

public pure function flux2d_LinearEuler2D_t(this, s, dsdx) result(flux)

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D_t), intent(in) :: this
real(kind=prec), intent(in) :: s(1:this%nvar)
real(kind=prec), intent(in) :: dsdx(1:this%nvar,1:2)

Return Value real(kind=prec)(1:this%nvar,1:2)

public pure function riemannflux2d_LinearEuler2D_t(this, sL, sR, dsdx, nhat) result(flux)

Characteristic-decomposition (impedance-matched) interface flux for linear acoustics with possibly discontinuous sound speed.

Read more…

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D_t), intent(in) :: this
real(kind=prec), intent(in) :: sL(1:this%nvar)
real(kind=prec), intent(in) :: sR(1:this%nvar)
real(kind=prec), intent(in) :: dsdx(1:this%nvar,1:2)
real(kind=prec), intent(in) :: nhat(1:2)

Return Value real(kind=prec)(1:this%nvar)


Subroutines

public subroutine AdditionalInit_LinearEuler2D_t(this)

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D_t), intent(inout) :: this

public subroutine SetMetadata_LinearEuler2D_t(this)

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D_t), intent(inout) :: this

public subroutine SetNumberOfVariables_LinearEuler2D_t(this)

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D_t), intent(inout) :: this

public subroutine SphericalSoundWave_LinearEuler2D_t(this, rhoprime, Lr, x0, y0, c)

This subroutine sets the initial condition for a weak blast wave problem. The initial condition is given by

Read more…

Arguments

TypeIntentOptionalAttributesName
class(LinearEuler2D_t), intent(inout) :: this
real(kind=prec), intent(in) :: rhoprime
real(kind=prec), intent(in) :: Lr
real(kind=prec), intent(in) :: x0
real(kind=prec), intent(in) :: y0
real(kind=prec), intent(in) :: c

public subroutine hbc2d_NoNormalFlow_LinearEuler2D(bc, mymodel)

No-normal-flow boundary condition for 2D linear Euler equations. Reflects the velocity vector about the boundary normal while preserving density, pressure, and sound speed.

Arguments

TypeIntentOptionalAttributesName
class(BoundaryCondition), intent(in) :: bc
class(Model), intent(inout) :: mymodel