Vanishing viscosity solution to a 2×2 system of conservation laws with linear damping

Basic Information

• Paper ID: 2408.10661
• Title: Vanishing viscosity solution to a 2×2 system of conservation laws with linear damping
• Author: Kayyunnapara Divya Joseph (Indian Institute of Science Education and Research Pune)
• Classification: math.AP (Analysis of PDEs)
• Submission Date: August 2024, Latest Version June 2025
• Paper Link: https://arxiv.org/abs/2408.10661v2

Abstract

This paper investigates a 2×2 non-strictly hyperbolic system of conservation laws with linear damping terms, arising from the Eulerian droplet model of air-particle flows. The author first studies the regularized viscous system with variable viscosity coefficients and obtains weak asymptotic solutions for general initial value conditions as well as solutions in Colombeau algebras. By studying the vanishing viscosity limit, it is proven that the limit is a solution in the distributional sense. Furthermore, the large-time asymptotic behavior of the viscous system is investigated. Due to analytical complexity, this system has not been adequately studied, with existing work limited to Riemann-type initial data. The important contribution of this paper lies in handling general initial value data rather than only Riemann-type data.

Research Background and Motivation

Problem Background

1. Core Problem: Study the initial value problem for the following 2×2 conservation law system: $\frac{\partial u}{\partial t} + \frac{\partial}{\partial x}\left(\frac{u^2}{2}\right) + \alpha u = 0$$\frac{\partial v}{\partial t} + \frac{\partial}{\partial x}(uv) = 0$ where α≠0 is a constant, with initial conditions u(x,0)=u₀(x), v(x,0)=v₀(x).
2. Physical Significance:
   • When α>0, the system arises from the Eulerian droplet model of air-particle flows
   • α represents the drag coefficient between air and particles
   • v≥0 and u represent the volume fraction and velocity of particles, respectively
   • Solutions develop delta-shocks
3. Mathematical Challenges:
   • Classical solutions of nonlinear first-order evolution equations break down in finite time
   • Extending solutions involves products of distributions, and Schwartz's impossibility result shows there exists no algebra containing the space of distributions and allowing such products
   • The system has a repeated eigenvalue u (with multiplicity 2) with only a one-dimensional eigenspace, making it a non-strictly hyperbolic system
   • Standard strictly hyperbolic system theory does not apply

Research Significance

1. Theoretical Value:
   • Non-conservative systems involving products of distributions naturally appear in multi-physics problems
   • Requires developing new generalized solution theory frameworks
   • Issues of non-smooth coefficients in continuity equations (at most BV regularity)
2. Application Background:
   • When α=0, related to modeling density inhomogeneity evolution in large-scale cosmic structure formation
   • Dynamics models of stochastic production flows (e.g., product flow in factories or supply chains)
   • Zero-pressure gas dynamics systems
   • Viscous particle dynamics

Limitations of Existing Methods

1. Existing Research:
   • Richard de la Cruz (2020) studied only Riemann problems (step-type initial data)
   • Work by Albeverio and Shelkovich mainly focused on the α=0 case
   • For the general initial value problem with α≠0, analysis is extremely complex and lacks systematic study
2. Technical Difficulties:
   • Solution space for u is the bounded variation function space (BV)
   • Solution space for v is the bounded Borel measure space
   • Characteristic lines are not straight (when α≠0), with du/dt = -αu along characteristics
   • Lack of uniqueness theory (uniqueness of the v component remains an open problem)

Research Motivation

This paper aims to construct explicit solutions for the system (1.1)-(1.2) under general initial value conditions, establish solution existence in different weak solution frameworks, and use the vanishing viscosity method as a physical selection criterion.

Core Contributions

1. Constructed explicit regularized solutions under general initial conditions: Obtained explicit solution formulas for the viscous system with variable viscosity coefficient ϵe^(-αt) (Theorem 3.1)
2. Established multiple weak solution frameworks:
   • Weak asymptotic solutions
   • Solutions in Colombeau algebras
   • Solutions in the distributional sense (based on Volpert products)
3. Proved existence of vanishing viscosity limits: Demonstrated that as ϵ→0, regularized solutions converge to distributional solutions of the original system (Theorems 4.1 and 4.2)
4. Verified entropy conditions: Proved that the velocity component u satisfies the Lax entropy inequality u(x⁺,t) ≤ u(x⁻,t)
5. Studied large-time asymptotic behavior: Analyzed the asymptotic form of viscous system solutions as t→∞, distinguishing between α>0 and α<0 cases (Theorem 3.2)
6. Handled general initial data: First to address general initial value data for this system (rather than only Riemann-type), significantly extending existing results

Detailed Methodology

Task Definition

Study the initial value problem for a 2×2 non-strictly hyperbolic conservation law system:

• Input: Initial value functions u₀∈L^∞(ℝ), v₀∈L¹(ℝ) (or more general distributions)
• Output: Solutions (u,v) where u is a BV function and v is a Borel measure
• Constraints: Satisfy the weak form of the conservation law equations

Core Method: Vanishing Viscosity Approximation

1. Regularized System

Consider parabolic regularization with variable viscosity: $\frac{\partial u}{\partial t} + \frac{\partial}{\partial x}\left(\frac{u^2}{2}\right) + \alpha u = \epsilon e^{-\alpha t}\frac{\partial^2 u}{\partial x^2}$$\frac{\partial v}{\partial t} + \frac{\partial}{\partial x}(uv) = \epsilon e^{-\alpha t}\frac{\partial^2 v}{\partial x^2}$

Key Design: The viscosity coefficient is ϵe^(-αt) rather than constant ϵ, enabling the use of the Hopf-Cole transformation.

2. Hopf-Cole Transformation

Introduce the transformation τ(t) = (1-e^(-αt))/α, and set:

• Û(x,t) = ∫₀ˣ u(s)ds, V̂(x,t) = ∫₀ˣ v(s)ds
• (u^ϵ, v^ϵ) = (e^(-αt)Û_x, V̂_x)

Further apply the Hopf-Cole transformation: $Û = -2\epsilon \log S^ϵ, \quad V̂ = \frac{C^ϵ}{S^ϵ}$

This transforms the nonlinear system into linear heat equations: $e^{\alpha t}\frac{\partial C^ϵ}{\partial t} = \epsilon\frac{\partial^2 C^ϵ}{\partial x^2}, \quad e^{\alpha t}\frac{\partial S^ϵ}{\partial t} = \epsilon\frac{\partial^2 S^ϵ}{\partial x^2}$

3. Explicit Solution Formula

By solving the heat equations, explicit solutions are obtained: $u^ϵ(x,t) = e^{-\alpha t}\frac{\int_{-\infty}^{\infty}\frac{x-y}{\tau(t)}e^{-\frac{\theta(x,y,\tau(t))}{2\epsilon}}dy}{\int_{-\infty}^{\infty}e^{-\frac{\theta(x,y,\tau(t))}{2\epsilon}}dy}$

$V^ϵ(x,t) = \frac{\int_{-\infty}^{\infty}(\int_0^y v_0(s)ds)e^{-\frac{\theta(x,y,\tau(t))}{2\epsilon}}dy}{\int_{-\infty}^{\infty}e^{-\frac{\theta(x,y,\tau(t))}{2\epsilon}}dy}$

where: $\theta(x,y,\tau) = \frac{(x-y)^2}{2\tau(t)} + \int_0^y u_0(s)ds$

This is the action functional, playing a central role in variational problems.

Weak Solution Theory Framework

1. Weak Asymptotic Solutions (Definition 2.1)

A family of functions (u^ϵ, v^ϵ) is called a weak asymptotic solution if:

• For all φ∈C_c^∞, the residual of the equation is o(1) (as ϵ→0)
• Initial value conditions are satisfied in the distributional sense

2. Volpert Product Theory

For a BV function u and v=V_x (V∈BV), define the non-conservative product ūv:

• At continuity points: ū(x,t) = u(x,t)
• At jump points: $\bar{u}(x,t) = \int_0^1 u((1-\alpha)u(x^-,t) + \alpha u(x^+,t))d\alpha$

This provides the distributional solution definition for the second equation v_t + (ūv)_x = 0.

3. Colombeau Algebra Solutions (Definition 2.5)

In the generalized function algebra G(Ω):

• Moderate elements: ∥∂_t^j∂x^ℓv^ϵ∥{L^∞(K)} = O(ϵ^(-N))
• Null elements: ∥∂_t^j∂x^ℓv^ϵ∥{L^∞(K)} = O(ϵ^M) for all M>0

Solutions (u,v) satisfy the equation in the associated sense (≈).

Technical Innovations

1. Variable Viscosity Coefficient Design: The choice of ϵe^(-αt) allows the time-transformed system to reduce to the standard heat equation, which is the key technique for handling the damping term
2. Minimization Problem: The limit solution is determined by the minimizer y(x,τ) of the variational problem min_y θ(x,y,τ), embodying a generalization of the Lax-Oleinik formula
3. Unified Framework: Establishes solution existence in three different generalized function theory frameworks, demonstrating the universality of the method
4. Large-Time Behavior Analysis: Distinguishes between α>0 (τ→1/α) and α<0 (τ→∞) cases, providing different asymptotic forms

Experimental Setup

This is pure mathematical theory research with no numerical experiments or datasets. Main results are established through rigorous mathematical proofs.

Theoretical Verification Framework

1. Initial Value Assumptions:
   • u₀∈L^∞(ℝ)∩L¹(ℝ)
   • v₀ is a bounded variation function with v₀∈L¹(ℝ)
2. Regularization Method: Use Friedrichs mollifier η^ϵ
   • u₀^ϵ = u₀ * η^ϵ
   • v₀^ϵ = v₀ * η^ϵ
3. Verification Steps:
   • Prove existence and smoothness of regularized solutions
   • Establish uniform estimates (independent of ϵ)
   • Prove existence of limits
   • Verify that limits satisfy the original equation

Main Results

Theorem 3.1: Explicit Formula for Regularized Solutions

For u₀∈L^∞(ℝ), v₀∈L¹(ℝ), the pair (u^ϵ, v^ϵ) defined by formula (3.4) is the unique C^∞ solution of the viscous system (3.1).

Proof Outline:

• Use Hopf-Cole transformation to reduce the problem to linear heat equations
• Obtain explicit solutions via convolution formulas
• Verify that all derivatives exist and satisfy the equation

Theorem 3.2: Large-Time Asymptotic Behavior

Case 1 (α>0): As t→∞, $\lim_{t\to\infty}-\frac{1}{2\epsilon}e^{\alpha t}u^ϵ(x,t) = \frac{dG/dx}{G}$$\lim_{t\to\infty}v^ϵ(x,t) = \frac{d}{dx}\left(\frac{H(x)}{G(x)}\right)$

where G(x) and H(x) are defined by formula (3.16), with limits holding uniformly on compact sets.

Case 2 (α<0): Introducing the scaled variable ξ=x/√(2ϵτ(t)), different asymptotic forms (3.19) are obtained.

Physical Meaning:

• α>0: Damping drives the system toward steady state
• α<0: Anti-damping leads to diffusive behavior

Theorem 4.1: Weak Asymptotic Solutions and Generalized Solutions

1. (ũ^ϵ, ṽ^ϵ) is a weak asymptotic solution: Using mollified initial values, residuals converge uniformly to zero
2. Existence of Limits: There exists a unique minimizer y(x,τ(t)) almost everywhere such that $u(x,t) = e^{-\alpha t}\frac{x-y(x,\tau(t))}{\tau(t)}$$V(x,t) = \int_0^{y(x,\tau(t))}v_0(s)ds$
3. Convergence: V^ϵ→V pointwise almost everywhere, v^ϵ→v in the distributional sense

Key Proof Points:

• Apply Hopf's variational principle and Jensen's inequality
• Use Evans' limit theorem (4.7) to handle exponential terms in integrals
• Establish uniform boundedness to apply dominated convergence theorem

Theorem 4.2: Distributional Solutions and Entropy Conditions

When v₀ is a bounded variation function:

1. Distributional Solution: The limit (u,v) satisfies the equation in the distributional sense, where uv is understood as the Volpert product
2. Entropy Condition: u satisfies the Lax entropy inequality u(x⁺,t) ≤ u(x⁻,t)

Proof Points:

• Verify ∂_tV + u∂_xV = 0 at continuity points
• Verify Rankine-Hugoniot condition at jump points with measure zero
• Use monotonicity of y(x,τ) to prove entropy condition

Theorem 4.3: Solutions in Colombeau Algebras

Under appropriate estimate conditions (4.16), (u^ϵ) and (v^ϵ) belong to G(Ω) and satisfy equation (2.6) in the associated sense.

Estimates: For all j,k≥0, $\|\partial_t^j\partial_x^k u^ϵ\|_{L^∞(\Omega_T)} = O(\epsilon^{-2(j+k)})$

These estimates are established through probabilistic measure representations (4.21)-(4.22) and Leibniz rules.

Related Work

Conservation Law System Theory

1. Classical Theory:
   • Lax (1957): Strictly hyperbolic system theory
   • Glimm (1965): Existence of large solutions
   • Kruzkov (1970): Entropy solution uniqueness theory
2. Non-Strictly Hyperbolic Systems:
   • LeFloch (1990): Existence and uniqueness for non-strictly hyperbolic systems
   • The system in this paper has repeated eigenvalues; standard theory does not apply

Distribution Product Theory

1. Volpert Theory (1967):
   • Definition of non-conservative products for BV functions
   • Handle jump points through averaging stacking
   • Application to non-conservative systems
2. Colombeau Algebras (1984):
   • Generalized function algebras containing distribution spaces
   • Allow products and nonlinear operations on distributions
   • Dal Maso, LeFloch, Murat (1995): Weak stability of non-conservative products
3. Weak Asymptotic Methods:
   • Maslov (1980): Method origin
   • Danilov, Shelkovich et al.: Interactions of nonlinear waves

Related Physical Models

1. Zero-Pressure Gas Dynamics (α=0 case):
   • Shandarin & Zeldovich (1989): Large-scale cosmic structure
   • Gurbatov & Saichev (1993): Viscous solutions
   • Albeverio & Shelkovich (2006): Delta-shocks
2. Eulerian Droplet Model:
   • Keita & Bourgault (2019): Delta-shocks for Riemann problems
   • Richard de la Cruz (2020): Riemann problems for system (1.1)
3. Related Applications:
   • Armbruster et al. (2004, 2019): Supply chain dynamics models
   • Forestier-Coste et al. (2015): Macroscopic production models

Advantages of This Paper

1. General Initial Data: First to address non-Riemann-type general initial data
2. Unified Framework: Establishes solutions in multiple weak solution theories
3. Explicit Formulas: Provides explicit expressions for regularized and limit solutions
4. Large-Time Behavior: Complete analysis of asymptotic behavior

Conclusions and Discussion

Main Conclusions

1. Constructive Results: Successfully constructed solutions for system (1.1) under general initial data through the vanishing viscosity method
2. Multiple Framework Verification:
   • Weak asymptotic solution framework
   • Colombeau algebra framework
   • Distributional solution framework (based on Volpert products)
   • Three frameworks mutually verify, enhancing result reliability
3. Physical Reasonableness:
   • Satisfies Lax entropy conditions
   • Vanishing viscosity limit provides physical selection criterion
   • Large-time behavior aligns with physical intuition
4. Theoretical Contributions:
   • Extends non-strictly hyperbolic system theory
   • Demonstrates effective methods for handling distribution products
   • Provides analytical paradigm for similar systems

Limitations

1. Uniqueness Problem:
   • Uniqueness of u component guaranteed by Kruzkov theory
   • Uniqueness of v component remains an open problem
   • v is a measure satisfying linear equations with discontinuous coefficients; standard L¹ estimates do not apply
2. Initial Value Restrictions:
   • Requires u₀∈L^∞∩L¹, v₀∈L¹ and bounded variation
   • More general distributional initial data (e.g., containing delta functions) not discussed
3. Missing Numerical Verification:
   • Pure theoretical research without numerical simulations
   • Explicit formulas involve complex integrals; numerical computation may be difficult
4. Multi-Dimensional Extension:
   • Only one-dimensional case considered
   • Multi-dimensional characteristic line structure more complex
5. Other Physical Effects:
   • Pressure terms not considered
   • External forces or source terms not included

Future Directions

1. Uniqueness Theory:
   • Develop uniqueness theory for v component
   • Explore new comparison principles in measure spaces
2. Numerical Methods:
   • Design numerical schemes preserving delta-shock structure
   • Verify theoretical predictions of large-time behavior
3. Generalizations:
   • Multi-dimensional cases
   • Systems with pressure terms
   • Variable coefficient damping α=α(x,t)
4. Applications:
   • Quantitative analysis of specific physical problems
   • Physical meaning and determination of parameter α
5. Related Systems:
   • Other non-strictly hyperbolic systems
   • Coupled systems (e.g., Euler-Poisson systems)

In-Depth Evaluation

Strengths

1. Theoretical Depth:
   • Systematically applies multiple modern PDE theory tools
   • Clever application of Hopf-Cole transformation demonstrates deep expertise
   • Elegant and powerful combination of variational principles and limit theory
2. Novelty:
   • First to address general initial value problem for this system, filling important gap
   • Variable viscosity coefficient ϵe^(-αt) design is original
   • Novel perspective of unifying three weak solution frameworks
3. Rigor:
   • All theorem proofs are complete and rigorous
   • Technical details carefully handled (e.g., distinction between pointwise and distributional convergence)
   • Comprehensive classification of special points (jump points, continuity points, irregular points)
4. Completeness:
   • Complete analytical chain from regularization to limit
   • Supplementary large-time behavior analysis enhances theoretical completeness
   • Multiple solution concepts mutually verify
5. Writing Quality:
   • Clear structure, rigorous logic
   • Sufficient background introduction, comprehensive related work review
   • Standard mathematical expression, consistent notation

Weaknesses

1. Practical Limitations:
   • Explicit formula (3.4) involves complex integrals, difficult for actual computation
   • Lacks concrete examples or numerical verification
   • Insufficient quantitative analysis of application problems
2. Theoretical Gaps:
   • Unresolved v component uniqueness is obvious deficiency
   • Insufficient discussion of properties of minimizer y(x,τ) (e.g., regularity)
   • Convergence rates not provided (only convergence proven)
3. Technical Limitations:
   • Hopf-Cole transformation depends on special viscosity coefficient form, limited generalizability
   • Strong one-dimensional restriction; multi-dimensional extension unclear
   • Relatively high regularity requirements on initial data (BV)
4. Insufficient Comparative Analysis:
   • Lacks detailed comparison with Richard de la Cruz's Riemann problem results
   • Does not discuss how general initial value solutions reduce to Riemann solutions
   • Connection to α=0 case could be deeper
5. Physical Interpretation:
   • Insufficient physical interpretation of delta-shocks
   • Physical meaning of large-time behavior could be more detailed
   • Physical range and actual values of parameter α not discussed

Impact Assessment

1. Theoretical Contribution:
   • Important supplement to non-strictly hyperbolic system theory
   • Demonstrates effective approach to handling distribution product problems
   • Provides methodological reference for related system research
   • Expected Citations: High citation potential in conservation laws, hyperbolic systems, and generalized function theory
2. Method Value:
   • Refined application of vanishing viscosity method
   • Comprehensive use of multiple weak solution frameworks has exemplary significance
   • Generalization of variational methods and Hopf formula
3. Application Prospects:
   • Provides theoretical foundation for droplet models, viscous particle dynamics
   • May inspire numerical method design
   • Provides mathematical support for multi-phase flow problems in engineering
4. Limitations:
   • Highly specialized; audience mainly limited to PDE theory researchers
   • Lack of numerical implementation reduces practical impact
   • Unresolved uniqueness problem limits theoretical completeness

Applicable Scenarios

1. Direct Application:
   • Eulerian droplet model of air-particle flows (α>0)
   • Damped zero-pressure gas dynamics
   • Deterministic moment equations of supply chain dynamics models
2. Method Reference:
   • Other non-strictly hyperbolic conservation law systems
   • Non-conservative systems involving distribution products
   • Singular limit problems requiring vanishing viscosity method
3. Theoretical Reference:
   • Application research of BV function theory
   • Concrete application cases of Colombeau algebras
   • Large-time asymptotic behavior analysis
4. Non-Applicable Cases:
   • Multi-dimensional problems (require new methods)
   • Strictly hyperbolic systems (simpler methods available)
   • Engineering problems requiring fast numerical solutions

Reproducibility

1. Theoretical Reproducibility: ★★★★★
   • All proofs complete and detailed
   • Key lemmas and theorems clearly marked
   • Mathematical derivations verifiable
2. Computational Reproducibility: ★★★☆☆
   • Explicit formula (3.4) in principle implementable
   • But involves high-dimensional numerical integration and optimization
   • Requires specialized numerical techniques for exponential integrals
3. Recommendations:
   • Provide reference code for numerical computation
   • Give analytical solutions for special cases
   • Discuss numerical stability issues

Overall Evaluation

This is a high-quality mathematical analysis paper making substantial contributions to non-strictly hyperbolic conservation law system theory. The author successfully synthesizes multiple modern PDE theory tools (Hopf-Cole transformation, Volpert products, Colombeau algebras, weak asymptotic methods) to address a complex system with important physical background, solving the general initial value problem for the first time.

Main Highlights:

• High theoretical depth and technical difficulty
• Original construction of explicit formulas
• Multiple framework verification enhances result reliability

Main Regrets:

• v component uniqueness unresolved
• Lack of numerical verification
• Insufficient application discussion

Recommendation Index: ★★★★☆ (4.5/5)

This paper is suitable for PDE theory researchers, applied mathematicians, and scholars interested in mathematical modeling of multi-phase flows. For researchers needing to handle similar non-conservative systems or distribution product problems, this paper provides valuable methodological reference.

Selected References

This paper cites 29 important references, with key references including:

• 29 Volpert (1967): Foundational theory of BV spaces and quasilinear equations
• 10 Dal Maso, LeFloch, Murat (1995): Definition and weak stability of non-conservative products
• 6-8 Colombeau (1984-1994): Generalized function algebra theory
• 16 Hopf (1950): Original work on Hopf-Cole transformation
• 23 Richard de la Cruz (2020): Latest research on Riemann problems for this system
• 1 Albeverio & Shelkovich (2006): Delta-shock theory
• 20 Kruzkov (1970): Entropy solution uniqueness theory

These references form the theoretical foundation and research context of this paper.