Advanced thermal analysis of microelectronics using spreading resistance models

Thermal analysis of electronic devices is one of the most important steps for designing of modern devices. Precise thermal analysis is essential for designing an effective thermal management system of modern electronic devices such as batteries, LEDs, microelectronics, ICs, circuit boards, semiconductors and heat spreaders. For having a precise thermal analysis, the temperature profile and thermal spreading resistance of the device should be calculated by considering the geometry, property and boundary conditions. Thermal spreading resistance occurs when heat enters through a portion of a surface and flows by conduction. It is the primary source of thermal resistance when heat flows from a tiny heat source to a thin and wide heat spreader. In this thesis, analytical models for modeling the temperature behavior and thermal resistance in some common geometries of microelectronic devices such as heat channels and heat tubes are investigated. Different boundary conditions for the system are considered. Along the source plane, a combination of discretely specified heat flux, specified temperatures and adiabatic condition are studied. Along the walls of the system, adiabatic or convective cooling boundary conditions are assumed. Along the sink plane, convective cooling with constant or variable heat transfer coefficient are considered. Also, the effect of orthotropic properties is discussed. This thesis contains nine chapters. Chapter one is the introduction and shows the concepts of thermal spreading resistance besides the originality and importance of the work. Chapter two reviews the literatures on the thermal spreading resistance in the past fifty years with a focus on the recent advances. In chapters three and four, thermal resistance of a twodimensional flux channel with non-uniform convection coefficient in the heat sink plane is studied. The non-uniform convection is modeled by using two functions than can simulate a wide variety of different heat sink configurations. In chapter five, a non-symmetrical flux channel with different heat transfer coefficient along the right and left edges and sink plane is analytically modeled. Due to the edge cooling and non-symmetry, the eigenvalues of the system are defined using the heat transfer coefficient on both edges and for satisfying the orthogonality condition, a normalized function is calculated. In chapter six, thermal behavior of two-dimensional rectangular flux channel with arbitrary boundary conditions on the source plane is presented. The boundary condition along the source plane can be a combination of the first kind boundary condition (Dirichlet or prescribed temperature) and the second kind boundary condition (Neumann or prescribed heat flux). The proposed solution can be used for modeling the flux channels with numerous different source plane boundary conditions without any limitations in the number and position of heat sources. In chapter seven, temperature profile of a circular flux tube with discretely specified boundary conditions along the source plane is presented. Also, the effect of orthotropic properties are discussed. In chapter 8, a three-dimensional rectangular flux channel with a non-uniform heat convection along the heat sink plane is analytically modeled. In chapter nine, a summary of the achievements is presented and some systems are proposed for the future studies. It is worth mentioning that all the models and case studies in the thesis are compared with the Finite Element Method (FEM).

Advanced Aerogel Composites for Oil Remediation and Recovery

Oil spills in marine environments often damage marine and coastal life if not remediated rapidly and efficiently. In spite of the strict enforcement of environmental legislations (i.e., Oil Pollution Act 1990) following the Exxon Valdez oil spill (June 1989; the second biggest oil spill in U.S. history), the Macondo well blowout disaster (April 2010) released 18 times more oil. Strikingly, the response methods used to contain and capture spilled oil after both accidents were nearly identical, note that more than two decades separate Exxon Valdez (1989) and Macondo well (2010) accidents.

The goal of this dissertation was to investigate new advanced materials (mechanically strong aerogel composite blankets-Cabot® Thermal Wrap™ (TW) and Aspen Aerogels® Spaceloft® (SL)), and their applications for oil capture and recovery to overcome the current material limitations in oil spill response methods. First, uptake of different solvents and oils were studied to answer the following question: do these blanket aerogel composites have competitive oil uptake compared to state-of-the-art oil sorbents (i.e., polyurethane foam-PUF)? In addition to their competitive mechanical strength (766, 380, 92 kPa for Spaceloft, Thermal Wrap, and PUF, respectively), our results showed that aerogel composites have three critical advantages over PUF: rapid (3-5 min.) and high (more than two times of PUF’s uptake) oil uptake, reusability (over 10 cycles), and oil recoverability (up to 60%) via mechanical extraction. Chemical-specific sorption experiments showed that the dominant uptake mechanism of aerogels is adsorption to the internal surface, with some contribution of absorption into the pore space.

Second, we investigated the potential environmental impacts (energy and chemical burdens) associated with manufacturing, use, and disposal of SL aerogel and PUF to remove the oil (i.e., 1 m3 oil) from a location (i.e., Macondo well). Different use (single and multiple use) and end of life (landfill, incinerator, and waste-to-energy) scenarios were assessed, and our results demonstrated that multiple use, and waste-to-energy choices minimize the energy and material use of SL aerogel. Nevertheless, using SL once and disposing via landfill still offers environmental and cost savings benefits relative to PUF, and so these benefits are preserved irrespective of the oil-spill-response operator choices.

To inform future aerogel manufacture, we investigated the different laboratory-scale aerogel fabrication technologies (rapid supercritical extraction (RSCE), CO2 supercritical extraction (CSCE), alcohol supercritical extraction (ASCE)). Our results from anticipatory LCA for laboratory-scaled aerogel fabrication demonstrated that RSCE method offers lower cumulative energy and ecotoxicity impacts compared to conventional aerogel fabrication methods (CSCE and ASCE).

The final objective of this study was to investigate different surface coating techniques to enhance oil recovery by modifying the existing aerogel surface chemistries to develop chemically responsive materials (switchable hydrophobicity in response to a CO2 stimulus). Our results showed that studied surface coating methods (drop casting, dip coating, and physical vapor deposition) were partially successful to modify surface with CO2 switchable chemical (tributylpentanamidine), likely because of the heterogeneous fiber structure of the aerogel blankets. A possible solution to these non-uniform coatings would be to include switchable chemical as a precursor during the gel preparation to chemically attach the switchable chemical to the pores of the aerogel.

Taken as a whole, the implications of this work are that mechanical deployment and recovery of aerogel composite blankets is a viable oil spill response strategy that can be deployed today. This will ultimately enable better oil uptake without the uptake of water, potential reuse of the collected oil, reduced material and energy burdens compared to competitive sorbents (e.g., PUF), and reduced occupational exposure to oiled sorbents. In addition, sorbent blankets and booms could be deployed in coastal and open-ocean settings, respectively, which was previously impossible.

Stable carbon isotope ratios of carbon dioxide from EDC and Berkner Island ice cores for 40-50 ka BP

The stable carbon isotopic signature of carbon dioxide (d13CO2) measured in the air occlusions of polar ice provides important constraints on the carbon cycle in past climates. In order to exploit this information for previous glacial periods, one must use deep, clathrated ice, where the occluded air is preserved not in bubbles but in the form of air hydrates. Therefore, it must be established whether the original atmospheric d13CO2 signature can be reconstructed from clathrated ice. We present a comparative study using coeval bubbly ice from Berkner Island and ice from the bubble-clathrate transformation zone (BCTZ) of EPICA Dome C (EDC). In the EDC samples the gas is partitioned into clathrates and remaining bubbles as shown by erroneously low and scattered CO2 concentration values, presenting a worst-case test for d13CO2 reconstructions. Even so, the reconstructed atmospheric d13CO2 values show only slightly larger scatter. The difference to data from coeval bubbly ice is statistically significant. However, the 0.16 per mil magnitude of the offset is small for practical purposes, especially in light of uncertainty from non-uniform corrections for diffusion related fractionation that could contribute to the discrepancy. Our results are promising for palaeo-atmospheric studies of d13CO2 using a ball mill dry extraction technique below the BCTZ of ice cores, where gas is not subject to fractionation into microfractures and between clathrate and bubble reservoirs.

Particle positions of numerical experiments 1 to 4

In mixed sediment beds, erosion resistance can change relative to that of beds composed of a uniform sediment because of varying textural and/or other grain-size parameters, with effects on pore water flow that are difficult to quantify by means of analogue techniques. To overcome this difficulty, a three-dimensional numerical model was developed using a finite difference method (FDM) flow model coupled with a distinct element method (DEM) particle model. The main aim was to investigate, at a high spatial resolution, the physical processes occurring during the initiation of motion of single grains at the sediment-water interface and in the shallow subsurface of simplified sediment beds under different flow velocities. Increasing proportions of very fine sand (D50=0.08 mm) were mixed into a coarse sand matrix (D50=0.6 mm) to simulate mixed sediment beds, starting with a pure coarse sand bed in experiment 1 (0 wt% fines), and proceeding through experiment 2 (6.5 wt% fines), experiment 3 (10.5 wt% fines), and experiment 4 (28.7 wt% fines). All mixed beds were tested for their erosion behavior at predefined flow velocities varying in the range of U 1-5=10-30 cm/s. The experiments show that, with increasing fine content, the smaller particles increasingly fill the spaces between the larger particles. As a consequence, pore water inflow into the sediment is increasingly blocked, i.e., there is a decrease in pore water flow velocity and, hence, in the flow momentum available to entrain particles. These findings are portrayed in a new conceptual model of enhanced sediment bed stabilization.

Multilayered graphene films prepared at moderate temperatures using energetic physical vapour deposition

Carbon films were energetically deposited onto copper and nickel foil using a filtered cathodic vacuum arc deposition system. Raman spectroscopy, scanning electron microscopy, transmission electron microscopy and UV–visible spectroscopy showed that graphene films of uniform thickness with up to 10 layers can be deposited onto copper foil at moderate temperatures of 750 C. The resulting films, which can be prepared at high deposition rates, were comparable to graphene films grown at 1050 C using chemical vapour deposition (CVD). This difference in growth temperature is attributed to dynamic annealing which occurs as the film grows from the energetic carbon flux. In the case of nickel substrates, it was found that graphene films can also be prepared at moderate substrate temperatures. However much higher carbon doses were required, indicating that the growth mode differs between substrates as observed in CVD grown graphene. The films deposited onto nickel were also highly non uniform in thickness, indicating that the grain structure of the nickel substrate influenced the growth of graphene layers. 

Design and development of a diagonal-airflow batch dryer for spatial drying homogeneity

Das Verfahren der Lebensmitteltrocknung wird häufig angewendet, um ein Produkt für längere Zeit haltbar zu machen. Obst und Gemüse sind aufgrund ihres hohen Wassergehalts leicht verderblich durch biochemische Vorgänge innerhalb des Produktes, nicht sachgemäße Lagerung und unzureichende Transportmöglichkeiten. Um solche Verluste zu vermeiden wird die direkte Trocknung eingesetzt, welche die älteste Methode zum langfristigen haltbarmachen ist. Diese Methode ist jedoch veraltet und kann den heutigen Herausforderungen nicht gerecht werden. In der vorliegenden Arbeit wurde ein neuer Chargentrockner, mit diagonalem Luftstömungskanal entlang der Länge des Trocknungsraumes und ohne Leitbleche entwickelt. Neben dem unbestreitbaren Nutzen der Verwendung von Leitblechen, erhöhen diese jedoch die Konstruktionskosten und führen auch zu einer Erhöhung des Druckverlustes. Dadurch wird im Trocknungsprozess mehr Energie verbraucht. Um eine räumlich gleichmäßige Trocknung ohne Leitbleche zu erreichen, wurden die Lebensmittelbehälter diagonal entlang der Länge des Trockners platziert. Das vorrangige Ziel des diagonalen Kanals war, die einströmende, warme Luft gleichmäßig auf das gesamte Produkt auszurichten. Die Simulation des Luftstroms wurde mit ANSYS-Fluent in der ANSYS Workbench Plattform durchgeführt. Zwei verschiedene Geometrien der Trocknungskammer, diagonal und nicht diagonal, wurden modelliert und die Ergebnisse für eine gleichmäßige Luftverteilung aus dem diagonalen Luftströmungsdesign erhalten. Es wurde eine Reihe von Experimenten durchgeführt, um das Design zu bewerten. Kartoffelscheiben dienten als Trocknungsgut. Die statistischen Ergebnisse zeigen einen guten Korrelationskoeffizienten für die Luftstromverteilung (87,09%) zwischen dem durchschnittlich vorhergesagten und der durchschnittlichen gemessenen Strömungsgeschwindigkeit. Um den Effekt der gleichmäßigen Luftverteilung auf die Veränderung der Qualität zu bewerten, wurde die Farbe des Produktes, entlang der gesamten Länge der Trocknungskammer kontaktfrei im on-line-Verfahren bestimmt. Zu diesem Zweck wurde eine Imaging-Box, bestehend aus Kamera und Beleuchtung entwickelt. Räumliche Unterschiede dieses Qualitätsparameters wurden als Kriterium gewählt, um die gleichmäßige Trocknungsqualität in der Trocknungskammer zu bewerten. Entscheidend beim Lebensmittel-Chargentrockner ist sein Energieverbrauch. Dafür wurden thermodynamische Analysen des Trockners durchgeführt. Die Energieeffizienz des Systems wurde unter den gewählten Trocknungsbedingungen mit 50,16% kalkuliert. Die durchschnittlich genutzten Energie in Form von Elektrizität zur Herstellung von 1kg getrockneter Kartoffeln wurde mit weniger als 16,24 MJ/kg und weniger als 4,78 MJ/kg Wasser zum verdampfen bei einer sehr hohen Temperatur von jeweils 65°C und Scheibendicken von 5mm kalkuliert. Die Energie- und Exergieanalysen für diagonale Chargentrockner wurden zudem mit denen anderer Chargentrockner verglichen. Die Auswahl von Trocknungstemperatur, Massenflussrate der Trocknungsluft, Trocknerkapazität und Heiztyp sind die wichtigen Parameter zur Bewertung der genutzten Energie von Chargentrocknern. Die Entwicklung des diagonalen Chargentrockners ist eine nützliche und effektive Möglichkeit um dei Trocknungshomogenität zu erhöhen. Das Design erlaubt es, das gesamte Produkt in der Trocknungskammer gleichmäßigen Luftverhältnissen auszusetzen, statt die Luft von einer Horde zur nächsten zu leiten.

Dimensionality hyper-reduction and machine learning for dynamical systems with varying parameters

Thesis (Ph.D.)--University of Washington, 2016-08

A Study on the liquefaction risk in seismic design of foundations

A fully coupled non-linear effective stress response finite difference (FD) model is built to survey the counter-intuitive recent findings on the reliance of pore water pressure ratio on foundation contact pressure. Two alternative design scenarios for a benchmark problem are explored and contrasted in the light of construction emission rates using the EFFC-DFI methodology. A strain-hardening effective stress plasticity model is adopted to simulate the dynamic loading. A combination of input motions, contact pressure, initial vertical total pressure and distance to foundation centreline are employed, as model variables, to further investigate the control of permanent and variable actions on the residual pore pressure ratio. The model is verified against the Ghosh and Madabhushi high acceleration field test database. The outputs of this work is aimed to improve the current computer-aided seismic foundation design that relies on ground’s packing state and consistency. The results confirm that on seismic excitation of shallow foundations, the likelihood of effective stress loss is greater in deeper depths and across free field. For the benchmark problem, adopting a shallow foundation system instead of piled foundation benefitted in a 75% less emission rate, a marked proportion of which is owed to reduced materials and haulage carbon cost.

Application of the latency insertion method (LIM) to the analysis of large circuit interconnect networks

The focus of this research is to explore the applications of the finite difference formulation based on the latency insertion method (LIM) to the analysis of circuit interconnects. Special attention is devoted to addressing the issues that arise in very large networks such as on-chip signal and power distribution networks. We demonstrate that the LIM has the power and flexibility to handle various types of analysis required at different stages of circuit design. The LIM is particularly suitable for simulations of very large scale linear networks and can significantly outperform conventional circuit solvers (such as SPICE).

Physical Design Methodologies for Low Power and Reliable 3D ICs

As the semiconductor industry struggles to maintain its momentum down the path following the Moore's Law, three dimensional integrated circuit (3D IC) technology has emerged as a promising solution to achieve higher integration density, better performance, and lower power consumption. However, despite its significant improvement in electrical performance, 3D IC presents several serious physical design challenges. In this dissertation, we investigate physical design methodologies for 3D ICs with primary focus on two areas: low power 3D clock tree design, and reliability degradation modeling and management. Clock trees are essential parts for digital system which dissipate a large amount of power due to high capacitive loads. The majority of existing 3D clock tree designs focus on minimizing the total wire length, which produces sub-optimal results for power optimization. In this dissertation, we formulate a 3D clock tree design flow which directly optimizes for clock power. Besides, we also investigate the design methodology for clock gating a 3D clock tree, which uses shutdown gates to selectively turn off unnecessary clock activities. Different from the common assumption in 2D ICs that shutdown gates are cheap thus can be applied at every clock node, shutdown gates in 3D ICs introduce additional control TSVs, which compete with clock TSVs for placement resources. We explore the design methodologies to produce the optimal allocation and placement for clock and control TSVs so that the clock power is minimized. We show that the proposed synthesis flow saves significant clock power while accounting for available TSV placement area. Vertical integration also brings new reliability challenges including TSV's electromigration (EM) and several other reliability loss mechanisms caused by TSV-induced stress. These reliability loss models involve complex inter-dependencies between electrical and thermal conditions, which have not been investigated in the past. In this dissertation we set up an electrical/thermal/reliability co-simulation framework to capture the transient of reliability loss in 3D ICs. We further derive and validate an analytical reliability objective function that can be integrated into the 3D placement design flow. The reliability aware placement scheme enables co-design and co-optimization of both the electrical and reliability property, thus improves both the circuit's performance and its lifetime. Our electrical/reliability co-design scheme avoids unnecessary design cycles or application of ad-hoc fixes that lead to sub-optimal performance. Vertical integration also enables stacking DRAM on top of CPU, providing high bandwidth and short latency. However, non-uniform voltage fluctuation and local thermal hotspot in CPU layers are coupled into DRAM layers, causing a non-uniform bit-cell leakage (thereby bit flip) distribution. We propose a performance-power-resilience simulation framework to capture DRAM soft error in 3D multi-core CPU systems. In addition, a dynamic resilience management (DRM) scheme is investigated, which adaptively tunes CPU's operating points to adjust DRAM's voltage noise and thermal condition during runtime. The DRM uses dynamic frequency scaling to achieve a resilience borrow-in strategy, which effectively enhances DRAM's resilience without sacrificing performance. The proposed physical design methodologies should act as important building blocks for 3D ICs and push 3D ICs toward mainstream acceptance in the near future.

POSITIVE FILTERED P_N METHOD FOR LINEAR TRANSPORT EQUATIONS AND THE ASSOCIATED OPTIMIZATION ALGORITHM

We propose a positive, accurate moment closure for linear kinetic transport equations based on a filtered spherical harmonic (FP_N) expansion in the angular variable. The FP_N moment equations are accurate approximations to linear kinetic equations, but they are known to suffer from the occurrence of unphysical, negative particle concentrations. The new positive filtered P_N (FP_N+) closure is developed to address this issue. The FP_N+ closure approximates the kinetic distribution by a spherical harmonic expansion that is non-negative on a finite, predetermined set of quadrature points. With an appropriate numerical PDE solver, the FP_N+ closure generates particle concentrations that are guaranteed to be non-negative. Under an additional, mild regularity assumption, we prove that as the moment order tends to infinity, the FP_N+ approximation converges, in the L2 sense, at the same rate as the FP_N approximation; numerical tests suggest that this assumption may not be necessary. By numerical experiments on the challenging line source benchmark problem, we confirm that the FP_N+ method indeed produces accurate and non-negative solutions. To apply the FP_N+ closure on problems at large temporal-spatial scales, we develop a positive asymptotic preserving (AP) numerical PDE solver. We prove that the propose AP scheme maintains stability and accuracy with standard mesh sizes at large temporal-spatial scales, while, for generic numerical schemes, excessive refinements on temporal-spatial meshes are required. We also show that the proposed scheme preserves positivity of the particle concentration, under some time step restriction. Numerical results confirm that the proposed AP scheme is capable for solving linear transport equations at large temporal-spatial scales, for which a generic scheme could fail. Constrained optimization problems are involved in the formulation of the FP_N+ closure to enforce non-negativity of the FP_N+ approximation on the set of quadrature points. These optimization problems can be written as strictly convex quadratic programs (CQPs) with a large number of inequality constraints. To efficiently solve the CQPs, we propose a constraint-reduced variant of a Mehrotra-predictor-corrector algorithm, with a novel constraint selection rule. We prove that, under appropriate assumptions, the proposed optimization algorithm converges globally to the solution at a locally q-quadratic rate. We test the algorithm on randomly generated problems, and the numerical results indicate that the combination of the proposed algorithm and the constraint selection rule outperforms other compared constraint-reduced algorithms, especially for problems with many more inequality constraints than variables.

Multiple solutions and numerical analysis to the dynamic and stationary models coupling a delayed energy balance model involving latent heat and discontinuous albedo with a deep ocean

We study a climatologically important interaction of two of the main components of the geophysical system by adding an energy balance model for the averaged atmospheric temperature as dynamic boundary condition to a diagnostic ocean model having an additional spatial dimension. In this work, we give deeper insight than previous papers in the literature, mainly with respect to the 1990 pioneering model by Watts and Morantine. We are taking into consideration the latent heat for the two phase ocean as well as a possible delayed term. Non-uniqueness for the initial boundary value problem, uniqueness under a non-degeneracy condition and the existence of multiple stationary solutions are proved here. These multiplicity results suggest that an S-shaped bifurcation diagram should be expected to occur in this class of models generalizing previous energy balance models. The numerical method applied to the model is based on a finite volume scheme with nonlinear weighted essentially non-oscillatory reconstruction and Runge–Kutta total variation diminishing for time integration.

Theoretical investigation of intra- and inter-cellular spatiotemporal calcium patterns in microcirculation

Microcirculatory vessels are lined by endothelial cells (ECs) which are surrounded by a single or multiple layer of smooth muscle cells (SMCs). Spontaneous and agonist induced spatiotemporal calcium (Ca2+) events are generated in ECs and SMCs, and regulated by complex bi-directional signaling between the two layers which ultimately determines the vessel tone. The contractile state of microcirculatory vessels is an important factor in the determination of vascular resistance, blood flow and blood pressure. This dissertation presents theoretical insights into some of the important and currently unresolved phenomena in microvascular tone regulation. Compartmental and continuum models of isolated EC and SMC, coupled EC-SMC and a multi-cellular vessel segment with deterministic and stochastic descriptions of the cellular components were developed, and the intra- and inter-cellular spatiotemporal Ca2+ mobilization was examined.^ Coupled EC-SMC model simulations captured the experimentally observed localized subcellular EC Ca2+ events arising from the opening of EC transient receptor vanilloid 4 (TRPV4) channels and inositol triphosphate receptors (IP3Rs). These localized EC Ca2+ events result in endothelium-derived hyperpolarization (EDH) and Nitric Oxide (NO) production which transmit to the adjacent SMCs to ultimately result in vasodilation. The model examined the effect of heterogeneous distribution of cellular components and channel gating kinetics in determination of the amplitude and spread of the Ca2+ events. The simulations suggested the necessity of co-localization of certain cellular components for modulation of EDH and NO responses. Isolated EC and SMC models captured intracellular Ca2+ wave like activity and predicted the necessity of non-uniform distribution of cellular components for the generation of Ca2+ waves. The simulations also suggested the role of membrane potential dynamics in regulating Ca2+ wave velocity. The multi-cellular vessel segment model examined the underlying mechanisms for the intercellular synchronization of spontaneous oscillatory Ca2+ waves in individual SMC. ^ From local subcellular events to integrated macro-scale behavior at the vessel level, the developed multi-scale models captured basic features of vascular Ca2+ signaling and provide insights for their physiological relevance. The models provide a theoretical framework for assisting investigations on the regulation of vascular tone in health and disease.^

numerical problems in the calculation of electric field and charge density profiles for mass impregnated non-draining hvdc cables subjected to fast polarity reversals

The goal of this simulation thesis is to present a tool for studying and eliminating various numerical problems observed while analyzing the behavior of the MIND cable during fast voltage polarity reversal. The tool is built on the MATLAB environment, where several simulations were run to achieve oscillation-free results. This thesis will add to earlier research on HVDC cables subjected to polarity reversals. Initially, the code does numerical simulations to analyze the electric field and charge density behavior of a MIND cable for certain scenarios such as before, during, and after polarity reversal. However, the primary goal is to reduce numerical oscillations from the charge density profile. The generated code is notable for its usage of the Arithmetic Mean Approach and the Non-Uniform Field Approach for filtering and minimizing oscillations even under time and temperature variations.

Models for reliability estimation of HVDC cable systems

The work carried out in this thesis aims at: - studying – in both simulative and experimental methods – the effect of electrical transients (i.e., Voltage Polarity Reversals VPRs, Temporary OverVoltages TOVs, and Superimposed Switching Impulses SSIs) on the aging phenomena in HVDC extruded cable insulations. Dielectric spectroscopy, conductivity measurements, Fourier Transform Infra-Red FTIR spectroscopy, and space charge measurements show variation in the insulating properties of the aged Cross-Linked Polyethylene XLPE specimens compared to non-aged ones. Scission in XLPE bonds and formation of aging chemical bonds is also noticed in aged insulations due to possible oxidation reactions. The aged materials show more ability to accumulate space charges compared to non-aged ones. An increase in both DC electrical conductivity and imaginary permittivity has been also noticed. - The development of life-based geometric design of HVDC cables in a detailed parametric analysis of all parameters that affect the design. Furthermore, the effect of both electrical and thermal transients on the design is also investigated. - The intrinsic thermal instability in HVDC cables and the effect of insulation characteristics on the thermal stability using a temperature and field iterative loop (using numerical methods – Finite Difference Method FDM). The dielectric loss coefficient is also calculated for DC cables and found to be less than that in AC cables. This emphasizes that the intrinsic thermal instability is critical in HVDC cables. - Fitting electrical conductivity models to the experimental measurements using both models found in the literature and modified models to find the best fit by considering the synergistic effect between field and temperature coefficients of electrical conductivity.