general momentum source/hydraulic pressure jump/positive displacement machine
I need some help in setting up a momentum source. I want to model an internal flow system that contains throttles and a positive
displacement machine. I want to see the effects of the positive displacement machine in the system, however I'm not interested on modeling the actual machine in 3D but a reduced mathematical model. I have the machine performance curves (speed, pressure ratio, volume flowrate, icentropic efficiency etc. Therefore I have created a general momentum source with a simple array that relates dp with flowrate for testing. As for BC inlet is relative pressure 0 [Pa] and outlet's velocity 5 [m/s] (will be massflow later to match the real problem setup but should be ok for now) FUNCTION: DpQ Argument Units = [Pa] Option = Interpolation Profile Function = Off Result Units = [m^3/s] INTERPOLATION DATA: Data Pairs = 0,0.025,500,0.05,1000,0.075 Extend Max = On Extend Min = On Option = One Dimensional END END #density = 1.284 [kg/m^3] #duct area = 0.05 [m^2] Pin = massFlowAve(p)@F91.89 2 Pout = massFlowAve(p)@F96.89 2 deltaP = abs(PoutPin) velocity = volflow/area@F96.89 2 volflow = DpQ(VardeltaP) #the variable VardeltaP is the algebraic expression deltaP scoef5 = 10e5 [kg m^3 s^1] svalue5x = scoef5 * (Velocity u  0[m/s]) svalue5y = scoef5 * (Velocity v  1[m/s]) svalue5z = scoef5 * (Velocity w  0[m/s]) When I use the above basic source definition I can get the thing to converge with the specified velocity in the defined subdomain but when I use the equations below there are many issues: scoef5b = Density * volflow / volume()@B89 2 svalue5bx = scoef5b * (Velocity u  0[m/s]) svalue5by = scoef5b * (Velocity v  velocity) svalue5bz = scoef5b * (Velocity w  0[m/s]) a) for example cfxpost wont let me evaluate the velocity and volflow scalars so I dont know if my coefficient is calculated correctly. b) the problem converges but velocity has the value from the outlet boundary and pressure's are..... well wrong IN = 1.655e+01 [Pa] OUT= 1.547e+02 [Pa] dp across subdomain = 1.381e+02 [Pa] c) am I doing something wrong with my problem approach? I would appreciate any comments especially if the source coefficient is properly defined and if the expression volflow = DpQ(VardeltaP) is correct. many thanks 
Hi,
I would not use "velocity" as a CEL variable. It might get corrupted with the internal velocity variable. Use a different variable name. Glenn Horrocks 
1 Attachment(s)
Quote:
However in cfxpost I do get the follwing error "The following unrecognised name was referenced: DpQ." thats the only bit that I'm unsure, the definition of the volflow expression as the function calculates at the pressure diferential at the inlet and outlet of the subdomain and then it should look up the pressure diferential from the array and output the volume flowrate. I attach the full ccl if anyone has the time to check and comment 
Hi,
Some issues I can see: You refer to "Density" in scoef5b. This is not defined. You define velocityx = volflow/area@F96.89 2. It should be area()@F96.89 2 You have no materials defined. Did you snip that out of the CCL? You need to have some materials defined. Glenn Horrocks 
materials is the stanadrd air @ 25 [C] mate (was a bit lazy :) ) I just copy paste the expression ccl code and forgot materials but use the standard ones for the example
Density is the variable name for density of the fluid m8 :confused: EDIT: yes stupid of me density but where its located thanks I got it now just checking the area()@ :o 
fixed the mistakes and cfx finds a solution but I still get still wrong results which means there is still something else wrong as I still get the "The following unrecognised name was referenced: DpQ." in cfx post
The 1D array is there but maybe there is another way to write the volflow = DpQ(VardeltaP) expression :confused: 
still bugging me this flippin expression :(

Hi,
Sorry, did not read your initial post fully  I don't think CEL interpolation expressions are sent to CFDPost. The best way to check where on your curve your function is lying is to set up monitor points for the input and output of the expression then you can export them from Solver Manager and check them externally. Glenn Horrocks 
thanks mate will try this tomo..
if I remember correcly fluent had this setup as ready made BC ... cant be that hard of a thing to do really... cant see any other way to model loads of fans/pumps inside a system using the machine characteristics... 
Hello Sir,
I'm going through the same problem at the moment. Just curious to see if your model worked in the end. Best regards, Phil P.S: CFX is a pain with fans !! 
Several post on this forum talk about this technique to model fans. It is quite commonly done and works fine when set up correctly.

All times are GMT 4. The time now is 19:49. 