Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat transfer coefficient”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Group-III nitride devices and systems on IBAD-textured substrates

A multilayer structure including a hexagonal epitaxial layer, such as GaN or other group III-nitride (III-N) semiconductors, a <111> oriented textured layer, and a non-single crystal substrate, and methods for making the same. The textured layer has a crystalline alignment preferably formed by the ion-beam assisted deposition (IBAD) texturing process and can be biaxially aligned. The in-plane crystalline texture of the textured layer is sufficiently low to allow growth of high quality hexagonal material, but can still be significantly greater than the required in-plane crystalline texture of the hexagonal material. The IBAD process enables low-cost, large-area, flexible metal foil substrates to be used as potential alternatives to single-crystal sapphire and silicon for manufacture of electronic devices, enabling scaled-up roll-to-roll, sheet-to-sheet, or similar fabrication processes to be used. The user is able to choose a substrate for its mechanical and thermal properties, such as how well its coefficient of thermal expansion matches that of the hexagonal epitaxial layer, while choosing a textured layer that more closely lattice matches that layer. Electronic devices such as LEDs can be manufactured from such structures. Because the substrate can act as both a reflector and a heat sink, transfer to other substrates, and use of external reflectors and heat sinks, is not required, greatly reducing costs. Large area devices such as light emitting strips or sheets may be fabricated using this technology.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Group-III nitride devices and systems on IBAD-textured substrates

A multilayer structure including a hexagonal epitaxial layer, such as GaN or other group III-nitride (III-N) semiconductors, a <111>oriented textured layer, and a non-single crystal substrate, and methods for making the same. The textured layer has a crystalline alignment preferably formed by the ion-beam assisted deposition (IBAD) texturing process and can be biaxially aligned. The in-plane crystalline texture of the textured layer is sufficiently low to allow growth of high quality hexagonal material, but can still be significantly greater than the required in-plane crystalline texture of the hexagonal material. The IBAD process enables low-cost, large-area, flexible metal foil substrates to be used as potential alternatives to single-crystal sapphire and silicon for manufacture of electronic devices, enabling scaled-up roll-to-roll, sheet-to-sheet, or similar fabrication processes to be used. The user is able to choose a substrate for its mechanical and thermal properties, such as how well its coefficient of thermal expansion matches that of the hexagonal epitaxial layer, while choosing a textured layer that more closely lattice matches that layer. Electronic devices such as LEDs can be manufactured from such structures. Because the substrate can act as both a reflector and a heat sink, transfer to other substrates, and use of external reflectors and heat sinks, is not required, greatly reducing costs. Large area devices such as light emitting strips or sheets may be fabricated using this technology.

Matias, Vladimir↗
DOE OSTI2021

Transformational challenge reactor analysis to inform preconceptual core design decisions: Sensitivity study of transient analysis in a hydride-moderated microreactor

The Transformational Challenge Reactor (TCR) program aims to demonstrate a revolutionary design approach enabled by advanced manufacturing and data analytics in the design of nuclear reactors. This article discusses scoping analyses of preconceptual designs to inform TCR design decisions and the evaluation of sensitivities and uncertainties on postulated transient scenarios. The applicability of the systems codes TRACE and RELAP5-3D to TCR transient analysis are examined, and RELAP5-3D models are used to examine the transient response of two candidate core designs at multiple power levels. Then, the uncertainty quantification code RAVEN is used to quantify the effect of several design parameters on reactivity-initiated accident (RIA) progression at hot zero power (HZP) and hot full power (HFP) as well as to assess the impact of uncertainties in transient parameters for the pressurized loss of forced cooling (PLOFC).When results were compared, TRACE and RELAP5-3D showed good agreement in their ability to predict system behavior, but RELAP5-3D calculations were closer to analytical predictions for the RIA. Furthermore, models for a PLOFC accident in two designs (a UO2 and TRISO core) at multiple power levels showed greater temperature margins for the TRISO core at all power levels. Using this information, along with other scoping analyses and constraints, the TCR design team selected a power level of 3 MWth and a TRISO-based core design. For this design RAVEN was applied to vary RIA parameters in RELAP5-3D models at HZP and HFP to understand the effect on figures of merit such as peak power, fuel and coolant temperature, and energy deposition. This sensitivity study found that the inserted reactivity worth was the most important parameter controlling all figures of merit, but for insertion up to 1.5$ no failure of TCR fuel is anticipated. For constant reactivity insertions, the magnitude of the fuel temperature coefficient was found to have the greatest effect on all figures of merit under most circumstances. These results not only demonstrate the anticipated robust safety of the proposed TCR fuel form but also provide a reference for future metal-hydride moderated systems to understand RIA behavior. In the PLOFC, the impact of heat transfer enhancement due to wavy flow channel effects dominated the variance in peak temperatures, and variations in heat exchanger elevation provided the greatest control on natural circulation flow rate. Additionally, no fuel particle failure is anticipated in the PLOFC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Performance Assessment of the Overall Building Envelope Thermal Performance—Building Envelope Performance (BEP) Metric

Today, to describe the thermal performance of the building envelope and its components we use a variation of metrics; such as, R-value, ACH (air exchange rate per hour), SHGC (solar heat gain coefficient) of windows, U-factor etc. None of these performance indicators is meant to represent the overall thermal performance. In this paper, such a metric is introduced, the BEP (building envelope performance) value. Unlike the thermal resistance, typically expressed as an R-value, the BEP-value considers additional elements of heat transfer that affect the energy demand of the building because of exterior and interior (solar) thermal loads: conductive and radiant heat transfer, and air infiltration. To demonstrate BEP’s utility, validation studies were carried out by comparing the BEP-value to theoretical results using whole building energy simulation tools such as EnergyPlus and WUFI Plus. Results show that BEP calculations are comparable to calculations made using these simulation tools and that unlike other similar metrics, the BEP-value accounts for all heat transfer mechanisms that are relevant for the overall energy performance of the building envelope. The BEP-value thus allows comparing envelopes of buildings with different use types in a fair and realistic manner.

42 ENGINEERING↗
DOE OSTI2024

Unraveling Interdiffusion Phenomena and the Role of Nanoscale Diffusion Barriers in the Copper–Gold System

Diffusion is one of the most fundamental concepts in materials science, playing a pivotal role in materials synthesis, forming, and degradation. Of particular importance is solid state interdiffusion of metals which defines the usable parameter space for material combinations in the form of alloys. This parameter space can be explored on the macroscopic scale by using diffusion couples. However, this method reaches its limit when going to low temperatures, small scales, and when testing ultrathin diffusion barriers. Therefore, this work transfers the principle of the diffusion couples to small scales by using core–shell nanowires and in situ heating. This allows us to delve into the interdiffusion dynamics of copper and gold, revealing the interplay between diffusion and the disorder–order phase transition. Our in situ TEM experiments in combination with chemical mapping reveal the interdiffusion coefficients of Cu and Au at low temperatures and highlight the impact of ordering processes on the diffusion behavior. The formation of ordered domains within the solid-solution is examined using high-resolution imaging and nanodiffraction including strain mapping. In addition, we examine the effectiveness of ultrathin Al 2 O 3 barrier layers to control interdiffusion of the diffusion couple. Our findings indicate that a 5 nm thick layer serves as an efficient diffusion barrier. Furthermore, this research provides valuable insights into the interdiffusion behavior of Cu and Au on the nanoscale, offering potential applications in the development of miniaturized integrated circuits and nanodevices.

alloys↗
DOE OSTI2023

Solar thermal aerogel receiver and materials therefor

A silica aerogel having a mean pore size less than 5 nm with a standard deviation of 3 nm. The silica aerogel may have greater than 95% solar-weighted transmittance at a thickness of 8 mm for wavelengths in the range of 250 nm to 2500 nm, and a 400° C. black-body weighted specific extinction coefficient of greater than 8 m2/kg for wavelengths of 1.5 μm to 15 μm. Silica aerogel synthesis methods are described. A solar thermal aerogel receiver (STAR) may include an opaque frame defining an opening, an aerogel layer disposed in the opaque frame, with at least a portion of the aerogel layer being proximate the opening, and a heat transfer fluid pipe in thermal contact with and proximate the aerogel layer. A concentrating solar energy system may include a STAR and at least one reflector to direct sunlight to an opening in the STAR.

Chen, Gang↗
DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Solar thermal aerogel receiver and materials therefor

A silica aerogel having a mean pore size less than 5 nm with a standard deviation of 3 nm. The silica aerogel may have greater than 95% solar-weighted transmittance at a thickness of 8 mm for wavelengths in the range of 250 nm to 2500 nm, and a 400° C. black-body weighted specific extinction coefficient of greater than 8 m 2 /kg for wavelengths of 1.5 μm to 15 μm. Silica aerogel synthesis methods are described. A solar thermal aerogel receiver (STAR) may include an opaque frame defining an opening, an aerogel layer disposed in the opaque frame, with at least a portion of the aerogel layer being proximate the opening, and a heat transfer fluid pipe in thermal contact with and proximate the aerogel layer. A concentrating solar energy system may include a STAR and at least one reflector to direct sunlight to an opening in the STAR.

Chen, Gang↗
DOE OSTI2025

Estimating Soil Thermal Inertia Profiles From the Passive Equilibration of a Temperature Probe

Knowledge of the distribution of soil thermal properties is important for understanding subsurface hydrological and biogeochemical processes. This study describes and evaluates quick thermal profiling (QTP), a new measurement technique aimed at providing rapid, depth-resolved measurements of soil thermal inertia at numerous locations across the landscape. A cylindrical probe with temperature sensors at multiple depths is quickly inserted into the ground, and soil thermal inertia is estimated from how quickly the probe temperature equilibrates with the soil. To this end, a finite volume heat transfer model is used to generate temperature equilibration time series across combinations of controlling factors, and a gridded search inversion approach is applied to infer soil thermal inertia. Field tests in the Arctic indicate that QTP measurements have a minimum uncertainty of 0.14 J m −2 K −1 s −1/2 and covary with dual-probe heat pulse thermal analyzer measurements (concordance correlation coefficient = 0.56) with a root-mean-square error of 0.40 J m −2 K −1 s −1/2 . Besides demonstrating the value of QTP for estimating thermal inertia, this study identifies various sources of measurement uncertainty, particularly probe-soil contact resistance and frictional heating. Further, analysis of soil samples indicates that thermal inertia can be used to estimate thermal conductivity and dry bulk density in the studied area, although such inferences are highly site-specific. Overall, the QTP method holds promise to generate thermal inertia data products and to complement other characterization approaches for advancing understanding of soil properties across far more locations than is currently possible.

Lamb, J. R. [Lawrence Berkeley National Laboratory↗
DOE OSTI2021

Saturation physics of threshold heat-flux reduction

The saturation physics of ion-temperature-gradient-driven turbulence is examined in relation to the temperature-gradient variation of the heat flux, which can exhibit an upshift of the critical gradient for significant flux relative to the linear instability threshold. Gyrokinetic measurements of saturation properties and spectral energy transfer, which will be defined in Sec. II, are presented, indicating that the physics of saturation is fundamentally unchanged on either side of the upshifted gradient. To analyze heat transport below and above the upshifted critical gradient, a fluid model for toroidal ion-temperature-gradient turbulence is modified to include the kinetic instability threshold. The model and the heat flux are rendered in the eigenmode decomposition to track the dominant mode-coupling channel of zonal-flow-catalyzed transfer to a conjugate stable mode. Given linear and nonlinear symmetries, the stable mode level and the cross-correlation of the unstable and stable mode amplitudes are related to the unstable mode level via linear physics. Here, the heat flux can then be written in terms of the unstable-mode level, which through a nonlinear balance depends on the eigenmode-dependent coupling coefficients and the triplet correlation time of the dominant coupled modes. Resonance in these quantities leads to suppressed heat flux above the linear threshold, with a nonlinear upshift of the critical gradient set by the resonance broadening of a finite perpendicular wavenumber and collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fouling modeling and prediction approach for heat exchangers using deep learning

In this article, we develop a generalized and scalable statistical model for accurate prediction of fouling resistance using commonly measured parameters of industrial heat exchangers. This prediction model is based on deep learning where a scalable algorithmic architecture learns non-linear functional relationships between a set of target and predictor variables from large number of training samples. Here, the efficacy of this modeling approach is demonstrated for predicting fouling in an analytically modeled cross-flow heat exchanger, designed for waste heat recovery from flue-gas using room temperature water. The performance results of the trained models demonstrate that the mean absolute prediction errors are under 10 –4 KW –1 for flue-gas side, water side and overall fouling resistances. The coefficients of determination (R 2 ), which characterize the goodness of fit between the predictions and observed data, are over 99%. Even under varying levels of measurement noise in the inputs, we demonstrate that predictions over an ensemble of multiple neural networks achieves better accuracy and robustness to noise. We find that the proposed deep-learning fouling prediction framework learns to follow heat exchanger flow and heat transfer physics, which we confirm using locally interpretable model agnostic explanations around randomly selected operating points. Overall, we provide a robust algorithmic framework for fouling prediction that can be generalized and scaled to various types of industrial heat exchangers.

42 ENGINEERING↗
DOE OSTI2022

Understanding the low-temperature chemistry of 1,2,4-trimethylbenzene

1,2,4-trimethylbenzene is an important representative aromatic component of gasoline/diesel/jet fuels and thus it is necessary to understand its low-temperature chemistry. In this paper, ignition delay times (IDTs) of both 1,2,4-trimethylbenzene (124TMB) and its blends with n -heptane were measured at engine-like conditions using both a high-pressure shock tube and a rapid compression machine for fuel in ‘air’ mixtures at pressures of 10 and 30 atm and at temperatures in the range 600 – 1100K. The experiments in this study show for the first time that 124TMB presents a two-stage ignition behavior at engine relevant conditions. Blending n -heptane with 124TMB can significantly increase mixture reactivity at temperatures below 1000K. A new detailed mechanism has been developed to simulate the experimentally measured IDT data. The mechanism can capture well the two-stage ignition behavior as well as the ignition delays at different pressures, equivalence ratios over a wide temperature range, for both pure fuels and their blended mixtures. Flux analyses show that the benzylic radicals (formed via H-atom abstraction from the methyl groups ortho-sites on 124TMB) can add to O 2 forming $\mathrm{R\dot{O}_2}$ radicals, which can isomerize to $\dot{Q}$OOH by intramolecular H-atom transfer from the ortho- methyl group and these $\dot{Q}$OOH radicals undergo a second addition to O 2 . This is analogous to the chain branching reaction pathways of alkanes. The chain branching reaction pathways are responsible for the first-stage heat release of 124TMB. In conclusion, the competitions between chain branching and both chain propagating and chain termination reaction pathways lead to a less pronounced negative temperature coefficient (NTC) behavior for 124TMB oxidation, compared to two-stage ignition behavior observed for alkanes and other fuels.

1,2,4-trimethylbenzene↗
DOE OSTI2022

An ab initio molecular dynamics investigation of the thermophysical properties of molten NaCl-MgCl 2

Molten salts have many applications in the nuclear and solar energy industries for thermal storage and heat transfer applications. However, there is a knowledge gap in molten salt thermophysical properties which hinders the technical readiness level of molten salt applications, especially in the nuclear industry. A common method of investigating new materials is through ab initio Molecular Dynamics (AIMD) simulations which is an effective tool to investigate structural and thermophysical properties at realistic temperatures. NaCl-MgCl 2 is an inexpensive salt that is a good candidate for use as a heat transfer medium in solar power applications or in the secondary loop of a nuclear reactor. In this article, the thermophysical properties of NaCl-MgCl 2 are computed via AIMD calculations to supplement the limited experimental studies in the literature. Here a wide range of compositions and temperatures for the pseudo-binary NaCl-MgCl 2 were used to calculate the density, heat capacity, compressibility, enthalpy of mixing, and volumetric thermal expansion coefficient. AIMD is shown to accurately model the densities of molten NaCl-MgCl 2 as there is good agreement with the available literature. This work observed a transition to a monotonic increase of the density with respect to MgCl 2 composition occurring above 1100 K. The heat capacity values increase uniformly with respect to concentration of MgCl 2 at a rate of 2.85 J/mol-K per 10 mol% of MgCl 2 . Select thermophysical properties are fit to a Redlich-Kister expansion for utilization in multiphysics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Influence of Reservoir Convection on Heat Extraction with Closed-Loop Geothermal Systems

Over the past few years, advancements in closed-loop geothermal systems (CLGS), also called advanced geothermal systems (AGS), have sparked a renewed interest in these types of designs. CLGS have certain advantages over traditional and enhanced geothermal systems (EGS), including not requiring in-situ reservoir permeability, conservation of the circulating fluid, and allowing for different fluids, including working fluids directly driving a turbine at the surface. CLGS may be attractive in environments where water resources are limited, rock contaminants must be avoided, and stimulation treatments are not available (e.g., due to regulatory or technical reasons). Despite these advantages, CLGS have some challenges, including limited surface area for heat transfer and requiring long wellbores and laterals to obtain multi-MW output in conduction-only reservoirs. CLGS have been investigated in conduction-only systems. In this paper, we explore the impact of both forced and natural convection on the levels of heat extraction with a CLGS deployed in a hot wet rock reservoir. We bound potential benefits of convection by investigating liquid reservoirs over a range of natural and forced convective coefficients. Additionally, we investigate the effects of permeability, porosity, and geothermal temperature gradient in the reservoir on CLGS outputs. Reservoir simulations indicate that reservoir permeabilities of at least ~100 mD are required for natural convection to increase the heat output with respect to a conduction-only scenario. The impact increases with increasing reservoir temperature. When subject to a forced convection flow field, Darcy velocities of at least 10-7 m/s are required to obtain an increase in heat output.

advanced geothermal systems↗
DOE OSTI2024

Significantly Promoting the Thermal Conductivity and Machinability of Negative Thermal Expansion Alloy via In Situ Precipitation of Copper Networks

Rapid advancements in electronic devices yield an urgent demand for high-performance electronic packaging materials with high thermal conductivity, low thermal expansion, and great mechanical properties. However, it is a great challenge for current design philosophies to fulfill all the requirements simultaneously. Here, an effective strategy is proposed for significantly promoting the thermal conductivity and machinability of negative thermal expansion alloy (Zr,Nb)Fe 2 through eutectic precipitation of copper networks. The eutectic dual-phase alloy exhibits an isotropic chips-matched thermal expansion coefficient and a thermal conductivity enhancement exceeding 200% compared with (Zr,Nb)Fe 2 , along with an ultimate compressive strength of 550 MPa. The addition of copper reorganizes the composition of (Zr,Nb)Fe 2 , which smooths the magnetic transition and shifts it toward higher temperature, resulting in linear low thermal expansion in a wide temperature range. The highly fine eutectic copper lamellae construct high thermal conductivity networks within (Zr,Nb)Fe 2 , serving as highways for heat transfer electrons and phonons. The in situ forming of eutectic copper lamellae also facilitates the mechanical properties by enhancing interfacial bonding and bearing additional stress after yielding of (Zr,Nb)Fe 2 . This work provides a novel strategy for promoting thermal conductivity and mechanical properties of negative thermal expansion alloys via eutectic precipitation of copper networks.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Conversion Kinetics and Ionic Conductivity in Na-β”-Alumina + YSZ (Naβ”AY) Sodium Solid Electrolyte via Vapor Phase Conversion Process

Sodium ion batteries have been receiving increasing attention and may see potential revival in the near future, particularly in large-scale grid energy storage coupling with wind and solar power generation, due to the abundant sodium resources, low cost, and sufficiently high energy density. Among the known sodium ion conductors, the Na-β”-alumina electrolyte remains highly attractive because of its high ionic conductivity. This study focuses on the vapor phase synthesis of a Na-β”-Alumina + YSZ (Naβ”AY) composite sodium electrolyte, which has higher mechanical strength and stability than conventional single phase β”-Alumina. The objectives are the measurement of conversion kinetics through a newly developed weight-gain based model and the determination of sodium ionic conductivity in the composite electrolyte. Starting samples contained ~70 vol% α-Alumina and ~30 vol% YSZ (3 mol% Y 2 O 3 stabilized Zirconia) with and without a thin alumina surface layer made by sintering in air at 1600 °C. The sintered samples were placed in a powder of Na-β”-alumina and heat-treated at 1250 °C for various periods. Sample dimensions and weight were measured as a function of heat treatment time. The conversion of α-Alumina in the α-Alumina + YSZ composite into Naβ”AY occurred by coupled diffusion of sodium ions through Na-β”-alumina and of oxygen ions through YSZ, effectively diffusing Na2O. From the analysis of the time dependence of sample mass and dimensions, the effective diffusion coefficient of Na 2 O through the sample, D eff , was estimated to be 1.74 x 10 -7 cm2 s -1 , and the effective interface transfer parameter, k eff , was estimated as 2.33 x 10 -6 cm s -1 . By depositing a thin alumina coating layer on top of the bulk composite, the chemical diffusion coefficient of oxygen through single phase Na-β”-alumina was estimated as 4.35 x 10 -10 cm 2 s -1 . An AC impedance measurement was performed on a fully converted Naβ”AY composite, and the conductivity of the composite electrolyte was 1.3 x 10 -1 S cm -1 at 300 °C and 1.6 x 10 -3 S cm -1 at 25 °C, indicating promising applications in solid state or molten salt batteries at low to intermediate temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synergistic learning with multi-task DeepONet for efficient PDE problem solving

Multi-task learning (MTL) is an inductive transfer mechanism designed to leverage useful information from multiple tasks to improve generalization performance compared to single-task learning. It has been extensively explored in traditional machine learning to address issues such as data sparsity and overfitting in neural networks. In this work, we apply MTL to problems in science and engineering governed by partial differential equations (PDEs). However, implementing MTL in this context is complex, as it requires task-specific modifications to accommodate various scenarios representing different physical processes. To this end, we present a multi-task deep operator network (MT-DeepONet) to learn solutions across various functional forms of source terms in a PDE and multiple geometries in a single concurrent training session. We introduce modifications in the branch network of the vanilla DeepONet to account for various functional forms of a parameterized coefficient in a PDE. Additionally, we handle parameterized geometries by introducing a binary mask in the branch network and incorporating it into the loss term to improve convergence and generalization to new geometry tasks. Our approach is demonstrated on three benchmark problems: (1) learning different functional forms of the source term in the Fisher equation; (2) learning multiple geometries in a 2D Darcy Flow problem and showcasing better transfer learning capabilities to new geometries; and (3) learning 3D parameterized geometries for a heat transfer problem and demonstrate the ability to predict on new but similar geometries. Finally, our MT-DeepONet framework offers a novel approach to solving PDE problems in engineering and science under a unified umbrella based on synergistic learning that reduces the overall training cost for neural operators.

42 ENGINEERING↗
DOE OSTI2026

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗