STConvectiveMatrix example 31

!!! note "" This example shows how to USE the SUBROUTINE called STConvectiveMatrix to create a space-time convective matrix. Triangle3 in space and Line2 in time.

Here, we want to DO the following.

M(I,J,a,b)=\int_{I_{n}}\int_{\Omega}c\frac{\partial N^{I}T_{a}}{\partial x}\frac{\partial N^{J}T_{b}}{\partial t}d\Omega dt

M(I,J,a,b)=\int_{I_{n}}\int_{\Omega}c\frac{\partial N^{I}T_{a}}{\partial y}\frac{\partial N^{J}T_{b}}{\partial t}d\Omega dt

M(I,J,a,b)=\int_{I_{n}}\int_{\Omega}c\frac{\partial N^{I}T_{a}}{\partial z}\frac{\partial N^{J}T_{b}}{\partial t}d\Omega dt

M(I,J,a,b)=\int_{I_{n}}\int_{\Omega}c\frac{\partial N^{I}T_{a}}{\partial t}\frac{\partial N^{J}T_{b}}{\partial x}d\Omega dt

M(I,J,a,b)=\int_{I_{n}}\int_{\Omega}c\frac{\partial N^{I}T_{a}}{\partial t}\frac{\partial N^{J}T_{b}}{\partial y}d\Omega dt

M(I,J,a,b)=\int_{I_{n}}\int_{\Omega}c\frac{\partial N^{I}T_{a}}{\partial t}\frac{\partial N^{J}T_{b}}{\partial z}d\Omega dt

!!! warning "" $c$ is scalar [[FEVariable_]], which can be a constant, or a FUNCTION of space-time, or some nonlinear FUNCTION.

In this example, convective matrix is formed for

[[ReferenceTriangle_]] Triangle3 element for space
[[ReferenceLine_]] Line2 element for time
[[QuadraturePoint_]] GaussLegendre
Space-time nodal values $c$

Modules and classes

[[ElemshapeData_]]
[[STElemshapeData_]]
[[QuadraturePoint_]]
[[ReferenceLine_]]
[[ReferenceTriangle_]]

Usage

PROGRAM main
    USE easifemBase
    IMPLICIT NONE
    TYPE(STElemshapeData_), ALLOCATABLE :: test(:)
    TYPE(ElemshapeData_) :: time_elemsd
    TYPE(Quadraturepoint_) :: quadFortime
    TYPE(Quadraturepoint_) :: quadForspace
    TYPE(ReferenceTriangle_):: refelemForSpace
    TYPE(ReferenceLine_) :: refelemForTime
    INTEGER(I4B) :: ii
    INTEGER(I4B), PARAMETER :: nsd=2, nns=3, nnt=2
    INTEGER(I4B), PARAMETER :: orderForTime=2, orderForSpace=1
    REAL(DFP), PARAMETER :: tij(1, 2) = RESHAPE([-1,1], [1,2])
    REAL(DFP), PARAMETER :: xij(2, 3)=RESHAPE([0,0,1,0,0,1], [nsd, nns])
    ! spatial nodal coordinates
    REAL(DFP), ALLOCATABLE :: xija(:, :, :), mat(:,:)
    ! spatial-temporal nodal coordinates
    REAL(DFP), parameter :: c(3,2)=reshape([1.0,1.0,1.0,1.0,1.0,1.0],[3,2])
    type(FEVariable_) :: cvar

!!! note "" First, we initiate a [[ReferenceLine_]] element for time domain. Note that nsd should be 1 when making reference element for time domain. Generate the quadrature points, and initiates an instance of [[ElemshapeData_]].

    refelemForTime= ReferenceLine(nsd=1)
    CALL Initiate(obj=quadFortime, &
		& refelem=refelemForTime,&
		& order=orderForTime, &
      	& quadratureType="GaussLegendre" )
    CALL Initiate( &
    	& obj=time_elemsd, &
	  	& quad=quadForTime, &
		& refelem=refelemForTime, &
      	& ContinuityType=typeH1,&
		& InterpolType=TypeLagrangeInterpolation)
    CALL Set(obj=time_elemsd, &
		& val=tiJ, N=time_elemsd%N, &
        & dNdXi=time_elemsd%dNdXi)

!!! note "" Initiate [[STElemshapeData_]].

    CALL Initiate(obj=test, elemsd=time_elemsd)

!!! note "" Generating shape functions for space-elements. Here, we are selecting a triangular element

    refelemForSpace = ReferenceTriangle(nsd=nsd)
    CALL Initiate(obj=quadForSpace, &
		& refelem=refelemForSpace, &
		& order=orderForSpace, &
		& quadratureType='GaussLegendre')

    DO ii = 1, SIZE(test)
      CALL Initiate( obj=test(ii), &
	  	& quad=quadForSpace, &
		& refelem=refelemForSpace, &
        & ContinuityType=typeH1, &
		& InterpolType=TypeLagrangeInterpolation)
    END DO

!!! note "" Setting the remaining DATA in obj. Here, xija are the space-time nodal coordinates.

	CALL Reallocate(xija, nsd, nns, nnt)
    DO ii = 1, nnt; xija(:, :, ii) = xij; END DO
    DO ii = 1, SIZE(test)
        CALL Set(obj=test(ii), &
            & val=xija, &
			& N=test(ii)%N, &
            & dNdXi=test(ii)%dNdXi, &
            & T=test(ii)%T)
    END DO

!!! note "" Let us now create the space-time convective matrix.

    cvar = NodalVariable(c, typeFEVariableScalar, typeFEVariableSpaceTime)

    mat=ConvectiveMatrix(test=test, trial=test, &
        & term1=del_t, term2=del_x, &
        & c=cvar)
    CALL Display(mat, "mat:")

??? example "Results"

                                mat:                                   
--------------------------------------------------------------------------
8.33333E-2  -8.33333E-2  0.00000E+0   8.33333E-2  -8.33333E-2  0.00000E+0
8.33333E-2  -8.33333E-2  0.00000E+0   8.33333E-2  -8.33333E-2  0.00000E+0
8.33333E-2  -8.33333E-2  0.00000E+0   8.33333E-2  -8.33333E-2  0.00000E+0
-8.33333E-2   8.33333E-2  0.00000E+0  -8.33333E-2   8.33333E-2  0.00000E+0
-8.33333E-2   8.33333E-2  0.00000E+0  -8.33333E-2   8.33333E-2  0.00000E+0
-8.33333E-2   8.33333E-2  0.00000E+0  -8.33333E-2   8.33333E-2  0.00000E+0

    mat=ConvectiveMatrix(test=test, trial=test, &
        & term1=del_t, term2=del_y, &
        & c=cvar)
    CALL Display(mat, "mat:")

??? example "Results"

                                   mat:                                   
--------------------------------------------------------------------------
 8.33333E-2  0.00000E+0  -8.33333E-2   8.33333E-2  0.00000E+0  -8.33333E-2
 8.33333E-2  0.00000E+0  -8.33333E-2   8.33333E-2  0.00000E+0  -8.33333E-2
 8.33333E-2  0.00000E+0  -8.33333E-2   8.33333E-2  0.00000E+0  -8.33333E-2
-8.33333E-2  0.00000E+0   8.33333E-2  -8.33333E-2  0.00000E+0   8.33333E-2
-8.33333E-2  0.00000E+0   8.33333E-2  -8.33333E-2  0.00000E+0   8.33333E-2
-8.33333E-2  0.00000E+0   8.33333E-2  -8.33333E-2  0.00000E+0   8.33333E-2

    mat=ConvectiveMatrix(test=test, trial=test, &
        & term1=del_x, term2=del_t, &
        & c=cvar)
    CALL Display(mat, "mat:")

??? example "Results"

                                    mat:                                    
----------------------------------------------------------------------------
8.33333E-2   8.33333E-2   8.33333E-2  -8.33333E-2  -8.33333E-2  -8.33333E-2
-8.33333E-2  -8.33333E-2  -8.33333E-2   8.33333E-2   8.33333E-2   8.33333E-2
0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0
8.33333E-2   8.33333E-2   8.33333E-2  -8.33333E-2  -8.33333E-2  -8.33333E-2
-8.33333E-2  -8.33333E-2  -8.33333E-2   8.33333E-2   8.33333E-2   8.33333E-2
0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0

    mat=ConvectiveMatrix(test=test, trial=test, &
        & term1=del_y, term2=del_t, &
        & c=cvar)
    CALL Display(mat, "mat:")

                                    mat:                                    
----------------------------------------------------------------------------
8.33333E-2   8.33333E-2   8.33333E-2  -8.33333E-2  -8.33333E-2  -8.33333E-2
0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0
-8.33333E-2  -8.33333E-2  -8.33333E-2   8.33333E-2   8.33333E-2   8.33333E-2
8.33333E-2   8.33333E-2   8.33333E-2  -8.33333E-2  -8.33333E-2  -8.33333E-2
0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0   0.00000E+0
-8.33333E-2  -8.33333E-2  -8.33333E-2   8.33333E-2   8.33333E-2   8.33333E-2

!!! settings "Cleanup"

END PROGRAM main

STConvectiveMatrix example 31

Modules and classes​

Usage​

Modules and classes

Usage