Define domain with multiple material functions

[TLWise]

I am currently using, type = ADGenericConstantMaterial, prop_names = 'rho mu cp k'. I would like to define each of the materials with an analytical function. is there something similar to ADGenericConstantMaterial that allows me to define rho mu cp k, with an analytical function for each?

[GiudGiud]

Hello

you can use GenericFunctionMaterial for functions of time and space.

For analytical functions of things like variables you may use the ParsedMaterial

Guillaume

[GiudGiud]

If you have a simple expression you can write in an input files, the ParsedMaterial

if not then you ll have to make your own material, using C++ code. The tutorial explains how

[TLWise]

It appears that these (Material property 'rho', requested by 'ins_mat' is not defined on block Coolant
Material property 'rho', requested by 'ins_mat_face' is not defined on block Coolant
Material property 'rho', requested by 'ins_mat_neighbor' is not defined on block Coolant
Material property 'rho', requested by 'mass_pspg' is not defined on block Coolant
) expect that material properties be defined in from ADGenericConstantMaterial is preventing my use of ParsedMaterial to specify rho as a function of temperature.

[GiudGiud]

You can only specify a property once.

You also need to be consistent with using AD. So you ll want to use an ADParsedMaterial most likely

[TLWise]

Yep. I figured that would probably be the cause for the issue. I tried the command using AD, and it is what resulted in the previous issue I shared. I removed rho as a definition in the ADConst defined it in the ADParsed. The other portions of AD did not recognize rho as being defined.

[GiudGiud]

can you please paste your former and current Materials block?

[TLWise]

[Materials] # Old
  [./const]
    type = ADGenericConstantMaterial
    prop_names = 'rho mu cp k'
    prop_values = '1000  0.001 4179.0  0.598' #[kg/m^3*1e-9] [Pa*sec*1e-6] [J/kg/K*1e6] [W/m/K*1e3] non exponential forms
  [../]
  [ins_mat]
    type = INSADStabilized3Eqn
    velocity = velocity
    pressure = p
    temperature = T_fluid
  []
[]

[Materials] # New
  [./const]
    type = ADGenericConstantMaterial
    prop_names = 'mu cp k'
#    prop_values = '0.001 4179.0  0.598' #[kg/m^3*1e-9] [Pa*sec*1e-6] [J/kg/K*1e6] [W/m/K*1e3] non exponential forms 
  [../]
  [./rho_density]
    type = ADParsedMaterial
    property_name = 'rho'
    coupled_variables = 'T_fluid'
    function = '99.98*((1.000088*(T_fluid-283.2))/(1-(p-1.0e5)/2.15e9))'
#    outputs = exodus
  [../]
  [ins_mat]
    type = INSADStabilized3Eqn
    velocity = velocity
    pressure = p
    temperature = T_fluid
  []
[]

[GiudGiud]

oh is your MOOSE old?

if so property_name is not the name of the parameter, it s still called f_name in an old version of MOOSE
property_name would return an unused parameter warning

[TLWise]

Yes, the version on this system is old, and I have to switch between expression and function when moving from one computer to the other.

[GiudGiud]

You ll have to use the old parameters. property-name is not recognized, so the default for f_name is being used.
so instead of declaring rho you are declaring a property as F and naturally getting an error

[TLWise]

This is the error that I get when I use T_fluid in the definition of rho:

*** ERROR ***
Invalid function
99.98*((1.000088*(T_fluid-283.2))/(1-(p-1.0e5)/2.15e9))

in ParsedMaterialHelper.
Syntax error: Unknown identifier

[GiudGiud]

you need to adapt all the parameters
All your neew syntax parameters are not recognized

you should have args = 'p T_fluid'

[TLWise]

My apologies upfront for not having mastered how to prevent the super large fonts.

[GiudGiud]

no worries. I use the triple backquotes to do that

[TLWise]

Had to remove the reference to T_fluid, and it got me further in the processing chain. I now get this error instead:

*** ERROR ***
/home/twise/projects/moose/modules/combined/7x7PWR_Assembly/Models/CoolantVelocityPWR.i:217: unused parameter 'Materials/rho_density/coupled_variables'

[GiudGiud]

should be args, not coupled_variables.
You can add T and p back once you have that fixed