---
title: Semidiscrete Inverse Source Problem
url: https://www.emergentmind.com/topics/semidiscrete-inverse-source-problem
type: topic
---

# Semidiscrete Inverse Source Problem

A semidiscrete inverse source problem is an inverse problem in which the source is reconstructed under a mixed continuous–discrete representation: one variable is discretized while the remaining variables remain continuous, or a continuous model is reduced to a finite-dimensional inverse system through a semidiscretization procedure. In the recent literature, this term covers several non-equivalent but structurally related settings: time-discretized reconstruction of a time-dependent source coefficient in an evolution equation, piecewise constant-in-time source identification in parabolic dynamics, space-semidiscrete stochastic parabolic inversion, and frequency-semidiscrete reformulations of multi-frequency Helmholtz source recovery [2502.04821], [2111.02285], [2509.03760], [1901.10047]. A central theme is that semidiscretization converts the inverse problem into a coupled system of algebraic, elliptic, or finite-dimensional subproblems while preserving enough analytic structure to establish existence, uniqueness, stability, or convergence.

## 1. Scope and defining formulations

The phrase “semidiscrete inverse source problem” is used in at least three precise senses in the cited literature.

First, it can mean a **time-discrete inverse source problem** for a continuous-in-space PDE. In the semilinear pseudo-parabolic setting, time is partitioned by Rothe’s method, derivatives are replaced by backward differences, and at each time level one reconstructs a discrete source coefficient together with the state [2502.04821].

Second, it can mean that the **unknown source itself is semi-discrete**, typically piecewise constant in time and continuous in space. For the linear parabolic equation,
\[
F(x,t)=\sum_{k=1}^K p_k(x)\,\chi_{t\in[t_{k-1},t_k)},
\]
the unknown consists of the switching times and the spatial source states \(p_k\in L^2(\Omega)\), while the governing PDE remains continuous [2111.02285].

Third, it can mean a **space semidiscrete inverse problem**, where the underlying stochastic parabolic equation is discretized in space but remains continuous in time. In that setting, the unknown is a random source term in the diffusion part of the semidiscrete SPDE, identified from terminal data and partial boundary flux traces [2509.03760].

A related usage appears in inverse potential and Helmholtz problems. There the forward model is continuous, but the inverse reconstruction is carried out over a finite-dimensional source space or a truncated frequency basis. The potential-theoretic formulation approximates a continuous boundary source by piecewise constant coefficients on boundary segments and solves a Nonnegative Least Squares problem [2510.04693]. The multi-frequency Helmholtz formulation expands the \(k\)-dependence in a finite orthonormal basis and reduces the inverse source problem to a coupled elliptic system for finitely many modal coefficients [1901.10047]. This suggests that “semidiscrete” is best understood structurally rather than purely by whether time or space is discretized.

## 2. Canonical PDE models and observation operators

A representative time-discrete model is the semilinear pseudo-parabolic equation on a bounded Lipschitz domain \(\Omega\subset\mathbb{R}^d\),
\[
\left\{
\begin{array}{rl}
\partial_t u(t,\mathbf{x}) - \nabla \cdot\left(\eta(t,\mathbf{x}) \nabla \partial_t u(t,\mathbf{x})\right) - \nabla \cdot \left( \kappa(t,\mathbf{x}) \nabla u (t,\mathbf{x})\right) & = f(u(t,\mathbf{x})) + p(t,\mathbf{x})\,h(t),\\
\kappa(t,\mathbf{x}) \nabla u(t,\mathbf{x}) \cdot \mathbf{n}(\mathbf{x}) = g(t,\mathbf{x}),\\
u(0,\mathbf{x}) = \tilde{u}_0(\mathbf{x}),
\end{array}
\right.
\]
where the source is separated into a known spatial factor \(p(t,\mathbf{x})\) and an unknown scalar coefficient \(h(t)\) [2502.04821]. The inverse data are nonlocal:
\[
\int_\Omega u(t,\mathbf{x})\,d\mathbf{x} = m(t),\quad t>0.
\]
Integrating the PDE over \(\Omega\) and using the Neumann boundary condition yields an explicit reconstruction formula for \(h(t)\), provided
\[
\omega(t):=\int_\Omega p(t,\mathbf{x})\,d\mathbf{x}\neq 0,
\quad
\omega_0:=\min_{t\in[0,T]}\omega(t)>0,
\]
so identifiability is tied directly to a nonvanishing spatial average of the source shape [2502.04821].

In the parabolic semi-discrete source model, the forward equation is
\[
(\partial_t + A)u(x,t) = \sum_{k=1}^K p_k(x)\,\chi_{t\in[t_{k-1},t_k)},
\]
with homogeneous Dirichlet boundary condition and zero initial condition. The measurements are sparse boundary fluxes
\[
\partial_{\nu_\kappa}u(x,t)\quad\text{for }(x,t)\in\Gamma\times(0,\infty),
\]
where \(\Gamma\subset\partial\Omega\) is an arbitrary nonempty open subset [2111.02285]. Here the discrete structure lies in time, but the profiles \(p_k\) remain continuous in space.

In the space-semidiscrete stochastic parabolic problem, the unknown is the random source term \(g\) in
\[
dw - \sum_{i=1}^n D_i(\gamma_i D_i w)\,dt
=
\Big(\sum_{i=1}^n a_{1i} A_i D_i(w) + a_2 w\Big)\,dt + g\,dB(t)
\]
on the spatial mesh \(\mathcal M\), with homogeneous discrete Dirichlet boundary condition [2509.03760]. The observation operator is
\[
\Lambda_1(g)
=
\big(\partial_\nu w|_{\Gamma^+\times(0,T)},\, w(T)\big),
\]
combining terminal data with the trace of the discrete normal derivative on a subset of the discrete boundary [2509.03760].

Other observation models illustrate further semidiscrete patterns. The inverse potential problem uses values of the potential on an observation curve \(\Gamma\) and reconstructs a nonnegative boundary density on a candidate enclosing boundary \(S=\partial\Omega\) [2510.04693]. The multi-frequency inverse source problem for the Helmholtz equation assumes boundary data for all wavenumbers \(k\in(\underline k,\overline k)\) and projects the frequency dependence onto a finite orthonormal basis \(\{\Psi_m\}_{m=1}^N\), producing finitely many spatial coefficient functions \(v_m(\mathbf x)\) [1901.10047].

## 3. Time semidiscretization, Rothe’s method, and stepwise reconstruction

The most explicit semidiscrete construction in the cited corpus is the Rothe scheme for the pseudo-parabolic inverse source problem. Let \(n\in\mathbb N\), \(\tau=T/n\), \(t_i=i\tau\), and
\[
\delta z_i:=\frac{z_i-z_{i-1}}{\tau}.
\]
The time-discrete inverse problem at level \(t_i\) is: find \(u_i\in H^1(\Omega)\) and \(h_i\in\mathbb R\) such that, for all \(\varphi\in H^1(\Omega)\),
\[
\langle \delta u_i,\varphi\rangle
+ \langle \eta_i\nabla \delta u_i,\nabla \varphi\rangle
+ \langle \kappa_i \nabla u_i,\nabla \varphi\rangle
=
h_i \langle p_i,\varphi\rangle
+ \langle f(u_{i-1}),\varphi\rangle
+ \langle \eta_i(\partial_t G)_i,\varphi\rangle_{\partial\Omega}
+ \langle g_i,\varphi\rangle_{\partial\Omega},
\]
with
\[
h_i =\frac{1}{\omega_i}\Big[ (m^\prime)_i - \int_{\partial\Omega} \eta_i(\partial_t G)_i\,d\gamma - \int_{\partial\Omega} g_i\,d\gamma - \int_\Omega f(u_{i-1})\,dx \Big],
\qquad
u_0=\tilde u_0.
\]
This produces a two-stage update: first reconstruct \(h_i\) explicitly from the previous state \(u_{i-1}\) and measured data, then solve a linear elliptic variational problem for \(u_i\) [2502.04821].

The analysis is stepwise. Once \(u_{i-1}\in H^1(\Omega)\) is known, \(h_i\) is explicit, and the bilinear form
\[
a(u,\varphi)
:=
\frac{1}{\tau}\langle u,\varphi\rangle
+
\frac{1}{\tau}\langle \eta_i\nabla u,\nabla\varphi\rangle
+
\langle \kappa_i\nabla u,\nabla\varphi\rangle
\]
is coercive:
\[
a(u,u)\ge \min\left\{\frac{1}{\tau},\,\frac{\eta_0}{\tau}+\kappa_0\right\}\|u\|_{H^1(\Omega)}^2.
\]
Hence the Lax–Milgram lemma yields a unique \(u_i\in H^1(\Omega)\), and induction gives a unique sequence \(\{(u_i,h_i)\}_{i=0}^n\) [2502.04821].

The same paper derives the discrete energy estimate
\[
\max_{1\le j\le n}\|u_j\|_{H^1(\Omega)}^2
+
\sum_{i=1}^n\|\delta u_i\|_{H^1(\Omega)}^2\tau
+
\sum_{i=1}^n\|u_i-u_{i-1}\|_{H^1(\Omega)}^2
+
\sum_{i=1}^n |h_i|^2\tau
\le C
\]
for \(\tau<\tau_0\), and then introduces piecewise linear and piecewise constant Rothe interpolants \(U_n,\overline U_n,\overline h_n\). Compactness and convergence arguments show that a subsequence converges to the unique weak solution of the continuous inverse problem [2502.04821].

A different notion of time semidiscretization appears in the parabolic source-sequence model. There the source is itself piecewise constant in time, so the unknown is not generated by discretizing a continuous source coefficient but prescribed as
\[
F(x,t)=\sum_{k=1}^K p_k(x)\chi_{t\in[t_{k-1},t_k)}.
\]
This structure enables interval-by-interval spectral analysis and a sequential Bayesian interpretation in which the spatial states \(p_k\) are estimated state by state [2111.02285]. The distinction is significant: in the Rothe formulation, the semidiscrete problem is an approximation method for a continuous inverse problem; in the parabolic sequence model, the semidiscrete structure is part of the model itself.

## 4. Spatial semidiscretization and finite-dimensional source representations

Space-semidiscrete inverse source problems replace differential operators by discrete operators while retaining continuous time. In the stochastic parabolic setting, the discrete difference operators \(D_i\), averaging operators \(A_i\), and discrete trace operators \(t_r^i\) are defined on a mesh \(\mathcal M\subset G\cap h\mathbb Z^n\), and the semidiscrete operator is
\[
\mathcal A_h w := dw - \sum_{i=1}^n D_i(\gamma_i D_i w)\,dt.
\]
The inverse problem is to recover the adapted random source \(g\in L^2_{\mathcal F}(0,T;H_h^1(\mathcal M))\) from the observation operator \(\Lambda_1(g)\) [2509.03760].

The principal analytic tool is a global Carleman estimate for the semidiscrete stochastic parabolic operator. For homogeneous boundary data, the estimate controls weighted norms of \(w\), discrete first and second derivatives, and the stochastic term \(g\), by the interior source terms, boundary flux terms on \(\Gamma^+\), and the terminal term at \(t=T\) [2509.03760]. Applied to the difference of two solutions, it yields the Lipschitz estimate
\[
\|g_1-g_2\|_{L^2_{\mathcal F}(0,T;L_h^2(\mathcal M))}
\le
C \|\Lambda_1(g_1)-\Lambda_1(g_2)\|_{X_1},
\]
with \(C\) independent of \(h\), under the structural condition
\[
|D_i(\tilde g)| \le C |A_i(\tilde g)|
\quad\text{in }Q_i^*,\ \mathbb P\text{-a.s.}
\]
for \(\tilde g=g_1-g_2\) [2509.03760]. This is a genuinely semidiscrete stability theorem rather than a continuous theorem subsequently discretized.

A different finite-dimensional source representation appears in inverse potential theory. There, a continuous volume potential is approximated by a single-layer potential on a candidate boundary \(S\), and \(S\) is partitioned into \(N\) boundary segments \(S_j\). Approximating the single-layer density by piecewise constants \(v_j\) yields
\[
A v = f,
\qquad
a_{ij}=|S_j|\,G(\mathbf x_i,\mathbf y_j),
\]
and the inverse problem becomes the discrete Nonnegative Least Squares problem
\[
\|A v - f\|_2^2 \rightarrow \min_{v\ge 0}.
\]
The paper explicitly describes this as a semidiscrete inverse source problem in which the forward model is continuous but the inverse reconstruction is performed over a finite-dimensional source space [2510.04693].

The multi-frequency Helmholtz formulation is semidiscrete in the auxiliary frequency variable. Writing
\[
v(\mathbf x,k)\approx \sum_{m=1}^N v_m(\mathbf x)\Psi_m(k)
\]
and projecting the differentiated Helmholtz equation onto the basis \(\{\Psi_m\}_{m=1}^N\) leads to the coupled elliptic system
\[
D_N\,L(V(\mathbf x)) + S_N\,\mathbf n^2(\mathbf x)\,V(\mathbf x)=0,
\qquad
V=(v_1,\dots,v_N)^T.
\]
The source term is eliminated before spatial discretization, and the inverse source problem is recast as an overdetermined boundary value problem for \(V\) [1901.10047]. This suggests that semidiscretization can act not only on physical coordinates but also on spectral or parameter variables.

## 5. Identifiability, uniqueness, and regularization mechanisms

Identifiability in semidiscrete inverse source problems is tied to the chosen representation.

For the pseudo-parabolic problem, uniqueness follows from the representation of \(h(t)\) through the measurement \(m(t)\) and the condition \(\omega(t)\neq0\). Theorem 3.1 states that there exists at most one couple \(\{u,h\}\) such that
\[
h\in L^2(0,T),\quad
u\in C([0,T];H^1(\Omega)),\quad
\partial_t u\in L^2(0,T;H^1(\Omega)),
\]
and Theorem 4.4 establishes existence and uniqueness of the weak solution pair \((u,h)\) under the stated assumptions on \(\eta,\kappa,f,\tilde u_0,G,p,m\) [2502.04821].

For the parabolic semi-discrete source model, the main theorem is an identifiability result from sparse boundary measurements: the flux data on any nonempty open subset \(\Gamma\subset\partial\Omega\) uniquely determine the semi-discrete source
\[
\sum_{k=1}^K p_k(x)\chi_{t\in[t_{k-1},t_k)}.
\]
More precisely, equality of the measured flux traces implies equality of the time grid, the number of intervals, and all spatial profiles \(p_k\) [2111.02285]. No stability estimate is proved there, so the result is uniqueness without quantitative continuity.

In the space-semidiscrete stochastic problem, the Carleman estimate yields a stability theorem rather than merely uniqueness. The source is not only identifiable but determined in a Lipschitz-stable way from terminal and boundary-flux observations [2509.03760].

Regularization strategies differ correspondingly. The pseudo-parabolic Rothe method regularizes indirectly through time discretization, coercivity, energy bounds, and compactness [2502.04821]. The inverse potential formulation uses a priori nonnegativity, implemented through NNLS, and also discusses Tikhonov regularization
\[
\|A v - \tilde f\|^2 + \alpha \|v\|^2 \rightarrow \min_v
\]
for noisy data [2510.04693]. The Helmholtz frequency-semidiscrete method applies the quasi-reversibility method to a coupled elliptic system and proves convergence of the minimizers by means of a new Carleman estimate [1901.10047]. The hyperbolic inverse source paper transforms the problem to a Volterra integro-differential equation and proves Lipschitz-like convergence of the quasi-reversibility method:
\[
\|p_{\delta,\alpha}-p^*\|_{L^2(\Omega)}^2
\le
C\Big(\delta^2 + \alpha \|w^*\|_{H^3(\Omega\times[0,T])}^2\Big)
\]
[1806.03921].

A recurrent misconception is that semidiscretization is purely numerical. The literature indicates otherwise. In some papers, semidiscreteness is an approximation device for a continuous inverse problem [2502.04821]; in others, it is the actual source model [2111.02285]; in others still, it is the natural analytic framework for a spatially discrete operator [2509.03760] or a finite-dimensional boundary source representation [2510.04693]. This suggests that the term denotes a class of inverse formulations rather than a single numerical technique.

## 6. Algorithms, numerical evidence, and broader research directions

The pseudo-parabolic Rothe algorithm is explicitly stepwise: at each time level, compute \(h_i\) from the discrete reconstruction formula and then solve the linear elliptic variational problem for \(u_i\) [2502.04821]. In one-dimensional experiments on \(\Omega=(0,1)\) with
\[
\eta(t,x)=0.5,\quad \kappa(t,x)=t+1,\quad f(s)=s,\quad p(t,x)=\sin(\pi x),
\]
Backward Euler with \(\tau=T/200=0.005\), P1 Lagrange finite elements on 200 elements, and implementation in FEniCSx (DOLFINx 0.9.0), the paper reports that the reconstructed \(h(t)\) closely matches the exact source for all tested noise levels and that, in the noise-free case,
\[
E_{\max}^u(\tau)=\mathcal O(\tau),\qquad
E_{\max}^h(\tau)=\mathcal O(\tau),
\]
indicating first-order convergence in time [2502.04821].

The parabolic semi-discrete source paper adopts a Bayesian sequential prediction viewpoint. In the numerical section, the source states are parameterized as moving square patches on \([0,1]^2\), the time grid is known, and a Sequential Importance Resampling particle filter is used with \(N_p=320\) particles. Across two multiscale permeability fields and two source trajectories, the estimated positions closely match the true trajectories [2111.02285]. This does not numerically reconstruct general \(p_k\in L^2(\Omega)\), but it demonstrates that the semidiscrete temporal state-sequence model is compatible with sequential Bayesian inversion.

In inverse potential theory, the NNLS formulation is embedded in a moving-window search over candidate source domains. The residual
\[
r(x_0,y_0)=\min_{v\ge0}\|A v-f\|
\]
and the total mass
\[
m(x_0,y_0)\approx \sum_j v_j |S_j|
\]
act as localization indicators. The paper reports that low residuals form a plateau when the window contains the true sources, whereas residuals increase when part of the source lies outside the candidate domain [2510.04693].

The frequency-semidiscrete Helmholtz method combines a basis truncation in \(k\), a coupled elliptic PDE system in \(\mathbf x\), and a discretized quasi-reversibility solver. In 2D tests with 5% noise, the reconstructed shapes and amplitudes of the source are reported to agree well with the truth both for Cauchy data and Dirichlet-only data, although convergence is proved only for the overdetermined formulation [1901.10047].

The hyperbolic quasi-reversibility paper gives a fully discrete implementation of the Volterra-based reformulation. Noisy boundary data are differentiated in time by a Tikhonov-regularized linear system, then a regularized normal equation is solved for the discrete approximation of the auxiliary field \(w\), and finally the source is reconstructed from the discrete analogue of the identity relating \(p\) to \(w(\cdot,0)\). Numerical tests in two dimensions with noise levels \(2\%\), \(5\%\), and \(10\%\) show visually accurate recovery of shapes and amplitudes [1806.03921].

Broader directions already present in the literature include stochastic semidiscrete inverse sources in arbitrary dimensions [2509.03760], time-fractional inverse source problems with discrete random noise and spectral cutoff regularization [1606.06371], discrete Helmholtz multi-frequency inverse sources with phased and phaseless data [2401.14103], and temporal-source recovery in coupled subdiffusion systems, where the paper explicitly states that the continuous framework is one “you would want to semidiscretize in space” [2603.11700]. A plausible implication is that semidiscrete inverse source problems now function as an interface between PDE identifiability theory, finite-dimensional regularization, and computational state-space inversion rather than as a narrowly numerical subtopic.

Source: https://www.emergentmind.com/topics/semidiscrete-inverse-source-problem