MATLAB Answers

Simscape 2P energy conservation bug ?

조회 수: 2(최근 30일)
Joris Naudin
Joris Naudin 2021년 3월 8일
댓글: Joris Naudin 2021년 5월 3일
Hello everyone,
I'd like to report something which I think could be a bug ?
I attached a very simple model showcasing the effect presented below.
I have a tank of steam discharging into a translation mechanical converter (so a piston), no heat losses/gains, no friction etc.. It is linked to a spring just so that it reaches an equilibrium after some time.
The piston moves just like expected but the enthalpy and temperature inside it does not: it actually goes higher than the supply !! Whereas initial conditions in the cylinder are colder than the supply.
I can't think of this being possible in real world physics or am I missing something ??
Thanks, have fun !
  댓글 수: 3
Joris Naudin
Joris Naudin 2021년 4월 29일
Impatient to know the reason behind
Thanks !

댓글을 달려면 로그인하십시오.

채택된 답변

Yifeng Tang
Yifeng Tang 2021년 5월 1일
Hi Joris,
After a closer look at the results, it looks like everything was physically reasonable. There was some very quick action that was happening during the first 2e-5 sec and 0.15 sec that created the initial spike in the T_I of the cylinder. In fact, it's more than a simple spike.
0 ~ 2.5e-5 sec: vapor rush into the cylinder (mdot_A), pressure quickly rises (p_I), T_I quick rises because the incoming fluid is doing work to the volume. At this time, there is a huge amount of force applied by the cylinder on the mass because the initial force from the spring is zero. The mass will accelerate quickly.
0 ~ 0.05sec: the mass is accelerated quickly to 30m/s, expanding the cylinder volume. The pressure and temperature drops quickly because mdot_A can't keep up. The net force on the mass crosses zero ~0.035sec as the spring is compressed enough.
0.05 ~ 0.15sec: mass still moving away but slowing down. mdot_A catches up, doing work to the low pressure volume, increasing both the pressure and temperature. Around this time, the spring force starts to match the cylinder force and the volume starts to stablize.
After that, things start to approach steady-state.
During the process, the temperature of the chamber drops a little bit, while also fluctuates according to the sign/direction of mdot_A. That's also where all the energy comes from.
Hope this makes sense.
  댓글 수: 6
Joris Naudin
Joris Naudin 2021년 5월 3일
Interesting ... something to get my head around finally!
I would love to see this phenomenon on a real life test bench.
Thank you very much !

댓글을 달려면 로그인하십시오.

추가 답변(0개)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by