Numerical solution of the transport equation for passive contaminants in three-dimensional complex terrains

This paper presents a numerical solution of the transport equation for heat and species in complex three-dimensional spaces. The solution domain of the equation is transformed into a cube, as also is the governing equation; the resultant equation is solved in the transformed space via a finite difference technique. The validity of the developed computer code is tested by predicting test cases for which either analytical or reliable experimental results exist. Results are also presented for the rate of convergence of the method and the computer storage requirements, from which the validity, the flexibility and the economy of the developed method are proved for flows in real three-dimensional complex terrains.     

Analytical approximation for solidification processes in PCM storage with internal fins: imposed heat flux

Uses of thermal energy storage systems have expanded notably in recent decades. In thermal energy systems, internal heat transfer enhancement techniques such as fins are often used because of the low thermal conductivity of the phase change materials (PCMs). In this paper, solidification of a PCM is studied in a rectangular storage with horizontal internal plate fins and an imposed constant heat flux on the vertical walls. A simplified analytical solution is presented and its results are compared to those for a numerical approach based on an enthalpy method. The fraction of solidified PCM in storage is calculated with the derived analytical model which determines how much of the storage is solidified after a certain time. The results show that the analytical model satisfactorily estimates the solid–liquid interface and the temperature distribution for the fin, which are useful in the design of PCM-based thermal energy storage or cooling systems.     

Thermal stresses in infinite circular cylinder subjected to rotation

The present investigation is concerned with the effect of rotation on an infinite circular cylinder subjected to certain boundary conditions.An analytical procedure for evaluation of thermal stresses,displacements,and temperature in rotating cylinder subjected to thermal load along the radius is presented.The dynamic thermal stresses in an infinite elastic cylinder of radius a due to a constant temperature applied to a variable portion of the curved surface while the rest of surface is maintained at zero temperature are discussed.Such situation can arise due to melting of insulating material deposited on the surface cylinder.A solution and numerical results are obtained for the stress components,displacement components,and temperature.The results obtained from the present semi-analytical method are in good agreement with those obtained by using the previously developed methods.     

Thermal stresses in an isotropic trimaterial interacted with a pair of point heat source and heat sink

A general analytical solution for an isotropic trimaterial interacted with a point heat source is provided in this paper. Based on the method of analytical continuation in conjunction with the alternating technique, the solutions to heat conduction and thermoelasticity problems for three dissimilar media are first derived. A rapidly convergent series solution for both the temperature and stress functions, which is expressed in terms of an explicit general term of the complex potential of the corresponding homogeneous problem, is obtained in an elegant form. As a numerical illustration, the distributions of thermal stresses along the interface are presented for various material combinations and for different positions of the applied heat source and heat sink.     

Wave motion in a thick cylindrical rod undergoing longitudinal impact

The paper presents a new formulation and the comprehensive analytical solution to longitudinal impact of thick elastic rods. In contrast to previously published works, the solution is derived based on the exact three-dimensional theory without using the classic Skalak’s decomposition. This new direct approach makes the analytical solution more transparent and much easier to obtain. The resulting formulas for basic mechanical quantities are derived using the residue theorem and their evaluation is made in such a way that the accuracy of presented results is significantly higher than those previously published. Based on these results, the transient wave phenomena occurring in the rods are discussed in detail. Additionally, the solution in time domain is obtained also by semi-analytical approach making use of numerical inverse Laplace transform. It is shown that the selected FFT based algorithm is accurate and robust enough, such that the analysis of wave motion in spatial and time domain can be done effectively preserving the results precision. Presented solution can be used as a benchmark for verification of numerical and experimental methods applied to elastodynamics problems.     

Flow and heat transfer over hyperbolic stretching sheets

The boundary layer flow and heat transfer analysis of an incompressible viscous fluid for a hyperbolically stretching sheet is presented. The analytical and numerical results are obtained by a series expansion method and a local non-similarity (LNS) method, respectively. The analytical and numerical results for the skin friction and the Nusselt number are calculated and compared with each other. The significant observation is that the momentum and the thermal boundary layer thickness decrease as the distance from the leading edge increases. The well-known solution of linear stretching is found as the leading order solution for the hyperbolic stretching.     

PRECISE INTEGRAL ALGORITHM BASED SOLUTION FOR TRANSIENT INVERSE HEAT CONDUCTION PROBLEMS WITH MULTI-VARIABLES

IntroductionIHCPs (InverseHeatConductionProblems)arecloselyassociatedwithmanyengineeringaspects,andwelldocumentedintheliteratures,coveringtheidentificationsofthermalparameters[1,2 ],boundaryshapes[3],boundaryconditions[4 ]andsource_relatedterms[5 ,6 ]etc .Howeveritseemsthatonlylittleworkisdirectlyconcernedwithmulti_variablesidentificationsbyauthors’knowledge.Tsengetal.presentedanapproachtodeterminingtwokindsofvariables[7],butonlygavefewnumericalexamplestodeterminethemsimultaneously .Thesol…     

Mechanism of unsteady aerodynamic heating with sudden change in surface temperature

The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.     

Exact analytical solution of the MHD Jeffery-Hamel flow problem

This article shows that the well known nonlinear boundary value problem namely MHD Jeffery-Hamel flow problem, investigated in recent years by many numerical and semi-analytical approximative methods, is exactly solvable and furthermore, gives analytical exact solution in the implicit form for further physical interpretation.     

A semi-analytical solution for linearized multicomponent cation exchange reactive transport in groundwater

Cation exchange in groundwater is one of the dominant surface reactions. Mass transfer of cation exchanging pollutants in groundwater is highly nonlinear due to the complex nonlinearities of exchange isotherms. This makes difficult to derive analytical solutions for transport equations. Available analytical solutions are valid only for binary cation exchange transport in 1-D and often disregard dispersion. Here we present a semi-analytical solution for linearized multication exchange reactive transport in steady 1-, 2- or 3-D groundwater flow. Nonlinear cation exchange mass–action–law equations are first linearized by means of a first-order Taylor expansion of log-concentrations around some selected reference concentrations and then substituted into transport equations. The resulting set of coupled partial differential equations (PDEs) are decoupled by means of a matrix similarity transformation which is applied also to boundary and initial concentrations. Uncoupled PDE’s are solved by standard analytical solutions. Concentrations of the original problem are obtained by back-transforming the solution of uncoupled PDEs. The semi-analytical solution compares well with nonlinear numerical solutions computed with a reactive transport code (CORE^2D) for several 1-D test cases involving two and three cations having moderate retardation factors. Deviations of the semi-analytical solution from numerical solutions increase with increasing cation exchange capacity (CEC), but do not depend on Peclet number. The semi-analytical solution captures the fronts of ternary systems in an approximate manner and tends to oversmooth sharp fronts for large retardation factors. The semi-analytical solution performs better with reference concentrations equal to the arithmetic average of boundary and initial concentrations than it does with reference concentrations derived from the arithmetic average of log-concentrations of boundary and initial waters.     

Fully-developed pipe and planar flows of multimode viscoelastic fluids

Two solutions are presented for fully-developed pipe and planar flows of multimode viscoelastic models. The fluids have a Newtonian solvent contribution and the polymer modes are described by the Phan-Thien—Tanner (PTT), the FENE-P or the Giesekus equation. The first solution is exact and can handle any number of modes, but is only semi-analytical. The second solution, which is presented only for the PTT model with a linear stress coefficient and the FENE-P model, can also handle any number of modes. It is based on a truncated series expansion and is completely analytical, but provides only an approximated solution. The complexity of the multimode solutions is investigated first with the exact semi-analytical method and it is shown that at high Deborah number flows the high-order stresses can become as important as the stress of the first mode. It is also under these conditions that the approximated analytical solution deviates from the exact semi-analytical solution. A criterion for the accurate use of the approximated solution is presented. Fortran codes are provided to obtain these solutions at the internet address at the end.     

A semi-analytical elastic stress–displacement solution for notched circular openings in rocks

A semi-analytical plane elasticity solution of the circular hole with diametrically opposite notches in a homogeneous and isotropic geomaterial is presented. This solution is based on: (i) the evaluation of the conformal mapping function of a hole of prescribed shape by an appropriate numerical scheme and (ii) the closed-form solutions of the Kolosov–Muskhelishvili complex potentials. For the particular case of circular notches––which resemble to the circular cavity breakout in rocks––it is demonstrated that numerical results pertaining to boundary stresses and displacements predicted by the finite differences model FLAC^2D, as well as previous analytical results referring to the stress-concentration-factor, are in agreement with analytical results. It is also illustrated that the solution may be easily applied to non-rounded diametrically opposite notch geometries, such as “dog-eared” breakouts by properly selecting the respective conformal mapping function via the methodology presented herein. By employing a stress-mean-value brittle failure criterion that takes into account the stress-gradient effect in the vicinity of the curved surfaces in rock as well as the present semi-analytical solution, it is found that a notched hole, e.g. borehole or tunnel breakout, may exhibit stable propagation. The practical significance of the proposed solution lies in the fact that it can be used as a quick-solver for back-analysis of borehole breakout images obtained in situ via a televiewer for the estimation of the orientation and magnitude of in situ stresses and of strain–stress measurements in laboratory tests.     

Analytical solution for a three-dimensional non-homogeneous bivariate population balance equation—a special case

There has been a dramatic increase in the number of research publications using the population balance equation (PBE). The PBE allows the prediction of the spatial distribution of the dispersed phase size for an accurate estimation of the flow fields in multiphase flows. A few recent studies have proposed new efficient numerical methods to solve non-homogeneous multivariate PBE and implemented the same in computational fluid dynamics (CFD) codes. However, these codes are generally benchmarked against other numerical methods and applied without verification. To address this gap, an analytical solution for a three-dimensional non-homogeneous bivariate PBE is presented here for the first time. The method of manufactured solutions (MMS) has been used to construct a solution of the non-homogeneous PBE containing breakage and coalescence terms, and an additional source term appearing as a result of this method. The analytical solution presented in this work can be used for the rigorous verification of computer codes written to solve the non-homogeneous bivariate PBE. Quantification of the errors due to different numerical schemes will also become possible with the availability of this analytical solution for the PBE.     

Thermal rhythms in composite ladle walls

Conclusions Analytical and numerical solutions have been obtained for the propagation of thermal rhythms in a composite ladle wall. In particular an analytical model simulating switch on-off thermal conditions at the inner ladle wall with Newton Cooling at the outer ladle wall is found to be in exact agreement with the full numerical solution for the composite wall. This analytical solution is seen to possess novel mathematical features.From a practical viewpoint, to obtain minimum heat loss from the ladle of molten steel it is clear that in early operations the insulating region is most effective if it were located close to the inner wall. As the thermal rhythms approach periodicity it is seen that minimum heat loss would be achieved if the insulating region were placed close to the outer ladle wall. Moreover in the later operational stages it was shown that the minimum heat loss is actually obtained if the insulating layer was omitted altogether. Since it is also shown that the time taken to achieve a periodic thermal rhythm is always least when the insulating region is located centrally in the composite ladle wall it can be argued that the evolution of the thermal rhythm is controlled by its inclusion. Inclusion of the insulating layer yields a marked reduction in the outside wall temperature and this could have considerable bearing on the life of the ladle in an integrated steel works.Preheating the ladle walls was found to reduce the ladle heat loss by more than 50 percent during the initial cycle.Finally it is hoped that the above theoretical results will be of considerable value in the understanding of thermal stratification in ladles of molten steel during standing and pouring.     

Analytical model for melting in a semi-infinite PCM storage with an internal fin

The most PCMs with high energy storage density have an unacceptably low heat conductivity and hence internal heat transfer enhancement techniques such as fins or other metal structures are required in latent heat thermal storage (LHTS) applications. Previous work has concentrated on numerical and experimental examination in determining the influence of the fins in melting phase change material. This paper presents a simplified analytical model based on a quasi-linear, transient, thin-fin equation which predicts the solid-liquid interface location and temperature distribution of the fin in the melting process with a constant imposed end-wall temperature. The analytical results are compared to the numerical results and they show good agreement. Due to the assumptions made in the model, the speed of the solid-liquid interface during the melting process is slightly too slow.     

Numerical solution of thermal entry length problem with variable viscosities and viscous dissipation

A numerical solution of the thermal entry length problem for circular pipes with constant temperature boundaries and with a prescribed dependence of viscosity on temperature and viscous dissipation of some specified order are presented. The results indicate a significant effect of these factors. The results for both heating and cooling processes are presented and compared with the available analytical results. They are useful in the design of heat exchangers with high Prandtl number fluids as the working media.     

Analytical solution to the problem of heat interaction of fluid flow and solid filling

The results of analytical solution to the problem of thermal interaction of a fluid stream with a media of solid particles found on the basis of the ODTPM model [1] are presented. A non-stationary behaviour of the volume-averaged temperature of solid filling, being a measure of the efficiency from heat transfer between fluid and solid phases, and of the heat-carrier temperature at the exit of the layer is studied. Time dependences computed according to the obtained relationships are in good agreement with the results of numerical calculations [1] for different sets of the problem parameters. The found solution expands the possibilities to analyze the efficiency of the “fluid-solid phase” systems and can be a basis to develop methods of experimental evaluation of different parameters of the system.     

等几何壳体分析与形状优化

首先基于Reissner-Mindlin理论进行了三维壳体等几何分析,而后基于此对三维壳体进行形状优化,提出了形状优化中灵敏度的全解析计算方法,包括位移应变阵、雅克比阵和刚度阵等相对控制顶点位置的灵敏度解析计算公式;通过实例验证了壳体等几何分析和灵敏度全解析计算方法的有效性。与传统的基于网格的灵敏度半解析计算方法相比,基于NURBS的灵敏度全解析计算具有精确、计算效率高的特点,且可以避免优化迭代中的网格畸变。     

Thermal analysis of composite superconductors subjected to time-dependent disturbances

A one-dimensional heat conduction equation with time- and temperature-dependent heat sources was employed to study the steady-state and transient response of a composite superconductor subjected to a thermal disturbance. An integral formulation was used to solve the steady-state problem of current redistribution and heat generation. The results of the integral formulation are compared with those of an analytical solution. The two solutions agree with each other except when the analytical solution fails as the temperature in the superconductor begins to exceed the critical temperature. Transient solutions were obtained by the finite-difference technique and the results are compared with a known analytical solution. Results of numerical calculations of the transient response of a composite superconductor subjected to an initial pulsed disturbance are presented. It is demonstrated that the superconductor can switch between the superconducting and the current-sharing state. The transient response and the stability of the composite conductor depend on the magnitude and duration of the disturbance, the dimensionless temperature θ^*, and the dimensionless parameter φ. Received on 18 November 1996     

反平面断裂问题的无单元伽辽金比例边界法

将比例边界法与无单元伽辽金法相结合,建立了反平面断裂分析的无单元伽辽金比例边界法。这是一种边界型无网格法,在环向方向上采用无单元伽辽金法进行离散,因此计算时仅需要边界上的节点信息,不需要边界元所要求的基本解。为了便于施加本质边界条件,通过建立节点值和虚拟节点值之间的关系给出了修正的移动最小二乘形函数。在径向方向上,该方法利用解析的方法求解,因此是一种半解析的数值方法。最后,给出了数值算例,并验证了所提方法后处理简单和计算精度高的特点,适合于求解反平面断裂问题。