Search NASA⌕ Search

Engineering topics

Roy, Pratanu

Publications and source records attributed to Roy, Pratanu.

Exact enforcement of temporal continuity in sequential physics-informed neural networks

The use of deep learning methods in scientific computing represents a potential paradigm shift in engineering problem solving. One of the most prominent developments is Physics-Informed Neural Networks (PINNs), in which neural networks are trained to satisfy partial differential equations (PDEs). While this method shows promise, the standard version has been shown to struggle in accurately predicting the dynamic behavior of time-dependent problems. To address this challenge, methods have been proposed that decompose the time domain into multiple segments, employing a distinct neural network in each segment and directly incorporating continuity between them in the loss function of the minimization problem. In this work we introduce a method to exactly enforce continuity between successive time segments via a solution ansatz. This hard constrained sequential PINN (HCS-PINN) method is simple to implement and eliminates the need for any loss terms associated with temporal continuity. The method is tested for a number of benchmark problems involving both linear and non-linear PDEs. Examples include various first order time dependent problems in which traditional PINNs struggle, namely advection, Allen–Cahn, and Korteweg–de Vries equations. Furthermore, second and third order time-dependent problems are demonstrated via wave and Jerky dynamics examples, respectively. Notably, the Jerky dynamics problem is chaotic, making the problem especially sensitive to temporal accuracy. Finally, the numerical experiments conducted with the proposed method demonstrated superior convergence and accuracy over both traditional PINNs and the soft-constrained counterparts.

42 ENGINEERING↗

Interface PINNs (I-PINNs): A physics-informed neural networks framework for interface problems

Here, we present a novel physics-informed neural networks (PINNs) framework for modeling interface problems, termed Interface PINNs (I-PINNs). I-PINNs uses different neural networks for any two subdomains separated by a sharp interface such that the neural networks differ only through their activation functions while the other parameters remain identical. The performance of I-PINNs, conventional PINNs, and other existing domain-decomposition PINNs methods such as extended PINNs (XPINNs) and multi-domain PINN (M-PINN) is compared through several one-dimensional, two-dimensional, and three-dimensional benchmark elliptic interface problems. The results demonstrate that I-PINNs provides a root-mean-square-error accuracy, at least two orders of magnitude better than conventional PINNs and XPINNs at approximately one-tenth of the computational cost of conventional PINNs and half the cost of XPINNs. Additionally, while I-PINNs and M-PINN provide comparable accuracies, M-PINN is found to be approximately 50% more expensive.

42 ENGINEERING↗

PhyML

Python scripts for physics informed machine learning for time-dependent problems. The code implements hard constrained sequential physics-informed neural networks (HCSPINNs) using JAX library. The details of the method and implementation can be found in the following paper: Roy, P., & Castonguay, S. (2024). Exact Enforcement of Temporal Continuity in Sequential Physics- Informed Neural Networks. arXiv preprint arXiv:2403.03223. (https://arxiv.org/abs/2403.03223)

Roy, Pratanu↗

Train small, model big: Scalable physics simulators via reduced order modeling and domain decomposition

Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the pilot-scale model, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. The ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. Here this method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about 15–40 times faster with only ~1% relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DynamicSimCCS

This code solves for the one-dimensional concentration and temperature distribution within a carbon capture column using a homogenized absorber model. These solutions are compared to experimental data.

Roy, Pratanu↗

A chemo-mechanical model for describing sorption hysteresis in a glassy polyurethane

Abstract Hysteretic sorption and desorption of water is observed from 0 to 95% relative humidity and 298–333 K on a glassy polyurethane foam. It is postulated that sorption-induced swelling of the glassy polyurethane increases the concentration of accessible hydrogen-bonding adsorption sites for water. The accessibility of sites is kinetically controlled due to the restricted thermal motions of chains in the glassy polymer, causing a difference in accessible site concentrations during sorption and desorption. This discrepancy leads to hysteresis in the sorbed concentrations of water. A coupled chemo-mechanical model relating volumetric strain, adsorption site concentration, and sorbed water concentration is employed to describe water sorption hysteresis in the glassy polyurethane. This model not only describes the final mass uptake for each relative humidity step, but also captures the dynamics of water uptake, which exhibit diffusion and relaxation rate-controlled regimes.

36 MATERIALS SCIENCE↗

Systems and methods for periodic nodal surface based reactors, distributors, contactors and heat exchangers

The present disclosure relates to a transport mechanism apparatus for transporting at least one of a gas or a fluid. The transport mechanism may have an inlet, an outlet and an engineered cellular structure forming a periodic nodal surface, which may include a triply periodic minimal surface (TPMS) structure. The structure is formed in a layer-by-layer three dimensional (3D) printing operation to include cells propagating in three dimensions, where the cells include non-intersecting, continuously curving wall portions having openings, and where the opening in the cells form a plurality of flow paths throughout the transport mechanism from the inlet to the outlet, and where portions of the cells form the inlet and the outlet.

Iyer, Jaisree Kannan↗

Water sorption equilibrium and dynamics on Epon 828-Epikure 3125

Epon 828-Epikure 3125 absorbs a significant amount of water (up to 4 weight percent), but diffusion of water into the material is slow. A one cm thick open slab may take 80 days to reach equilibrium at 70 °C and five years at 20 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable Reduced Order Model with Discontinuous Galerkin Domain Decomposition

scaleupROM is a scalable, physics-constrained reduced order model (ROM). It aims to provide robust, accelerated physics predictions at extrapolated scales, based on the small, component-level data. This is implemented by combining projection-based ROM with discontinuous Galerkin domain decomposition, in the framework of MFEM and libROM. It currently supports the Poisson equation and Stokes flow equation, and more work is in progress toward general, nonlinear physics systems.

Chung, Seung Whan↗

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hierarchical triply periodic minimal surface structures as heat exchangers and reactors

The present disclosure relates to a transport mechanism apparatus for transporting at least one of a gas or a fluid. The transport mechanism may have an inlet, an outlet and a triply periodic minimal surface (TPMS) structure. The TPMS structure is formed in a layer-by-layer three dimensional (3D) printing operation to include cells propagating in three dimensions, where the cells include wall portions having openings, and where the cells form a plurality of flow paths throughout the transport mechanism from the inlet to the outlet, and where the cells form the inlet and the outlet.

Roy, Pratanu↗

Hydrodynamics of countercurrent flow in an additive-manufactured column with triply periodic minimal surfaces for carbon dioxide capture

Computational investigation of the countercurrent flows in an additively manufactured (AM) column packed with triply periodic minimal surface (TPMS) is challenging. Multiphase flow simulations using the volume of fluid (VOF) method can explain the column’s hydrodynamic characteristics, including interfacial area, liquid holdup, pressure drop, flow regime, etc. The interfacial area is a key factor dictating the mass transfer via the absorption process, and thereby the overall efficiency of the column. The impacts of a solvent's physical properties, liquid and gas loads, and contact angle on the hydrodynamics are extensively explored. The solvents include traditional solvent (monoethanolamine [MEA]) used in industry as well as a water-lean solvent (EEMPA) as prospective solvents for carbon capture. As expected, interfacial area and liquid holdup increase with increasing liquid loads (q L ). The EEMPA exhibits higher values for the interfacial area and liquid holdup compared to those for the MEA. The dry and wet pressure drops per unit length in the column at different gas loads (F G ) are calculated and further compared with the empirical correlations for traditional packings. At low gas loads, i.e., in the preloading regime, the predicted values of the wet and dry pressure drops match well with the corresponding values calculated from available correlations. The gas load has a marginal impact on the interfacial area at its lower value, but the higher value of F G (>2.18 $\sqrt{Pa}$) causes flooding in the column. The effect of varying contact angles on the interfacial area and liquid holdup in the preloading regime is also studied. Both interfacial area and liquid holdup decrease with the increase in the value of the contact angle. Altogether, the effect of the contact angle on interfacial area is more pronounced than other parameters influencing the interfacial area.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced absorber heat integration via heat exchange packings

A rate-based model of an absorption column was developed and used to analyze several intercooling strategies utilizing “heat exchange packings.” These packings are capable of removing heat from the column and transferring it to a cooling fluid within the packing. For absorption of CO 2 into aqueous monoethanolamine under industrial conditions, intercooling via heat exchange packings placed along 10–20% of the column could reduce the column height by ~15%. The height of these columns was close to the minimum theoretical value, calculated by numerically optimizing the temperature profile. Effective intercooling could also be achieved by using the cool, rich solvent as the cooling fluid. This reduces the cooling load and facilitates recovery of waste heat. Heat exchange packings could also be used to redistribute heat within the column, reducing the column height with no net cooling load. However, this approach requires larger heat transfer coefficients than have been experimentally observed.

20 FOSSIL-FUELED POWER PLANTS↗