SEARCH · Search NASA

Results for “THERMAL SIMULATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Open-source simulation platform for air source heat pump integrated with thermal energy storage

Here, this article introduces a modular simulation platform for assessing thermal energy storage (TES) integrated with air source heat pumps (ASHP). The Python platform is an open-source library that includes classes for modeling air-air and air-water heat pumps, TES devices, and the heating load of residential buildings. To validate the ASHP model, the study utilized experimental data obtained from a commercial heat pump evaluated at the National Renewable Energy Laboratory (NREL). The results indicate a mean deviation of 0.7% for COP across the operating range, with a maximum relative deviation of 12.6%. In terms of system heating capacity, the model had an average deviation of 4.3% compared to experimental results, with a maximum deviation of 8.2%. Three classes were implemented for modeling distinct types of TES devices: a generic TES based on energy balances, a sensible isothermal water tank, and a stratified water tank. Details of the mathematical models are provided, along with their respective strengths and limitations. An example is provided showcasing the integration of a residential 10 kWh thermal storage unit with an ASHP operating in Denver, CO. The comparison of two different discharge criteria for the TES unit highlights the importance of control strategies in the system performance.

25 ENERGY STORAGE↗

DOE OSTI2026

Verification of global gyrokinetic simulation of low frequency mode excited by thermal plasma in spherical tokamak

A low frequency Alfven eigenmode in the ST40 spherical tokamak (ST) is found to be resonantly excited by thermal ions from global gyrokinetic GTC simulations. The analyses of phase-space resonances and wave-particle energy exchanges show that the Beta-induced Alfvén-Acoustic Eigenmode (BAAE) is primarily driven unstable by thermal ions, but is stabilized by energetic particles (EPs). When all kinetic effects are artificially suppressed, GTC simulations in the ideal magnetohydrodynamic (MHD) limit confirm the existence of the MHD eigenmode with a frequency in the BAAEs gap, as previously found by the ideal MHD code NOVA. The kinetic effects of thermal ions induce a frequency downshift compared to the stable BAAE in the MHD simulations. These results support the verification for global gyrokinetic simulations of the BAAE in the ST. However, the simulations do not explain the experimental observations of the BAAE in the ST40, possibly owing to the limitation of the collisionless gyrokinetic Monte Carlo model for the EPs.

Alfvén eigenmode↗

DOE OSTI2025

Thermalization and criticality on an analogue–digital quantum simulator

Abstract Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators 1,2 . Unlocking the full potential of such systems towards this goal requires flexible initial state preparation, precise time evolution and extensive probes for final state characterization. Here we present a quantum simulator comprising 69 superconducting qubits that supports both universal quantum gates and high-fidelity analogue evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. This hybrid platform features more versatile measurement capabilities compared with analogue-only simulators, which we leverage here to reveal a coarsening-induced breakdown of Kibble–Zurek scaling predictions 3 in theXYmodel, as well as signatures of the classical Kosterlitz–Thouless phase transition 4 . Moreover, the digital gates enable precise energy control, allowing us to study the effects of the eigenstate thermalization hypothesis 5–7 in targeted parts of the eigenspectrum. We also demonstrate digital preparation of pairwise-entangled dimer states, and image the transport of energy and vorticity during subsequent thermalization in analogue evolution. These results establish the efficacy of superconducting analogue–digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics.

Science & Technology - Other Topics↗

DOE OSTI2026

A time-parallel method for scalable heat transfer simulations of additive manufacturing

Here, a major challenge in simulating the thermal behavior in additive manufacturing processes is the disparate length and time scales between transport phenomena occurring in the melt pool and the component. A common simulation approach relies on spatial decomposition for parallel computing, but due to the nature of heat transfer in AM, where most of the computational expenditure is localized near the melt pool, the computational speedup from spatial parallelization saturates quickly. Therefore, additional parallelism by means of time-domain decomposition is needed to fully take advantage of high-performance computing (HPC) resources. This work introduces a time-parallel method to improve the computational scalability of additive manufacturing simulations on HPC systems, while maintaining high temporal resolution of heat transfer near the melt pool. The method, inspired by the nonlinear paraexp formalism, performs an iterative superposition of nonlinear solutions to the initial value problem, integrating the heat equation across overlapping time-parallel intervals. For a single layer of the NIST AMB2018–01 L7 benchmark problem, the method achieves a 38.51x speedup in wall-clock time with a maximum error in the global temperature solution of 0.99%. This reduces the total solution time from 196.72 min to 5.11 min on 128 nodes of the ORNL Frontier supercomputer. The tradeoff between accuracy and total wall-clock time is investigated and recommendations for time-parallel deployment for AM problems are made.

Additive manufacturing↗

DOE OSTI2025

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗

DOE OSTI2026

Reservoir Characterization and Techno-Economic Analysis of Enhanced Geothermal and Closed-Loop Systems in the Wattenberg Field of the Denver-Julesburg Basin, Colorado

This report provides details of reservoir characterization and techno-economic analysis for enhanced geothermal systems and closed-loop geothermal systems in the Denver-Julesburg Basin Wattenberg area. The analysis contributes to the Geothermal Limitless Approach to Drilling Efficiencies (GLADE) project, a multi-institutional research initiative focused on advancing geothermal energy development by reducing drilling costs and improving the rate of penetration. The GLADE project brings together national laboratories, academic institutions, and industry partners, including the U.S. Department of Energy's Geothermal Technologies Office, Occidental Petroleum, the National Laboratory of the Rockies (NLR), Los Alamos National Laboratory, Colorado School of Mines, Louisiana State University, Texas A&M University, and several drilling technology companies. As part of this effort, NLR developed and applied slender-body theory modeling tools and the GEOPHIRES simulator to support techno-economic analysis. The slender-body theory wellbore simulator evaluated thermal performance in U-loop, Eavor-type multilateral, and enhanced geothermal system well configurations. Results indicated that an optimized single U-loop design delivers an outlet temperature of 168.8 degrees Celsius,stabilizing at a final outlet temperature of 135.3 degrees Celsius. Heat production also improved significantly, with an average thermal output of 16.9 MWth and a final thermal output of 11.0 MWth. These results provided critical baseline data for selecting optimized systems for techno-economic analysis. Building on these outputs, NLR applied the GEOPHIRES techno-economic simulator to estimate the levelized cost of electricity for both the Eavor multilateral and enhanced geothermal systems. The sensitivity analysis showed that drilling costs exert a particularly strong impact on the levelized cost of electricity of multilateral systems. Therefore, advances in drilling technology and improved drilling rates could yield significant gains in the economic performance of these systems. In this context, the GLADE project directly supports the U.S. Department of Energy's mission to reduce the cost of geothermal energy and accelerate the deployment of next-generation geothermal technologies.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Effects of Galactic Irradiation on Thermal and Electronic Transport in Tungsten

The impact of irradiation on the thermal and electronic properties of materials is a persistent puzzle, particularly defect formation at the atomic and nanoscales. This work examines the nanoscale effects of low-energy irradiation on tungsten (W), focusing on defect-induced modifications to thermal and electronic transport. Using the site-projected thermal conductivity (SPTC) method [A. Gautam et al. PSS-RRL, 2400306, 2024], bulk and twin-grain boundary W with vacancy defects based on the Norgett–Robinson–Torrens displacements per atom model is analyzed. SPTC provides a detailed prediction of postcascade spatial thermal conductivity distribution. Electronic conductivity activity is estimated using the “N 2 method” [K. Nepal et al. Carbon, 119711, 2025] to explore the consequences of vacancies and grain boundaries, highlighting the defect-dependent nature of charge transport behavior. These findings offer high-resolution insights into irradiation-driven transport phenomena, with implications for space-exposed materials and nanoscale thermal/electronic management.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

RaDIATE Collaboration Thermal shock studies

As next-generation accelerator target facilities (High Energy Physics, Spallation Sources, ... ) become increasingly more powerful and intense, high power target systems face key technical challenges, such as radiation damage and thermal shock. Those ultimately degrade the performance and lifetime of targets and have been identified as the leading cross-cutting challenges of high-power target facilities. In order to operate reliable beam-intercepting devices in the framework of energy and intensity increase for next generation accelerators, the RaDIATE Collaboration (Radiation Damage In Accelerator Target Environment), established in 2012 and managed by Fermilab, brings together existing expertise in nuclear material and accelerator targets from 20 international institutions, including CERN, to execute a coordinated strategy for high power targetry R&D. In this context, several thermal shock studies were performed at CERN's HiRadMat (High-Radiation to Materials) facility, that took a key step towards improving our knowledge on target damage tolerance. The first experiment HRMT-24, supported by EURCARD2 and completed in 2015, successfully validated the Johnson-Cook strength model developed at SwRI on beryllium S200FH (used for beam window material), providing a better confidence in simulating the thermal shock response of current and future S200FH beryllium components. HRMT-43, supported by ARIES and completed in 2018, tested various materials (Be, C, SiC, Si, Ti and ceramic nanofiber). It was a first and unique test which included pre-irradiated specimens from high energy proton beam irradiation to identify thermal shock response differences between non-irradiated and previously irradiated materials. Real-time measurement of dynamic thermomechanical response of graphite helped to benchmark numerical simulations.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Integrated Steady-State System Package for Nuclear Thermal Propulsion Analysis Using Multi-Dimensional Thermal Hydraulics and Dimensionless Turbopump Treatment

Nuclear thermal propulsion is an evolving technology that can be utilized for long-distance space travel. This technology yields the advantage of a high thrust and specific impulse, but requires an examination of the potential design adjustments necessary to enhance its feasibility. The development of nuclear thermal propulsion requires a comprehensive understanding of the system-level behavior during transient and steady-state operation. This paper extends our previous research by including the proper handling of turbomachinery with multi-channel thermal hydraulic simulations only for steady-state solutions. The system-level approach presented here enables the treatment of the turbopump components through non-dimensional analysis that eliminates the assumption of constant efficiencies. All the other components within the system (e.g., reflector and core) can be discretized to multiple channels and layers, in which the full thermal hydraulic solution is established. The approach chosen here enables the realistic modeling of the propellant flow within the expander cycle by capturing the pressure losses, mass flow rate splits, and enthalpy gain for various operational conditions. The verification of the package is completed through point comparisons of previous investigations into similar system designs. Furthermore, sensitivity studies are used to benchmark the capabilities of the package and investigate solution variations due to the perturbation of operational conditions and regimes. The sensitivity studies performed here are important to capture variation in flow characteristics (e.g., temperature, pressure, mass flow rates) for different design objectives such as the thrust and specific impulse. This work demonstrates that system-level simulations lacking multi-channel capability and proper turbomachinery treatment may yield higher uncertainties in understanding the engine’s response and characteristics to changing various requirements. This is extremely important when screening the design space of such propulsion systems and when transient simulations are required.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Design and Study of Inductively Coupled Plasma Chamber Components Using the SupRISE Test Device at DIII-D

The DIII-D National Fusion Facility aims to increase the auxiliary heating power for the tokamak by upgrading the Neutral Beam Injection (NBI) system. In collaboration with North Carolina State University, the conventional ‘arc and filament’ NBI ion sources will be converted to inductively coupled plasma (ICP) sources which utilize radio frequency (RF) coupling to maximize reliability for high power operation. In support of this initiative, a full scale test device (SupRISE, Superior Radiofrequency Ion Source Experiment) is currently under construction at the DIII-D Facility. In preparation for the construction of a full scale prototype that can be installed on the DIII-D NBI system, experiments on SupRISE will be conducted to determine the optimal RF frequency for high power coupling, Faraday shield slit configuration, and ICP chamber length. SupRISE is comprised of an approximately 30 x 70 cm quartz dielectric vessel with an internal Faraday shield enclosed in a secondary vacuum chamber to ensure structural stability of the dielectric. Actively cooled front and back plates are designed to reduce the thermal stresses on the plasma facing components and mate with the existing accelerator used by the NBI system at DIII-D. 50 kW of RF power will be coupled to the plasma through the quartz over a variable frequency range of 4-8 MHz to sustain a plasma density of ∼ 10 18 m −3 for a 10 s ON, 210 s OFF duty cycle. Various modelling efforts have been employed to simulate the thermal and stress profiles over the primary components of the SupRISE device as well as the inductance behavior of the RF antenna. These simulation results and the final design for SupRISE will be presented. An additional reduced-scale predecessor ICP source (called RISE) has been used to develop a predictive match model that will be applied to frequency optimization studies on SupRISE. Furthermore, the outcomes of this research and complementary efforts at North Carolina State University are essential for the incorporation of ICP NBI positive ion sources at the DIII-D facility.

DIII-D↗

DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

The formation of an radial edge electric field due to finite ion orbit width effects is the possible root cause of the H-mode edge

Abstract Full orbit-following simulations of thermal ions show that finite ion-orbit width effects create charge separation near the last closed flux surface (LCFS) which generates a localized radial electric field. Experimentally, edge electric fields are observed in H-mode plasmas and they are necessary for the edge turbulence suppression via the E × B flow shear mechanism. Confined trapped (and to a lesser extent co-passing) ions near the plasma edge form a positive charge distribution outside the LCFS, while thermal electrons are tied more tightly to field lines owing to their small mass and are poorly confined outside the LCFS, hence charge neutrality is violated outside the LCFS. A large number of reported observations from spherical and conventional tokamaks support the results from the simulations although the simulations were not performed fully self consistently. The results suggest ways to lower the H-mode power threshold and optimize the H-mode plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Which way does the dendrite grow? Competition among epitaxy, preferred growth direction, and thermal gradients in powder bed fusion additive manufacturing

The as-processed microstructure of metal alloy parts manufactured through laser powder bed fusion (LPBF) is heavily derived from the cellular dendritic solidification. The growth direction of dendrites within the melt pool is determined through competition among epitaxial growth, preferred growth directions, and maximum thermal gradients. However, the dominant factor and the specific role of each in developing melt pool microstructures remain unknown. Here, in this study, we performed single laser track scans on an SS316L single crystal substrate and combined experimental characterization of microstructure and crystal orientations with Computational Fluid Dynamics simulations of thermal gradients to evaluate the role of each factor in determining dendritic growth direction and evolution. Our results reveal that epitaxial growth dominates microstructure development by preferentially growing along a single 〈100〉 variant of the single crystal substrate adjacent to the melt pool boundary. Under LPBF’s highly curved and rapidly evolving thermal field, this preferential dendrite variant selection and its continued growth from the melt pool boundary to the centerline are governed by the local temperature gradient magnitude at the solid-liquid interface, rather than by the instantaneous maximum temperature gradient direction alone. Using these findings, we successfully predict changes in the dendrite growth direction with changing laser scan direction on a single crystal substrate, and show that the geometric melt pool centerline can deviate from the microstructural centerline because asymmetric local temperature gradient magnitudes transiently limit growth, resulting in different dendrite travel distances on each side of the melt pool.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Quantifying the thermal effect and methyl radical production in nanosecond repetitively pulsed glow discharges applied to a methane-air flame

In this work, we investigated non-equilibrium plasma produced by nanosecond repetitively pulsed glow discharges applied across a lean premixed methane-air flame. The flame is stationary, axisymmetric, and laminar. The discharges are applied on the symmetry axis crossing the reactant gases, flame front, and product gases, allowing phase-locked averaged measurements and comparisons with axisymmetric numerical simulations. The thermal effect and methyl radical production are quantified in the discharge in the reactant gas region. One-dimensional, two-beam, hybrid, femtosecond-picosecond, coherent anti-Stokes Raman scattering is used to acquire spatial and temporal profiles of temperature and oxygen-to-nitrogen concentration ratio. Photo-fragmentation laser-induced fluorescence is used to acquire quantitative two-dimensional profiles of methyl radicals in the discharge providing the first quantitative imaging of methyl produced ahead of a flame by plasma-induced methane dissociation. The spatial profiles of temperature and oxygen-to-nitrogen concentration ratio are in steady state, indicating that individual discharges have an insignificant heating effect. Upper and lower bounds of the produced mole fraction of methyl radicals in the plasma are obtained due to uncertainties in the collisional quenching rates of excited state methylidyne radicals in the plasma. The discharges produce a maximum of 600–1100 ppm of methyl radicals upstream of the flame front within 25 ns. This amount is similar to the predicted methyl mole fraction for the flame without plasma and thus represents a significant chemical perturbation to the reactants upstream of the flame front. The produced methyl follows an exponential decay in the first microsecond after the discharge with a decay constant of 8 µs close to the flame, and 0.8 µs further from the flame. The decay then deviates from the exponential curve and the methyl persists for tens of microseconds. The results suggest that for the tested configuration, the thermal effect of individual discharges through fast gas heating is negligible, while active chemical species are produced in large quantities in the reactant gases, upstream of the flame front.

hybrid fs-ps CARS↗

DOE OSTI2026

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

WEC Irradiation Can Thermal Testing Report

With renewed interest for the next generation of small modular reactors and space-based nuclear reactors, metal hydrides have been a focus of research for moderating materials to reduce the need for highly enriched fuel to maintain small core sizes. Yttrium hydride (YH) is an attractive candidate for these metal hydrides due to its ability to retain hydrogen at elevated temperatures as well as the low neutron absorption cross section of yttrium metal. Given the necessity for reactors to operate at elevated temperatures for extended periods of time, understanding hydrogen retention in moderator channels in simulated operating thermal conditions is paramount in predicting reactor performance and the long-term stability of YH. Additionally, coating strategies on the moderator encapsulating will be explored to ascertain if it aids in retention of hydrogen within the capsule and moderator.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Fabrication and Performance Evaluation of Double-Sided Copper Nanowire Arrays as Thermal and Electrical Interfacial Layers

Reducing contact interface thermal and electrical resistances is in great demand across various industries, particularly in the semiconductor industry. This study introduces an approach using double-sided copper nanowire (Cu NW) arrays on copper sheets as both thermal and electrical interfacial layers, designed to effectively accommodate the topographical inconsistencies between contact surfaces. Further, experimental outcomes reveal a significant reduction in thermal contact resistance (TCR), with a value of 2.5 mm 2 K W –1 , thereby exceeding the efficiency of reported nanostructural thermal interface materials (TIMs). Additionally, when utilized as an electrical interfacial layer, these double-sided Cu NWs arrays dramatically reduced electrical contact resistance (ECR), outperforming traditional conductive grease in applications necessitating separable bonding, though showing comparable performance to costly silver-based conductive epoxies required for permanent, inseparable bonds. The promising results of the double-sided Cu NWs arrays in reducing both TCR and ECR, confirmed by finite element simulation, highlight their substantial potential in advancing TIMs and electrical interconnection applications across various sectors.

42 ENGINEERING↗