model.f90 File Reference

Physical model of bubble growth. More...

Go to the source code of this file.

Functions/Subroutines
subroutine, public	model::odesystem (neq, t, y, ydot)
	Model supplied to integrator. More...
subroutine, public	model::molar_balance (y)
	Calculates molar amount of blowing agents in bubble and shell. More...
subroutine, public	model::dim_var (t, y)
	Calculates dimensional variables. More...
subroutine, public	model::kinmodel (y)
	Evaluates kinetic source terms. More...

Detailed Description

Physical model of bubble growth.

Author

Pavel Ferkl

Model is a collection of momentum, mass and enthalpy balances.

Definition in file model.f90.

Function/Subroutine Documentation

dim_var()

subroutine, public model::dim_var	(	real(dp), intent(in)	t,
		real(dp), dimension(:), intent(in)	y
	)

Calculates dimensional variables.

Transforms integrated variables to human readable format with easy to remember names.

Parameters

[in]	t	time
[in]	y	integrated variables

Definition at line 239 of file model.f90.

References globals::avconc, globals::conv, globals::eqconc, globals::fpeq, globals::geometry, globals::kh, globals::laplace_pres, globals::mb, globals::mbl, globals::ngas, globals::pair, globals::pair0, globals::pamb, globals::porosity, globals::pressure, globals::r0, globals::radius, globals::req, globals::rhofoam, globals::rhop, globals::s0, globals::sigma, globals::st, globals::temp, globals::temp0, globals::teq, globals::time, globals::vsh, globals::wblpol, and globals::xoheq.

kinmodel()

subroutine, public model::kinmodel

(

real(dp), dimension(:), intent(in)

y

)

Evaluates kinetic source terms.

Uses Modena model. Other kinetic models are implemented directly in the model subroutine.

Parameters

[in]

y

integrated variables

Definition at line 283 of file model.f90.

References globals::kin_model, globals::kineq, globals::kinsource, globals::teq, and globals::time.

molar_balance()

subroutine, public model::molar_balance

(

real(dp), dimension(:), intent(in)

y

)

Calculates molar amount of blowing agents in bubble and shell.

Can be used to check that we don't lose physical blowing agent.

Parameters

[in]

y

integrated variables

Definition at line 206 of file model.f90.

References globals::dz, globals::fceq, globals::geometry, globals::mb, globals::mb2, globals::mb3, globals::ngas, globals::p, globals::pressure, globals::radius, and globals::temp.

odesystem()

subroutine public model::odesystem	(	integer, intent(in)	neq,
		real(dp), intent(in)	t,
		real(dp), dimension(neq), intent(in)	y,
		real(dp), dimension(neq), intent(out)	ydot
	)

Model supplied to integrator.

Implementation of the ODEs.

Bubble growth model. Discretized by FVM, nonequidistant mesh.

Parameters

[in]	neq	number of equations
[in]	t	time
[in]	y	integrated variables
[out]	ydot	derivatives of integrated variables

FVM, equidistant mesh, cylindrical geometry, sin(alpha)=(dh/dr)/(1+(dh/dr)**2)

Parameters

[in]	neq	number of equations
[in]	t	time
[in]	y	integrated variables
[out]	ydot	derivatives of integrated variables

Definition at line 22 of file model.f90.

References globals::aoh, globals::aw, globals::co2_pos, globals::cp, globals::d, globals::dhoh, globals::dhv, globals::dhw, globals::dilution, globals::dz, globals::eoh, globals::eta, globals::ew, globals::fceq, globals::fpeq, globals::geometry, globals::inertial_term, globals::kh, globals::kin_model, globals::kineq, globals::kinsource, globals::lpeq, globals::mb, globals::mbl, globals::nco0, globals::ngas, globals::oh0, globals::p, globals::pair, globals::pamb, globals::radius, globals::req, globals::rhobl, globals::rhop, globals::sigma, globals::temp, globals::teq, globals::time, globals::vsh, globals::w0, globals::xoheq, and globals::xweq.