This class implements a coupled heat and mass transfer material model. More...

#include <hemotkmat.h>

Inheritance diagram for oofem::HeMoTKMaterial:

Collaboration diagram for oofem::HeMoTKMaterial:

Public Member Functions
	HeMoTKMaterial (int n, Domain *d)
	Constructor. More...

virtual	~HeMoTKMaterial ()
	Destructor. More...

virtual void	giveFluxVector (FloatArray &answer, GaussPoint gp, const FloatArray &grad, const FloatArray &field, TimeStep tStep)
	Returns the flux for the field and its gradient. More...

virtual void	giveCharacteristicMatrix (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes characteristic matrix of receiver in given integration point. More...

virtual double	giveCharacteristicValue (MatResponseMode mode, GaussPoint gp, TimeStep tStep)
	Computes the characteristic value of receiver in given integration point, respecting its history. More...

virtual bool	isCharacteristicMtrxSymmetric (MatResponseMode rMode)
	Returns true if stiffness matrix of receiver is symmetric Default implementation returns true. More...

virtual IRResultType	initializeFrom (InputRecord *ir)
	Initializes receiver according to object description stored in input record. More...

virtual double	give (int aProperty, GaussPoint *gp)
	Returns the value of material property 'aProperty'. More...

virtual int	hasMaterialModeCapability (MaterialMode mode)
	Tests if material supports material mode. More...

virtual const char *	giveInputRecordName () const

virtual const char *	giveClassName () const

double	sorption_isotherm (double phi)

double	inverse_sorption_isotherm (double w)

double	give_dphi_dw (double w)

Public Member Functions inherited from oofem::TransportMaterial
	TransportMaterial (int n, Domain *d)
	Constructor. More...

virtual	~TransportMaterial ()
	Destructor. More...

virtual void	updateInternalState (const FloatArray &state, GaussPoint gp, TimeStep tStep)
	Updates internal state of material according to new state vector. More...

virtual int	hasInternalSource ()
	Returns nonzero if receiver generates internal source of state variable(s), zero otherwise. More...

virtual void	computeInternalSourceVector (FloatArray &val, GaussPoint gp, TimeStep tStep, ValueModeType mode)
	Computes the internal source vector of receiver. More...

virtual MaterialStatus *	CreateStatus (GaussPoint *gp) const
	Creates new copy of associated status and inserts it into given integration point. More...

Public Member Functions inherited from oofem::Material
	Material (int n, Domain *d)
	Constructor. More...

virtual	~Material ()
	Destructor. More...

virtual bool	hasProperty (int aProperty, GaussPoint *gp)
	Returns true if 'aProperty' exists on material. More...

virtual void	modifyProperty (int aProperty, double value, GaussPoint *gp)
	Modify 'aProperty', which already exists on material. More...

double	giveCastingTime ()

virtual bool	isActivated (TimeStep *tStep)

virtual int	hasNonLinearBehaviour ()
	Returns nonzero if receiver is non linear. More...

virtual int	hasCastingTimeSupport ()
	Tests if material supports casting time. More...

virtual void	giveInputRecord (DynamicInputRecord &input)
	Setups the input record string of receiver. More...

virtual void	printYourself ()
	Prints receiver state on stdout. Useful for debugging. More...

virtual contextIOResultType	saveIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
	Stores integration point state to output stream. More...

virtual contextIOResultType	restoreIPContext (DataStream &stream, ContextMode mode, GaussPoint *gp)
	Reads integration point state to output stream. More...

virtual int	checkConsistency ()
	Allows programmer to test some internal data, before computation begins. More...

virtual int	initMaterial (Element *element)
	Optional function to call specific procedures when initializing a material. More...

virtual MaterialStatus *	giveStatus (GaussPoint *gp) const
	Returns material status of receiver in given integration point. More...

virtual int	packUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
	Pack all necessary data of integration point (according to element parallel_mode) into given communication buffer. More...

virtual int	unpackAndUpdateUnknowns (DataStream &buff, TimeStep tStep, GaussPoint ip)
	Unpack and updates all necessary data of given integration point (according to element parallel_mode) into given communication buffer. More...

virtual int	estimatePackSize (DataStream &buff, GaussPoint *ip)
	Estimates the necessary pack size to hold all packed data of receiver. More...

virtual double	predictRelativeComputationalCost (GaussPoint *gp)
	Returns the weight representing relative computational cost of receiver The reference material model is linear isotropic material - its weight is set to 1.0 The other material models should compare to this reference model. More...

virtual double	predictRelativeRedistributionCost (GaussPoint *gp)
	Returns the relative redistribution cost of the receiver. More...

virtual void	initTempStatus (GaussPoint *gp)
	Initializes temporary variables stored in integration point status at the beginning of new time step. More...

virtual int	setIPValue (const FloatArray &value, GaussPoint *gp, InternalStateType type)
	Sets the value of a certain variable at a given integration point to the given value. More...

Public Member Functions inherited from oofem::FEMComponent
	FEMComponent (int n, Domain *d)
	Regular constructor, creates component with given number and belonging to given domain. More...

virtual	~FEMComponent ()
	Virtual destructor. More...

Domain *	giveDomain () const

virtual void	setDomain (Domain *d)
	Sets associated Domain. More...

int	giveNumber () const

void	setNumber (int num)
	Sets number of receiver. More...

virtual void	updateLocalNumbering (EntityRenumberingFunctor &f)
	Local renumbering support. More...

virtual contextIOResultType	saveContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Stores receiver state to output stream. More...

virtual contextIOResultType	restoreContext (DataStream &stream, ContextMode mode, void *obj=NULL)
	Restores the receiver state previously written in stream. More...

virtual void	printOutputAt (FILE file, TimeStep tStep)
	Prints output of receiver to stream, for given time step. More...

virtual Interface *	giveInterface (InterfaceType t)
	Interface requesting service. More...

std::string	errorInfo (const char *func) const
	Returns string for prepending output (used by error reporting macros). More...

Protected Member Functions
void	computeConductivityMtrx (FloatMatrix &answer, MatResponseMode mode, GaussPoint gp, TimeStep tStep)

void	matcond1d (FloatMatrix &d, GaussPoint gp, MatResponseMode mode, TimeStep tStep)

void	matcond2d (FloatMatrix &d, GaussPoint gp, MatResponseMode mode, TimeStep tStep)

void	matcond3d (FloatMatrix &d, GaussPoint gp, MatResponseMode mode, TimeStep tStep)

double	computeCapacityCoeff (MatResponseMode mode, GaussPoint gp, TimeStep tStep)

double	giveHumidity (GaussPoint *gp, ValueModeType mode)
	Returns positive value of humidity, use VM_Velocity for previous (equilibrated) value. More...

double	get_latent (double w, double t)

double	get_ceff (double w, double t)

double	get_chi (double w, double t)

double	perm_wt (double w, double t)

double	perm_ww (double w, double t)

double	give_delta_gw (double phi)

double	give_dpgw_dt (double t, double phi)

double	get_b (double w, double t)

double	get_sat (double w, double t)

double	give_p_gws (double t)

virtual int	giveIPValue (FloatArray &answer, GaussPoint gp, InternalStateType type, TimeStep tStep)
	Returns the integration point corresponding value in Reduced form. More...

Protected Attributes
double	a_0
	Constant (obtained from experiments) [Bazant and Najjar, 1972]. More...

double	nn
	Constant-exponent (obtained from experiments) [Bazant and Najjar, 1972]. More...

double	phi_c
	Constant-relative humidity (obtained from experiments) [Bazant and Najjar, 1972]. More...

double	delta_wet
	Constant-water vapor permeability (obtained from experiments) [Bazant and Najjar, 1972]. More...

double	w_h
	Constant water content (obtained from experiments) [Pedersen, 1990]. More...

double	n
	Constant-exponent (obtained from experiments) [Pedersen, 1990]. More...

double	a
	Constant (obtained from experiments) [Pedersen, 1990]. More...

double	latent
	Latent heat of evaporation. More...

double	c
	Thermal capacity. More...

double	rho
	Volume density. More...

double	chi_eff
	Effective thermal conductivity. More...

double	por
	Porosity. More...

double	rho_gws
	Saturation volume density. More...

Protected Attributes inherited from oofem::Material
Dictionary	propertyDictionary
	Property dictionary. More...

double	castingTime
	Casting time. More...

Protected Attributes inherited from oofem::FEMComponent
int	number
	Component number. More...

Domain *	domain
	Link to domain object, useful for communicating with other FEM components. More...

Detailed Description

This class implements a coupled heat and mass transfer material model.

It computes conductivity and capacity matrices for coupled heat and moisture transfer. Assumptions: Water vapor is the only driving mechanism; relative humidity is from range 0.2 - 0.98 (I and II region).

Author: Tomas Krejci Source: T.K. Doctoral Thesis; Bazant and Najjar, 1972; Pedersen, 1990.

Definition at line 66 of file hemotkmat.h.

Constructor & Destructor Documentation

oofem::HeMoTKMaterial::HeMoTKMaterial	(	int	n,
		Domain *	d
	)

inline

Constructor.

Creates material with given number, belonging to given domain.

Parameters

n	Material number.
d	Domain to which new material will belong.

Definition at line 93 of file hemotkmat.h.

virtual oofem::HeMoTKMaterial::~HeMoTKMaterial ( )

inlinevirtual

Destructor.

Definition at line 95 of file hemotkmat.h.

References give(), giveCharacteristicMatrix(), giveCharacteristicValue(), giveFluxVector(), hasMaterialModeCapability(), initializeFrom(), and isCharacteristicMtrxSymmetric().

Member Function Documentation

double oofem::HeMoTKMaterial::computeCapacityCoeff	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protected

Definition at line 298 of file hemotkmat.C.

References oofem::FloatArray::at(), get_b(), get_ceff(), get_latent(), oofem::Material::giveStatus(), oofem::TransportMaterialStatus::giveTempField(), oofem::FloatArray::isEmpty(), OOFEM_ERROR, and rho.

Referenced by giveCharacteristicValue(), and giveClassName().

void oofem::HeMoTKMaterial::computeConductivityMtrx	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

protected

Definition at line 151 of file hemotkmat.C.

References oofem::GaussPoint::giveMaterialMode(), matcond2d(), matcond3d(), and OOFEM_ERROR.

Referenced by giveCharacteristicMatrix(), and giveClassName().

double oofem::HeMoTKMaterial::get_b	(	double	w,
		double	t
	)

protected

Definition at line 536 of file hemotkmat.C.

References get_sat(), give_dphi_dw(), inverse_sorption_isotherm(), por, and rho_gws.

Referenced by computeCapacityCoeff(), and giveClassName().

double oofem::HeMoTKMaterial::get_ceff	(	double	w,
		double	t
	)

protected

Definition at line 567 of file hemotkmat.C.

References c, and rho.

Referenced by computeCapacityCoeff(), and giveClassName().

double oofem::HeMoTKMaterial::get_chi	(	double	w,
		double	t
	)

protected

Definition at line 580 of file hemotkmat.C.

References chi_eff.

Referenced by giveClassName(), giveFluxVector(), matcond1d(), matcond2d(), and matcond3d().

double oofem::HeMoTKMaterial::get_latent	(	double	w,
		double	t
	)

protected

Definition at line 524 of file hemotkmat.C.

References latent.

Referenced by computeCapacityCoeff(), giveClassName(), giveFluxVector(), matcond1d(), matcond2d(), and matcond3d().

double oofem::HeMoTKMaterial::get_sat	(	double	w,
		double	t
	)

protected

Definition at line 554 of file hemotkmat.C.

Referenced by get_b(), and giveClassName().

double oofem::HeMoTKMaterial::give	(	int	aProperty,
		GaussPoint *	gp
	)

virtual

Returns the value of material property 'aProperty'.

Property must be identified by unique int id. Integration point also passed to allow for materials with spatially varying properties

Parameters

aProperty	ID of property requested.
gp	Integration point,

Returns: Property value.

Reimplemented from oofem::Material.

Definition at line 82 of file hemotkmat.C.

References oofem::Material::give().

Referenced by oofem::HydratingHeMoMaterial::updateInternalState(), and ~HeMoTKMaterial().

double oofem::HeMoTKMaterial::give_delta_gw ( double phi )

protected

Definition at line 402 of file hemotkmat.C.

References a_0, delta_wet, nn, OOFEM_ERROR, and phi_c.

Referenced by giveClassName(), perm_wt(), and perm_ww().

double oofem::HeMoTKMaterial::give_dpgw_dt	(	double	t,
		double	phi
	)

protected

Definition at line 489 of file hemotkmat.C.

Referenced by giveClassName(), and perm_wt().

double oofem::HeMoTKMaterial::give_dphi_dw ( double w )

Definition at line 471 of file hemotkmat.C.

References a, n, and w_h.

Referenced by get_b(), oofem::HydratingHeMoMaterial::giveCharacteristicValue(), giveClassName(), and perm_ww().

double oofem::HeMoTKMaterial::give_p_gws ( double t )

protected

Definition at line 508 of file hemotkmat.C.

Referenced by giveClassName(), and perm_ww().

void oofem::HeMoTKMaterial::giveCharacteristicMatrix	(	FloatMatrix &	answer,
		MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes characteristic matrix of receiver in given integration point.

The algorithm should use temporary or equilibrium history variables stored in integration point status to compute and return required result.

Parameters

answer	Contains result.
mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::TransportMaterial.

Definition at line 126 of file hemotkmat.C.

References oofem::__MatResponseModeToString(), computeConductivityMtrx(), and OOFEM_ERROR.

Referenced by ~HeMoTKMaterial().

double oofem::HeMoTKMaterial::giveCharacteristicValue	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	tStep
	)

virtual

Computes the characteristic value of receiver in given integration point, respecting its history.

The algorithm should use temporary or equilibrium history variables stored in integration point status to compute and return required result.

Parameters

mode	Material response mode.
gp	Integration point.
tStep	Time step (most models are able to respond only when tStep is current time step).

Implements oofem::TransportMaterial.

Reimplemented in oofem::HydratingHeMoMaterial.

Definition at line 143 of file hemotkmat.C.

References computeCapacityCoeff().

Referenced by oofem::HydratingHeMoMaterial::giveCharacteristicValue(), and ~HeMoTKMaterial().

virtual const char* oofem::HeMoTKMaterial::giveClassName ( ) const

inlinevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Reimplemented in oofem::HydratingHeMoMaterial.

Definition at line 118 of file hemotkmat.h.

References computeCapacityCoeff(), computeConductivityMtrx(), get_b(), get_ceff(), get_chi(), get_latent(), get_sat(), give_delta_gw(), give_dpgw_dt(), give_dphi_dw(), give_p_gws(), giveHumidity(), giveIPValue(), inverse_sorption_isotherm(), matcond1d(), matcond2d(), matcond3d(), perm_wt(), perm_ww(), and sorption_isotherm().

void oofem::HeMoTKMaterial::giveFluxVector	(	FloatArray &	answer,
		GaussPoint *	gp,
		const FloatArray &	grad,
		const FloatArray &	field,
		TimeStep *	tStep
	)

virtual

Returns the flux for the field and its gradient.

Todo:: { Should the field variable just be a scalar? This might change when we rethink the coupled-fields approach. Now its either just [temperature], or [temperature, concentration] so to cover both cases there is a floatarray. }

Parameters

answer	The flux.
gp	Gauss point.
grad	Gradient of the primary field, usually the main input.
field	The value of the field itself.
tStep	Active time step.

Implements oofem::TransportMaterial.

Definition at line 93 of file hemotkmat.C.

References oofem::FloatArray::addSubVector(), oofem::FloatArray::at(), oofem::FloatArray::beScaled(), get_chi(), get_latent(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), perm_wt(), perm_ww(), oofem::FloatArray::resize(), oofem::TransportMaterialStatus::setTempField(), oofem::TransportMaterialStatus::setTempFlux(), oofem::TransportMaterialStatus::setTempGradient(), and oofem::FloatArray::zero().

Referenced by ~HeMoTKMaterial().

double oofem::HeMoTKMaterial::giveHumidity	(	GaussPoint *	gp,
		ValueModeType	mode
	)

protectedvirtual

Returns positive value of humidity, use VM_Velocity for previous (equilibrated) value.

Reimplemented from oofem::TransportMaterial.

Definition at line 339 of file hemotkmat.C.

References oofem::FloatArray::at(), oofem::TransportMaterialStatus::giveField(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), oofem::TransportMaterialStatus::giveTempField(), inverse_sorption_isotherm(), and OOFEM_ERROR.

Referenced by giveClassName(), and giveIPValue().

virtual const char* oofem::HeMoTKMaterial::giveInputRecordName ( ) const

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Reimplemented in oofem::HydratingHeMoMaterial.

Definition at line 117 of file hemotkmat.h.

References _IFT_HeMoTKMaterial_Name.

int oofem::HeMoTKMaterial::giveIPValue	(	FloatArray &	answer,
		GaussPoint *	gp,
		InternalStateType	type,
		TimeStep *	tStep
	)

protectedvirtual

Returns the integration point corresponding value in Reduced form.

Parameters

answer	Contain corresponding ip value, zero sized if not available.
gp	Integration point to which the value refers.
type	Determines the type of internal variable.
tStep	Determines the time step.

Returns: Nonzero if the assignment can be done, zero if this type of variable is not supported.

Reimplemented from oofem::TransportMaterial.

Reimplemented in oofem::HydratingHeMoMaterial.

Definition at line 604 of file hemotkmat.C.

References oofem::FloatArray::at(), giveHumidity(), oofem::TransportMaterial::giveIPValue(), and oofem::FloatArray::resize().

Referenced by giveClassName(), and oofem::HydratingHeMoMaterial::giveIPValue().

int oofem::HeMoTKMaterial::hasMaterialModeCapability ( MaterialMode mode )

virtual

Tests if material supports material mode.

Parameters

mode	Required material mode.

Returns: Nonzero if supported, zero otherwise.

Reimplemented from oofem::Material.

Definition at line 45 of file hemotkmat.C.

Referenced by ~HeMoTKMaterial().

IRResultType oofem::HeMoTKMaterial::initializeFrom ( InputRecord * ir )

virtual

Initializes receiver according to object description stored in input record.

This function is called immediately after creating object using constructor. Input record can be imagined as data record in component database belonging to receiver. Receiver may use value-name extracting functions to extract particular field from record.