Optimal shape design for packaging containing heating elements by inverse heat transfer method

The present study is concerned with the determination of the optimal shape for a package containing multiple heating elements. The optimization tool has been developed based on the inverse heat transfer (IHT) approach, incorporating a direct problem solver, a numerical grid generator, a direct-differentiation sensitivity analyzer, and the conjugate gradient method. Shape design that leads to a specified outer surface temperature distribution is predicted by the approach. In this study, the effects of internal heat generation on optimal shapes of the packagings have also been evaluated. Several practical cases with various imbedded heating elements and thermal conditions are studied. Results show that the approach provides an efficient computer-aided design scheme for the shape profile determination.     

Unsteady conjugate problem of a dissipative fluid in a horizontal channel with a periodic variation temperature

The unsteady two-dimensional transient heat transfer problem referring to a fully laminar flow developing in a parallel-plane channel exposed to a periodic variation surface temperature with distance is numerically studied. The effects of channel thickness, Péclet number, wall-to-fluid conductivity ratio, thermal diffusivity ratio, angular frequency and the viscous dissipation parameter are determined in the solutions. The non-linear equations are discretized by means an implicit finite difference scheme and the electric analogy to the resulting system is applied to convert these equations into a network-electrical model that was solved using a computer code (electric circuits simulator). In this scheme, only spatial discretization is necessary, while time remains as a real continuous variable, and its programming does not require manipulation of the sophisticated mathematical software that is inherent in other numerical methods. The network simulation method, which satisfies the conservation law for the heat flux variable and the uniqueness law for temperature, also permits the direct visualization of the local and/or integrated transport variables at any point or section of the medium.     

A least-squares procedure for the solution of transport problems

In this paper a least-squares formulation associated with a conjugate gradient algorithm is proposed for the solution of transport problems. In this procedure the advection–diffusion equation is first discretized in time using an implicit scheme. At each time step the resulting partial differential equation is replaced by an optimal control problem. This minimization problem involves the minimization of a functional defined via a state equation. This functional is chosen in order to force the numerical solution of the advection–diffusion equation to be equal to the hyperbolic advective part of this equation. The effectiveness of the method is shown through a one-dimensional example involving advective and diffusive transport. No oscillation and high accuracy have been obtained for the entire range of Peclet numbers with a Courant number well in excess of unity.     

Fully coupled finite volume solutions of the incompressible Navier–Stokes and energy equations using an inexact Newton method

An inexact Newton method is used to solve the steady, incompressible Navier–Stokes and energy equation. Finite volume differencing is employed on a staggered grid using the power law scheme of Patankar. Natural convection in an enclosed cavity is studied as the model problem. Two conjugate-gradient -like algorithms based upon the Lanczos biorthogonalization procedure are used to solve the linear systems arising on each Newton iteration. The first conjugate-gradient-like algorithm is the transpose-free quasi-minimal residual algorithm (TFQMR) and the second is the conjugate gradients squared algorithm (CGS). Incomplete lower-upper (ILU) factorization of the Jacobian matrix is used as a right preconditioner. The performance of the Newton- TFQMR algorithm is studied with regard to different choices for the TFQMR convergence criteria and the amount of fill-in allowed in the ILU factorization. Performance data are compared with results using the Newton-CGS algorithm and previous results using LINPACK banded Gaussian elimination (direct-Newton). The inexact Newton algorithms were found to be CPU competetive with the direct-Newton algorithm for the model problem considered. Among the inexact Newton algorithms, Newton-CGS outperformed Newton- TFQMR with regard to CPU time but was less robust because of the sometimes erratic CGS convergence behaviour.     

Streamline upwind finite element method for conjugate heat transfer problems

This paper presents a combined finite element method for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow. The streamline upwind finite element method is used for the analysis of thermal viscous flow in the fluid region, whereas the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the proposed method is to consistently couple heat transfer along the fluid-solid interface. Three test cases, i.e. conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the efficiency of the present method. The English text was polished byYunming Chen.     

Simulation of two- and three-dimensional internal subsonic flows using a finite element method

In this paper a finite element method is presented to predict internal subsonic flows. Using a low-Mach-number approximation, the pressure is decomposed into a mean thermodynamic contribution and a dynamic fluctuation to deal with the complex role of the pressure in internal aerodynamics. A semi-implicit time integration and a finite element method with a moving mesh are described to take into account complex geometries and moving boundaries. An Uzawa algorithm accelerated by a preconditioned residual method is introduced to solve the coupled non-symmetric linear system for the velocity components and the pressure. An efficient conjugate gradient method combined with an incomplete LU preconditioning is used to solve the non-symmetric linear systems arising from the discretization. The implementation of the numerical scheme on parallel supercomputers is also discussed. Efficient algorithms for the finite element assembly phase and for the solution of linear systems are described which take advantage of the parallel architecture of the new generation of supercomputers. With this technique a global speed-up of 10 is achieved on a supercomputer with eight processors. To illustrate the capabilities of the numerical method, 2D and 3D simulations of flows in the combustion chamber of a reciprocating engine and around the combustor dome of a gas turbine engine are presented.     

Non-symmetric CG-like schemes and the finite element solution of the advection–dispersion equation

Seven leading iterative methods for non-symmetric linear systems (GMRES, BCG, QMR, CGS, Bi-CGSTAB, TFQMR and CGNR) are compared in the specific context of solving the advection–dispersion equation by a classic approach: The space derivatives are approximated by linear finite elements while an implicit scheme is used to integrate the time derivatives. Convergence formulas that predict the behaviour of the iterative methods as a function of the discretization parameters are developed and validated by experiments. It is shown that all methods converge nicely when the coefficent matrix of the linear system is close to normal and the finite element approximation of the advection–dispersion equation yields accurate results.     

Iterative solution techniques for the stokes and Navier-Stokes equations

The linear system arising from a Lagrange-Galerkin mixed finite element approximation of the Navier–Stokes and continuity equations is symmetric indefinite and has the same block structure as a system arising from a mixed finite element discretization of a Stokes problem. This paper considers the iterative solution of such a system, comparing the performance of the one-level preconditioned conjugate residual method for indefinite matrices with that of a more traditional two-level pressure correction approach. Asymptotic estimates for the amount of work involved in each method are given together with the results of related numerical experiments.     

Fast finite difference solution for steady-state Navier-Stokes equations using the BID method

A first biharmonic boundary value problem is obtained by combining the coupled steady-state Navier-Stokes equations in their stream-function-vorticity formulation. This biharmonic boundary value problem is solved by a fast biharmonic solver developed by the authors wherein the idea of preconditioned conjugate gradient method is used. The biharmonic driver (BID) method using this solver has been found fast converging, and produces accurate results up to moderately large Reynolds numbers. Also, the mesh size does not affect the convergence rate.     

Transient finite element method for calculating steady state three-dimensional free surfaces

In order to reduce the cost of large three-dimensional calculations of steady state free surfaces, we have combined a time-dependent approach, a decoupling algorithm and a conjugate gradient solver along the lines introduced earlier by Gresho and Chan. The free surface is calculated separately by applying the kinematic condition to a number of faces defined on the undeformed surface. For the pseudo-time-marching technique we show that it is economical to adopt different time steps for the free surface calculation and the other fields. The accuracy of the method is tested on the well-known circular die problem; the method is then used to reveal the effects of inertia and shear thinning on square and rectangular dies.     

Selective lumping finite element method for shallow water flow

A finite element method for solving shallow water flow problems is presented. The standard Galerkin method is employed for spatial discretization. The numerical integration scheme for the time variation is the explicit two step scheme, which was originated by the authors and their co-workers. However, the original scheme has been improved to remove the erroneous artifical damping effect. Since the improved scheme employs a combination of lumped and unlumped coefficients, the scheme is referred to as a selective lumping scheme. Stability conditions and accuracy are investigated by considering several numerical examples. The method has been applied to the tidal flow in Osaka Bay and Yatsushiro Bay.     

Inertia dominated flow and heat transfer in liquid drop spreading on a hot substrate

The present work deals with computational modeling of the fluid flow and heat transfer taking place in the process of impact of a cold liquid drop (T_d = 20-25 °C) onto a dry heated substrate characterized by different thermophysical properties. The computational model, based on the volume-of-fluid method for the free-surface capturing, is validated by simulating the configurations accounting for the conjugate heat transfer. The simulations were performed in a range of impact Reynolds numbers (Re = 2000-4500), Weber numbers (We = 27-110) and substrate temperatures (T_s = 100-120 °C). The considered temperature range of the drop-surface, i.e. liquid-solid system does not account for the phase change, that is boiling and evaporation. The model performances are assessed by contrasting the results to the reference database originating from the experimental and complementary numerical investigations by Pasandideh-Fard et al. [Pasandideh-Fard, M., Aziz, S., Chandra, S., Mostaghimi, J., 2001. Cooling effectiveness of a water drop impinging on a hot surface. International Journal of Heat and Fluid Flow, 22, 201-210] and Healy et al. [Healy, W., Hartley, J., Abdel-Khalik, S., 2001. On the validity of the adiabatic spreading assumption in droplet impact cooling. International Journal of Heat and Mass Transfer, 44, 3869-3881]. In addition, the thermal field obtained is analyzed along with the corresponding asymptotic analytical solution proposed by Roisman [Roisman, I.V., 2010. Fast forced liquid film spreading on a substrate: flow, heat transfer and phase transition. Journal of Fluid Mechanics, 656, 189-204]. Contrary to some previous numerical studies, the present computational model accounts for the air flow surrounding the liquid drop. This model feature enables a small air bubble to be resolved in the region of the impact point. The reported results agree reasonably well with experimental and theoretical findings with respect to the drop spreading pattern and associated heat flux and temperature distribution.     

A hybrid scheme for the numerical simulation of shock/discontinuity problems

To address accuracy issues for direct numerical simulation, a hybrid scheme based on the weighted compact scheme (WCS) and weighted essentially non-oscillatory (WENO) scheme is developed. The new hybrid method incorporates the advantages of both schemes. Time integration is performed using the fourth-order total variation diminishing Runge–Kutta method with a characteristic filter. The accuracy of the scheme is assessed using several benchmark problems. Results show that the proposed scheme produces a more accurate solution for problems involving shocks and discontinuities in comparison with the traditional shock-capturing methods.     

随机结构静力反应概率密度演化方程的差分方法

随机结构分析的概率密度演化方法是分析随机结构静力反应的一种具有良好前景的方法。本文研究了求解随机结构静力反应概率密度演化方程的差分方法，分别探讨了单边差分格式和Lax-Wendroff格式的计算性态。二者均能满足概率相容性条件并且能够保证均值线性增长。以八层框架结构的静力随机反应为例，对两种差分格式的结果及精确解答进行了具体的比较分析。研究表明，两种差分格式均是收敛和稳定的，在不连续点处存在角点效应．单边差分格式能够保证概率非负性，而Lax-Wendroff格式具有往往更快的收敛速度。就变异系数而言，通常单边差分格式的变异系数随着区间离散数的增长而趋于稳定值，Lax-Wendroff格式则一开始就可得到恒定的值。     

单点子域积分与多点子域积分

基于单点子域精细积分的思想,针对抛物线型热传导方程初边值问题,提出了多点子域积分的概念,推出了一种多点子域积分的FTCS格式。该格式为显格式,并证明其为无条件稳定。数值算例表明,多点子域积分的FTCS格式具有比单点子域积分的FTCS格式收敛速度快的特点。     

AUFS 格式在无网格方法中的应用简

将计算量小,激波分辨率高的AUFS (artificially upstream flux vector splitting) 格式应用于无网格方法. 所发展算法基于多项式基函数最小二乘无网格方法,采用线性基函数曲面拟合及AUFS 格式计算各离散点的空间导数,应用四阶Runge-Kutta 法进行时间显式推进. 为验证算法健壮性、精度以及计算效率,对Riemann 问题、超音速平面流动,以及不同攻角NACA0012 翼型跨音速流场进行了数值模拟,其结果同采用HLLC (Harten-Lax-van Leer-contact) 格式的无网格方法以及文献报道结果吻合较好,并且计算量较形式简单HLLC 格式减少约15%.     

压电材料反平面应变状态的椭圆夹杂及界面裂纹问题

本文采用共法求解了压电材料反平面变形的椭圆夹杂及界面裂纹问题，前者的解答表明当远场外力均匀分布对夹杂内的应力场及电位移场是常量，后者解答表明在界面裂纹的裂尖处，应力及电位移都具有γ＾－１／２的奇异性。     

Parallel algorithms for panel methods

A parallel algorithm for the solution of potential flow problems using the panel method of Hess and Smith and conjugate and bi-conjugate gradient techniques is presented. Analysis of the parallelism for the matrix. solvers shows the algorithms to have scalable properties as the problem size grows indefinitely large. Speed-up and efficiency values are presented along with experimental and theoretical values for the optimum number of processors for maximum speed-up. It is envisaged that the parallel techniques presented here have applications using other boundary integral methods for solving engineering problems of a more complex nature.     

An accurate computation for rapidly varied flow in an open channel

A new method combining the Preissmann four-point scheme and the Holly–Preissmann reach-back scheme is introduced to solve the rapidly varied flow problem in an open channel. The Preissmann four-point scheme is well known for the computation of one-dimensional unsteady flow. The Holly–Preissmann reach-back scheme integrates the Holly-Preissmann two-point scheme with the concept of reach-back characteristics, which allows the characteristics to project several time steps beyond the current time level. A spontaneous surge formation case is used to demonstrate and evaluate the applicability of the new method. It has been found that the results from this method are quite compatible with those of Preissmann four-point scheme. In addition, with the appropriate choice of the number of reach-back time steps, this new method can always avoid the numerical oscillation which usually exists when one uses the Preissmann four-point scheme for the condition of Courant number not close to unity.