Skip to article frontmatterSkip to article content

Function

The learned function, fθ\boldsymbol{f_\theta}, will map the spatial coordinates, xϕRDϕ\mathbf{x}_\phi \in \mathbb{R}^{D\phi}, and time coordinate, tRt \in \mathbb{R}, to sea surface height, uRu \in \mathbb{R}.

u=fθ(xϕ,t)u = \boldsymbol{f_\theta}(\mathbf{x}_\phi, t)

Loss

The standard loss term is data-driven

Ldata=MSE(u,u^)=1Nn=1N(ufϕ(xϕ,t))2\mathcal{L}_{data} = \text{MSE}(u, \hat{u}) = \frac{1}{N} \sum_{n=1}^N \left(u - \boldsymbol{f_\phi}(\mathbf{x}_\phi, t) \right)^2

However, there is no penalization to make the field behave the way we would expect. We also want a regularization which makes the field, uu, behave how we would expect. This can be achieved by adding a physics-informed loss regularization term to the total loss.

L=Ldata+λLphy\mathcal{L} = \mathcal{L}_{data} + \lambda \mathcal{L}_{phy}

This loss term can be minimized by effectively minimizing a PDE function. For example:

fphy(x,t):=tu(x,t)+N[u(x,t)]=0\boldsymbol{f}_{phy}(\mathbf{x},t):= \partial_t u(\mathbf{x},t) + \mathcal{N}[u(\mathbf{x},t)] = 0

where t\partial_t is the derivative of the field, uu, wrt to time and N[]\mathcal{N}[\cdot] are some partial differential equations. We are interested in minimizing the full PDE, which we denote fphy\boldsymbol{f}_{phy}, st it is 0. So the standard loss function applies.


Examples

Below are some examples of how I have used the PINNs loss/regularization formulation in my own research.


QG Equation

We have the following PDE for the QG dynamics:

tq+detJ(ψ,q)=0\partial_t q + \det J(\psi, q) = 0

where q(x,t)R2×RRq(x,t) \in \mathbb{R}^2 \times \mathbb{R} \rightarrow \mathbb{R} is the potential vorticity (PV), ψ(x,t)R2×RR\psi(x,t) \in \mathbb{R}^2 \times \mathbb{R} \rightarrow \mathbb{R} is the stream function, t\partial_t is the partial derivative wrt tt, J\boldsymbol{J}, is the Jacobian operator and detJ(,)\det \boldsymbol{J}(\cdot,\cdot) is the determinant of the Jacobian.

Objective: We want to convert this PDE in terms of sea surface height (SSH) instead of PV and the stream function.


QG Equation 4 SSH (TLDR)

Note: For the ease of notation, let’s denote uu as the SSH. The above PDE can be written in terms of uu

t2u+c2tu+c1detJ(u,2u)+c3detJ(u,u)=0\partial_t \nabla^2 u + c_2 \partial_t u + c_1\det \boldsymbol{J}(u, \nabla^2 u) + c_3 \det \boldsymbol{J}(u, u) = 0

See the following page for more and how this equation was derived.

Example I: Divergence + Curl Free

Example II: Quasi-Geostrophic Equations