Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermal boundary conductance”

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 145 records · Page 8

Selective laser sintering and spark plasma sintering of (Zr,Nb,Ta,Ti,W)C compositionally complex carbide ceramics

Abstract Two advanced manufacturing processes, spark plasma sintering (SPS) and selective laser sintering (SLS), have been developed for synthesis of (Zr,Nb,Ta,Ti,W)C compositionally complex carbide (CCC) via reactive sintering of a powder mixture of constitute monocarbides. X‐ray diffraction analysis confirmed that the single‐phase CCC can be formed by both SPS and SLS. While a homogenous microstructure with uniform metal element distributions was developed during SPS, three‐layer microstructures with a thin TiC‐rich layer and two TaC‐rich layers along with a TiO 2 ‐rich surface layer containing W nanoparticles were formed during SLS. In addition, cellular structures with W, Zr, and Ti element segregation and dislocations on cell boundaries were observed in the SLS‐CCC sample, indicating the effect of nonequilibrium conditions on microstructure formation during laser melting followed by rapid cooling and solidification process. Compared to the SPS‐CCC sample, the SLS‐CCC showed enhanced hardness and reduced thermal conductivity, which may be related to their unique cellular structures.

36 MATERIALS SCIENCE↗

Numerical investigation of thermomechanical response of multiscale porous Ultra-High Temperature Ceramics

Recent advances in Ultra-High Temperature Ceramics (UHTC) manufacturing have permitted the development of multiscale porous UHTC microstructures. Within the target application of hypersonic vehicles, dense UHTCs are suitable for thermal protection on leading edges, whereas porous UHTCs may find a role in providing thermally insulated interfaces for temperature-sensitive interior components. Designing vehicles incorporating porous UHTCs requires a characterization of their thermomechanical properties across the full range of expected operating temperatures spanning –20 °C–2500 °C. This research represents a preliminary study in performing this characterization. Several numerical experiments are performed using a coupled thermomechanical implementation of the Material Point Method to determine the temperature dependence of effective material properties both with and without damage. Furthermore, complex time-dependent boundary conditions derived from known hypersonic flight profiles are simulated in order to probe the various couplings between deformation, damage and heat transfer. It is shown that the model reveals the importance of micro-buckling in determining effective material stiffness and thermal conductivity.

36 MATERIALS SCIENCE↗

Development of process parameters and post-build conditions for qualification of LPBF 316 SS

To harness the potential of laser powder bed fusion (LPBF) 316H stainless steel (SS) for use in advanced nuclear reactors, extensive research efforts are needed to develop optimal laser processing parameters and appropriate heat treatments, taking into account various manufacturing platforms and locations. Furthermore, it is crucial to assess the material properties in environments relevant to reactor operating conditions. In collaboration with Oak Ridge National Laboratory (ORNL) and Los Alamos National Laboratory (LANL), Argonne National Laboratory (ANL) is dedicated to gaining insights into how the manufacturing process and post-build treatments impact the performance of LPBF 316H SS. This report provides an overview of ANL's research findings for FY23, focusing on LPBF 316H SS produced using a Renishaw AM400 laser system. Specifically, our research has concentrated on three key areas: 1) Optimum laser processing parameters: Leveraging the high-throughput printing technology, we systematically varied the laser power, exposure time and point distance within the same build. By measuring the as-printed porosity, we identified an optimum printing parameter window that produced materials with near-complete density; 2) Post-built treatment development: Materials printed with the optimum laser parameters underwent treatments including stress relief, solution annealing, and hot isostatic pressing. Subsequently, we fabricated tension, creep, and fatigue specimens from materials subjected to these different conditions. Over the next year, a series of scoping tests will be conducted to facilitate the selection of the most suitable post-build treatment condition; 3) Thermal aging effects on microstructure and properties: thermal aging at 550°C, 650°C and 750°C was conducted on LPBF 316H SS specimens up to 2500 h. Microstructural characterization was performed with electron microscopy, and the evolution of the dislocation cell structures and secondary phases was studied. Microhardness tests and tension tests were conducted to quantify changes in mechanical properties. The results revealed that the aged LPBF 316H SS exhibited a high density of Cr-rich and Mo-rich fine precipitates within grains due to the presence of a high density of dislocations and the globally distributed dislocation cell structures, which served as heterogeneous nucleation sites for secondary phases. Along grain boundaries, a Mo-Cr-rich phase, observed in 750°C-aged conditions, displayed significant thermal coarsening. Throughout the aging process, various mechanisms, including stress relaxation, dislocation cell structure recovery, solution hardening, and precipitation hardening, competed to influence materials strength. These results contribute towards establishing the technical basis for qualifying LPBF 316H SS for advanced nuclear reactor applications.

36 MATERIALS SCIENCE↗

Dual-Fuel Ammonia Equivalence Ratio Sweep Data

Ammonia (NH3) has garnered significant interest as an alternative fuel for meeting international emissions reduction mandates in sectors with high weight and distance requirements, such as shipping. Technical barriers and unanswered questions remain on the combustion strategies that can maximize ammonia utilization and minimize emissions. Prior research studies at the US Department of Energy’s Oak Ridge National Laboratory have shown strong performance with NH3 under dual-fuel mode using conventional diesel combustion (CDC) manifold air pressure settings. Diesel airflow was initially used to simplify retrofitting (no turbocharger modification), which resulted in air-fuel equivalence ratios (λ) greater than 1.5. To characterize potential improvements in dual-fuel NH3 combustion performance at richer in-cylinder conditions, a global λ sweep compared the use of early (E-pilot) and late (L-pilot) single diesel injections. The experiments were conducted at 1200 RPM and 12.8 ± 0.2 bar (75 % load), and λ was varied by decreasing the commanded air flow to the engine at greater than 90 % ammonia energy substitution level. A diesel injection timing sweep was conducted for both the injection strategies at fixed λ, and the timing with the lowest engine-out N2O emissions was identified. The results indicated an optimal balance between CO2,eq and thermal efficiency benefits both E-pilot and L-pilot injection strategy cases compared with CDC at a λ of 1.4. The indicated nitrogen-based emissions exhibited a strong correlation to the ratio of CA5–50 and ignition delay for L-pilot, but no apparent trend emerged for the E-pilot injection strategy at the tested boundary conditions. This dataset includes the raw experimental data and documentation of the experiment conditions and methods.

30 DIRECT ENERGY CONVERSION↗

Numerical investigation of the influence of shear and thermal stratification on the wind turbine tip‐vortex stability

Summary The interaction between wind turbine wakes and atmospheric turbulence is characterised by complex dynamics. In this study, two major components of the atmospheric boundary layer dynamics have been isolated, namely, the mean velocity profile shear and the thermal stratification, to examine their impact on the near‐wake development by undertaking a series of highly resolved large‐eddy simulations. Subsequently, instantaneous flow fields are extracted from the simulations and used to conduct Fourier analysis and proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip‐vortex instability mechanisms. Our findings indicate that shear can significantly affect the breakup of the wind turbine tip‐vortices and the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near‐wake field. Finally, our analysis shows that the applied perturbation frequency determines the tip‐vortex breakup location as it controls the onset of the mutual inductance instability.

17 WIND ENERGY↗

Arctic Coastal Erosion: Modeling and Experimentation

Increasing Arctic coastal erosion rates have put critical infrastructure and native communities at risk while also mobilizing ancient organic carbon into modern carbon cycles. Although the Arctic comprises one-third of the global coastline and has some of the fastest eroding coasts, current tools for quantifying permafrost erosion are unable to explain the episodic, storm-driven erosion events. Our approach, mechanistically coupling oceanographic predictions with a terrestrial model to capture the thermo-mechanical dynamics of erosion, enables this much needed treatment of transient erosion events. The Arctic Coastal Erosion Model consists of oceanographic and atmospheric boundary conditions that force a coastal terrestrial permafrost environment in Albany (a multi-physics based finite element model). An oceanographic modeling suite (consisting of WAVEWATCH III, Delft3D-FLOW, and Delft3D-WAVE) produced time-dependent surge and run-up boundary conditions for the terrestrial model. In the terrestrial model, a coupling framework unites the mechanical and thermal aspects of erosion. 3D stress/strain fields develop in response to a plasticity model of the permafrost that is controlled by the frozen water content determined by modeling 3D heat conduction and solid-liquid phase change. This modeling approach enables failure from any allowable deformation (block failure, slumping, etc.). Extensive experimental work has underpinned the ACE Model development including field campaigns to measure in situ ocean and erosion processes, strength properties derived from thermally driven geomechanical experiments, as well as extensive physical composition and geochemical analyses. Combined, this work offers the most comprehensive and physically grounded treatment of Arctic coastal erosion available in the literature. The ACE model and experimental results can be used to inform scientific understanding of coastal erosion processes, contribute to estimates of geochemical and sediment land-to-ocean fluxes, and facilitate infrastructure susceptibility assessments.

54 ENVIRONMENTAL SCIENCES↗

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗

The influence of cooling rates on strain phase diagrams and domain structures of ferroelectric thin films: A case study of PbTiO 3

Strain engineering has been established as an effective approach to control phase equilibria, domain configurations, and functional properties of ferroelectric thin films. Temperature-strain phase diagrams have been used as powerful tools for providing insights into strain engineering. However, almost all existing phase diagrams established using the phase-field approach assume quenching conditions without considering actual cooling rates during the post-deposition annealing process of ferroelectric thin films. Within this work, we systematically investigate the influence of cooling rates on domain structures and the strain-phase diagram of ferroelectric thin films using phase-field simulations, taking PbTiO 3 thin films as a model system. We found that both the position of phase boundaries in the strain phase diagrams and the domain morphology are significantly influenced by the cooling rates. It is revealed that while the paraelectric-ferroelectric phase boundary remains invariant, the phase boundaries between single-phase and multi-phase regions tend to shift toward the corresponding multi-phase region as the cool rate reduces. Slow cooling generally leads to more ordered domain structures with increased domain size. Using the obtained equilibrium domain structures, we calculated effective thermal conductivities and found significant variations that can be tuned by the cooling rates. In conclusion, this work reveals an underexplored yet critical impact of cooling rates on phase equilibria and domain structures in ferroelectric thin films, which may inspire further fine-tuning of domains and domain walls in low-dimensional ferroelectrics for multifunctional applications.

Domain structures↗

Morphological-informed Thermal Property Prediction in the Engineering Domain

Using first-principle, atomic scale methods to predict thermal carrier (electron and phonon) behavior at nanoscale yields highly resolved thermal properties. Density functional theory simulations can model myriad effects on the transport of phonons and electrons (e.g., phonon and electron coupling, electron correlation, defect presence, carrier scattering), and yield thermal properties for small geometric domains such as thermal and electrical conductivity, or heat capacity. While these thermal properties are accurate in atomistic and nanoscopic systems, they must be scaled up to the microstructural domain to be a useful predictor for experimental basis, as the microscale is where myriad changes and physical phenomenon occur within a material (e.g., grain boundaries, precipitate aggregations, interfaces, defect clusters). My talk will discuss the necessity of thermal property predictions for nuclear fuels applications and more. Additionally I will discuss the necessity in developing multi-physics, multi-scale methods for future applications in predictions for not only nuclear fuels, but other materials applications.

36 MATERIALS SCIENCE↗

Towards High-Speed friction stir welding of 25 mm Thick AA2139-T8: tool innovation and process development

Joining thick plates (≥ 12 mm) of aluminum (Al) alloys is challenging due to high tool forces, uneven material flow, and non-uniform heat distribution through the material’s thickness. Recent advancements in tool design and welding parameters for friction stir welding (FSW) in thick plates butt joining have encountered a developmental plateau, highlighting the need for innovative approaches to overcome existing limitations. Here, this study focuses on the systematic development of single-pass and double-pass FSW processes for high-strength aluminum alloy AA2139-T8, to improve joint efficiency and enable high-speed welding capabilities. Experimental evidence is presented for a novel tool design with opposing pin threads, enabling high-speed (178 mm/min) single-pass friction stir butt welding of 25 mm thick AA2139-T8. A series of FSW trials was conducted both in air and with a trailing water spray, using steel backing plates (BPs) to investigate the impact of quenching and cooling rates on process response and joint performance. A joint efficiency of 83% was attained using the novel tool features and effectively controlling process forces and thermal boundary conditions.

AA2139-T8↗

Microstructural Changes to Thermally Sprayed Materials Subjected to Dynamic Compression

Dynamic compression of materials can induce a variety of microstructural changes. As thermally-sprayed materials have highly complex microstructures, the expected pressure at which changes occur cannot be predicted a priori. In addition, typical in-situ measurements such as velocimetry are unable to adequately diagnose microstructural changes such as failure or pore collapse. Quasi-isentropic compression experiments with sample recovery were conducted to examine microstructural changes in thermally sprayed tantalum and tantalum-niobium blends up to 8 GPa pressure. Spall fracture was observed in all tests, and post-shot pore volume decreased relative to the initial state. The blended material exhibited larger spall planes with fracture occurring at interphase boundaries. An estimate of the pressure at which pore collapse is complete was determined to be ~26 GPa for pure tantalum and ~19 GPa for the tantalumniobium blend under these loading conditions.

36 MATERIALS SCIENCE↗

An open source fast fluid dynamics model for data center thermal management

Although computational fluid dynamics (CFD) has been widely adopted to improve data center thermal management, the high computational demand limits its applications, such as multivariate optimal design and operation. Fast fluid dynamics (FFD), which has been applied for fast airflow simulation, shows great potential. However, few research applied FFD for optimal design and operation of data center thermal management. This research improves the FFD model for data centers and conducts a comprehensive evaluation and demonstration. First, the FFD model is improved by solving the advection and diffusion equations together using an upwind scheme instead of a semi-Lagrangian advection solver in the conventional FFD model. Second, new features for data centers are added, such as a pressure correction method to simulate plenum airflow and dynamic boundary conditions for IT racks. The new FFD model is first validated with two indoor environment cases and the results show that the new FFD model has slightly better overall prediction accuracy and faster speed compared to the conventional FFD model. It is also observed that both FFD models achieve acceptable accuracy, except for a few localized disparities with experimental data, which might be due to simplified handling of turbulence viscosity near the boundaries. Furthermore, validation with a real data center shows that the FFD model achieves a similar level of accuracy as CFD when compared to the experimental measurements with some level of uncertainties. It is then demonstrated for data center optimal design and operation, which saves 53.4–58.8% of annual energy while still meeting the thermal requirements. In conclusion, with a much faster speed and comparable accuracy compared to CFD, the FFD model parallelized on a graphics processing unit is promising for practical model-based data center early design and operation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The Effect of Membrane Casting Irregularities on Initial Fuel Cell Performance

Membrane electrode assemblies (MEA) for polymer electrolyte membrane fuel cells (PEMFC) may possess as-manufactured non-uniformities in any of its constituent components. This article studies casting irregularities located within the PEM to understand their potential impact on MEA initial performance. Membrane material was cast on a lab-scale polymer casting line operating either within or at the boundaries of the process window. The resulting membrane material was either pristine in the former case, or, in the latter case, contained air bubbles, cracks, and other irregularities, ranging in size from about 0.6 to 3 mm. Spatial polarization experiments were conducted using a 121-channel segmented fuel cell system, and thermal imaging was performed subsequently to map hydrogen crossover to spatial performance. While total-cell polarization data was minimally impacted by the irregularities, the spatial diagnostics showed local performance impacts that, in operation over time, could cause degradation in performance or earlier failure of the MEA. Such impacts could lead to these irregularities being classified as defects, i.e., manufacturing variations that should be identified and not included in a fuel cell stack. Classification of irregularities as defects will ultimately assist the industry by contributing to the development of threshold detection limits for in-line quality control diagnostics.

30 DIRECT ENERGY CONVERSION↗

The radiation instability of thermally stable nanocrystalline platinum gold

Here, recent experimentally validated alloy design theories have demonstrated nanocrystalline binary alloys that are stable against thermally induced grain growth. An open question is whether such thermal stability also translates to stability under irradiation. In this study, we investigate the response to heavy ion irradiation of a nanocrystalline platinum gold alloy that is known to be thermally stable from previous studies. Heavy ion irradiation was conducted at both room temperature and elevated temperatures on films of nanocrystalline platinum and platinum gold. Using scanning/transmission electron microscopy equipped with energy-dispersive spectroscopy and automated crystallographic orientation mapping, we observe substantial grain growth in the irradiated area compared to the controlled area beyond the range of heavy ions, as well as compositional redistribution under these conditions, and discuss mechanisms underpinning this instability. These findings highlight that grain boundary stability against one external stimulus, such as heat, does not always translate into grain boundary stability under other stimuli, such as displacement damage.

36 MATERIALS SCIENCE↗

Dynamical System Scaling of a Thermocline Thermal Storage System in the Thermal Energy Distribution System (TEDS) Facility

The purpose of this study was to develop a process to convert input signals from one facility into another by reflecting geometric and environmental settings. The Dynamic Energy Transport and Integration Laboratory (DETAIL) is one facility in development that aims to emulate the daily interactions among power production industry systems and be capable of receiving real-time data from those systems as inputs. To convert signals and ensure that the temporal sequences and magnitudes reflect the laboratory settings, the ability to scale and project data is essential. To demonstrate this ability, Dynamical System Scaling (a methodology that enables systems to scale and project or extrapolate datasets to desired environments while conserving the observed transient behavior based on first principles) was applied to DETAIL’s thermocline thermal storage system in the Thermal Energy Distribution System. The thermocline system was successfully scaled and a test case was conducted to generate a doubly accelerated energy charge and discharge in reference to past experimental data from the facility. The accelerated data were determined as able to conserve the amount of energy stored and the associated test boundary conditions were charge line maximum temperature, charge line velocity, and thermocline maximum temperature at 354 °C, 0.458 m/s, and 418 °C, respectively. The research results represented a case that required signals to be accelerated without altering the stored energy.

Yoshiura, Ramon (ORCID:0000000319781313)↗

Validation of Oregon State University High Temperature Test Facility Experiments Using Pronghorn

The OSU High Temperature Test Facility is a quarter-scale diameter, 1/64 scale volume test facility meant to replicate thermophysical phenomena in the prototypical General Atomics Modular High Temperature Gas Reactor. Tests pertaining to conduction cooldown events were performed from 2016-2019, providing a large database by which computational methods that are applicable to different length scales can be validated. One of these codes is Pronghorn, which is a coarse-meshed, porous-based subchannel thermal hydraulics code based on the MOOSE application, with the intention of better capturing the physics of both conduction and convection heat transfer within the OSU HTTF core. The goal of this summer project is to develop the framework by which Pronghorn can perform validation exercises of the HTTF core for benchmarking, by generating a mesh appropriate to the geometry of the HTTF core, developing input decks that accurately capture the initial and boundary conditions, materials, and relevant equations to the physics seen in the HTTF core, and using a postprocessor to compare simulation results to various experimental data. While validation of codes is a multi-year project, a mesh has been generated and tested in Pronghorn that meets mass conservation and basic heat transfer principles. The next step is to accurate depict the fluid inlet and outlet boundary conditions, which will be performed using computational fluid dynamics software.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Magnon orbital angular momentum of ferromagnetic honeycomb and zigzag lattice models

By expanding the gauge 𝜆 𝑛 ⁡(𝐤) for magnon band 𝑛 in harmonics of momentum 𝐤=(𝑘,𝜙), we demonstrate that the only observable component of the magnon orbital angular momentum 𝑂 𝑛 ⁡(𝐤) is its angular average over all angles 𝜙, denoted by 𝐹 𝑛 ⁡(𝑘). Although 𝐹 𝑛 ⁡(𝑘) vanishes for antiferromagnetic honeycomb and zigzag (0<𝐽 1 <𝐽 2 ) lattices, it is nonzero for the ferromagnetic (FM) versions of those lattices in the presence of Dzyaloshinskii-Moriya interactions. For a FM zigzag model with equal exchange interactions 𝐽 1⁢𝑥 and 𝐽 1⁢𝑦 along the 𝑥 and 𝑦 axes, the magnon bands are degenerate along the boundaries of the Brillouin zone with 𝑘 𝑥 −𝑘 𝑦 =±𝜋/𝑎 and the Chern numbers 𝐶 𝑛 are not well defined. However, a revised model with 𝐽 1⁢𝑦 ≠𝐽 1⁢𝑥 lifts those degeneracies and produces well-defined Chern numbers of 𝐶𝑛=±1 for the two magnon bands. When 𝐽 1⁢𝑦 =𝐽 1⁢𝑥 , the thermal conductivity 𝜅 𝑥⁢𝑦⁡ (𝑇) of the FM zigzag lattice is largest for 𝐽2/𝐽1>6 but is still about four times smaller than that of the FM honeycomb lattice at high temperatures. Due to the removal of band degeneracies, 𝜅𝑥⁢𝑦⁡(𝑇) is slightly enhanced when 𝐽 1⁢𝑦 ≠𝐽 1⁢𝑥 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Physics-Informed Neural Networks for Heat Transfer Problems

Abstract Physics-informed neural networks (PINNs) have gained popularity across different engineering fields due to their effectiveness in solving realistic problems with noisy data and often partially missing physics. In PINNs, automatic differentiation is leveraged to evaluate differential operators without discretization errors, and a multitask learning problem is defined in order to simultaneously fit observed data while respecting the underlying governing laws of physics. Here, we present applications of PINNs to various prototype heat transfer problems, targeting in particular realistic conditions not readily tackled with traditional computational methods. To this end, we first consider forced and mixed convection with unknown thermal boundary conditions on the heated surfaces and aim to obtain the temperature and velocity fields everywhere in the domain, including the boundaries, given some sparse temperature measurements. We also consider the prototype Stefan problem for two-phase flow, aiming to infer the moving interface, the velocity and temperature fields everywhere as well as the different conductivities of a solid and a liquid phase, given a few temperature measurements inside the domain. Finally, we present some realistic industrial applications related to power electronics to highlight the practicality of PINNs as well as the effective use of neural networks in solving general heat transfer problems of industrial complexity. Taken together, the results presented herein demonstrate that PINNs not only can solve ill-posed problems, which are beyond the reach of traditional computational methods, but they can also bridge the gap between computational and experimental heat transfer.

Engineering↗