Hi,(adsbygoogle = window.adsbygoogle || []).push({});

For my Bachelor's thesis I've been working on a finite time Carnot cycle.

I've finished my numerical analysis using the differential equations governing the time evolution.

My next step should be a simulation.

First I should stick to a 1 dimensional system.

This system consists of a piston and a thermalising wall.

This thermalising wall acts as a heat reservoir. Every particle colliding with this wall is absorbed.

The wall than ejects a 'new' particle with a certain velocity.

This velocity is governed by the stochastic distribution

[itex]f(\vec{v},T_i)=C v \exp\left(-\frac{m v^2}{2kT_i}\right)[/itex] with [itex]i=h[/itex] while expanding and [itex]i=c[/itex] while compressing. C is the normalisation constant.

Since I use a 1D system in this first approximating step, the Maxwell-Boltzmann distribution isn't necessary.

The piston has a constant velocity u. This is chosen because of the fact that the article I base my calculations on is targeting a system that is easy to control.

The article is "Molecular Kinetic analysis of a finite-time Carnot cycle" by Y. Izumida and K. Okuda published in september 2008.

I reckon I have to use some sort of Monte-Carlo method because of the stochastic nature of the reservoir. I have however not a clue on how to start.

But the paper talks about Molecular Dynamics. Am I thinking about it in a wrong way?

My previous experience with computational physics is small.

I've only worked with a driven pendulum using the GSL library and the Ising model using the metropolis algorithm.

Sincerely,

Joris

**Physics Forums | Science Articles, Homework Help, Discussion**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Comp. Physics: Finite time Carnot cycle

**Physics Forums | Science Articles, Homework Help, Discussion**