I am trying to implement a Casson viscosity model for fluid flow in an artery. The Casson fluid is non-Newtonian in nature and thus depends on the shear rate.
In the material models module I created another case with string casson and compiled the materialmodels.mod file.
Code: Select all
CASE ('casson')
c1n = ListGetReal( Material, 'Shear Yield Strength',n,Element % NodeIndexes )
c1 = SUM( Basis(1:n) * c1n(1:n) )
c2 = ListGetConstReal( Material, 'M Value',gotIt)
IF(gotIt) THEN
s= SQRT(ss)
mu= (SQRT(Viscosity)+SQRT(c1/s)*(1-2.718**(-1*SQRT(c2*s))))**2
END IF
Code: Select all
Material 1
Name = "Material 1"
Procedure "materialmodels" "EffectiveViscosity"
Viscosity Model = casson
M Value = 2
Viscosity = 0.0035
Porosity Model = Always saturated
Shear Yield Strength = 100
End
Code: Select all
MAIN: Reading Model: case.sif
LoadInputFile: Scanning input file: case.sif
ERROR:: Model Input:
ERROR:: Model Input: Unknown specifier: [viscosity model]
ERROR:: Model Input: In section: [material 1]
ERROR:: Model Input: For property name:[procedure materialmodels EffectiveViscosity]