Estimation of thermal breakdown voltage of HVDC cables - A theoretical framework

The insulation in a dc cable is subjected to both thermal and electric stress at the same time. While the electric stress is generic to the cable, the temperature rise in the insulation is, by and large, due to the Ohmic losses in the conductor. The consequence of this synergic effect is to reduce the maximum operating voltage and causes a premature failure of the cable. The authors examine this subject in some detail and propose a comprehensive theoretical formulation relating the maximum thermal voltage (MTV) to the physical and geometrical parameters of the insulation. The heat flow patterns and boundary conditions considered by the authors here and those found in earlier literature are provided. The MTV of a dc cable is shown to be a function of the load current apart from the resistance of the insulation. The results obtained using the expressions, developed by the authors, are compared with relevant results published in the literature and found to be in close conformity.

Thermal stresses in rectangular plates

A method of determining the thermal stresses in a flat rectangular isotropic plate of constant thickness with arbitrary temperature distribution in the plane of the plate and with no variation in temperature through the thickness is presented. The thermal stress have been obtained in terms of Fourier series and integrals that satisfy the differential equation and the boundary conditions. Several examples have been presented to show the application of the method.

Thermal stresses in a spherical shell with a conical nozzle

A thermal stress problem of a spherical shell with a conical nozzle is solved using a continuum approach. The thermal loading consists of a steady temperature which is uniform on the inner and outer surfaces of the shell and the conical nozzle but may vary linearly across the thickness. The thermal stress problem is converted to an equivalent boundary value problem and boundary conditions are specified at the junction of the spherical shell and conical nozzle. The stresses are obtained for a uniform increase in temperature and for a linear variation of temperature across the thickness of the shell, and are presented in graphical form for ready use.

Analysis of continuum using boundary compatibility conditions of integrated force method

Static and vibration problems of an indeterminate continuum are traditionally analyzed by the stiffness method. The force method is more or less non-existent for such problems. This situation is primarily due to the incomplete state of development of the compatibility conditions which are essential for the analysis of indeterminate structures by the flexibility method. The understanding of the Compatibility Conditions (CC) has been substantially augmented. Based on the understanding of CC, a novel formulation termed the Integrated Force Method (IFM) has been established. In this paper IFM has been extended for the static and vibration analyses of a continuum. The IFM analysis is illustrated taking three examples: 1. (1) rectangular plate in flexure 2. (2) analysis of a cantilevered dam 3. (3) free vibration analysis of a beam. From the examples solved it is observed that the force response of an indeterminate continuum with mixed boundary conditions can be generated by IFM without any reference to displacements in the field or on the boundary. Displacements if required can be calculated by back substitution.

Intensity Distribution Design for Laser-Induced Thermal Loading Based on Numerical Simulation

To accomplish laser-induced thermal loading simulation tests for pistons,the Gaussian beam was modulated into multi-circular beam with specific intensity distribution.A reverse method was proposed to design the intensity distribution for the laser-induced thermal loading based on finite element(FE) analysis.Firstly,the FE model is improved by alternating parameters of boundary conditions and thermal-physical properties of piston material in a reasonable range,therefore it can simulate the experimental resul...

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Thermal fatigue on pistons induced by shaped high power laser. Part II: Design of spatial intensity distribution via numerical simulation

In the laser induced thermal fatigue simulation test on pistons, the high power laser was transformed from the incident Gaussian beam into a concentric multi-circular pattern with specific intensity ratio. The spatial intensity distribution of the shaped beam, which determines the temperature field in the piston, must be designed before a diffractive optical element (DOE) can be manufactured. In this paper, a reverse method based on finite element model (FEM) was proposed to design the intensity distribution in order to simulate the thermal loadings on pistons. Temperature fields were obtained by solving a transient three-dimensional heat conduction equation with convective boundary conditions at the surfaces of the piston workpiece. The numerical model then was validated by approaching the computational results to the experimental data. During the process, some important parameters including laser absorptivity, convective heat transfer coefficient, thermal conductivity and Biot number were also validated. Then, optimization procedure was processed to find favorable spatial intensity distribution for the shaped beam, with the aid of the validated FEM. The analysis shows that the reverse method incorporated with numerical simulation can reduce design cycle and design expense efficiently. This method can serve as a kind of virtual experimental vehicle as well, which makes the thermal fatigue simulation test more controllable and predictable. (C) 2007 Elsevier Ltd. All rights reserved.

Finite element analysis of compressor disc heat transfer using the transient conduction method with sparse boundary collocation points

Surface temperature measurements from two discs of a gas turbine compressor rig are used as boundary conditions for the transient conduction solution (inverse heat transfer analysis). The disc geometry is complex, and so the finite element method is used. There are often large radial temperature gradients on the discs, and the equations are therefore solved taking into account the dependence of thermal conductivity on temperature. The solution technique also makes use of a multigrid algorithm to reduce the solution time. This is particularly important since a large amount of data must be analyzed to obtain correlations of the heat transfer. The finite element grid is also used for a network analysis to calculate the radiant heat transfer in the cavity formed between the two compressor discs. The work discussed here proved particularly challenging as the disc temperatures were only measured at four different radial locations. Four methods of surface temperature interpolation are examined, together with their effect on the local heat fluxes. It is found that the choice of interpolation method depends on the available number of data points. Bessel interpolation gives the best results for four data points, whereas cubic splines are preferred when there are considerably more data points. The results from the analysis of the compressor rig data show that the heat transfer near the disc inner radius appears to be influenced by the central throughflow. However, for larger radii, the heat transfer from the discs and peripheral shroud is found to be consistent with that of a buoyancy-induced flow.

Simulations of thermal upgrading methods for oil shale reservoirs

This paper analyzes reaction and thermal front development in porous reservoirs with reacting flows, such as those encountered in shale oil extraction. A set of dimensionless parameters and a 3D code are developed in order to investigate the important physical and chemical variables of such reservoirs when heated by in situ methods. This contribution builds on a 1D model developed for the precursor study to this work. Theory necessary for this study is presented, namely shale decomposition chemical mechanisms, governing equations for multiphase flow in porous media and necessary closure models. Plotting the ratio of the thermal wave speed to the fluid speed allows one to infer that the reaction wave front ends where this ratio is at a minimum. The reaction front follows the thermal front closely, thus allowing assumptions to be made about the extent of decomposition solely by looking at thermal wave progression. Furthermore, this sensitivity analysis showed that a certain minimum permeability is required in order to ensure the formation of a traveling thermal wave. It was found that by studying the non-dimensional governing parameters of the system one can ascribe characteristic values for these parameters for given initial and boundary conditions. This allows one to roughly predict the performance of a particular method on a particular reservoir given approximate values for initial and boundary conditions. Channelling and flow blockage due to carbon residue buildup impeded each method's performance. Blockage was found to be a result of imbalanced heating. Copyright 2012, Society of Petroleum Engineers.

Thermal Deformation Analysis of Micro satellite Structure

Mounting accuracy of satellite payload and ADCS (attitude determination and control subsystem) seats is one of the requirements to achieve the satellite mission with acceptable performance. Components of mounting inaccuracy are technological inaccuracies, residual plastic deformations after loading (during transportation and orbital insertion), elastic deformations, and thermal deformations during orbital operation. This paper focuses on estimation of thermal deformations of satellite structure. Thermal analysis is executed by applying finite-difference method (IDEAS) and temperature profile for satellite components case is evaluated. Then, Perform thermal finite-element analysis applying the finite-difference model results as boundary conditions; and calculate the resultant thermal strain. Next, applying the resultant thermal strain, perform finite-element structure analysis to evaluate structure deformations at the payload and ADCS equipments seats.

Development and Validation of a Compact Thermal Model for an Aircraft Compartment

The development of accurate structural/thermal numerical models of complex systems, such as aircraft fuselage barrels, is often limited and determined by the smallest scales that need to be modelled. The development of reduced order models of the smallest scales and consequently their integration with higher level models can be a way to minimise the bottle neck present, while still having efficient, robust and accurate numerical models. In this paper a methodology on how to develop compact thermal fluid models (CTFMs) for compartments where mixed convection regimes are present is demonstrated. Detailed numerical simulations (CFD) have been developed for an aircraft crown compartment and validated against experimental data obtained from a 1:1 scale compartment rig. The crown compartment is defined as the confined area between the upper fuselage and the passenger cabin in a single aisle commercial aircraft. CFD results were utilised to extract average quantities (temperature and heat fluxes) and characteristic parameters (heat transfer coefficients) to generate CTFMs. The CTFMs have then been compared with the results obtained from the detailed models showing average errors for temperature predictions lower than 5%. This error can be deemed acceptable when compared to the nominal experimental error associated with the thermocouple measurements.

The CTFMs methodology developed allows to generate accurate reduced order models where accuracy is restricted to the region of Boundary Conditions applied. This limitation arises from the sensitivity of the internal flow structures to the applied boundary condition set. CTFMs thus generated can be then integrated in complex numerical modelling of whole fuselage sections.

Further steps in the development of an exhaustive methodology would be the implementation of a logic ruled based approach to extract directly from the CFD simulations numbers and positions of the nodes for the CTFM.

Generation of inertia-gravity waves in the rotating thermal annulus by a localised boundary layer instability

Waves with periods shorter than the inertial period exist in the atmosphere (as inertia-gravity waves) and in the oceans (as Poincaré and internal gravity waves). Such waves owe their origin to various mechanisms, but of particular interest are those arising either from local secondary instabilities or spontaneous emission due to loss of balance. These phenomena have been studied in the laboratory, both in the mechanically-forced and the thermally-forced rotating annulus. Their generation mechanisms, especially in the latter system, have not yet been fully understood, however. Here we examine short period waves in a numerical model of the rotating thermal annulus, and show how the results are consistent with those from earlier laboratory experiments. We then show how these waves are consistent with being inertia-gravity waves generated by a localised instability within the thermal boundary layer, the location of which is determined by regions of strong shear and downwelling at certain points within a large-scale baroclinic wave flow. The resulting instability launches small-scale inertia-gravity waves into the geostrophic interior of the flow. Their behaviour is captured in fully nonlinear numerical simulations in a finite-difference, 3D Boussinesq Navier-Stokes model. Such a mechanism has many similarities with those responsible for launching small- and meso-scale inertia-gravity waves in the atmosphere from fronts and local convection.

Ventilation in a street canyon under diurnal heating conditions

To study the thermal effects on airflow in a street canyon under real heating conditions (due to diurnal solar radiation), a one-way static approach combining an urban canopy model and CFD is proposed in this paper. An urban canopy model was developed to calculate the individual temperatures of surfaces in the street canyon. The calculated surface temperature may be used as a thermal boundary for CFD simulation. The reliability of this model was validated against a field experiment in Harbin, China. Using the coupling calculation method, the wind flow and air exchange process inside an idealized street canyon was studied. The simulation results show that the thermal effect has significant impacts on the transfer process in the street canyon, especially when the approaching wind is weak. Under a real diurnal thermal forcing, the flow structure within the street canyon changes from one primary vortex to two counter-rotating vortices. The change of transfer process, induced by the buoyancy force, was determined by the thermal condition of all surfaces rather than a single one. Key words: thermal effect, street canyon, numerical simulation, transfer process, diurnal heating.

Boundary and internal conditions for adjoint fluid-flow problems

The ever-increasing robustness and reliability of flow-simulation methods have consolidated CFD as a major tool in virtually all branches of fluid mechanics. Traditionally, those methods have played a crucial role in the analysis of flow physics. In more recent years, though, the subject has broadened considerably, with the development of optimization and inverse design applications. Since then, the search for efficient ways to evaluate flow-sensitivity gradients has received the attention of numerous researchers. In this scenario, the adjoint method has emerged as, quite possibly, the most powerful tool for the job, which heightens the need for a clear understanding of its conceptual basis. Yet, some of its underlying aspects are still subject to debate in the literature, despite all the research that has been carried out on the method. Such is the case with the adjoint boundary and internal conditions, in particular. The present work aims to shed more light on that topic, with emphasis on the need for an internal shock condition. By following the path of previous authors, the quasi-1D Euler problem is used as a vehicle to explore those concepts. The results clearly indicate that the behavior of the adjoint solution through a shock wave ultimately depends upon the nature of the objective functional.

Influence of hydraulics and control of thermal storage in solar assisted heat pump combisystems

This paper studies the influence of hydraulics and control of thermal storage in systems combined with solar thermal and heat pump for the production of warm water and space heating in dwellings. A reference air source heat pump system with flat plate collectors connected to a combistore was defined and modeled together with the IEA SHC Task 44 / HPP Annex 38 (T44A38) “Solar and Heat Pump Systems” boundary conditions of Strasbourg climate and SFH45 building. Three and four pipe connections as well as use of internal and external heat exchangers for DHW preparation were investigated as well as sensor height for charging of the DHW zone in the store. The temperature in this zone was varied to ensure the same DHW comfort was achieved in all cases. The results show that the four pipe connection results in 9% improvement in SPF compared to three pipe and that the external heat exchanger for DHW preparation leads to a 2% improvement compared to the reference case. Additionally the sensor height for charging the DHW zone of the store should not be too low, otherwise system performance is adversely affected