纳米流体强化传热问题的高效算法设计及其高精度数值模拟

对周期性温度边界条件下方腔内Copper-水纳米流体的对流传热问题建立了完全四阶高精度紧致差分格式,并探讨了该格式高效算法设计,同时引入伪时间导数项采用Von Neumann方法证明了格式的绝对稳定性.研究了不同Ri数下纳米颗粒体积分数、左右边界温度的正弦振幅比和相位偏差对纳米流体传热效率的影响.数值结果表明,对于给定的Ri数,随着纳米颗粒体积分数的增大和边界温度的正弦振幅比的增加,纳米流体的传热效果会显著增强,而相位偏差只影响方腔右侧的温度分布.     

A Mixed Analytical/Numerical Method for Velocity and Heat Transfer of Laminar Power-Law Fluids

This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudo-similarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer.     

Numerical Computation of Heat and Mass Transfer in Fluidized Beds with Liquid Injection

The goal of this work is the simulation of concentration and temperature distributions insides fluidized bed with a uniform liquid distribution. Further, a physically based 2D model is developed for the heat and mass transfer processes in fluidized beds with a spray nozzle. The model is a coupled and semi-linear system of convection-reaction-diffusion equations. We considered the numerical solution of these semi-linear partial differential equations with discrete boundary conditions using linear finite elements on an adaptive triangular grid in space and implicit methods in time. We present calculations using, semi implicit and implicit methods in time, and different solvers for solving the linear systems. The complex correlations of mass and liquid flow rates, mass transfer, heat transfer, drying, and transient two dimensional air humidity, air temperature, degree of wetness, liquid film temperature and particle temperature were simulated. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Numerical Treatment of One-dimensional Models of Radiative Heat Transfer in Co-current Tube Furnaces

A model of a top-fired furnace, such as is used in chemicalindustrial plants, is described. The model gives rise to a systemof first-order differential equations, with boundary conditionsat both the top and bottom of the furnace. This also involvesa non-linear coupling of the radiation field to the energy ofthe gases flowing. Published accounts of this problem make useof shooting methods which converge extremely slowly. An alternativemethod will be described in which the radiative field equationsare solved by an algorithm of Grant and, to deal with the non-linearitiesof the system, a generalized Newton-Raphson method is employed.     

Inverse Problem for Equations of Complex Heat Transfer

Computational Mathematics and Mathematical Physics - The inverse problem with integral overdetermination for the equations of complex heat transfer, including the $${{P}_{1}}$$ approximation for...     

细菌堆浸中流动反应传热耦合过程数值模拟分析

针对辉铜矿微生物浸出过程进行数学建模,利用ComsolMultiphysics软件对浸堆中热量传递、氧气流动、目的金属离子分布以及氧化转化率等进行数值模拟分析．结果表明:沿堆体斜坡边界处氧的浓度较大,而在堆的中央部分,氧气的浓度非常小,导致这部分区域浸出反应缓慢;沿浸堆边坡处的温度偏低,温度最高的部分在底部区域附近靠近边坡位置处,且温度升高值超过6℃;矿堆底部靠近边坡区域目的金属离子浓度最高;靠近底部和斜坡部分,浸出反应速率快、氧化转化率高．矿堆的其余部分,因氧浓度低,导致浸出反应速率缓慢、氧化转化率低．     

Mesoscopic Numerical Study on Flow Boiling Heat Transfer Performance in Channels With Multiple Rectangular Heaters北大核心CSCD

The flow boiling phenomenon in a channel with multiple rectangular heaters under a constant wall temperature was numerically studied with the lattice Boltzmann method. The effects of spacings between heaters, heater lengths and heater surface wettabilities on the bubble morphology, the bubble area and the heat flux on the heater surface, were studied. The results show that, the bubble growth rate increases with the spacing between heaters. The larger the bubble area is, the earlier the nucleated bubbles will leave the heater surface. The corresponding boiling heat transfer performance increases by 12% with the spacing between heaters growing from 250 lattices to 1 000 lattices. On the other hand, the longer the heater length is, the earlier the bubble will nucleate and leave the heater surface, and the better the boiling heat transfer performance will be. The boiling heat transfer performance increases by 13% with the heater length rising from 16 lattices to 22 lattices. In addition, the bubble nucleates later on the hydrophilic surface than on the hydrophobic surface. Compared with the hydrophilic surface, the hydrophobic surface retains residual bubbles after the leaving of bubbles from the heater. The average heat flux and the bubble area of the hydrophilic surface are less than those of the hydrophobic surface. With the contact angle changing from 77° to 120°, the heat transfer performance increases by 26%. Finally, the orthogonal test results indicate that, the wettability of the heat exchanger surface has the greatest influence on the flow boiling heat transfer performance, while the heater length has the least influence. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

矩形微通道热沉内单相稳态层流流体的流动与传热分析

采用解析方法分析了矩形微通道热沉内单相稳态层流流体的流动与传热．基于y方向流速和导热不变的假设,建立流体在矩形微通道内流动的流速方程和传热的温度方程,进而推导出Nusselt数和Poiseuille数的理论表达式．通过计算结果可以看出,推导的Nusselt数和Poiseuille数的解析解与其他文献的结果吻合较好,而且当宽高比趋于无穷大时,Nusselt数和Poiseuille数分别趋近于8.235和96,这与其他文献结果完全相同．在Reynolds数相同时,摩擦因数随着宽高比的增加而增加,而在相同宽高比时,摩擦因数随Reynolds数的增加而减小．     

带有球凹/球凸的旋转通道内流动结构和强化换热的数值研究

对带有凹坑和凸包的内流通道在不同旋转数下的对流换热特性进行了数值分析,探讨了Coriolis力对通道中流场和换热特征的影响．研究发现,随着旋转数增加,通道前缘呈现出较弱的流动冲击,但存在较大的尾迹和延迟的流动再附着,后缘凹坑内部有一较小旋涡和较强射流使得后缘传热得到强化,最高可达60%．总体Nusselt数随着旋转数的增加先减小而后增大.     

电渗流中传热传质过程与熵的分析

在壁面存在恒定热通量条件下,分析微通道内电渗流中传热传质过程与熵的生成.建立数值计算模型,分别采用Poisson-Boltzmann方程、Navier-Stokes方程、Nernst-Planck方程和能量方程来描述微通道内双电层电势、流场、离子浓度和温度的分布情况.引入熵产生,进一步研究不同流动参数对流体传热过程的作用,讨论不同流动参数下各热效应的变化规律,并具体分析热效应参数对流体总熵增加及各部分热效应对总熵比重的影响.结果表明,动电参数与Joule(焦耳)热系数的增大会使得传热性能减弱,动电参数对传热性能影响更为明显;流体的总熵为动电参数、传质系数和质量弥散系数的增函数.     

Numerical Approximation of First Eigenvalues for a Heat Conduction Problem

An eigenvalue problem for a differential operator connected with the heat conduction of a non-steady flow is considered. By using an iterative method for computing the eigenvalues of second kind Fredholm operators (see [1]-[11]) we derive approximation of first eigenvalues closer with respect to results appearing in preceding literature.     

Heat Transfer in Rotor/Rotor and Rotor/Stator Cavity: Physics,Correlations and Numerical Modeling

The transitional and turbulent flow in the near wall sublayer is now mostly modeled based on the existing knowledge of simple 2D flows. To determine the effect of three dimensionality on the turbulent flow structures and turbulent heat transfer in the near wall areas the authors investigate numerically (SVV) turbulent flow in rotor/stator and rotor/rotor flows (with and without axial throughflow). These simple model flows contain most of the phenomena that are needed to understand more complex, 3D transitional and turbulent flows. Attention is focused on the turbulent characteristics which should have more universal character. To stabilize calculations for high Reynolds numbers (up to Re=800 000) the SVV operator is introduced into the Navier-Stokes and energy equations solver for cylindrical coordinate system without using complex numbers. Code optimization and parallelization have speeded up computations 20 times. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Unsteady Slip Flow and Heat Transfer Analysis of Oldroyd-B Fluid Over the Stretching Wedge北大核心CSCD

研究了在速度滑移现象存在下,上随体Oldroyd-B流体绕加热的楔形体的非稳态流动。采用松弛-延迟热通量模型,模拟了传热过程和热延迟时间对传热的影响,通过考虑浮升力、热辐射和对流换热边界条件,进一步研究了流动及传热特性。利用同伦分析方法获得常微分方程组的近似解析解,发现滑移参数的增大可以促进流体的流动,以及流体的温度随热辐射参数增大而升高。此外还发现,温度场在热松弛时间和热延迟时间中出现相反的变化趋势。     

Numerical Computation of Optimal Trajectories for Coplanar, Aeroassisted Orbital Transfer

This paper is concerned with the problem of the optimal coplanaraeroassisted orbital transfer of a spacecraft from a high Earth orbitto a low Earth orbit. It is assumed that the initial and final orbits arecircular and that the gravitational field is central and is governed by theinverse square law. The whole trajectory is assumed to consist of twoimpulsive velocity changes at the begin and end of one interior atmosphericsubarc, where the vehicle is controlled via the lift coefficient.The problem is reduced to the atmospheric part of the trajectory, thusarriving at an optimal control problem with free final time and liftcoefficient as the only (bounded) control variable. For this problem,the necessary conditions of optimal control theory are derived. Applyingmultiple shooting techniques, two trajectories with different controlstructures are computed. The first trajectory is characterized by a liftcoefficient at its minimum value during the whole atmospheric pass. For thesecond trajectory, an optimal control history with a boundary subarcfollowed by a free subarc is chosen. It turns out, that this secondtrajectory satisfies the minimum principle, whereas the first one fails tosatisfy this necessary condition; nevertheless, the characteristicvelocities of the two trajectories differ only in the sixth significantdigit.In the second part of the paper, the assumption of impulsive velocitychanges is dropped. Instead, a more realistic modeling with twofinite-thrust subarcs in the nonatmospheric part of the trajectory isconsidered. The resulting optimal control problem now describes the wholemaneuver including the nonatmospheric parts. It contains as controlvariables the thrust, thrust angle, and lift coefficient. Further,the mass of the vehicle is treated as an additional state variable. For thisoptimal control problem, numerical solutions are presented. They are comparedwith the solutions of the impulsive model.     

Convective Heat Transfer past Small Cylindrical Bodies

Convective heat transfer from an array of small, cylindrical bodies of arbitrary shape in an unbounded, two-dimensional domain is a singular perturbation problem involving an infinite logarithmic expansion in the small parameter ε, representing the order of magnitude of the size of the bodies. Using the method of matched asymptotic expansions, we formulate a hybrid asymptotic-numerical method to solve for the dimensionless, steady-state temperature. We assume that the velocity field of the fluid surrounding the bodies is arbitrary but known. From our asymptotic solution for an arbitrary velocity field, we present the results for two special cases: a uniform flow field and a simple shear flow field. We demonstrate the asymptotic results of the hybrid method through a number of examples and, in a particular case, we compare these results to an exact analytical solution.     

Solvability of Inverse Extremal Problems for Stationary Heat and Mass Transfer Equations

We consider inverse extremal problems for the stationary system of heat and mass transfer equations describing the propagation of a substance in a viscous incompressible heat conducting fluid in a bounded domain with Lipschitz boundary. The problems consist in finding some unknown parameters of a medium or source densities from a certain information of a solution. We study solvability of the direct boundary value problem and the inverse extremal problem, justify application of the Lagrange principle, introduce and analyze the optimality systems, and establish sufficient conditions for uniqueness of solutions.     

Theoretical and Numerical Analysis of an Initial-Boundary Value Problem for the Radiative Transfer Equation with Fresnel Matching Conditions

A Cauchy problem for the time-dependent radiative transfer equation in a three-dimensional multicomponent medium with generalized matching conditions describing Fresnel reflection and refraction at the interface of the media is considered. The unique solvability of the problem is proven, a Monte Carlo method for solving the initial-boundary value problem is developed, and computational experiments for different implementations of the algorithm are conducted.