Search results

Present study aims to extend the lattice Boltzmann method (LBM) to simulate radiation in geometries with curved boundaries, as the first step to simulate radiation in complex porous media. In recent years, researchers have increasingly explored the use of porous media to improve the heat transfer processes. The lattice Boltzmann method (LBM) is one of the most effective techniques for simulating heat transfer in such media. However, the application of the LBM to study radiation in complex geometries that contain curved boundaries, as found in many porous media, has been limited.

The numerical evaluation of the effect of the radiation-conduction parameter and extinction coefficient on temperature and incident radiation distributions demonstrates that the proposed LBM algorithm provides highly accurate results across all cases, compared to those found in the literature or those obtained using the finite volume method (FVM) with the discrete ordinates method (DOM) for radiative information.

For the case with a conduction-radiation parameter equal to 0.01, the maximum relative error is 1.9% in predicting temperature along vertical central line. The accuracy improves with an increase in the conduction-radiation parameter. Furthermore, the comparison between computational performances of two approaches reveals that the LBM-LBM approach performs significantly faster than the FVM-DOM solver.

The difficulty of radiative modeling in combined problems involving irregular boundaries has led to alternative approaches that generally increase the computational expense to obtain necessary radiative details. To address the limitations of existing methods, this study presents a new approach involving a coupled lattice Boltzmann and first-order blocked-off technique to efficiently model conductive-radiative heat transfer in complex geometries with participating media. This algorithm has been developed using the parallel lattice Boltzmann solver.

3D steady heat conduction analysis considering heat source is conducted on the fundamental of the fast multipole method (FMM)-accelerated line integration boundary element method (LIBEM).

Due to considering the heat source, domain integral is generated in the traditional heat conduction boundary integral equation (BIE), which will counteract the well-known merit of the BEM, namely, boundary-only discretization. To avoid volume discretization, the enhanced BEM, the LIBEM with dimension reduction property is introduced to transfer the domain integral into line integrals. Besides, owing to the unsatisfactory performance of the LIBEM when it comes to large-scale structures requiring massive computation, the FMM-accelerated LIBEM (FM-LIBEM) is proposed to improve the computation efficiency further.

Assuming N and M are the numbers of nodes and integral lines, respectively, the FM-LIBEM can reduce the time complexity from O(NM) to about O(N+ M), and a full discussion and verification of the advantage are done based on numerical examples under heat conduction.

(1) The LIBEM is applied to 3D heat conduction analysis with heat source. (2) The domain integrals can be transformed into boundary integrals with straight line integrals by the LIM. (3) A FM-LIBEM is proposed and can reduce the time complexity from O(NM) to O(N+ M). (4) The FM-LIBEM with high computational efficiency is exerted to solve 3D heat conduction analysis with heat source in massive computation successfully.

The purpose of this paper is to address unsteady heat conduction in two subsets of ordinary bodies. One subset consists of a large plane wall, a long cylinder and a sphere in one dimension. The other subset consists of a short cylinder and a large rectangular bar in two dimensions. The prevalent assumptions in the two subsets are: constant initial temperature, uniform surface heat flux and thermo-physical properties invariant with temperature. The engineering applications of the unsteady heat conduction deal with the determination of temperature–time histories in the two subsets using electric resistance heating, radiative heating and fire pool heating.

To this end, a novel numerical procedure named the enhanced method of discretization in time (EMDT) transforms the linear one-dimensional unsteady, heat conduction equations with non-homogeneous boundary conditions into equivalent nonlinear “quasi–steady” heat conduction equations having the time variable embedded as a time parameter. The equivalent nonlinear “quasi–steady” heat conduction equations are solved with a finite difference method.

Based on the numerical computations, it is demonstrated that the approximate temperature–time histories in the simple subset of ordinary bodies (large plane wall, long cylinder and sphere) exhibit a perfect matching over the entire time domain 0 < t < ∞ when compared against the rigorous exact temperature–time histories expressed by classical infinite series. Furthermore, using the method of superposition of solutions in the convoluted subset (short cylinder and large rectangular crossbar), the same level of agreement in the approximate temperature–time histories in the simple subset of ordinary bodies is evident.

The performance of the proposed EMDT coupled with a finite difference method is exhaustively assessed in the solution of the unsteady, one-dimensional heat conduction equations with prescribed surface heat flux for: a subset of one-dimensional bodies (plane wall, long cylinder and spheres) and a subset of two-dimensional bodies (short cylinder and large rectangular bar).

Gas insulated systems, such as gas insulated lines (GIL), use insulating gas, mostly sulfur hexalfluoride (SF₆), to enable a higher dielectric strength compared to e.g. air. However, under high voltage direct current conditions, charge accumulation and electric field stress may occur, which may lead to partial discharge or system failure. Therefore, numerical simulations are used to design the system and determine the electric field and charge distribution. Although the gas conduction shows a more complex current–voltage characteristic compared to solid insulation, the electric conductivity of the SF₆ gas is set as constant in most works. The purpose of this study is to investigate different approaches to address the conduction in the gas properly for numerical simulations.

In this work, two approaches are investigated to address the conduction in the insulating gas and are compared to each other. One method is an ion-drift-diffusion model, where the conduction in the gas is described by the ion motion in the SF₆ gas. However, this method is computationally expensive. Alternatively, a less complex approach is an electro-thermal model with the application of an electric conductivity model for the SF₆ gas. Measurements show that the electric conductivity in the SF₆ gas has a nonlinear dependency on temperature, electric field and gas pressure. From these measurements, an electric conductivity model was developed. Both methods are compared by simulation results, where different parameters and conditions are considered, to investigate the potential of the electric conductivity model as a computationally less expensive alternative.

The simulation results of both simulation approaches show similar results, proving the electric conductivity for the SF₆ gas as a valid alternative. Using the electro-thermal model approach with the application of the electric conductivity model enables a solution time up to six times faster compared to the ion-drift-diffusion model. The application of the model allows to examine the influence of different parameters such as temperature and gas pressure on the electric field distribution in the GIL, whereas the ion-drift-diffusion model enables to investigate the distribution of homo- and heteropolar charges in the insulation gas.

This work presents numerical simulation models for high voltage direct current GIL, where the conduction in the SF₆ gas is described more precisely compared to a definition of a constant electric conductivity value for the insulation gas. The electric conductivity model for the SF₆ gas allows for consideration of the current–voltage characteristics of the gas, is computationally less expensive compared to an ion-drift diffusion model and needs considerably less solution time.

In the present work, we focus on developing an in-house parallel meshless local Petrov-Galerkin (MLPG) code for the analysis of heat conduction in two-dimensional and three-dimensional regular as well as complex geometries.

The parallel MLPG code has been implemented using open multi-processing (OpenMP) application programming interface (API) on the shared memory multicore CPU architecture. Numerical simulations have been performed to find the critical regions of the serial code, and an OpenMP-based parallel MLPG code is developed, considering the critical regions of the sequential code.

Based on performance parameters such as speed-up and parallel efficiency, the credibility of the parallelization procedure has been established. Maximum speed-up and parallel efficiency are 10.94 and 0.92 for regular three-dimensional geometry (343,000 nodes). Results demonstrate the suitability of parallelization for larger nodes as parallel efficiency and speed-up are more for the larger nodes.

Few attempts have been made in parallel implementation of the MLPG method for solving large-scale industrial problems. Although the literature suggests that message-passing interface (MPI) based parallel MLPG codes have been developed, the OpenMP model has rarely been touched. This work is an attempt at the development of OpenMP-based parallel MLPG code for the very first time.

A few earlier studies presented infeasible heatline trajectories for natural convection within annular domains involving an inner circular cylinder and outer square/circular enclosure. The purpose of this paper is to revisit and illustrate the correct heatline trajectories for various test cases.

Galerkin finite element based methodology and space adaptive grid have been used to simulate natural convective flows within the annular domains. The prediction of heatlines involves derivatives at the nodes, which are evaluated based on finite element basis functions and contributions from neighboring elements.

The heatlines in the earlier work indicate infeasible heat flow paths such as heat flow from one portion to the other of isothermal hot walls and heat flow across the adiabatic walls. Current results illustrate physically consistent heat flow paths involving perpendicularly emerging heatlines from hot to cold walls for conductive transport, long heat flow paths around the closed-loop heatline cells for convective transport and parallel layout of heatlines to the adiabatic walls. Results also demonstrate complex heatlines involving multiple flow vortices and complex flow structures.

Current work translates heatfunctions from energy flux vectors, which are determined by using basis sets. This work demonstrates the expected heatline trajectories for various scenarios involving conductive and convective heat transport within enclosures with an inner hot object as a first attempt, and the results are precursors for the understanding of energy flow estimates.

The purpose of this study is to clarify the mechanism of ambient and transient temperature effects on piezoelectric pressure sensors, and to propose corresponding compensation measures. The temperature of the explosion field has a significant influence on the piezoelectric sensor used to measure the shock wave pressure. For accurate shock wave pressure measurement, based on the actual piezoelectric pressure sensors used in the explosion field, the effects of ambient and transient temperatures on the sensor should be studied.

The compensation method of the ambient temperature is discussed according to the sensor size and material. The theoretical analysis method of the transient temperature is proposed. For the transient temperature conduction problem of the sensor, the finite element simulation method of structure-temperature coupling is used to solve the temperature distribution of the sensor and the change in the contact force on the quartz crystal surface under the step and triangle temperatures. The simulation results are highly consistent with the theory.

Based on the analysis results, a transient temperature control method is proposed, in which 0.5 mm thick lubricating silicone grease is applied to the sensor diaphragm, and 0.2 mm thick fiberglass cloth is wrapped around the sensor side. Simulation experiments are carried out to verify the feasibility of the control method, and the results show that the control method effectively suppresses the output of the thermal parasitic.

The above thermal protection methods can effectively improve the measurement accuracy of shock wave pressure and provide technical support for the evaluation of the power of explosion damage.

This study aims to investigate a new circuitous minichannel cold plate (MCP) design involving flow fragmentation. The overall thermal performance and the temperature uniformity analysis are performed and compared with the traditional serpentine design. The substrate thickness and its thermal conductivity are varied to analyse the effect of axial-back conduction due to the conjugate nature of heat transfer.

The traditional serpentine minichannel is modified into five new fragmented designs with two inlets and two outlets. A three-dimensional numerical model involving the effect of conjugate heat transfer with a single-phase laminar fluid flow subjected to constant heat flux is solved using a finite volume-based computational fluid dynamics solver.

The minimum and maximum temperature differences are observed for the two branch fragmented flow designs. The two-branch and middle channel fragmented design shows better temperature uniformity over other designs while the three-branch fragmented designs exhibited better hydrodynamic performance.

MCPs could be used as an indirect liquid cooling method for battery thermal management of pouch and prismatic cells. Coupling the modified cold plates with a battery module and investigating the effect of different battery parameters and environmental effects in a transient state are the prospects for further research.

The study involves several aspects of evaluation for a conclusive decision on optimum channel design by analysing the performance plot between the temperature uniformity index, average base temperature and overall thermal performance. The new fragmented channels are designed in a way to facilitate the fluid towards the outlet in the minimum possible path thereby reducing the pressure drop, also maximizing the heat transfer and temperature uniformity from the substrate due to two inlets and a reversed-flow pattern. Simplified minichannel designs are proposed in this study for practical deployment and ease of manufacturability.

Multilevel inverters play a major role in the development of high-power industrial applications. In traditional low-level inverters (e.g. 2-level), the switching frequency is restricted and the harmonic spectrum of the system is hard to meet power requirements. Similarly, high-level inverters consist of a large number of switches, complex modulation techniques and complex hardware architecture, which results in high power loss and a significant amount of harmonic distortion. Furthermore, it is a must to ensure that every switch experiences the same stress of voltage and current. The purpose of this paper is to present an inverter topology with lower conduction and switching losses via reduced number of switches and equal voltage source-sharing technique.

Herein, the authors present a cascaded multilevel inverter having less power switches, a simple modulation technique and an equal voltage source-sharing phenomenon implementation.

The modulation technique becomes more complex when equal voltage source-sharing is to be implemented. In this study, a novel topology for the multilevel inverter with fewer switches, novel modulation technique, equal voltage source-sharing and Inductor-Capacitor-Inductor filter implementation is demonstrated to the reduce harmonic spectrum and power losses of the proposed system.

The nine-level inverter design is validated using software simulations and hardware prototype testing; the power losses of the proposed inverter design are elaborated and compared with the traditional approach.

For a thermal protection system (TPS) of long endurance hypersonic flight vehicle (HFV), its thermal insulation property not only determines by the manufactured morphology but also changes along time. A thermal conductivity prediction model for aerogel considering heat treatment effect is carried out and applied to solve the heat conduction problem of a TPS. The aim of this study is to provide theoretical and numerical references for further development of aerogels applying to TPSs.

A thermal conductivity prediction model for aerogel is established considering treatment effect. The heat conduction problem of a TPS is derived and solved by combining the differential quadrature method and the Runge–Kutta method. The prediction results of aerogel thermal conductivities are verified by comparing with those in literature, while the calculated temperature field of TPS is verified by comparing with that by ABAQUS.

Numerical results show that when applying the current prediction model, the calculated high temperature area in the aerogel layer is narrowed due to the decrease of the thermal conductivity during heat treatment process.

This study will be beneficial to carry out the precise design of TPS for long endurance HFVs.