Topology Optimization of Two Fluid Heat Transfer Problems for Heat Exchanger Design

Topology optimization of thermal-fluid coupling problems has received widespread attention.This article proposes a novel topology optimization method for laminar two-fluid heat exchanger design.The proposed method utilizes an artificial density field to create two permeability interpolation functions that exhibit opposing trends,ensuring separation between the two fluid domains.Additionally,a Gaussian function is employed to construct an interpolation function for the thermal conductivity coefficient.Furthermore,a computational program has been developed on the OpenFOAM platform for the topology optimization of two-fluid heat exchangers.This program leverages parallel computing,significantly reducing the time required for the topology optimization process.To enhance computational speed and reduce the number of constraint conditions,we replaced the conventional pressure drop constraint condition in the optimization problem with a pressure inlet/outlet boundary condition.The 3D optimization results demonstrate the characteristic features of a surface structure,providing valuable guidance for designing heat exchangers that achieve high heat exchange efficiency while minimizing excessive pressure loss.At the same time,a new structure appears in large-scale topology optimization,which proves the effectiveness and stability of the topology optimization program written in this paper in large-scale calculation.

Comparative study of two fluid model and dense discrete phase model for simulations of gas-solid hydrodynamics in circulating fluidized beds

Computational fluid dynamics(CFD)has become a valuable tool to study the complex gas-solid hydrodynamics in the circulating fluidized bed(CFB).Based on the two fluid model(TFM)under the Eulerian-Eulerian framework and the dense discrete phase model(DDPM)under the Eulerian-Lagrangian framework,this work conducts the comparative study of the gas-solid hydrodynamics in a CFB riser by these two different models.Results show that DDPM could be used to predict gas-solid hydrodynamics in the circulating fluidized bed,and there are differences between TFM and DDPM,especially in the radial distribution profiles of solid phase.Sensitivity analysis results show that the gas-solid drag model exhibits significant effects on the results for both the two models.The specularity coefficient and the restitution coefficient in the TFM,as well as the reflection coefficient and the parcel number in the DDPM,exhibit less impact on the simulated results.

Numerical simulation of gas-solid flow with two fluid model in a spouted-fluid bed

The flow characteristics in a spouted-fluid bed differ from those in spouted or fluidized beds because of the injection of the spouting gas and the introduction of a fluidizing gas. The flow behavior of gas-solid phases was predicted using the Eulerian-Eulerian two-fluid model （TFM） approach with kinetic theory for granular flow to obtain the flow patterns in spouted-fluid beds. The gas flux and gas incident angle have a significant influence on the porosity and particle concentration in gas-solid spouted-fluid beds. The fluidizing gas flux affects the flow behavior of particles in the fountain. In the spouted-fluid bed, the solids volume fraction is low in the spout and high in the annulus. However, the solids volume fraction is reduced near the wall.

An Interacting Two-Fluid Scenario for Quintom Dark Energy

The Quintom dark energy is a proposal that explains the recent observations that mildly favor the equation of state of dark energy w crossing -1 near the past. The Quintom model is often constructed by two scalar fields, where one is the quintessence field and another is the phantom field. The cosmological implication of the coupling of the two fields of the dark energy is out of question worth investigating. However, the consideration of the coupling in the field scenario is somewhat complex thus we propose an interacting two-fluid Quintom scenario for simplicity. The interaction between the two components is parametrized by a constant η in this scenario. The cosmological implications of this parametrization are investigated in detail in this paper. Also, a diagnostic for this model is performed by using the statefinder pairs {s, r} and {q, r}.

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 interracial 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.

Validity ranges of interfacial wave theories in a two-layer fluid system

In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness dl, and lower layer thick-ness d2, instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehaute＇s plot for free surface waves if water depth ratio r= d1/d2 approaches to infinity and the upper layer water density p1 to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of σ=（P2 - Pl）/P2 → 1.0 and r 〉 1.0. In the end, several figures of the validity ranges for various interfacial wavetheories in the two-layer fluid are given and compared with the results for surface waves.

Green Functions with Pulsating Sources in A Two-Layer Fluid of Finite Depth

The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.

A Coupled Variable Coefficient Modified KdV Equation Arising from a Two-Layer Fluid System

A quite general coupled variable coefficient modified KdV (VCmKdV) equation in a two-layer fluid systemis derived by means of the reductive perturbation method.Making use of the CK's direct method,some similarityreductions of the coupled VCmKdV equation are obtained and their corresponding group explanations are discussed.Some exact solutions of the coupled equations are also presented.

Diffraction of Water Waves by A Vertically Floating Cylinder in A Two-Layer Fluid

In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder are obtained by use of the method of eigenfunction expansions. The hydrodynamic loads on the vertically floating cylinder in a two-layer fluid inelude not only the surge, heave and pitch exciting forces due to the incident wave of the surface-wave mode, but also those due to the incident wave of the internal-wave mode. This is different from the case of a homogenous fluid. Some given examples show that, for a two-layer fluid system with a small density difference, the hydrodynamic loads for the surface-wave mode do not differ significantly from those due to surface waves in a single-layer fluid, but the hydrodynamic loads for the internal-wave mode are important over a wide range of frequencies. Moreover, also considered are the free surface and interface elevations generated by the diffraction wave due to the incident wave of the surface-wave and interhal-wave modes, and transfer of energy between modes.

Second-order random interfacial wave solutions for two-layer fluid with a free surface

A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song （2004. Geophys Res Lett, 31 （15）：L15302）. The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves （Thorpe. 1968a.Trans R Soc London, 263A：563～614） and standing interracial waves （Thorpe. 1968b. J Fluid Mech, 32：489-528） for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean （1979.Ocean Engineering Rep 20） if the density of the upper layer is much smaller than that of the lower layer.

Hydroelastic Response of A Circular Plate in Waves on A Two-Layer Fluid of Finite Depth

The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.

Oblique Water Wave Scattering by Bottom Undulation in a Two-layer Fluid Flowing Through a Channel

在在由二层组成的隧道流动的小变丑上的繁殖被考虑的倾斜的波浪(内部波浪) 的问题。上面的液体被假定被一个僵硬的盖上面跳，它是为免费表面的近似，并且更低的被有小变丑的透不过的底部表面下面围住；隧道在水平方向是无界的。假定无旋的运动，不安技术被采用由与适当草地功能的介绍合适使用草地不可分的定理在二液体计算速度潜力的一阶的修正。那些功能在以包含代表底部变丑的形状功能 c (x) 的积分计算思考和传播系数帮助。三维的线性水波浪理论为提出相关边界价值问题被利用。底部变丑的二个特殊例子被认为验证结果。正弦曲线涟漪的一个补丁的考虑(让一样的波浪数) 思考系数是比率的摆动的功能的表演两次到涟漪波浪数字的波浪数字的 x 部件。当这比率接近一个时，理论预言在床和接口之间的一个反响的相互作用，并且思考系数变得一涟漪的数字多重。如果这个数字大，事件波浪精力的高思考发生。类似的结果为有不同波浪数字的正弦曲线涟漪的一个补丁被观察。为小入射角，反映的精力与直到 / 的另外的入射角相比是更大的，这也被观察 4。这些理论观察被图形的结果支持。

Influence of the earth rotation on the interfacial wave solutions and wave-induced tangential stress in a two-layer fluid system

Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a fiat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected, the solutions derived from the present work include as special cases those obtained by Sun et al. （2004. Science in China, Ser. D, 47（12）： 1147-1154） for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if the density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth＇s rotation both on the surface wave solutions and the interfacial wave solutions should be considered.

Higher-order Boussinesq-type equations for interfacial waves in a two-fluid system

Interracial waves propagating along the interface between a three-dimensional two-fluid system with a rigid upper boundary and an uneven bottom are considered. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. A set of higher-order Boussinesq-type equations in terms of the depth-averaged velocities accounting for stronger nonlinearity are derived. When the small parameter measuring frequency dispersion keeping up to lower-order and full nonlinearity are considered, the equations include the Choi and Camassa＇s results （1999）. The enhanced equations in terms of the depth-averaged velocities are obtained by applying the enhancement technique introduced by Madsen et al. （1991） and Schaffer and Madsen （1995a）. It is noted that the equations derived from the present study include, as special cases, those obtained by Madsen and Schaffer （1998）. By comparison with the dispersion relation of the linear Stokes waves, we found that the dispersion relation is more improved than Choi and Camassa＇s （1999） results, and the applicable scope of water depth is deeper.

Three-dimensional Computational Fluid Dynamics Modeling of Two-phase Flow in a Structured Packing Column

Characterizing the complex two-phase hydrodynamics in structured packed columns requires a powerful modeling tool.The traditional two-dimensional model exhibits limitations when one attempts to model the detailed two-phase flow inside the columns.The present paper presents a three-dimensional computational fluid dynamics(CFD) model to simulate the two-phase flow in a representative unit of the column.The unit consists of an entire corrugation channel and describes well the real liquid flow conditions.The detailed unsteady two-phase 3D CFD calculations on column packed with Flexipak 1Y were implemented within the volume of fluid(VOF) mathematical framework.The CFD model was validated by comparing the calculated thickness of liquid film with the available experimental data.Special attention was given to quantitative analysis of the effects of gravity on the hydrodynamics.Fluctuations in the liquid mass flow rate and the calculated pressure drop loss were found to be qualitatively in agreement with the experimental observations.

Construction of Wave-free Potentials and Multipoles in a Two-layer Fluid Having Free-surface Boundary Condition with Higher-order Derivatives

有高顺序的边界条件为秒顺序产生的液体结构相互作用的域里的问题的一个大类部分微分方程。各种各样的方法正在被使用处理与 Laplaces 方程联系的这些种混合边界值问题(或 Helmholtz 方程) 在通过固体或液体宣传的波浪的学习产生。在波浪结构相互作用的广泛地使用的方法之一是多极扩大方法。这扩大包含一个常规波浪，波浪来源，波浪偶极子和常规没有波浪的部分的一般联合。没有波浪的部分能进一步以作为没有波浪的潜力被称为的没有波浪的多极被扩展。这些是 Laplaces 方程或二维的 Helmholz 方程的单个答案。这些没有波浪的潜力和多极的建设为象二维的非倾斜、倾斜的波浪那样的很多种状况以一种系统的方式这里被介绍，在有有部分衍生物被考虑的更高的顺序的免费表面状况的二层的液体的三个维的波浪。特别地，这些被获得面对作为一个薄有弹性的盘子建模的一个冰封面在免费表面并且也表面紧张的存在的效果考虑。另外为限制大小写，多极和没有波浪的潜在的功能走向单个层多极和没有波浪的潜力，这能被显示出。

Transient electro-osmotic and pressure driven flows of two-layer fluids through a slit microchannel

By method of the Laplace transform, this arti- cle presents semi-analytical solutions for transient electro- osmotic and pressure-driven flows （EOF/PDF） of two-layer fluids between microparallel plates. The linearized Poisson- Boltzmann equation and the Cauchy momentum equation have been solved in this article. At the interface, the Maxwell stress is included as the boundary condition. By numerical computations of the inverse Laplace transform, the effects of dielectric constant ratio e, density ratio p, pressure ratio p, viscosity ratioμ of layer II to layer I, interface zeta potential difference △ψ, interface charge density jump Q, the ratios of maximum electro-osmotic velocity to pressure velocity , and the normalized pressure gradient B on transient veloc- ity amplitude are presented.We find the velocity amplitude becomes large with the interface zeta potential difference and becomes small with the increase of the viscosity. The ve- locity will be large with the increases of dielectric constant ratio; the density ratio almost does not influence the EOF ve- locity. Larger interface charge density jump leads to a strong jump of velocity at the interface. Additionally, the effects of the thickness of fluid layers （hi and h2） and pressure gradient on the velocity are also investigated.

Fluid Flow Behavior of Liquid in Cylindrical Vessels Stirred by One or Two Air Jets

Based on the two-phase model (Eulerian-Eulerian model), the three dimensional fluid flow in water and that liquid steel systems stirred by one or two multiple gas jets are simulated. In the Eulerian-Eulerian two-phase model, the gas and the liquid phase are considered to be two different continuous fluids interacting with each other through the finite inter-phase areas. The exchange between the phases is represented by source terms in conversation equations. Turbulence is assumed to be a property of the liquid phase. A new turbulence modification - model is introduced to consider the bubbles movement contribution to and . The dispersion of phases due to turbulence is represented by introducing a diffusion term in mass conservation equation. The mathematical simulation agrees well with the experiment results. The study results indicate that the distance of two nozzles has big effect on fluid flow behavior in the vessel. Using two gas injection nozzles at the half radii of one diameter of the bottom generates a much better mixing than with one nozzle under the condition of the same total gas flow rate.

Exact Solution for Pressure Driven Flow of Two Immiscible Phan-Thien-Tanner Fluids in a Pipe

This paper provides the exact solutions for the fully developed two layer pressure driven flows of incompressible Phan-Thien-Tanner fluids in a horizontal cylindrical pipe. Exact equations are formulated and solved for important kinematic properties, such as, velocity profiles, normal and shear stresses, total volume fluxes through a circular cross-section and average velocities. Graphical results are provided and discussed for the different flow parameters. A comparison of Upper Convected Maxwell (UCM), Linear Phan-Thien-Tanner (LPTT) and Exponential Phan-Thien-Tanner (EPTT) shows that UCM is a low viscosity fluid as compared to LPTT, and EPTT and LPTT is lighter than EPTT.

Simulation and Analysis on the Two-Phase Flow Fields in a Rotating-Stream-Tray Absorber by Using Computational Fluid Dynamics

The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The simulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, mad a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall;the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.