Fact-checked by Grok 2 weeks ago

Volterra integral equation

A Volterra integral equation is a functional equation in which an unknown function appears linearly under an integral sign whose upper limit is the independent variable, distinguishing it from Fredholm integral equations with fixed limits.^[1] Named after the Italian mathematician Vito Volterra, who developed the foundational theory in a series of papers published in 1896 on the inversion of definite integrals, these equations form a core part of integral equation theory and are closely related to initial value problems for ordinary differential equations.^[2]^[3] Volterra's work on these equations began in the 1880s during his tenure at the University of Pisa, building on earlier studies in mathematical physics, but his 1896 publications provided the first systematic treatment, including iterated kernels and resolvent methods that influenced later developments by Ivar Fredholm and David Hilbert.^[2]^[4] The equations were further advanced by Traian Lalescu in his 1908 doctoral thesis, Sur les équations de Volterra, which established key existence and uniqueness results.^[3] Over the subsequent century, the theory expanded to include nonlinear variants, stochastic forms, fractional variants, and generalizations on time scales, with ongoing research in numerical methods and regularization techniques for ill-posed problems.^[3]^[5]^[6] Mathematically, Volterra integral equations are classified into two main types: those of the first kind, where the unknown function y appears only inside the integral, as in

\int_a^x K(x, t) y(t) \, dt = f(x),

and those of the second kind, where y also appears outside, as in

y(x) = g(x) + \int_a^x K(x, t) y(t) \, dt.

Here, K(x, t) is the kernel, f and g are given functions, and a is typically the initial point.^[1] These can be homogeneous or nonhomogeneous, linear or nonlinear, and may include delay terms or integro-differential components. Solutions to linear second-kind equations with continuous kernels are often unique and Lipschitz continuous, provable via the Banach fixed-point theorem in appropriate function spaces.^[1] For first-kind equations, which are more challenging due to potential ill-posedness, algebraic formulas and regularization methods have been developed, particularly in topological algebras.^[7] The theory of Volterra integral equations draws from functional analysis and operator theory, with resolvent kernels enabling successive approximations akin to Picard iteration for differential equations.^[8] Existence and uniqueness theorems, such as those by Corduneanu, extend to weakly singular kernels and abstract spaces, while numerical solutions employ product integration or variational methods for nonlinear cases.^[9]^[10] Volterra integral equations arise naturally in diverse scientific and engineering contexts, modeling phenomena where history-dependent processes are key. In biology, they describe population growth dynamics, extending Volterra's own predator-prey models.^[11] Applications in physics include viscoelasticity for material stress-strain relations and potential theory in electrostatics.^[1]^[12] In engineering, they appear in nuclear reactor kinetics, particle transport, and conformal mapping; in finance, they model ruin probabilities in risk theory.^[1]^[13] These equations thus bridge pure mathematics with applied sciences, underscoring their enduring relevance.

Definition and Fundamentals

General Form

A Volterra integral equation is defined as an equation in which the unknown function appears inside an integral whose upper limit is the independent variable t and lower limit is a fixed value a, typically with the integral depending on values of the function up to t.^[14] Volterra integral equations are classified into first kind and second kind. For the first kind, the unknown function u appears only under the integral:

\int_a^t K(t,s) u(s) \, ds = f(t), \quad t \in [a, T],

where f is a given continuous function and K(t,s) is the kernel. These are often ill-posed and require special treatment.^[1] The general form of a linear Volterra integral equation of the second kind is given by

u(t) = f(t) + \int_a^t K(t,s) u(s) \, ds, \quad t \in [a, T],

where f is a given continuous function, K(t,s) is the kernel defined on the domain \{(t,s) : a \leq s \leq t \leq T\}, and u is the unknown function to be determined.^[14] A special case is the convolution-type equation, where the kernel depends only on the difference t - s, i.e., K(t,s) = \tilde{K}(t - s).^[14] These equations often arise from converting initial value problems for ordinary differential equations (ODEs) into integral form by successive integration. For instance, starting from a first-order ODE y'(t) = g(t) + \lambda y(t) with y(a) = 0, integration from a to t yields the equivalent Volterra equation y(t) = \int_a^t g(s) \, ds + \lambda \int_a^t y(s) \, ds.^[15] Key properties of Volterra integral equations include causality, meaning the value of the solution at t depends only on the forcing function and the solution for s \leq t, reflecting a "memory" limited to past values. Under mild conditions, such as continuity of the kernel K and the function f, the equation is well-posed, possessing a unique continuous solution on [a, T].^[14] As an illustrative example, consider the simple convolution equation

f(t) = g(t) + \lambda \int_0^t f(s) \, ds, \quad t \geq 0,

where \lambda is a constant and g is given; this admits an explicit solution via successive approximation or Laplace transforms, demonstrating the solvability of such forms.^[14] Unlike Fredholm integral equations, which feature fixed integration limits, Volterra equations have variable limits tied to the independent variable.^[14]

Historical Background

The Volterra integral equation was introduced by Italian mathematician Vito Volterra in the 1890s as part of his investigations into functional equations arising in physical and astronomical problems.^[2] Volterra's seminal contributions appeared in a series of papers published in 1896, where he systematically treated integral equations of the second kind with variable upper limits, focusing on their inversion and iterated kernels.^[3] These works laid the foundational framework for what would later be named after him, emphasizing equations that model systems where the present state depends on past history.^[3] Volterra's 1896 work stemmed from studies in mathematical physics and astronomy during his tenure at the University of Pisa and Rome. His framework for integral equations with memory effects was later applied to modeling hereditary phenomena, such as viscoelasticity in materials (around 1909) and population dynamics in ecology (1926 predator-prey model).^[2]^[16] A key advancement came in 1908 with the doctoral thesis of Romanian mathematician Traian Lalescu, titled Sur les équations de Volterra, supervised by Émile Picard at the Sorbonne.^[17] Lalescu's work provided the first comprehensive theory for these equations, including the development of the resolvent kernel concept, which enabled explicit solutions through series expansions and established existence and uniqueness results under suitable conditions.^[17] In the early 20th century, further developments connected Volterra equations to ordinary differential equations through explorations of conversion techniques and singular kernels to bridge integrodifferential systems.^[18] These contributions expanded applications in boundary value problems and stability analysis.^[18] By the mid-20th century, Volterra integral equations had evolved into a distinct branch of analysis, with numerical methods emerging in the late 1950s–1960s to address practical computations.^[19] The theory was solidified in the late 20th century through comprehensive monographs, such as Volterra Integral and Functional Equations by Gustaf Gripenberg, Stig-Olof Londen, and Olof Staffans (1990), which unified linear and nonlinear aspects with modern functional analytic tools.^[20]

Classification and Properties

First-Kind Equations

Volterra integral equations of the first kind take the standard form

\int_a^t K(t,s) f(s) \, ds = g(t),

where f is the unknown function, K(t,s) is the kernel, and the integral is over the interval from a to t.^[21] This equation differs from those of the second kind by the absence of an explicit term involving f(t) outside the integral, which often results in ill-posedness; specifically, for non-degenerate continuous kernels, the range of the associated integral operator is non-closed, leading to a potential loss of information and sensitivity to perturbations in the right-hand side g(t).^[21]^[22] Existence and uniqueness of solutions require certain conditions on the data and kernel. Typically, g must be differentiable with g(a) = 0, and the kernel K should be continuous on the domain with K(t,t) \neq 0 for all t \in [a,b] to ensure the solution is unique in appropriate function spaces, such as L^2(a,b).^[21]^[23] Under these conditions, particularly when g is sufficiently smooth and the kernel is differentiable in its first argument, the first-kind equation can be converted to an equivalent second-kind form by differentiation:

f(t) = \frac{1}{K(t,t)} \left[ g'(t) - \int_a^t \frac{\partial K}{\partial t}(t,s) f(s) \, ds \right].

This transformation highlights the inherent relation to differentiation equations but preserves the triangular structure of the Volterra operator.^[21] A prominent example is Abel's integral equation of the first kind, a special case with a singular (weakly singular) kernel K(t,s) = \frac{(t-s)^{\alpha-1}}{\Gamma(\alpha)} for $0 < \alpha < 1:

\int_a^t \frac{(t-s)^{\alpha-1}}{\Gamma(\alpha)} f(s) \, ds = g(t).

This arises in applications like fractional calculus and tautochrone problems, where the singularity at s = t exacerbates the ill-posedness.^[21]^[24] Due to the ill-posed nature, solutions to first-kind Volterra equations are unstable in the presence of noisy data, necessitating regularization techniques—such as Tikhonov or Lavrent'ev methods—to achieve stable approximations in practical computations.^[21]

Second-Kind Equations

Second-kind Volterra integral equations are characterized by the presence of the unknown function both outside and inside the integral, making them generally more tractable than their first-kind counterparts. The standard form is given by

f(t) = g(t) + \lambda \int_a^t K(t,s) f(s) \, ds,

where f(t) is the unknown function, g(t) is a given forcing function, \lambda is a complex parameter, K(t,s) is the kernel, and the integration is over the interval from a to t with a \leq t \leq b.^[25] When g(t) = 0, the equation reduces to the homogeneous case

f(t) = \lambda \int_a^t K(t,s) f(s) \, ds,

which typically admits only the trivial solution f \equiv 0 under suitable conditions on the kernel.^[25] Vito Volterra introduced this classification in the late 19th century, distinguishing second-kind equations from first-kind forms and highlighting their solvability through iterative methods, treating them as limits of finite systems of linear equations.^[2] These equations possess favorable properties arising from the Volterra integral operator Vf(t) = \int_a^t K(t,s) f(s) \, ds, which is compact on spaces like the continuous functions C[a,b] when K is continuous.^[26] This compactness implies that the operator I - \lambda V leads to results analogous to the Fredholm alternative, but simplified by the triangular structure of the kernel (where K(t,s) = 0 for s > t), ensuring unique solvability for all \lambda \neq 0 without nontrivial homogeneous solutions, as the spectrum of V is \{[0](/page/0)\}.^[26] A prominent subclass consists of convolution-type equations, where the kernel depends only on the difference t - s:

f(t) = g(t) + \lambda \int_0^t k(t-s) f(s) \, ds.

Here, the solution can be expressed using the resolvent kernel r(u), satisfying r = \lambda k + \lambda k * r (with * denoting convolution), yielding f(t) = g(t) + \int_0^t r(t-s) g(s) \, ds.^[1] For example, with k(u) = e^{-u}, the resolvent is r(u) = 1, so the solution simplifies to f(t) = g(t) + \int_0^t g(s) \, ds.^[1] An illustrative case with an exponential kernel is the equation f(t) = 1 + \int_0^t e^{-(t-s)} f(s) \, ds, solvable via the Neumann series expansion f(t) = \sum_{n=0}^\infty V^n(1)(t), where V^0(1) = 1 and V^{n+1}(1)(t) = \int_0^t e^{-(t-s)} V^n(1)(s) \, ds. This series converges to the explicit solution f(t) = t + 1, demonstrating the iterative solvability emphasized by Volterra.^[27]

Analytical Solutions

Resolvent Kernel Method

The resolvent kernel method provides an analytical approach to solving linear Volterra integral equations of the second kind, expressed as u(t) = f(t) + \lambda \int_a^t K(t, s) u(s) \, ds, where \lambda is a complex parameter, K(t, s) is the kernel, and f(t) is the forcing function. The resolvent kernel R(t, s; \lambda) is defined such that the unique solution can be written in the variation-of-constants form u(t) = f(t) + \lambda \int_a^t R(t, s; \lambda) f(s) \, ds. This kernel satisfies the resolvent equation R(t, s; \lambda) = K(t, s) + \lambda \int_s^t K(t, v) R(v, s; \lambda) \, dv, or equivalently the symmetric form R(t, s; \lambda) = K(t, s) + \lambda \int_s^t R(t, v; \lambda) K(v, s) \, dv.^[28] The resolvent kernel is derived using the Neumann series expansion, where the iterated kernels are defined recursively by K_1(t, s) = K(t, s) and K_n(t, s) = \lambda \int_s^t K(t, v) K_{n-1}(v, s) \, dv for n \geq 2. Under suitable conditions on the kernel, such as |K(t, s)| \leq M for some constant M > 0 on the domain D = \{(t, s) \mid a \leq s \leq t \leq b\}, the series R(t, s; \lambda) = \sum_{n=1}^\infty \lambda^{n-1} K_n(t, s) converges absolutely and uniformly on D when |\lambda| < 1/M(b-a). This convergence ensures the existence and uniqueness of the continuous solution u(t) for continuous f(t) and K(t, s), as guaranteed by the contraction mapping principle or Gronwall's inequality.^[28] For convolution-type kernels, where K(t, s) = k(t - s), the resolvent kernel can be computed explicitly using Laplace transforms. The Laplace transform of the resolvent satisfies \hat{R}(p; \lambda) = \frac{\hat{k}(p)}{1 - \lambda \hat{k}(p)}, where \hat{k}(p) and \hat{R}(p; \lambda) denote the transforms with respect to the first argument, assuming the transforms exist and the denominator does not vanish in the region of convergence. The inverse Laplace transform then yields R(t, s; \lambda). This approach simplifies computations for specific kernels like exponential or polynomial forms.^[29] The method extends to Volterra equations of the first kind, f(t) = \lambda \int_a^t K(t, s) u(s) \, ds, by differentiating both sides to obtain a second-kind equation f'(t) = \lambda K(t, t) u(t) + \lambda \int_a^t \frac{\partial}{\partial t} K(t, s) u(s) \, ds, assuming f and K are differentiable and K(t, t) \neq 0. The solution is then found using the resolvent of the derived second-kind equation, but care must be taken with singularities arising if K(t, t) = 0, which may require additional regularization or alternative techniques to ensure well-posedness.^[30]

Laplace Transform Approach

The Laplace transform approach provides an analytical method for solving Volterra integral equations, particularly those of convolution type, by converting the integral equation into an algebraic equation in the transform domain. For a second-kind equation of the form f(t) = g(t) + \lambda \int_0^t k(t - s) f(s) \, ds, applying the Laplace transform utilizes the convolution theorem, yielding \hat{f}(p) = \frac{\hat{g}(p)}{1 - \lambda \hat{k}(p)}, where \hat{\cdot}(p) denotes the Laplace transform with respect to the variable p > 0. The solution f(t) is recovered via the inverse Laplace transform \mathcal{L}^{-1} \{ \hat{f}(p) \}. This transformation is effective when the transforms of g(t) and k(t) exist and the denominator does not vanish in the region of convergence.^[31] The inversion step often results in closed-form expressions, especially for kernels leading to rational functions in the p-domain. For polynomial kernels, solutions frequently involve exponentials or hyperbolic functions, as the poles of \hat{f}(p) determine the form through partial fraction decomposition. Consider the convolution-type equation f(t) = 1 + \int_0^t (t - s) f(s) \, ds; its Laplace transform gives \hat{f}(p) = \frac{p}{p^2 - 1}, with inverse f(t) = \cosh t. Even for non-convolution cases like f(t) = 1 + \int_0^t s f(s) \, ds, the Laplace transform leads to a differential equation in the p-domain, \frac{d}{dp} \hat{f}(p) + p \hat{f}(p) = 1, whose solution inverts to f(t) = e^{t^2/2}. These examples illustrate how the method yields explicit solutions when the inverse transform is tractable.^[31]^[32] For first-kind Volterra equations of convolution type, g(t) = \lambda \int_0^t k(t - s) f(s) \, ds, the transform simplifies to \hat{g}(p) = \lambda \hat{k}(p) \hat{f}(p), so \hat{f}(p) = \frac{\hat{g}(p)}{\lambda \hat{k}(p)}, assuming \hat{k}(p) \neq 0. In more general first-kind forms without strict convolution structure, the transform may introduce boundary terms, yielding \hat{g}(p) = p \hat{K}(p) \hat{f}(p) - boundary contributions from initial conditions or kernel behavior at boundaries, necessitating careful handling during inversion to ensure consistency, such as requiring g(0) = 0. Inversion proceeds similarly but can be ill-posed if \hat{k}(p) has zeros near the origin.^[32]^[33] This approach is limited to equations with convolution structure, where the kernel depends only on t - s; for general kernels K(t, s), additional techniques like expansion in Fourier series may be required to apply the transform effectively. The method assumes the functions are Laplace-transformable, with sufficient decay for convergence, and is most powerful for linear equations where exact inversion is possible.^[32]

Numerical Methods

Trapezoidal Rule

The trapezoidal rule provides a straightforward quadrature-based numerical method for approximating solutions to Volterra integral equations of the second kind, given by f(t) = g(t) + \int_a^t K(t,s) f(s) \, ds. To apply this method, the interval [a, b] is partitioned into n equal subintervals of width h = (b - a)/n, yielding grid points t_i = a + i h for i = 0, 1, \dots, n. The integral is then approximated using the composite trapezoidal rule:

\int_a^{t_i} K(t_i, s) f(s) \, ds \approx \frac{h}{2} \left[ K(t_i, t_0) f(t_0) + 2 \sum_{j=1}^{i-1} K(t_i, t_j) f(t_j) + K(t_i, t_i) f(t_i) \right].

This approximation leads to a system of linear equations that can be solved recursively for the values f(t_i).^[34]^[35] Substituting the trapezoidal approximation into the second-kind equation yields the recurrence relation

f(t_i) \approx \frac{g(t_i) + \frac{h}{2} K(t_i, t_0) f(t_0) + h \sum_{j=1}^{i-1} K(t_i, t_j) f(t_j)}{1 - \frac{h}{2} K(t_i, t_i)},

with the initial value f(t_0) = g(t_0). For i \geq 1, each f(t_i) depends only on previously computed values, forming a lower triangular system that ensures computational efficiency.^[34]^[35] Under the assumption of smooth kernel K(t,s) and forcing function g(t) (sufficiently differentiable on [a,b]), the global error of this method is O(h^2), consistent with the local truncation error of the trapezoidal quadrature. The lower triangular structure of the resulting system guarantees stability for the forward recursion, avoiding ill-conditioned matrices common in other integral equation discretizations. Implementation involves forward substitution, where values are computed sequentially from i=0 to i=n, with evaluations of K and g at grid points.^[34]^[35] A representative example is the equation f(t) = 1 + \int_0^t f(s) \, ds on [0,1], where K(t,s) = 1 and g(t) = 1, with exact solution f(t) = e^t. For n=2 (h=0.5), the approximations are f(0) = 1, f(0.5) \approx 1.667, and f(1) \approx 2.778, compared to exact values f(0.5) \approx 1.649 and f(1) = [e](/page/E!) \approx 2.718, illustrating the O(h^2) convergence as n increases.^[34]

Collocation Methods

Collocation methods approximate the solution of Volterra integral equations by representing it with piecewise polynomial functions on a partition of the integration interval. The domain [a, b] is divided into subintervals [t_{n-1}, t_n] for n = 1, \dots, N, where t_0 = a and t_N = b, with uniform step size h = (b - a)/N. On each subinterval, the approximate solution u_n(t) is a polynomial of degree m-1, constructed such that it interpolates previously computed values at t_{n-1} and satisfies the collocation equation at m-1 interior points t_{n j} = t_{n-1} + c_j h, where 0 < c_1 < \dots < c_{m-1} < 1 are fixed collocation parameters, often chosen as Gauss-Legendre nodes for optimal accuracy. For m=2 (linear polynomials, degree 1), there is one interior collocation point.^[36]^[37] For second-kind Volterra equations of the form u(t) = g(t) + \int_a^t K(t, s) u(s) , ds, the method leads to a system of algebraic equations. At each collocation point t_{n i}, the equation becomes u(t_{n i}) - g(t_{n i}) = \int_a^{t_{n i}} K(t_{n i}, s) P_n(s) , ds, where P_n(s) denotes the piecewise polynomial interpolant of the approximate solution up to t_n. This integral is evaluated using quadrature rules adapted to the polynomial degree, resulting in a recursive scheme solved sequentially over subintervals. The trapezoidal rule can be viewed as a special case of collocation methods using endpoint evaluation with linear approximation, though it relies on simpler quadrature without full interior collocation.^[36]^[37] These methods offer spectral accuracy for smooth solutions, achieving convergence rates that improve with polynomial degree, and can handle weakly singular kernels through graded meshes that refine near the singularity to restore optimal order. Error analysis relies on projection operators onto polynomial spaces, yielding uniform bounds of O(h^m) at collocation points, with potential superconvergence to O(h^{2m}) using appropriate nodes. For linear collocation (m=2), the global error is O(h^2); for quadratic (m=3), it reaches O(h^3).^[36]^[37]^[38] A representative example is the linear collocation method applied to the weakly singular second-kind Abel equation u(t) = 1 + \int_0^t (t - s)^{-1/2} u(s) , ds / \Gamma(1/2), whose exact solution is u(t) = \cosh(\sqrt{2 t}). Using piecewise linear polynomials on a uniform mesh, numerical experiments demonstrate a convergence rate of O(h^2) in the maximum norm, with errors reducing from approximately 10^{-2} for h=0.1 to 10^{-4} for h=0.01, confirming the theoretical bounds. Graded meshes further enhance accuracy near t=0, achieving near-optimal rates despite the singularity.^[36]^[37]^[38]

Applications

Ruin Theory

In the classical ruin problem within the Cramér-Lundberg model, the surplus of an insurance company at time t is given by U(t) = u + c t - \sum_{i=1}^{N(t)} X_i, where u \geq 0 is the initial surplus, c > 0 is the constant premium income rate, N(t) is a Poisson process with intensity \lambda > 0 counting the number of claims, and the X_i are independent and identically distributed claim sizes with distribution function F and finite mean \mu. The non-ruin (survival) probability \psi(u) = P(\inf_{t \geq 0} U(t) \geq 0 \mid U(0) = u) satisfies an integral equation obtained by conditioning on the occurrence time T and size X of the first claim.^[39] This conditioning yields the defective renewal equation

\psi(u) = \int_0^\infty \left[ \int_0^{u + c t} \psi(u + c t - x) \, dF(x) \right] \lambda e^{-\lambda t} \, dt,

where the inner integral accounts for cases where the first claim does not cause ruin (and the process restarts with adjusted surplus), while cases where the claim exceeds the accumulated surplus contribute 0 (implicitly). After a change of integration variable to express the equation in terms of the surplus level, this takes the form of a Volterra integral equation of the second kind.^[39] Solutions to this equation can be obtained using Laplace transforms or the resolvent kernel method, leading to explicit closed-form expressions when claims follow an exponential distribution F(x) = 1 - e^{-x/\mu} for x \geq 0. In this case, assuming the safety loading condition c > \lambda \mu, the non-ruin probability is

\psi(u) = 1 - \frac{\lambda \mu}{c} \exp\left( -\left( \frac{1}{\mu} - \frac{\lambda}{c} \right) u \right).

This formula highlights the exponential decay of the ruin risk as initial capital increases.^[40] For general claim distributions where closed forms are unavailable, numerical computation of \psi(u) relies on solving the associated Volterra equation via methods such as the trapezoidal rule or collocation techniques, which discretize the integral and iteratively approximate the solution over a grid of surplus levels. These approaches provide accurate estimates for practical ruin assessments in insurance.^[41] The Cramér-Lundberg model, which formalized this application of integral equations to ruin probabilities, originated from Filip Lundberg's 1903 doctoral thesis and was rigorously developed by Harald Cramér in the 1930s.^[42]

Viscoelasticity Models

Viscoelastic materials combine viscous and elastic properties, resulting in time-dependent stress-strain responses that depend on the material's deformation history. These behaviors are commonly modeled using Volterra integral equations, which capture the hereditary nature of the material through convolution integrals. Specifically, the stress \sigma(t) at time t can be expressed as

\sigma(t) = \int_{-\infty}^t G(t - s) \frac{d\varepsilon(s)}{ds} \, ds,

where \varepsilon(t) is the strain, and G(t - s) is the relaxation modulus representing the material's stress response to a unit step strain.^[43] This formulation arises from the Boltzmann superposition principle, assuming linearity and causality in the material response.^[44] In standard linear viscoelastic models, the relaxation modulus G(t) defines the stress relaxation under constant strain, while the creep compliance J(t) describes the strain accumulation under constant stress. These functions are interrelated through a Volterra integral equation of the second kind, such as

G(t) = \delta(t) - \int_0^t \dot{J}(t - \tau) G(\tau) \, d\tau,

where \delta(t) is the Dirac delta function, and the dot denotes differentiation. This equation allows computation of one function from the other, facilitating the modeling of either relaxation or creep scenarios.^[43]^[44] Solution techniques for these models often leverage the Laplace transform to simplify the convolution structure. Applying the transform yields the frequency-domain relation

\hat{\sigma}(p) = p \hat{G}(p) \hat{\varepsilon}(p),

where hats denote Laplace transforms, enabling straightforward computation of frequency responses for dynamic loading or harmonic excitations. Inverse transforms or numerical inversion can then recover time-domain solutions.^[43] Representative examples include the Maxwell and Kelvin-Voigt models, which reduce to specific kernel forms in the Volterra framework. The Maxwell model, consisting of a spring and dashpot in series, features an exponential relaxation modulus G(t) = E e^{-t/\lambda}, where E is the modulus and \lambda is the relaxation time, leading to unbounded creep under sustained stress. In contrast, the Kelvin-Voigt model, with a spring and dashpot in parallel, has an instantaneous response combined with viscous damping, expressed as \sigma(t) = E \varepsilon(t) + \eta \dot{\varepsilon}(t), which integrates to a Volterra form with a delta function kernel component for relaxation.^[43] Hereditary integrals in these models often require numerical simulation, such as multistep methods or collocation schemes, to evaluate the convolution efficiently for complex loading histories.^[45] Extensions of these Volterra-based models apply to quasi-static engineering problems, such as the deflection of viscoelastic beams under slow loading. For instance, the governing equation for a Timoshenko beam reduces to a second-kind Volterra integral equation incorporating the relaxation kernel, allowing analysis of time-dependent deformations in structures like polymer composites or biological tissues.^[46]

Extensions

Nonlinear Variants

Nonlinear Volterra integral equations extend the linear case by incorporating nonlinearity in the kernel or the unknown function, leading to forms such as f(t) = g(t) + \int_a^t K(t, s, f(s)) \, ds, where the kernel K depends on f(s).^[47]^[48] Existence and uniqueness of solutions for these equations are established using fixed-point theorems in appropriate Banach spaces, such as Banach's contraction mapping principle under Lipschitz continuity conditions on the nonlinearity, or Schauder's fixed-point theorem for compact operators when growth conditions are satisfied.^[49]^[50] These results require the nonlinearity to satisfy boundedness or sublinear growth to ensure the operator maps a suitable ball into itself and is contractive or compact.^[51] For perturbations of linear Volterra equations, where the nonlinearity is small, successive approximations—starting from the linear solution and iteratively incorporating the nonlinear term—converge to the unique solution under suitable smallness assumptions on the perturbation parameter.^[52] This method leverages the resolvent kernel from the linear theory as an initial guess, providing a constructive path to solutions when direct fixed-point application is challenging.^[53] In population dynamics, quadratic nonlinearities appear in models derived from logistic growth, capturing intraspecific competition effects.^[54] These models allow qualitative analysis through substitution in some cases. Nonlinear variants pose theoretical challenges, including the potential for multiple solutions when Lipschitz conditions fail, and finite-time blow-up under superlinear growth in the kernel, where solutions become unbounded in finite time despite continuous data.^[55]^[56] Blow-up occurs when the nonlinearity amplifies growth uncontrollably, necessitating careful analysis of continuation properties beyond local existence intervals.^[57] Recent advances as of 2025 include studies on Ulam stabilities for nonlinear Volterra equations and new numerical schemes like Bernstein polynomial methods and general linear methods for solving them.^[58]^[59]

Stochastic Versions

Stochastic Volterra integral equations (SVIEs) extend the classical deterministic Volterra framework by incorporating random noise, typically through Itô integrals with respect to Brownian motion or more general semimartingales, to model systems with hereditary effects under uncertainty. These equations capture memory-dependent dynamics in stochastic environments, such as diffusion processes with time-lagged influences. Seminal early work on SVIEs with Itô integrals established foundational existence and uniqueness results for linear and nonlinear forms driven by continuous martingales and processes of bounded variation.^[60] A standard forward SVIE takes the form

X(t) = \phi(t) + \int_0^t \mu(t,s,X(s)) \, ds + \int_0^t \sigma(t,s,X(s^-)) \, dL(s),

where X evolves in a Hilbert space, \phi is a given initial function, \mu and \sigma are coefficient functions satisfying Lipschitz and linear growth conditions, and L is a Hilbert space-valued Lévy process. Under these assumptions, a unique adapted càdlàg solution exists, obtained by embedding the equation into a Hilbert space of functions and relating it to a boundary value problem for a first-order stochastic partial differential equation (SPDE). This connection facilitates analysis of limiting laws and stability properties, with applications in energy markets and epidemic modeling where past states influence future volatility. Pathwise comparison theorems hold for equations with monotone coefficients, ensuring that solutions preserve ordering from initial conditions almost surely.^[61] Backward stochastic Volterra integral equations (BSVIEs) provide a dual framework, solving backward in time with terminal conditions and non-local generators to handle anticipation or recursive utilities in stochastic settings. The general Type-II BSVIE is

Y(t) = \psi(t) + \int_t^T g(t,s,Y(s),Z(t,s),Z(s,t)) \, ds - \int_t^T Z(t,s) \, dW(s),

where \psi is the terminal process, g is the generator with Lipschitz continuity, and W is Brownian motion. Unique adapted M-solutions exist under mild measurability and growth conditions, with L^2-estimates bounding the solution norms by those of \psi and g. These solutions admit representations via partial differential equations (PDEs) driven by forward diffusions, enabling decoupling into forward-backward systems for control and pricing problems. For instance, Type-I BSVIEs (without the Z(s,t) term) link directly to nonlinear PDEs for utility maximization.^[62] Extensions of SVIEs include versions with jumps via Lévy-Itô decompositions or singular kernels, preserving well-posedness under relaxed regularity assumptions on coefficients. In linear-quadratic (LQ) control settings, SVIEs model controlled state processes with quadratic costs, yielding optimal feedback controls through coupled forward-backward systems and Riccati-like equations derived from spike variations. Necessary and sufficient conditions for optimality involve positive semidefiniteness of weighting operators, with applications in stochastic portfolio optimization and risk management where memory kernels represent delayed market impacts. Numerical schemes, such as Euler-Maruyama discretizations adapted for Volterra structures, achieve strong convergence rates of order 0.5 under Hölder continuity of solutions.^[63]^[64] Recent developments as of 2025 encompass singular backward SVIEs in infinite dimensions, improved ϑ-methods for numerical solutions, and well-posedness for equations with weakly singular kernels driven by fractional Brownian motion.^[65]^[66]

References

[1]
Volterra Integral Equation - an overview | ScienceDirect Topics
Volterra integral equations refer to a class of integral equations that involve a function defined by an integral with respect to a variable limit, ...
[2]
Vito Volterra - Biography - MacTutor - University of St Andrews
Vito Volterra published papers on partial differential equations, particularly the equation of cylindrical waves. His most famous work was done on integral ...
[3]
1896–1996: One hundred years of Volterra integral equations of the ...
We review Vito Volterra's seminal papers (on the inversion of definite integrals) of 1896, with regard to their mathematical results and within the context ...
[4]
The Scientific Work of Vito Volterra - jstor
But it was in 1896 that his systematic treatment of integral equations was first printed. Half a dozen papers of that year discuss iterated kernels, the resol-.
[5]
[PDF] Volterra Integral Equations on Time Scales
This article introduces the basic qualitative and basic quantitative theory of Volterra integral equations on time scales and thus may be considered as a ...
[6]
Volterra integral equations of the first kind and applications to linear ...
Jul 29, 2020 · An algebraic formula for the solution of a Volterra integral equa- tion of the first kind is given in the topological algebra of locally ...
[7]
Volterra Integral Equations
This book offers a comprehensive introduction to the theory of linear and nonlinear Volterra integral equations (VIEs), ranging from Volterra's fundamental ...
[8]
https://www.cambridge.org/core/books/volterra-integral-equations/7046C767D498B31F59FE600A2F1892DD
[9]
[PDF] A Novel Approach for Solving Volterra Integral Equations Involving ...
The paper presents an approximation method called local fractional variational iteration method (LFVIM) for solving the linear and nonlinear Volterra integral ...
[10]
[PDF] Approximation of Volterra and Fredholm Integral Equations by ...
Integral equations arise in various fields of science and engineering and find their applications in conformal mapping, Volterra's population growth model ...
[11]
[PDF] Analytical and Numerical Solutions of Volterra Integral Equation of ...
This equation has wide range of applications in physics and engineering such as potential theory,. Dirichlet problems, electrostatics, the particle transport ...
[12]
[PDF] CALIFORNIA STATE UNIVERSITY, NORTHRIDGE VOLTERRA ...
Integral equations arise in a wide variety of mathematical, scientific and engineering problems. In this chapter we will point out some places where Volterra.
[13]
[PDF] Linear Volterra Integral Equations
Dec 12, 2016 · This chapter presents an introduction to the history and the classical theory of lin- ear Volterra integral equations of the first and second ...
[14]
None
### Summary: Converting Initial Value Problems (IVPs) of ODEs to Volterra Integral Equations
[15]
[PDF] Traian Lalescu - an outstanding romanian mathematician - Aosr.ro
Traian Lalescu stated that ”Picard was my master”. In 1908 he obtained his PhD with the thesis ”Sur l'`equation de. Volterra”. In 1908-1909 he participated ...
[16]
Integral equations between theory and practice: the cases of Italy ...
Tricomi, Francesco. 1957. Integral equations. New York: Interscience. Page references here from the 1985 reprint by Dover, New York. Volterra, Vito. 1954 ...
[17]
[PDF] On the History of Numerical Methods - for Volterra Integral Equations
As far as the general integral equation (1.1) is concerned, one finds the approach used by VOLTERRA already in the thesis by LE ROUX in 1894 (published as [56] ...
[18]
Volterra Integral and Functional Equations
Chapter 1 - Introduction and Overview pp 1-34 Access Part I - Linear Theory Access Chapter 2 - Linear Convolution Integral Equations pp 35-75 Access
[19]
[PDF] A Survey of Regularization Methods for First-Kind Volterra Equations
Abstract. We survey continuous and discrete regularization methods for first-kind Volterra problems with continuous kernels. Classical regu-.
[20]
[PDF] Ill-Posed Problems and Integral Equations. - DTIC
The severeness or mildness of ill-posedness of integral equations of tile first kind Af = g depends upon properties of the kernel and the regularity of the data.
[21]
[PDF] On the Uniqueness of Solutions of First Kind Volterra Integral ...
One widely used assumption to obtain uniqueness of the solution function is that the kernel function does not vanish on the diagonal, i.e. k(x, x) 6= 0 for all ...
[22]
[PDF] Solutions of the Generalized Abel's Integral Equation using ...
Taking a(x) = 0, b(x) = x, k(x, t) = (x − t)−α, Equations (1)-(2) can be regarded as Volterra integral equations of first and second kind, respectively. 2.2.
[23]
Volterra Integral Equation of the Second Kind - Wolfram MathWorld
Volterra Integral Equation of the Second Kind. An integral equation of the form. phi(x)=f(x)+int_a^xK(x,. where K(x,t) is the integral kernel, f(x) is a ...
[24]
[PDF] 2. Bounded and Compact Operators - LSU Math
In terms of spectral theory we can formulate the last result as follows: A Volterra integral operator with continuous kernel has no spectral values different ...
[25]
[PDF] 1 Introduction 2 Solution of Volterra Equations - Faculty
The second example has a quadratic polyno- mial solution x(t) = t2 + 2t + 1 and the same decaying exponential kernel. The right hand side is then g(t) = e;t.
[26]
None
### Summary of Resolvent Kernel for Linear Volterra Integral Equations (Second Kind)
[27]
None
### Summary of Resolvent Kernel Using Laplace Transform for Convolution Type Volterra Equations
[28]
(PDF) Resolution of system of Volterra integral equations of the first ...
A solution method for various systems of integral equations of the first kind is presented in this paper. The method starts off by transforming the systems ...
[29]
Using the Laplace Transform to Solve the Volterra Integral Equation
Aug 25, 2025 · This paper presents study the Laplace transform of the convolution integral and use this concept to study the solution of the Volterra ...
[30]
[PDF] MODULE 2 UNIT 1 S2 VOLTERRA INTEGRAL EQUATIONS 1.0 ...
Laplace transformation serves as a powerful tool in the solving of integral equations. Convolution, the inverse of Laplace transformation, is necessary to ...
[31]
Laplace transforms for approximation of highly oscillatory Volterra ...
Apr 1, 2014 · This paper focuses on Laplace and inverse Laplace transforms for approximation of Volterra integral equations of the first kind with highly ...
[32]
Analytical and Numerical Methods for Volterra Equations
Presents an aspect of activity in integral equations methods for the solution of Volterra equations for those who need to solve real-world problems.
[33]
Computational Methods for Integral Equations
L. M. Delves, University of Liverpool, J. L. Mohamed, University of Liverpool. Publisher: Cambridge University Press. Online publication date: December 2009.
[34]
Collocation Methods for Volterra Integral Equations
$$217.00On the Divergence of Collocation Solutions in Smooth Piecewise Polynomial Spaces for Volterra Integral Equations. BIT Numerical Mathematics, Vol. 44, Issue. 4, ...
[35]
Collocation Methods for Volterra Integral Equations Review
Oct 16, 2025 · We consider both exact and discretized one-step and multistep collocation methods, and illustrate main convergence results, making some ...
[36]
Collocation Methods for Weakly Singular Second-kind Volterra ...
Abstract. Collocation type methods are studied for the numerical solution of the weakly singular Volterra integral equation of the second kind: f(t)=g(t)+∫
[37]
[PDF] ull7b.tex Week 7: pm, 7.11.2018 6. The ruin problem and the ...
Combining, we obtain the integral equation for the non-ruin probability φ(u): ... For the Cramér-Lundberg ... the now-standard Cramér-Lundberg model in insurance, ...
[38]
A survey of some recent results on Risk Theory
The exact formula for exponential claim amounts becomes both a bound and an asymptotic formula for the ruin probability: one has ψ(u) ≤ e−Ru for all u ...<|control11|><|separator|>
[39]
[PDF] Numerical Ultimate Ruin Probabilities under Interest Force
We derive a linear Volterra integral equation of the second kind and apply an order four Block-by-block method in conjuction with the Simpson rule to solve the ...
[40]
[PDF] Stochastic Processes in Non-Life Insurance (SkadeStok) 2023/2024 ...
We start by giving a brief overview of the Cramér-Lundberg model. The model was first introduced informally in the doctoral thesis of Filip Lundberg in 1903.
[41]
None
### Summary of Volterra Integral Equations in Linear Viscoelasticity
[42]
On the Volterra integral equation relating creep and relaxation
Apr 8, 2008 · For viscoelastic materials, when the response is linear, the dual linear Boltzmann causal integral equations are the appropriate model. The ...
[43]
On the Numerical Solution of Volterra Integral Equations Arising in ...
On the numerical solution of Volterra integral equations arising in linear viscoelasticity. Available. A. Wineman.
[44]
A locking-free finite element formulation for a non-uniform linear ...
Oct 1, 2021 · Then, we have the quasi-static version of the viscoelastic Timoshenko beam model, stated as a Volterra integral equation of second kind ...
[45]
Nonlinear Volterra Integral Equations - SpringerLink
It is well known that linear and nonlinear Volterra integral equations arise in many scientific fields such as the population dynamics, spread of epidemics, ...Missing: general | Show results with:general
[46]
Numerical solution for a generalized form of nonlinear cordial ...
In this article, we propose a numerical method for a general form of nonlinear cordial Volterra integral equations.
[47]
The existence and uniqueness of the solution for nonlinear ...
In this subsection, we prove the existence and uniqueness of a solution for nonlinear. Fredholm integral equations and nonlinear Volterra integral equations by ...Missing: form | Show results with:form
[48]
[PDF] Existence and approximate solutions of nonlinear integral equations
Oct 10, 2003 · The existence of such solutions is based on some well-known fixed point theorems in Banach spaces such as Schaefer fixed point theorem, Schauder.
[49]
An existence theorem for nonlinear functional Volterra integral ...
Using the method of Petryshyn's fixed point theorem in Banach algebra, we investigate the existence of solutions for functional integral equations.
[50]
(PDF) The Method of Successive Approximations (Neumann's ...
Aug 6, 2025 · The method of successive approximations (Neumann's series) is applied to solve linear and nonlinear Volterra integral equation of the second kind.Missing: perturbation | Show results with:perturbation
[51]
[PDF] The Successive Approximation Method for Solving Nonlinear ...
Sep 30, 2019 · In this paper, we will use the successive approximation method for solving. Fredholm integral equation of the second kind using Maple18.
[52]
Approximations of the nonlinear Volterra's population model by an ...
Aug 8, 2011 · In this paper, we apply the new homotopy perturbation method to solve the Volterra's model for population growth of a species in a closed ...
[53]
Analytical and Numerical Solutions of the Riccati Equation Using the ...
Aug 6, 2025 · This work presents new approximate analytical solutions for the Riccati equation (RE) resulting from the application of the method of ...
[54]
Conditions for blow-up of solutions of some nonlinear Volterra ...
In this paper we give necessary and sufficient conditions for blow-up of solutions for a particular class of nonlinear Volterra equations. We also give some ...Missing: multiple | Show results with:multiple
[55]
A system of nonlinear Volterra equations with blow-up solutions
Aug 7, 2025 · A pair of coupled nonlinear Volterra equations are examined for solutions that can have either global or blow-up behavior.Missing: multiple | Show results with:multiple
[56]
Blow-up and superexponential growth in superlinear Volterra ...
This paper concerns the finite-time blow-up and asymptotic behaviour of solutions to nonlinear Volterra integro-differential equations.Missing: multiple | Show results with:multiple
[57]
Comparison of stochastic Volterra equations - Project Euclid
Berger, M. and Mizel, V. (1980) Volterra equations with Ito├ integrals. J. Integral Equations, 2, 187±. 245, 319±337. Cochran, W.G., Lee, J.-S. and Potthoff ...Missing: Vinter | Show results with:Vinter
[58]
None
### Summary of Stochastic Volterra Integral Equations (SVIEs) from arXiv:1903.05045
[59]
None
### Definition of Backward Stochastic Volterra Integral Equations (BSVIEs)
[60]
[PDF] Linear quadratic control problems of stochastic Volterra integral ...
For the study of SVIEs with deterministic coefficients, a class of stochastic Fredholm−Volterra integral equations is introduced to replace conventional forward ...<|control11|><|separator|>
[61]
Numerical methods for stochastic Volterra integral equations with ...
In this paper we first establish the existence, uniqueness and Hölder continuity of the solution to stochastic Volterra integral equations (SVIEs) with weakly ...