Dimensionless Compressibility Research Articles

Theoretical analysis of the structure, thermodynamics, and shear elasticity of deeply metastable hard sphere fluids.

The structure, thermodynamics, and slow activated dynamics of the equilibrated metastable regime of glass-forming fluids remain a poorly understood problem of high theoretical and experimental interest. We apply a highly accurate microscopic equilibrium liquid state integral equation theory, in conjunction with naïve mode coupling theory of particle localization, to study in a unified manner the structural correlations, thermodynamic properties, and dynamic elastic shear modulus in deeply metastable hard sphere fluids. Distinctive behaviors are predicted including divergent inverse critical power laws for the contact value of the pair correlation function, pressure, and inverse dimensionless compressibility, and a splitting of the second peak and large suppression of interstitial configurations of the pair correlation function. The dynamic elastic modulus is predicted to exhibit two distinct exponential growth regimes with packing fraction that have strongly different slopes. These thermodynamic, structural, and elastic modulus results are consistent with simulations and experiments. Perhaps most unexpectedly, connections between the amplitude of long wavelength density fluctuations, dimensionless compressibility, local structure, and the dynamic elastic shear modulus have been theoretically elucidated. These connections are more broadly relevant to understanding the slow activated relaxation and mechanical response of colloidal suspensions in the ultradense metastable region and deeply supercooled thermal liquids in equilibrium.

Compressibility Effect on Darcy Porous Convection

Perfectly incompressible materials do not exist in nature but are a useful approximation of several media which can be deformed in non-isothermal processes but undergo very small volume variations. In this paper, the linear analysis of the Darcy-Bénard problem is performed in the class of extended-quasi-thermal-incompressible fluids, introducing a factor beta which describes the compressibility of the fluid and plays an essential role in the instability results. In particular, in the Oberbeck-Boussinesq approximation, a more realistic constitutive equation for the fluid density is employed in order to obtain more thermodynamically consistent instability results. The critical Rayleigh-Darcy number for the onset of convection is determined, via linear instability analysis of the conduction solution, as a function of a dimensionless parameter widehat{beta } proportional to the compressibility factor beta, proving that widehat{beta } enhances the onset of convective motions.Article HighlightsThe onset of convection in fluid-saturated porous media is analyzed, taking into account fluid compressibility effect.The critical Rayleigh-Darcy number is determined in a closed algebraic form via linear instability analysis.The critical Rayleigh-Darcy number is shown to be a decreasing function of the dimensionless compressibility factor.

Dynamics of Compression-Driven Gas-Liquid Displacement in a Capillary Tube.

We study two-phase displacement via the steady compression of an air reservoir connected to an oil-filled capillary tube. Our experiments and modeling reveal complex displacement dynamics depending on compression rate and reservoir volume that, for large reservoirs, depend on a single dimensionless compressibility number. We identify two distinct displacement regimes, separated by a critical value of the compressibility number. While the subcritical regime exhibits quasisteady displacement after an initial transient, the supercritical regime exhibits burstlike expulsion.

Understanding dynamics in coarse-grained models. II. Coarse-grained diffusion modeled using hard sphere theory.

The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S0(T) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S0(T) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S0(T) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Long Wavelength Thermal Density Fluctuations in Molecular and Polymer Glass-Forming Liquids: Experimental and Theoretical Analysis under Isobaric Conditions.

We establish via an in-depth analysis of experimental data that the dimensionless compressibility (proportional to the dimensionless amplitude of long wavelength thermal density fluctuations) of one-component normal and supercooled liquids of chemically complex nonpolar and weakly polar molecules and polymers follows extremely well a surprisingly simple and general temperature dependence over an exceptionally wide range of pressures and temperatures. A theoretical basis for this behavior is shown to exist in the venerable van der Waals model and its more modern interpretations. Although associated hydrogen-bonding (and to a lesser degree strongly polar) liquids display modestly more complex behavior, rather simple temperature and pressure dependences are also discovered. A new approach to collapse the temperature- and pressure-dependent dimensionless compressibility data onto a master curve is formulated that differs from the empirical thermodynamic scaling approach. As a practical matter, we also find that the dimensionless compressibility scales well as an inverse power law with temperature with an exponent that is system dependent and decreases with pressure. At very high pressures and low temperatures, the thermal liquid behavior appears to approach (but not reach) a repulsion-dominated random close packing limit. All these findings are relevant to our recent theoretical work on the problem of activated relaxation and vitrification of supercooled molecular and polymeric liquids.

Theory of Transient Localization, Activated Dynamics, and a Macromolecular Glass Transition in Ring Polymer Liquids.

We construct a segmental scale force level theory for the center-of-mass diffusion constant and corresponding relaxation time for globally compact unconcatenated ring polymer solutions and melts (degree of polymerization N). The approach is based on slowly decaying macromolecular scale intermolecular force dynamic correlations as the origin of their unusual dynamics. Unentangled Rouse, weakly caged, and activated regimes are predicted. The barrier of the activated regime scales linearly with N and as a power law of concentration, which drives a kinetic glass transition on the radius-of-gyration scale. The values of N at the two dynamic crossovers (Rouse to weakly caged, weakly caged to activated) are proportional, with nonuniversality entering mainly via macromolecular volume fraction and dimensionless compressibility. Quantitative comparisons with simulation data reveal good agreement. Aspects of intermediate time dynamics are analyzed, and predictions are made for the conditions required to observe a macromolecular glass transition in the laboratory and on the computer.

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses.

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical "slaving onset" stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

Experimental test of a predicted dynamics–structure–thermodynamics connection in molecularly complex glass-forming liquids

Understanding in a unified manner the generic and chemically specific aspects of activated dynamics in diverse glass-forming liquids over 14 or more decades in time is a grand challenge in condensed matter physics, physical chemistry, and materials science and engineering. Large families of conceptually distinct models have postulated a causal connection with qualitatively different "order parameters" including various measures of structure, free volume, thermodynamic properties, short or intermediate time dynamics, and mechanical properties. Construction of a predictive theory that covers both the noncooperative and cooperative activated relaxation regimes remains elusive. Here, we test using solely experimental data a recent microscopic dynamical theory prediction that although activated relaxation is a spatially coupled local-nonlocal event with barriers quantified by local pair structure, it can also be understood based on the dimensionless compressibility via an equilibrium statistical mechanics connection between thermodynamics and structure. This prediction is found to be consistent with observations on diverse fragile molecular liquids under isobaric and isochoric conditions and provides a different conceptual view of the global relaxation map. As a corollary, a theoretical basis is established for the structural relaxation time scale growing exponentially with inverse temperature to a high power, consistent with experiments in the deeply supercooled regime. A criterion for the irrelevance of collective elasticity effects is deduced and shown to be consistent with viscous flow in low-fragility inorganic network-forming melts. Finally, implications for relaxation in the equilibrated deep glass state are briefly considered.

Parametric studies on the behaviour of reinforced elastomeric bearings considering the compressibility and extensibility index, shape factor, and lift-off

Reinforced elastomeric bearings (REBs) are widely used in bridges as well as for seismic isolation applications. Analytical solutions for the critical design properties of REBs (e.g. compression and bending modulus, stress distributions) have been derived based on the assumptions of the pressure solution. In all cases, the analytical solutions can be expressed as functions of the dimensionless compressibility and extensibility index. This parameter accounts for the influence of the compressibility of the elastomer and the extensibility of the reinforcement. Together with the shape factor, these variables are highly influential on REB design. The objective of this study is to investigate mechanical behaviour of bonded and unbonded REBs with the same compressibility and extensibility index but different other parameters and the influence of the shape factor by finite element analysis (FEA) including lift-off. For unbonded REBs, lift-off (i.e. the separation between the superstructure and the bearings) can occur under certain combinations of large rotations and relatively low axial stress. Three different loading cases (pure compression, pure bending and the combined case) are considered. The results show that REBs with the same compressibility and extensibility index but different other parameters have similar behaviour under the same loading condition. Available analytical solutions have good agreement with numerical results provided an appropriate shear modulus is selected and the level of rotation is reasonable for structural applications.

Thermodynamics-Structure-Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids.

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus-Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics-structure-dynamics "correlations" based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. The predicted thermodynamics-structure-dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

Integral equation theory of thermodynamics, pair structure, and growing static length scale in metastable hard sphere and Weeks-Chandler-Andersen fluids.

We employ the Ornstein-Zernike integral equation theory with the Percus-Yevick (PY) and modified-Verlet (MV) closures to study the equilibrium structural and thermodynamic properties of metastable monodisperse hard sphere and continuous repulsion Weeks-Chandler-Andersen (WCA) fluids under density and temperature conditions where the system is strongly overcompressed or supercooled, respectively. The theoretical results are compared to crystal-avoiding simulations of these dense monodisperse model one-component fluids. The equation of state (EOS) and dimensionless compressibility are computed using both the virial and compressibility routes. For hard spheres, the MV-based virial route EOS and dimensionless compressibility are in very good agreement with simulation for all packing fractions, much better than the PY analogs. The corresponding MV-based predictions for the static structure factor are also very good. The amplitude of density fluctuations on the local cage scale and in the long wavelength limit, and three technically different measures of the density correlation length, are studied with both closures. All five properties grow in a roughly exponential manner with density in the metastable regime up to packing fractions of 58% with no sign of saturation. The MV-based results are in good agreement with our crystal-avoiding simulations. Interestingly, the density dependences of long and short wavelength quantities are closely related. The MV-based theory is also quite accurate for the thermodynamics and structure of supercooled monodisperse WCA fluids. Overall our findings are also relevant as critical input to microscopic theories that relate the equilibrium pair correlation function or static structure factor to dynamical constraints, barriers, and activated relaxation in glass-forming liquids.

Intermolecular structural correlations in model globular and unconcatenated ring polymer liquids.

We employ the field theoretic polymer integral equation theory to construct a segment-level theory for the thermodynamics and pair structure of dense liquids of interpenetrating ring polymers and a simple globule model. The latter is defined by a fractal mass distribution on all internal length scales with an exponent equal to the spatial dimension (dF = ds = 3). In an isochoric ensemble the dimensionless compressibility and pressure is predicted to vary exponentially with macromolecular volume fraction. An intermolecular correlation hole exists down to small length scales. This model appears to be useful for a recently studied experimental soft nanoparticle suspension, and also serves as a reference system for our analysis of ring liquids. Motivated by simulations, a two-fractal exponent ring model is adopted for the intramolecular structure factor. At smaller lengths it describes chain-like macromolecules, while on larger scales it corresponds to a space-filling object in the sense that dF = ds = 3. The crossover between these two regimes is of order the entanglement length of the linear chain analog. Based on a constant compressibility ensemble, the effective volume fraction grows at intermediate values of degree of polymerization (N), and crosses over to a very slow logarithmic growth at large N. A weaker intermolecular correlation hole is predicted. The number of nearest neighbor rings increases dramatically at small N, akin to linear chain melts, but then tends to saturate at large N, in accord with simulations. The tools developed may be relevant for other partially interpenetrating soft objects such as core-shell nanogels or microgels.

New analytical solutions for weakly compressible Newtonian Poiseuille flows with pressure-dependent viscosity

Steady-state, isothermal, Poiseuille flows in straight channels and circular tubes of weakly compressible Newtonian fluids are considered. The major assumption is that both the mass density and the shear viscosity of the fluid vary linearly with pressure. The non-zero velocity components, the pressure, the mass density and viscosity of the fluid are represented over the flow domain as asymptotic expansions in which the dimensionless isothermal compressibility coefficient ɛ is taken as small parameter. A perturbation analysis is performed and asymptotic solutions for all variables are obtained up to first order in ɛ. The derived solutions, which hold for not necessarily small values of the dimensionless pressure-dependence coefficient, extend previous regular perturbation results and analytical works in the literature for weakly compressible fluids with constant viscosity (solved with a regular perturbation scheme), for incompressible flows with pressure-dependent viscosity (solved analytically), as well as for compressible fluids with pressure-dependent viscosity (solved with double regular perturbation schemes). In contrast to the previous analytical studies in the literature, a non-zero wall-normal velocity is predicted at first order in ε, even at zero Reynolds number. A severe reduction of the volumetric flow-rate at the entrance of the tube/channel and multiplicity of solutions in the flow curves (volumetric flow-rate versus pressure drop) are also predicted. Last, it is shown that weak compressibility of the fluid and the viscosity pressure-dependence have competing effects on the mean friction factor and the average pressure difference required to drive the flow.

Multi-scale entropic depletion phenomena in polymer liquids.

We apply numerical polymer integral equation theory to study the entropic depletion problem for hard spheres dissolved in flexible chain polymer melts and concentrated solutions over an exceptionally wide range of polymer radius of gyration to particle diameter ratios (Rg/D), particle-monomer diameter ratios (D/d), and chain lengths (N) including the monomer and oligomer regimes. Calculations are performed based on a calibration of the effective melt packing fraction that reproduces the isobaric dimensionless isothermal compressibility of real polymer liquids. Three regimes of the polymer-mediated interparticle potential of mean force (PMF) are identified and analyzed in depth. (i) The magnitude of the contact attraction that dominates thermodynamic stability scales linearly with D/d and exhibits a monotonic and nonperturbative logarithmic increase with N ultimately saturating in the long chain limit. (ii) A close to contact repulsive barrier emerges that grows linearly with D/d and can attain values far in excess of thermal energy for experimentally relevant particle sizes and chain lengths. This raises the possibility of kinetic stabilization of particles in nanocomposites. The barrier grows initially logarithmically with N, attains a maximum when 2Rg ∼ D/2, and then decreases towards its asymptotic long chain limit as 2Rg ≫ D. (iii) A long range (of order Rg) repulsive, exponentially decaying component of the depletion potential emerges when polymer coils are smaller than, or of order, the nanoparticle diameter. Its amplitude is effectively constant for 2Rg ≤ D. As the polymer becomes larger than the particle, the amplitude of this feature decreases extremely rapidly and becomes negligible. A weak long range and N-dependent component of the monomer-particle pair correlation function is found which is suggested to be the origin of the long range repulsive PMF. Implications of our results for thermodynamics and miscibility are discussed.

Elastically cooperative activated barrier hopping theory of relaxation in viscous fluids. II. Thermal liquids.

Building on the elastically collective nonlinear Langevin equation theory developed for hard spheres in Paper I, we propose and implement a quasi-universal theory for the alpha relaxation of thermal liquids based on mapping them to an effective hard sphere fluid via the dimensionless compressibility. The result is a zero adjustable parameter theory that can quantitatively address in a unified manner the alpha relaxation time over 14 or more decades. The theory has no singularities above zero Kelvin, and relaxation in the equilibrium low temperature limit is predicted to be of a roughly Arrhenius form. The two-barrier (local cage and long range collective elastic) description results in a rich dynamic behavior including apparent Arrhenius, narrow crossover, and deeply supercooled regimes, and multiple characteristic or crossover times and temperatures of clear physical meaning. Application of the theory to nonpolar molecules, alcohols, rare gases, and liquids metals is carried out. Overall, the agreement with experiment is quite good for the temperature dependence of the alpha time, plateau shear modulus, and Boson-like peak frequency for van der Waals liquids, though less so for hydrogen-bonding molecules. The theory predicts multiple growing length scales upon cooling, which reflect distinct aspects of the coupled local hopping and cooperative elastic physics. Calculations of the growth with cooling of an activation volume, which is strongly correlated with a measure of dynamic cooperativity, agree quantitatively with experiment. Comparisons with elastic, entropy crisis, dynamic facilitation, and other approaches are performed, and a fundamental basis for empirically extracted crossover temperatures is established. The present work sets the stage for addressing distinctive glassy phenomena in polymer melts, and diverse liquids under strong confinement.

Compressible Poiseuille flows with exponential type pressure-dependent mass density

The effect of fluid compressibility in steady, Newtonian and viscoelastic, isothermal, laminar flows is studied. The mass density of the fluid is assumed to vary exponentially with pressure, and viscoelasticity is taken into account by employing the Oldroyd-B constitutive model. The flow problems are solved using a regular perturbation scheme in terms of the dimensionless isothermal compressibility parameter, ε. The sequence of partial differential equations resulting from the perturbation procedure is solved analytically up to third-order in ε. So far, this is the highest order analytical solution derived for Newtonian and viscoelastic compressible flows. Special emphasis on the validity and convergence of the perturbation series is given. It is shown that the inclusion of the third-order term in the perturbation series has important consequences for the preservation of the physical properties of the solution, which is composed up to third order, as well as its convergence. Furthermore, the high accuracy of the perturbation solution in conjunction with the requirement for a positive definite conformation tensor results in monotonic flow curves, as opposed to non-monotonic ones derived in previous works. Cross-sectional averaged expressions for the pressure drop, volumetric flow rate, skin friction factor, and the first- and second-normal stress differences are derived and discussed. For the Newtonian case, the theoretical results for the pressure drop versus the Reynolds number follow the same trend with experimental data from the literature for a much wider range of the dimensionless numbers and parameters to those reported in previous works, which is mainly due to the fact that a more realistic (exponentially, pressure-dependent) equation of state has been utilized.

Strictly two-dimensional self-avoiding walks: Thermodynamic properties revisited

The density crossover scaling of various thermodynamic properties of solutions and melts of self-avoiding and highly flexible polymer chains without chain intersections confined to strictly two dimensions is investigated by means of molecular dynamics and Monte Carlo simulations of a standard coarse-grained bead-spring model. In the semidilute regime we confirm over an order of magnitude of the monomer density ρ the expected power law scaling for the interaction energy between different chains e(int) ~ ρ(21/8), the total pressure P ~ ρ(3) and the dimensionless compressibility g(T) = lim(q→0)S(q) ~ 1/ρ(2). Various elastic contributions associated to the affine and non-affine response to an infinitesimal strain are analyzed as functions of density and sampling time. We show how the size ξ(ρ) of the semidilute blob may be determined experimentally from the total monomer structure factor S(q) characterizing the compressibility of the solution at a given wave vector q. We comment briefly on finite persistence length effects.

Perturbation solutions of weakly compressible Newtonian Poiseuille flows with Navier slip at the wall

We consider both the planar and axisymmetric steady, laminar Poiseuille flows of a weakly compressible Newtonian fluid assuming that slip occurs along the wall following Navier’s slip equation and that the density obeys a linear equation of state. A perturbation analysis is performed in terms of the primary flow variables using the dimensionless isothermal compressibility as the perturbation parameter. Solutions up to the second order are derived and compared with available analytical results. The combined effects of slip, compressibility, and inertia are discussed with emphasis on the required pressure drop and the average Darcy friction factor.

Laminar axisymmetric flow of a weakly compressible viscoelastic fluid

The combined effects of weak compressibility and viscoelasticity in steady, isothermal, laminar axisymmetric Poiseuille flow are investigated. Viscoelasticity is taken into account by employing the Oldroyd-B constitutive model. The fluid is assumed to be weakly compressible with a density that varies linearly with pressure. The flow problem is solved using a regular perturbation scheme in terms of the dimensionless isothermal compressibility parameter. The sequence of partial differential equations resulting from the perturbation procedure is solved analytically up to second order. The two-dimensional solution reveals the effects of compressibility and the other dimensionless numbers and parameters in the flow. Expressions for the average pressure drop, the volumetric flow rate, the total axial stress, as well as for the skin friction factor are also derived and discussed. The validity of other techniques used to obtain approximate solutions of weakly compressible flows is also discussed in conjunction with the present results.