tem_miescatter_module Module

  1. tem_miescatter_module.f90
  2. tem_miescatter_module

This module gathers the definition of the analytical solution by scattering of an electromagnetic wave at a dielectric sphere/cylinder. The solution is given by means of Mie-series. The solution is given in: Cai, W., & Deng, S. (2003). An upwinding embedded boundary method for Maxwell's equations in media with material interfaces: 2D case. Journal of Computational Physics, 190(1), 159-183. doi:10.1016/S0021-9991(03)00269-9


Uses

  • module~~tem_miescatter_module~~UsesGraph module~tem_miescatter_module tem_miescatter_module module~tem_aux_module tem_aux_module module~tem_miescatter_module->module~tem_aux_module module~env_module env_module module~tem_miescatter_module->module~env_module module~tem_logging_module tem_logging_module module~tem_miescatter_module->module~tem_logging_module module~flu_binding flu_binding module~tem_miescatter_module->module~flu_binding module~aot_table_module aot_table_module module~tem_miescatter_module->module~aot_table_module module~aotus_module aotus_module module~tem_miescatter_module->module~aotus_module module~tem_param_module tem_param_module module~tem_miescatter_module->module~tem_param_module module~tem_aux_module->module~env_module module~tem_aux_module->module~tem_logging_module module~tem_aux_module->module~flu_binding module~tem_aux_module->module~aot_table_module module~tem_aux_module->module~aotus_module module~tem_revision_module tem_revision_module module~tem_aux_module->module~tem_revision_module module~tem_comm_env_module tem_comm_env_module module~tem_aux_module->module~tem_comm_env_module module~tem_tools_module tem_tools_module module~tem_aux_module->module~tem_tools_module module~tem_lua_requires_module tem_lua_requires_module module~tem_aux_module->module~tem_lua_requires_module mpi mpi module~tem_aux_module->mpi module~env_module->module~flu_binding module~env_module->module~aotus_module iso_fortran_env iso_fortran_env module~env_module->iso_fortran_env module~env_module->mpi module~tem_logging_module->module~env_module module~tem_logging_module->module~aot_table_module module~tem_logging_module->module~aotus_module module~tem_param_module->module~env_module module~tem_comm_env_module->mpi module~tem_tools_module->module~env_module module~tem_lua_requires_module->module~env_module module~tem_lua_requires_module->module~flu_binding module~tem_lua_requires_module->module~aot_table_module module~tem_lua_requires_module->module~aotus_module iso_c_binding iso_c_binding module~tem_lua_requires_module->iso_c_binding

Used by

  • module~~tem_miescatter_module~~UsedByGraph module~tem_miescatter_module tem_miescatter_module module~tem_spatial_module tem_spatial_module module~tem_spatial_module->module~tem_miescatter_module module~tem_spacetime_fun_module tem_spacetime_fun_module module~tem_spacetime_fun_module->module~tem_miescatter_module module~tem_spacetime_fun_module->module~tem_spatial_module module~tem_derived_module tem_derived_module module~tem_derived_module->module~tem_spacetime_fun_module module~tem_spacetime_var_module tem_spacetime_var_module module~tem_derived_module->module~tem_spacetime_var_module module~tem_variable_module tem_variable_module module~tem_derived_module->module~tem_variable_module program~tem_varsys_stfunvar_test tem_varSys_stfunVar_test program~tem_varsys_stfunvar_test->module~tem_spacetime_fun_module program~tem_varsys_stfunvar_test->module~tem_derived_module program~tem_varsys_stfunvar_test->module~tem_variable_module program~tem_varsys_derivevar_test tem_varSys_deriveVar_test program~tem_varsys_derivevar_test->module~tem_spacetime_fun_module module~tem_spacetime_var_module->module~tem_spatial_module module~tem_spacetime_var_module->module~tem_spacetime_fun_module module~tem_spacetime_var_module->module~tem_variable_module program~tem_variable_evaltype_test tem_variable_evaltype_test program~tem_variable_evaltype_test->module~tem_spacetime_fun_module program~tem_variable_evaltype_test->module~tem_derived_module program~tem_variable_evaltype_test->module~tem_variable_module program~tem_varsys_test tem_varSys_test program~tem_varsys_test->module~tem_spacetime_fun_module program~tem_varsys_test->module~tem_spacetime_var_module program~tem_varsys_test->module~tem_variable_module module~tem_varmap_module tem_varMap_module module~tem_varmap_module->module~tem_spacetime_fun_module module~tem_varmap_module->module~tem_spacetime_var_module module~tem_surfacedata_module tem_surfaceData_module module~tem_surfacedata_module->module~tem_spacetime_fun_module program~tem_spacetime_fun_test tem_spacetime_fun_test program~tem_spacetime_fun_test->module~tem_spacetime_fun_module program~tem_logical_opertor_test tem_logical_opertor_test program~tem_logical_opertor_test->module~tem_spacetime_fun_module program~tem_logical_opertor_test->module~tem_derived_module program~tem_logical_opertor_test->module~tem_variable_module program~tem_variable_extract_test tem_variable_extract_test program~tem_variable_extract_test->module~tem_spacetime_fun_module program~tem_variable_extract_test->module~tem_derived_module program~tem_variable_extract_test->module~tem_variable_module program~tem_variable_combine_test tem_variable_combine_Test program~tem_variable_combine_test->module~tem_spacetime_fun_module program~tem_variable_combine_test->module~tem_derived_module program~tem_variable_combine_test->module~tem_variable_module module~tem_face_module tem_face_module module~tem_face_module->module~tem_spacetime_fun_module module~tem_ini_condition_module tem_ini_condition_module module~tem_ini_condition_module->module~tem_spatial_module module~tem_variable_module->module~tem_spacetime_fun_module module~tem_bc_module tem_bc_module module~tem_bc_module->module~tem_spacetime_fun_module module~tem_bc_module->module~tem_varmap_module program~tem_varsys_opvar_test tem_varSys_opVar_test program~tem_varsys_opvar_test->module~tem_spacetime_fun_module program~tem_varsys_opvar_test->module~tem_derived_module program~tem_varsys_opvar_test->module~tem_variable_module program~tem_face_test~3 tem_face_test program~tem_face_test~3->module~tem_face_module module~tem_depend_module tem_depend_module module~tem_depend_module->module~tem_varmap_module module~tem_restart_module tem_restart_module module~tem_restart_module->module~tem_varmap_module module~hvs_output_module hvs_output_module module~hvs_output_module->module~tem_varmap_module program~tem_face_test tem_face_test program~tem_face_test->module~tem_face_module program~tem_face_test~2 tem_face_test program~tem_face_test~2->module~tem_face_module module~tem_tracking_module tem_tracking_module module~tem_tracking_module->module~tem_varmap_module module~tem_operation_var_module tem_operation_var_module module~tem_operation_var_module->module~tem_varmap_module module~tem_operation_var_module->module~tem_variable_module module~tem_convergence_module tem_convergence_module module~tem_convergence_module->module~tem_varmap_module

Contents

Interfaces

Derived Types

Functions

Subroutines


Interfaces

public interface tem_load_miescatter_displacementfieldz

  • private subroutine tem_load_miescatter(conf, thandle, me)

    load gauss pulse variables to set initial condition

    Arguments

    TypeIntentOptionalAttributesName
    type(flu_State) :: conf

    lua state type
    integer, intent(in) :: thandle

    aotus parent handle
    type(tem_miescatter_field_type), intent(out) :: me

    Global gauss pulse data type

public interface tem_load_miescatter_magneticfieldx

  • private subroutine tem_load_miescatter(conf, thandle, me)

    load gauss pulse variables to set initial condition

    Arguments

    TypeIntentOptionalAttributesName
    type(flu_State) :: conf

    lua state type
    integer, intent(in) :: thandle

    aotus parent handle
    type(tem_miescatter_field_type), intent(out) :: me

    Global gauss pulse data type

public interface tem_load_miescatter_magneticfieldy

  • private subroutine tem_load_miescatter(conf, thandle, me)

    load gauss pulse variables to set initial condition

    Arguments

    TypeIntentOptionalAttributesName
    type(flu_State) :: conf

    lua state type
    integer, intent(in) :: thandle

    aotus parent handle
    type(tem_miescatter_field_type), intent(out) :: me

    Global gauss pulse data type

Derived Types

type, public :: tem_miescatter_field_type

Parameters of the solution for Mie-Scatter at dielectric cylinder (infinite height in z direction).

Components

TypeVisibilityAttributesNameInitial
type(tem_miescatter_type), private :: miescatter

Parameter of the geometrical and material setup
type(tem_mieexpansion_type), private :: mieexpansion

The expansion parameters for the scattered solution

type, private :: tem_miescatter_type

Parameters of the solution for Mie-Scatter at dielectric cylinder (infinite height in z direction)

Components

TypeVisibilityAttributesNameInitial
real(kind=rk), private :: center(2) =[0.0_rk, 0.0_rk]

The center of the cylinder (in the x-y plane)
real(kind=rk), private :: radius =1.0_rk

The radius of the cylinder
real(kind=rk), private :: permeability_background =1.0_rk

Permeability (mu) of the background.
real(kind=rk), private :: permitivity_background =1.0_rk

Permitivity (epsilon) of the background.
real(kind=rk), private :: permeability_cylinder =1.0_rk

Permeability (mu) of the cylinder.
real(kind=rk), private :: permitivity_cylinder =1.0_rk

Permitivity (epsilon) of the cylinder.
real(kind=rk), private :: wavenumber_background

Wave number for the background
real(kind=rk), private :: wavenumber_cylinder

Wave number for the cylinder
real(kind=rk), private :: omega =2.0_rk*PI*1.0_rk

Angular frequency ( = 2 * pi * f , where f is the original frequency in Hertz)

type, private :: tem_mieexpansion_type

Expansion data for the Mier scatter solution.

Components

TypeVisibilityAttributesNameInitial
integer, private :: nCoeffs

The number of expansion coefficients
complex(kind=rk), private, allocatable:: c_tot(:)

Expansion coefficients inside the cylinder
complex(kind=rk), private, allocatable:: c_scat(:)

Expansion coefficients outside the cylinder (for the scattered field)

Functions

public function tem_eval_miescatter_displz(me, coord, time, n) result(res)

Evaluate displacement field (z component) for Mie-Scattering of electromagnetic wave at dielectric cylinder.

Arguments

TypeIntentOptionalAttributesName
type(tem_miescatter_field_type), intent(in) :: me

The function to evaluate
real(kind=rk), intent(in) :: coord(n,3)

barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates
real(kind=rk), intent(in) :: time

The time to evaluate the function at.
integer, intent(in) :: n

Number of points to evaluate the function for.

Return Value real(kind=rk)(n)

return value of the function

public function tem_eval_miescatter_magnx(me, coord, time, n) result(res)

Evaluate magnetic field (x-component) for Mie-Scattering of electromagnetic wave at dielectric cylinder.

Arguments

TypeIntentOptionalAttributesName
type(tem_miescatter_field_type), intent(in) :: me

The function to evaluate
real(kind=rk), intent(in) :: coord(n,3)

barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates
real(kind=rk), intent(in) :: time

The time to evaluate the function at.
integer, intent(in) :: n

Number of points to evaluate the function for.

Return Value real(kind=rk)(n)

return value of the function

public function tem_eval_miescatter_magny(me, coord, time, n) result(res)

Evaluate magnetic field (y-component) for Mie-Scattering of electromagnetic wave at dielectric cylinder.

Arguments

TypeIntentOptionalAttributesName
type(tem_miescatter_field_type), intent(in) :: me

The function to evaluate
real(kind=rk), intent(in) :: coord(n,3)

barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates
real(kind=rk), intent(in) :: time

The time to evaluate the function at.
integer, intent(in) :: n

Number of points to evaluate the function for.

Return Value real(kind=rk)(n)

return value of the function

private function tem_eval_miescatter_magnangular(me, coord, time, n) result(res)

Evaluate magnetizing field (angular-component) for Mie-Scattering of electromagnetic wave at dielectric cylinder.

Arguments

TypeIntentOptionalAttributesName
type(tem_miescatter_field_type), intent(in) :: me

The function to evaluate
real(kind=rk), intent(in) :: coord(n,3)

barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates
real(kind=rk), intent(in) :: time

The time to evaluate the function at.
integer, intent(in) :: n

Number of points to evaluate the function for.

Return Value real(kind=rk)(n)

return value of the function

private function tem_eval_miescatter_magnradial(me, coord, time, n) result(res)

Evaluate magnetizing field (radial-component) for Mie-Scattering of electromagnetic wave at dielectric cylinder.

Arguments

TypeIntentOptionalAttributesName
type(tem_miescatter_field_type), intent(in) :: me

The function to evaluate
real(kind=rk), intent(in) :: coord(n,3)

barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates
real(kind=rk), intent(in) :: time

The time to evaluate the function at.
integer, intent(in) :: n

Number of points to evaluate the function for.

Return Value real(kind=rk)(n)

return value of the function

private function cart2polar(x, y) result(polar)

Convert from cartesian coordinates (in the x-y plane) to polar coordinates (radius,angle)

Arguments

TypeIntentOptionalAttributesName
real(kind=rk) :: x

X coordinate
real(kind=rk) :: y

Y coordinate

Return Value real(kind=rk)(2)

Polar coordinates, radius (first entry) and angle (second entry)

private function bessel_jn_derivative(n, x) result(der)

Compute derivative for Bessel function of first kind of order n.

Arguments

TypeIntentOptionalAttributesName
integer :: n

The order of the function.
real(kind=rk), intent(in) :: x

The evaluation point.

Return Value real(kind=rk)

The function value.

private function hankel2_n_derivative(n, x) result(der)

Compute derivative for Hankel function of second kind of order n.

Arguments

TypeIntentOptionalAttributesName
integer :: n

The order of the function.
real(kind=rk), intent(in) :: x

The evaluation point.

Return Value complex(kind=rk)

The function value.

private function hankel2_n(n, x) result(val)

Compute Hankel function of second kind of order n.

Arguments

TypeIntentOptionalAttributesName
integer :: n

The order of the function.
real(kind=rk), intent(in) :: x

The evaluation point.

Return Value complex(kind=rk)

The function value.

private function tem_get_ncoeffs_miescat(conf, thandle) result(nCoeffs)

Read more…

Arguments

TypeIntentOptionalAttributesName
type(flu_State) :: conf

lua state type
integer, intent(in) :: thandle

aotus parent handle

Return Value integer

The number of coefficients necessary for the Fourier series


Subroutines

private subroutine tem_load_miescatter(conf, thandle, me)

load gauss pulse variables to set initial condition

Arguments

TypeIntentOptionalAttributesName
type(flu_State) :: conf

lua state type
integer, intent(in) :: thandle

aotus parent handle
type(tem_miescatter_field_type), intent(out) :: me

Global gauss pulse data type

private subroutine tem_init_data(conf, thandle, miescatter, mieexpansion)

Init the expansion coefficients for the Mie-Scattering.

Arguments

TypeIntentOptionalAttributesName
type(flu_State) :: conf

lua state type
integer, intent(in) :: thandle

aotus parent handle
type(tem_miescatter_type), intent(in) :: miescatter

The description of the Mie-Scattering setup.
type(tem_mieexpansion_type), intent(out) :: mieexpansion

The expansion coefficients to be filled