Methods to Calculate Entropy Generation

Entropy generation, formulated by combining the first and second laws of thermodynamics with an appropriate thermodynamic potential, emerges as the difference between a phenomenological entropy function and a reversible entropy function. The phenomenological entropy function is evaluated over an irreversible path through thermodynamic state space via real-time measurements of thermodynamic states. The reversible entropy function is calculated along an ideal reversible path through the same state space. Entropy generation models for various classes of systems—thermal, externally loaded, internally reactive, open and closed—are developed via selection of suitable thermodynamic potentials. Here we simplify thermodynamic principles to specify convenient and consistently accurate system governing equations and characterization models. The formulations introduce a new and universal Phenomenological Entropy Generation (PEG) theorem. The systems and methods presented—and demonstrated on frictional wear, grease degradation, battery charging and discharging, metal fatigue and pump flow—can be used for design, analysis, and support of diagnostic monitoring and optimization.


Introduction
System analyses based on energy conservation alone have been shown inadequate for consistent characterization of real, often nonlinear, system transformation.The introduction of the second law via the works of Carnot, Clapeyron [1], Clausius, Massieu [2] and others established thermodynamics as a field for consistent description of changes in a system undergoing any form of energy conversion.For centuries, classical thermodynamics has been restricted to equilibrium and near-equilibrium transformations.Defining a minimum condition for a system to exist or process to occur in nature, Rayleigh's dissipation function [3], Onsager's least energy dissipation [4] (an application of his reciprocal relations of microscopic reversibility) and Prigogine's minimum entropy generation [5,6]-similar statements expressed as δS ′ = YdX/T-set the stage for thermodynamic characterization of real system transformation, a field commonly known as irreversible thermodynamics.Here, δS ′ is entropy generation, Y is generalized force, dX is generalized displacement and T is temperature.The correlation between energy dissipation (or entropy generation) and system degradation-advanced and permanent disorganization of material structure-has been theorized and experimentally verified.Recently, several multi-disciplinary system characterizations have emerged, presenting experimental results consistent with entropy-based formulations.These works show high accuracies in analyzing active system transformations, with inconsistencies attributable to the often-used "steady state" assumption.While some long-running systems operate predominantly in pseudo-steady state, e.g., very high-cycle fatigue of steels, most loaded systems often transform unsteadily, limiting Entropy 2024, 26, 237 2 of 27 the applicability of energy conservation and steady-state entropy characterizations.Recent works by Osara and Bryant [7][8][9][10] using the Degradation-Entropy Generation (DEG) methodology [11] to assess battery degradation, grease degradation and metal fatigue, showed a near 100% accurate and consistent characterization of these systems undergoing severely abusive loads.DEG methods relate increments of degradation to increments of entropy generation.Instead of a "steady state" assumption, Osara and Bryant [7][8][9][10] augmented laws of thermodynamics with the thermodynamic potentials to formulate entropy generation.Methods to calculate entropy generation under general conditions for all systems are needed to enable real-time assessment of system transformation.This article will develop entropy generation δS ′ for open and closed systems as the difference δS ′ = δS phen − dS rev ≥ 0 (1a) between a phenomenological entropy generation function δS phen , evaluated via suitable measurements of variables over an irreversible path in a thermodynamic state space between initial and final (or current) thermodynamic states; and a reversible entropy function dS rev , evaluated over a reversible path in the same state space between the same initial and final states.The thermodynamic state space consists of any independent variables that characterize the thermodynamic state and all active irreversible dissipative processes.The reversible path for dS rev , which is the projection of the irreversible phenomenological path of δS phen onto the reversible subspace, consists of the thermodynamic state variables of the phenomenological path sans the irreversible dissipative process variables.This approach will eliminate many tenuous "steady state" assumptions and loopholes in existing approaches.To obtain entropy generation S ′ , Equation (1a) must be integrated between the initial and final/current states, with entropy functions δS phen and dS rev integrated over their respective paths.True reversibility excludes process rates and time effects, hence, the reversible path defined herein marks the theoretical limit of a real-time-based process via a linear function joining the initial and final states.

Local Equilibrium
Prigogine posited: given that true equilibrium is asymptotic for all real systems, every continuous macroscopic system is made up of elements for which observable state properties (such as temperature and pressure) can be instantaneously determined or measured, thereby rendering equilibrium formulations describing these properties valid for each element in the macrosystem.Each element is, therefore, in local equilibrium [5,6,12,13].This theorem allows the extension of reversibility-based formulations to real systems, with entropy generation representation.While the system variables are spatial and temporal functions, many real systems operate with near-uniform internal properties, primarily changing with time.

States, Paths and Path (Line) Integrals
The following principles will be judicious to Section 2: 1.
Thermodynamics often involves changes in variables between two states.Variables include a set of the independent thermodynamic state variables Z chosen and measured for a particular system, and state dependent system properties which are functions of the states of Z.The Z characterize a system's thermodynamic state and can include temperature T, pressure P, and number of moles N, among others.Changes in system properties such as energy E, entropy S, temperature T and S rev use the exact differential d; are path independent, wherein changes in properties over an irreversible path (irr) are identical to changes in properties over a reversible path (rev), e.g., dE = dE irr = dE rev ; and the line integral Path-dependent variables such as work W, heat transfer Q, entropy generation S ′ and S phen depend on what occurs along the path between states o and f, and use the inexact differential δ such that f o δW = W must be accumulated over all instants of time t along the (assumed known) transformation path between times t o and t f .The path-dependent parameters will depend on a set of variables Z = {Z, ζ k } assumed to be time dependent, observable, and measurable.The Z characterize the thermodynamic state, whereas the ζ k characterize any active irreversible dissipative processes.Via a suitable numerical integration such as the trapezoid rule, with the Z (t) = {Z(t), ζ k (t)} measured as points Z t j = {Z(t j ), ζ k (t j ) suitably spaced at time instants t o < t j < t f in accord with the sampling theorem [14], the increments δW can be accumulated into the line integral where Z(t f ) and Z(t o ) must be measured or known at the beginning and final times t o and t f .With this, the line integral dZ dt, where sum index Z denotes a sum over all the components of Z, and dE/dZ must be evaluated at each time instant t o < t j < t f along the reversible path.

4.
For reversible processes dS rev with initial and final states dS revo and dS rev f , ) + dS revo satisfies the reversible path approximations of item 3.

5.
A phenomenological path (phen) through the thermodynamic state space enclosing Z = {Z, ζ k } defined in item 2 includes nonzero ζ k (t).A reversible path (rev) in the reversible subspace {Z} of Z [15] involves only the thermodynamic states Z rev , not the ζ k .The projection of the set of points Z t j = {Z t j , ζ k (t j ) that comprise an irreversible phenomenological path onto the reversible subspace is the set of points { Z rev t j [15].The system inputs, active mechanisms and dynamics are determined by the phenomenological path.
These principles are based on classical approaches [2,16,17].While classical theories are limited to equilibrium and near-equilibrium states, we here introduce concurrent physically observable states to characterize real, often highly dissipative, far-from-equilibrium systems.

Combining Internal Energy and Entropy Balances
For a stationary thermodynamic system (open and closed), excluding gravitational effects, the first law of thermodynamics balances dU, the change in overall system "internal" energy.Here ∑ δQ is the sum of all heat exchanges across the system boundary, ∑ δW is the sum of all work transfers across the system boundary, ∑ (u + Pv)dN e is the sum of energy transfers by matter flows dN e across the system boundary (for open systems only), and ∑ µ k dN k is the sum of all compositional energy changes within the system boundary.Embedded in the compositional change term ∑ µ k dN k are, in terms of mole number N k of species k, changes in the quantity of matter due to chemical reactions dN r k and ionic mass diffusions within the system boundary dN d k , i.e., dN k = dN r k + dN d k , where N r k and N d k are reactive and internally diffusive species, respectively.Also, u is molar internal energy, P is pressure, v is molar volume and µ is chemical potential.The product Pv is termed the flow work in open systems.Via the molar mass, Equation ( 2) can be re-written in terms of mass m.Inexact differential δ indicates path-dependent variables.
A statement of the second law-the Clausius inequality-gives the change in entropy of a closed system as dS ≥ δQ T where δQ/T is entropy flow by heat transfer, which can be positive or negative depending on the transfer direction, and T is the temperature of the boundary where the energy/entropy transfer takes place.Via the thermodynamic state principle-item 1 of Section 1.2-and the entropy balance [5,6,12], entropy change accompanying an open system process can be evaluated along a real and often nonlinear irreversible path (irr) as where δS Q = δQ T is entropy transfer via heat, dS N e = (u+Pv)dN e T is entropy transfer via mass, and δS ′ is internal entropy generation which always accompanies permanent change and structural disorganization of the evolving system.While dS, δS Q and dS N e can be positive or negative, the second law asserts δS ′ ≥ 0. Without a priori knowledge of the entropy generation, change along the irreversible path, Equation (3), cannot be determined.Along a reversible rev (ideal and linear) path with δS ′ = 0, Substituting δW = ∑ Y l dX l [18] for external/boundary works including compression PdV, strain work Vσdε, electrical work νdq and others into Equation ( 2) and combining with entropy Equation (3) gives the combined first and second laws as where Y l are intensive variables such as pressure P, strength/stress σ, voltage ν, etc.; X l are the system's extensive variables such as volume V, strain ε, charge q; and TδS U (S, X l , N k ), a path-dependent inexact differential, is defined by and equal to the middle expression of Equation (5).The script notation distinguishes S U as an entropy related function of the independent variables in parenthesis.Subscript phen indicates evaluation along the phenomenological path, the observable path where the independent states and dissipative process variables are available and measured at each instant.Since the product of temperature and entropy change (TdS) in Equation ( 5) subsumes the heat and flow transfer terms in Equations ( 2) and (3), Equation ( 5) is valid for all systems, open and closed [7][8][9][10]12].Equation (5), which governs the system along an irreversible path, has a pair of unknowns δS ′ > 0 and dU; all other terms are observable and can be measured, as shown later.
To solve, a second independent equation will arise via δS ′ = 0 [7-10] applied to Equation (5).Recalling the thermodynamic state principle, all process paths between the same initial and end states must have the same difference between states: dU = dU irr = dU rev , dS = dS irr = dS rev , etc., whether the path be reversible or irreversible.Substituting difference values into Equation (5) with δS ′ = 0, and solving yields dU = dU rev = TdS rev − ∑ Y l,rev dX l + ∑ µ k,rev dN k = TdS U,rev (S, X l , N k ), (6) where subscript rev indicates a quantity evaluated under reversible conditions.(Classical thermodynamics calculates entropy change along a reversible path, Equations ( 4) and (6), with reversible energy changes conveniently calculated via ideal response to loads, e.g., elastic deformation or non-saturating magnetics).Substituting dU = dU rev from Equation (6) into Equation ( 5) and rearranging yields entropy generation Entropy 2024, 26, 237 5 of 27 for all systems.The middle equalities of Equations ( 5) and ( 6) serve to evaluate the rightside terms of Equation (7).
In Equations ( 3) and ( 4) and forthcoming entropy formulations, summation signs ∑ indicating multiple heat and mass flow terms are omitted for convenience.

Entropy Generation of Various System Classes via Thermodynamic Potentials
To evaluate entropy generation, Equation (7) requires the system's changes of entropy dS and internal energy dU rev .These are difficult to determine, especially for non-thermal and multi-component systems [16], necessitating the widely used "steady state" assumption (dS = dU rev = 0) that led to Prigogine's stationary nonequilibrium transformation or minimum entropy generation δS ′ min = PdV T + ∑ µ k dN k T for a reacting compressible system [5,6,12].Via Legendre transforms, equivalent forms of Equations ( 5) and (7) are derived in terms of observable, measurable and more easily controllable system properties such as temperature and pressure [16] for all systems, open or closed, steady or unsteady.The thermodynamic potentials (with PV work replaced by generalized YX work)--enthalpy H = U + YX, Helmholtz free energy A = U-TS and Gibbs free energy G = U + YX-TS--are differentiated and solved for dU, the result of which is then inserted into Equation (5) to get: Analogous to the operations performed on Equation ( 5) that led to Equation (7), similar operations performed on Equations ( 8) and solved for δS ′ render Equations ( 7) and ( 9), of the form of Equation (1a), require: • δS phen : evaluated via the expressions of the middle equalities of Equations ( 5) and ( 8) at points over the irreversible phenomenological path defined in items 2 and 5 of Section 1.2, where salient quantities can be measured.• dS rev : evaluated as an exact differential over a reversible path between the beginning and final values of the irreversible path for δS phen , as discussed in list items 3, 4 and 5 of Section 1.2.Here, the independent states in the parentheses of Equations ( 5) and ( 8) are linear in time, as per items 3 and 4 of Section 1.2.

Phenomenological Entropy Generation Theorem
The results of Equations ( 7) and ( 9) can be summarized as the Phenomenological Entropy Generation (PEG) Theorem: Spontaneous entropy (and energy) changes along the (observable) phenomenological path is the sum of the ideal (linear, reversible transformation) entropy and the internally generated entropy, respectively: which can be solved for entropy generation, as in Equation (1a), repeated here as where, for energy extraction/decomposition or system loading, dS rev ≤ δS phen < 0, and for energy addition/formation, 0 < dS rev ≤ δS phen .Only end state measurements of system variables (before and after process interaction) are required for evaluation of dS rev (see items 3 and 4 of Section 1.2), unlike δS phen which requires an instantaneous account of all active processes.Note that a system's energy change and entropy change during a process can be negative or positive, depending on the direction of energy or entropy flow across system boundaries.Thermal systems include heat engines (thermal energy to mechanical work), heat pumps (mechanical work to thermal energy) and hydrocarbon fuels (chemical energy to thermal energy).Convenient for fuels and open systems is enthalpy H [16], which replaces volume with pressure as independent variable and measures the amount of thermal energy in a system.In a chemical reaction, change in enthalpy sums the heat absorbed or released by the reaction, by non-boundary deforming interactions, and by change in internal compositional energies.A fuel source's heating value is its enthalpy of combustion.For heat energy, the external work terms ∑ Y l dX l T in Equations (8a) and ( 7) are neglected and dH rev and dU rev , the maximum/minimum theoretical thermal energies available are often specified in tables, e.g., standard enthalpy of formation of pure substances [19], standard enthalpy of reaction, or standard enthalpy of combustion or heating value of fuels [20].The heat entropy change dS in Equations ( 7) and (9a) is evaluated via Equation (4) using the thermal energy balance δQ = CdT-for heat flow to or from a thermal system at uniform temperature with no heat generation-rendering dS rev = C rev dT T + (u+Pv) rev dN e T

Phenomenological Entropy Generation
. For a given transformation and flow, dH rev and dU rev are constants, hence the middle set of terms of Equations ( 7) and (9a) are extracted to obtain the phenomenological entropy generation functions: and C rev is the heat capacity (which can be obtained at standard room temperature and pressure) and (u + Pv) rev is the open-system flow enthalpy which, in the case of evaporation, is the latent heat of vaporization.Equation (10a) applies to a thermal system doing external work or receiving non-thermal energy across its boundary; Equation (10b) to a thermal system undergoing non-boundary-deforming internal transformation.Here, , the sum of the energy change due to combustion, nuclear or other exothermic/endothermic chemical reaction r, and diffusion d (as in a flame).

Boundary-Loaded (Work-Capable) Systems: Helmholtz Potential
Most electrical, structural and mechanical systems do not use or produce useful thermal energy but output work and/or use an external power supply.With Equations (10a) and (10b) inadequate for these systems, the Helmholtz free energy A [16] replaces entropy S with temperature T as independent variable, offers an adequate, consistent and convenient characterization, and measures maximum/minimum boundary work from/to a thermodynamic system.Referring to Equation (9b), dA rev is the maximum/minimum (theoretical) work possible.During work output, energy extraction or system loading, dT ≥ 0, dX l ≥ 0, dN k ≤ 0 and dA rev ≤ 0, rendering δS ′ ≥ 0. For work input, energy addition or product formation, dT ≤ 0, dX l ≤ 0, dN k ≥ 0 and dA rev ≥ 0, reversing the respective signs of terms in Equation (9b) to accord with the second law δS ′ ≥ 0 [5,8,9,19].A comparison of Equations ( 6) and (9b) shows that the Helmholtz relation conveniently absorbs dU rev and dS into dA rev and −S A dT, removing the need to measure heat or mass transfer across the system boundary and the need to determine dU rev , which is ambiguous for non-thermal systems.Most work-capable systems have a standardized maximum work dA rev obtainable, e.g., the elastic energy function for deformable solids.With a specified dA rev , δS ′ in Equation (9b) measures the irreversible entropy generation pertaining to dissipation of useful energy via work across the thermodynamic boundary, which requires the instantaneous evaluation of the Helmholtz phenomenological entropy generation [8,9] terms in Equation (8b), where Both δS H,phen and δA phen have consistent interpretations in all boundary-loaded systems, with terms composed of conjugate pairs involving physically observable and readily measurable changes in intensive and extensive system variables dT, dX l and dN k .For a non-reactive non-diffusive (dN r k = dN d k = 0) system, such as a lubricated mechanical interface, a fatigue-loaded component, and others, the last term of Equation (10c) can be neglected.Note that the very slow and/or (typically laboratory-controlled) isothermal case gives minimum entropy generation δS min = −∑ Y l dX l T .Various forms of work ∑ Y l dX l include frictional F f dx, electrical vdq, shear Vτdγ, compression PdV, and magnetic BdM, among others.

Internally Reactive Systems and Energy Storage Systems: Gibbs Potential
Reactive systems (chemical, nuclear, among others) undergo energy transformations via changes in composition.Energy storage and power sources such as batteries, nuclear power plants, and super capacitors involve changes in active species.The Gibbs free energy G in Equation (8c) replaces entropy S with temperature T and generalized position X with generalized force Y as independent variables [16], and measures maximum internal work or compositional (reactive) energy obtainable from a thermodynamic system.Equations (8c) and (9c) apply to all reactions, such as chemical formation/decomposition of substances, phase transitions, radioactive decay, etc.During active species consumption or system decomposition, dT ≥ 0, dY l ≤ 0, dN k ≤ 0 and dG rev ≤ 0, rendering δS ′ ≥ 0. For active species production or system formation, dT ≤ 0, dY l ≥ 0, dN k ≥ 0 and dG rev ≥ 0, reversing the signs of the respective terms in Equation (9c) to preserve δS ′ ≥ 0. With a specified dG rev (most energy systems have rated capacities and specific energies), Equation (9c) measures the actual irreversible entropy generation pertaining to the dissipation of useful energy via compositional changes.From Equation (8c), the Gibbs phenomenological entropy generation [7,10] where Both δS G,phen and δG phen have consistent interpretations in all systems undergoing active compositional reactions (charge, discharge or combinations) and are composed of conjugate pairs involving physically observable and readily measurable system variables dT, dY and dN k .For constant-pressure reactive system-process interactions such as cycling of electrochemical energy systems [7,10], the term XdY/T can be neglected.For non-reactive (dN k = 0) energy systems such as hydraulic/pressure accumulators, the last term in Equation (10d) can be dropped.

Generalized Material Properties, Entropy Content S, Internal Free Energy Dissipation "-SdT"
The irreversible Helmholtz and Gibbs fundamental relations, Equation (8b,c), introduced "-SdT", the portion of the free energy dissipated and accumulated internally by a loaded system, typically observed as the rise in temperature of the system under non-thermal loading.This can include effects of plastic work, friction, resistive Ohmic or Joule heat, chemical reaction heat generation, and sometimes heat from an external source.Temperature change dT is driven by the system's entropy content S. Without an entropy measurement device, S is often neglected or dT = 0, which requires experiments to be isothermal and/or significant temperature corrections for real-world applications.With δW = YdX, the energy-based Helmholtz and Gibbs equations suggest A = A(T, X, N) and G = G(T, Y, N).The entropy-based Massieu functions suggest S A = S A (T, X, N) and S G = S G (T, Y, N) wherein the entropy of a system depends on temperature T, generalized position X (for Helmholtz potential), generalized force Y (for Gibbs) and number of moles N, all of which are experimentally and instantaneously measurable.Via partial derivatives, total Helmholtz and Gibbs entropy changes Here, dN = dN e + dN k includes effects of mass flow and internal compositional changes/ chemical reactions, respectively.From Maxwell relations and Callen's derivatives reduction [16], Equations ( 11) can be re-stated using derived measurable system parameters [16,21,22], in terms of generalized work variables X, Y, as where C X > 0 and C Y > 0 are heat capacities (for solids, > 0 is generalized "isothermal loadability" [9], obtained via a reduction of the isothermal Gibbs derivative [16] = −Xκ T , a system/material property whose inverse defines the load modulus E' also derived from the second partial isothermal Helmholtz derivative (Appendix A).The equation α κ T = αE ′ = β > 0 is the thermal coefficient of generalized force (pressure, stress, voltage, etc.); λ X > 0 and λ Y > 0 are the coefficients of thermal chemico-transport decay (for solids λ X ≈ λ Y = λ) including the combined effects of internal reaction and mass flow on entropy content.
Heat capacity C measures the system's thermal response to heat transfer, retaining consistent meanings in all systems, and α measures non-thermal and non-chemical response (e.g., mechanical, electrical, etc.) to heat and temperature changes, obtained by defining YX for the specific system-process interaction.Generalized κ T -the inverse of the load-specific modulus E'-represents isothermal loadability, a measure of a material/system's "cold" response to boundary loading: loadability is compressibility (inverse of bulk modulus) for a compressible system [16], bendability for a beam under bending [9], shearability (inverse of shear modulus) for shearing or torsional loading [8], and conductance (inverse of resistance) for electrical work.Derivations and detailed discussions of the properties in Equations ( 12) are presented in Appendix A and in reference [22] specifically for a shear-loaded system.These formulations can be used to define new system-and processspecific material properties for assessing system/material performance and behavior.
By substituting Equations (12) into Equations ( 11), we have Integrating with an initial condition S 0 = 0 on entropy (valid in degradation analysis) gives Helmholtz and Gibbs entropy contents as functions of observable and measurable system-process phenomenological variables T, X, Y, N and material properties C, α, κ T , µ. Internal free energy dissipation via Helmholtz and Gibbs potentials are then While the system's state response and process variables (T,X,Y,N) are directly dependent on prevalent interaction rates (e.g., strain rates, electric currents, loads, etc.) and conditions (e.g., external source heating or cooling), the material properties (C, α, κ T , λ) can be assumed steady over a wide range of values of the state variables.Although the material properties may vary with changing state variables, the effects of such changes are minimal in a stable system with no discontinuities such as phase changes, severe chemical reactions, etc.For compositionally changing systems via chemical, nuclear or other reactions (conveniently characterized by the Gibbs potential), an alternate formulation is the Gibbs-Duhem equation [7,10,12,16]: Similar expressions can be established for other compositionally changing systems.For energy systems, the choice of Equation ( 16) or the second of Equations ( 15) depends on convenience and desired analysis output."-SdT" is termed MicroStructuroThermal (MST) energy dissipation [8,9] to suggest the dissipated energies as the source of heat.For electrochemical energy systems such as batteries and capacitors-where the last terms in Equation ( 16) and the second of Equations ( 15) are expressed via Faraday's electrolysis law in terms of cell charge capacity q and potential v-these equations are more specifically named ElectroChemicoThermal (ECT) energy dissipation [7,10].An application is presented in Section 4.3.
For reactive boundary-loaded (YdX) systems, substitute the first of Equations ( 15) into Equation (10c) to give and for internally reactive systems under displacement-controlled or non-boundary-deforming loading XdY, substitute the second of Equations ( 15) into Equation (10d) to give Equations ( 17) are posed in terms of the system's phenomenological variables, which are instantaneously measurable intensive and extensive system properties and process parameters that characterize the active phenomena along irreversible and reversible paths.

Stress vs. Strength Sign Conventions
For material property definitions, generalized force is interpreted as generalized useful force, strength or potential in line with the definitions of the free energies as maximum useful work obtainable from a system (Gibbs: internal work or compositional change; Helmholtz: external or boundary work).All the energy and entropy balances here accord with the IUPAC convention of representing energy leaving the system via work (and heat) as negative.In such a loaded system, dY is the decrease in strength, and hence is negative and dX is the increase in displacement.This accords with an expanding gas, for which pressure drops with increasing volume.However, in mechanics and other science/engineering fields that deal with solid materials, it is common to observe and use increase in stress as the system is loaded.In such cases, dY is positive.The derivations in this article and appendix consider dY < 0 the decrease in strength in a loaded system.

Helmholtz-Gibbs Coupling
Energy storage systems that provide boundary (external) work via direct interaction, e.g., batteries, capacitors, and pressure tanks, undergo internal changes driven by active external interaction.As such, internal phenomenological transformations can be monitored via boundary work measurements at the work transfer interface/terminal.The pressure stored in a hydraulic accumulator reduces as the accumulator provides external work, e.g., moves a weight over a distance.The phenomenological free energy change of an operational hydraulic accumulator can be expressed via the Gibbs potential as δG phen = −S G dT + VdP or via the Helmholtz potential as δA phen = −S A dT − PdV.For electrochemical energy systems, Osara and Bryant [7,10] presented a coupling of the internal chemical/diffusion kinetics with externally measured discharge/charge energy, µdN = −Vdq, to replace the phenomenological Gibbs relation δG phen = −S G dT + µdN with the more convenient phenomenological Helmholtz relation δA phen = −S A dT − Vdq.The Helmholtz and Gibbs relations for systems with interdependent internal and external interactions provides convenient characterization methods.

Rates
For application to time-based measurements, the phenomenological entropy generation functions of Equations (10)  . . . .
The dot notation represents time rate of change d( )/dt.
where the last term can be output turbine power or input compressor/pump power.Equation ( 19) can be used to monitor all boundary-loaded open systems in operation, including chemically reactive systems such as fuel cells-which have coupled boundary work, as discussed previously-simply by measuring inlet and exit flow rates, temperatures, generalized forces and velocities.Similar formulations for enthalpy can be derived, based on Equation (18b).

Evaluating Total Entropy Generation and Path (Line) Integrals
To estimate total entropy generation S', δS ′ in Equation (1a) must be integrated from the initial state o to the final state f in thermodynamic state space: where t o and t f are the times of the initial and final states, and the entropy generation functions δS phen and dS rev and their time rates .

S phen and
.
S rev are defined by the middle equality functions in Equations ( 5) and (8).The phenomenological entropy generation function δS phen (or .S phen ) is evaluated along the irreversible phenomenological path where states { Z, ζ k } are measured.The reversible entropy change function dS rev (or .S rev ) is evaluated along the reversible path with states { Z rev }.Proper selection of δS phen -the focus of Section 2-depends on system internal conditions and boundary loads.
To estimate S ′ = S ′ (t) at time t, t o < t < t f , replace t f in Equation ( 20) with t and let { Z rev (t)} be the projection of { Z(t), ζ k (t)}, i.e., the thermodynamic states { Z rev (t)} of the reversible path are related to their counterparts { Z(t), ζ k (t)} from the phenomenological path.With this, the integrals can be estimated via the methods of item 2 of Section 1.2.

Sample (Phenomenological) Entropy Generation Calculations
Sliding of copper against steel, shearing of grease, discharge and recharge of a lithiumion battery, fatigue of a steel rod, and flow through a pump will illustrate application of the entropy generation theory.

Friction Sliding of Copper against Steel at Steady Speed-(Steady State)
During a series of friction and wear tests, a copper rider pressed by 9.7 kg dead weight against a steel countersurface slid at steady speed .x = 3.3 ms −1 under carefully maintained thermal and lubricated boundary conditions [23].Measured were friction force F and temperatures at three locations in the copper, to estimate friction heat generation F .
x, heat flow and surface temperature T, all of which were steady during sliding to render .T = 0, which when substituted into Equation (18c) yields S ′ phen = − x T dt.Using measured values, the integral was evaluated over the test time interval to obtain the entropy generation plot in Figure 1.

Mechanical Shearing of Grease-Shear Stress and Shear Strain (Helmholtz Potential)
0.25 kg of Aeroshell 14 aircraft NLGI 4 lithium grease in a cup was sheared by a rotating impeller.The grease-in-cup system was treated as closed and non-reacting.Tests and procedures [8] at impeller speed 3 Hz measured impeller power M T ω (the product of torque M T and rotational speed ω) and temperatures (via thermocouples) of grease T and ambient, which are plotted versus time in Figure 2a.In terms of native grease internal variables, M T ω = Vτ .γ gives the shear power as the product of volume V, shear stress τ and shear strain rate

Mechanical Shearing of Grease-Shear Stress and Shear Strain (Helmholtz Potential)
0.25 kg of Aeroshell 14 aircraft NLGI 4 lithium grease in a cup was sheared by a rotating impeller.The grease-in-cup system was treated as closed and non-reacting.Tests and procedures [8] at impeller speed 3 Hz measured impeller power MTω (the product of torque MT and rotational speed ω) and temperatures (via thermocouples) of grease T and ambient, which are plotted versus time in Figure 2a.In terms of native grease internal variables,    gives the shear power as the product of volume , shear stress  and shear strain rate  .Helmholtz entropy content density SA-the first of Equations ( 14) divided by V-was substituted into Equation (18c) with  ˙ 0,  ,  ˙  , to yield phenomenological entropy density With the data of Figure 2a, the integrals in Equation ( 21) were estimated numerically via the methods of Section 1.2 to yield the entropy generation versus time plots in Figure 2b, where the MST entropy density and shear entropy density are the first and second integrals of Equation ( 21

Mechanical Shearing of Grease-Shear Stress and Shear Strain (Helmholtz Potential)
0.25 kg of Aeroshell 14 aircraft NLGI 4 lithium grease in a cup was sheared by a rotating impeller.The grease-in-cup system was treated as closed and non-reacting.Tests and procedures [8] at impeller speed 3 Hz measured impeller power MTω (the product of torque MT and rotational speed ω) and temperatures (via thermocouples) of grease T and ambient, which are plotted versus time in Figure 2a.In terms of native grease internal variables,    gives the shear power as the product of volume , shear stress  and shear strain rate  .Helmholtz entropy content density SA-the first of Equations ( 14) divided by V-was substituted into Equation (18c) with  ˙ 0,  ,  ˙  , to yield phenomenological entropy density With the data of Figure 2a, the integrals in Equation ( 21) were estimated numerically via the methods of Section 1.2 to yield the entropy generation versus time plots in Figure 2b, where the MST entropy density and shear entropy density are the first and second integrals of Equation ( 21), respectively.With the data of Figure 2a, the integrals in Equation ( 21) were estimated numerically via the methods of Section 1.2 to yield the entropy generation versus time plots in Figure 2b, where the MST entropy density and shear entropy density are the first and second integrals of Equation ( 21), respectively.v (applying Faraday's electrolysis law, see paragraph after Equation ( 16)), substituted into Equation (22) and integrated over time rendered, for discharge and charge, Here, q = t f t 0 Idt is the charge content.With the data of Figure 3a, the integrals in Equation ( 23) were estimated via the methods of Section 1.2 to yield the entropy generation versus time plots in Figure 3b, where the ElectroChemicoThermal ECT entropy and Ohmic entropy are the first and second integrals of Equation ( 23), respectively.Four 3.7 V, 11.5 Ah single-cell lithium-ion batteries were discharged at a variable discharge current of 5 A and recharged at a constant current of 3 A [7]. Figure 3a plots voltage , current I, temperatures of battery T and ambient, measured during the battery cycle.For a battery, the Helmholtz-Gibbs coupling, Section 2.7, yields ∑   ˙ , the Ohmic power.Substituting into Equation (19c) yielded Via the Gibbs-Duhem formulation, Equation ( 16), at constant pressure,  ˙  ˙ (applying Faraday's electrolysis law, see paragraph after Equation ( 16)), substituted into Equation (22) and integrated over time rendered, for discharge and charge, Here,    is the charge content.With the data of Figure 3a, the integrals in Equation (23) were estimated via the methods of Section 1.2 to yield the entropy generation versus time plots in Figure 3b, where the ElectroChemicoThermal ECT entropy and Ohmic entropy are the first and second integrals of Equation ( 23), respectively.

Fatigue of Metals-Stress and Strain (Helmholtz Potential)
A high-resolution infra-red camera monitored the temperature profile of an SS 304 stainless steel rod subjected to a 10 Hz displacement-controlled cyclic bending load until fatigue failure [24].See Figure 4a.Here, boundary work Y .X = Vσ : .ε, where σ is the stress tensor and ε is the elemental strain rate tensor, both having elastic and plastic components, viz σ = σ e + σ p , ε = ε e + ε p .Existing models [25] estimated the stress and strain.As in the case of grease shearing, Helmholtz entropy content density S A -the first of Equations ( 14) divided by volume V-was substituted into (18c) with .
Substituting stress, strain, temperature and material property values, the integrals in Equation ( 24) were estimated via the methods of Section 1.2 to render the entropy generation plots in Figure 4b, where the MST entropy density S ′ µT (red plot) and load entropy density S' W (blue plot) are the first and second integrals of Equation ( 24), respectively.For low-cycle fatigue, with significant plastic deformation, the modulus defined in Equations ( 12) and ( 13), substituted into entropy content density S A in Equation ( 24), is the hardness modulus.

Fatigue of Metals-Stress and Strain (Helmholtz Potential)
A high-resolution infra-red camera monitored the temperature profile of an SS 304 stainless steel rod subjected to a 10 Hz displacement-controlled cyclic bending load until fatigue failure [24].See Figure 4a.Here, boundary work  :  , where  is the stress tensor and  is the elemental strain rate tensor, both having elastic and plastic components, viz    ,    .Existing models [25] estimated the stress and strain.
As in the case of grease shearing, Helmholtz entropy content density SA-the first of Equations ( 14) divided by volume V-was substituted into (18c) with  ˙ 0 to obtain Substituting stress, strain, temperature and material property values, the integrals in Equation (24) were estimated via the methods of Section 1.2 to render the entropy generation plots in Figure 4b, where the MST entropy density ′ (red plot) and load entropy density S'W (blue plot) are the first and second integrals of Equation ( 24), respectively.For low-cycle fatigue, with significant plastic deformation, the modulus defined in Equations ( 12) and ( 13), substituted into entropy content density SA in Equation ( 24), is the hardness modulus.

Pump Flow-Pressure and Flow Rate (Internal Energy)
Water flowed through a three-phase 15-hp 260-gpm centrifugal motor-pump, instrumented to measure inlet and exit pump pressures and flow rate [26], which are plotted versus time in Figure 5a.A valve adjusted the flow rate.Measured volumetric flow rate was converted to mass flow rate using the density of water (  ), and pump power    , the product of torque  and rotational speed .Substituting into Equation (19), with positive input power, yielded With no external heat source and assuming negligible rise in flow temperature which was not measured during this test-hence T is constant ambient temperature-the first rightside term in Equation ( 25) vanished.Integrating with respect to time,  .
With no external heat source and assuming negligible rise in flow temperature which was not measured during this test-hence T is constant ambient temperature-the first right-side term in Equation ( 25) vanished.Integrating with respect to time, Substituting the data of Figure 5a, the integrals in Equation ( 26) yielded Figure 5b plots, where the flow entropy (measuring the change in the flow between inlet and exit, green plot) and load entropy (measuring the effect of power input into the pump, blue plot) are the first and second integrals of Equation ( 26), respectively.
Substituting the data of Figure 5a, the integrals in Equation ( 26) yielded Figure 5b plots, where the flow entropy (measuring the change in the flow between inlet and exit, green plot) and load entropy (measuring the effect of power input into the pump, blue plot) are the first and second integrals of Equation ( 26), respectively. (

Discussion
Here, reversible implies thermodynamic reversibility: any real system undergoing a spontaneous process cannot "revert" back to its original state without work from an external source, hence is thermodynamically irreversible.The reversible forms of Equations (9) (wherein  0 ) derived directly from Legendre transforms of entropy by Francois Massieu are called the Massieu functions [16,27].Here, by using the irreversible form of internal energy change (Equation ( 5)), Equations ( 9) are termed the irreversible Massieu functions.

Phenomenology
In open thermal systems,  in the internal energy and enthalpy Equations ( 7) and (9a), including the heat and mass transfer entropies, is analogous to   ⁄ in the Helmholtz and Gibbs energy Equation (9a,b).
Comparing the entropy Equations ( 3) and ( 4) and the corresponding energy counterparts in Equations ( 5), (6) ( 8) and ( 9), via the state principle, showed that changes in entropy and energy between two states are path-independent, whether the process path is reversible or irreversible, i.e.,     ′ ;     ′, where E is any of internal energy U, enthalpy H, Helmholtz potential A or

Gibbs potential G.
With a transformation representable along the linear and nonlinear paths, most thermodynamic characterizations employ energy change   and entropy change   : only end state measurements of system variables (before and after process interaction) are required for evaluation; unlike    ′ and     which require instantaneous account of all active processes.Note that a system's energy change dE and entropy change dS during a process can be negative or positive, depending on the direction of energy or entropy flow across system boundaries, hence neither dE nor dS measures the permanent changes in the system.On the other

Discussion
Here, reversible implies thermodynamic reversibility: any real system undergoing a spontaneous process cannot "revert" back to its original state without work from an external source, hence is thermodynamically irreversible.The reversible forms of Equations ( 9) (wherein δS ′ = 0) derived directly from Legendre transforms of entropy by Francois Massieu are called the Massieu functions [16,27].Here, by using the irreversible form of internal energy change (Equation ( 5)), Equations ( 9) are termed the irreversible Massieu functions.

Phenomenology
In open thermal systems, dS in the internal energy and enthalpy Equations ( 7) and (9a), including the heat and mass transfer entropies, is analogous to SdT/T in the Helmholtz and Gibbs energy Equation (9a,b).
Comparing the entropy Equations ( 3) and ( 4) and the corresponding energy counterparts in Equations ( 5), (6) ( 8) and ( 9), via the state principle, showed that changes in entropy and energy between two states are path-independent, whether the process path is reversible or irreversible, i.e., dS = dS rev = dS irr = δS phen − δS ′ ; dE = dE rev = dE irr = δE phen − TδS ′ , where E is any of internal energy U, enthalpy H, Helmholtz potential A or Gibbs potential G.
With a transformation representable along the linear and nonlinear paths, most thermodynamic characterizations employ energy change dE = dE rev and entropy change dS = dS rev : only end state measurements of system variables (before and after process interaction) are required for evaluation; unlike dE = dE irr = δE phen − TδS ′ and dS = dS irr = δS phen − δS ′ which require instantaneous account of all active processes.Note that a system's energy change dE and entropy change dS during a process can be negative or positive, depending on the direction of energy or entropy flow across system boundaries, hence neither dE nor dS measures the permanent changes in the system.On the other hand, entropy generation, Equations (1a), ( 7) and ( 9), evolves monotonically as stipulated by the second law.Figure 6 depicts the Phenomenological Entropy Generation theorem for a loaded or spontaneously transforming system.Equation (1b), which measures the entropy generated by the system's internal irreversibilities alone, is in accordance with experience, appearing similar to the Gouy-Stodola theorem of availability (exergy) analysis [21,[28][29][30].Rearranging Equation (1b) renders another statement of the second law or entropy balance as which replaces entropy transfer in Prigogine's entropy balance, Equation ( 3), with phenomenological entropy, from which entropy generation is subtracted.As discussed previously, the quasi-reversible terms in the foregoing formulations proceed at constant, often standardized or predetermined rates for a given transformation, making them negligible in instantaneous energy dissipation/degradation monitoring for which the phenomenological terms are used.Table 1 summarizes phenomenological entropy generations derived in this article for various classes of open and closed systems.According to Onsager and Prigogine, the minimum entropy generation rate for a system to exist (minimally active) is the quotient of its primary work interaction (or energy transfer) and boundary temperature, i.e., .

FIGURE CORRECTIONS FOR METHODS TO CALCULATE ENTROPY GENERATION
T [5,6,12,13].The minimum entropy generation rates for internally reactive non-thermal open systems, including terms that characterize every significant interaction, can be derived from Equations (18a) and (19) .
Equations ( 28) only apply to steady interactions where temperature is controlled or assumed to be constant (dT ≈ 0).Equations (28) present the minimum conditions for simultaneous real interactions to occur (i.e., minimum perturbation from equilibrium).The steady-state frictional wear in Section 4.1 demonstrates minimum entropy generation.
The MST entropy of Equation ( 29) accompanies the primary interactions defined by Minimum Entropy Generation (MEG) in Equations ( 28).This accords with the thermodynamic state postulate [12,16,31] which requires r + 1 independent, intensive properties to fully specify the state of a simple system undergoing r primary work interactions.In non-thermal systems, the MST entropy measures the effects of energy dissipated as heat in the system, which is therefore unavailable for work.As such, the MST entropy must be minimized to decelerate degradation of non-thermal systems, the limit of which is the MEG (where .
In the sample demonstrations in Figures 3b and 4b, the ECT/MST entropy (red plots) is significantly lower than the primary interaction entropy or MEG (blue plots), which may have justified prior approaches neglecting the former via an order of magnitude analysis.However, recent works [7][8][9][10] have shown that the ECT/MST entropy, measuring the free energy dissipation, contributes significantly to degradation.

Dissipation Factor and Entropic Efficiency
To measure a non-thermal system's dissipation tendencies relative to useful work output, define the dissipation factor the ratio of the MicroStructuroThermal MST entropy S ′ µT to the sum of boundary work S ′ W and compositional change S ′ N entropies, which can consistently characterize system response to dissipative mechanisms as well as multiple systems undergoing the same output/input work.A low J (minimal dissipation relative to available work) is desired for optimum performance and durability.
An examination of entropy Equations ( 7) and (9) indicates that low entropy generation ( δS ′ → 0 ) would make more useful system energy available, the limit of which is the reversible system (δS ′ = 0).Define entropic efficiency the ratio of the sum of boundary work and compositional change entropies to reversible entropy S rev .With an ideal (reversible or perfect) system-for which S ′ W + S ′ N = S rev -establishing 100% efficiency, a high η S ′ is preferred for slow degradation and optimum performance.Equation (31) appears similar to the exergy-based second-law efficiency which uses reversible work and boundary work at end states (i.e., not instantaneously determined during process interaction) [31,32].

Thermal vs. Non-Thermal Systems-Internal Energy and Enthalpy vs. the Free Energies
Internal energy-and enthalpy-based entropy generations, Equations ( 7) and (9a), have an entropy change term dS, which includes heat transfer entropy, accentuating the already known suitability and convenience of internal energy and enthalpy (heat content) in characterizing primarily heat-based (or thermal) systems.The Helmholtz-and Gibbs-based entropies of Equation (9b,c) encapsulate the effects of entropy transfers in the observable evolution of the system's internal variables (material properties and temperature).This makes the free energies and the consequent "free entropies" particularly convenient for non-thermal systems and interactions where heat transfer is not readily measurable and for which thermal mechanisms primarily emanate from energy dissipation.The free entropies are more consistent for assessing system transformations from observable and measurable system variables, irrespective of surrounding conditions.
Significant increase in temperature reduces available free energy via the MicroStructur-oThermal MST term, see Equation (8b,c).For non-thermal (mechanical, structural, chemical, electronic, magnetic, electrical, etc.) systems, this accords with experience, further making the free entropy formulations derived from the free energies suitable for non-thermal system characterizations.Thermal systems are utilized for heat content and typically have high temperature increase rate.This indicates that a high MST component-which accompanies high thermal energy-is favorable for thermal systems.Hence, the free energies are subject to misinterpretation and not recommended for characterizing thermal systems.

Summary and Conclusions
Building on the first-principles foundations of modern irreversible thermodynamics laid by Clausius, Rayleigh, Onsager and Prigogine, this article formulated universally consistent, instantaneous entropy generations for diverse macroscopic system categories (see summary in Table 1).Presented was thermodynamic resolution of the active processes and unsteady system responses during loading using readily evaluated entropy generation.A Phenomenological Entropy Generation (PEG) theorem was derived and proposed, expounding the significances of the newly introduced MicroStructuroThermal MST entropy (or ElectroChemicoThermal ECT entropy for electrochemical power systems) and previously neglected reversible entropy to characteristic entropy generation.Extending and generalizing Gibbs theory of thermodynamic stability, the Clausius inequality, Rayleigh's energy dissipation principle, Onsager's reciprocity, and Prigogine's entropy balance-the hallmarks of classical and modern irreversible thermodynamics [4][5][6]12,13,33]-this article, using results from recently published experimental works [7][8][9][10] • along the phenomenological path, thermodynamic and other states germane to evaluating S ′ phen and S rev are observable (measurable), permitting evaluation of entropy generation δS ′ ; • entropy generation is always non-negative in accordance with the second law, while its constituent terms S ′ phen and S rev are directional, negative for a loaded system and positive for an energized system.This implies S ′ phen ≤ |S rev | during load application/work output or active species decomposition, and S ′ phen ≥ |S rev | during energization/work input or active species formation, in accordance with experience and thermodynamic laws.In plain words, measurable energy obtained from a real system is always less than the theoretical maximum/reversible energy, and measurable energy added to a real system is always more than the minimum/reversible energy.(Modulus signs indicate magnitudes only).

Data Availability Statement:
No new data were created or analyzed in this study.Data analyzed here are available in the referenced works from which they were reproduced.Data sharing is not applicable to this article.

Conflicts of Interest:
Author Michael D. Bryant is employed by Machine Essence.The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.A13) and (A15) shows that E ′ = 1 κ T , i.e., the load modulus is the inverse of the isothermal loadability (as the bulk modulus is the inverse of the isothermal compressibility in compressible systems).

Abbreviations
According with the sign convention for a loaded system discussed in Section 2.6, the negative signs in Equations (A13) and (A15) indicate that an increase in generalized displacement X corresponds to a decrease in generalized strength/force Y, a stability condition for all loadable/work-capable systems.For example, as a compressible system expands (energy extraction or system loading), its accumulated/internal pressure decreases, i.e., ∂V > 0 ⇒ ∂P < 0, and vice versa for compression (energy addition), to give ∂P ∂V < 0. In mechanical systems, loading reduces component strength or toughness, ∂ε > 0 ⇒ ∂σ < 0 to give ∂σ ∂ε < 0. In electrochemical systems such as batteries, discharge (energy extraction via charge transfer) reduces electrochemical potential or voltage, dq > 0 ⇒ dV < 0 yielding ∂V ∂q < 0. The loadability or modulus of most solids and liquids (the latter, to much less extent) is not significantly affected by temperature over a wide temperature range, yielding κ T ≈ κ for most practical purposes.defined for a system as the displacement induced in the system by an increase in temperature at constant force/strength.The thermal displacement coefficient α measures the system's natural and characteristic mechanical response to temperature change via displacement change or deformation, and has universally consistent meaning in all systems.Equation (A16) implies that the system's Gibbs entropy content S G must increase in response to decrease in generalized force Y for process continuity and system stability.It also establishes that the position/displacement X in the system spontaneously increases (e.g., expansion, strain, etc.) with increase in temperature for stability and process continuity.
The second partial derivative of the Helmholtz free energy with respect to temperature T and generalized displacement X defines the thermal force (strength) coefficient β.Substituting the first and second equalities of Equation (A6) into the mixed second partial Helmholtz derivative yields rearranged to give the (solid) system's thermal force (strength) coefficient the ratio of the thermal displacement coefficient to the isothermal loadability or the product of the thermal displacement coefficient and the load modulus.β is the force/strength/stress induced by a degree rise in temperature at constant displacement.For stress-strain systems, β is also called the thermal stress coefficient [22].Equation (A18) implies that the system's Helmholtz entropy content must increase in response to increase in displacement (or deformation, expansion, etc.) for process continuity.In general, a solid material's strength spontaneously decreases with increase in temperature, ∂Y/∂T < 0, whereas for fluids, ∂Y/∂T > 0 due to the increase in the useful kinetic and thermal energy of fluids with temperature.
f o dE = ∆E = E f − E odepends only on the property values at the beginning and end states (o and f ).This is the thermodynamic state principle.
will be expressed in rate forms.With transport of active species into and out of an open system, flow rate .N e replacing dN e , Equations(10) in rate forms become: .

For
one-dimensional flow in open systems such as pumps, compressors, turbines and heat exchangers, among others, the second right-side term in Equation (18a,b) can be apthe control volume must equal the amounts in and out).Here, h = u + Pv is the specific standard enthalpy of the flowing fluid, and the subscripts inlet and exit denote pertinent quantities at those ports.The molar flow rate .N e can be converted to mass flow rate by multiplying by molar mass.Then, the internal energy-based phenomenological entropy rate, Equation (18a), in the absence of chemical or other reactions, becomes Entropy 2024,

Figure 2 .
Figure 2. (a) Monitored parameters during grease shearing and (b) phenomenological entropy density terms: shear entropy density and MST entropy density over time [8].Figure 2. (a) Monitored parameters during grease shearing and (b) phenomenological entropy density terms: shear entropy density and MST entropy density over time [8].

Figure 2 .
Figure 2. (a) Monitored parameters during grease shearing and (b) phenomenological entropy density terms: shear entropy density and MST entropy density over time [8].Figure 2. (a) Monitored parameters during grease shearing and (b) phenomenological entropy density terms: shear entropy density and MST entropy density over time [8].

4. 3 .
Discharge of Lithium-Ion Battery-Voltage and Charge (Helmholtz-Gibbs Coupling) Four 3.7 V, 11.5 Ah single-cell lithium-ion batteries were discharged at a variable discharge current of 5 A and recharged at a constant current of 3 A [7]. Figure 3a plots voltage v, current I, temperatures of battery T and ambient, measured during the battery Entropy 2024, 26, 237 13 of 27 cycle.For a battery, the Helmholtz-Gibbs coupling, Section 2.7, yields ∑ µ k .N k = vI, the Ohmic power.Substituting into Equation (19c) yielded .-Duhem formulation, Equation (16), at constant pressure, −S .T = q .

Figure 3 .
Figure 3. (a) Monitored parameters during Li-ion battery cycling showing 1.5 h discharge starting at ~5 A, followed by 1.5 h charge at 3 A. (b) phenomenological entropy terms: Ohmic entropy and ECT entropy over time [7].

Figure 3 .
Figure 3. (a) Monitored parameters during Li-ion battery cycling showing 1.5 h discharge starting at ~5 A, followed by 1.5 h charge at 3 A. (b) phenomenological entropy terms: Ohmic entropy and ECT entropy over time [7].

Figure 4 .
Figure 4. (a) Monitored parameters during steel rod fatigue loading and (b) phenomenological entropy terms: load entropy density and MST entropy density over time [9].

Figure 4 .
Figure 4. (a) Monitored parameters during steel rod fatigue loading and (b) phenomenological entropy terms: load entropy density and MST entropy density over time [9].

4. 5 .
Pump Flow-Pressure and Flow Rate (Internal Energy) Water flowed through a three-phase 15-hp 260-gpm centrifugal motor-pump, instrumented to measure inlet and exit pump pressures and flow rate [26], which are plotted versus time in Figure 5a.A valve adjusted the flow rate.Measured volumetric flow rate was converted to mass flow rate using the density of water ( .m = ρ .V), and pump power Y .X = M T ω, the product of torque M T and rotational speed ω.Substituting into Equation (19), with positive input power, yielded .

Figure 5 .
Figure 5. (a) Monitored inlet/exit pressures and flow rate during a centrifugal motor pump operation, and (b) phenomenological entropy terms: flow entropy and load entropy over time.

Figure 5 .
Figure 5. (a) Monitored inlet/exit pressures and flow rate during a centrifugal motor pump operation, and (b) phenomenological entropy terms: flow entropy and load entropy over time.

Table 1 .
Summary of entropy generation formulations for various system categories, with examples.Here, subscript e, and superscripts r, p and d represent flow, reactants, products and diffusion, respectively.For thermal systems, µ is the molar enthalpy ∂H ∂N .For other (chemical, etc.) reactive/energy systems, µ is the molar Gibbs energy or chemical potential ∂G ∂N .More information on c, β and λ is in Section 2.5 and Appendix A.

Figure 7 .
Figure 7. Illustrations of the Minimum Entropy Generation theorem, showing reversible (green) and minimum (purple) paths; the difference between the paths defines minimum entropy generation (orange).(a) Rates dS rev , δS min , δS ′ and (b) Accumulations S rev , S min , S ′ .
3.Exact differentials dE, state functions of the independent state variables Z(t), if intractable, can be numerically integrated over the (ideal) reversible path per methods of the prior paragraph.The reversible path must transit states o to f in quasiequilibrium and be continuous and maximally smooth over time, which can be approximated by linear functions with slope determined by the end states, for exam- ple, if dE(t) = dE(Z(t)), components of Z(t) with slope = δS phen − dS rev ≥ 0,

the Phenomenological Entropy Generation (PEG) Theorem;
, demonstrated that: • a combination of the thermodynamic potentials and the irreversible form of the TdS equation yields the irreversible Massieu functions; • steady-state systems generate entropy at a minimum rate which, for a boundaryloaded, internally reactive open system, is the sum of boundary work/load entropy S ′ partial derivative of the Helmholtz free energy with respect to generalized position/displacement X defines the load (or strength) modulus E ' of a material.Substituting the second equality of Equation (A6) into the second partial Helmholtz derivative yields the generalized force/strength response to generalized position/displacement, commonly used in analysis of systems under external boundary load, such as stress-loaded (elastic, shear, hardness moduli) and compressible systems (bulk modulus).A comparison of Equations ( second