dada-engine.org - User contributions [en]

Thermodynamic and Mechanical Study/fr

2026-09-07T14:59:58Z

FuzzyBot: Updating to match new version of source page

<span id="1._Scope_and_method"></span>
== 1. Périmètre et méthode ==

<div class="mw-translate-fuzzy">
Cette étude décrit le cycle récepteur complet de la machine Dada sous forme symbolique. Le modèle couple une cinématique imposée des deux pistons, des volumes de gaz bien mélangés, des échangeurs thermiques 0D, des clapets passifs et une fermeture hydraulique compressible générique.
</div>

L'origine du cycle est fixée au volume maximal du grand cylindre, avec :

<math>t=0,\qquad \theta(0)=0,\qquad V_L(0)=V_{L,\max}.</math>

Les quatre volumes de gaz sont :

<div class="mw-translate-fuzzy">
* <code>S</code> : petit cylindre, directement associé à l'échangeur froid <code>C</code> ;
* <code>L</code> : grand cylindre, directement associé à l'échangeur chaud <code>H</code> ;
* <code>C</code> : échangeur thermique froid, qui prélève de la chaleur au réservoir froid ;
* <code>H</code> : échangeur thermique chaud, qui rejette de la chaleur vers le réservoir chaud.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The hydraulic topology is fixed:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>S\to H_i\to L,</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>L\to H_o\to S.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The passive check valves therefore always allow:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>H_i\to L,\qquad H_o\to S.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The physical function of each heat exchanger is independent of the operating mode. What changes between refrigeration and motor operation is the external thermal reservoir connected to each exchanger.
The cycle comprises four hydraulic/thermodynamic phases:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
* '''Phase I''' — check valves closed, nominally adiabatic: compression on the <math display="inline">L+H_o</math> side;
* '''Phase II''' — heat exchange, nominally isothermal: transfer <math display="inline">L \to H_o \to S</math>, with heat removed from the gas;
* '''Phase III''' — check valves closed, nominally adiabatic: expansion on the <math display="inline">S+H_i</math> side;
* '''Phase IV''' — heat exchange, nominally isothermal: transfer <math display="inline">S \to H_i \to L</math>, with heat supplied to the gas.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
These phase definitions describe the permanent physical topology of the machine. The external reservoir associated with each heat exchanger depends on the operating mode. The phase boundaries and their angular positions are determined independently in each operating mode by the actual check-valve events.
</div>

Les qualificatifs ''nominalement adiabatique'' et ''nominalement isotherme'' décrivent l'objectif du cycle. Ils ne constituent pas des contraintes thermodynamiques exactes : les températures sont calculées par les bilans, les échangeurs restent couplés à leurs réservoirs, et les zones de faible déplacement des pistons ne sont pas supposées parfaitement immobiles.

<span id="1.1_First-level_assumptions"></span>
=== 1.1 Hypothèses de premier niveau ===

Le modèle repose sur les hypothèses suivantes :

* gaz monophasique, parfait et calorifiquement parfait ;
* propriétés constantes <math display="inline">R</math>, <math display="inline">C_p</math>, <math display="inline">C_v</math>, <math display="inline">\gamma</math> avec <math display="inline">R=C_p-C_v</math> et <math display="inline">\gamma=C_p/C_v</math> ;
* chaque volume de gaz est uniforme et bien mélangé ;
* énergie cinétique et énergie potentielle du gaz négligées dans les bilans 0D ;
* volumes des échangeurs fixes ;
* cinématique des pistons imposée ;
* pas de frottement mécanique dans le modèle thermodynamique ;
* échanges thermiques représentés par une conductance globale <math display="inline">UA</math> ;
* clapets passifs commandés par la différence de pression.

<span id="2._Notation_and_conventions"></span>
== 2. Notations et conventions ==

<span id="2.1_Geometry_and_kinematics"></span>
=== 2.1 Géométrie et cinématique ===

<div class="mw-translate-fuzzy">
Pour <math display="inline">i\in\{S,L\}</math> :
</div>

<div class="mw-translate-fuzzy">
<math>V_{i,\min}>0,\qquad V_{i,\max}>V_{i,\min}.</math>
</div>

Le volume balayé est :

<div class="mw-translate-fuzzy">
<math>V_{i,\mathrm{swept}}=V_{i,\max}-V_{i,\min}.</math>
</div>

Les volumes thermodynamiques imposés par la mécanique sont :

<math>V_S=V_S(t),\qquad V_L=V_L(t),</math>

avec leurs dérivées '''signées''' :

<math>\dot V_S=\frac{dV_S}{dt},\qquad \dot V_L=\frac{dV_L}{dt}.</math>

<div lang="en" dir="ltr" class="mw-content-ltr">
The crank angular velocity is signed:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>\theta(t)=\omega t,\qquad
\dot V_k=\omega\frac{dV_k}{d\theta},
\qquad k\in\{S,L\}.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The driven refrigeration direction is chosen as positive:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>\omega>0.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
Motor operation uses the opposite crank direction:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>\omega<0.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The geometric origin is identical in both modes:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>t=0,\qquad \theta(0)=0,\qquad V_L(0)=V_{L,\max}.</math>
</div>

Lorsque seule la valeur absolue de la vitesse volumique est utile :

<div class="mw-translate-fuzzy">
<math>\nu_{V,i}=|\dot V_i|.</math>
</div>

Une '''zone quasi immobile''' désigne un intervalle où le déplacement ou <math display="inline">|\dot V|</math> reste faible par rapport aux phases de transfert. Cette zone correspond au « plateau » de l'optimisation cinématique, sans supposer <math display="inline">\dot V=0</math> exactement.

<span id="2.2_Kinematic_closure_fraction"></span>
=== 2.2 Fraction de fermeture cinématique ===

Pour décrire la fermeture normalisée d'un cylindre :

<div class="mw-translate-fuzzy">
<math>\Lambda_i(t)=\frac{V_{i,\max}-V_i(t)}{V_{i,\max}-V_{i,\min}},\qquad i\in\{S,L\}.</math>
</div>

<div class="mw-translate-fuzzy">
Ainsi <math display="inline">\Lambda_i=0</math> correspond au volume maximal et <math display="inline">\Lambda_i=1</math> au volume minimal. Les valeurs <math display="inline">\Lambda_L^*</math> et <math display="inline">\Lambda_S^*</math> sont des cibles cinématiques nominales aux transitions. La valeur réellement atteinte lors d'un événement de clapet est :
</div>

<math>\Lambda_{\mathrm{real}}=\Lambda(t_{\mathrm{event}}).</math>

Les transitions physiques restent déterminées par les pressions ; <math display="inline">\Lambda^*</math> n'est donc pas une condition d'ouverture imposée.

<span id="2.3_Thermodynamic_variables"></span>
=== 2.3 Variables thermodynamiques ===

<div class="mw-translate-fuzzy">
Pour chaque volume <math display="inline">i\in\{S,L,C,H\}</math> :
</div>

<div class="mw-translate-fuzzy">
<math>m_i,\qquad U_i,\qquad T_i,\qquad P_i,\qquad V_i.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The subscripts <math display="inline">i</math> and <math display="inline">o</math> denote respectively the gas contained in <math display="inline">H_i</math> and <math display="inline">H_o</math>.
</div>

Le vecteur d'état complet est choisi sous la forme :

<div class="mw-translate-fuzzy">
<math>\mathbf X=(m_S,U_S,m_L,U_L,m_C,U_C,m_H,U_H).</math>
</div>

<div class="mw-translate-fuzzy">
Les températures et pressions en sont déduites :
</div>

<div class="mw-translate-fuzzy">
<math>T_i=\frac{U_i}{m_iC_v},\qquad P_i=\frac{m_iRT_i}{V_i}.</math>
</div>

<div class="mw-translate-fuzzy">
Les volumes <math display="inline">V_S(t)</math> et <math display="inline">V_L(t)</math> sont imposés par la cinématique et ne sont pas des variables d'état thermodynamiques indépendantes. Les volumes <math display="inline">V_C</math> et <math display="inline">V_H</math> sont constants.
</div>

<span id="2.4_Energy_sign_convention"></span>
=== 2.4 Convention de signe énergétique ===

La chaleur est positive lorsqu'elle est '''reçue par le gaz'''. Le travail est positif lorsque le '''gaz fournit du travail''' :

<math>\dot W=P\dot V.</math>

<div lang="en" dir="ltr" class="mw-content-ltr">
By definition of the two physical heat exchangers:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>Q_i>0,\qquad Q_o<0</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
in the intended operating regime of both refrigeration and motor operation.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The net cycle work distinguishes the two modes:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>W_{\mathrm{cycle}}<0</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
for driven refrigeration operation, whereas:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>W_{\mathrm{cycle}}>0</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
for motor operation.
</div>

<span id="3._Thermal_closure_and_validity_domain"></span>
== 3. Fermeture thermique et domaine de validité ==

<span id="3.1_Exchange_with_the_thermal_reservoirs"></span>
=== 3.1 Échange avec les réservoirs thermiques ===

<div lang="en" dir="ltr" class="mw-content-ltr">
For the heat-in exchanger:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>\dot Q_i=(UA)_i(T_{i,\mathrm{res}}-T_i).</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
In the intended operating regime:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>T_i<T_{i,\mathrm{res}}
\quad\Rightarrow\quad
\dot Q_i>0.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
For the heat-out exchanger:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>\dot Q_o=(UA)_o(T_{o,\mathrm{res}}-T_o).</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
In the intended operating regime:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>T_o>T_{o,\mathrm{res}}
\quad\Rightarrow\quad
\dot Q_o<0.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The reservoir temperatures depend on the operating mode.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
For refrigeration operation:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>T_{i,\mathrm{res}}=T_{\mathrm{cold}},\qquad
T_{o,\mathrm{res}}=T_{\mathrm{hot}}.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
For motor operation:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>T_{i,\mathrm{res}}=T_{\mathrm{hot}},\qquad
T_{o,\mathrm{res}}=T_{\mathrm{cold}}.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
Thus the heat-transfer equations themselves are identical in both modes.
</div>

Une phase nominalement isotherme ne signifie donc pas <math display="inline">T_{HX}=T_{HX,\mathrm{res}}</math> : un écart de température fini est nécessaire pour transférer une puissance thermique finie lorsque <math display="inline">UA</math> est fini.

Un indicateur de qualité d'isothermie peut être défini sur une phase donnée par :

<math>\varepsilon_T=\frac{T_{\max}-T_{\min}}{T_{\mathrm{ref}}}.</math>

<span id="3.2_Thermophysical_validity_domain"></span>
=== 3.2 Domaine de validité thermophysique ===

Le modèle de base suppose :

<math>PV=mRT,\qquad Z=1,</math>

<math>C_p=\mathrm{const},\qquad C_v=\mathrm{const},\qquad \gamma=\mathrm{const}.</math>

La validité doit être contrôlée a posteriori sur tout le cycle, notamment par :

<math>|Z-1|\ll1,</math>

et par une faible variation des propriétés thermophysiques, par exemple :

<math>\varepsilon_{C_p}=\frac{C_{p,\max}-C_{p,\min}}{C_{p,\mathrm{ref}}}\ll1.</math>

Le fluide doit rester monophasique gazeux et suffisamment éloigné de toute condensation ou transition de phase dans l'ensemble du domaine <math display="inline">(P,T)</math> parcouru.

Si ces critères deviennent insuffisants, une extension pourra employer <math display="inline">Z(P,T)</math>, <math display="inline">C_p(T)</math>, <math display="inline">C_v(T)</math> ou une équation d'état réelle sans modifier l'architecture générale des bilans de masse et d'énergie.

<span id="4._Reduced_formulation_of_a_quasi-pressure-equalized_pair"></span>
== 4. Formulation réduite d'une paire quasi égalisée en pression ==

Lorsqu'un cylindre et son échangeur sont reliés par une liaison interne très peu résistive, on peut utiliser l'approximation :

<math>P_{\mathrm{cyl}}\approx P_{HX}=P_{\mathrm{pair}}.</math>

Elle est acceptable si :

<math>\varepsilon_P=\frac{|\Delta P_{\mathrm{int}}|}{P_{\mathrm{pair}}}\ll1,\qquad \Delta P_{\mathrm{int}}=P_{\mathrm{cyl}}-P_{HX}.</math>

Un faible nombre de Mach interne constitue un contrôle complémentaire :

<math>Ma_{\mathrm{int}}=\frac{|u_{\mathrm{int}}|}{a}\ll1,</math>

mais il ne suffit pas, à lui seul, à garantir l'égalisation des pressions.

<span id="4.1_Closed-pair_case"></span>
=== 4.1 Cas fermé de la paire ===

Pour une paire fermée, avec <math display="inline">V_{HX}</math> fixe et <math display="inline">V=V_{\mathrm{cyl}}+V_{HX}</math> :

<math>
\frac{dP}{dt}=
\frac{(\gamma-1)(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX})-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

Lorsque <math display="inline">(UA)_{HX}=0</math> :

<math>P(V_{\mathrm{cyl}}+V_{HX})^\gamma=\mathrm{const}.</math>

<span id="4.2_Internal_redistribution_flow_rate"></span>
=== 4.2 Débit interne de redistribution ===

<div class="mw-translate-fuzzy">
On définit positif le débit interne cylindre → échangeur. Posons :
</div>

<math>N=V_{HX}\dot P+(\gamma-1)(UA)_{HX}(T_{HX}-T_{HX,\mathrm{res}}).</math>

Le débit transporte l'enthalpie de l'état amont :

<math>
\dot m_{\mathrm{int}}=
\begin{cases}
\dfrac{N}{\gamma R T_{\mathrm{cyl}}}, & N\ge0 \quad (\mathrm{cyl}\to HX),\\[6pt]
\dfrac{N}{\gamma R T_{HX}}, & N<0 \quad (HX\to\mathrm{cyl}).
\end{cases}
</math>

La température de l'échangeur évolue selon :

<math>
\dot T_{HX}=\frac{T_{HX}}{P}\dot P-
\frac{RT_{HX}^2}{PV_{HX}}\dot m_{\mathrm{int}}
.</math>

<span id="4.3_Open_receiving_pair"></span>
=== 4.3 Paire réceptrice ouverte ===

Lors d'une phase active, le débit externe entre physiquement dans le '''cylindre récepteur''', et non directement dans son échangeur associé. Si <math display="inline">\dot m_{\mathrm{ext}}>0</math> entre dans le cylindre à la température <math display="inline">T_{\mathrm{ext}}</math> :

<math>
\dot P=
\frac{\gamma RT_{\mathrm{ext}}\dot m_{\mathrm{ext}}
+(\gamma-1)(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX})
-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

Les masses vérifient :

<math>\dot m_{\mathrm{cyl}}=\dot m_{\mathrm{ext}}-\dot m_{\mathrm{int}},\qquad
\dot m_{HX}=\dot m_{\mathrm{int}},</math>

et donc :

<math>\dot m_{\mathrm{pair}}=\dot m_{\mathrm{ext}}.</math>

Le cas fermé est obtenu immédiatement avec <math display="inline">\dot m_{\mathrm{ext}}=0</math>.

<span id="5._Hydraulic_closure"></span>
== 5. Fermeture hydraulique ==

<span id="5.1_Generic_formulation"></span>
=== 5.1 Formulation générique ===

Toute liaison hydraulique est décrite par une loi générique :

<math>\dot m=\Phi(P_u,P_d,T_u,\mathcal G,\mathcal F)</math>

<math display="inline">u</math> et <math display="inline">d</math> désignent respectivement l'amont (''upstream'') et l'aval (''downstream''). L'enthalpie transportée est celle de l'amont :

<math>\dot H_{\mathrm{mass}}=\dot m C_pT_u.</math>

Pour une liaison bidirectionnelle, l'état amont est déterminé par le sens réel du gradient de pression. Pour un clapet, le débit inverse est interdit.

<span id="5.2_First-level_closure_using_a_compressible_orifice"></span>
=== 5.2 Fermeture de premier niveau par orifice compressible ===

Une première approximation consiste à employer une aire hydraulique effective :

<math>(C_dA)_{\mathrm{eff}},</math>

qui représente la facilité globale de passage du gaz dans la liaison réelle.

Avec :

<math>r=\frac{P_d}{P_u},\qquad
r_{\mathrm{crit}}=\left(\frac{2}{\gamma+1}\right)^{\gamma/(\gamma-1)},</math>

le débit non étranglé, pour <math display="inline">r>r_{\mathrm{crit}}</math>, est :

<math>
\dot m=(C_dA)_{\mathrm{eff}}P_u
\sqrt{\frac{2\gamma}{RT_u(\gamma-1)}
\left(r^{2/\gamma}-r^{(\gamma+1)/\gamma}\right)}
.</math>

Pour <math display="inline">r\le r_{\mathrm{crit}}</math> :

<math>
\dot m=(C_dA)_{\mathrm{eff}}P_u
\sqrt{\frac{\gamma}{RT_u}}
\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}}
.</math>

Cette fermeture n'est pas l'essence du modèle : elle pourra être remplacée par une loi plus fidèle aux pertes de charge d'un échangeur réel, d'une conduite ou d'un clapet.

<span id="6._Complete_thermodynamic_cycle"></span>
== 6. Cycle thermodynamique complet ==

<span id="6.1_Phase_I_—_check_valves_closed,_nominally_adiabatic:_compression_on_the_L_side"></span>
=== 6.1 Phase I — clapets fermés, nominalement adiabatique : compression côté L ===

<div class="mw-translate-fuzzy">

Les deux clapets sont fermés. Les paires <code>L+H</code> et <code>S+C</code> sont fermées. Le côté <code>L+H</code> est nominalement comprimé ; le mouvement du petit piston reste celui fourni par la cinématique réelle.
</div>

<div class="mw-translate-fuzzy">
Pour <code>L+H</code> :
</div>

<div class="mw-translate-fuzzy">
<math>\dot P_{LH}=
\frac{(\gamma-1)(UA)_H(T_{H,\mathrm{res}}-T_H)-\gamma P_{LH}\dot V_L}
{V_L+V_H}.</math>
</div>

<div class="mw-translate-fuzzy">
Pour <code>S+C</code> :
</div>

<div class="mw-translate-fuzzy">
<math>\dot P_{SC}=
\frac{(\gamma-1)(UA)_C(T_{C,\mathrm{res}}-T_C)-\gamma P_{SC}\dot V_S}
{V_S+V_C}.</math>
</div>

Les équations de redistribution interne et de température du §4 s'appliquent aux deux paires.

<div class="mw-translate-fuzzy">
La transition vers la Phase II survient lorsque le clapet <code>H → S</code> satisfait sa condition d'ouverture.
</div>

<span id="6.2_Phase_II_—_heat_exchange,_nominally_isothermal:_L_→_Ho_→_S"></span>
<div class="mw-translate-fuzzy">
=== 6.2 Phase II — échange thermique, nominalement isotherme : L → H → S ===
</div>

<div class="mw-translate-fuzzy">
Le gaz quitte <code>L</code>, traverse l'échangeur chaud <code>H</code>, y rejette de la chaleur, franchit le clapet <code>H → S</code>, puis entre dans le cylindre récepteur <code>S</code>. La paire <code>S+C</code> reste quasi égalisée si le critère <math display="inline">\varepsilon_P\ll1</math> est satisfait.
</div>

<span id="6.2.1_Donor_cylinder_L"></span>
==== 6.2.1 Cylindre donneur L ====

Le bilan fondamental est :

<math>
\frac{d(m_LC_vT_L)}{dt}
=-P_L\dot V_L-\dot m_{L,\mathrm{out}}C_pT_L
.</math>

La cinématique est conçue pour maintenir un débit sortant du cylindre donneur. Dans ce cas, la solution analytique est :

<math>
\frac{T_L}{T_{L,\mathrm{ref}}}=
\left[\frac{m_L}{m_{L,\mathrm{ref}}}\frac{V_{L,\mathrm{ref}}}{V_L}\right]^{\gamma-1}
</math>

et :

<math>
P_L=P_{L,\mathrm{ref}}
\left[\frac{m_L}{m_{L,\mathrm{ref}}}\frac{V_{L,\mathrm{ref}}}{V_L}\right]^\gamma
.</math>

On obtient également :

<math>\frac{T_L}{T_{L,\mathrm{ref}}}=
\left(\frac{P_L}{P_{L,\mathrm{ref}}}\right)^{(\gamma-1)/\gamma},</math>

et l'entropie spécifique du gaz restant vérifie <math display="inline">ds=0</math> sous ces hypothèses. Si un débit rétrograde apparaît malgré la conception cinématique, la solution analytique n'est plus applicable et il faut utiliser le bilan ouvert complet en <math display="inline">(m_L,U_L)</math>.

<span id="6.2.2_Heat-out_exchanger_Ho"></span>
<div class="mw-translate-fuzzy">
==== 6.2.2 Échangeur chaud H ====
</div>

Conservation de la masse :

<div class="mw-translate-fuzzy">
<math>\dot m_H=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}}.</math>
</div>

Bilan énergétique fondamental :

<div class="mw-translate-fuzzy">
<math>
\frac{d(m_HC_vT_H)}{dt}
=\dot m_{\mathrm{in}}C_pT_L
-\dot m_{\mathrm{out}}C_pT_H
+(UA)_H(T_{H,\mathrm{res}}-T_H)
.</math>
</div>

Sous forme développée :

<div class="mw-translate-fuzzy">
<math>
\dot T_H=
\frac{\dot m_{\mathrm{in}}(C_pT_L-C_vT_H)
-\dot m_{\mathrm{out}}RT_H
+(UA)_H(T_{H,\mathrm{res}}-T_H)}
{C_vm_H}
.</math>
</div>

Les débits sont déterminés par les lois hydrauliques :

<div class="mw-translate-fuzzy">
<math>\dot m_{\mathrm{in}}=\Phi_{\mathrm{HX,H}}(P_L,P_H,T_L,\ldots),</math>
</div>

<div class="mw-translate-fuzzy">
<math>\dot m_{\mathrm{out}}=\Phi_{\mathrm{valve,H}}(P_H,P_{SC},T_H,\ldots).</math>
</div>

<span id="6.2.3_Receiving_pair_S+Hi"></span>
<div class="mw-translate-fuzzy">
==== 6.2.3 Paire réceptrice S+C ====
</div>

<div class="mw-translate-fuzzy">
Le débit issu de <code>H</code> entre dans '''S'''. La pression de la paire vérifie :
</div>

<div class="mw-translate-fuzzy">
<math>
\dot P_{SC}=
\frac{\gamma R\dot m_{\mathrm{out}}T_H
+(\gamma-1)(UA)_C(T_{C,\mathrm{res}}-T_C)
-\gamma P_{SC}\dot V_S}
{V_S+V_C}
.</math>
</div>

<div class="mw-translate-fuzzy">
Les équations de redistribution interne du §4 restent inchangées : le débit externe agit sur l'échangeur <code>C</code> indirectement par l'évolution de la pression de paire.
</div>

<div class="mw-translate-fuzzy">
La transition vers la Phase III est l'événement de fermeture du clapet <code>H → S</code>.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
=== 6.3 Phase III — check valves closed, nominally adiabatic: expansion on the S side ===
</div>

<div class="mw-translate-fuzzy">
Les deux clapets sont fermés. La paire <code>S+C</code> se détend nominalement ; la paire <code>L+H</code> reste également fermée. Aucun piston n'est supposé rigoureusement immobile.
</div>

<div class="mw-translate-fuzzy">
Pour <code>S+C</code> :
</div>

<div class="mw-translate-fuzzy">
<math>\dot P_{SC}=
\frac{(\gamma-1)(UA)_C(T_{C,\mathrm{res}}-T_C)-\gamma P_{SC}\dot V_S}
{V_S+V_C}.</math>
</div>

<div class="mw-translate-fuzzy">
Pour <code>L+H</code> :
</div>

<div class="mw-translate-fuzzy">
<math>\dot P_{LH}=
\frac{(\gamma-1)(UA)_H(T_{H,\mathrm{res}}-T_H)-\gamma P_{LH}\dot V_L}
{V_L+V_H}.</math>
</div>

<div class="mw-translate-fuzzy">
La transition vers la Phase IV survient lorsque le clapet <code>C → L</code> satisfait sa condition d'ouverture.
</div>

<span id="6.4_Phase_IV_—_heat_exchange,_nominally_isothermal:_S_→_Hi_→_L"></span>
<div class="mw-translate-fuzzy">
=== 6.4 Phase IV — échange thermique, nominalement isotherme : S → C → L ===
</div>

<div class="mw-translate-fuzzy">
Le gaz quitte <code>S</code>, traverse l'échangeur froid <code>C</code>, y reçoit de la chaleur du réservoir froid, franchit le clapet <code>C → L</code>, puis entre dans le cylindre récepteur <code>L</code>. La paire <code>L+H</code> reste quasi égalisée si le critère <math display="inline">\varepsilon_P\ll1</math> est satisfait.
</div>

<span id="6.4.1_Donor_cylinder_S"></span>
==== 6.4.1 Cylindre donneur S ====

Le bilan fondamental est :

<math>
\frac{d(m_SC_vT_S)}{dt}
=-P_S\dot V_S-\dot m_{S,\mathrm{out}}C_pT_S
.</math>

Sous débit sortant garanti par la conception cinématique :

<math>
\frac{T_S}{T_{S,\mathrm{ref}}}=
\left[\frac{m_S}{m_{S,\mathrm{ref}}}\frac{V_{S,\mathrm{ref}}}{V_S}\right]^{\gamma-1}
</math>

et :

<math>
P_S=P_{S,\mathrm{ref}}
\left[\frac{m_S}{m_{S,\mathrm{ref}}}\frac{V_{S,\mathrm{ref}}}{V_S}\right]^\gamma
.</math>

Un débit rétrograde impose de revenir au bilan ouvert complet en <math display="inline">(m_S,U_S)</math>.

<div lang="en" dir="ltr" class="mw-content-ltr">
==== 6.4.2 Heat-in exchanger Hi ====
</div>

<div class="mw-translate-fuzzy">
<math>\dot m_C=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}}.</math>
</div>

Le bilan énergétique fondamental est :

<div class="mw-translate-fuzzy">
<math>
\frac{d(m_CC_vT_C)}{dt}
=\dot m_{\mathrm{in}}C_pT_S
-\dot m_{\mathrm{out}}C_pT_C
+(UA)_C(T_{C,\mathrm{res}}-T_C)
.</math>
</div>

Sous forme développée :

<div class="mw-translate-fuzzy">
<math>
\dot T_C=
\frac{\dot m_{\mathrm{in}}(C_pT_S-C_vT_C)
-\dot m_{\mathrm{out}}RT_C
+(UA)_C(T_{C,\mathrm{res}}-T_C)}
{C_vm_C}
.</math>
</div>

Les débits sont déterminés par :

<div class="mw-translate-fuzzy">
<math>\dot m_{\mathrm{in}}=\Phi_{\mathrm{HX,C}}(P_S,P_C,T_S,\ldots),</math>
</div>

<div class="mw-translate-fuzzy">
<math>\dot m_{\mathrm{out}}=\Phi_{\mathrm{valve,C}}(P_C,P_{LH},T_C,\ldots).</math>
</div>

<span id="6.4.3_Receiving_pair_L+Ho"></span>
<div class="mw-translate-fuzzy">
==== 6.4.3 Paire réceptrice L+H ====
</div>

<div class="mw-translate-fuzzy">
Le débit issu de <code>C</code> entre dans '''L'''. La pression de la paire vérifie :
</div>

<div class="mw-translate-fuzzy">
<math>
\dot P_{LH}=
\frac{\gamma R\dot m_{\mathrm{out}}T_C
+(\gamma-1)(UA)_H(T_{H,\mathrm{res}}-T_H)
-\gamma P_{LH}\dot V_L}
{V_L+V_H}
.</math>
</div>

<div class="mw-translate-fuzzy">
La fermeture du clapet <code>C → L</code> ramène à la Phase I. La fermeture géométrique des pistons ne suffit pas à garantir la fermeture thermodynamique du cycle.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
=== 6.5 Motor operation ===
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
No second set of thermodynamic balance equations is required for motor operation.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The physical hydraulic topology remains:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>L\to H_o\to S,\qquad
S\to H_i\to L.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The check-valve directions remain:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>H_o\to S,\qquad H_i\to L.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
Motor operation is obtained by:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
# reversing the crank direction,
# exchanging the external reservoirs connected to <math display="inline">H_i</math> and <math display="inline">H_o</math>.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
Thus:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>\omega<0,</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
with the same cycle origin:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>V_L(0)=V_{L,\max}.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The thermal-reservoir assignment becomes:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>T_{i,\mathrm{res}}=T_{\mathrm{hot}},\qquad
T_{o,\mathrm{res}}=T_{\mathrm{cold}}.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The heat-in exchanger therefore absorbs heat from the hot reservoir:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>Q_i>0,</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
while the heat-out exchanger rejects heat to the cold reservoir:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>Q_o<0.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The pressure histories, mass flow rates, check-valve events, and periodic thermodynamic state must be recalculated with the reversed kinematics.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The motor regime is obtained when:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>W_{\mathrm{cycle}}>0.</math>
</div>

<span id="7._Physical_transitions_of_the_check_valves"></span>
== 7. Transitions physiques des clapets ==

Pour un clapet orienté de l'amont <math display="inline">u</math> vers l'aval <math display="inline">d</math> :

<math>P_u-P_d\ge\Delta P_{\mathrm{open}}
\quad\Rightarrow\quad \text{opening},</math>

<math>P_u-P_d\le\Delta P_{\mathrm{close}}
\quad\Rightarrow\quad \text{closing},</math>

avec, en présence d'hystérésis :

<math>\Delta P_{\mathrm{close}}\le\Delta P_{\mathrm{open}}.</math>

Une ouverture ou une fermeture modifie la topologie hydraulique et donc les équations actives ; elle ne provoque aucun saut instantané de l'état thermodynamique. Pour chaque volume :

<div class="mw-translate-fuzzy">
<math>m_i^+=m_i^-,\qquad U_i^+=U_i^-,\qquad V_i^+=V_i^-.</math>
</div>

Pour le gaz parfait :

<div class="mw-translate-fuzzy">
<math>T_i^+=T_i^-,\qquad P_i^+=P_i^-.</math>
</div>

Il n'y a donc pas d'égalisation instantanée des pressions à l'ouverture d'un clapet.

Les valeurs <math display="inline">\Lambda^*</math> servent de cibles de conception cinématique ; les valeurs réellement observées lors des transitions sont <math display="inline">\Lambda_{\mathrm{real}}=\Lambda(t_{\mathrm{event}})</math>.

<div lang="en" dir="ltr" class="mw-content-ltr">
The two passive check valves have permanent physical orientations:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>H_o\to S,\qquad H_i\to L.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
Their opening and closing conditions are identical in both operating modes. Only the pressure histories change because the crank kinematics are reversed in motor operation.
</div>

<span id="8._Work,_heat,_and_performance"></span>
<div class="mw-translate-fuzzy">
== 8. Travail, chaleur et coefficient de performance ==
</div>

Le travail instantané du gaz sur les deux pistons est calculé pendant toutes les phases :

<math>
\dot W=P_S^\star\dot V_S+P_L^\star\dot V_L
.</math>

<math display="inline">P_S^\star</math> et <math display="inline">P_L^\star</math> désignent la pression thermodynamique effectivement appliquée au gaz dans chaque cylindre selon la topologie de la phase : pression de paire lorsque le cylindre appartient à une paire quasi égalisée, pression propre lorsqu'il est donneur isolé hydrauliquement.

Le travail net du cycle est :

<math>
W_{\mathrm{cycle}}=
\int_0^\tau
\left(P_S^\star\dot V_S+P_L^\star\dot V_L\right)dt
.</math>

Les chaleurs échangées sont :

<div class="mw-translate-fuzzy">
<math>Q_C=\int_0^\tau (UA)_C(T_{C,\mathrm{res}}-T_C)\,dt,</math>
</div>

<div class="mw-translate-fuzzy">
<math>Q_H=\int_0^\tau (UA)_H(T_{H,\mathrm{res}}-T_H)\,dt.</math>
</div>

En régime périodique :

<math>\Delta U_{\mathrm{cycle}}=0,</math>

et le premier principe donne :

<div class="mw-translate-fuzzy">
<math>Q_C+Q_H=W_{\mathrm{cycle}}.</math>
</div>

Le COP frigorifique est :

<div class="mw-translate-fuzzy">
<math>COP_c=\frac{Q_C}{-W_{\mathrm{cycle}}}.</math>
</div>

Le COP en pompe à chaleur est :

<div class="mw-translate-fuzzy">
<math>COP_h=\frac{-Q_H}{-W_{\mathrm{cycle}}}=COP_c+1.</math>
</div>

<div class="mw-translate-fuzzy">
Les signes <math display="inline">Q_C>0</math>, <math display="inline">Q_H<0</math> et <math display="inline">W_{\mathrm{cycle}}<0</math> constituent des contrôles du régime frigorifique recherché.
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
For motor operation:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>Q_i>0,\qquad Q_o<0,\qquad W_{\mathrm{cycle}}>0.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The thermal efficiency is:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>
\eta_{\mathrm{th}}
=
\frac{W_{\mathrm{cycle}}}{Q_i}
=
1+\frac{Q_o}{Q_i}.
</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
The mean thermodynamic motor power is:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>
\overline{\dot W}
=
\frac{W_{\mathrm{cycle}}}{\tau}.
</math>
</div>

La force thermodynamique exercée par le gaz sur une face de piston peut être écrite <math display="inline">F_{\mathrm{gas}}=PS</math>. La force mécanique nette, les inerties et les frottements relèvent du dimensionnement mécanique ultérieur.

<span id="9._Gas_charge_and_periodic_regime"></span>
== 9. Charge de gaz et régime périodique ==

La quantité totale de gaz enfermée est un paramètre physique :

<div class="mw-translate-fuzzy">
<math>M_{\mathrm{tot}}=m_S+m_L+m_C+m_H=\mathrm{const}.</math>
</div>

Elle peut être imposée directement ou définie expérimentalement par une pression et une température de charge. La configuration de référence choisie est celle de <math display="inline">t=0</math>, avec le grand cylindre au volume maximal. Si tous les volumes communiquent et sont en équilibre uniforme pendant le remplissage :

<div class="mw-translate-fuzzy">
<math>
M_{\mathrm{tot}}=
\frac{P_{\mathrm{charge}}
\left[V_S(0)+V_{L,\max}+V_C+V_H\right]}
{RT_{\mathrm{charge}}}
.</math>
</div>

<math display="inline">P_{\mathrm{charge}}</math> et <math display="inline">T_{\mathrm{charge}}</math> définissent la quantité de gaz chargée ; ils ne constituent pas une condition que le cycle périodique doit retrouver.

Le régime périodique établi est une solution du système telle que, entre deux passages successifs au volume maximal du grand cylindre avec la même direction cinématique :

<math>\mathbf X(t+\tau)=\mathbf X(t).</math>

La périodicité géométrique seule :

<div class="mw-translate-fuzzy">
<math>V_i(t+\tau)=V_i(t)</math>
</div>

ne suffit pas à garantir la périodicité thermodynamique.

L'état numérique utilisé pour lancer un calcul peut être approximatif ; il ne doit pas être confondu avec un paramètre physique de la machine. Le futur solveur pourra rechercher le point fixe périodique par cycles successifs, méthode de tir ou méthode de Newton.

<span id="10._Parameters,_design_data,_and_results"></span>
== 10. Paramètres, données de conception et résultats ==

<span id="10.1_Prescribed_data"></span>
=== 10.1 Données imposées ===

<div class="mw-translate-fuzzy">
* fluide de travail et propriétés de référence <math display="inline">R</math>, <math display="inline">C_p</math>, <math display="inline">C_v</math>, <math display="inline">\gamma</math> ;
* températures des réservoirs <math display="inline">T_{C,\mathrm{res}}</math>, <math display="inline">T_{H,\mathrm{res}}</math> ;
* charge totale <math display="inline">M_{\mathrm{tot}}</math>, ou équivalent <math display="inline">(P_{\mathrm{charge}},T_{\mathrm{charge}})</math> dans la configuration de remplissage définie au §9 ;
* cinématique <math display="inline">V_S(t)</math>, <math display="inline">V_L(t)</math> et, si pertinent, <math display="inline">\omega</math>.
</div>

<span id="10.2_Design_parameters"></span>
=== 10.2 Paramètres de conception ===

<div class="mw-translate-fuzzy">
* <math display="inline">V_{S,\min}</math>, <math display="inline">V_{S,\max}</math>, <math display="inline">V_{L,\min}</math>, <math display="inline">V_{L,\max}</math> ;
* <math display="inline">V_C</math>, <math display="inline">V_H</math> ;
* <math display="inline">(UA)_C</math>, <math display="inline">(UA)_H</math> ;
* paramètres hydrauliques des échangeurs, conduites et clapets, représentés au premier niveau par <math display="inline">(C_dA)_{\mathrm{eff}}</math> ;
* seuils <math display="inline">\Delta P_{\mathrm{open}}</math>, <math display="inline">\Delta P_{\mathrm{close}}</math> ;
* cibles cinématiques <math display="inline">\Lambda_L^*</math>, <math display="inline">\Lambda_S^*</math>.
</div>

<span id="10.3_Calculated_variables_and_results"></span>
=== 10.3 Variables et résultats calculés ===

<div class="mw-translate-fuzzy">
* <math display="inline">m_i</math>, <math display="inline">U_i</math>, <math display="inline">T_i</math>, <math display="inline">P_i</math> ;
* débits massiques internes et externes ;
* <math display="inline">\dot Q_C</math>, <math display="inline">\dot Q_H</math>, <math display="inline">Q_C</math>, <math display="inline">Q_H</math> ;
* <math display="inline">W_{\mathrm{cycle}}</math>, <math display="inline">COP_c</math>, <math display="inline">COP_h</math> ;
* extrema de pression, température et débit ;
* événements réels des clapets et <math display="inline">\Lambda_{\mathrm{real}}</math> ;
* qualité d'isothermie <math display="inline">\varepsilon_T</math> ;
* critères de validité <math display="inline">\varepsilon_P</math>, <math display="inline">Ma</math>, <math display="inline">Z</math> et variation des propriétés.
</div>

<span id="11._Global_conservation_checks"></span>
== 11. Contrôles globaux de conservation ==

<span id="11.1_Mass_conservation"></span>
=== 11.1 Conservation de la masse ===

Le solveur doit satisfaire :

<math>\frac{dM_{\mathrm{tot}}}{dt}=0.</math>

Un résidu numérique utile est :

<math>\varepsilon_M(t)=M_{\mathrm{tot}}(t)-M_{\mathrm{tot}}(0).</math>

<span id="11.2_Global_energy_conservation"></span>
=== 11.2 Conservation globale de l'énergie ===

Quelle que soit la phase, les flux massiques et enthalpiques internes doivent s'annuler lorsque les bilans de tous les volumes sont additionnés. Le bilan global doit se réduire à :

<div class="mw-translate-fuzzy">
<math>
\frac{dU_{\mathrm{tot}}}{dt}
=\dot Q_C+\dot Q_H
-P_S^\star\dot V_S
-P_L^\star\dot V_L
.</math>
</div>

Un résidu énergétique cumulatif peut être défini par :

<div class="mw-translate-fuzzy">
<math>
\varepsilon_E(t)=
U_{\mathrm{tot}}(t)-U_{\mathrm{tot}}(0)
-Q_C(0,t)-Q_H(0,t)+W(0,t)
.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
Over a periodic cycle:
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>Q_i+Q_o=W_{\mathrm{cycle}}.</math>
</div>

Le solveur doit maintenir <math display="inline">\varepsilon_M</math> et <math display="inline">\varepsilon_E</math> proches de zéro à la précision numérique attendue.

----

<span id="Appendix_A_—_Symbolic_derivations_and_validated_checks"></span>
<div class="mw-translate-fuzzy">
= Annexe A — Dérivations symboliques et contrôles validés =
</div>

<span id="A.1_Pressure_equation_for_a_closed_pair"></span>
<div class="mw-translate-fuzzy">
== A.1 Équation de pression d'une paire fermée ==
</div>

Pour une paire cylindre + échangeur à pression quasi uniforme :

<math>U=\frac{P(V_{\mathrm{cyl}}+V_{HX})}{\gamma-1}.</math>

Le premier principe donne :

<math>\frac{dU}{dt}=\dot Q-P\dot V_{\mathrm{cyl}},</math>

avec :

<math>\dot Q=(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX}).</math>

En dérivant <math display="inline">U</math> :

<math>\frac{1}{\gamma-1}
\left[(V_{\mathrm{cyl}}+V_{HX})\dot P+P\dot V_{\mathrm{cyl}}\right]
=\dot Q-P\dot V_{\mathrm{cyl}}.</math>

D'où :

<math>
\dot P=
\frac{(\gamma-1)\dot Q-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

Si <math display="inline">\dot Q=0</math> :

<math>\frac{\dot P}{P}=-\gamma\frac{\dot V}{V},</math>

puis :

<math>PV^\gamma=\mathrm{const}.</math>

<span id="A.2_Internal_flow_rate_of_the_pair"></span>
<div class="mw-translate-fuzzy">

== A.2 Débit interne de la paire ==
</div>

Pour l'échangeur seul, volume fixe :

<math>U_{HX}=\frac{PV_{HX}}{\gamma-1}.</math>

Donc :

<math>\frac{V_{HX}}{\gamma-1}\dot P
=(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX})
+\dot m_{\mathrm{int}}C_pT_{\mathrm{up}}.</math>

En utilisant <math display="inline">C_p=\gamma R/(\gamma-1)</math> :

<math>V_{HX}\dot P+(\gamma-1)(UA)_{HX}(T_{HX}-T_{HX,\mathrm{res}})
=\gamma R T_{\mathrm{up}}\dot m_{\mathrm{int}}.</math>

On retrouve la définition du numérateur <math display="inline">N</math> et la sélection de la température amont selon le signe du débit.

<span id="A.3_Evolution_of_the_heat-exchanger_temperature_within_a_pair"></span>
<div class="mw-translate-fuzzy">
== A.3 Évolution de la température de l'échangeur dans une paire ==
</div>

Pour un volume fixe :

<math>m_{HX}=\frac{PV_{HX}}{RT_{HX}}.</math>

En dérivant :

<math>\frac{\dot m_{HX}}{m_{HX}}
=\frac{\dot P}{P}-\frac{\dot T_{HX}}{T_{HX}}.</math>

Avec <math display="inline">\dot m_{HX}=\dot m_{\mathrm{int}}</math> :

<math>
\dot T_{HX}=\frac{T_{HX}}P\dot P-
\frac{RT_{HX}^2}{PV_{HX}}\dot m_{\mathrm{int}}
.</math>

Cette relation découle uniquement de la conservation de masse et de l'équation d'état ; elle est valable pour les deux sens de débit.

<span id="A.4_Open_receiving_pair"></span>
<div class="mw-translate-fuzzy">
== A.4 Paire réceptrice ouverte ==
</div>

Pour une paire recevant <math display="inline">\dot m_{\mathrm{ext}}</math> dans son cylindre :

<math>\frac{dU}{dt}=\dot m_{\mathrm{ext}}C_pT_{\mathrm{ext}}
+\dot Q-P\dot V_{\mathrm{cyl}}.</math>

Avec <math display="inline">U=P(V_{\mathrm{cyl}}+V_{HX})/(\gamma-1)</math> :

<math>
\dot P=
\frac{\gamma RT_{\mathrm{ext}}\dot m_{\mathrm{ext}}
+(\gamma-1)\dot Q
-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

Le bilan de masse local :

<math>\dot m_{\mathrm{cyl}}=\dot m_{\mathrm{ext}}-\dot m_{\mathrm{int}},</math>

<math>\dot m_{HX}=\dot m_{\mathrm{int}},</math>

redonne immédiatement :

<math>\dot m_{\mathrm{pair}}=\dot m_{\mathrm{ext}}.</math>

<span id="A.5_Analytical_solution_for_the_adiabatic_donor_cylinder"></span>
<div class="mw-translate-fuzzy">
== A.5 Solution analytique du cylindre donneur adiabatique ==
</div>

Pour un cylindre adiabatique, bien mélangé, avec débit uniquement sortant :

<math>\frac{d(mC_vT)}{dt}=-P\dot V-\dot m_{\mathrm{out}}C_pT,</math>

et :

<math>\dot m=-\dot m_{\mathrm{out}}.</math>

En développant :

<math>C_vm\dot T+C_vT\dot m=-P\dot V+C_pT\dot m,</math>

soit :

<math>C_vm\dot T=-P\dot V+RT\dot m.</math>

Avec <math display="inline">P=mRT/V</math> et <math display="inline">R/C_v=\gamma-1</math> :

<math>\frac{dT}{T}=(\gamma-1)
\left(\frac{dm}{m}-\frac{dV}{V}\right).</math>

Après intégration :

<math>
\frac{T}{T_0}=
\left[\frac{m}{m_0}\frac{V_0}{V}\right]^{\gamma-1}
.</math>

Puis, avec <math display="inline">PV=mRT</math> :

<math>
P=P_0\left[\frac{m}{m_0}\frac{V_0}{V}\right]^\gamma
.</math>

et :

<math>
\frac{T}{T_0}=
\left(\frac{P}{P_0}\right)^{(\gamma-1)/\gamma}
.</math>

Sous ces hypothèses, l'entropie spécifique du gaz restant est constante : <math display="inline">ds=0</math>. L'entropie totale du contenu du cylindre n'est pas constante puisque sa masse varie.

<span id="A.6_Active_heat_exchanger:_expanded_balance"></span>
<div class="mw-translate-fuzzy">
== A.6 Échangeur actif : bilan développé ==
</div>

Bilan fondamental :

<math>\frac{d(mC_vT)}{dt}
=\dot m_{\mathrm{in}}C_pT_{\mathrm{in}}
-\dot m_{\mathrm{out}}C_pT
+\dot Q.</math>

En développant le membre de gauche et en utilisant :

<math>\dot m=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}},</math>

on obtient :

<math>
\dot T=
\frac{\dot m_{\mathrm{in}}(C_pT_{\mathrm{in}}-C_vT)
-\dot m_{\mathrm{out}}RT+\dot Q}
{C_vm}
.</math>

Contrôles limites :

* sans débit, l'équation redonne la relaxation thermique d'un volume fermé ;
* à débit stationnaire égal entrant/sortant, elle redonne <math display="inline">mC_v\dot T=\dot mC_p(T_{\mathrm{in}}-T)+\dot Q</math> ;
* si <math display="inline">T_{\mathrm{in}}=T</math> et <math display="inline">\dot m_{\mathrm{in}}=\dot m_{\mathrm{out}}</math>, la contribution nette du débit à <math display="inline">\dot T</math> s'annule.

<span id="A.7_Global_mass_test_during_an_active_phase"></span>
<div class="mw-translate-fuzzy">
== A.7 Test global de masse en phase active ==
</div>

Pour la Phase II :

<math>\dot m_L=-\dot m_{\mathrm{in}},</math>

<div class="mw-translate-fuzzy">
<math>\dot m_H=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}},</math>
</div>

<div class="mw-translate-fuzzy">
<math>\dot m_{SC}=\dot m_{\mathrm{out}}.</math>
</div>

En sommant :

<math>\dot M_{\mathrm{tot}}=0.</math>

La Phase IV donne exactement le même résultat par symétrie.

<span id="A.8_Global_energy_test_during_an_active_phase"></span>
<div class="mw-translate-fuzzy">
== A.8 Test global d'énergie en phase active ==
</div>

Pour la Phase II :

<math>\dot U_L=-P_L\dot V_L-\dot m_{\mathrm{in}}h_L,</math>

<div class="mw-translate-fuzzy">
<math>\dot U_H=\dot m_{\mathrm{in}}h_L-\dot m_{\mathrm{out}}h_H+\dot Q_H,</math>
</div>

<div class="mw-translate-fuzzy">
<math>\dot U_{SC}=\dot m_{\mathrm{out}}h_H+\dot Q_C-P_{SC}\dot V_S.</math>
</div>

Les flux enthalpiques internes s'annulent exactement :

<math>-\dot m_{\mathrm{in}}h_L+\dot m_{\mathrm{in}}h_L=0,</math>

<div class="mw-translate-fuzzy">
<math>-\dot m_{\mathrm{out}}h_H+\dot m_{\mathrm{out}}h_H=0.</math>
</div>

Il reste :

<div class="mw-translate-fuzzy">
<math>
\dot U_{\mathrm{tot}}
=\dot Q_C+\dot Q_H
-P_L\dot V_L-P_{SC}\dot V_S
.</math>
</div>

La Phase IV fournit la relation symétrique. Dans les phases fermées, la même structure résulte directement de la somme des bilans des deux paires. D'où, pour toute phase :

<div class="mw-translate-fuzzy">
<math>
\dot U_{\mathrm{tot}}
=\dot Q_C+\dot Q_H
-P_S^\star\dot V_S-P_L^\star\dot V_L
.</math>
</div>

Intégrée sur un cycle périodique, cette relation redonne :

<div class="mw-translate-fuzzy">
<math>Q_C+Q_H=W_{\mathrm{cycle}}.</math>
</div>

<span id="A.9_Continuity_at_transitions"></span>
<div class="mw-translate-fuzzy">
== A.9 Continuité aux transitions ==
</div>

À l'instant d'un événement de clapet, aucune masse ni énergie finie ne peut être transférée en un temps nul. Les variables conservées et la géométrie sont donc continues :

<div class="mw-translate-fuzzy">
<math>m_i^+=m_i^-,\qquad U_i^+=U_i^-,\qquad V_i^+=V_i^-.</math>
</div>

<div lang="en" dir="ltr" class="mw-content-ltr">
<math>m_j^+=m_j^-,\qquad U_j^+=U_j^-,\qquad V_j^+=V_j^-.</math>
</div>

Avec le gaz parfait :

<div class="mw-translate-fuzzy">
<math>T_i=\frac{U_i}{m_iC_v},\qquad P_i=\frac{m_iRT_i}{V_i},</math>
</div>

ce qui implique :

<div class="mw-translate-fuzzy">
<math>T_i^+=T_i^-,\qquad P_i^+=P_i^-.</math>
</div>

L'événement ne crée qu'un changement de topologie hydraulique et de système d'équations actif.

Thermodynamic and Mechanical Study/en

2026-09-07T14:59:46Z

FuzzyBot: Updating to match new version of source page

== 1. Scope and method ==

This study describes the complete thermodynamic cycle of the Dada engine in symbolic form, in both refrigeration and motor operation. The model couples prescribed kinematics of the two pistons, well-mixed gas volumes, 0D heat exchangers, passive check valves, and a generic compressible hydraulic closure model.

The cycle origin is set with the large cylinder at maximum volume, with:

<math>t=0,\qquad \theta(0)=0,\qquad V_L(0)=V_{L,\max}.</math>

The four gas volumes are:

* <math display="inline">S</math>: small cylinder;
* <math display="inline">L</math>: large cylinder;
* <math display="inline">H_i</math>: heat-in exchanger, located on the <math display="inline">S \to L</math> hydraulic path and transferring heat into the gas;
* <math display="inline">H_o</math>: heat-out exchanger, located on the <math display="inline">L \to S</math> hydraulic path and transferring heat out of the gas.

The hydraulic topology is fixed:

<math>S\to H_i\to L,</math>

<math>L\to H_o\to S.</math>

The passive check valves therefore always allow:

<math>H_i\to L,\qquad H_o\to S.</math>

The physical function of each heat exchanger is independent of the operating mode. What changes between refrigeration and motor operation is the external thermal reservoir connected to each exchanger.
The cycle comprises four hydraulic/thermodynamic phases:

* '''Phase I''' — check valves closed, nominally adiabatic: compression on the <math display="inline">L+H_o</math> side;
* '''Phase II''' — heat exchange, nominally isothermal: transfer <math display="inline">L \to H_o \to S</math>, with heat removed from the gas;
* '''Phase III''' — check valves closed, nominally adiabatic: expansion on the <math display="inline">S+H_i</math> side;
* '''Phase IV''' — heat exchange, nominally isothermal: transfer <math display="inline">S \to H_i \to L</math>, with heat supplied to the gas.

These phase definitions describe the permanent physical topology of the machine. The external reservoir associated with each heat exchanger depends on the operating mode. The phase boundaries and their angular positions are determined independently in each operating mode by the actual check-valve events.

The terms ''nominally adiabatic'' and ''nominally isothermal'' describe the objective of the cycle. They do not constitute exact thermodynamic constraints: temperatures are calculated from the balances, the heat exchangers remain coupled to their reservoirs, and the low-displacement regions of the pistons are not assumed to be perfectly stationary.

=== 1.1 First-level assumptions ===

The model is based on the following assumptions:

* single-phase, ideal and calorically perfect gas;
* constant properties <math display="inline">R</math>, <math display="inline">C_p</math>, <math display="inline">C_v</math>, <math display="inline">\gamma</math>, with <math display="inline">R=C_p-C_v</math> and <math display="inline">\gamma=C_p/C_v</math>;
* each gas volume is uniform and well mixed;
* gas kinetic and potential energies are neglected in the 0D balances;
* fixed heat-exchanger volumes;
* prescribed piston kinematics;
* no mechanical friction in the thermodynamic model;
* heat exchange represented by an overall conductance <math display="inline">UA</math>;
* passive check valves controlled by the pressure difference.

== 2. Notation and conventions ==

=== 2.1 Geometry and kinematics ===

For <math display="inline">k\in\{S,L\}</math>:

<math>V_{k,\min}>0,\qquad V_{k,\max}>V_{k,\min}.</math>

The swept volume is:

<math>V_{k,\mathrm{swept}}=V_{k,\max}-V_{k,\min}.</math>

The thermodynamic volumes prescribed by the mechanism are:

<math>V_S=V_S(t),\qquad V_L=V_L(t),</math>

with their '''signed''' derivatives:

<math>\dot V_S=\frac{dV_S}{dt},\qquad \dot V_L=\frac{dV_L}{dt}.</math>

The crank angular velocity is signed:

<math>\theta(t)=\omega t,\qquad
\dot V_k=\omega\frac{dV_k}{d\theta},
\qquad k\in\{S,L\}.</math>

The driven refrigeration direction is chosen as positive:

<math>\omega>0.</math>

Motor operation uses the opposite crank direction:

<math>\omega<0.</math>

The geometric origin is identical in both modes:

<math>t=0,\qquad \theta(0)=0,\qquad V_L(0)=V_{L,\max}.</math>

When only the absolute value of the volumetric speed is useful:

<math>\nu_{V,k}=|\dot V_k|.</math>

A '''quasi-stationary region''' denotes an interval in which the displacement or <math display="inline">|\dot V|</math> remains small compared with the transfer phases. This region corresponds to the “plateau” of the kinematic optimization, without assuming <math display="inline">\dot V=0</math> exactly.

=== 2.2 Kinematic closure fraction ===

To describe the normalized closure of a cylinder:

<math>\Lambda_k(t)=\frac{V_{k,\max}-V_k(t)}{V_{k,\max}-V_{k,\min}},\qquad k\in\{S,L\}.</math>

Thus <math display="inline">\Lambda_k=0</math> corresponds to maximum volume and <math display="inline">\Lambda_k=1</math> to minimum volume. The values <math display="inline">\Lambda_L^*</math> and <math display="inline">\Lambda_S^*</math> are nominal kinematic targets at the transitions. The value actually reached at a check-valve event is:

<math>\Lambda_{\mathrm{real}}=\Lambda(t_{\mathrm{event}}).</math>

The physical transitions remain determined by the pressures; <math display="inline">\Lambda^*</math> is therefore not an imposed opening condition.

=== 2.3 Thermodynamic variables ===

For each volume <math display="inline">j\in\{S,L,i,o\}</math>:

<math>m_j,\qquad U_j,\qquad T_j,\qquad P_j,\qquad V_j.</math>

The subscripts <math display="inline">i</math> and <math display="inline">o</math> denote respectively the gas contained in <math display="inline">H_i</math> and <math display="inline">H_o</math>.

The complete state vector is chosen as:

<math>\mathbf X=(m_S,U_S,m_L,U_L,m_i,U_i,m_o,U_o).</math>

Temperatures and pressures are derived from:

<math>T_j=\frac{U_j}{m_jC_v},\qquad
P_j=\frac{m_jRT_j}{V_j}.</math>

The volumes <math display="inline">V_S(t)</math> and <math display="inline">V_L(t)</math> are prescribed by the kinematics and are not independent thermodynamic state variables. The volumes <math display="inline">V_i</math> and <math display="inline">V_o</math> are constant.

=== 2.4 Energy sign convention ===

Heat is positive when it is '''received by the gas'''. Work is positive when the '''gas delivers work''':

<math>\dot W=P\dot V.</math>

By definition of the two physical heat exchangers:

<math>Q_i>0,\qquad Q_o<0</math>

in the intended operating regime of both refrigeration and motor operation.

The net cycle work distinguishes the two modes:

<math>W_{\mathrm{cycle}}<0</math>

for driven refrigeration operation, whereas:

<math>W_{\mathrm{cycle}}>0</math>

for motor operation.

== 3. Thermal closure and validity domain ==

=== 3.1 Exchange with the thermal reservoirs ===

For the heat-in exchanger:

<math>\dot Q_i=(UA)_i(T_{i,\mathrm{res}}-T_i).</math>

In the intended operating regime:

<math>T_i<T_{i,\mathrm{res}}
\quad\Rightarrow\quad
\dot Q_i>0.</math>

For the heat-out exchanger:

<math>\dot Q_o=(UA)_o(T_{o,\mathrm{res}}-T_o).</math>

In the intended operating regime:

<math>T_o>T_{o,\mathrm{res}}
\quad\Rightarrow\quad
\dot Q_o<0.</math>

The reservoir temperatures depend on the operating mode.

For refrigeration operation:

<math>T_{i,\mathrm{res}}=T_{\mathrm{cold}},\qquad
T_{o,\mathrm{res}}=T_{\mathrm{hot}}.</math>

For motor operation:

<math>T_{i,\mathrm{res}}=T_{\mathrm{hot}},\qquad
T_{o,\mathrm{res}}=T_{\mathrm{cold}}.</math>

Thus the heat-transfer equations themselves are identical in both modes.

A nominally isothermal phase therefore does not mean <math display="inline">T_{HX}=T_{HX,\mathrm{res}}</math>: a finite temperature difference is required to transfer finite thermal power when <math display="inline">UA</math> is finite.

An indicator of isothermal quality may be defined over a given phase by:

<math>\varepsilon_T=\frac{T_{\max}-T_{\min}}{T_{\mathrm{ref}}}.</math>

=== 3.2 Thermophysical validity domain ===

The base model assumes:

<math>PV=mRT,\qquad Z=1,</math>

<math>C_p=\mathrm{const},\qquad C_v=\mathrm{const},\qquad \gamma=\mathrm{const}.</math>

Validity must be checked a posteriori over the entire cycle, notably through:

<math>|Z-1|\ll1,</math>

and through small variations of the thermophysical properties, for example:

<math>\varepsilon_{C_p}=\frac{C_{p,\max}-C_{p,\min}}{C_{p,\mathrm{ref}}}\ll1.</math>

The working fluid must remain single-phase and gaseous, and sufficiently far from any condensation or phase transition throughout the <math display="inline">(P,T)</math> domain traversed.

If these criteria become insufficient, an extension may use <math display="inline">Z(P,T)</math>, <math display="inline">C_p(T)</math>, <math display="inline">C_v(T)</math>, or a real-gas equation of state without changing the general architecture of the mass and energy balances.

== 4. Reduced formulation of a quasi-pressure-equalized pair ==

When a cylinder and its heat exchanger are connected by a very low-resistance internal path, the approximation

<math>P_{\mathrm{cyl}}\approx P_{HX}=P_{\mathrm{pair}}.</math>

may be used. It is acceptable if:

<math>\varepsilon_P=\frac{|\Delta P_{\mathrm{int}}|}{P_{\mathrm{pair}}}\ll1,\qquad \Delta P_{\mathrm{int}}=P_{\mathrm{cyl}}-P_{HX}.</math>

A low internal Mach number provides an additional check:

<math>Ma_{\mathrm{int}}=\frac{|u_{\mathrm{int}}|}{a}\ll1,</math>

but it is not sufficient on its own to guarantee pressure equalization.

=== 4.1 Closed-pair case ===

For a closed pair, with fixed <math display="inline">V_{HX}</math> and <math display="inline">V=V_{\mathrm{cyl}}+V_{HX}</math>:

<math>
\frac{dP}{dt}=
\frac{(\gamma-1)(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX})-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

When <math display="inline">(UA)_{HX}=0</math>:

<math>P(V_{\mathrm{cyl}}+V_{HX})^\gamma=\mathrm{const}.</math>

=== 4.2 Internal redistribution flow rate ===

The internal flow rate is defined as positive from cylinder \to heat exchanger. Let:

<math>N=V_{HX}\dot P+(\gamma-1)(UA)_{HX}(T_{HX}-T_{HX,\mathrm{res}}).</math>

The flow carries the enthalpy of the upstream state:

<math>
\dot m_{\mathrm{int}}=
\begin{cases}
\dfrac{N}{\gamma R T_{\mathrm{cyl}}}, & N\ge0 \quad (\mathrm{cyl}\to HX),\\[6pt]
\dfrac{N}{\gamma R T_{HX}}, & N<0 \quad (HX\to\mathrm{cyl}).
\end{cases}
</math>

The heat-exchanger temperature evolves according to:

<math>
\dot T_{HX}=\frac{T_{HX}}{P}\dot P-
\frac{RT_{HX}^2}{PV_{HX}}\dot m_{\mathrm{int}}
.</math>

=== 4.3 Open receiving pair ===

During an active phase, the external flow physically enters the '''receiving cylinder''', not directly its associated heat exchanger. If <math display="inline">\dot m_{\mathrm{ext}}>0</math> enters the cylinder at temperature <math display="inline">T_{\mathrm{ext}}</math>:

<math>
\dot P=
\frac{\gamma RT_{\mathrm{ext}}\dot m_{\mathrm{ext}}
+(\gamma-1)(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX})
-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

The masses satisfy:

<math>\dot m_{\mathrm{cyl}}=\dot m_{\mathrm{ext}}-\dot m_{\mathrm{int}},\qquad
\dot m_{HX}=\dot m_{\mathrm{int}},</math>

and therefore:

<math>\dot m_{\mathrm{pair}}=\dot m_{\mathrm{ext}}.</math>

The closed case is obtained immediately with <math display="inline">\dot m_{\mathrm{ext}}=0</math>.

== 5. Hydraulic closure ==

=== 5.1 Generic formulation ===

Any hydraulic connection is described by a generic law:

<math>\dot m=\Phi(P_u,P_d,T_u,\mathcal G,\mathcal F)</math>

<math display="inline">u</math> and <math display="inline">d</math> respectively denote the upstream and downstream states. The transported enthalpy is that of the upstream state:

<math>\dot H_{\mathrm{mass}}=\dot m C_pT_u.</math>

For a bidirectional connection, the upstream state is determined by the actual direction of the pressure gradient. For a check valve, reverse flow is prohibited.

=== 5.2 First-level closure using a compressible orifice ===

A first approximation consists in using an effective hydraulic area:

<math>(C_dA)_{\mathrm{eff}},</math>

which represents the overall ease of gas flow through the actual connection.

With:

<math>r=\frac{P_d}{P_u},\qquad
r_{\mathrm{crit}}=\left(\frac{2}{\gamma+1}\right)^{\gamma/(\gamma-1)},</math>

the unchoked flow rate, for <math display="inline">r>r_{\mathrm{crit}}</math>, is:

<math>
\dot m=(C_dA)_{\mathrm{eff}}P_u
\sqrt{\frac{2\gamma}{RT_u(\gamma-1)}
\left(r^{2/\gamma}-r^{(\gamma+1)/\gamma}\right)}
.</math>

For <math display="inline">r\le r_{\mathrm{crit}}</math>:

<math>
\dot m=(C_dA)_{\mathrm{eff}}P_u
\sqrt{\frac{\gamma}{RT_u}}
\left(\frac{2}{\gamma+1}\right)^{\frac{\gamma+1}{2(\gamma-1)}}
.</math>

This closure is not essential to the model: it may later be replaced by a law that more faithfully represents the pressure losses of a real heat exchanger, pipe, or check valve.

== 6. Complete thermodynamic cycle ==

=== 6.1 Phase I — check valves closed, nominally adiabatic: compression on the L side ===

Both check valves are closed. The <math display="inline">L+H_o</math> and <math display="inline">S+H_i</math> pairs are closed. The <math display="inline">L+H_o</math> side is nominally compressed; the motion of the small piston remains that provided by the actual kinematics.

For <math display="inline">L+H_o</math>:

<math>\dot P_{Lo}=
\frac{(\gamma-1)(UA)_o(T_{o,\mathrm{res}}-T_o)-\gamma P_{Lo}\dot V_L}
{V_L+V_o}.</math>

For <math display="inline">S+H_i</math>:

<math>\dot P_{Si}=
\frac{(\gamma-1)(UA)_i(T_{i,\mathrm{res}}-T_i)-\gamma P_{Si}\dot V_S}
{V_S+V_i}.</math>

The internal redistribution and temperature equations of §4 apply to both pairs.

The transition to Phase II occurs when the <math display="inline">H_o \to S</math> check valve satisfies its opening condition.

=== 6.2 Phase II — heat exchange, nominally isothermal: L → Ho → S ===

The gas leaves <math display="inline">L</math>, passes through the heat-out exchanger <math display="inline">H_o</math>, where it rejects heat, crosses the <math display="inline">H_o \to S</math> check valve, and then enters the receiving cylinder <math display="inline">S</math>. The <math display="inline">S+H_i</math> pair remains quasi-pressure-equalized if the criterion <math display="inline">\varepsilon_P\ll1</math> is satisfied.

==== 6.2.1 Donor cylinder L ====

The fundamental balance is:

<math>
\frac{d(m_LC_vT_L)}{dt}
=-P_L\dot V_L-\dot m_{L,\mathrm{out}}C_pT_L
.</math>

The kinematics are designed to maintain an outflow from the donor cylinder. In this case, the analytical solution is:

<math>
\frac{T_L}{T_{L,\mathrm{ref}}}=
\left[\frac{m_L}{m_{L,\mathrm{ref}}}\frac{V_{L,\mathrm{ref}}}{V_L}\right]^{\gamma-1}
</math>

and:

<math>
P_L=P_{L,\mathrm{ref}}
\left[\frac{m_L}{m_{L,\mathrm{ref}}}\frac{V_{L,\mathrm{ref}}}{V_L}\right]^\gamma
.</math>

One also obtains:

<math>\frac{T_L}{T_{L,\mathrm{ref}}}=
\left(\frac{P_L}{P_{L,\mathrm{ref}}}\right)^{(\gamma-1)/\gamma},</math>

and the specific entropy of the remaining gas satisfies <math display="inline">ds=0</math> under these assumptions. If reverse flow occurs despite the kinematic design, the analytical solution is no longer applicable and the complete open-system balance in <math display="inline">(m_L,U_L)</math> must be used.

==== 6.2.2 Heat-out exchanger Ho ====

Mass conservation:

<math>\dot m_o=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}}.</math>

Fundamental energy balance:

<math>
\frac{d(m_oC_vT_o)}{dt}
=\dot m_{\mathrm{in}}C_pT_L
-\dot m_{\mathrm{out}}C_pT_o
+(UA)_o(T_{o,\mathrm{res}}-T_o)
.</math>

In expanded form:

<math>
\dot T_o=
\frac{\dot m_{\mathrm{in}}(C_pT_L-C_vT_o)
-\dot m_{\mathrm{out}}RT_o
+(UA)_o(T_{o,\mathrm{res}}-T_o)}
{C_vm_o}
.</math>

The flow rates are determined by the hydraulic laws:

<math>\dot m_{\mathrm{in}}=\Phi_{\mathrm{HX,o}}(P_L,P_o,T_L,\ldots),</math>

<math>\dot m_{\mathrm{out}}=\Phi_{\mathrm{valve,o}}(P_o,P_{Si},T_o,\ldots).</math>

==== 6.2.3 Receiving pair S+Hi ====

The flow from <math display="inline">H_o</math> enters '''S'''. The pair pressure satisfies:

<math>
\dot P_{Si}=
\frac{\gamma R\dot m_{\mathrm{out}}T_o
+(\gamma-1)(UA)_i(T_{i,\mathrm{res}}-T_i)
-\gamma P_{Si}\dot V_S}
{V_S+V_i}
.</math>

The internal redistribution equations of §4 remain unchanged: the external flow acts on heat exchanger <math display="inline">H_i</math> indirectly through the evolution of the pair pressure.

The transition to Phase III is the closing event of the <math display="inline">H_o \to S</math> check valve.

=== 6.3 Phase III — check valves closed, nominally adiabatic: expansion on the S side ===

Both check valves are closed. The <math display="inline">S+H_i</math> pair expands nominally; the <math display="inline">L+H_o</math> pair also remains closed. Neither piston is assumed to be strictly stationary.

For <math display="inline">S+H_i</math>:

<math>\dot P_{Si}=
\frac{(\gamma-1)(UA)_i(T_{i,\mathrm{res}}-T_i)-\gamma P_{Si}\dot V_S}
{V_S+V_i}.</math>

For <math display="inline">L+H_o</math>:

<math>\dot P_{Lo}=
\frac{(\gamma-1)(UA)_o(T_{o,\mathrm{res}}-T_o)-\gamma P_{Lo}\dot V_L}
{V_L+V_o}.</math>

The transition to Phase IV occurs when the <math display="inline">H_i \to L</math> check valve satisfies its opening condition.

=== 6.4 Phase IV — heat exchange, nominally isothermal: S → Hi → L ===

The gas leaves <math display="inline">S</math>, passes through the heat-in exchanger <math display="inline">H_i</math>, where it receives heat from its external reservoir, crosses the <math display="inline">H_i \to L</math> check valve, and then enters the receiving cylinder <math display="inline">L</math>. The <math display="inline">L+H_o</math> pair remains quasi-pressure-equalized if the criterion <math display="inline">\varepsilon_P\ll1</math> is satisfied.

==== 6.4.1 Donor cylinder S ====

The fundamental balance is:

<math>
\frac{d(m_SC_vT_S)}{dt}
=-P_S\dot V_S-\dot m_{S,\mathrm{out}}C_pT_S
.</math>

For outflow guaranteed by the kinematic design:

<math>
\frac{T_S}{T_{S,\mathrm{ref}}}=
\left[\frac{m_S}{m_{S,\mathrm{ref}}}\frac{V_{S,\mathrm{ref}}}{V_S}\right]^{\gamma-1}
</math>

and:

<math>
P_S=P_{S,\mathrm{ref}}
\left[\frac{m_S}{m_{S,\mathrm{ref}}}\frac{V_{S,\mathrm{ref}}}{V_S}\right]^\gamma
.</math>

Reverse flow requires returning to the complete open-system balance in <math display="inline">(m_S,U_S)</math>.

==== 6.4.2 Heat-in exchanger Hi ====

<math>\dot m_i=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}}.</math>

The fundamental energy balance is:

<math>
\frac{d(m_iC_vT_i)}{dt}
=\dot m_{\mathrm{in}}C_pT_S
-\dot m_{\mathrm{out}}C_pT_i
+(UA)_i(T_{i,\mathrm{res}}-T_i)
.</math>

In expanded form:

<math>
\dot T_i=
\frac{\dot m_{\mathrm{in}}(C_pT_S-C_vT_i)
-\dot m_{\mathrm{out}}RT_i
+(UA)_i(T_{i,\mathrm{res}}-T_i)}
{C_vm_i}
.</math>

The flow rates are determined by:

<math>\dot m_{\mathrm{in}}=\Phi_{\mathrm{HX,i}}(P_S,P_i,T_S,\ldots),</math>

<math>\dot m_{\mathrm{out}}=\Phi_{\mathrm{valve,i}}(P_i,P_{Lo},T_i,\ldots).</math>

==== 6.4.3 Receiving pair L+Ho ====

The flow from <math display="inline">H_i</math> enters '''L'''. The pair pressure satisfies:

<math>
\dot P_{Lo}=
\frac{\gamma R\dot m_{\mathrm{out}}T_i
+(\gamma-1)(UA)_o(T_{o,\mathrm{res}}-T_o)
-\gamma P_{Lo}\dot V_L}
{V_L+V_o}
.</math>

Closing the <math display="inline">H_i \to L</math> check valve returns the system to Phase I. Geometric closure of the pistons alone is not sufficient to guarantee thermodynamic closure of the cycle.

=== 6.5 Motor operation ===

No second set of thermodynamic balance equations is required for motor operation.

The physical hydraulic topology remains:

<math>L\to H_o\to S,\qquad
S\to H_i\to L.</math>

The check-valve directions remain:

<math>H_o\to S,\qquad H_i\to L.</math>

Motor operation is obtained by:

# reversing the crank direction,
# exchanging the external reservoirs connected to <math display="inline">H_i</math> and <math display="inline">H_o</math>.

Thus:

<math>\omega<0,</math>

with the same cycle origin:

<math>V_L(0)=V_{L,\max}.</math>

The thermal-reservoir assignment becomes:

<math>T_{i,\mathrm{res}}=T_{\mathrm{hot}},\qquad
T_{o,\mathrm{res}}=T_{\mathrm{cold}}.</math>

The heat-in exchanger therefore absorbs heat from the hot reservoir:

<math>Q_i>0,</math>

while the heat-out exchanger rejects heat to the cold reservoir:

<math>Q_o<0.</math>

The pressure histories, mass flow rates, check-valve events, and periodic thermodynamic state must be recalculated with the reversed kinematics.

The motor regime is obtained when:

<math>W_{\mathrm{cycle}}>0.</math>

== 7. Physical transitions of the check valves ==

For a check valve oriented from upstream <math display="inline">u</math> to downstream <math display="inline">d</math>:

<math>P_u-P_d\ge\Delta P_{\mathrm{open}}
\quad\Rightarrow\quad \text{opening},</math>

<math>P_u-P_d\le\Delta P_{\mathrm{close}}
\quad\Rightarrow\quad \text{closing},</math>

with hysteresis, if present:

<math>\Delta P_{\mathrm{close}}\le\Delta P_{\mathrm{open}}.</math>

An opening or closing event changes the hydraulic topology and therefore the active equations; it does not cause any instantaneous jump in the thermodynamic state. For each volume:

<math>m_j^+=m_j^-,\qquad U_j^+=U_j^-,\qquad V_j^+=V_j^-.</math>

For an ideal gas:

<math>T_j^+=T_j^-,\qquad P_j^+=P_j^-.</math>

There is therefore no instantaneous pressure equalization when a check valve opens.

The values <math display="inline">\Lambda^*</math> serve as kinematic design targets; the values actually observed at the transitions are <math display="inline">\Lambda_{\mathrm{real}}=\Lambda(t_{\mathrm{event}})</math>.

The two passive check valves have permanent physical orientations:

<math>H_o\to S,\qquad H_i\to L.</math>

Their opening and closing conditions are identical in both operating modes. Only the pressure histories change because the crank kinematics are reversed in motor operation.

== 8. Work, heat, and performance ==

The instantaneous work delivered by the gas on the two pistons is calculated during all phases:

<math>
\dot W=P_S^\star\dot V_S+P_L^\star\dot V_L
.</math>

<math display="inline">P_S^\star</math> and <math display="inline">P_L^\star</math> denote the thermodynamic pressure effectively applied to the gas in each cylinder according to the phase topology: pair pressure when the cylinder belongs to a quasi-pressure-equalized pair, and its own pressure when it is a hydraulically isolated donor.

The net work over the cycle is:

<math>
W_{\mathrm{cycle}}=
\int_0^\tau
\left(P_S^\star\dot V_S+P_L^\star\dot V_L\right)dt
.</math>

The exchanged heats are:

<math>Q_i=\int_0^\tau (UA)_i(T_{i,\mathrm{res}}-T_i)\,dt,</math>

<math>Q_o=\int_0^\tau (UA)_o(T_{o,\mathrm{res}}-T_o)\,dt.</math>

In periodic steady operation:

<math>\Delta U_{\mathrm{cycle}}=0,</math>

and the first law gives:

<math>Q_i+Q_o=W_{\mathrm{cycle}}.</math>

The refrigeration COP is:

<math>COP_c=\frac{Q_i}{-W_{\mathrm{cycle}}}.</math>

The heat-pump COP is:

<math>COP_h=\frac{-Q_o}{-W_{\mathrm{cycle}}}=COP_c+1.</math>

The signs <math display="inline">Q_i>0</math>, <math display="inline">Q_o<0</math>, and <math display="inline">W_{\mathrm{cycle}}<0</math> provide checks of the intended refrigeration regime.

For motor operation:

<math>Q_i>0,\qquad Q_o<0,\qquad W_{\mathrm{cycle}}>0.</math>

The thermal efficiency is:

<math>
\eta_{\mathrm{th}}
=
\frac{W_{\mathrm{cycle}}}{Q_i}
=
1+\frac{Q_o}{Q_i}.
</math>

The mean thermodynamic motor power is:

<math>
\overline{\dot W}
=
\frac{W_{\mathrm{cycle}}}{\tau}.
</math>

The thermodynamic force exerted by the gas on a piston face may be written <math display="inline">F_{\mathrm{gas}}=PS</math>. Net mechanical force, inertia, and friction belong to the subsequent mechanical sizing stage.

== 9. Gas charge and periodic regime ==

The total amount of enclosed gas is a physical parameter:

<math>M_{\mathrm{tot}}=m_S+m_L+m_i+m_o=\mathrm{const}.</math>

It may be imposed directly or defined experimentally by a charging pressure and temperature. The chosen reference configuration is that at <math display="inline">t=0</math>, with the large cylinder at maximum volume. If all volumes communicate and are in uniform equilibrium during charging:

<math>
M_{\mathrm{tot}}=
\frac{P_{\mathrm{charge}}
\left[V_S(0)+V_{L,\max}+V_i+V_o\right]}
{RT_{\mathrm{charge}}}
.</math>

<math display="inline">P_{\mathrm{charge}}</math> and <math display="inline">T_{\mathrm{charge}}</math> define the amount of gas charged; they are not conditions that the periodic cycle must recover.

The established periodic regime is a solution of the system such that, between two successive passages through the maximum volume of the large cylinder with the same kinematic direction:

<math>\mathbf X(t+\tau)=\mathbf X(t).</math>

Geometric periodicity alone:

<math>V_k(t+\tau)=V_k(t),\qquad k\in\{S,L\}.</math>

is not sufficient to guarantee thermodynamic periodicity.

The numerical state used to initialize a calculation may be approximate; it must not be confused with a physical parameter of the machine. The future solver may search for the periodic fixed point by successive cycles, a shooting method, or a Newton method.

== 10. Parameters, design data, and results ==

=== 10.1 Prescribed data ===

* operating mode;
* reservoir temperatures <math display="inline">T_{\mathrm{cold}}</math> and <math display="inline">T_{\mathrm{hot}}</math>;
* signed crank angular velocity <math display="inline">\omega</math>.
* working fluid and reference properties <math display="inline">R</math>, <math display="inline">C_p</math>, <math display="inline">C_v</math>, <math display="inline">\gamma</math>;
* total charge <math display="inline">M_{\mathrm{tot}}</math>, or equivalently <math display="inline">(P_{\mathrm{charge}},T_{\mathrm{charge}})</math> in the charging configuration defined in §9;
* kinematics <math display="inline">V_S(t)</math>, <math display="inline">V_L(t)</math>.

=== 10.2 Design parameters ===

* <math display="inline">V_{S,\min}</math>, <math display="inline">V_{S,\max}</math>, <math display="inline">V_{L,\min}</math>, <math display="inline">V_{L,\max}</math>;
* <math display="inline">V_i</math>, <math display="inline">V_o</math>;
* <math display="inline">(UA)_i</math>, <math display="inline">(UA)_o</math>;
* hydraulic parameters of the heat exchangers, pipes, and check valves, represented at first level by <math display="inline">(C_dA)_{\mathrm{eff}}</math>;
* thresholds <math display="inline">\Delta P_{\mathrm{open}}</math>, <math display="inline">\Delta P_{\mathrm{close}}</math>;
* kinematic targets <math display="inline">\Lambda_L^*</math>, <math display="inline">\Lambda_S^*</math>.

=== 10.3 Calculated variables and results ===

* <math display="inline">m_j,U_j,T_j,P_j,\qquad j\in\{S,L,i,o\}</math>;
* internal and external mass flow rates;
* <math display="inline">\dot Q_i</math>, <math display="inline">\dot Q_o</math>, <math display="inline">Q_i</math>, <math display="inline">Q_o</math>;
* <math display="inline">W_{\mathrm{cycle}}</math>, <math display="inline">COP_c</math>, <math display="inline">COP_h</math>;
* <math display="inline">\eta_{\mathrm{th}}</math> in motor operation.
* pressure, temperature, and flow-rate extrema;
* actual check-valve events and <math display="inline">\Lambda_{\mathrm{real}}</math>;
* isothermal quality <math display="inline">\varepsilon_T</math>;
* validity criteria <math display="inline">\varepsilon_P</math>, <math display="inline">Ma</math>, <math display="inline">Z</math>, and property variations.

== 11. Global conservation checks ==

=== 11.1 Mass conservation ===

The solver must satisfy:

<math>\frac{dM_{\mathrm{tot}}}{dt}=0.</math>

A useful numerical residual is:

<math>\varepsilon_M(t)=M_{\mathrm{tot}}(t)-M_{\mathrm{tot}}(0).</math>

=== 11.2 Global energy conservation ===

Whatever the phase, the internal mass and enthalpy fluxes must cancel when the balances of all volumes are summed. The global balance must reduce to:

<math>
\frac{dU_{\mathrm{tot}}}{dt}
=\dot Q_i+\dot Q_o
-P_S^\star\dot V_S
-P_L^\star\dot V_L
.</math>

A cumulative energy residual may be defined by:

<math>
\varepsilon_E(t)=
U_{\mathrm{tot}}(t)-U_{\mathrm{tot}}(0)
-Q_i(0,t)-Q_o(0,t)+W(0,t)
.</math>

Over a periodic cycle:

<math>Q_i+Q_o=W_{\mathrm{cycle}}.</math>

The solver must keep <math display="inline">\varepsilon_M</math> and <math display="inline">\varepsilon_E</math> close to zero to the expected numerical accuracy.

----

== Appendix A — Symbolic derivations and validated checks ==

=== A.1 Pressure equation for a closed pair ===

For a cylinder + heat-exchanger pair at quasi-uniform pressure:

<math>U=\frac{P(V_{\mathrm{cyl}}+V_{HX})}{\gamma-1}.</math>

The first law gives:

<math>\frac{dU}{dt}=\dot Q-P\dot V_{\mathrm{cyl}},</math>

with:

<math>\dot Q=(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX}).</math>

Differentiating <math display="inline">U</math>:

<math>\frac{1}{\gamma-1}
\left[(V_{\mathrm{cyl}}+V_{HX})\dot P+P\dot V_{\mathrm{cyl}}\right]
=\dot Q-P\dot V_{\mathrm{cyl}}.</math>

Hence:

<math>
\dot P=
\frac{(\gamma-1)\dot Q-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

If <math display="inline">\dot Q=0</math>:

<math>\frac{\dot P}{P}=-\gamma\frac{\dot V}{V},</math>

then:

<math>PV^\gamma=\mathrm{const}.</math>

=== A.2 Internal flow rate of the pair ===

For the heat exchanger alone, at fixed volume:

<math>U_{HX}=\frac{PV_{HX}}{\gamma-1}.</math>

Therefore:

<math>\frac{V_{HX}}{\gamma-1}\dot P
=(UA)_{HX}(T_{HX,\mathrm{res}}-T_{HX})
+\dot m_{\mathrm{int}}C_pT_{\mathrm{up}}.</math>

Using <math display="inline">C_p=\gamma R/(\gamma-1)</math>:

<math>V_{HX}\dot P+(\gamma-1)(UA)_{HX}(T_{HX}-T_{HX,\mathrm{res}})
=\gamma R T_{\mathrm{up}}\dot m_{\mathrm{int}}.</math>

This recovers the definition of the numerator <math display="inline">N</math> and the selection of the upstream temperature according to the sign of the flow rate.

=== A.3 Evolution of the heat-exchanger temperature within a pair ===

For a fixed volume:

<math>m_{HX}=\frac{PV_{HX}}{RT_{HX}}.</math>

Differentiating:

<math>\frac{\dot m_{HX}}{m_{HX}}
=\frac{\dot P}{P}-\frac{\dot T_{HX}}{T_{HX}}.</math>

With <math display="inline">\dot m_{HX}=\dot m_{\mathrm{int}}</math>:

<math>
\dot T_{HX}=\frac{T_{HX}}P\dot P-
\frac{RT_{HX}^2}{PV_{HX}}\dot m_{\mathrm{int}}
.</math>

This relation follows solely from mass conservation and the equation of state; it is valid for both flow directions.

=== A.4 Open receiving pair ===

For a pair receiving <math display="inline">\dot m_{\mathrm{ext}}</math> into its cylinder:

<math>\frac{dU}{dt}=\dot m_{\mathrm{ext}}C_pT_{\mathrm{ext}}
+\dot Q-P\dot V_{\mathrm{cyl}}.</math>

With <math display="inline">U=P(V_{\mathrm{cyl}}+V_{HX})/(\gamma-1)</math>:

<math>
\dot P=
\frac{\gamma RT_{\mathrm{ext}}\dot m_{\mathrm{ext}}
+(\gamma-1)\dot Q
-\gamma P\dot V_{\mathrm{cyl}}}
{V_{\mathrm{cyl}}+V_{HX}}
.</math>

The local mass balance:

<math>\dot m_{\mathrm{cyl}}=\dot m_{\mathrm{ext}}-\dot m_{\mathrm{int}},</math>

<math>\dot m_{HX}=\dot m_{\mathrm{int}},</math>

immediately gives:

<math>\dot m_{\mathrm{pair}}=\dot m_{\mathrm{ext}}.</math>

=== A.5 Analytical solution for the adiabatic donor cylinder ===

For an adiabatic, well-mixed cylinder with outflow only:

<math>\frac{d(mC_vT)}{dt}=-P\dot V-\dot m_{\mathrm{out}}C_pT,</math>

and:

<math>\dot m=-\dot m_{\mathrm{out}}.</math>

Expanding:

<math>C_vm\dot T+C_vT\dot m=-P\dot V+C_pT\dot m,</math>

thus:

<math>C_vm\dot T=-P\dot V+RT\dot m.</math>

With <math display="inline">P=mRT/V</math> and <math display="inline">R/C_v=\gamma-1</math>:

<math>\frac{dT}{T}=(\gamma-1)
\left(\frac{dm}{m}-\frac{dV}{V}\right).</math>

After integration:

<math>
\frac{T}{T_0}=
\left[\frac{m}{m_0}\frac{V_0}{V}\right]^{\gamma-1}
.</math>

Then, using <math display="inline">PV=mRT</math>:

<math>
P=P_0\left[\frac{m}{m_0}\frac{V_0}{V}\right]^\gamma
.</math>

and:

<math>
\frac{T}{T_0}=
\left(\frac{P}{P_0}\right)^{(\gamma-1)/\gamma}
.</math>

Under these assumptions, the specific entropy of the remaining gas is constant: <math display="inline">ds=0</math>. The total entropy of the gas contained in the cylinder is not constant because its mass varies.

=== A.6 Active heat exchanger: expanded balance ===

Fundamental balance:

<math>\frac{d(mC_vT)}{dt}
=\dot m_{\mathrm{in}}C_pT_{\mathrm{in}}
-\dot m_{\mathrm{out}}C_pT
+\dot Q.</math>

Expanding the left-hand side and using:

<math>\dot m=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}},</math>

one obtains:

<math>
\dot T=
\frac{\dot m_{\mathrm{in}}(C_pT_{\mathrm{in}}-C_vT)
-\dot m_{\mathrm{out}}RT+\dot Q}
{C_vm}
.</math>

Limiting checks:

* with no flow, the equation recovers the thermal relaxation of a closed volume;
* with equal steady inlet/outlet flow rates, it recovers <math display="inline">mC_v\dot T=\dot mC_p(T_{\mathrm{in}}-T)+\dot Q</math>;
* if <math display="inline">T_{\mathrm{in}}=T</math> and <math display="inline">\dot m_{\mathrm{in}}=\dot m_{\mathrm{out}}</math>, the net contribution of the flow to <math display="inline">\dot T</math> vanishes.

=== A.7 Global mass test during an active phase ===

For Phase II:

<math>\dot m_L=-\dot m_{\mathrm{in}},</math>

<math>\dot m_o=\dot m_{\mathrm{in}}-\dot m_{\mathrm{out}},</math>

<math>\dot m_{Si}=\dot m_{\mathrm{out}}.</math>

Summing:

<math>\dot M_{\mathrm{tot}}=0.</math>

Phase IV gives exactly the same result by symmetry.

=== A.8 Global energy test during an active phase ===

For Phase II:

<math>\dot U_L=-P_L\dot V_L-\dot m_{\mathrm{in}}h_L,</math>

<math>\dot U_o=\dot m_{\mathrm{in}}h_L-\dot m_{\mathrm{out}}h_o+\dot Q_o,</math>

<math>\dot U_{Si}=\dot m_{\mathrm{out}}h_o+\dot Q_i-P_{Si}\dot V_S.</math>

The internal enthalpy fluxes cancel exactly:

<math>-\dot m_{\mathrm{in}}h_L+\dot m_{\mathrm{in}}h_L=0,</math>

<math>-\dot m_{\mathrm{out}}h_o+\dot m_{\mathrm{out}}h_o=0.</math>

What remains is:

<math>
\dot U_{\mathrm{tot}}
=\dot Q_i+\dot Q_o
-P_L\dot V_L-P_{Si}\dot V_S
.</math>

Phase IV provides the symmetric relation. In the closed phases, the same structure follows directly from summing the balances of the two pairs. Hence, for any phase:

<math>
\dot U_{\mathrm{tot}}
=\dot Q_i+\dot Q_o
-P_S^\star\dot V_S-P_L^\star\dot V_L
.</math>

Integrated over a periodic cycle, this relation gives:

<math>Q_i+Q_o=W_{\mathrm{cycle}}.</math>

=== A.9 Continuity at transitions ===

At the instant of a check-valve event, no finite mass or energy can be transferred in zero time. The conserved variables and geometry are therefore continuous:

<math>\qquad j\in\{S,L,i,o\}</math>

<math>m_j^+=m_j^-,\qquad U_j^+=U_j^-,\qquad V_j^+=V_j^-.</math>

For an ideal gas:

<math>T_j=\frac{U_j}{m_jC_v},\qquad P_j=\frac{m_jRT_j}{V_j},</math>

which implies:

<math>T_j^+=T_j^-,\qquad P_j^+=P_j^-.</math>

The event only creates a change in hydraulic topology and in the active system of equations.

Translations:Thermodynamic and Mechanical Study/344/en

2026-09-07T14:59:19Z