Thermodynamics of the Stockmayer fluid in an applied field

The thermodynamic properties of the Stockmayer fluid in an applied field are studied using theory and computer simulation. Theoretical expressions for the second and third virial coefficients are obtained in terms of the dipolar coupling constant (λ, measuring the strength of dipolar interactions as compared to thermal energy) and dipole–field interaction energy (α, being proportional to the applied field strength). These expressions are tested against numerical results obtained by Mayer sampling calculations. The expression for the second virial coefficient contains terms up to λ⁴, and is found to be accurate over realistic ranges of dipole moment and temperature, and over the entire range of the applied field strength (from zero to infinity). The corresponding expression for the third virial coefficient is truncated at λ³, and is not very accurate: higher order terms are very difficult to calculate. The virial coefficients are incorporated in to a thermodynamic theory based on a logarithmic representation of the Helmholtz free energy. This theory is designed to retain the input virial coefficients, and account for some higher order terms in the sense of a resummation. The compressibility factor is obtained from the theory and compared to results from molecular dynamics simulations with a typical value λ = 1. Despite the mathematical approximations of the virial coefficients, the theory captures the effects of the applied field very well. Finally, the vapour–liquid critical parameters are determined from the theory, and compared to published simulation results; the agreement between the theory and simulations is good.     

Thermodynamics of a long-range triangle-well fluid

The long-range triangle-well fluid has been studied using three different approaches: firstly, by an analytical equation of state obtained by a perturbation theory, secondly via a self-consistent integral equation theory, the so-called self-consistent Ornstein–Zernike approach (SCOZA) which is presently one of the most accurate liquid-state theories, and finally by Monte Carlo simulations. We present vapour–liquid phase diagrams and thermodynamic properties such as the internal energy and the pressure as a function of the density at different temperatures and for several values of the potential range. We assess the accuracy of the theoretical approaches by comparison with Monte Carlo simulations: the SCOZA method accurately predicts the thermodynamics of these systems and the first-order perturbation theory reproduces the overall thermodynamic behaviour for ranges greater than two molecular diameters except that it overestimates the critical point. The simplicity of the equation of state and the fact that it is analytical in the potential range makes it a good candidate to be used for calculating other thermodynamic properties and as an ingredient in more complex theoretical approaches.     

The green function of an infinite,fluid loaded membrane

In this paper the response of a fluid loaded plane structure (a membrane) to a concentrated line force excitation is considered in great detail. The normalized velocity response—here called the Green function G—depends upon a dimensionless range x₀=k_m|x|, where k_m is the free wavenumber on the membrane in a vacuum, on the Mach number $\text{M=}\text{k}{}_0\text{k}{}_{\mathrm{m}}$, the ratio of wave phase speed ω/k_m on the unloaded membrane to the sound speed ω/k₀, and on a parameter ? which can be regarded as a measure of fluid loading at the “coincidence” condition M=1. In the analogous problem involving a thin elastic plate, the corresponding parameter is independent of frequency and plate thickness and may be regarded as an intrinsic measure of fluid loading; moreover, in cases of common interest (steel in water, aluminium in air) that parameter is small. In the present paper, the asymptotic structure of G(x₀, M, ?) is therefore sought in the limit ? → 0. Naturally, no single asymptotic expansion can be expected to be valid throughout the (x₀, M) plane, and the programme therefore involves the delineation of regions of that plane in which distinct asymptotic results apply, the construction and discussion of those results, and the asymptotic matching (according to the procedures of the method of matched asymptotic expansions) of results holding in adjoining regions. The Fourier integral for G is broken into surface wave and acoustic components, and the asymptotic structure obtained for each. Previously obtained results for the behaviour at large distances are recovered, with a demonstration that very large distances indeed (x₀ ? ?^?2) may be needed for their validity for some ranges of M; and the drive point behaviour, of G(x₀=0, M, ?) as ? → 0 qua function of M, is shown to correspond to that already discussed in the literature. Elsewhere, in the covering of the whole (x₀, M) plane by different asymptotic expressions, a wide variety of analytical results is found, reflecting the achievement in different regions of different balances among the five competing physical mechanisms represented in the model: namely, structural stiffness, structural inertia, fluid pressures, fluid compressibility and fluid inertia. These different balances give rise to a wide variety of expressions for the phase and amplitude of the surface wave and acoustic components which can now be used to isolate the dominant structural and acoustic mechanisms at any point in the (x₀, M) plane.     

The thermoelectric working fluid: Thermodynamics and transport

Thermoelectric devices are heat engines, which operate as generators or refrigerators using the conduction electrons as a working fluid. The thermoelectric heat-to-work conversion efficiency has always been typically quite low, but much effort continues to be devoted to the design of new materials boasting improved transport properties that would make them of the electron crystal–phonon glass type of systems. On the other hand, there are comparatively few studies where a proper thermodynamic treatment of the electronic working fluid is proposed. The present article aims at contributing to bridge this gap by addressing both the thermodynamic and transport properties of the thermoelectric working fluid covering a variety of models, including interacting systems.     

Thermodynamics of the λ line in binary fluid mixtures

Thermodynamics of the λ line, variations on a theme by Wheeler and Griffiths, have been presented.     

高分子锚定的流体膜

将高分子自洽平均场理论和膜弹性理论结合起来研究高分子锚定流体膜体系，得到了高分子链热力学平衡浓度的分布和膜的形状.由于膜对高分子的不可穿透性，造成高分子可以实现的构象受到限制，导致高分子对膜施加不均匀的熵压，从而锚定点附近的膜弯离高分子，这些结果符合前人的理论分析和蒙特卡罗模拟的结果.此外，考察了膜和高分子链段之间的相互作用，高分子链长以及膜的弯曲刚性、张力对膜形变的影响. 关键词： 高分子 流体膜 自洽场理论 Helfrich膜弹性理论     

Stochastic Thermodynamics of an Electromagnetic Energy Harvester

We study the power extracted by an electromagnetic energy harvester driven by broadband vibrations. We describe the system with a linear model, featuring an underdamped stochastic differential equation for an effective mass in a harmonic potential, coupled electromechanically with the current in the circuit. We compare the characteristic curve (power vs. load resistance) obtained in experiments for several values of the vibration amplitude with the analytical results computed from the model. Then, we focus on a more refined analysis, taking into account the temporal correlations of the current signal and the fluctuations of the extracted power over finite times. We find a very good agreement between the analytical predictions and the experimental data, showing that the linear model with effective parameters can describe the real system, even at the fine level of fluctuations. Our results could be useful in the framework of stochastic thermodynamics applied to energy harvesting systems.     

Third-order Stokes wave solutions for interfacial internalwaves in three-layer dendity-stratified fluid

<正>Interfacial internal waves in a three-layer density-stratified fluid are investigated using a singular perturbation method,and third-order asymptotic solutions of the velocity potentials and third-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory.As expected,the third-order solutions describe the third-order nonlinear modification and the third-order nonlinear interactions between the interfacial waves.The wave velocity depends on not only the wave number and the depth of each layer but also on the wave amplitude.     

Thermodynamics predicts how confinement modifies the dynamics of the equilibrium hard-sphere fluid

We study how confining the equilibrium hard-sphere fluid to restrictive one- and two-dimensional channels with smooth interacting walls modifies its structure, dynamics, and entropy using molecular dynamics and transition-matrix Monte Carlo simulations. Although confinement strongly affects local structuring, the relationships between self-diffusivity, excess entropy, and average fluid density are, to an excellent approximation, independent of channel width or particle-wall interactions. Thus, thermodynamics can be used to predict how confinement impacts dynamics.     

耦合界面张力的三维流体界面不稳定性的格子Boltzmann模拟

采用介观格子Boltzmann方法模拟界面张力作用下三维流体界面的Rayleigh-Taylor (RT)不稳定性的增长过程,主要分析表面张力对流体界面动力学行为及尖钉和气泡后期增长的影响机制.首先发现三维RT不稳定性的发生存在临界表面张力(σ_c),其值随着流体Atwood数的增大而增大,且数值预测值与理论分析结果σ_c=(ρ_h-ρ₁)g/k~2一致.另外,随着表面张力的增大,不稳定性演化过程中界面卷吸程度和结构复杂性逐渐减弱,系统中界面破裂形成离散液滴的数目也显著减少.相界面的后期动力学行为也从非对称发展转向始终保持关于中轴线对称.尖钉与气泡振幅在表面张力较小时对其变化不显著,当表面张力增大到一定值后,可以有效地抑制尖钉与气泡振幅的增长.进一步发现,高雷诺数三维RT不稳定性在不同表面张力下均经历4个不同的发展阶段:线性阶段、饱和速度阶段、重加速和混沌混合阶段.尖钉与气泡在饱和速度阶段以近似恒定的速度增长,其渐进速度的值与修正的势流理论模型结果一致.受非线性Kelvin-Helmholtz旋涡的剪切作...     

A hybrid numerical method for interfacial fluid flow with soluble surfactant

We address a significant difficulty in the numerical computation of fluid interfaces with soluble surfactant that occurs in the physically representative limit of large bulk Peclet number Pe. At the high values of Pe in typical fluid-surfactant systems, there is a transition layer near the interface in which the surfactant concentration varies rapidly, and large gradients at the interface must be resolved accurately to evaluate the exchange of surfactant between the interface and bulk flow. We use the slenderness of the layer to develop a fast and accurate ‘hybrid’ numerical method that incorporates a separate, singular perturbation analysis of the dynamics in the transition layer into a full numerical solution of the interfacial free boundary problem. The accuracy and efficiency of the method is assessed by comparison with a more ‘traditional’ numerical approach that uses finite differences on a curvilinear coordinate system exterior to the bubble, without the separate transition layer reduction. The traditional method implemented here features a novel fast calculation of fluid velocity off the interface.     

Stability of fluid flow past a membrane

The stability of the flow of a fluid past a solid membrane of infinitesimal thickness is investigated using a linear stability analysis. The system consists of two fluids of thicknesses R and H R and bounded by rigid walls moving with velocities and , and separated by a membrane of infinitesimal thickness which is flat in the unperturbed state. The fluids are described by the Navier-Stokes equations, while the constitutive equation for the membrane incorporates the surface tension, and the effect of curvature elasticity is also examined for a membrane with no surface tension. The stability of the system depends on the dimensionless strain rates and in the two fluids, which are defined as and for a membrane with surface tension , and and for a membrane with zero surface tension and curvature elasticity K. In the absence of fluid inertia, the perturbations are always stable. In the limit , the decay rate of the perturbations is O(k ³ ) smaller than the frequency of the fluctuations. The effect of fluid inertia in this limit is incorporated using a small wave number asymptotic analysis, and it is found that there is a correction of smaller than the leading order frequency due to inertial effects. This correction causes long wave fluctuations to be unstable for certain values of the ratio of strain rates and ratio of thicknesses H. The stability of the system at finite Reynolds number was calculated using numerical techniques for the case where the strain rate in one of the fluids is zero. The stability depends on the Reynolds number for the fluid with the non-zero strain rate, and the parameter , where is the surface tension of the membrane. It is found that the Reynolds number for the transition from stable to unstable modes, , first increases with , undergoes a turning point and a further increase in the results in a decrease in . This indicates that there are unstable perturbations only in a finite domain in the plane, and perturbations are always stable outside this domain. Received: 29 May 1997 / Revised: 9 October 1997 / Accepted: 26 November 1997     

Critical fluctuations of tense fluid membrane tubules

We show that, contrary to planar membranes under tension, tense membrane tubules exhibit important critical fluctuations originating from unidimensional Goldstone modes. The latter yield unexpected behavior, such as correlations extending over the whole tube length and the increase of the fluctuating area over the projected area with increasing tension.     

Visualization of transient interfacial waves induced by spin-up of two immiscible fluid layers

Abstract  

Interfacial waves of two immiscible layers in a spin-up container were investigated using experimental visualization. While the interface near the central part rose up, instability waves propagated in an azimuthal direction on the interface. These waves were mainly caused by Kelvin–Helmholtz instability for the velocity difference between two layers during spin-up, but had complicated transient characteristics owing to the rotation in a closed system. We visualized the structure of the interfacial waves by the use of three types of optical characteristics of the interface. Image processing provided the detailed factors of the interfacial waves that were classified in four life stages from their generation to disappearance. The initial generation process involved many frequency modes due to a large velocity difference, and then a low mode stood out during the growth, and disappeared with an ellipsoidal sloshing mode to achieve the rigid rotation in both layers.     

A set of Boussinesq-type equations for interfacial internal waves in two-layer stratified fluid

Many new forms of Boussinesq-type equations have been developed to extend the range of applicability of the classical Boussinesq equations to deeper water in the study of the surface waves. One approach was used by Nwogu (1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618--638) to improve the linear dispersion characteristics of the classical Boussinesq equations by using the velocity at an arbitrary level as the velocity variable in derived equations and obtain a new form of Boussinesq-type equations, in which the dispersion property can be optimized by choosing the velocity variable at an adequate level. In this paper, a set of Boussinesq-type equations describing the motions of the interfacial waves propagating alone the interface between two homogeneous incompressible and inviscid fluids of different densities with a free surface and a variable water depth were derived using a method similar to that used by Nwogu (1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618--638) for surface waves. The equations were expressed in terms of the displacements of free surface and density-interface, and the velocity vectors at arbitrary vertical locations in the upper layer and the lower layer (or depth-averaged velocity vector across each layer) of a two-layer fluid. As expected, the equations derived in the present work include as special cases those obtained by Nwogu (1993, J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638) and Peregrine (1967, J. Fluid Mech. 27, 815-827) for surface waves when the density of the upper fluid is taken as zero.