In this contribution, we present a discrete explicit-in-time LB extension of the non-equilibrium extrapolation method for the use of ODE-based boundary conditions (previous usages of non-equilibrium extrapolation have been mostly confined to given Dirichlet-based pressure or velocity boundary conditions of the fluid system [1]}, [2]}, [3]}). In order to update the distribution functions from a timestep \(t_n\)  to a timestep \(t_{n+1}=t_n+\Delta t\)  of a fluid point \(x_f\in \Omega \)  that is streamed from the corresponding 0D-3D interface (boundary) node \(\mathbf {x}\in \partial \Omega _2\)  of the lattice (i.e., \(x_f = x+e_i\Delta t\) ), such a formulation can be derived by separating collision and streaming operations of Equation (REF ) into
\(f_i^*(x)=f_i(x,t_n)-\frac{1}{\tau }[f_i(x,t_n)-f_i^{eq}(x,t_n)]+\Delta t F_i(x,t_n), \quad \mathbf {x} \in \partial \Omega _2,\) 
