- #1
JorisL
- 492
- 189
Hi,
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.
Joris
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.
Joris
Last edited: