It is not in general possible (at this point) to call WL functions from SystemModeler components.
I agree that it would be a great feature to have.
What you can currently do is use functions that are simple enough to be converted to Modelica (the language used by SystemModeler models) in CreateSystemModel:
model = CreateSystemModel[{Unevaluated[Function[#1 + Sin[#2]]][1,
u[t]] == u'[t]}, t]
While Malte is correct that there is no general way to do this directly, with a little extra effort you can do this.
You can call Mathematica from C using WSTP. See the documentation here. Also the tutorials referenced from that page will describe how to create a C wrapper function that will send any data, run any function, and return results from MMA using C.
The C wrapper function can be called from SystemModeler using an external C block in WSM. The documentation for calling C from WSM can be found in the WSM help system and also in Modelica By Example at this link within the WSM documentation center:
The overall process would be to write a small (few lines of code) C function based on the examples in the WSTP documentation that executes whatever MMA function you want using data sent in from C. Call the function from SystemModeler using the examples above. I hope this helps.
Malte hints in 2018 that it will eventually be possible to call MMA functions from SystemModeler, does anyone know how this is progressing. To date the focus has been on adding SystemModeler’s capabilities into Wolfram Language.