---
title: 'PI-GANO: Physics-Informed Neural Operator'
url: https://www.emergentmind.com/topics/pi-gano
type: topic
---

# PI-GANO: Physics-Informed Neural Operator

PI-GANO (Physics-Informed Geometry-Aware Neural Operator) is a mesh-free neural operator architecture designed to model steady-state flows simultaneously through and around porous structures of arbitrary geometry, boundary condition, and material parameter, without retraining for each new configuration. It systematically integrates the governing physics of the Navier–Stokes and Darcy–Forchheimer equations within a unified loss formulation, supports variable boundary conditions and material properties, and achieves rapid, geometry-agnostic inference suitable for design acceleration in engineering and environmental flow scenarios [2602.14108].

## 1. Problem Definition and Motivation

Many problems in engineering, environmental, or atmospheric science feature bodies with both permeable (porous) interiors and free-fluid exteriors—for example, windbreaks, catalysts, or submerged gabions. Standard computational fluid dynamics (CFD) utilizes the Navier–Stokes equations in fluid regions (Ωₚ) and the Darcy–Forchheimer extension in porous regions (Ω_f), with appropriate coupling on the fluid–porous interface (Γ). Classical approaches demand meshing bespoke to each geometry and repeated PDE solves whenever boundary conditions or material parameters (e.g., Darcy D, Forchheimer F, viscosity μ) change.

PI-GANO circumvents these bottlenecks by learning the mapping
$${ \{\partial \Omega_p,\ D,\ F,\ \mu,\ \mathrm{BC}\}\ \longrightarrow\ \{\mathbf{u}(x),\ p(x)\} }$$
in a single, end-to-end physics-informed neural operator. This architecture generalizes to unseen geometries, boundary conditions, and parameter settings without retraining [2602.14108].

## 2. Mathematical and Physical Model

### 2.1. Governing Equations

PI-GANO partitions the domain Ω into:
- $\Omega_p$: porous region
- $\Omega_f$: free-fluid region
- Γ = ∂Ωₚ: interface
- ∂Ωₚ, ∂Ω_f: exterior boundaries

On $\Omega_f$, it enforces the incompressible steady Navier–Stokes equations:
\[
\nabla\cdot \mathbf{u} = 0, \quad \rho\,(\mathbf{u}\cdot\nabla)\mathbf{u} = -\nabla p + \mu\,\nabla^2\mathbf{u}
\]
On $\Omega_p$, the Darcy–Forchheimer extension applies:
\[
\rho(\mathbf{u}\cdot\nabla)\mathbf{u} = -\nabla p + \mu \nabla^2 \mathbf{u} - (\mu D + \tfrac{1}{2} \rho F |\mathbf{u}|)\mathbf{u}
\]
Here, $D$ and $F$ are calculated from the porosity $\phi$ and particle diameter $d_p$ via Kozeny–Carman relations:
\[
D = \frac{180}{d_p^2}\frac{(1-\phi)^2}{\phi^3}, \qquad F = \frac{1.8}{d_p}\frac{1-\phi}{\phi^3}
\]
Boundary and interface conditions require continuity of velocity and normal stress at the interface Γ, and prescribed Dirichlet/Neumann conditions at the domain boundary ($\partial \Omega_f$).

### 2.2. Unified Physics-Informed Loss

Collocation points $x_i$ throughout Ω are annotated with an indicator function
\[
\chi_p(x) = \begin{cases} 1, & x \in \Omega_p \\ 0, & x \in \Omega_f \end{cases}
\]
At each point, the residual for the full steady PDE is:
\[
r_m(x_i) = \rho (\mathbf{u}^\theta_i \cdot \nabla) \mathbf{u}^\theta_i + \nabla p^\theta_i - \mu \nabla^2 \mathbf{u}^\theta_i + \chi_p(x_i) (\mu D + \tfrac12 \rho F |\mathbf{u}^\theta_i|)\mathbf{u}^\theta_i
\]
The incompressibility residual is:
\[
r_c(x_i) = \nabla \cdot \mathbf{u}^\theta_i
\]
The total loss minimized by PI-GANO is:
\[
\mathcal{L}(\theta) = \lambda_m\,\mathcal{L}_m + \lambda_c\,\mathcal{L}_c + \lambda_b\,\mathcal{L}_b + \lambda_d\,\mathcal{L}_d
\]
with
\[
\begin{aligned}
&\mathcal{L}_m = \frac{1}{N_f}\sum_{i=1}^{N_f}\|r_m(x_i)\|^2, \qquad \mathcal{L}_c = \frac{1}{N_f}\sum_{i=1}^{N_f}\|r_c(x_i)\|^2 \\
&\mathcal{L}_b = \frac{1}{N_b}\sum_{j=1}^{N_b} \left(\|\mathbf{u}^\theta(x_j)-\mathbf{u}(x_j)\|^2 + \|p^\theta(x_j)-p(x_j)\|^2\right) \\
&\mathcal{L}_d = \frac{1}{N_d}\sum_{k=1}^{N_d} \left(\|\mathbf{u}^\theta(x_k)-\mathbf{u}(x_k)\|^2 + \|p^\theta(x_k)-p(x_k)\|^2\right)
\end{aligned}
\]
Here, $\lambda$ hyperparameters control the enforcement strength, and $N$ denotes the number of various collocation/data points [2602.14108].

## 3. Architecture and Training

### 3.1. Input Encoding and Latent Representations

PI-GANO extends the PointNet backbone of PIPN, forming a full neural operator through:
- **Per-point features:** coordinates $x \in \mathbb{R}^d$, signed-distance to interface $\mathrm{SDF}(x)$, porous-region indicator $\chi_p(x)$, and boundary-type vector $I_b(x)$ (e.g., inlet, outlet, wall).
- **Global parameter inputs:** inlet velocity $U$, Darcy $D$, Forchheimer $F$ (with absent BCs set to zero).

A shared MLP encoder $\phi_e$ produces a per-point embedding, with max-pooling aggregation yielding a global geometry latent $g\in\mathbb{R}^{512}$. Boundary/parameter information is embedded by a second MLP $\psi_e$ acting on boundary samples, generating a latent $b\in\mathbb{R}^{d_b}$.

For each query (or collocation) point $x_i$, the triple $[$ per-point encoding, $g$, $b$ $]$ is fed to a shared decoder MLP $\rho$ which outputs $(u_x, u_y, u_z)^\theta(x_i)$ and $p^\theta(x_i)$. This shared decoder design reduces parameter count and achieves better field consistency than separate output branches [2602.14108].

### 3.2. Training Protocol

- **Optimizer:** Adam with initial learning rate $10^{-3}$; exponential decay factor $\alpha$ set per task.
- **Activation:** SiLU; $\tanh$ for manufactured solutions only.
- **Dropout:** $p=0.05$ on last two layers (disabled for manufactured solution).
- **Epochs:** $3\,000$.
- **Loss weights:** see Table 1.

| Experiment | $\lambda_b$ | $\lambda_m=\lambda_c$ | $\lambda_d$ | $\alpha$ (lr–decay) | #int | #bnd | #obs |
|------------|-------------|-----------------------|-------------|---------------------|------|------|------|
| MMS        |      1      |          1            |    —        |    0.9995           |  667 | 168  |  —   |
| Duct (fixed) |    1      |          1            |   100       |    0.9990           | 1000 | 200  | 500  |
| Windbreak (3D) | 1      |         10            |    1        |    0.9990           | 2000 |1000  |1000  |

*Explanation of terms: MMS—manufactured solutions; Duct—2D; Windbreak—3D scenarios. #int/#bnd/#obs: counts of interior, boundary, and observation points.*

## 4. Dataset Generation and Test Cases

### 4.1. CFD Dataset Creation

Simulation data is generated with OpenFOAM's `simpleFoam` solver. Domains include both 2D ducts with porous obstacles and 3D windbreak environments (trees, buildings). Boundary conditions and meshes are precisely matched to the requirements of steady, laminar, incompressible, Newtonian flow, with accurate convergence.

- **2D shapes:** star, octagon, ellipse; composite geometries held out for generalization tests.
- **3D:** Extracted from 3D tree species (oak, pine, cypress, willow, acacia, eucalyptus) and house models, with porosity estimated by projected leaf count.

Features per-point include velocity $\mathbf{u}$, pressure $p$, region indicator $\chi_p$, porosity $\phi$, SDF, and boundary-type. Uniform sampling with increased density near the interface ensures balanced loss enforcement.

### 4.2. Data Preprocessing

Point features are standardized (binary fields to [0,1], $(\mathbf{u},p)$ via $Z$-score). PDE coefficients are scaled accordingly.

- **Splits:** $60\%$ train, $20\%$ validation, $20\%$ test for all cases.

## 5. Performance and Generalization Results

### 5.1. Method of Manufactured Solutions

Verification with manufactured 2D fields,
\[
u_x = \cos x \sin y,\quad u_y = -\cos y \sin x,\quad p \equiv 0
\]
demonstrates PI-GANO can reproduce known analytical solutions to within MAEs of order $10^{-2}$ on both training and unseen geometries.

### 5.2. 2D Fixed-BC, Variable-BC, and 3D Cases

**Case 2 (2D, fixed-BC duct, PIPN):**
- Test MAE: $u_x\approx 3.1\times 10^{-3}$, $u_y\approx 1.2\times 10^{-3}$, $p\approx 9\times 10^{-4}$
- Inference: $0.02$ s per case, OpenFOAM baseline $1.17$ s

**Case 3 (2D, variable-BC, PI-GANO):**
- Test MAE: $u_x\approx 1.95\times 10^{-3}$, $u_y\approx 8.5\times 10^{-4}$, $p\approx 4.0\times 10^{-4}$
- Unseen BC/geometry: $u_x\approx 1.57\times 10^{-2}$, $u_y\approx 5.12\times 10^{-3}$, $p\approx 6.86\times 10^{-3}$
- Inference: $0.010$ s, OpenFOAM baseline $2.08$ s

**Case 4 (3D windbreak, PI-GANO):**
- Global MAE: $u_x\approx 1.15\times 10^{-7}$, $u_y\approx 4.14\times 10^{-8}$, $u_z\approx 3.18\times 10^{-8}$, $p\approx 2.39\times 10^{-12}$
- Inference: $0.014$ s, OpenFOAM baseline $20$ s

In all cases, maximum errors concentrate near sharp corners, region interfaces, and steep-gradient zones; the streamwise velocity $u_x$ is particularly affected.

## 6. Comparative Analysis and Limitations

### 6.1. Comparative Assessment

PI-GANO, in contrast to PINNs and PIPN:
- Handles novel geometry, boundary condition, and material-parameter combinations in a single forward pass, by virtue of its operator structure.
- Achieves test and generalization errors below those of PIPN when boundary/parameter variations are present.
- Reduces inference latency by one to three orders of magnitude relative to full CFD simulations.
- Provides mesh-free prediction and does not require re-meshing or solver reruns for new cases.

### 6.2. Limitations and Prospective Advancements

Dominant sources of error remain at flow–porous interfaces and in regions with steep gradients. Remedies discussed include:
- Employing PointNet++ for enhanced capture of fine-scale geometry,
- Adoption of spectral (FNO-style) kernels,
- Adaptive collocation-point sampling targeting high-error zones,
- Incorporating additional PDE physics, such as $k$–$\epsilon$ turbulence modeling and automatic multi-objective loss weighting.

These enhancements are suggested to further improve the robustness and accuracy of the PI-GANO operator [2602.14108].

## 7. Conclusion

PI-GANO establishes a systematic, physics-informed, mesh-free neural operator framework for steady-state flows through and around porous geometries. It enables real-time, accurate inference across unseen spatial configurations, boundary settings, and material properties, potentially accelerating parametric and geometric design processes without repeated training or meshing. It unifies geometric awareness, PDE constraints, and operator learning within a single model and demonstrates generalization performance and computational speed superior to conventional CFD approaches and prior physics-informed neural architectures [2602.14108].

Source: https://www.emergentmind.com/topics/pi-gano