Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Flux”

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.

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

At least 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

High-fidelity pool boiling simulations on multiple nucleation sites using interface capturing method

Boiling has proved to be one of the most efficient means for heat transfer and is a very important phenomenon during severe accident scenarios in light water reactors. High-fidelity pool boiling simulations can provide a numerical database for improving mechanistic boiling models by allowing for specific evaluation of interactions among bubbles. Previously published pool boiling simulations investigated two nucleation sites in which bubble growth at one site suppressed nucleation at the other site. Based on previous study results, more complicated interface-capturing simulations on pool boiling were conducted using PHASTA code with locally refined unstructured mesh. First, different boundary conditions (BCs) were assessed to support robustness and reproducibility of the boiling model. Then, a scale study was conducted at a larger domain with nine nucleation sites where either nine or four nucleation sites are activated. Involving more nucleation sites increased the complexity of bubble interactions from surrounding sites. Finally, bubble departure behavior influenced by wall heat flux was investigated. In this work, when heat flux was increased, the order of bubble departure changed, but diagonal bubbles always departed after one another. The departure time interval between the first and second bubble reduced as heat flux increased. The corresponding frequency was almost linearly proportional to the heat flux. In addition, bubble departure behavior was found to be greatly influenced by the nucleation site pattern. Multiple nucleation sites resulted in superimposed inhibitive effects from surrounding sites to each bubble, which extensively delayed the departure. This new observation was not discussed in previously published works. The work presented here provides new insight on the fundamental understanding of boiling phenomena, contributes to the development of a 3D multiphase computational fluid dynamics (M-CFD) model, and provides a more comprehensive database for data-driven pool boiling studies.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2012

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗
DOE OSTI2024

ESM data downscaling: a comparison of super-resolution deep learning models

Abstract Climate projections at fine spatial resolutions are required to conduct accurate risk assessment for critical infrastructure and design adaptation planning. Generating these projections using advanced Earth system models (ESM) requires significant computational resources. To address this issue, various statistical downscaling techniques have been introduced to generate fine-resolution data from coarse-resolution simulations. In this study, we evaluate and compare five deep learning-based downscaling techniques, namely, super-resolution convolutional neural networks, fast super-resolution convolutional neural network ESM, efficient sub-pixel convolutional neural network, enhanced deep residual network (EDRN), and super-resolution generative adversarial network (SRGAN). These techniques are applied to a dataset generated by the Energy Exascale Earth System Model (E3SM), focusing on key surface variables such as surface temperature, shortwave heat flux, and longwave heat flux. Models are trained and validated using paired fine-resolution (0.25 $$^{\circ }$$ ∘ ) and coarse-resolution (1 $$^{\circ }$$ ∘ ) monthly data obtained from a 9-year simulation. Next, blind testing is performed using monthly data obtained from two different years outside of the training and validation set. To evaluate the efficiency of each technique, different statistical metrics are used, including mean squared error (MSE), peak signal-to-noise ratio (PSNR), structural similarity index measure (SSIM), and learned perceptual image patch similarity (LPIPS). The results show that EDRN outperforms other algorithms in terms of PSNR, SSIM, and MSE, but struggles to capture fine-scale features in the data. In contrast, SRGAN, a generative model that uses perceptual loss, excels in capturing fine details at boundaries and internal structures, resulting in lower LPIPS than other methods.

Pawar, Nikhil M. (ORCID:0000000211613289)↗
NASA NTRS1972

Thermal elastoplastic structural analysis of non-metallic thermal protection systems

An incremental theory and numerical procedure to analyze a three-dimensional thermoelastoplastic structure subjected to high temperature, surface heat flux, and volume heat supply as well as mechanical loadings are presented. Heat conduction equations and equilibrium equations are derived by assuming a specific form of incremental free energy, entropy, stresses and heat flux together with the first and second laws of thermodynamics, von Mises yield criteria and Prandtl-Reuss flow rule. The finite element discretization using the linear isotropic three-dimensional element for the space domain and a difference operator corresponding to a linear variation of temperature within a small time increment for the time domain lead to systematic solutions of temperature distribution and displacement and stress fields. Various boundary conditions such as insulated surfaces and convection through uninsulated surface can be easily treated. To demonstrate effectiveness of the present formulation a number of example problems are presented.

Chung, T. J.↗
NASA NTRS1985

Methods for heat transfer and temperature field analysis of the insulated diesel phase 2 progress report

This report describes work done during Phase 2 of a 3 year program aimed at developing a comprehensive heat transfer and thermal analysis methodology for design analysis of insulated diesel engines. The overall program addresses all the key heat transfer issues: (1) spatially and time-resolved convective and radiative in-cylinder heat transfer, (2) steady-state conduction in the overall structure, and (3) cyclical and load/speed temperature transients in the engine structure. During Phase 2, radiation heat transfer model was developed, which accounts for soot formation and burn up. A methodology was developed for carrying out the multi-dimensional finite-element heat conduction calculations within the framework of thermodynamic cycle codes. Studies were carried out using the integrated methodology to address key issues in low heat rejection engines. A wide ranging design analysis matrix was covered, including a variety of insulation strategies, recovery devices and base engine configurations. A single cylinder Cummins engine was installed at Purdue University, and it was brought to a full operational status. The development of instrumentation was continued, concentrating on radiation heat flux detector, total heat flux probe, and accurate pressure-crank angle data acquisition.

Morel, T.↗
NASA NTRS2001

Thermal Protection Materials for Reentry and Planetary Applications

Thermal protection materials and systems (TPS) are used to protect spacecraft during reentry into Earth's atmosphere or entry into planetary atmospheres. As such, these materials are subject to severe environments with high heat fluxes and rapid heating. Catalytic effects can increase the temperatures substantially. These materials are also subject to impact damage from micrometeorites or other debris during ascent, orbit, and descent, and thus must be able to withstand damage and function following damage. Thermal protection materials and coatings used in reusable launch vehicles will be reviewed, including the needs and directions for new materials to enable new missions that require faster turnaround and much greater reusability. The role of ablative materials for use in high heat flux environments, especially for non-reusable applications and upcoming planetary missions, will be discussed. New thermal protection system materials may enable the use of sharp nose caps and leading edges on future reusable space transportation vehicles. Vehicles employing this new technology would have significant increases in maneuverability and out-of-orbit cross range compared to current vehicles, leading to increased mission safety in the event of the need to abort during ascent or from orbit. Ultrahigh temperature ceramics, a family of materials based on HfB2and ZrB2 with SiC, will be discussed. The development, mechanical and thermal properties, and uses of these materials will be reviewed.

Johnson, Sylvia M.↗
NASA NTRS2001

Thermal Protection Materials for Reentry Applications

Thermal protection materials and systems (IRS) are used to protect spacecraft during reentry into Earth's atmosphere or entry into planetary atmospheres. As such, these materials are subject to severe environments with high heat fluxes and rapid heating. Catalytic effects can increase the temperatures substantially. These materials are also subject to impact damage from micrometeorites or other debris during ascent, orbit, and descent, and thus must be able to withstand damage and to function following damage. Thermal protection materials and coatings used in reusable launch vehicles will be reviewed, including the needs and directions for new materials to enable new missions that require faster turnaround and much greater reusability. The role of ablative materials for use in high heat flux environments, especially for non-reusable applications and upcoming planetary missions, will be discussed. New thermal protection system materials may enable the use of sharp nose caps and leading edges on future reusable space transportation vehicles. Vehicles employing this new technology would have significant increases in maneuverability and out-of-orbit cross range compared to current vehicles, leading to increased mission safety in the event of the need to abort during ascent or from orbit. Ultrahigh temperature ceramics, a family of materials based on HfB2 and ZrB2 with SiC, will be discussed. The development, mechanical and thermal properties, and uses of these materials will be reviewed.

Johnson, Sylvia M.↗
NASA NTRS2017

Validation of a Thermo-Ablative Model of Elastomeric Internal Insulation Materials

In thermo-ablative material modeling, as in many fields of analysis, the quality of the existing models significantly exceeds that of the experimental data required for their validation. In an effort to narrow this gap, a laboratory-scale internal insulation test bed was developed that exposes insulation samples to realistic solid rocket motor (SRM) internal environments while being instrumented to record real-time rates of both model inputs (i.e., chamber pressure, total surface heat flux, and radiative heat flux) as well as model outputs (i.e., material decomposition depths (MDDs) and in-depth material temperatures). In this work, the measured SRM internal environment parameters were used in conjunction with equilibrium thermochemistry codes as inputs to one-dimensional thermo-ablative models of the PBINBR and CFEPDM insulation samples used in the lab-scale test firings. The computed MDD histories were then compared with those deduced from real-time X-ray radiography of the insulation samples, and the calculated in-depth temperatures were compared with those measured by embedded thermocouples. The results of this exercise emphasize the challenges of modeling and testing elastomeric materials in SRM environments while illuminating the path forward to improved fidelity.

Martin, Heath T.↗
NASA NTRS2022

Experimental Heat Transfer Results and Flow Visualization of Vertical Upflow Boiling in Earth Gravity with Subcooled Inlet Conditions – In Preparation for Experiments Onboard the International Space Station

Since 2012, researchers at the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and NASA Glenn Research Center have been collaborating on a long-term effort to study flow boiling and condensation in microgravity. The ultimate goal has been to develop the Flow Boiling and Condensation Experiment (FBCE) for the International Space Station (ISS). Based on the findings from prior flow boiling experiments both at different orientations in Earth gravity and onboard parabolic flights simulating short durations of microgravity, a final refined experiment design, construction, and operating procedure have been arrived at for long-duration microgravity flow boiling experiments onboard the ISS. This study investigates flow boiling of n-Perfluorohexane with subcooled inlet in a rectangular channel of dimensions 114.6 mm heated length, 2.5 mm width, and 5 mm height. These pre-launch experiments (Mission Sequence Testing) were conducted in vertical upflow orientation in Earth gravity using the same experimental rig that was launched to the ISS in August 2021. The various operating parameters varied are heating configuration (single- and double-sided), mass velocity (180 – 3200 kg/m2s), inlet subcooling (+0 – 32°C, encompassing both highly subcooled and near-saturated inlet conditions), and inlet pressure (119 – 191 kPa). High-speed video flow visualization images are presented to explain the two-phase interfacial physics within the channel’s heated section. Heat transfer results in terms of flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and averaged heat transfer coefficients are analyzed and parametric effects elucidated. Severe temporal thermodynamic equilibrium is observed for near-saturated inlet at very low velocities. Nucleate boiling degradation starts at larger heat fluxes for single-sided heating than double sided at low mass velocities with highly subcooled inlet, and conversely at high mass velocities with near-saturated inlet. Nucleate boiling degradation can be delayed to higher heat fluxes by highly subcooling the inlet and increasing mass velocity. The entire local heat transfer coefficient profiles are degraded at higher heat fluxes for near-saturated inlet, but only the downstream part for highly subcooled inlet. This study also confirmed reliability of the upcoming ISS experimental data for subcooled inlet conditions and the collected Earth-gravity data will be used for comparison against the ISS data.

flow boiling↗
DOE OSTI2021

Experimental Results of 2-m Heptane, Bakken Crude Oil, and Dilbit Crude Oil Pool Fire Tests Performed for the National Research Council of Canada

This report provides results from a series of 2-m pool fire experiments performed in the Thermal Test Complex at Sandia National Laboratories testing heptane, Bakken crude oil, and dilbit crude oil. The effect of the presence and placement of a calorimeter, fuel supply temperature, and maintaining a constant fuel level were assessed. Measurements include burn rate, surface emissive power, flame height, heat flux to an engulfed calorimeter, heat flux to external instruments, thermocouple temperatures within the fuel and fire plume, and heat release rate. The results indicate that the presence and placement of the calorimeter has the most effect on the measured quantities for the Bakken crude oil and indicated no effect for the Dilbit crude oil. The fuel feed temperature had a slight effect for the heptane fuel, but not for the crude oils. Allowing the fuel to burn down did not have a significant effect on any of the fuels. The Bakken crude oil resulted in the highest average total heat flux to the calorimeter by a factor of about 1.5 and 1.3 higher compared to heptane and the dilbit crude oil, respectively.

02 PETROLEUM↗
NASA NTRS2010

Effects of Solar Photovoltaic Panels on Roof Heat Transfer

Building Heating, Ventilation and Air Conditioning (HVAC) is a major contributor to urban energy use. In single story buildings with large surface area such as warehouses most of the heat enters through the roof. A rooftop modification that has not been examined experimentally is solar photovoltaic (PV) arrays. In California alone, several GW in residential and commercial rooftop PV are approved or in the planning stages. With the PV solar conversion efficiency ranging from 5-20% and a typical installed PV solar reflectance of 16-27%, 53-79% of the solar energy heats the panel. Most of this heat is then either transferred to the atmosphere or the building underneath. Consequently solar PV has indirect effects on roof heat transfer. The effect of rooftop PV systems on the building roof and indoor energy balance as well as their economic impacts on building HVAC costs have not been investigated. Roof calculator models currently do not account for rooftop modifications such as PV arrays. In this study, we report extensive measurements of a building containing a flush mount and a tilted solar PV array as well as exposed reference roof. Exterior air and surface temperature, wind speed, and solar radiation were measured and thermal infrared (TIR) images of the interior ceiling were taken. We found that in daytime the ceiling surface temperature under the PV arrays was significantly cooler than under the exposed roof. The maximum difference of 2.5 C was observed at around 1800h, close to typical time of peak energy demand. Conversely at night, the ceiling temperature under the PV arrays was warmer, especially for the array mounted flat onto the roof. A one dimensional conductive heat flux model was used to calculate the temperature profile through the roof. The heat flux into the bottom layer was used as an estimate of the heat flux into the building. The mean daytime heat flux (1200-2000 PST) under the exposed roof in the model was 14.0 Watts per square meter larger than under the tilted PV array. The maximum downward heat flux was 18.7 Watts per square meters for the exposed roof and 7.0 Watts per square meters under the tilted PV array, a 63% reduction due to the PV array. This study is unique as the impact of tilted and flush PV arrays could be compared against a typical exposed roof at the same roof for a commercial uninhabited building with exposed ceiling and consisting only of the building envelope. Our results indicate a more comfortable indoor environment in PV covered buildings without HVAC both in hotter and cooler seasons.

Dominguez, A.↗
NASA NTRS2023

Assessing the Performance of the South American Land Data Assimilation System Version 2 (SALDAS-2) Energy Balance across Diverse Biomes

Understanding the exchange of energy between the surface and the atmosphere is important in view of the climate scenario. However, it becomes a challenging task due to a sparse network of observations. This study aims to improve the energy balance estimates for the Amazon, Cerrado, and Pampa biomes located in South America using the radiation and precipitation forcing obtained from the Clouds and the Earth’s Radiant Energy System (CERES) and the precipitation CPTEC/MERGE datasets. We employed three surface models—Noah-MP, Community Land Model (CLSM), and Integrated Biosphere Simulator (IBIS)—and conducted modeling experiments, termed South America Land Data Assimilation System (SALDAS-2). The results showed that SALDAS-2 radiation estimates had the smallest errors. Moreover, SALDAS-2 precipitation estimates were better than the Global Land Data Assimilation System (GLDAS) in the Cerrado (MBE = −0.16) and Pampa (MBE = −0.19). Noah-MP presented improvements compared with CLSM and IBIS in 100% of towers located in the Amazon. CLSM tends to overestimate the latent heat flux and underestimate the sensible heat flux in the Amazon. Noah-MP and Ensemble outperformed GLDAS in terms latent and sensible heat fluxes. The potential of SALDAS-2 should be emphasized to provide more accurate estimates of surface energy balance.

modeling↗
NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗
DOE OSTI2020

Nonlinear δf particle in cell gyrokinetic simulations of the microtearing mode

Nonlinear simulations are carried out for the microtearing mode using particle-based δf gyrokinetic simulations for parameters relevant to spherical tokamaks. The present study finds that the microtearing mode can generate significant electron heat flux, which is predominantly carried out by the electromagnetic component of the heat flux with a negligible contribution from the electrostatic component. The mode sustains without the electrostatic component. We observe that the electron heat flux increases with the electron temperature gradient. The heat flux exhibits a very weak dependence on the collisions. It increases with electron β initially; however, at very high β, the electron heat flux is reduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Near-Field Radiative Heat Transfer Between Two α-MoO3 Biaxial Crystals

The near-field radiative heat transfer (NFRHT) between two semi-infinite α-MoO3 biaxial crystals is investigated numerically based on the fluctuation–dissipation theorem combined with the modified 4 × 4 transfer matrix method in this paper. In the calculations, the near-field radiative heat flux (NFRHF) along each of the crystalline directions of α-MoO3 is obtained by controlling the orientation of the biaxial crystals. The results show that much larger heat flux than that between two semi-infinite hexagonal boron nitride can be achieved in the near-field regime, and the maximum heat flux is along the [001] crystalline direction. The mechanisms for the large radiative heat flux are explained as due to existence of hyperbolic phonon polaritons (HPPs) inside α-MoO3 and excitation of hyperbolic surface phonon polaritons (HSPhPs) at the vacuum/α-MoO3 interfaces. The effect of relative rotation between the emitter and the receiver on the heat flux is also investigated. It is found that the heat flux varies significantly with the relative rotation angle. The modulation contrast can be as large as two when the heat flux is along the [010] direction. We attribute the large modulation contrast mainly to the misalignment of HSPhPs and HPPs between the emitter and the receiver. Hence, the results obtained in this work may provide a promising way for manipulating near-field radiative heat transfer between anisotropic materials.

Engineering↗
DOE OSTI2023

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗
NASA NTRS1998

Decadal Air-Sea Interaction in the North Atlantic Based on Observations and Modeling Results

The decadal, 12-14 year, cycle observed in the North Atlantic SST and tide gauge data was examined using the NCEP/NCAR reanalyses, COADS data and an ocean model simulation. Besides this decadal mode, a shorter, subdecadal period of about 8 years exists in tide gauge data north of 40N, in the subpolar SST and in the winter North Atlantic Oscillation (NAO) index and in subpolar winter heat flux values. The decadal cycle is a well separated mode in a singular spectrum analysis (SSA) for a time series of SST EOF mode 1 with a center over the Gulf Stream extension. Tide gauge and SST data are consistent in that both show a significant subdecadal periodicity exclusively in the subpolar gyre, but in subtropics the 12-14 year period is the prominent, but nonstationary, decadal signal. The main finding of this study is that this 12-14 year cycle can be constructed based on the leading mode of the surface heat flux. This connection to the surface heat flux implicates the participation of the thermohaline circulation in the decadal cycle. During the cycle starting from the positive index phase of NAO, SST and oceanic heat content anomalies are created in subtropics due to local heat flux and intensification of the thermohaline circulation. The anomalies advect to the subpolar gyre where they are amplified by local heat flux and are part of the negative feedback of thermohaline circulation on itself. Consequently the oceanic thermohaline circulation slows down and the opposite cycle starts. The oscillatory nature would not be possible without the active atmospheric participation in the cycle, because it provides the unstable interaction through heat flux, without it, the oceanic mode would be damped. This analysis suggests that the two principal modes of heat flux variability, corresponding to patterns similar to North Atlantic Oscillation (NAO) and Western Atlantic (WA), are part of the same decadal cycle and an indirect measure of the north-south movement of the storm tracks.

Hakkinen, Sirpa↗
DOE OSTI2020

Budgets for Decadal Variability in Pacific Ocean Heat Content

A slowdown in the rate of surface warming in the early 2000s led to renewed interest in the redistribution of ocean heat content (OHC) and its relationship with internal climate variability. In this work we use the Community Earth System Model version 1 to study the relationship between OHC and the interdecadal Pacific oscillation (IPO), a major mode of decadal sea surface temperature variability in the Pacific Ocean. By comparing the relative contributions of surface heat flux and ocean dynamics to changes in OHC for different phases of the IPO, we try to identify the underlying physical processes involved. Our results suggest that during IPO phase transitions, changes of 0–300-m OHC across the northern extratropical Pacific are positively contributed by both surface heat flux and oceanic heat transport. By contrast, oceanic heat transport appears to drive the OHC changes in equatorial Pacific whereas surface heat flux acts as a damping term. During a positive IPO phase, weakened wind-driven circulation acts to increase the OHC in the equatorial Pacific while the enhanced evaporation acts to damp OHC anomalies. In the Kuroshio–Oyashio Extension region, a dipole anomaly of zonal heat advection amplifies an OHC dipole anomaly that moves eastward, while strong turbulent heat fluxes act to dampen this OHC anomaly. In the northern subtropical Pacific, both the wind-driven evaporation change and the change of zonal heat advection along Kuroshio Extension contribute to the OHC change during phase transition. For the northern subpolar Pacific, both surface heat flux and enhanced meridional advection contribute to the positive OHC anomalies during the positive IPO phase.

54 ENVIRONMENTAL SCIENCES↗