Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential geometry”

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

Characterization Strategy of a Limited Number of Radioactive Waste Effluent Tanks in View of Decommissioning a Large Number of Tanks - 20286

In order to organize the D and D of the former effluent treatment facility at the CEA Cadarache Center (France), the key data to gather through the facility characterization are the following:: volumes and masses, liquid/solid distribution and specific activity. By studying the treatment processes and the history of the facility, families of tanks have been defined and a reference tank chosen for investigation in each family. The intrusive investigations (video, residue thickness, dose rate..), the gamma spectrometry measurements carried out on site on the residue samples and radiochemical results collected in laboratories, led us to conclude that there was a relative radiological homogeneity (for specific activity and Am-241/Cs-137 ratio) within each tank. The next step consisted in extrapolating the specific activity evaluated for the reference tanks, to the other tanks of each family by using data such as geometry, dose rate or the function of the tanks within the process. Thus, a cartography of the 130 tanks of the facility is available. Finally, the conclusions drawn from the radiological homogeneity of the tanks and their affiliation to homogeneous, well-identified tank clusters lead to a simple and efficient operational approach for the D and D of the 130 tanks of the facility. Indeed, because of the residue characteristics (homogeneous activity and strong differentiation of Cs/Am ratio), a sample collected when opening the tank, and measured on site by gamma spectrometry will allow us to validate the operational design and the emptying techniques defined during the D and D study phase thanks to the characterization of a limited number of reference tanks. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Genetic algorithm optimization of nuclear criticality experiment for reduction of intermediate-energy 239 Pu nuclear data uncertainties

Nuclear criticality experiments are conducted to investigate specific nuclear data important for safe handling and storage of fissile materials, reactor design and operation, and the validation of radiation transport codes. Incorrect or uncertain nuclear data can prohibitively impact operational safety limits, reactor licensing, and predictive simulation capability; therefore, integral measurements from criticality experiments are necessary and should be performed frequently. To maximize the impact of the integral measurements, it is important to consider experiment geometry, material selection, and component dimensions. When taking these considerations into account, the experiment design process becomes iterative and very time intensive. This work utilizes a genetic algorithm to efficiently explore potential nuclear criticality experiment designs for the Laboratory Directed Research & Development project PARADIGM (PARallel Approach of Differential and InteGral Measurements) at Los Alamos National Laboratory. In this paper, the building blocks of the genetic algorithm are discussed in detail, the genetic algorithm methodology is verified, and the genetic algorithm is used to produce three candidate experiment models for the final PARADIGM design. The three candidate models produced by the genetic algorithm consist of copper-reflected assemblies containing 14 repeating units of alumina, graphite, boron, and plutonium plates. Furthermore, in addition to the optimization results, final design considerations are also discussed for designs with a height and/or weight very close to or slightly above assembly machine operational limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Ceramic encapsulated metal phase change material for high temperature thermal energy storage

Thermal energy storage (TES) is a broad-based technology for reducing CO₂ emissions and advancing concentrating solar, fossil, and nuclear power through improvements in efficiency and economics. Phase change materials (PCMs) are of interest as TES media because of their ability to store large amounts of heat in relatively small volumes. The volume expansion during a phase change, typically between a solid and liquid, can cause breakage of protective coatings. This paper reports on the fabrication of a ceramic encapsulated metal (CEM) high temperature TES technology using a rotary calcining furnace and a fluidized bed chemical vapor deposition coating technique. Aluminum beads were chosen as the PCM because Al has a high melting point (660 °C), low cost, high heat of fusion, and an ability to form a thin, strong alumina layer capable of supporting the Al melt for subsequent processing. Quite remarkably, this study shows that 1 mm diameter Al can be fluidized up to at least 1500 °C in an appropriate atmosphere while maintaining a spheroid geometry. This allowed for producing a first of a kind CEM whereby Al particles were encapsulated in pyro-carbon (PyC) and high purity, dense chemical vapor deposited SiC. The CEM with a PyC only coating was exposed to thermal cycling to test the performance with a differential scanning calorimeter; the melting point and latent heat were measured to be 648.4 ± 2.8 °C and 293.3 ± 6.2 J/g respectively. It was demonstrated that the CEM design is possible to produce, laying the foundation for manufacturing of high temperature, tunable, TES media.

25 ENERGY STORAGE↗
DOE OSTI2021

Sr and Co Vapor-Phase Transport from LSCF Cathodes

A series of tests were performed to differentiate between surface diffusion and vapor-phase diffusion of chemical species from LSCF cathode material at typical sintering temperatures (up to 1100°C) in SOFC production. A GDC source substrate was printed with LSCF and separated by an air gap from a YSZ target substrate. Various geometries with long surface paths were employed to reduce the possibility of surface transport. Sr and Co were detected on the target substrates via energy dispersive spectroscopy (EDS) and x-ray photoelectron spectroscopy (XPS), including spatially resolved XPS. These results support a vapor-phase transport mechanism for Sr and Co. Sintering conditions and barrier layer requirements may need to be revisited to prevent the formation of undesired strontium zirconate at cathode/electrolyte interface due to vapor-phase transport.

vapor phase transport, SOFC cathode, Sr transport↗
DOE OSTI2025

Fragmentation dynamics through geometrical distortion in low-energy electron attachment to carbon disulfide

The dissociation dynamics of dissociative electron attachment to carbon disulfide is studied using the velocity map imaging technique. Ion yields of the S− and CS− fragments are measured as a function of incident electron energy between 2 and 12 eV. Energy and angle-differential yields of the S− fragments are analyzed to investigate the dynamics of the CS2− anion around the 3.8-eV peak. The resonances involved in this fragmentation are identified and characterized in calculations using the well-known R-matrix method. Based on our measurements and the geometry dependence of the calculated resonances, we propose that a geometrical distortion of the transient negative-ion state of CS2 is necessary for S-CS bond dissociation to occur.

Paul, Anirban↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2022

Simulation of CEFR neutronic start-up tests with FENNECS

This paper presents simulation results of selected Neutronic Start-up Tests of the China Experimental Fast Reactor (CEFR) obtained by the neutronics code FENNECS that have been performed within the frame of the IAEA Coordinated Research Program I31032. The Finite Element Neutronics code FENNECS is developed at GRS and solves the few-group steady-state and transient diffusion equation using a Galerkin-based finite element approach. Its main purpose is the safety assessment of Small Modular Reactors and Micro Reactors with complex geometry (e.g., rotating control drums) which gain increased interest internationally. Serpent has been applied to create reference models of the CEFR and cross-section libraries in 10 energy groups for FENNECS. Using these libraries, the following experiments have been simulated with FENNECS: net criticality, control rod integral and differential worth, void reactivity effect, subassembly exchange reactivity effects and reaction rate distribution. The obtained satisfactory agreements with the measurements represent a valuable contribution to the validation of FENNECS. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Whole core analysis of OECD/NEA TVA-WB1 benchmark exercise 1 using continuous-energy Monte Carlo code MVP with JENDL-4.0u and ENDF/B-VII.0

The TVA WB1 benchmark problem was precisely modeled in accordance with the geometry and material specifications, and benchmark exercise 1 was solved with the MVP code and the nuclear data libraries based on JENDL-4.0u (JENDL-4.0 updated file) and ENDF/B-VII.0. The criticality, the control rod bank worth, the isothermal temperature coefficient (ITC), the integral rod worth of bank D, and differential soluble boron worth (DBW) were obtained. The criticality and the control rod bank worth were compared with the measured data and published calculation results with other codes, Serpent and KENO. The criticality and the control rod bank worth calculated by MVP and the libraries agree with the measured and the published calculation results within the criteria specified in the ANSI/ANS standard. From these results, the calculation model for the MVP and the libraries was verified. The other parameters such as ITC, the integral rod worth, and DBW were also compared between JENDL-4.0u and ENDF/B-VII.0. The comparison of the results shows JENDL-4.0u and ENDF/B-VII.0 give comparable results within the standard deviation, the difference between libraries is small. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Correlating Broadband Photoluminescence with Structural Dynamics in Layered Hybrid Halide Perovskites

The emission of white light from a single material is atypical and is of interest for solid-state lighting applications. Broadband light emission has been observed in some layered perovskite derivatives, A 2 PbBr 4 (A = R-NH 3 + ), and correlates with static structural distortions corresponding to out-of-plane tilting of the lead bromide octahedra. While materials with different organic cations can yield distinct out-of-plane tilts, the underlying origin of the octahedral tilting remains poorly understood. Using high energy resolution (e.g., quasi-elastic) neutron scattering, this contribution details the rotational dynamics of the organic cations in A 2 PbBr 4 materials where A = n-butylammonium (nBA), 1,8-diaminooctammonium (ODA), and 4-aminobutyric acid (GABA). The organic cation dynamics differentiate (nBA) 2 PbBr 4 from (ODA)PbBr 4 or (GABA) 2 PbBr 4 in that the larger spatial extent of dynamics of nBA yields a larger effective cation radius. The larger effective volume of the nBA cation in (nBA) 2 PbBr 4 yields a closer to ideal A-site geometry, preventing the out-of-plane tilt and broadband luminescence. In all three compounds, we observe hydrogen dynamics attributed to rotation of the ammonium headgroup and at a time scale faster than the white light photoluminescence studied by time-correlated single photon counting spectroscopy. This supports a previous assignment of the broadband emission as resulting from a single ensemble, such that the emissive excited state experiences many local structures faster than the emissive decay. Furthermore, the findings presented here highlight the role of the organic cation and its dynamics in hybrid organic–inorganic perovskites and white light emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Context-Aware Learning for Inverse Design in Photovoltaics

This document describes progress in the ARPA-E DIFFERENTIATE project titled “Context-Aware Learning for Inverse Design in Photovoltaics” during the period of May 2019 to May 2022. This project is being performed at Iowa State University, New York University, Stanford University, and National Renewable Energy Laboratory. The project aims to develop a new machine learning (ML) framework to significantly accelerate the design of organic microstructures for improved organic photovoltaic performance. In this project, we had developed an inverse design framework using Deep Learning called InvNets for generating microstructures with desired physics-driven properties. As a preliminary product, in Milestone 3, we demonstrated how InvNets show 20% improvement in the performance of the microstructures and over 100X speedup in the performance compared to traditional processes for physics-driven inverse design. Later, in Milestone 6, we demonstrated that InvNets work for more complex physics properties, specifically, generating microstructures for organic photovoltaic cells with desired current-voltage characteristics. Further, in Milestone 4, we explored the idea of using physics-aware surrogates for obtaining solutions of partial differential equations(PDE) called as DiffNets(now called as NeuFENets to avoid ambiguity of names). The connection between both frameworks is that DiffNet surrogates form the physics-aware surrogate in the InvNet framework. Finally in Milestone 8, we extend our framework for other physics domains. Specifically, we explore building geometry-aware NeuFENets by developing physics surrogates that exploit ideas from traditional immersed boundary finite element methods. With these updates, we are able to achieve all the Milestones.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Orbital-selective Dirac fermions and extremely flat bands in frustrated kagome-lattice metal CoSn

Layered kagome-lattice 3d transition metals are emerging as an exciting platform to explore the frustrated lattice geometry and quantum topology. However, the typical kagome electronic bands, characterized by sets of the Dirac-like band capped by a phase-destructive flat band, have not been clearly observed, and their orbital physics are even less well investigated. Here, we present close-to-textbook kagome bands with orbital differentiation physics in CoSn, which can be well described by a minimal tight-binding model with single-orbital hopping in Co kagome lattice. The capping flat bands with bandwidth less than 0.2 eV run through the whole Brillouin zone, especially the bandwidth of the flat band of out-of-plane orbitals is less than 0.02 eV along Γ-M. The energy gap induced by spin-orbit interaction at the Dirac cone of out-of-plane orbitals is much smaller than that of in-plane orbitals, suggesting orbital-selective character of the Dirac fermions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Detection, monitoring, and determination of location of changes in metallic structures using multimode acoustic signals

Acoustic transducers generate and receive acoustic signals at multiple locations along a surface of rigid structure, wherein longitudinal spacing between transducer locations define measurement zones. Acoustic signals with chosen amplitude-time-frequency characteristics excite multiple vibration modes in the structure within each zone. Small mechanical changes in inspection zones lead to scattering and attenuation of broadband acoustic signals, which are detectable as changes in received signal characteristics as part of a through-transmission technique. Additional use of short, narrowband pulse acoustic signals as part of a pulse-echo technique allows determination of the relative location of the mechanical change within each zone based on the differential delay profiles. For accurate acoustic modeling and simulation, the mesh size, time step, time delay, and time-window size are optimized. Frequency normalization of the Short-Time Fourier Transform of acoustic response output improves experiment-simulation cross-validation. Applications of the method to structures with arbitrarily complex geometries are also demonstrated.

Findikoglu, Alp Tugrul↗
DOE OSTI2020

BISON Robustness and Performance Improvements

BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Quantum Gravity and Laser Interferometry: Towards Observable Predictions

Understanding quantum gravity remains one of the deepest challenges in modern physics, as direct experimental access to Planck-scale effects is beyond current technological reach. However, recent theoretical advances indicate that quantum fluctuations of spacetime may produce measurable effects in precision experiments, particularly near causal horizons. This opens new avenues for testing quantum gravity phenomena through high-precision measurement techniques. This dissertation develops multiple theoretical models to characterize these effects and examines their potential observational signatures in future gravitational wave interferometers. We begin by investigating the role of quantum fluctuations in near-horizon geometries through the lens of the AdS/CFT correspondence, which provides a powerful framework for understanding the interplay between quantum field theory and general relativity via holographic principles. By modeling stochastic energy-momentum sources in Rindler-AdS spacetime, we demonstrate that vacuum fluctuations transform the Einstein equations into a Langevin-type stochastic differential equation, leading to potentially observable fluctuations in photon traversal times. Extending this approach to Minkowski spacetime, we establish a correspondence between gravitational shockwaves and fluid dynamics, showing that near-horizon perturbations satisfy an equation analogous to that governing incompressible fluids, thereby reinforcing the membrane paradigm and hydrodynamic analogies in the context of the fluid/gravity duality. Furthermore, we construct the covariant phase space of a spherically symmetric causal diamond in Minkowski spacetime, identifying two fundamental charges that govern its evolution. These results provide a foundation for quantizing causal horizons and understanding their microscopic degrees of freedom. Building upon these theoretical developments, we further examine a related stochastic phenomenon: the gravitational wave memory background arising from the cumulative memory steps produced by supermassive black hole mergers. After reviewing the standard stochastic gravitational wave background, gravitational memory effects, and BMS symmetries, we model the stochastic memory background using a Brownian motion framework. We show that while the cumulative memory background initially appears above the sensitivity curve of space-based interferometers like LISA, the realistic subtraction of individually resolvable merger events substantially suppresses the residual signal, making its detection more challenging. This highlights the critical importance of source subtraction when evaluating the detectability of gravitational memory effects. By bridging fundamental theory with experimental prospects, this dissertation contributes to the ongoing effort to uncover the quantum nature of spacetime through precision measurement techniques. Whether through detecting quantum spacetime fluctuations, gravitational memory backgrounds, or probing the symmetries of causal horizons, the pursuit of observable quantum gravity phenomena continues to expand the frontiers of both theory and experiment.

Zhang, Yiwen [Caltech] (ORCID:0000000323559416)↗
DOE OSTI2020

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗
DOE OSTI2025

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗
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↗