#include <hemokunzelmat.h>

Inheritance diagram for oofem::HeMoKunzelMaterial:

Collaboration diagram for oofem::HeMoKunzelMaterial:

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

virtual	~HeMoKunzelMaterial ()
	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 atTime)
	Computes characteristic matrix of receiver in given integration point. More...

virtual double	giveCharacteristicValue (MatResponseMode mode, GaussPoint gp, TimeStep atTime)
	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	sorptionIsotherm (double h)
	returns water content (in kg/m^3) More...

double	sorptionIsothermDerivative (double h)
	computes derivative of the moisture storage function (sorption isotherm) with respect to relative humidity More...

double	computeWaterVaporPerm (double T)

double	computeSatVaporPressure (double T)

double	computeSatVaporPressureDerivative (double T)

double	computeDw (double h)

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 Types
enum	isothermType { Hansen, Kunzeliso }

enum	permeabilityType { Multilin_h, Multilin_wV, Kunzelperm }

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

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

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

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

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

double	perm_mm (double h, double T)

double	perm_mh (double h, double T)

double	perm_hm (double h, double T)

double	perm_hh (double h, double T)

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

virtual double	giveHumidity (GaussPoint *gp, ValueModeType mode)
	Returns positive value of humidity if implemented and enabled in derived material, -1 otherwise. More...

Protected Attributes
enum oofem::HeMoKunzelMaterial::isothermType	Isotherm

enum oofem::HeMoKunzelMaterial::permeabilityType	Permeability

FloatArray	perm_h
	values of the multilinear permeability More...

FloatArray	perm_Dwh

FloatArray	perm_wV

FloatArray	perm_DwwV

double	A
	water absorption coefficient [kg m^-2 s^-0.5] More...

double	iso_wh
	Parameter of Hansen's/Kunzel's isotherm - max. adsorbed water content [kg/m^3]. More...

double	iso_n
	Parameter of Hansen's isotherm - exponent. More...

double	iso_a
	Parameter of Hansen's isotherm. More...

double	iso_b
	Parameter of Kunzel's isotherm. More...

double	mu
	water vapor diffusion resistance [-] More...

double	PL
	ambient atmospheric pressure [Pa] More...

double	lambda0
	thermal conductivity [W m^-1 K^-1] More...

double	b
	thermal conductivity supplement [-] More...

double	rho
	bulk density of dry building material [kg m^-3] More...

double	rhoH2O
	water density [kg m^-3] More...

double	cs
	specific heat capacity of the building material [J kg^-1 K^-1] More...

double	cw
	specific heat capacity of liquid water [J kg^-1 K^-1] More...

double	hv
	latent heat of phase change/evaporation enthalpy of pure water [J/kg] 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

Definition at line 67 of file hemokunzelmat.h.

Member Enumeration Documentation

enum oofem::HeMoKunzelMaterial::isothermType

protected

Enumerator
Hansen
Kunzeliso

Definition at line 71 of file hemokunzelmat.h.

enum oofem::HeMoKunzelMaterial::permeabilityType

protected

Enumerator
Multilin_h
Multilin_wV
Kunzelperm

Definition at line 73 of file hemokunzelmat.h.

Constructor & Destructor Documentation

oofem::HeMoKunzelMaterial::HeMoKunzelMaterial	(	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 109 of file hemokunzelmat.h.

virtual oofem::HeMoKunzelMaterial::~HeMoKunzelMaterial ( )

inlinevirtual

Destructor.

Definition at line 111 of file hemokunzelmat.h.

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

Member Function Documentation

double oofem::HeMoKunzelMaterial::computeCapacityCoeff	(	MatResponseMode	mode,
		GaussPoint *	gp,
		TimeStep *	atTime
	)

protected

Definition at line 344 of file hemokunzelmat.C.

References oofem::FloatArray::at(), cs, cw, give(), oofem::Material::giveStatus(), oofem::TransportMaterialStatus::giveTempField(), oofem::FloatArray::isEmpty(), OOFEM_ERROR, sorptionIsotherm(), and sorptionIsothermDerivative().

Referenced by giveCharacteristicValue(), and giveClassName().

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

protected

Definition at line 214 of file hemokunzelmat.C.

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

Referenced by giveCharacteristicMatrix(), and giveClassName().

double oofem::HeMoKunzelMaterial::computeDw ( double h )

Definition at line 499 of file hemokunzelmat.C.

References A, oofem::FloatArray::at(), oofem::FloatArray::giveSize(), iso_wh, Kunzelperm, Multilin_h, Multilin_wV, OOFEM_ERROR, perm_Dwh, perm_DwwV, perm_h, perm_wV, Permeability, rhoH2O, and sorptionIsotherm().

Referenced by giveClassName(), and perm_mm().

double oofem::HeMoKunzelMaterial::computeSatVaporPressure ( double T )

Definition at line 458 of file hemokunzelmat.C.

Referenced by giveClassName(), perm_hm(), and perm_mm().

double oofem::HeMoKunzelMaterial::computeSatVaporPressureDerivative ( double T )

Definition at line 478 of file hemokunzelmat.C.

Referenced by giveClassName(), perm_hh(), and perm_mh().

double oofem::HeMoKunzelMaterial::computeWaterVaporPerm ( double T )

vapor diffusion coefficient in air [kg m^-1 s^-1 Pa^-1]

Definition at line 445 of file hemokunzelmat.C.

References mu, and PL.

Referenced by giveClassName(), perm_hh(), perm_hm(), perm_mh(), and perm_mm().

double oofem::HeMoKunzelMaterial::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 143 of file hemokunzelmat.C.

References oofem::Material::give().

Referenced by computeCapacityCoeff(), perm_hh(), and ~HeMoKunzelMaterial().

void oofem::HeMoKunzelMaterial::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 189 of file hemokunzelmat.C.

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

Referenced by ~HeMoKunzelMaterial().

double oofem::HeMoKunzelMaterial::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.

Definition at line 206 of file hemokunzelmat.C.

References computeCapacityCoeff().

Referenced by ~HeMoKunzelMaterial().

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

inlinevirtual

Returns: Class name of the receiver.

Implements oofem::FEMComponent.

Definition at line 134 of file hemokunzelmat.h.

References computeCapacityCoeff(), computeConductivityMtrx(), computeDw(), computeSatVaporPressure(), computeSatVaporPressureDerivative(), computeWaterVaporPerm(), giveHumidity(), giveIPValue(), matcond1d(), matcond2d(), matcond3d(), perm_hh(), perm_hm(), perm_mh(), perm_mm(), sorptionIsotherm(), and sorptionIsothermDerivative().

void oofem::HeMoKunzelMaterial::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 154 of file hemokunzelmat.C.

References oofem::FloatArray::add(), oofem::FloatArray::addSubVector(), oofem::FloatArray::at(), oofem::FloatArray::giveSize(), oofem::Material::giveStatus(), perm_hh(), perm_hm(), perm_mh(), perm_mm(), oofem::FloatArray::resize(), oofem::TransportMaterialStatus::setTempField(), oofem::TransportMaterialStatus::setTempFlux(), oofem::TransportMaterialStatus::setTempGradient(), and oofem::FloatArray::zero().

Referenced by ~HeMoKunzelMaterial().

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

protectedvirtual

Returns positive value of humidity if implemented and enabled in derived material, -1 otherwise.

Reimplemented from oofem::TransportMaterial.

Definition at line 645 of file hemokunzelmat.C.

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

Referenced by giveClassName(), and giveIPValue().

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

inlinevirtual

Returns: Input record name of the receiver.

Implements oofem::FEMComponent.

Definition at line 133 of file hemokunzelmat.h.

References _IFT_HeMoKunzelMaterial_Name.

int oofem::HeMoKunzelMaterial::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.

Definition at line 632 of file hemokunzelmat.C.

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

Referenced by giveClassName().

int oofem::HeMoKunzelMaterial::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 hemokunzelmat.C.

Referenced by ~HeMoKunzelMaterial().

IRResultType oofem::HeMoKunzelMaterial::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.