Showing 1–15 of 15 results for author: Makuch, K

Tailoring the first law of thermodynamics for convective flows

Authors: Karol Makuch

Abstract: Equilibrium thermodynamics is grounded in the law of energy conservation, with a specific focus on how systems exchange energy with their environment during transitions between equilibrium states. These transitions are typically characterized by quantities such as heat absorption and the work needed to alter the system's volume. This study is inspired by the potential to develop an analogous, stra… ▽ More Equilibrium thermodynamics is grounded in the law of energy conservation, with a specific focus on how systems exchange energy with their environment during transitions between equilibrium states. These transitions are typically characterized by quantities such as heat absorption and the work needed to alter the system's volume. This study is inspired by the potential to develop an analogous, straightforward thermodynamic description for systems that are out of equilibrium. Here, we explore the global energy exchanges that occur during transitions between these nonequilibrium states. We study a system with heat flow and an external (gravity) field that exhibits macroscopic motion, such as Rayleigh-Bénard convection. We show that the formula for system's energy exchange has the same form as in equilibrium. It opens the possibility of describing out-of-equilibrium systems using a few simple laws similar to equilibrium thermodynamics. △ Less

Submitted 16 August, 2024; originally announced August 2024.

Comments: 6 pages, 1 figure

Direction of spontaneous processes in non-equilibrium systems with movable/permeable internal walls

Authors: Robert Hołyst, Paweł Jan Żuk, Anna Maciołek, Karol Makuch, Konrad Giżyński

Abstract: The second law of equilibrium thermodynamics explains the direction of spontaneous processes in a system after removing internal constraints. When the system only exchanges energy with the environment as heat, the second law states that spontaneous processes at constant temperature satisfy: $\textrm{d} U - δQ \leq 0$. Here, $\textrm{d} U$ is the infinitesimal change of the internal energy, and… ▽ More The second law of equilibrium thermodynamics explains the direction of spontaneous processes in a system after removing internal constraints. When the system only exchanges energy with the environment as heat, the second law states that spontaneous processes at constant temperature satisfy: $\textrm{d} U - δQ \leq 0$. Here, $\textrm{d} U$ is the infinitesimal change of the internal energy, and $δQ$ is the infinitesimal heat exchanged in the process. We will consider three different systems in a heat flow: ideal gas, van der Waals gas, and a binary mixture of ideal gases. We will also study ideal gas and van der Waals gas in the heat flow and gravitational field. We will divide each system internally into two subsystems by a movable wall. We will show that the direction of the motion of the wall, after release, at constant boundary conditions is determined by the same inequality as in equilibrium thermodynamics. The only difference between equilibrium and non-equilibrium law is the dependence of the net heat change, $δQ$, on the state parameters of the system. We will also consider a wall thick and permeable to gas particles and derive Archimedes' principle in the heat flow. Finally, we will study the ideal gas's Couette flow, where the direction of the motion of the internal wall follows from the inequality $\textrm{d} E - δQ - δW_s \leq 0$, with $\textrm{d} E$ being the infinitesimal change of the total energy (internal and kinetic) and $δW_s$ the infinitesimal work exchanged with the environment due to shear force. Ultimately, we will synthesize all these cases in a framework of the second law of non-equilibrium thermodynamics. △ Less

Submitted 29 April, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

The first law of thermodynamics in hydrodynamic steady and unsteady flows

Authors: Konrad Giżyński, Karol Makuch, Jan Paczesny, Paweł Żuk, Anna Maciołek, Robert Hołyst

Abstract: We studied planar compressible flows of ideal gas as models of a non-equilibrium thermodynamic system. We demonstrate that internal energy $U(S^{*},V,N)$ of such systems in stationary and non-stationary states is the function of only three parameters of state, i.e. non-equilibrium entropy $S^{*}$, volume $V$ and number of particles $N$ in the system. Upon transition between different states, the s… ▽ More We studied planar compressible flows of ideal gas as models of a non-equilibrium thermodynamic system. We demonstrate that internal energy $U(S^{*},V,N)$ of such systems in stationary and non-stationary states is the function of only three parameters of state, i.e. non-equilibrium entropy $S^{*}$, volume $V$ and number of particles $N$ in the system. Upon transition between different states, the system obeys the first thermodynamic law, i.e. $dU=T^{*}dS^{*}-p^{*}dV+μ^{*}dN$, where $U=3/2 NRT^{*}$ and $p^{*}V=NRT^{*}$. Placing a cylinder inside the channel, we find that U depends on the location of the cylinder $y_{c}$ only via the parameters of state, i.e. $U(S^{*}(y_{c}),V,N(y_{c}))$ at V=const. Moreover, when the flow around the cylinder becomes unstable, and velocity, pressure, and density start to oscillate as a function of time, t, U depends on t only via the parameters of state, i.e. $U(S^{*}(t),V,N(t))$ for V=const. These examples show that such a form of internal energy is robust and does not depend on the particular boundary conditions even in the unsteady flow. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 17 pages, 9 figures

Parameters of state in the global thermodynamics of binary ideal gas mixtures in a stationary heat flow

Authors: Anna Maciolek, Robert Holyst, Karol Makuch, Konrad Giżynski, Pawel J. Żuk

Abstract: We formulate the first law of global thermodynamics for stationary states of the binary ideal gas mixture subjected to heat flow. We map the non-uniform system onto the uniform one and show that the internal energy $U(S^*,V,N_1,N_2,f_1^*,f_2^*)$ is the function of the following parameters of state: a non-equilibrium entropy $S^*$, volume $V$, number of particles of the first component, $N_1$, numb… ▽ More We formulate the first law of global thermodynamics for stationary states of the binary ideal gas mixture subjected to heat flow. We map the non-uniform system onto the uniform one and show that the internal energy $U(S^*,V,N_1,N_2,f_1^*,f_2^*)$ is the function of the following parameters of state: a non-equilibrium entropy $S^*$, volume $V$, number of particles of the first component, $N_1$, number of particles of the second component $N_2$ and the renormalized degrees of freedom. The parameters $f_1^*,f_2^*$, $N_1, N_2$ satisfy the relation $x_1f_1^*/f_1+x_2f_2^*/f_2=1$ ($f_1$, where $x_i$ is the fraction of $i$ component, and $f_2$ are the degrees of freedom for each component respectively). Thus only 5 parameters of state describe the non-equilibrium state of the binary mixture in the heat flow. We calculate the non-equilibrium entropy $S^{*}$ and new thermodynamic parameters of state $f_1^*, f_2^*$ explicitly. The latter are responsible for heat generation due to the concentration gradients. The theory reduces to equilibrium thermodynamics, when the heat flux goes to zero. As in equilibrium thermodynamics, the steady-state fundamental equation also leads to the thermodynamic Maxwell relations for measurable steady-state properties. △ Less

Submitted 27 September, 2023; originally announced November 2023.

Comments: 8 pages, 1 figure

Steady state thermodynamics of ideal gas in shear flow

Authors: Karol Makuch, Konrad Giżyński, Robert Hołyst, Anna Maciołek, Paweł J. Żuk

Abstract: Equilibrium thermodynamics describes the energy exchange of a body with its environment. Here, we describe the global energy exchange of an ideal gas in the Coutte flow in a thermodynamic-like manner. We derive a fundamental relation between internal energy as a function of parameters of state. We analyze a non-equilibrium transition in the system and postulate the extremum principle, which determ… ▽ More Equilibrium thermodynamics describes the energy exchange of a body with its environment. Here, we describe the global energy exchange of an ideal gas in the Coutte flow in a thermodynamic-like manner. We derive a fundamental relation between internal energy as a function of parameters of state. We analyze a non-equilibrium transition in the system and postulate the extremum principle, which determines stable stationary states in the system. The steady-state thermodynamic framework resembles equilibrium thermodynamics. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Steady thermodynamic fundamental relation for the interacting system in a heat flow

Authors: Robert Hołyst, Karol Makuch, Konrad Giżyński, Anna Maciołek, Paweł J. Żuk

Abstract: There is a long-standing question of whether it is possible to extend the formalism of equilibrium thermodynamics to the case of non-equilibrium systems in steady states. We have made such an extension for an ideal gas in a heat flow [Hołyst \emph{et al.}, J. Chem. Phys. 157, 194108 (2022)]. Here we investigate whether such a description exists for the system with interactions: the Van der Waals g… ▽ More There is a long-standing question of whether it is possible to extend the formalism of equilibrium thermodynamics to the case of non-equilibrium systems in steady states. We have made such an extension for an ideal gas in a heat flow [Hołyst \emph{et al.}, J. Chem. Phys. 157, 194108 (2022)]. Here we investigate whether such a description exists for the system with interactions: the Van der Waals gas in a heat flow. We introduce the parameters of state, each associated with a single way of changing energy. The first law of non-equilibrium thermodynamics follows from these parameters. The internal energy $U$ for the non-equilibrium states has the same form as in equilibrium thermodynamics. For the Van der Waals gas, $U(S^*, V, N, a^*,b^* )$ is a function of only 5 parameters of state (irrespective of the number of parameters characterizing the boundary conditions): the entropy $S^*$, volume $V$, number of particles $N$, and the rescaled Van der Waals parameters $a^*$, $b^*$. The state parameters, $a^*$, $b^*$, together with $S^*$, determine the net heat exchange with the environment. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 7 pages, 1 figure

Thermodynamics of stationary states of the ideal gas in a heat flow

Authors: Robert Hołyst, Karol Makuch, Anna Maciołek, Paweł J. Żuk

Abstract: There is a long-standing question as to whether and to what extent it is possible to describe nonequilibrium systems in stationary states in terms of global thermodynamic functions. The positive answers have been obtained only for isothermal systems or systems with small temperature differences. We formulate thermodynamics of the stationary states of the ideal gas subjected to heat flow in the for… ▽ More There is a long-standing question as to whether and to what extent it is possible to describe nonequilibrium systems in stationary states in terms of global thermodynamic functions. The positive answers have been obtained only for isothermal systems or systems with small temperature differences. We formulate thermodynamics of the stationary states of the ideal gas subjected to heat flow in the form of the zeroth, first, and second law. Surprisingly, the formal structure of steady state thermodynamics is the same as in equilibrium thermodynamics. We rigorously show that $U$ satisfies the following equation $dU=T^{*}dS^{*}-pdV$ for a constant number of particles, irrespective of the shape of the container, boundary conditions, size of the system, or mode of heat transfer into the system. We calculate $S^{*}$ and $T^{*}$ explicitly. The theory selects stable nonequilibrium steady states in a multistable system of ideal gas subjected to volumetric heating. It reduces to equilibrium thermodynamics when heat flux goes to zero. △ Less

Submitted 21 October, 2022; v1 submitted 25 November, 2021; originally announced November 2021.

Internal energy in compressible Poiseuille flow

Authors: Konrad Gizynski, Karol Makuch, Jan Paczesny, Yirui Zhang, Anna Maciołek, Robert Holyst

Abstract: We analyse a compressible Poiseuille flow of ideal gas in a plane channel. We provide the form of internal energy U for a non-equilibrium stationary state (NESS) that includes viscous dissipation and pressure work. We demonstrate that U depends strongly on the ratio Δp/p_0, where Δp is the pressure difference between inlet and outlet and p_0 is the outlet's pressure. In addition, U depends on two… ▽ More We analyse a compressible Poiseuille flow of ideal gas in a plane channel. We provide the form of internal energy U for a non-equilibrium stationary state (NESS) that includes viscous dissipation and pressure work. We demonstrate that U depends strongly on the ratio Δp/p_0, where Δp is the pressure difference between inlet and outlet and p_0 is the outlet's pressure. In addition, U depends on two other variables: the channel aspect ratio and the parameter equivalent to Reynolds number. The stored internal energy, ΔU=U-U0, is small compared to the internal energy U0 of the equilibrium state (ES) for a moderate range of values of Δp/p_0. However, ΔU can become large for big Δp or close to vacuum conditions at the outlet (p_0~0 Pa). △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 17 pages, 9 figures

Continuous non-equilibrium transition driven by the heat flow

Authors: Yirui Zhang, Marek Litniewski, Karol Makuch, Pawel J. Zuk, Anna Maciolek, Robert Holyst

Abstract: We discovered an out-of-equilibrium transition in the ideal gas between two walls, divided by an inner, adiabatic, movable wall. The system is driven out-of-equilibrium by supplying energy directly into the volume of the gas. At critical heat flux, we have found a continuous transition to the state with a low-density, hot gas on one side of the movable wall and a dense, cold gas on the other side.… ▽ More We discovered an out-of-equilibrium transition in the ideal gas between two walls, divided by an inner, adiabatic, movable wall. The system is driven out-of-equilibrium by supplying energy directly into the volume of the gas. At critical heat flux, we have found a continuous transition to the state with a low-density, hot gas on one side of the movable wall and a dense, cold gas on the other side. Molecular dynamic simulations of the soft-sphere fluid confirm the existence of the transition in the interacting system. We introduce a stationary state Helmholtz-like function whose minimum determines the stable positions of the internal wall. This transition can be used as a paradigm of transitions in stationary states and the Helmholtz-like function as a paradigm of the thermodynamic description of these states. △ Less

Submitted 24 May, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: 8 pages, 6 figures; paragraph one rewritten; comments added regarding maxEP, and figures 4, 5 and 6 are added; references added and typos corrected

Journal ref: Phys. Rev. E 104, 024102 (2021)

Real-space dynamical mean-field theory of Friedel oscillations in strongly correlated electron systems

Authors: B. Chatterjee, J. Skolimowski, K. Makuch, K. Byczuk

Abstract: We study Friedel oscillations and screening effects of the impurity potential in the Hubbard model. Electronic correlations are accounted for by solving the real-space dynamical mean-field theory equations using the continuous time quantum Monte-Carlo simulations at finite temperatures and using a homogeneous self-energy approximation with the numerical renormalization group at zero temperature. W… ▽ More We study Friedel oscillations and screening effects of the impurity potential in the Hubbard model. Electronic correlations are accounted for by solving the real-space dynamical mean-field theory equations using the continuous time quantum Monte-Carlo simulations at finite temperatures and using a homogeneous self-energy approximation with the numerical renormalization group at zero temperature. We find that in the Fermi liquid phase both the amplitudes of Friedel oscillations and the screening charge decrease with increasing the interaction and follow the behavior of the Fermi liquid renormalization factor. Inside the Mott insulator regime the Friedel oscillations are absent but the residual screening charge remains finite. △ Less

Submitted 23 July, 2018; originally announced July 2018.

Comments: 13 pages, 17 figures

Journal ref: Phys. Rev. B 100, 115118 (2019)

Generalization of Clausius-Mossotti approximation in application to short-time transport properties of suspensions

Authors: Karol Makuch

Abstract: In 1983 Felderhof, Ford and Cohen gave microscopic explanation of the famous Clausius-Mossotti formula for the dielectric constant of nonpolar dielectric. They based their considerations on the cluster expansion of the dielectric constant, which relates this macroscopic property with the microscopic characteristics of the system. In this article, we analyze the cluster expansion of Felderhof, Ford… ▽ More In 1983 Felderhof, Ford and Cohen gave microscopic explanation of the famous Clausius-Mossotti formula for the dielectric constant of nonpolar dielectric. They based their considerations on the cluster expansion of the dielectric constant, which relates this macroscopic property with the microscopic characteristics of the system. In this article, we analyze the cluster expansion of Felderhof, Ford and Cohen by performing its resummation (renormalization). Our analysis leads to the ring expansion for the macroscopic characteristic of the system, which is an expression alternative to the cluster expansion. Using similarity of structures of the cluster expansion and the ring expansion, we generalize (renormalize) the Clausius-Mossotti approximation. We apply our renormalized Clausius-Mossotti approximation to the case of the short-time transport properties of suspensions, calculating the effective viscosity and the hydrodynamic function with the translational self-diffusion and the collective diffusion coefficient. We perform calculations for monodisperse hard-sphere suspensions in equilibrium with volume fraction up to 45%. To assess the renormalized Clausius-Mossotti approximation, it is compared with numerical simulations and the Beenakker-Mazur method. The results of our renormalized Clausius-Mossotti approximation lead to comparable or much less error (with respect to the numerical simulations), than the Beenakker-Mazur method for the volume fractions below $ φ\approx 30\% $ (apart from a small range of wave vectors in hydrodynamic function). For volume fractions above $φ\approx 30 \%$, the Beenakker-Mazur method gives in most cases lower error, than the renormalized Clausius-Mossotti approximation. △ Less

Submitted 29 September, 2015; v1 submitted 21 July, 2015; originally announced July 2015.

Journal ref: Phys. Rev. E 92, 042317 (2015)

Rotational self-diffusion in suspensions of charged particles: Revised Beenakker-Mazur and Pairwise Additivity methods versus numerical simulations

Authors: Karol Makuch, Marco Heinen, Gustavo C. Abade, Gerhard Nägele

Abstract: To the present day, the Beenakker-Mazur (BM) method is the most comprehensive statistical physics approach to the calculation of short-time transport properties of colloidal suspensions. A revised version of the BM method with an improved treatment of hydrodynamic interactions is presented and evaluated regarding the rotational short-time self-diffusion coefficient, $D^r$ , of suspensions of charg… ▽ More To the present day, the Beenakker-Mazur (BM) method is the most comprehensive statistical physics approach to the calculation of short-time transport properties of colloidal suspensions. A revised version of the BM method with an improved treatment of hydrodynamic interactions is presented and evaluated regarding the rotational short-time self-diffusion coefficient, $D^r$ , of suspensions of charged particles interacting by a hard-sphere plus screened Coulomb (Yukawa) pair potential. To assess the accuracy of the method, elaborate simulations of $D^r$ have been performed, covering a broad range of interaction parameters and particle concentrations. The revised BM method is compared in addition with results by a simplifying pairwise additivity (PA) method in which the hydrodynamic interactions are treated on a two-body level. The static pair correlation functions re- quired as input to both theoretical methods are calculated using the Rogers-Young integral equation scheme. While the revised BM method reproduces the general trends of the simulation results, it systematically and significantly underestimates the rotational diffusion coefficient. The PA method agrees well with the simulation data at lower volume fractions, but at higher concentrations $D^r$ is likewise underestimated. For a fixed value of the pair potential at mean particle distance comparable to the thermal energy, $D^r$ increases strongly with increasing Yukawa potential screening parameter. △ Less

Submitted 7 January, 2015; originally announced January 2015.

Comments: 24 pages, 13 figures

Journal ref: Soft Matter, 2015, 11, 5313-5326

Multipole matrix elements of Green function of Laplace equation

Authors: Karol Makuch, Przemysław Górka

Abstract: Multipole matrix elements of Green function of Laplace equation are calculated. The multipole matrix elements of Green function in electrostatics describe potential on a sphere which is produced by a charge distributed on the surface of a different (possibly overlapping) sphere of the same radius. The matrix elements are defined by double convolution of two spherical harmonics with the Green funct… ▽ More Multipole matrix elements of Green function of Laplace equation are calculated. The multipole matrix elements of Green function in electrostatics describe potential on a sphere which is produced by a charge distributed on the surface of a different (possibly overlapping) sphere of the same radius. The matrix elements are defined by double convolution of two spherical harmonics with the Green function of Laplace equation. The method we use relies on the fact that in the Fourier space the double convolution has simple form. Therefore we calculate the multipole matrix from its Fourier transform. An important part of our considerations is simplification of the three dimensional Fourier transformation of general multipole matrix by its rotational symmetry to the one-dimensional Hankel transformation. △ Less

Submitted 2 January, 2015; originally announced January 2015.

Journal ref: Acta Physica Polonica Series B 46(8), January 2015

Thermodynamic properties of correlated fermions in lattices with spin-dependent disorder

Authors: K. Makuch, J. Skolimowski, P. B. Chakraborty, K. Byczuk, D. Vollhardt

Abstract: Motivated by the rapidly growing possibilities for experiments with ultracold atoms in optical lattices we investigate the thermodynamic properties of correlated lattice fermions in the presence of an external spin-dependent random potential. The corresponding model, a Hubbard model with spin-dependent local random potentials, is solved within dynamical mean-field theory. This allows us to present… ▽ More Motivated by the rapidly growing possibilities for experiments with ultracold atoms in optical lattices we investigate the thermodynamic properties of correlated lattice fermions in the presence of an external spin-dependent random potential. The corresponding model, a Hubbard model with spin-dependent local random potentials, is solved within dynamical mean-field theory. This allows us to present a comprehensive picture of the thermodynamic properties of this system. In particular, we show that for a fixed total number of fermions spin-dependent disorder induces a magnetic polarization. The magnetic response of the polarized system differs from that of a system with conventional disorder. △ Less

Submitted 7 May, 2013; v1 submitted 14 February, 2013; originally announced February 2013.

Comments: 20 pages, 11 figures, published version

Journal ref: New J. Phys. 15, 045031 (2013)

Scattering series in mobility problem for suspensions

Authors: Karol Makuch

Abstract: The mobility problem for suspension of spherical particles immersed in an arbitrary flow of a viscous, incompressible fluid is considered in the regime of low Reynolds numbers. The scattering series which appears in the mobility problem is simplified. The simplification relies on the reduction of the number of types of single-particle scattering operators appearing in the scattering series. In our… ▽ More The mobility problem for suspension of spherical particles immersed in an arbitrary flow of a viscous, incompressible fluid is considered in the regime of low Reynolds numbers. The scattering series which appears in the mobility problem is simplified. The simplification relies on the reduction of the number of types of single-particle scattering operators appearing in the scattering series. In our formulation there is only one type of single-particle scattering operator. △ Less

Submitted 9 October, 2012; v1 submitted 21 August, 2012; originally announced August 2012.

Comments: 11 pages