1. Introduction

The age-old Riemann Hypothesis is a conjecture that the Riemann Zeta function has its zeros only at the negative even integers and complex numbers with real part $\frac{1}{2}$ and is firstly stated in the essay by Bernard Riemann [1].

The fact that the (warped) Riemann Zeta function $\zeta \left(\frac{1}{2}+iv\right)$ under certain conditions will arise as the solution of a second order ordinary differential equation on the time domain while the Riemann Zeta function $\zeta \left(1-v\right)$ may arise as the same solution but on the warped frequency domain (domain of the Dirichlet transform (93) image functions defined originally on a time domain) gave rise to possible presumptions that the two mentioned functions may be a Dirichlet Transform of each other. Since the extended form of the Riemann Zeta function is found to be an eigenfunction of the Fourier Transform operator, a composite Mellin Transform is found that projects $\zeta \left(\frac{1}{2}+2\pi iv\right)$ into $\zeta \left(1-v\right)$ preserving scale properties and is de facto a line-invariant transform between $\zeta \left(\frac{1}{2}+2\pi iv\right)$ and $\zeta \left(v\right)$,since there exists a well-known functional relationship between $\zeta \left(1-v\right)$ and $\zeta \left(v\right)$,the so-called functional equation of the extended Riemann Zeta function

$\varsigma \left(s\right)=2^s{\pi }^{s-1}\sin \left(\frac{\pi s}{2}\right)\Gamma \left(1-s\right)\varsigma \left(1-s\right)$ (1)

A fact that confirms the age-old conjecture of Bernard Riemann.

The $C^{\infty }$ -form of the eigenfunctions of the Fourier transform and the Dilation operator are thoroughly studied in [2]. The discrete and collapsed forms of these hyperbolic (eigen)functions are used as analytical input signals in [3] [4] [5] [6] [7] for an analytical signal analyzing affine model.

In this manuscript, the discrete super-position of the mentioned (warped) eigenfunctions are used to find a new eigenfunction of the Fourier Transform and to build solutions of the known and studied second order Ordinary Differential Equations (ODEs) on both time and warped frequency domain.

The solutions of the momentum equation of the Navier Stokes for both 2d and 3d versions are related to Dirichlet series (95), on the warped frequency domain. An appropriate inverse Dirichlet transform may be used to retrieve the solution on the time domain.

A general solution of the Navier Stokes 3d equation for the viscid incompressible flows is constructed from the begin conditions in the warped frequency domain, using geometric properties of the 2d version. This solution is called a minimal solution to the 3d Navier Stokes equation for the viscid incompressible flows. Those types of functions may contain domain singularities, which will cause the blow-up time when the numerical pivots strike a singularity. Please do see Appendix A for a study of the Navier Stokes d.e.

2. The Euler Differential Equation and a Think Experiment

Consider the solutions for the Euler differential equation

$t^2{y}^{″}+at{y}^{\prime }+by=0$ (2)

that may be deducted from the second order polynomial characteristic equation

$r^2+\left(a-1\right)r+b=0$ (3)

[With substitution of $t\frac{\text{d}y}{\text{d}t}=\frac{\text{d}y}{\text{d}x}$ to eliminate the non-constant terms does $t^2{y}^{″}+at{y}^{\prime }+by=0$ transpose to $y^{″}-\left(a-1\right)y^{\prime }+b=0$ and followed by the substitution of $y\left(t\right)={\text{e}}^{rx}$ delivers the found characteristic polynomial Equation (3).]

The solutions of the Euler differential equation [8] is given by

$y=\{\begin{array}{l}{c}_1{\left|t\right|}^{r_1}+c_2{\left|t\right|}^{r_2}\text{for}D>0\\ \left(c_1+c_2\ln \left|t\right|\right){\left|t\right|}^{\alpha }\text{for}D=0\\ \left[c_1\cos \left(\beta \ln \left|t\right|\right)+c_2\sin \left(\beta \ln \left|t\right|\right)\right]{\left|t\right|}^{\alpha }\text{for}D<0\end{array}$ (4)

With

$\{\begin{array}{l}\alpha =\frac{1}{2}\left(1-a\right)\\ \beta =\frac{1}{2}\sqrt{D}\end{array}$ (5)

and D the discriminant equivalent with

${\left(a-1\right)}^2-4b$ (6)

If we think the variable t as restricted to the set of natural numbers $\mathbb{N}$,and $c_1$ and $c_2$ as constants and equal to one and take improper super-position of all identified solutions, then the Riemann Zeta Function will arise as the family bundle of solutions of (2) as follows:

$y=\{\begin{array}{l}{c}_1\zeta \left(-\alpha -\beta \right)+c_2\zeta \left(-\alpha +\beta \right)\text{for}D>0\\ c_1\zeta \left(-\alpha \right)-c_2{\zeta }^{\prime }\left(-\beta \right)\text{for}D=0\\ \zeta \left(-\alpha +i\beta \right)\text{for}D<0\end{array}$ (7)

with the zeros of the characteristic polynomial (3) defined as $r_{1,2}=\alpha \pm \beta$ for

$D\ge 0$ (8)

and as $r_{1,2}=\alpha \pm i\beta$ when

$D<0$ (9)

In our think experiment we have got both the Riemann Zeta function $\zeta \left(-\alpha -\beta \right)$ and the warped Riemann Zeta function $\zeta \left(-\alpha +i\beta \right)$ and the derivative of the Riemann Zeta function ${\zeta }^{\prime }\left(-\beta \right)$ as solutions to the Equation (2).

The Riemann Zeta function arises also as a special case of the super-position of the general solutions (see TableB1 for examples in Appendix B) of other second order ordinary differential equations like the confluent hyper geometric differential equation called hyper geometric function of the first and second kinds, generalized Laguerre d.e., the Bessel d.e. etc. and the incompressible Navier-stokes equation

$\frac{\text{d}u}{\text{d}t}+\left(u\cdot \nabla \right)u-\nu \Delta u=\nabla p+g$ (10)

[9],but in the warped frequency domain. We would have expected then that the Riemann Zeta functions $\zeta \left(\frac{1}{2}+iv\right)$, $\zeta \left(\pm v\right)$ would therefore reside on the

warped frequency domain (as an image of the Dirichlet Transform operator). This is very fundamental in the sense that with proper assumptions and a proper time-frequency transformation these types of functions can co-exist both on the time and warped frequency domain as we are going to see next. Therefore it was very important for us to analyze the behavior of such functions on the warped frequency domain. Appendix B consists of a study of the relationship between discrete solutions of selected 2^nd order ODEs and the Riemann Zeta Function. Eigenfunctions of the Fourier Transform and the Dilation operators are thoroughly discussed in the thesis of Garas in [2] and the denomination “warped frequency domain” that is used in this manuscript is adapted from [2].

The Mellin transform operator [10] [11] does consist of the warping operator $U_x$,[2] p. 122 and further, followed by a Fourier transform. Its image will therefore reside on the warped frequency domain, when the initial function does reside on the time domain.

3. The Extended Riemann Zeta Function as an Eigenfunction of the Fourier Transform

Consider the monomial represented by ${\left|x\right|}^{\alpha }$ and its defined unitary Fourier transform (with ordinary frequency) $-2\sin \left(\frac{\pi \alpha }{2}\right)\Gamma \left(1+\alpha \right){\left|2\pi \upsilon \right|}^{-\alpha -1}$. If we take the improper sum of both ${\left|x\right|}^{\alpha }$ and its Fourier Transform, a procedure also known as the Poisson Summation formula, then will the extended form of the Riemann Zeta function in the complex plane emerge as follows,

$\mathcal{F}{\displaystyle {\sum }_{\left|x\right|=1}^{\left|x\right|=\infty }{\left|x\right|}^{\alpha }}={\displaystyle {\sum }_{\left|\upsilon \right|=1}^{\left|\upsilon \right|=\infty }-2\sin \left(\frac{\pi \alpha }{2}\right)\Gamma \left(1+\alpha \right){\left|2\pi \upsilon \right|}^{-\alpha -1}}=\varsigma \left(-\alpha \right)$ (11)

using the functional equation for the extended Riemann Zeta function (1). In other words the extended Riemann Zeta function is invariant under the Fourier Transform.

This is a very important find because we can now derive a general formula for the Mellin Transform with fractional derivatives of the Fourier Transform and use it to calculate the Mellin Transform of the Riemann Zeta function.

Hence

${\mathcal{M}}_f^{\xi }\left[\mathcal{X}\right]\left(\beta \right)={\displaystyle {\int }_0^{\infty }\mathcal{X}\left(f\right){\text{e}}^{2\pi i\xi f}{f}^{2\pi i\beta +r}\text{d}f}={\displaystyle {\int }_{-\infty }^{\infty }\mathcal{X}\left(f\right){\text{e}}^{2\pi i\xi f}{f}^{2\pi i\beta +r}\text{d}f}$ (12)

[6] if all frequency components of $\mathcal{X}\left(f\right)$ are positive (or if $\mathcal{X}\left(f\right)$ is an analytical signal), can be seen as a Fourier Transform with transform variable $-\left|\xi \right|$ and with

${\mathcal{M}}_f^{\xi =0}\left[\mathcal{X}\right]\left(\beta \right)={\mathcal{M}}_f\left[\mathcal{X}\left(f\right)\right]\left(\beta \right)$ (13)

Then is

${\mathcal{M}}_f^{-\left|\xi \right|}\left[\mathcal{X}\right]\left(\beta \right)=\frac{i^{\mu }}{2\pi }\frac{{\partial }^{\mu }\mathcal{F}\left[\mathcal{X}\right]\left(f\right)\left(\left|\xi \right|\right)}{\partial {\left|\xi \right|}^{\mu }}$ (14)

with $\mathcal{F}\left[\mathcal{X}\right]\left(f\right)\left(\xi \right)={\displaystyle {\int }_{-\infty }^{\infty }\mathcal{X}\left(f\right){\text{e}}^{-2\pi i\xi f}\text{d}f}$ (15)

and using the formula for the Fourier Transform of the μ-th derivative of a Fourier integrable function f,

$\mathcal{F}{t}^{\mu }f\left(t\right)\left(\xi \right)=i^{\mu }\frac{{\partial }^{\mu }\mathcal{F}\left[f\right]\left(t\right)\left(\xi \right)}{\partial {\xi }^{\mu }}$ (16)

with i the imaginary unit and $\mu \in ℂ$ extended from the set of natural numbers $\mathbb{N}$ to the set of complex numbers $ℂ$ and with

$\mu =2\pi i\beta +r$ (17)

If we take the Riemann Zeta function as the analytical signal

$\mathcal{X}\left(f\right)={\displaystyle {\sum }_{\left|f\right|=1}^{\left|f\right|=\infty }{\left|f\right|}^{-\alpha }}$ (18)

will above formula yield for

$\begin{array}{c}{\mathcal{M}}_f^{-\xi }\left[\mathcal{X}\right]\left(\beta \right)={\displaystyle {\int }_{-\infty }^{\infty }{\displaystyle {\sum }_{\left|f\right|=1}^{\left|f\right|=\infty }{\left|f\right|}^{-\alpha }{\text{e}}^{-2\pi i\xi f}{f}^{2\pi i\beta +r}\text{d}f}}\\ =\frac{i^{\mu }}{2\pi }\frac{{\partial }^{\mu }\mathcal{F}{\displaystyle {\sum }_{\left|f\right|=1}^{\left|f\right|=\infty }{\left|f\right|}^{-\alpha }\left(\xi \right)}}{\partial {\xi }^{\mu }}\\ =\frac{i^{\mu }}{2\pi }\mathcal{F}\frac{{\partial }^{\mu }{\displaystyle {\sum }_{\left|f\right|=1}^{\left|f\right|=\infty }{\left|f\right|}^{-\alpha }\left(\xi \right)}}{\partial f^{\mu }}\\ =\mathcal{F}{\displaystyle {\sum }_{\left|f\right|=1}^{\left|f\right|=\infty }\stackrel{\text{scalingeffect}}{\stackrel{︷}{2\pi i^{-\mu }\frac{\Gamma \left(\mu +\alpha \right)}{\Gamma \left(\alpha \right)}}}{\xi }^{\stackrel{\text{dilativeeffect}}{\stackrel{︷}{-\mu -\alpha }}}}\\ =2\pi i^{-\mu }\frac{\Gamma \left(\mu +\alpha \right)}{\Gamma \left(\alpha \right)}\varsigma \left(\mu +\alpha \right)\end{array}$ (19)

with $\mu =2\pi i\beta +r$ using the fact that the fractional derivative of

$\frac{{\partial }^{\mu }{x}^{-\alpha }}{\partial x^{\mu }}=-1^{\mu }\frac{\Gamma \left(\mu +\alpha \right)}{\Gamma \left(\alpha \right)}{x}^{-\left(\alpha +\mu \right)}$ (20)

for $\alpha \ge 0$ and $\xi \ge 0$ [7]. The identified scaling and dilative effect (or Doppler effect) makes this transform ideal for analytical signal reconstruction and signal analysis purposes applied in radar technology.

4. The Proof of the Age-Old Riemann Hypothesis

For normalization purposes and to be coherent with the used Fourier transform we will use $\zeta \left(\frac{1}{2}+2\pi iv\right)$ instead of $\zeta \left(\frac{1}{2}+iv\right)$.

An affine or a line-preserving map between $\zeta \left(\frac{1}{2}+2\pi iv\right)$ and $\zeta \left(1-s\right)$ may be formed by the composite Mellin transform operator ${\mathcal{M}}_v{\mathcal{M}}_n$,with

${\mathcal{M}}_{n}f=\mathcal{F}{U}_{n}f=\mathcal{F}{\text{e}}^{-\frac{n}{2}}f\left({\text{e}}^n\right)$ (21)

and $U_n$ the warping operator on variable n and $\mathcal{F}$ the Fourier transform. ( ${\mathcal{M}}_n$,similarly with the discrete Mellin Transform used in [6],entails a manipulation with the summation index variable to calculate the result). Using the found composite Mellin transform on $\zeta \left(\frac{1}{2}+2\pi iv\right)$ will yield

$\begin{array}{l}{\mathcal{M}}_v{\mathcal{M}}_n\zeta \left(\frac{1}{2}+2\pi iv\right)\\ ={\mathcal{M}}_v{U}_n\mathcal{F}\zeta \left(\frac{1}{2}+2\pi iv\right)={\mathcal{M}}_v\mathcal{F}{\displaystyle {\sum }_{n=1}^{n=\infty }{\text{e}}^{-\frac{n}{2}}\frac{{\text{e}}^{-2\pi iv\log \left({\text{e}}^{-n}\right)}}{\sqrt{{\text{e}}^{-n}}}}\\ ={\mathcal{M}}_v\mathcal{F}{\displaystyle {\sum }_{n=1}^{n=\infty }{\text{e}}^{2\pi ivn}}={\mathcal{M}}_v{\displaystyle {\sum }_{n=1}^{n=\infty }\delta \left(v-n\right)}=\zeta \left(1-s\right)\end{array}$ (22)

Q.e.d. The “quod erat demonstrandum” abbreviation refers to the confirmation of the age-old Riemann Hypothesis, knowing that there exists a functional equation between $\zeta \left(1-s\right)$ and $\zeta \left(s\right)$.

The popularized expression for the Riemann Zeta function in the complex plane on the so-called critical line, $\zeta \left(\frac{1}{2}+is\right)$,is also a solution of higher order

ODEs when the pair of complex conjugate solutions of the corresponding characteristic polynomials are repeating as we have seen in our thought experiment. Its warped version (thus warping of a warped function with ${\mathcal{M}}^2$ ) can be represented as an even characteristic polynomial, since the Riemann Zeta function is an eigenfunction of the Fourier transform. Hence

$\begin{array}{l}{\mathcal{M}}_n\zeta \left(\frac{1}{2}+2\pi iv\right)\stackrel{\text{using}(22)}{\stackrel{︷}{=}}{\mathcal{M}}_s^{-1}\left\{\zeta \left(1-s\right)\right\}=2{\displaystyle {\sum }_{v=1}^{v=\infty }\cos \left(v2\pi x\right)}\\ =2{\displaystyle {\sum }_{v=1}^{v=\infty }{T}_v\left(\cos \left(2\pi x\right)\right)}={\lim }_{v\to \infty }{U}_{2v}\left(\cos \pi x\right)-1\\ =a_{2v}{\tilde{x}}^{2v}+a_{2v-2}{\tilde{x}}^{2v-2}+\cdots +a_2{\tilde{x}}^2-1\end{array}$ (23)

and with $\tilde{x}=\cos \pi x$,with $a_{2v}$ defined as in the On-Line Encyclopedia of Integer Sequences as OAEIS A053117 [13],that is comparable with a characteristic polynomial of a higher order ODE for some $v\to \infty$. Here are $T_v$ and respectively $U_v$ the Chebyshev polynomials of resp. the first and second kinds [14]. Since the index v has to be taken even in the cumulative $U_v$ polynomial gives the opportunity to form repeating pairs of conjugate complex roots as the resulting characteristic polynomial may then be decomposed in product of lower order parabolic expressions which is a conditio sine qua non for the warped Riemann Zeta function to exist as a solution of the discussed ODE.

The expression $-2{\displaystyle {\sum }_{v=1}^{v=\infty }\cos \left(2\pi xv\right)}$ is equivalent with the Fourier series of the zero-th Bernoulli polynomials,

$B_n\left(x\right)=\frac{n!}{2\pi i^n}{\displaystyle {\sum }_{v\ne 0}\frac{{\text{e}}^{2\pi ivx}}{v^n}}=-2n!{\displaystyle {\sum }_{v\in \mathbb{N}/\left\{0\right\}}\frac{\cos \left(2\pi xv-\frac{n\pi }{2}\right)}{2v{\pi }^n}}$ (24)

The last definition can be found in the public accessible on-line Wikipedia.

If we do interchange the variable x and n with each other, a logic way to calculate the Dirichlet Transform of an Ordinary Generating Function resulting in a Dirichlet series, yields:

$B_x\left(n\right)=-2^{1-x}{\pi }^{-x}x!{\displaystyle {\sum }_{v=1}^{v=\infty }{k}^{-x}\sin \left(\frac{\pi }{2}\left[1-4vn+x\right]\right)}$ for

$0\le x\le 1$ (25)

what equivalent is with $-\zeta \left(1-s\right)$ for $n=0$ and is called the Hurwitz formula for the generalized extended zeta function $\zeta \left(x,n\right)$ and with $n\in \mathbb{Z}$ the well-known functional equation for the Riemann Zeta function will emerge.

The extended generalized zeta function is invariant under the Dirichlet Transform (and so under the Mellin Transform), when transforming only the domain of its arguments will undergo a variable swap as we have seen. The Hurwitz Zeta function is therefore an eigenfunction under both the Dirichlet and the Mellin transforms.

$\stackrel{{\mathcal{M}}_x}{\stackrel{︷}{\underset{}{\zeta \left(\frac{1}{2}+2\pi iv\right)}\to 2{\displaystyle {\sum }_{v=1}^{v=\infty }{T}_v\left(\cos \left(2\pi y\right)\right)}\equiv \underset{={\displaystyle {\sum }_{v>0}{a}_{2v}{\tilde{x}}^{2v}}-1}{\underset{︸}{{\lim }_{v\to \infty }{U}_{2v}\left(\cos \pi y\right)-1}}}}$ (26)

$\underset{{\mathcal{M}}_y}{\underset{︸}{2{\displaystyle {\sum }_{v=1}^{v=\infty }{T}_v\left(\cos \left(2\pi y\right)\right)}\equiv {\lim }_{v\to \infty }{U}_{2v}\left(\cos \pi y\right)-1\to \stackrel{\begin{array}{c}\text{basicDirichletseries}\\ \text{withnegativeargument}\end{array}}{\stackrel{︷}{\zeta \left(1-s\right)}}}}$ (27)

$\underset{\text{DirichletTransform}{D}_{s}f\left(x\right)=\frac{1}{\Gamma \left(s\right)}{\displaystyle {\int }_0^{\infty }{z}^{s-1}{f\left(x{\text{e}}^{-z}\right)|}_{x=1}\text{d}z}(93)\text{usedincombinatoricstheory}}{\underset{︸}{\underset{\text{OGF}}{\underset{︸}{{\displaystyle {\sum }_{v>0}{a}_v{x}^v}}}\to \underset{\text{Dirichletseries}}{\underset{︸}{{\displaystyle {\sum }_{v>0}{a}_v{v}^{-s}}}}}}$ (28)

[15] p. 177.

5. Conclusions and Future Work

Using the fact that the extended Riemann Zeta function is an eigenfunction under the Fourier Transform, it is proved that there exists a line-invariant composite Mellin transform operator that projects $\zeta \left(\frac{1}{2}+2\pi ix\right)$ into $\zeta \left(1-s\right)$ namely

${\mathcal{M}}_x{\mathcal{M}}_n\zeta \left(\frac{1}{2}+2\pi ix\right)=\zeta \left(1-s\right)$ (29)

what confirms the validation of the age-old Riemann Hypothesis, knowing that there exists a well-known functional equation between $\zeta \left(1-s\right)$ and $\zeta \left(s\right)$.

De facto they are conjugate pairs of each other using repetitively a fractional derivative Fourier Transform operator or the so-called Mellin Transform operator ${\mathcal{M}}_f^{\xi }$ (12). The Mellin transform will expand the transform variable for ${\mathcal{M}}_n$ from the set of natural numbers $\mathbb{N}$ to the set of real numbers $ℝ$ and ${\mathcal{M}}_x$ will expand the transform variable from the set of real numbers to the set of (negative) complex numbers $ℂ$.

Our future work will entail the study of the full spectrum of the eigenspaces belonging to the (fractional derivative) Fourier Transform operator and the behavior of the so-called conjugate pairs on their respectively domains with regard to known solution procedures for partial differential equations.

Appendix A

Minimal solution for the 3d viscid incompressible Navier Stokes equation: an application to physics

The solution of the Laplace equation in cylindrical coordinates is of paramount importance in many applications and this is not an exception for finding a minimal solution of the 3d Navier Stokes equations for the viscid incompressible media. The viscid flow arises when the viscous forces are prominent present in the equations. The solution is called minimal because a 2d geometric solution is extended to the 3d case, where domain singularities might be inherited.

If we extend the solutions to the complex space and take the super-positioned result of all positive discretized solutions will yield a polynomial or hyperbolic progression, which is not surprising as we are going to see.

Hence all second order ordinary differential equations that we have studied in this manuscript do have a polynomial progression as a solution as we are going the see in the next paragraphs.

The solution of the Mellin Transformed Laplace equation in cylindrical coordinates on the radius variable

$\{\begin{array}{l}{r}^2{f}_{rr}+r{f}_r+r{f}_{\theta \theta }=0\\ \frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\text{d}f}{\text{d}r}\right)+\frac{1}{r^2}\frac{{\partial }^{2}f}{\partial {\theta }^2}=0\\ F_{\theta \theta }+s^{2}F=0\end{array}$ (30)

is $A\left(s\right)\cos \left(\theta s\right)+B\left(s\right)\sin \left(\theta s\right)$ (31)

If we discretize $\theta =n\in \mathbb{N}$,then take the superposition of all discrete solutions, is still a general solution of above equation, namely

${\displaystyle {\sum }_1^{\infty }F\left(s,n\right)}={\displaystyle {\sum }_1^{\infty }{C}_1\left(s\right)\cos \left(sn\right)+C_2\left(s\right)\sin \left(sn\right)}$ (32)

and using $2\pi s$ instead of s and taking $A\left(s\right)=1$ and $B\left(s\right)=-i$ (to simplify the expression), then will the solution become

${\mathcal{M}}_s^{-1}{\displaystyle {\sum }_{n\in \mathbb{Z}/\left\{0\right\}}{\text{e}}^{-2\pi ins}}={\displaystyle {\sum }_{\begin{array}{l}n=-\infty \\ n\ne 0\end{array}}^{n=\infty }{x}^{2\pi in-\frac{1}{2}}}$ (33)

since ${\displaystyle {\sum }_{n\in \mathbb{Z}/\left\{0\right\}}{\text{e}}^{-2\pi ins}}$ is invariant under the Fourier Transform or

${\displaystyle {\sum }_{\left|x\right|}{\left|x\right|}^{2\pi in-\frac{1}{2}}}$ (34)

if we would have taken summation over the variable x.

If we take now the Dirichlet Transform of the solution (33) for $n>0$,will this result in $\phi \left(s\right)$ with $\phi \left(s\right)=f_4\left(s\right)={\displaystyle {\sum }_{n=1}^{n=\infty }{a}_n{n}^{-s}2\pi i^{-s}}$ if we would have used a general term $a_n$ [15]. Please note that $\phi \left(s\right)=f_4\left(s\right)={\displaystyle {\sum }_{n=1}^{n=\infty }{a}_n{n}^{-s}2\pi i^{-s}\Gamma \left(s\right)}$ for $a_n=1$,yields

$2\pi i^{-s}\Gamma \left(s\right)\varsigma \left(s\right)=\frac{\frac{i^{-s}}{2}\varsigma \left(1-s\right)}{\cos \left(\frac{\pi s}{2}\right)}=\frac{1}{2}\varsigma \left(1-s\right)-\frac{i}{2}\varsigma \left(1-s\right)\tan \left(\frac{\pi s}{2}\right)$

$=\varsigma \left(1-s\right)\left(\frac{1}{2}-\frac{i}{2}\tan \left(\frac{\pi s}{2}\right)\right)\cong \theta \left(-is\right)\varsigma \left(1-s\right)$

with $\theta \left(s\right)$ the Heaviside step function [16].

Above, the so-called functional equation of the Riemann Zeta function (1) is used, and

$\{\begin{array}{l}\Gamma \left(s\right)\Gamma \left(1-s\right)=\frac{\pi }{\sin \left(\pi s\right)}\\ \sin \left(\pi s\right)=\frac{1}{2}\cos \left(\frac{\pi s}{2}\right)\sin \left(\frac{\pi s}{2}\right)\\ \theta \left(s\right)\cong \frac{1}{2}+\frac{1}{2}\tanh \left(\frac{\pi s}{2}\right)=\{\begin{array}{l}1\text{for}s\to \infty \\ \frac{1}{2}\text{for}s=0\\ 0\text{for}s\to -\infty \end{array}\end{array}$ (35)

The Mellin Transform of (32), for $C_1\left(s\right)=f\left(s\right)$ and $C_2\left(s\right)=if\left(s\right)$,may be represented as

${\mathcal{M}}_x^{-\left|n\right|}\left[{\displaystyle {\sum }_1^{\infty }F\left(x,n\right)}\right]\left(\mu \right)={\mathcal{M}}_x{\displaystyle {\sum }_1^{\infty }{\text{e}}^{-i\left|n\right|x}f\left(\mu \right)}$ (36)

${\mathcal{M}}_x{\displaystyle {\sum }_1^{\infty }{\text{e}}^{-i\left|n\right|x}f\left(x\right)\left(\mu \right)}={\displaystyle {\sum }_1^{\infty }{\displaystyle {\int }_0^{\infty }{\text{e}}^{-i\left|n\right|x}f\left(x\right)x^{\mu }\text{d}x}}$ (37)

$={\displaystyle {\int }_0^{\infty }U\zeta \left(\frac{1}{2}+ix\right)\underset{\text{taylorseriesexpansion:}}{\underset{︸}{f\left(x\right)}}{x}^{\mu }\text{d}x}$

$={\displaystyle {\sum }_{k=0}^{k=\infty }{\displaystyle {\int }_0^{\infty }U\zeta \left(\frac{1}{2}+ix\right)x^{\mu }{f}^{\left(k\right)}\left(0\right)\frac{x^k}{k!}\text{d}x}}$

$\underset{\begin{array}{c}{\mathcal{M}}_x^0\left[U_n^{0,r}\mathcal{X}\right]\left(\mu +k-1\right):\\ \text{UnitaryFourierTransform}\\ \text{withordinaryfrequency}\end{array}}{\underset{︸}{=}}{\displaystyle {\sum }_{k=0}^{k=\infty }\frac{f^{\left(k\right)}\left(0\right)}{k!}{\displaystyle {\int }_0^{\infty }U\zeta \left(\frac{1}{2}+ix\right)x^{\mu +k}\text{d}x}}$

$={\displaystyle {\sum }_{k=0}^{k=\infty }\frac{f^{\left(k\right)}\left(0\right)}{k!}{\displaystyle {\int }_0^{\infty }{\text{e}}^{-2\pi in\mu }\zeta \left(\frac{1}{2}+ix\right)x^{\mu +k}\text{d}x}}$

$={\displaystyle {\sum }_{k=0}^{k=\infty }\frac{i^{\mu +k}}{2\pi }\frac{f^{\left(k\right)}\left(0\right)}{k!}\frac{{\partial }^{\mu +k}}{\partial n^{\mu +k}}\mathcal{F}\zeta \left(\frac{1}{2}+ix\right)}$ (38)

$={\displaystyle {\sum }_{k=0}^{k=\infty }\frac{i^{\mu +k}}{2\pi }\frac{f^{\left(k\right)}\left(0\right)}{k!}\frac{{\partial }^{\mu +k}}{\partial n^{\mu +k}}\zeta \left(\frac{1}{2}+ix\right)}$

$={\displaystyle {\sum }_{k=0}^{k=\infty }2\pi i^{-\mu -k}{f}^{\left(k\right)}\left(0\right)\frac{\Gamma \left(ix+\mu +k+\frac{1}{2}\right)}{\Gamma \left(k+1\right)\Gamma \left(ix+\frac{1}{2}\right)}\zeta \left(ix+\mu +k+\frac{1}{2}\right)}$

with ${\displaystyle {\sum }_0^{\infty }\frac{f^{\left(k\right)}\left(0\right)}{k!}}{x}^k$ (39)

the series expansion of the begin condition f around the value 0. Equation (38) is de facto a solution on the warped frequency domain to the eigenvalue problem Equation (30) or the Helmholtz equation that typically arises in physical problems involving partial differential equations in both space and time (e.g. the viscid Navier Stokes incompressible d.e.).

When the begin conditions differ per independent solution legs, the same method may be applied but now instead of using ${\text{e}}^{-i\left|n\right|x}$,its real part (with regard to the term with $C_1\left(s\right)$ ) must be used and respectively its imaginary part (with regard to the term with $C_2\left(s\right)$ ) must be used to calculate (38).

The time domain equivalent of (38) may be calculated by the inverse Dirichlet transform and is defined by

${\displaystyle {\sum }_{n=1}^{n=\infty }{\displaystyle {\sum }_{k=0}^{k=\infty }2\pi i^{-\mu -k}{f}^{\left(k\right)}\left(0\right)\frac{\Gamma \left(2n+k\right)}{\Gamma \left(k+1\right)\Gamma \left(n\right)\Gamma \left(n\right)}{x}^{\left(-2n-k\right)}}}$ (40)

taking $\mu =ix+\frac{1}{2}$ (41)

The Dirichlet Transform of a function f is defined as

$\frac{1}{\Gamma \left(\mu \right)}\mathcal{M}f\left({\text{e}}^x\right)=\frac{{\text{e}}^{-2\pi i\left|\xi \right|x}}{\Gamma \left(\mu \right)}\mathcal{M}f\left\{\mu \right\}=\frac{1}{\Gamma \left(\mu \right)}\frac{i^{\mu }}{2\pi }\frac{{\partial }^{\mu }}{\partial {\left|\xi \right|}^{\mu }}\mathcal{F}\left[f\right]\left(\mu \right)\left(\left|\xi \right|\right)$ (42)

what is equivalent with

$\frac{1}{\Gamma \left(\mu \right)}{\mathcal{M}}_x^{-\left|\xi \right|}f\left(x\right)\left(\mu \right)$ (43)

The above solution (38) can also be used for the reconstruction of the general solution (for the velocity profile, and thus also for the pressure warped frequency distribution) of the incompressible 2d-Navier Stokes equation in cylinder coordinates, taking

$\theta =n$ (44)

and $r=s$ (45)

since the Laplace equation in cylindrical coordinates is congruent with the 2d momentum equation for an incompressible flow. The solution for the (warped) velocity distribution function is therefore:

$u\left(\mu \right)={\displaystyle {\sum }_{k=0}^{k=\infty }2\pi i^{-\mu -k}{f}^{\left(k\right)}\left(0\right)\frac{\Gamma \left(ix+\mu +k+\frac{1}{2}\right)}{\Gamma \left(k+1\right)\Gamma \left(ix+\frac{1}{2}\right)}\zeta \left(ix+\mu +k+\frac{1}{2}\right)}$ (46)

With the begin function f equivalent with 1 does Equation (46) yield the Dirichlet series $2\pi i^{-\mu }\frac{\Gamma \left(2\mu \right)}{\Gamma \left(\mu \right)}\varsigma \left(2\mu \right)$ taking $\mu =ix+\frac{1}{2}$ :

$\{\begin{array}{l}\mu =ix+\frac{1}{2}\\ 2\pi i^{-\mu }\frac{\Gamma \left(2\mu \right)}{\Gamma \left(\mu \right)}\varsigma \left(2\mu \right)\text{warped frequency domain}\left(\text{see Figure}\text{A1}\right)\\ {\displaystyle {\sum }_{n=1}^{n=\infty }2\pi i^{-n}\frac{\Gamma \left(2n\right)}{\Gamma \left(n\right)\Gamma \left(n\right)}{x}^{-2n}}=-\frac{i\sqrt{\pi }}{2x^2{\left(\frac{\pi x^2+2i}{x^2}\right)}^{\frac{3}{2}}},\left|x\right|>\sqrt{\frac{2}{\pi }}\text{time}\text{domain}\left(\text{see Figure A2}\right)\end{array}$ (47)

FigureA1 confirms the volatile characteristics of the Dirichlet series. FigureA2 is the inverse Dirichlet transform of FigureA1, using [15].

The 3d Navier Stokes system of equations for the incompressible viscid flow [9] is reducible to an extent that its general and exact solutions are determined by the (general and exact) solutions of the 2d Laplace equation in cylindrical coordinates (as above) as follows. We will call this solution a minimal solution.

Consider the NVS 3d systems of equations for the incompressible viscid flow:

$\rho \frac{\partial u}{\partial t}+\rho \left(u\cdot \nabla \right)u-\mu {\nabla }^{2}u=-\nabla p+\rho g$ (48)

and $\{\begin{array}{l}\mu =\text{viscid constant}\\ \rho =\text{medium density}\end{array}$ (49)

Equation (48) per directional component in cylindrical coordinates:

$\begin{array}{l}\rho \left(\frac{\partial u_r}{\partial t}+u_r\frac{\partial u_r}{\partial r}+\frac{u_{\theta }}{r}\frac{\partial u_r}{\partial \theta }-\frac{u_{\theta }^2}{r}+u_z\frac{\partial u_r}{\partial z}\right)\\ =-\frac{\partial P}{\partial r}+\rho g_r+\mu \left(\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial u_r}{\partial r}\right)-\frac{u_r}{r^2}+\frac{1}{r^2}\frac{{\partial }^2{u}_r}{\partial {\theta }^2}-\frac{2}{r^2}\frac{\partial u_{\theta }}{\partial \theta }+\frac{{\partial }^2{u}_r}{\partial z^2}\right)\end{array}$ (50)

$\begin{array}{l}\rho \left(\frac{\partial u_{\theta }}{\partial t}+u_r\frac{\partial u_{\theta }}{\partial r}+\frac{u_{\theta }}{r}\frac{\partial u_{\theta }}{\partial \theta }+\frac{u_{\theta }{u}_r}{r}+u_z\frac{\partial u_{\theta }}{\partial z}\right)\\ =-\frac{\partial P}{\partial \theta }+\rho g_{\theta }+\mu \left(\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial u_{\theta }}{\partial r}\right)-\frac{u_{\theta }}{r^2}+\frac{1}{r^2}\frac{{\partial }^2{u}_{\theta }}{\partial {\theta }^2}+\frac{2}{r^2}\frac{\partial u_r}{\partial \theta }+\frac{{\partial }^2{u}_{\theta }}{\partial z^2}\right)\end{array}$ (51)

$\begin{array}{l}\rho \left(\frac{\partial u_z}{\partial t}+u_r\frac{\partial u_z}{\partial r}+\frac{u_{\theta }}{r}\frac{\partial u_z}{\partial \theta }+u_z\frac{\partial u_z}{\partial z}\right)\\ =-\frac{\partial P}{\partial z}+\rho g_z+\mu \left(\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial u_z}{\partial r}\right)+\frac{1}{r^2}\frac{{\partial }^2{u}_z}{\partial {\theta }^2}+\frac{{\partial }^2{u}_z}{\partial z^2}\right)\end{array}$ (52)

To solve the above 3d NVS system of equation we will have to reduce the r.h.s. of each equation to the sum of the pressure term and the external force and equating the Laplace expression in cylindrical coordinates to zero (solve up the Laplacian to get the velocity distribution), and then expand this using time convolution to $u\left(x,t\right)$ with the appropriate Gaussian diffusion Green’s function to solve the enclosed diffusion equation

$\rho \frac{\partial u}{\partial t}=\mu {\nabla }^{2}u$ (53)

and then use it to calculate the pressure distribution as follows

$\nabla p=\rho g-\rho \left(u\cdot \nabla \right)u$ (54)

To simplify this we will expand the super-position of discrete solutions Equation (32) over the whole set of integers $\mathbb{Z}$,then

$u\left(x,t\right)=\underset{\text{convolution integral}}{\underset{︸}{{\displaystyle {\int }_0^{\infty }{\displaystyle {\sum }_{k=0}^{k=\infty }\underset{{\displaystyle {\sum }_{k=-\infty }^{k=\infty }{\text{e}}^{i\pi ky}}}{\underset{︸}{\stackrel{\text{integrand of eq}\text{.(37)}}{\stackrel{︷}{2\pi \delta \left(y\right)}}}}\stackrel{k^{\text{th}}\text{Taylor term of}f\text{around}a}{\stackrel{︷}{\frac{{\left(y-a\right)}^k{f}^{\left(k\right)}\left(a\right)}{k!}}}\stackrel{\text{Green}’\text{s function of diffusion eq}\text{.}}{\stackrel{︷}{\theta \left(t\right){\frac{1}{4\pi \mu t}}^{\frac{1}{2}}{\text{e}}^{-\frac{{\left(x-y\right)}^2}{4\mu t}}}}\text{d}y}}}}$ (55)

$\rho =1={\displaystyle {\sum }_{k=0}^{k=\infty }\frac{1}{2}\sqrt{\frac{\pi }{\mu t}}{\text{e}}^{-\frac{{\left|x\right|}^2}{4\mu t}}\frac{-a^k{f}^{\left(k\right)}\left(a\right)}{k!}\theta \left(t\right)}=\frac{1}{2}\sqrt{\frac{\pi }{\mu t}}{\text{e}}^{-\frac{{\left|x\right|}^2}{4\mu t}}{\text{e}}^{-a\left(Df\right)a}\theta \left(t\right)$ (56)

The last function is rescaled along both axes compared to its time-domain version, and is called a “compressed” eigenfunction of the Mellin Transform operator, the corresponding transform pair on the time domain is equivalent with

$\pi \sqrt{\frac{1}{4\pi \mu t}}{\text{e}}^{\mu t{\left|x\right|}^2}{\text{e}}^{-a\left(Df\right)a}\theta \left(t\right)$ (57)

In two dimensions it will be congruent with

$\frac{\pi }{4\mu t}{\text{e}}^{\mu t\left({\left|x\right|}^2+{\left|y\right|}^2\right)}{\text{e}}^{-a\left(D{f}_1\right)a}{\text{e}}^{-a\left(D{f}_2\right)a}{\theta }^2\left(t\right)$ (58)

The interested reader may try to derived it, the warped frequency 3d pressure distribution in Cartesian coordinates would be then, taking

$g=0$ (59)

out of Equation (54) and disregarding the terms ${\text{e}}^{-a\left(D{f}_j\right)a}$ for $j=1,2,3$ and the Heaviside time restriction in the different directions,

$\left(\begin{array}{c}{P}_x\\ P_y\\ P_z\end{array}\right)=-\rho \left(\begin{array}{c}\frac{1}{2}\frac{\left|y\right|\sqrt{\pi }}{{\mu }^{2\frac{1}{2}}{t}^{2\frac{1}{2}}}k\text{Erf}\left(\frac{\left|x\right|}{\sqrt{4\mu t}}\right)\exp \left(-\frac{{\left|y\right|}^2}{4\mu t}\right)\exp \left(-\frac{{\left|z\right|}^2}{2\mu t}\right)-\frac{1}{2}\frac{\left|z\right|\sqrt{\pi }}{{\mu }^{2\frac{1}{2}}{t}^{2\frac{1}{2}}}j\text{Erf}\left(\frac{\left|x\right|}{\sqrt{4\mu t}}\right)\exp \left(-\frac{{\left|y\right|}^2}{2\mu t}\right)\exp \left(-\frac{{\left|z\right|}^2}{2\mu t}\right)(60)\\ \frac{1}{2}\frac{\left|x\right|\sqrt{\pi }}{{\mu }^{2\frac{1}{2}}{t}^{2\frac{1}{2}}}k\text{Erf}\left(\frac{\left|y\right|}{\sqrt{4\mu t}}\right)\exp \left(-\frac{{\left|z\right|}^2}{4\mu t}\right)\exp \left(-\frac{{\left|x\right|}^2}{2\mu t}\right)+\frac{1}{2}\frac{\left|z\right|\sqrt{\pi }}{{\mu }^{2\frac{1}{2}}{t}^{2\frac{1}{2}}}i\text{Erf}\left(\frac{\left|y\right|}{\sqrt{4\mu t}}\right)\exp \left(-\frac{{\left|z\right|}^2}{4\mu t}\right)\exp \left(-\frac{{\left|x\right|}^2}{2\mu t}\right)(61)\\ \frac{1}{2}\frac{\left|x\right|\sqrt{\pi }}{{\mu }^{2\frac{1}{2}}{t}^{2\frac{1}{2}}}j\text{Erf}\left(\frac{\left|z\right|}{\sqrt{4\mu t}}\right)\exp \left(-\frac{{\left|x\right|}^2}{4\mu t}\right)\exp \left(-\frac{{\left|y\right|}^2}{2\mu t}\right)-\frac{1}{2}\frac{\left|y\right|\sqrt{\pi }}{{\mu }^{2\frac{1}{2}}{t}^{2\frac{1}{2}}}i\text{Erf}\left(\frac{\left|z\right|}{\sqrt{4\mu t}}\right)\exp \left(-\frac{{\left|x\right|}^2}{2\mu t}\right)\exp \left(-\frac{{\left|y\right|}^2}{4\mu t}\right)(62)\end{array}\right\}$

using $\begin{array}{l}{u}_k\left(x,y,t\right)=\frac{k}{\mu t}\exp \left(-\frac{{\left|x\right|}^2}{4\mu t}\right)\exp \left(-\frac{{\left|y\right|}^2}{4\mu t}\right),\\ u_i\left(y,z,t\right)=\frac{i}{\mu t}\exp \left(-\frac{{\left|y\right|}^2}{4\mu t}\right)\exp \left(-\frac{{\left|z\right|}^2}{4\mu t}\right)\text{etc}\end{array}$ (63)

based on an i, j, k indexed classical Hamilton’ian unit vectors with magnitude $\sqrt{-1}$ where $i^2=j^2=k^2=ijk=-1$ and

$ij=k$, $jk=i$, $ki=j$, $ij=-ji$, $jk=-kj$, $ki=-ik$ and where Erf(x) the Error function is, defined as

$\text{Erf}\left(z\right)=\frac{2}{\sqrt{\pi }}{\displaystyle {\int }_0^z{\text{e}}^{-t^2}\text{d}t}$ (64)

On the edge of the domain the following is true for the normal vectors

$\frac{\partial u}{\partial n}=\frac{\partial u_k}{\partial z}\cong u\left(x+k\right)-u\left(x\right)=0$ etc. (65)

Please do see FigureA3. FigureA3 does show the contour graph of the warped frequency distribution in 3 dimensions.

The total energy of each wave system is equal on both the time and frequency domain [6] [7]. In symbols,

${\displaystyle \int {\displaystyle {\int }_{-\infty }^{\infty }P\left(t,f\right)\partial t\partial f}}={\displaystyle {\int }_{-\infty }^{\infty }{\left|\mathcal{X}\left(f\right)\right|}^2\partial f}={\displaystyle {\int }_{-\infty }^{\infty }{\left|x\left(t\right)\right|}^2\partial t}={\mathcal{E}}^2$ (66)

with above the expression for the Wigner-Ville Distribution function ( a class of dual energy distributions), and below, the integrand [is equivalent to the Ambiguity Function] as defined in [3] [6] [7]

$P\left(t,f\right)=f{\displaystyle \underset{ℝ}{\int }{\text{e}}^{-2\pi iuft}\frac{u^{2r+2}}{2\sinh \frac{u}{2}}\mathcal{X}\left(\frac{fu{\text{e}}^{\frac{-u}{2}}}{2\sinh \frac{u}{2}}\right){\mathcal{X}}^{*}\left(\frac{fu{\text{e}}^{\frac{u}{2}}}{2\sinh \frac{u}{2}}\right)\text{d}u}$ (67)

$\mathcal{X}\left(f\right)$ is the analytical signal on the frequency domain and its Fourier Transform inverse $x\left(t\right)$ on the time domain and ${\mathcal{E}}^2$ the total (kinetic) energy of the system and r is a warp parameter that may be set $-\frac{1}{2}$ to yield a factor

$\frac{u{\text{e}}^{\frac{u}{2}}}{{\text{e}}^u-1}$ (68)

in the integrand and what may be calculated with the Residue theorem. Expression (65) will yield an unbounded wave energy-level of ${\mathcal{E}}^2={\displaystyle {\int }_{-\infty }^{\infty }{\left|u\left(\mu \right)\right|}^2\partial \mu }$ for every component with $\zeta \left(1\right)=\infty$.

The Wigner-Ville time-frequency distribution function is a unitary affine function that suffices the following equation

${\displaystyle {\int }_{-\infty }^{\infty }\text{d}t}{\displaystyle {\int }_0^{\infty }P\left(t,f\right)\delta \left(t-\xi -\frac{\beta }{f}\right)f^{2r}\text{d}f}={\left|{\mathcal{M}}_f^{\xi }\mathcal{X}\left(\beta \right)\right|}^2$ (69)

with ${\mathcal{M}}_f^{\xi }\left[\mathcal{X}\right]\left(\beta \right)={\displaystyle {\int }_0^{\infty }\mathcal{X}\left(f\right){\text{e}}^{2\pi i\xi f}{f}^{2\pi i\beta +r}\text{d}f}$ (12).

The total wave energy does inherit the known singularities within the Dirichlet series from the components of the distribution function of $u\left(x\right)$ what

makes it also unbounded.

The last maximum value of the total energy is also supported by the fact that the Wigner-Ville Distribution function definition also does have a factor of the hyperbolic cosecant with half argument in its integrand (defined on $C^{\infty }\left(ℂ\backslash \left\{2\pi \mathbb{Z}\right\}\right)$,with poles at $\left\{2\pi i{\left\{k\right\}}_{k\in \mathbb{Z}}\right\}$ what may lead to variants of the improper Riemann contour integral [1]

$2\sin \left(\pi s\right)\Gamma \left(s\right)\varsigma \left(s\right)=i{\displaystyle {\int }_{\infty }^{\infty }\frac{{\left(-x\right)}^{s-1}}{{\text{e}}^x-1}\text{d}x}=i{\displaystyle \oint \frac{{\left(-x\right)}^{s-1}}{{\text{e}}^x-1}\text{d}x}$ (70)

Below schematically, the application of the presented formulas for solving up the Laplace equation in cylindrical coordinates in an analytic signal analyzing model used by [6]. The constructed transformation is scale invariance:

An analytic signal analyzing model

Out of the general solution is it clear that the distribution functions for $u\left(x\right)$ and $p\left(x\right)$ on the warped frequency domain do inherit the known singularities of the Dirichlet series, naturally depending on the begin-conditions which may cause shifts in singularities at $\left\{1\right\}$,removable singularities at

$-\frac{2\pi \mathbb{Z}}{\log p}i+\frac{\log a\left(p\right)}{\log p}$ (72)

with p a prime number as a consequence of the Euler denomination of the Dirichlet series ( $a\left({\mathbb{N}}^{+}\right)=1$ for the Riemann Zeta function) in the form of the product formula

${\displaystyle {\sum }_{n=1}^{n=\infty }a\left(n\right)n^{-s}}={\displaystyle {\prod }_{p\text{prime}}\frac{1}{1-a\left(p\right)p^{-s}}}$ (73)

only for $a\left(n\right)$ is totally multiplicative [17]. Solving up the zeros of the denominator in the last expression

${\text{e}}^{2\pi i\mathbb{Z}}-a\left(p\right)p^{-s}=0$ (74)

does result in the mentioned complex removable singularities. So may be said that the distribution functions $u\left(x\right)$ and $p\left(x\right)$ are defined on $C^{\infty }\left(ℂ\backslash \left\{\mathbb{Z}\right\}\right)$ for each component.

The begin-conditions may also dictate the periodicity of the general solution. There exists no general solution that is independent of the begin-solutions, since the solution is constructed out of the begin-conditions. The begin-conditions may cause singularities or may annihilate them, depending on the initial problem and its basic solution. The basic solution of the Laplace equation in cylindrical coordinates, which is the solution core for both the velocity and pressure frequency distributions (in the incompressible case) for the Navier Stokes equations in 2d and to certain extent also for the 3d version, does contain and consequently inherits the singularities observed in the cited improper Riemann contour integral in particular or from the Dirichlet series in general. The input signal

$\mathcal{X}\left(f\right)={\displaystyle {\sum }_1^{\infty }{f}^{\alpha }}$ (75)

of the used model may manifest or cause singularities on both time and warped frequency domains with both the summation parameter indexed by f (hyperbolas) and by $\alpha$ (polynomials). On the time domain the signal may experience singularities due to the vast powers that are used that may be unpractical for most commercial available micro-processors to compute [for example when f is chosen as the series of orthogonal polynomials like the Hermite polynomials or due to the known singularities of the Zeta function. On the warped frequency domain the singularities are due to the properties of the formed Dirichlet series output signals for both cases.

Interesting to note is that the so-called blow-up time within the numeric schemes may manifest whenever the numerical pivoting method strikes a known singularity of the general solution.

When the analytical signal condition is dropped or when negative discrete solutions of the Laplace equation (read the momentum equation) are allowed, the motion solution of the 3d Navier Stokes equations for the viscid incompressible flow will become a (multi-dimensional) Gaussian function with no singularities on neither the time domain nor the warped frequency domain.

Appendix B

The Riemann Zeta Function represented as a Dirichlet Transform of discrete super-position of second order ODEs solutions.

The confluent hyper geometric functions are solutions to the Weber differential equations, obtained after separation of variables when decomposing the Laplacian in cylindrical coordinates;

$\{\begin{array}{l}\frac{\partial u^2}{\partial x^2}-\left(c+k^2{x}^2\right)u=0\\ \frac{\partial v^2}{\partial y^2}-\left(c-k^2{x}^2\right)v=0\end{array}$ (76)

The substitution $u=w\cdot \exp \left(-\frac{x^2}{4}\right)$ converts the weber equation to

$\frac{\partial w^2}{\partial x^2}-x\frac{\partial w}{\partial x}+\upsilon w=0$ (77)

with solution $u=H_{\upsilon }\left(x\right)\cdot \exp \left(-\frac{x^2}{4}\right)$,where $H_{\upsilon }\left(x\right)$ is the Hermite polynomial.

In general the solutions of the second order ordinary confluent hypergeometric d.e.

$x\frac{\partial w^2}{\partial x^2}+\left(c-x\right)\frac{\partial w}{\partial x}-aw=0$ (78)

are of the form

$w=b_1{}_{1}F{}_1\left(a;c;x\right)+b_{2}U\left(a,c,x\right)$ (79)

and are called the confluent hypergeometric function of the first and second kind respectively. The first kind version is also denoted as $M\left(a,c,x\right)$ or $\Phi \left(a;c;x\right)$ with

${}_{1}F{}_1\left(a;c;x\right)={\displaystyle {\sum }_{n=0}^{n=\infty }\frac{a_{\left(n\right)}{z}^n}{c_{\left(n\right)}n!}}$ (80)

and $a_{\left(n\right)}$ the Pochhammer symbol defined as

$a_{\left(n\right)}=\frac{\Gamma \left(a+n\right)}{\Gamma \left(a\right)}=a\left(a+1\right)\cdots \left(a+n-1\right)$ (81).

An integral representation for $U\left(a,b,z\right)$ is

$U\left(a,b,z\right)=\frac{1}{\Gamma \left(a\right)}{\displaystyle {\int }_0^{\infty }{\text{e}}^{-z\cdot t}{t}^{a-1}{\left(1+t\right)}^{b-a-1}\text{d}t}$ (82)

[18] [19].

The confluent hypergeometric function of the second kind, can be transposed to the improper contour integral

$i{\displaystyle \oint \frac{{\left(-x\right)}^{s-1}}{{\text{e}}^x-1}\text{d}x}$

by choosing $b=a+1$,and $a=s$. Namely

$U\left(s,s+1,n\right)=\frac{1}{\Gamma \left(s\right)}{\displaystyle {\int }_0^{\infty }{\text{e}}^{-n\cdot t}{t}^{s-1}\text{d}t}$

$U\left(s,s+1,n\right)=n^{-s}$ (83)