A nonlinear inverse problem for the prediction of local thermal contact conductance in plate finned-tube heat exchangers

A nonlinear inverse problem utilizing the Conjugate Gradient Method (CGM) of minimization is used successfully to estimate the temporally and circumferentially varying thermal contact conductance of a plate finned-tube heat exchanger by reading the simulated transient temperature measurement data from the thermocouples located on the plate. The thermal properties of the fin and tube are assumed to be functions of temperature, and this makes the problem nonlinear. It is assumed that no prior information is available on the functional form of the unknown thermal contact conductance in the present study, thus, it is classified as the function estimation in the inverse calculation. The accuracy of the inverse analysis is examined by using the simulated temperature measurements. Finally the inverse solutions with and without the consideration of temperature-dependent thermal properties are compared. Results show that when the nonlinear inverse calculations are performed an excellent estimation on the thermal contact conductance can be obtained with any arbitrary initial guesses within a couple of minute's CPU time on a HP-730 workstation.     

Numerical nonlinear inverse problem of determining wall heat flux

The inverse problem of determining time-variable surface heat flux in a plane wall, with constant or temperature dependent thermal properties, is numerically studied. Different kinds of incident heat flux, including rectangular waveform, are assumed. The solution is numerically solved as a function estimation problem, so that no a priori information for the functional waveforms of the unknown heat flux is needed. In all cases, a solution in the form of a piece-wise function is used to approach the incident flux. Transient temperature measurements at the boundary, from the solution of the direct problem, served as the simulated experimental data needed as input for the inverse analysis. Both direct and inverse heat conduction problems are solved using the network simulation method. The solution is obtained step-by-step by minimising the classical functional that compares the above input data with those obtained from the solution of the inverse problem. A straight line of variable slope and length is used for each one of the stretches of the desired solution. The influence of random error, number of functional terms and the effect of sensor location are studied. In all cases, the results closely agree with the solution.     

Inverse thermal analysis of the drying zone of the evaporator of an axially grooved heat pipe

An inverse approach is performed to characterize the thermal behaviour of an axially grooved heat pipe, in steady state, for various operating conditions. For this purpose, an experimental set up, as well as a network conduction model, are developed to simulate the heat transfer in the wall at the evaporator section. The minimization of an objective function, taking into account the discrepancy between measured temperatures and computed ones, allows then the estimation of a heat transfer coefficient as well as the drying out front positions for all the axial grooves. Hence, at the burnout point, the significant temperature increase in the evaporator extremity is considered to be a direct consequence of the restriction of the evaporative zone. Therefore, the distribution of liquid phase in the capillary structure of the heat pipe can be obtained through the analysis of the measured temperature gradient in the evaporator section where the dry out front was expected to occur. Furthermore, the dry out front expansion can be observed when the input heat load is increased or when the adiabatic temperature is decreased. Introducing an adverse tilt angle also shows the effect of the puddle.     

Retrieving the heat transfer coefficient for jet impingement from transient temperature measurements

Algorithm of retrieving the heat transfer coefficient (HTC) from transient temperature measurements is presented. The unknown distributions of two types of boundary conditions: the temperature and heat flux are parameterized using a small number of user defined functions. The solutions of the direct heat conduction problems with known boundary temperature and flux are expressed as a superposition of auxiliary temperature fields multiplied by unknown parameters. Inverse problem is formulated as a least squares fit of calculated and measured temperatures and is cast in a form of a sum of two objective functions. The first results originates from an inverse problem for retrieving the boundary temperature the second comes from the inverse problem for reproducing the boundary heat flux. The final form of the objective function is obtained by enforcing constant in time value of the heat transfer coefficient. This approach leads to substantial regularization of the results, when compared with the standard technique, where HTC is calculated from separately reconstructed temperature and heat flux on the boundary. The validation of the numerical procedure is carried out by reconstructing a known distribution of the HTC using simulated measurements laden by stochastic error. The proposed approach is also used to reconstruct the distribution of the HTC in a physical experiment of heating a cylindrical sample using an impinging jet.     

基于数据驱动模型求解热传导反问题

热传导反问题求解在工程领域具有重要的应用价值.本文发展数据驱动模型识别了管道内壁几何形状和皮肤肿瘤生长参数等热传导反问题.在管道内壁几何形状识别问题中,首先采用随机生成模型结合有限元法求解热传导正问题,并采用有效导热系数转化的思想,建立机器学习模型,求解了测点温度与有效导热系数之间的抽象映射关系,进而实现管道内壁几何形状的识别.然后,应用数据驱动模型识别了皮肤肿瘤的生长参数,分别讨论了不同测量误差对计算结果的影响.数值算例表明,本文提出的数据驱动模型能够准确估算肿瘤的生热率和血液灌注率.这些工作显示了数据驱动模型在求解热传导反问题方面具有广阔的应用前景.     

Dynamic analysis of functionally graded plates using the hybrid Fourier-Laplace transform under thermomechanical loading

In this study, the analytical solution is presented for dynamic response of a simply supported functionally graded rectangular plate subjected to a lateral thermomechanical loading. The first-order and third-order shear deformation theories and the hybrid Fourier-Laplace transform method are used. The material properties of the plate, except Poisson’s ratio, are assumed to be graded in the thickness direction according to a power-law distribution in terms of the volume fractions of the constituents. The plate is subjected to a heat flux on the bottom surface and convection on the upper surface. A third-order polynomial temperature profile is considered across the plate thickness with unknown constants. The constants are obtained by substituting the profile into the energy equation and applying the Galerkin method. The obtained temperature profile is considered along with the equations of motion. The governing partial differential equations are solved using the finite Fourier transformation method. Using the Laplace transform, the unknown variables are obtained in the Laplace domain. Applying the analytical Laplace inverse method, the solution in the time domain is derived. The computed results for static, free vibration, and dynamic problems are presented for different power law indices for a plate with simply supported boundary conditions. The results are validated with the known data reported in the literature. Furthermore, the results calculated by the analytical Laplace inversion method are compared with those obtained by the numerical Newmark method.     

Efficient reconstruction of local heat fluxes in pool boiling experiments by goal-oriented adaptive mesh refinement

In this paper, we consider the efficient estimation of local boiling heat fluxes from transient temperature measurements in the heater close to the heater surface. For accurate prediction, heat flux estimation is formulated as a transient three-dimensional (3D) inverse heat conduction problem (IHCP). This inverse problem is ill-posed and cannot be treated straightforwardly by established numerical methods. In order to obtain a regularized stable solution, a large-scale time-dependent PDE-constrained optimization problem has to be solved and an appropriate stopping criterion for the termination of the iterative solution process has to be chosen. Since the boiling heat flux is non-uniformly distributed on the heater surface due to the strong local activity of the boiling process, the use of a fixed uniform spatial discretization is not efficient. Instead, an adaptive mesh refinement strategy can be used to obtain an appropriate discretization which significantly reduces the total computational effort. In this work, we present an automatic algorithm incorporating an adaptive mesh refinement via a heat flux-based a-posteriori error estimation technique. The suggested algorithm can cope with both spatially point-wise or highly resolved temperature observations efficiently. It is applied to real measurement data obtained from two different types of pool boiling experiments. The numerical results show that the computational effort can be reduced significantly for given estimation quality. This adaptive IHCP solution technique can be also viewed as an efficient soft sensor to deduce unmeasurable local boiling heat fluxes.     

Effect of measurement volume size on turbulent flow measurement using ultrasonic Doppler method

This paper presents the results of an investigation on the effects of measurement volume size on the mean velocity profile and the Reynolds stress for fully developed turbulent pipe flows. The study employs the ultrasonic velocity profile method, which is based on the ultrasonic Doppler method. The ultrasonic Doppler method offers many advantages over conventional methods for flow rate measurement in the nuclear power plant piping system. This method is capable of measuring the instantaneous velocity profile along the measuring line and is applicable for opaque liquids and opaque pipe wall materials. Furthermore, the method has the characteristic of being non-intrusive. Although it is applicable to various flow conditions, it requires a relatively large measurement volume. The measurement volume of the present method has a disk-shape determined by the effective diameter of the piezoelectric element and the number of the wave cycles of the ultrasonic pulse. Considering this disk-shaped measurement volume and expressing the time-averaged velocity in a truncated Taylor series expansion around the value at the center of the measuring control volume, the value of the velocity can be obtained. The results are then compared with the data obtained from DNS and LDA measurements. The result shows that the effect of the measurement volume size appears in the buffer region and viscous sublayer.     

A method for estimating wall friction in turbulent wall-bounded flows

We describe a simple method for estimating turbulent boundary layer wall friction using the fit of measured velocity data to a boundary layer model profile that extends the logarithmic profile all the way to the wall. Two models for the boundary layer profile are examined, the power-series interpolation scheme of Spalding and the Musker profile which is based on the eddy viscosity concept. The performance of the method is quantified using recent experimental data in zero pressure gradient flat-plate turbulent boundary layers, and favorable pressure gradient turbulent boundary layers in a pipe, for which independent measurements of wall shear are also available. Between the two model profiles tested, the Musker profile performs much better than the Spalding profile. Results show that the new procedure can provide highly accurate estimates of wall shear with a mean error of about 0.5% in friction velocity, or 1% in shear stress, an accuracy that is comparable to that from independent direct measurements of wall shear stress. An important advantage of the method is its ability to provide accurate estimates of wall shear not only based on many data points in a velocity profile but also very sparse data points in the velocity profile, including only a single data point such as that originating from a near-wall probe.     

Rheometry using velocity measurements

Traditional methods of rheometry employ simple flows such as viscometric flows and measure stress or volumetric flow rate to determine the rheological parameters in the constitutive equation. One can find analytic solutions for stress and volumetric flow rates for these simple flows, and comparison of them with experimental data determines rheological parameter values. In the present investigation, rheological parameters are estimated by measuring velocity at certain locations. A pulsatile flow in a circular pipe, which can be implemented easily, is adopted to estimate rheological parameters in a general constitutive equation. The inverse problem of determining the rheological parameters from velocity measurements is solved using a conjugate gradient method. The present method is found to yield a reasonably accurate estimation of rheological parameters even with noisy velocity measurements.     

Analysis of laminar flow forced convection heat transfer in the entrance region of a circular pipe

A numerical solution, for incompressible, steady-state, laminar flow heat transfer in the combined entrance region of a circular tube is presented for the case of constant wall heat flux and constant wall temperature. The development of velocity profile is obtained from Sparrow's entrance region solution. This velocity distribution is used in solving the energy equation numerically to obtain temperature profiles. Variation of the heat transfer coefficient for these two different boundary conditions for the early stages of boundary layer formation on the pipe wall is obtained. Local Nusselt numbers are calculated and the results are compared with those given byUlrichson andSchmitz. The effect of the thermal boundary conditions is studied by comparing the uniform wall heat flux results with uniform wall temperature.     

Inverse approach and sensitivity analysis for identification of ingot-mould thermal resistance in continuous casting of metals

The problem of appropriate location of the sensors for identification of ingot – mould thermal resistance during continuous casting of metals is the subject of the paper. Analysed problem belongs to the group of inverse problems. The present work shows also the method of identification of unknown thermal resistance using the temperature measurements at the number of sensors located in the wall of the mould. The influence of the location of the sensors on the results of identification is analysed. The best location of the sensors results from the sensitivity analysis for the steady-state inverse heat conduction problem. Validation of the proposed inverse method is realized by comparison of the results taken from solution of inverse and direct problems. Several numerical examples are presented and analysed.     

Conjugate heat transfer in thermally developing laminar flow with viscous dissipation effects

In this study, thermally developing laminar forced convection in a pipe including viscous dissipation and wall conductance is investigated numerically. The constant heat flux is assumed to be imposed at the outer surface of the pipe wall. The finite volume method is used. The distributions for the developing temperature and local Nusselt number in the entrance region are obtained. The dependence of the results on the Brinkman number and the dimensionless thermal conductivity are shown. The viscous heating effect on the wall is shown. Significant viscous dissipation effects have been observed for large Br.     

Heat flux sensor with minimal impact on boundary conditions

A technique for determining the heat transfer on the far surface of a wall based on measuring the heat flux and temperature on the near wall is presented. Although heat transfer measurements have previously been used to augment temperature measurements in inverse heat conduction methods, the sensors used alter the heat flow through the surface, disturbing the very quantity that is desired to be measured. The ideal sensor would not alter the boundary condition that would exist were the sensor not present. The innovation of this technique in that it has minimal impact on the wall boundary condition. Since the sensor is placed on the surface of the wall, no alteration of the wall is needed. The theoretical basis for the experimental technique as well as experimental results showing the heat flux sensor performance is presented.     

Application of an Inverse Simulator of Single-Phase Flow Analysis for Leakage Pathway Estimation in a Multiphase System for Geologic $$\hbox {CO}_{2}$$ Storage

When \(\hbox {CO}_{2}\) is injected in a brine reservoir, brine or \(\hbox {CO}_{2}\) can be discharged into a permeable formation saturated with brine above the storage reservoir along a leakage pathway, if present. In most situations, the overlying formation can act as a single-phase aquifer with only brine leakage before \(\hbox {CO}_{2}\) leaks. This study examines the applicability of a developed inverse code for single-phase fluids to detect leakage pathway locations in view of the overlying formation using pressure anomalies induced by leaks. Before applying inverse analysis, forward modeling is performed using the TOUGH2 model to determine brine and \(\hbox {CO}_{2}\) leakage in a homogeneous conceptual model, and the simulated pressure profiles at monitoring wells are used as measurements in the inverse model. In the inverse code, an important consideration is that the vertical hydraulic conductivity and cross-sectional area of a leakage pathway that are inherent to a leakage term in the mass balance equation are integrated as a single parameter to estimate the leakage pathway locations. This method eliminates the impact of the uncertainty of the leakage pathway size on the accuracy of leakage pathway estimation. The inverse model examines the effect of the number of monitoring wells, monitoring periods and \(\hbox {CO}_{2}\) leakage into the overlying formation on the accuracy of leakage pathway estimation according to eleven application examples. The comparison between the results of the single-phase inverse code and iTOUGH2 code illustrates that the single-phase inverse model can be applicable to the leakage pathway estimation in a multiphase flow system.     

Drag reduction and heat transfer in surfactant solutions with excess counterion

We present the results of a study of turbulent drag reduction in a small circulating loop using surfactant solutions with excess counterion. In addition, these solutions were used in measurements of heat transfer, both in pipe flow and in an impinging jet. Both frictional drag and heat transfer were reduced in the pipe flow experiments. Measurements of heat transfer in the impinging jet revealed a dependence on the molar concentration ratio of the counterion. When the counterion was added at a molar concentration 30 times higher than that of the surfactant, the resulting surfactant solution did not reduce the rate of heat transfer in the impinging jet. By using this surfactant system in an impinging jet, we show both a reduction in pipe friction and normal heat transfer potential in a circulating heat exchange system. In order to investigate this difference in heat transfer between pipe flows and impinging jet flows, a comparison was made of the wall shear stress between these two flow regimes. The estimated wall shear stress was of the same order in both flows, and thus was not considered to be the primary cause of the difference in heat transfer. It is instead suggested that the micellar structure of the surfactant is influenced by a compressive deformation of the impinging flow in a manner that is different from the shear deformation observed in pipe flow.     

Slip ratio in dispersed viscous oil-water pipe flow

In this article, dispersed flow of viscous oil and water is investigated. The experimental work was performed in a 26.2-mm-i.d. 12-m-long horizontal glass pipe using water and oil (viscosity of 100 mPa s and density of 860 kg/m³) as test fluids. High-speed video recording and a new wire-mesh sensor based on capacitance (permittivity) measurements were used to characterize the flow. Furthermore, holdup data were obtained using quick-closing-valves technique (QCV). An interesting finding was the oil-water slip ratio greater than one for dispersed flow at high Reynolds number. Chordal phase fraction distribution diagrams and images of the holdup distribution over the pipe cross-section obtained via wire-mesh sensor indicated a significant amount of water near to the pipe wall for the three different dispersed flow patterns identified in this study: oil-in-water homogeneous dispersion (o/w H), oil-in-water non-homogeneous dispersion (o/w NH) and Dual continuous (Do/w & Dw/o). The phase slip might be explained by the existence of a water film surrounding the homogeneous mixture of oil-in-water in a hidrofilic-oilfobic pipe.     

Local heat transfer from a hot plate to a water jet

Jet impingement boiling is very efficient in cooling of hot surfaces as a part of the impinging liquid evaporates. Because of its importance to many cooling procedures, investigations on basic mechanisms of jet impingement boiling heat transfer are needed. Until now, most of the experimental studies, carried out under steady-state conditions, used a heat flux controlled system and were limited by the critical heat flux (CHF). The present study focuses on steady-state experiments along the entire boiling curve for hot plate temperatures of up to 700°C. A test section has been built up simulating a hot plate. It is divided into 8 independently heated modules of 10 mm length to enable local heat transfer measurements. By means of temperature controlled systems for each module local steady-state experiments in the whole range between single phase heat transfer and film boiling are possible. By solving the two dimensional inverse heat conduction problem, the local heat flux and the corresponding wall temperature on the surface of each module can be computed. The measurements show important differences between boiling curves measured at the stagnation line and those obtained in the parallel flow region. At the stagnation line, the transition boiling regime is characterised by very high heat fluxes, extended to large wall superheats. Inversely, boiling curves in the parallel flow region are very near to classical ones obtained for forced convection boiling. The analysis of temperature fluctuations measured at a depth of 0.8 mm from the boiling surface enables some conclusions on the boiling mechanism in the different boiling regimes.     

Reduction of beam refraction in optical pipe flow experiments by use of sheet-fabricated pipe walls

A measurement section is presented that allows for the study of liquid flows in straight pipes with noninvasive optical techniques. The amount of refraction of light rays is minimized by the use of short pipe sections manufactured of transparent film. For the situation of water, it is shown that walls as thin as 85?μm can be used in turbulent flows at high velocities. In our specific case, a flow of water at Reynolds numbers upto Re _D =300?000 in a pipe of Ø 70?mm diameter (average velocity 4?m/s) was realized without unacceptable vibrations of the pipe wall. In this situation laser Doppler measurements can be carried out without the need to correct for the position of the measurement volume or changes in the velocity response of the system.     

Cross-flow heat transfer in fixed bed

Radial flow reactor operated at cross-flow heat transfer is focused for large scale methanol synthesis. The effects of operating conditions including the reactor inlet air temperature, the heating pipe temperature and the air flow rate on the cross-flow heat transfer were investigated and results show that the temperature profile of the area in front of the heating pipe is slightly affected by all the operating conditions. The main area whose temperature profile is influenced is located behind the heating pipe. The heat transfer direction is related to the direction of the flow. In order to obtain the basic parameters for radial flow reactor designing calculation, the dimensionless number group method was used for data fitting of the bed effective thermal conductivity and the wall heat transfer coefficient which were calculated by the mathematical model with the product of Reynolds number and Prandtl number. The comparison of experimental data and calculated values shows that the calculated values fit the experimental data satisfactorily and the formulas can be used for reactor designing calculation.