Blockchain Papers

Follow blockchain research across journals, conferences, and preprint repositories.

5 papersLast indexed Aug 31, 2026
Search papers

Paper index

5 results · page 1 of 1

Clear filters
Mar 23, 2026·Zenodo (CERN European Organization for Nuclear Research)
0 cites
Symmetry, Triadic Sparsity, and Global Regularity for Kida-Pelz Navier-Stokes Flows

Andrea Cavazzini

This paper proves unconditional global regularity with quantitative exponential decay for the three-dimensional incompressible Navier–Stokes equations on the periodic box, restricted to velocity fields invariant under the Kida–Pelz symmetry group of order 48, for viscosities above an explicit threshold. The entire proof reduces, through a chain of six independently verifiable steps, to a single integer arithmetic fact: 20,625 < 31,104. This is, to the author's knowledge, the first time a Navier–Stokes regularity result has been distilled to a verifiable inequality between two five-digit integers, with every constant computed exactly and no numerical approximation entering the argument at any stage. The Kida–Pelz flow and why it matters. The Kida–Pelz initial datum, introduced by Kida (1985) and studied extensively by Pelz (2001), has occupied a special place in the blow-up literature for decades. It was originally proposed as a candidate for finite-time singularity formation precisely because its high octahedral symmetry concentrates vortex stretching into a small number of interacting structures, producing some of the most intense enstrophy growth observed in direct numerical simulations. The fact that the same symmetry that was expected to promote blow-up turns out to prevent it is itself a significant finding: it demonstrates that vortex stretching intensity and blow-up potential are fundamentally different quantities, a distinction that is often blurred in heuristic discussions of turbulence. The proof architecture. The argument has a deliberately transparent two-layer structure separating analysis from arithmetic, so that each layer can be checked independently by specialists in different fields. The analytic layer establishes three quantitative inputs. First, a spectral gap: representation-theoretic analysis of the octahedral group acting on Fourier space shows that the first two shells of the Laplacian spectrum are entirely killed by symmetry, tripling the effective Poincaré constant from 1 to at least 3. This means the KP symmetry forces vorticity to reside at higher wavenumbers where viscous dissipation is three times stronger than for generic flows. Second, a triadic density bound: the GKP equivariance constrains the Fourier support so severely that the number of resonant triads contributing to the nonlinear stretching term is reduced by a factor involving the group order, yielding a geometric density bound of at most 2. Third, an exact initial enstrophy: the KP datum is monochromatic, with all Fourier modes sitting at a single shell of squared wavenumber 11, giving the exact rational value 33/4 for the initial enstrophy. No floating-point computation, truncation, or discretisation enters this calculation. These three inputs feed into a Bernoulli differential inequality for the enstrophy whose separatrix is computed in closed form. The arithmetic layer then verifies that the initial enstrophy lies below this separatrix, which reduces to the integer comparison 33 times 625 equals 20,625, which is less than 31,104 equals 4 times 7,776. The safety margin is 50.8 percent, meaning the result would survive even if the analytic constants were degraded by up to 20 percent. Bounded enstrophy then gives global existence via the standard H1-continuation criterion, and exponential decay in all Sobolev norms follows by a Gronwall bootstrap. Beyond the core result. The paper establishes several extensions that go beyond mere regularity. Exponential decay is proved not only for the enstrophy but for all Sobolev norms simultaneously, with explicit prefactors and rates. The decay is shown to hold in all Lebesgue spaces from L2 to L-infinity and for all derivative orders, meaning that every physically measurable quantity associated with the flow decays exponentially. The pressure decays at double the velocity rate, a consequence of the quadratic structure of the pressure Poisson equation. Time analyticity is established for all positive times, meaning the solution extends to a holomorphic function in a strip around the real time axis. A shell-by-shell energy spectrum analysis shows that higher Fourier shells decay faster, with rates proportional to the squared wavenumber — a quantitative version of the physical intuition that small-scale structures are dissipated more rapidly. A Reynolds number characterisation shows that the Bernoulli closure holds if and only if the KP Reynolds number is below approximately 235, giving a concrete, physically interpretable criterion. The stability result deserves particular emphasis: global regularity is shown to persist under small perturbations that need not respect the KP symmetry. This means the result is not a fragile artifact of exact symmetry but a robust property of a neighbourhood in function space around the KP datum. The self-frustration connection. This paper is designed as a companion to the author's monograph "Self-Frustration of Vortex Stretching and the Architecture of the Navier–Stokes Blow-Up Barrier" (Cavazzini, 2026), which identifies a twelve-link chain of structural obstructions to finite-time blow-up for general three-dimensional Navier–Stokes. Three of those twelve links have concrete, quantitative realisations in the Kida–Pelz class. The enhanced spectral gap is a realisation of Link 5 (the spectral gap threshold that governs alignment stability). The triadic density reduction is a realisation of Link 3 (the oscillation bound that controls the pressure Hessian for tube-like vorticity). The identically vanishing helicity — proved here as a consequence of the parity inversion in the octahedral group — is a realisation of Link 6 (the gap–alignment complementarity), because it eliminates the eigenframe injection mechanism entirely: with zero helicity budget, the pressure Hessian cannot rotate the strain eigenframe to sustain the dangerous compressive component identified in the companion paper as the sole variable separating regularity-compatible from blow-up-compatible configurations. When all three mechanisms act simultaneously, as enforced by the octahedral symmetry, the self-frustration chain that remains open for general flows closes completely and unconditionally. The arithmetic inequality 20,625 < 31,104 is the quantitative expression of this closure. This provides the first concrete validation of the self-frustration framework as a genuine regularity tool rather than merely a classification scheme: the structural architecture described in the companion monograph is not an abstract taxonomy but a machinery that produces theorems when supplied with sufficient quantitative input. The minimal symmetry result strengthens this connection further: the octahedral group of order 48 is proved to be the smallest finite subgroup of O(3) for which the Bernoulli method closes. This characterises the precise boundary between symmetry groups where the self-frustration mechanisms are strong enough to guarantee regularity and those where they are not, providing a sharp answer to the question of how much geometric structure is needed to resolve the regularity problem within this framework. Context within the broader landscape. The Navier–Stokes regularity problem has a long history of partial results exploiting symmetry, from the classical two-dimensional theory (where regularity is known unconditionally due to the absence of vortex stretching) to various axisymmetric and helical reductions. The present work differs from these in a fundamental respect: the Kida–Pelz flow is fully three-dimensional with active, sustained vortex stretching — the mechanism responsible for the supercritical character of the equations is present and operative, not eliminated by dimensional reduction. What the symmetry does is not remove the stretching but quantitatively constrain it, tilting the balance between stretching and dissipation in favour of dissipation by a computable margin. This is a qualitatively different use of symmetry from the classical approach, and it suggests that the boundary between regularity and potential blow-up may be more accessible than previously thought — not through eliminating the dangerous mechanism, but through measuring and constraining it. The paper also contributes to the broader programme of understanding which structural properties of the Navier–Stokes equations are responsible for regularity. The identification of three independent mechanisms (enhanced dissipation, triadic depletion, topological obstruction) that close the regularity chain when acting together, combined with the companion monograph's demonstration that these same mechanisms are present but quantitatively insufficient for general flows, suggests a precise research programme: strengthen the quantitative estimates on any one of the three mechanisms sufficiently to close the chain without symmetry. The open problems listed in the paper — removal of the viscosity threshold, full Gevrey bootstrap, exact spectral gap computation — are formulated with this programme in mind. Methodological note. Every result in the paper carries an explicit epistemic label. All constants are computed exactly as rational numbers or algebraic expressions. The paper makes no claim regarding the Clay Millennium Prize and explicitly discusses the four gaps separating the present result from the Prize requirements: symmetry restriction, periodic domain, viscosity threshold, and partial Gevrey bootstrap. The distance from each gap to a resolution is assessed individually, with the viscosity threshold identified as an artifact of the Bernoulli method rather than a physical phase transition. MSC 2020 Classification: 35Q30 (primary — Navier–Stokes equations); 76D03 (existence, uniqueness, and regularity for incompressible viscous fluids); 42B25 (maximal functions and Littlewood–Paley theory); 20C15 (ordinary representations and characters of finite groups); 35B65 (smoothness and regularity of solutions to PDE

Open access
2 source records
Navier-Stokes equation solutions
Fluid Dynamics and Turbulent Flows
Advanced Numerical Methods in Computational Mathematics
Original source
Feb 5, 2023·Journal of the Mechanics and Physics of Solids
157 cites
Neural networks meet hyperelasticity: A guide to enforcing physics

Lennart Linden, Dominik K. Klein, Karl A. Kalina, Jörg Brummund · 6 authors

In the present work, a hyperelastic constitutive model based on neural networks is proposed which fulfills all common constitutive conditions by construction, and in particular, is applicable to compressible material behavior. Using different sets of invariants as inputs, a hyperelastic potential is formulated as a convex neural network, thus fulfilling symmetry of the stress tensor, objectivity, material symmetry, polyconvexity, and thermodynamic consistency. In addition, a physically sensible stress behavior of the model is ensured by using analytical growth terms, as well as normalization terms which ensure the undeformed state to be stress free and with zero energy. In particular, polyconvex, invariant-based stress normalization terms are formulated for both isotropic and transversely isotropic material behavior. By fulfilling all of these conditions in an exact way, the proposed physics-augmented model combines a sound mechanical basis with the extraordinary flexibility that neural networks offer. Thus, it harmonizes the theory of hyperelasticity developed in the last decades with the up-to-date techniques of machine learning. Furthermore, the non-negativity of the hyperelastic neural network-based potentials is numerically examined by sampling the space of admissible deformations states, which, to the best of the authors' knowledge, is the only possibility for the considered nonlinear compressible models. For the isotropic neural network model, the sampling space required for that is reduced by analytical considerations. In addition, a proof for the non-negativity of the compressible Neo-Hooke potential is presented. The applicability of the model is demonstrated by calibrating it on data generated with analytical potentials, which is followed by an application of the model to finite element simulations. In addition, an adaption of the model to noisy data is shown and its [...]

Open access
2 source records
Elasticity and Material Modeling
Model Reduction and Neural Networks
Advanced Numerical Methods in Computational Mathematics
Original source
Jan 1, 2018·TUScholarShare (Temple University)
1 cites
Asynchronous Optimized Schwarz Methods for Partial Differential Equations in Rectangular Domains

José C. Garay

Asynchronous iterative algorithms are parallel iterative algorithms in which communications and iterations are not synchronized among processors. Thus, as soon as a processing unit finishes its own calculations, it starts the next cycle with the latest data received during a previous cycle, without waiting for any other processing unit to complete its own calculation. These algorithms increase the number of updates in some processors (as compared to the synchronous case) but suppress most idle times. This usually results in a reduction of the (execution) time to achieve convergence. Optimized Schwarz methods (OSM) are domain decomposition methods in which the transmission conditions between subdomains contain operators of the form \linebreak $\partial/\partial \nu +\Lambda$, where $\partial/\partial \nu$ is the outward normal derivative and $\Lambda$ is an optimized local approximation of the global Steklov-Poincar\'e operator. There is more than one family of transmission conditions that can be used for a given partial differential equation (e.g., the $OO0$ and $OO2$ families), each of these families containing a particular approximation of the Steklov-Poincar\'e operator. These transmission conditions have some parameters that are tuned to obtain a fast convergence rate. Optimized Schwarz methods are fast in terms of iteration count and can be implemented asynchronously. In this thesis we analyze the convergence behavior of the synchronous and asynchronous implementation of OSM applied to solve partial differential equations with a shifted Laplacian operator in bounded rectangular domains. We analyze two cases. In the first case we have a shift that can be either positive, negative or zero, a one-way domain decomposition and transmission conditions of the $OO2$ family. In the second case we have Poisson's equation, a domain decomposition with cross-points and $OO0$ transmission conditions. In both cases we reformulate the equations defining the problem into a fixed point iteration that is suitable for our analysis, then derive convergence proofs and analyze how the convergence rate varies with the number of subdomains, the amount of overlap, and the values of the parameters introduced in the transmission conditions. Additionally, we find the optimal values of the parameters and present some numerical experiments for the second case illustrating our theoretical results. To our knowledge this is the first time that a convergence analysis of optimized Schwarz is presented for bounded subdomains with multiple subdomains and arbitrary overlap. The analysis presented in this thesis also applies to problems with more general domains which can be decomposed as a union of rectangles.

Open access
Advanced Numerical Methods in Computational Mathematics
Matrix Theory and Algorithms
Differential Equations and Numerical Methods
Original source
Jan 1, 2013·SIAM Journal on Numerical Analysis
34 cites
Study of Full Implicit Petroleum Engineering Finite-Volume Scheme for Compressible Two-Phase Flow in Porous Media

Bilal Saad, Mazen Saad

An industrial scheme, to simulate the compressible two-phase flow in porous media, consists of a finite volume method together with a phase-by-phase upstream scheme. The implicit finite volume scheme satisfies industrial constraints of robustness since the proposed scheme discretizes the equations with gravity and capillary terms. We show that the proposed scheme satisfies the maximum principle for the saturation, a discrete-energy estimate on the pressures, and a function of the saturation that denotes capillary terms. These stability results allow us to derive the convergence of a subsequence to a weak solution of the continuous equations as the size of the discretization tends to zero. To our knowledge, this is the first convergence result of a finite volume scheme in the case of two-phase compressible flow in several space dimensions. The proof is given for the complete system when the density of each phase depends on its own pressure.

Advanced Numerical Methods in Computational Mathematics
Enhanced Oil Recovery Techniques
Computational Fluid Dynamics and Aerodynamics
Original source
May 4, 1976·Proceedings of the Royal Society of London A Mathematical and Physical Sciences
37 cites
On the free surface of a viscous fluid motion

D. H. Sattinger

Abstract We consider a container of fluid with a rod inserted in the centre. As the rod rotates the surface of the fluid forms a curved surface whose shape correctly balances the forces of gravity, internal stress, atmospheric pressure, and surface tension. The surface of the fluid is depressed in the neighbourhood of the rod in the case of a Newtonian fluid, but the fluid may climb along the rod in the case of non-Newtonian fluids. In recent work by D. D. Joseph and R. S. Fosdick this free surface problem has been treated quantitatively by virtue of a formal perturbation series in which the solution is developed in powers of the angular velocity of the rod. The purpose of the present paper is to give a rigorous proof of convergence in the case of a Newtonian fluid. The method here may possibly be applied to other free surface problems - for example, surface waves of a viscous fluid. The convergence proof provides, to our knowledge, the first rigorous existence theorem for a free surface problem in the theory of viscous fluids. The proof also raises a novel problem in the theory of elliptic systems of partial differential equations. By means of the implicit function theorem, the question of convergence is reduced to that of obtaining a priori estimates for an elliptic boundary value problem. That problem is formulated on a domain with a ridge where the fluid surface meets the rod. In addition, the type of boundary conditions prescribed differ on the free surface and on the rod. In general, the solution to such a mixed boundary value problem is not smooth at the ridge, even if the boundary data is smooth. The solution will be smooth, however, if the boundary data satisfy certain consistency conditions which make it compatible with the given set of partial differential equations. The consistency conditions in question here are a set of linear relations between various derivatives of the boundary and inhomogeneous data. It is shown that if one makes the assumption that the wetting angle of the surface is zero, the consistency conditions are invariant under the full nonlinear equations of the free surface problem. This makes it possible to consider the original problem on a smaller function space - namely the subclass of functions satisfying the appropriate consistency conditions - and in this subclass one can apply the implicit function theorem and obtain the required a priori estimates. The solution thus obtained is regular up to the ridge.

Navier-Stokes equation solutions
Elasticity and Material Modeling
Advanced Numerical Methods in Computational Mathematics
Original source