Macroscopic traffic flow is stochastic, but the physics-informed deep learning methods currently used in transportation literature embed deterministic PDEs and produce point-valued outputs; the stochasticity of the governing dynamics plays no role in the learned representation. This work develops a framework in which the physics constraint itself is distributional and directly derived from stochastic traffic-flow dynamics. Starting from an Ito-type Lighthill-Whitham-Richards model with Brownian forcing, we derive a one-point forward equation for the marginal traffic density at each spatial location. The spatial coupling induced by the conservation law appears as an explicit conditional drift term, which makes the closure requirement transparent. Based on this formulation, we derive an equivalent deterministic Probability Flow ODE that is pointwise evaluable and differentiable once a closure is specified. Incorporating this as a physics constraint, we then propose a score network with an advection-closure module, trainable by denoising score matching together with a Fokker-Planck residual loss. The resulting model targets a data-conditioned density distribution, from which point estimates, credible intervals, and congestion-risk measures can be computed. The framework provides a basis for distributional traffic-state estimation and for stochastic fundamental-diagram analysis in a physics-informed generative setting.
Figures & tables
Fan–Du SLWR CACIE 2023; Trans. B 2024
Shi–Di PIDL AAAI 2021
This paper
Stochastic model
Random-parameter PDE
None (det. ARZ)
Itô SPDE
Distributional eq.
None
N/A
Fokker–Planck
Physics in learning
Forward solver
Det. PDE residual
PF-ODE residual
Output
Mean & std
Point estimate
Full pθ(ρ^;x,t)
UQ mechanism
Parametric
Post-hoc only
Intrinsic (score)
Dynamic noise
No
No
Yes ( dWt )
Table 1: Positioning of the proposed framework relative to the closest prior works.
Symbol
Description
Domain and probability space
D=[0,L]
Highway segment
QT=D×(0,T]
Spatio-temporal domain
(Ω,F,{Ft},P)
Filtered probability space
Wtk(k=1,…,K)
Independent standard Brownian motions
Traffic variables and flux
Table 2: Summary of notation.
Figure 1: Architectural comparison between PIDL and the proposed framework. Both share a two-component structure (estimation network + physics evaluator), but the proposed framework operates at the distributional level: the network learns a score function, the physics evaluator regularises via the Fokker–Planck equation, and the output is a data-conditioned distributional estimate.
Component
PIDL
This paper
Network input
(x,t)
(ρ^,x,t)
Network output
(ρ^,u^) point values
sθ (score, scalar)
Physics evaluator
ARZ residual (f^1,f^2)
FPE residual Rθ,ϕ
Data loss
∑∣ρ^−ρobs∣2
DSM: ∑∣sθ+ϵ/σ∣2
Physics loss
∑∣f^1∣2+∣f^2∣2
∑∣Rθ,ϕ∣2
Collocation domain
(x,t) , 2D
(ρ^,x,t) , 3D
Table 3: Component-level comparison between PIDL and the proposed framework.