We consider the PDE flow associated to Riemann zeta and general DirichletLfunctions. These are models characterized by nonlinearities appearing in classical number theory problems, and generalizing the classical holomorphic Riemann flow studied by Broughan and Barnett. Each zero of a Dirichlet Lfunction is an exact solution of the model. In this paper, we first show local existence of bounded continuous solutions in the Duhamel sense to any DirichletLfunction flow with initial condition far from the pole (as long as this exists). In a second result, we prove global existence in the case of nonlinearities of the form DirichletLfunctions and data initially on the right of a possible pole at s = 1. Additional global wellposedness and convergence results are proved in the case of the defocussing Riemann zeta nonlinearity and initial data located on the real line and close to the trivial zeros of the zeta. The asymptotic stability of any stable zero is also proved. Finally, in the Riemann zeta case, we consider the "focusing" model, and prove blow-up of real-valued solutions near the pole s = 1. (c) 2025 Elsevier Inc. All rights are reserved, including those for text and data mining, AI training, and similar technologies.