Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.
Solving Linearized Euler Equation in frequency space using Weak Form PDE module
Posted Jul 26, 2019, 10:23 p.m. EDT Fluid & Heat, Acoustics & Vibrations, Equation-Based Modeling Version 5.4 0 Replies
Please login with a confirmed email address before reporting spam
Hi,
I am trying to solve the Linearized Euler Equation (LEE) in frequency space using 'Weak form PDE' module. For 1-D, uniform mean quantity assumptions, the LEE is given by:
Momentum: rho_mean(iomegau+u_meanux)+px = 0 Energy: iomegap+u_meanpx+gammap_mean*ux = 0
If the BCs are rigid at both ends, the first theoretical eigen-frequency is given by:
f1 = c/(2L(1-M^2))
This works perfect when there is no mean flow, i.e, M=0, (eigenfrequency = 1075.7 rad/s with c = 343m/s, L = 1m). However, when I include the mean flow effect (let's say, M=sqrt(10)), it gives a different answer (968.17 rad/s) from the theoretical one (1195.2 rad/s). And their difference from the solution without mean flow is the same:
1075.7/968.17 = 1.1111 1195.2/1075.7 = 1.1111
I guess COMSOL is subtracting something instead of adding it, but not sure exactly.
The attached files contain the ppt that summarized the procedure I took as well as the COMSOL file.
Thanks,
Attachments:
Hello JW Kim
Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.
If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.