PCM thermal insulation in buildings

This paper presents the results of a numerical and experimental study of phase change material (PCM) filled walls and roofs under real operational conditions to achieve passive thermal comfort. The numerical part of the study was based on a one-dimensional model for the phase change problem controlled by pure conduction. Real radiation data was used to determine the external face temperature. The numerical treatment was based upon using finite difference approximations and the ADI scheme. The results obtained were compared with field measurements. The experimental set-up consisted of a small room with movable roof and side wall. The roof was constructed in the traditional way but with the phase change material enclosed. Thermocouples were distributed across the cross section of the roof. Another roof, identical but without the PCM, was also used during comparative tests. The movable wall was also constructed as is done traditionally but with the PCM enclosed. Again, thermocouples were distributed across the wall thickness to enable measurement of the local temperatures. Another wall, identical but without the PCM, was also used during comparative tests. The PCM used in the numerical and experimental tests was composed of a mixture of two commercial grades of glycol in order to obtain the required fusion temperature range. Comparison between the simulation results and the experiments indicated good agreement. Field tests also indicated that the PCM used was adequate and that the concept was effective in maintaining the indoor temperature very close to the established comfort limits. Further economical analysis indicated that the concept could effectively help in reducing the electric energy consumption and improving the energy demand pattern. © 1997 by John Wiley & Sons, Ltd.

Numerical and experimental analysis of natural convection in a cavity with flush mounted heat sources on a side wall

This work considers a problem of interest in several technological applications such as the thermal control of electronic equipment. It is also important to study the heat transfer performance of these components under off-normal conditions, such as during failure of cooling fans. The effect of natural convection on the flow and heat transfer in a cavity with two flush mounted heat sources on the left vertical wall, simulating electronic components, is studied numerically and experimentally. The influence of the power distribution, spacing between the heat sources and cavity aspect ratio have been investigated. An analysis of the average Nusselt number of the two heat sources was performed to investigate the behavior of the heat transfer coefficients. The results obtained numerically and experimentally, after an error analysis, showed a good agreement.

Simulação numérica de trocadores de calor multitubulares equipados com fitas helicoidais externamente aos tubos

Classical shell-and-tube heat exchangers are usually equipped with segmental baffles. These baffles serve two basic functions: (a) they provide tube supports, thereby preventing or reducing mechanical problems, such as sagging or vibration; (b) they direct the fluid flow over the tubes so as to introduce a cross-flow component, thereby increasing the heat transfer. Segmented baffles have several sources of performance loss, some due to various leakage flows and others caused by stagnation zones. A new concept of longitudinal flow heat exchanger - based on placing twisted tapes along the tube bundle subchannels - was developed to mitigate drawbacks of other types of tubular heat exchangers. In this paper, a numerical model has been implemented in order to simulate the thermal-hydraulic feature of tubular heat exchangers equipped either with segmental baffles or with subchannel twisted tapes. The tube bundle has been described by means of an equivalent porous medium type model, allowing a macroscopic description of the shell-side flow. The basic equations - continuity, momentum and energy - have been solved by using the finite volume method. Typical numerical results have been compared with experimental data, reaching a very good agreement. A comparative analysis of different types of heat exchangers has been carried out, revealing the satisfactory thermal-hydraulic efficiency level of the twisted tapes heat exchangers.

Study on transient ice formation of laminar flow inside externally supercooled rectangular duct

The transient process of solidification of laminar liquid flow (water) submitted to super-cooling was investigated both theoretically and experimentally. In this study an alternative analytical formulation and numerical approach were adopted resulting in the unsteady model with temperature dependent thermophysical properties in the solid region. The proposed model is based upon the fundamental equations of energy balance in the solid and liquid regions as well as across the solidification front. The basic equations and the associated boundary and initial conditions were made dimensionless by using the Landau transformation to immobilize the moving front and render the problem to a fixed plane type problem. A laminar velocity profile is admitted in the liquid domain and the resulting equations were discretized using the finite difference approach. The numerical predictions obtained were compared with the available results based on other models and concepts such as Neumann analytical model, the apparent thermal capacity model due to Bonacina and the conventional fixed grid energy model due to Goodrich. To obtain further comparisons and more validation of the model and the numerical solution, an experimental rig was constructed and instrumented permitting very well controlled experimental measurements. The numerical predictions were compared with the experimental results and the agreement was found satisfactory.

Optimal design of smart structures using bonded piezoelectrics for vibration control

Smart material technology has become an area of increasing interest for the development of lighter and stronger structures which are able to incorporate actuator and sensor capabilities for collocated control. In the design of actively controlled structures, the determination of the actuator locations and the controller gains, is a very important issue. For that purpose, smart material modelling, modal analysis methods, control and optimization techniques are the most important ingredients to be taken into account. The optimization problem to be solved in this context presents two interdependent aspects. The first one is related to the discrete optimal actuator location selection problem, which is solved in this paper using genetic algorithms. The second is represented by a continuous variable optimization problem, through which the control gains are determined using classical techniques. A cantilever Euler-Bernoulli beam is used to illustrate the presented methodology.

A FEM approach to the equations of magneto-aerodynamics

In this work, a Finite Element Method treatment is outlined for the equations of Magnetoaerodynamics. In order to provide a good basis for numerical treatment of Magneto-aerodynamics, a full version of the complete equations is presented and FEM contribution matrices are deduced, as well as further terms of stabilization for the compressible flow case.

Pressure field in the cathode-anode region of a high-power klystron amplifier

The pressure field of a high-power klystron amplifier in the cathode and anode region was investigated. The investigation was performed using a 1.3 GHz, 100 A and 240 kV high-power klystron with five reentrant coaxial cavities, assembled in cylindrical drift tube 1.2 m long. The diffusion equation in mathematical model was also solved by using a 3-D finite element method code, in order to obtain pressure profile in region of interest. The results show that density profile of molecules between cathode-anode region was determined, where cathode pressure is approximately 10% higher than anode pressure.

A 2-D delaunay refinement algorithm using an initial prerefinement from the boundary mesh

An improvement to the quality bidimensional Delaunay mesh generation algorithm, which combines the mesh refinement algorithms strategy of Ruppert and Shewchuk is proposed in this research. The developed technique uses diametral lenses criterion, introduced by L. P. Chew, with the purpose of eliminating the extremely obtuse triangles in the boundary mesh. This method splits the boundary segment and obtains an initial prerefinement, and thus reducing the number of necessary iterations to generate a high quality sequential triangulation. Moreover, it decreases the intensity of the communication and synchronization between subdomains in parallel mesh refinement. © 2008 IEEE.

Análise da execução de um algoritmo numérico cbs para simulação de escoamento utilizando elementos finitos

This work intends to analyze the application and execution time of a numerical algorithm that simulates incompressible and isothermal flows. It was used the explicit scheme of the Characteristic Based Split (CBS) algorithm and the Artificial Compressibility (AC) scheme for coupling pressure-velocity equations. The discretization was done with the finite element method using a bilinear elements grid. The free software GNU-Octave was used for implementation and execution of routines. The results were analyzed using the classic lid-driven cavity problem. This work shows results for tests with several Reynolds' number. The results for these tests show a good agreement when compared with previous ones obtained from bibliography. The code runtime's analysis shows yet that the matrix's assembly is the part of greater consumption time in the implementation.

Coupling rate measurement of a novel multi-actuated piezoelectric device using optical interferometry

The class of piezoelectric actuators considered in this paper consists of a multi-flexible structure actuated by two or more piezoceramic devices that must generate different output displacements and forces at different specified points of the domain and in different directions. The devices were modeled by finite element using the software ANSYS and the topology optimization method. The following XY actuators were build to achieve maximum displacement in the X and Y directions with a minimum crosstalk between them. The actuator prototypes are composed of an aluminum structure, manufactured by using a wire Electrical Discharge Machining, which are bonded to rectangular PZT5A piezoceramic blocks by using epoxy resin. Multi-actuator piezoelectric device displacements can be measured by using optical interferometry, since it allows dynamic measurements in the kHz range, which is of the order of the first resonance frequency of these piezomechanisms. A Michelson-type interferometer, with a He-Ne laser source, is used to measure the displacement amplitudes in nanometric range. A new optical phase demodulation technique is applied, based on the properties of the triangular waveform drive voltage applied to the XY piezoelectric nanopositioner. This is a low-phase-modulation-depth-like technique that allows the rapid interferometer auto-calibration. The measurements were performed at 100 Hz frequency, and revealed that the device is linear voltage range utilized in this work. The ratio between the generated and coupled output displacements and the drive voltages is equal to 10.97 nm/V and 1.76 nm/V, respectively, which corresponds to a 16% coupling rate. © 2010 IEEE.

Influence of high insertion torque on implant placement - an anisotropic bone stress analysis

The aim of this study was to evaluate the influence of the high values of insertion torques on the stress and strain distribution in cortical and cancellous bones. Based on tomography imaging, a representative mathematical model of a partial maxilla was built using Mimics 11.11 and Solid Works 2010 softwares. Six models were built and each of them received an implant with one of the following insertion torques: 30, 40, 50, 60, 70 or 80 Ncm on the external hexagon. The cortical and cancellous bones were considered anisotropic. The bone/implant interface was considered perfectly bonded. The numerical analysis was carried out using Ansys Workbench 10.0. The convergence of analysis (6%) drove the mesh refinement. Maximum principal stress (σ max) and maximum principal strain (ε max) were obtained for cortical and cancellous bones around to implant. Pearson's correlation test was used to determine the correlation between insertion torque and stress concentration in the periimplant bone tissue, considering the significance level at 5%. The increase in the insertion torque generated an increase in the σ max and ε max values for cortical and cancellous bone. The σmax was smaller for the cancellous bone, with greater stress variation among the insertion torques. The ε max was higher in the cancellous bone in comparison to the cortical bone. According to the methodology used and the limits of this study, it can be concluded that higher insertion torques increased tensile and compressive stress concentrations in the periimplant bone tissue.

Dependência estrutural da capacidade de rotação plástic em vigas de concreto armado

Summary In this work the structural dependence of plastic rotation capacity in RC beams is evaluated using the Finite Element Method. The objective is to achieve a better understanding of the non-linear behavior of reinforced concrete members and perform extensive parameter studies, using a rational model developed by Bigaj [1] to analyze the phenomenon of plastic rotation capacity in reinforced concrete members. It is assumed that only bending failure is relevant due to sufficient member resistance against shear and torsion. The paper begins with the physical and theoretical background of the phenomenon of plastic hinge development in RC structures. Special emphasis is laid on the issue of structural dependence of deformation capacity of plastic hinges in RC members. Member size dependence and influence of properties of construction materials were emphasized as well. The essential components of the Bigajs model for calculating the plastic rotation capacity are discussed. The behaviour of the plastic hinge is analysed taking into account the strain localisation in the damage zones of the hinge region. The Fictitious Crack Model (FCM) and the Compressive Damage Zone Model (CDZ) are adopted in a Fracture Mechanics approach to model the behaviour of concrete in tension and compression, respectively. The approach is implemented in FEMOOP, a FEM in-house solver under development, and applied to evaluate ductility in 2D beams. The models were generated with GiD, a pre-processor and post-processor developed by CIMNE, and analyzed with the capabilities implemented in FEMOOP. © Universitat Politècnica de Catalunya, Barcelona, España 2010.

Integrating open source codes for solid modeling and automatic mesh generation in FEM applications

The applications of the Finite Element Method (FEM) for three-dimensional domains are already well documented in the framework of Computational Electromagnetics. However, despite the power and reliability of this technique for solving partial differential equations, there are only a few examples of open source codes available and dedicated to the solid modeling and automatic constrained tetrahedralization, which are the most time consuming steps in a typical three-dimensional FEM simulation. Besides, these open source codes are usually developed separately by distinct software teams, and even under conflicting specifications. In this paper, we describe an experiment of open source code integration for solid modeling and automatic mesh generation. The integration strategy and techniques are discussed, and examples and performance results are given, specially for complicated and irregular volumes which are not simply connected. © 2011 IEEE.

Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

Influence of convergence angle of tooth preparation on the fracture resistance of Y-TZP-based all-ceramic restorations

Objective: To investigate the influence of the convergence angle of tooth preparation on the fracture load of Y-TZP-based ceramic (YZ-Vita YZ) substructure (SB) veneered with a feldspathic porcelain (VM9-Vita VM9). Methods: Finite element stress analysis (FEA) was performed to examine the stress distribution of the system. Eighty YZ SB were fabricated using a CAD-CAM system and divided into four groups (n = 20), according to the total occlusal convergence (TOC) angle: G6-6° TOC; G12-12° TOC; G20-20° TOC; and G20MOD-20° TOC with modified SB. All SB were veneered with VM9, cemented in a fiber reinforced epoxy resin die, and loaded to failure. Half of the specimens from each group (n = 10) were cyclic fatigued (106 cycles) before testing. Failure analysis was performed to determine the fracture origin. Data were statistically analyzed using Anova and Tukey's tests (α = 0.05). Results: The greatest mean load to fracture value was found for the G20MOD, which was predicted by the FEA. Cyclic fatigue did not significantly affect the load of fracture. Catastrophic failure originating from the internal occlusal surface of the SB was the predominant failure mode, except for G20MOD. Significance: The YZ-VM9 restorations resisted greater compression load than the usual physiological occlusal load, regardless of the TOC angle of preparations. Yet, the G20MOD design produced the best performance among the experimental conditions evaluated. © 2013 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.