physicalProperties.f90

Go to the documentation of this file.

module physicalproperties
    use constants
    use fmodena
    implicit none
    private
    public polymerconductivity,gasconductivity,strutcontent
contains
!********************************BEGINNING*************************************
subroutine polymerconductivity(ksol,temp)
    real(dp), intent(out) :: ksol
    real(dp), intent(in) :: temp
    !modena variables
    integer(c_size_t) :: ksolTemppos

    integer(c_int) :: ret

    type(c_ptr) :: ksolModena = c_null_ptr
    type(c_ptr) :: ksolInputs = c_null_ptr
    type(c_ptr) :: ksolOutputs = c_null_ptr
    ksolmodena = modena_model_new(&
        c_char_"polymer_thermal_conductivity"//c_null_char)
    if (modena_error_occurred()) then
        call exit(modena_error())
    endif
    ksolinputs = modena_inputs_new(ksolmodena)
    ksoloutputs = modena_outputs_new(ksolmodena)
    ksoltemppos = modena_model_inputs_argpos(&
        ksolmodena, c_char_"T"//c_null_char)
    call modena_model_argpos_check(ksolmodena)
    call modena_inputs_set(ksolinputs, ksoltemppos, temp)
    ret = modena_model_call(ksolmodena, ksolinputs, ksoloutputs)
    if(ret /= 0) then
        call exit(ret)
    endif
    ksol=modena_outputs_get(ksoloutputs, 0_c_size_t)
    call modena_inputs_destroy (ksolinputs)
    call modena_outputs_destroy (ksoloutputs)
    call modena_model_destroy (ksolmodena)
end subroutine polymerconductivity
!***********************************END****************************************


!********************************BEGINNING*************************************
subroutine gasconductivity(kgas,temp,xO2,xN2,xCO2,xCyP)
    real(dp), intent(out) :: kgas
    real(dp), intent(in) :: temp
    real(dp), intent(in) :: xO2
    real(dp), intent(in) :: xN2
    real(dp), intent(in) :: xCO2
    real(dp), intent(in) :: xCyP
    !modena variables
    integer(c_size_t) :: kgasTemppos
    integer(c_size_t) :: kgasXCO2pos
    integer(c_size_t) :: kgasXCyPpos
    integer(c_size_t) :: kgasXO2pos
    integer(c_size_t) :: kgasXN2pos

    integer(c_int) :: ret

    type(c_ptr) :: kgasModena = c_null_ptr
    type(c_ptr) :: kgasInputs = c_null_ptr
    type(c_ptr) :: kgasOutputs = c_null_ptr
    kgasmodena = modena_model_new(&
        c_char_"gasMixtureConductivity"//c_null_char)
    if (modena_error_occurred()) then
        call exit(modena_error())
    endif
    kgasinputs = modena_inputs_new(kgasmodena)
    kgasoutputs = modena_outputs_new(kgasmodena)
    kgastemppos = modena_model_inputs_argpos(&
        kgasmodena, c_char_"T"//c_null_char)
    kgasxco2pos = modena_model_inputs_argpos(&
        kgasmodena, c_char_"x[A=CO2]"//c_null_char)
    kgasxcyppos = modena_model_inputs_argpos(&
        kgasmodena, c_char_"x[A=CyP]"//c_null_char)
    kgasxo2pos = modena_model_inputs_argpos(&
        kgasmodena, c_char_"x[A=O2]"//c_null_char)
    kgasxn2pos = modena_model_inputs_argpos(&
        kgasmodena, c_char_"x[A=N2]"//c_null_char)
    call modena_model_argpos_check(kgasmodena)
    call modena_inputs_set(kgasinputs, kgastemppos, temp)
    call modena_inputs_set(kgasinputs, kgasxco2pos, xco2)
    call modena_inputs_set(kgasinputs, kgasxcyppos, xcyp)
    call modena_inputs_set(kgasinputs, kgasxo2pos, xo2)
    call modena_inputs_set(kgasinputs, kgasxn2pos, xn2)
    ret = modena_model_call(kgasmodena, kgasinputs, kgasoutputs)
    if(ret /= 0) then
        call exit(ret)
    endif
    kgas=modena_outputs_get(kgasoutputs, 0_c_size_t)
    call modena_inputs_destroy (kgasinputs)
    call modena_outputs_destroy (kgasoutputs)
    call modena_model_destroy (kgasmodena)
end subroutine gasconductivity
!***********************************END****************************************


!********************************BEGINNING*************************************
subroutine strutcontent(strut_content,foam_density)
    real(dp), intent(out) :: strut_content
    real(dp), intent(in) :: foam_density
    !modena variables
    integer(c_size_t) :: fspos
    integer(c_size_t) :: rhopos

    integer(c_int) :: ret

    type(c_ptr) :: fsModena = c_null_ptr
    type(c_ptr) :: fsInputs = c_null_ptr
    type(c_ptr) :: fsOutputs = c_null_ptr
    fsmodena = modena_model_new(c_char_"strutContent"//c_null_char)
    if (modena_error_occurred()) then
        call exit(modena_error())
    endif
    fsinputs = modena_inputs_new(fsmodena)
    fsoutputs = modena_outputs_new(fsmodena)
    rhopos = modena_model_inputs_argpos(&
        fsmodena, c_char_"rho"//c_null_char)
    call modena_model_argpos_check(fsmodena)
    call modena_inputs_set(fsinputs, rhopos, foam_density)
    ret = modena_model_call(fsmodena, fsinputs, fsoutputs)
    if(ret /= 0) then
        call exit(ret)
    endif
    strut_content=modena_outputs_get(fsoutputs, 0_c_size_t)
    call modena_inputs_destroy (fsinputs)
    call modena_outputs_destroy (fsoutputs)
    call modena_model_destroy (fsmodena)
end subroutine strutcontent
!***********************************END****************************************
end module physicalproperties