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.

Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.

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.

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.

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.

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.

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.

Two-fluid oscillatory flow in a channel

The validity of Navier’s partial slip condition is investigated by studying the oscillatory flow in a coated channel.The two-fluid model is used to solve the unsteady viscous equations exactly.Partial slip is experienced by the core fluid.It is found that Naviers condition does not hold for an unsteady core fluid.

A kind of extended Korteweg-de Vries equation and solitary wave solutions for interfacial waves in a two-fluid system

This paper considers interfacial waves propagating along the interface between a two-dimensional two-fluid with a flat bottom and a rigid upper boundary. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. It just focuses on the weakly non-linear small amplitude waves by introducing two small independent parameters： the nonlinearity ratio ε, represented by the ratio of amplitude to depth, and the dispersion ratio μ, represented by the square of the ratio of depth to wave length, which quantify the relative importance of nonlinearity and dispersion. It derives an extended KdV equation of the interfacial waves using the method adopted by Dullin et al in the study of the surface waves when considering the order up to O（μ^2）. As expected, the equation derived from the present work includes, as special cases, those obtained by Dullin et al for surface waves when the surface tension is neglected. The equation derived using an alternative method here is the same as the equation presented by Choi and Camassa. Also it solves the equation by borrowing the method presented by Marchant used for surface waves, and obtains its asymptotic solitary wave solutions when the weakly nonlinear and weakly dispersive terms are balanced in the extended KdV equation.

Optimizing control of a two-span rotor system with magnetorheological fluid dampers

A control system aims at vibration reduction in a two-span rotor system with two shear mode magnetorheological （MRF） dampers is designed. A finite element model of the MRF damper- rotor system is built and used to analyze the rotor vibration characteristics. Based on Hooke and Jeeves algorithm and the numerical simulation analysis, an optimal appropriate controller is proposed and designed. Experimental results show that rotor vibration caused by unbalance is well controlled （ first critical speed region 37% , second critical speed region 42% ）. To reflect advantages of optimi- zing strategy presented and validate the intelligent optimization control technology, detailed experi- ments were developed on a two-span rotor-vibration-control platform. The influence on accuracy, rapidity and stability of optimizing control for rotor vibration are analyzed. It provides a powerful technical support for the extension and application in target and control for shafting vibration.

MATHEMATICAL MODEL OF TWO-PHASE FLUID NONLINEAR FLOW IN LOW-PERMEABILITY POROUS MEDIA WITH APPLICATIONS

A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.

Linearization and analytic solution of fluid dynamics of cell two-phase flow

Linearized equations of fluid dynamics of cell two phase flow for one dimensional case are proposed. Based on the equations, an analytic solution is derived, in which the frequency of wave is observed. The frequency formula consists of all important parameters of the fluid dynamics. In our observation, the group velocity and phase velocity of the motion of wave propagation are explicitly exhibited as well.

Interaction of Water Waves with Small Undulations on a Porous Bed in a Two-layer Ice-covered Fluid

The scattering problem involving water waves by small undulation on the porous ocean-bed in a two-layer fluid,is investigated within the framework of the two-dimensional linear water wave theory where the upper layer is covered by a thin uniform sheet of ice modeled as a thin elastic plate.In such a two-layer fluid there exist waves with two different modes,one with a lower wave number propagate along the ice-cover whilst those with a higher wave number propagate along the interface.An incident wave of a particular wave number gets reflected and transmitted over the bottom undulation into waves of both modes.Perturbation analysis in conjunction with the Fourier transform technique is used to derive the first-order corrections of reflection and transmission coefficients for both the modes due to incident waves of two different modes.One special type of bottom topography is considered as an example to evaluate the related coefficients in detail.These coefficients are depicted in graphical forms to demonstrate the transformation of wave energy between the two modes and also to illustrate the effects of the ice sheet and the porosity of the undulating bed.