Gravity-Induced Modulation of Negative Differential Thermal Resistance in Fluids

Basic Information

• Paper ID: 2510.08909
• Title: Gravity-Induced Modulation of Negative Differential Thermal Resistance in Fluids
• Authors: Qiyuan Zhang, Juncheng Guo, Juchang Zou, Rongxiang Luo
• Classification: cond-mat.stat-mech (Condensed Matter Physics - Statistical Mechanics)
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.08909

Abstract

This study investigates the effects of gravity on negative differential thermal resistance (NDTR) in fluids through multi-particle collision dynamics (MPC) modeling. In the integrable limit, the authors derive an analytical heat flux formula exhibiting the NDTR effect and demonstrate that when gravity is introduced along the thermodynamic force direction, the critical temperature difference required for NDTR occurrence can be substantially reduced. The research further shows that the thermal bath-induced NDTR mechanism, originally applicable only to weakly interacting systems, can operate in systems with stronger interactions under gravity, while maintaining robustness in binary mixtures. These findings provide new insights into heat transport and establish a theoretical foundation for designing fluid thermal devices that exploit the NDTR effect under gravity.

Research Background and Motivation

1. Core Problem

Negative differential thermal resistance (NDTR) refers to the phenomenon where heat flux decreases as the temperature difference between thermal baths increases. This phenomenon is crucial for understanding the fundamental mechanisms of heat transport and for realizing novel thermal management functionalities. While NDTR has been extensively studied in lattice models, research on fluid systems remains relatively limited.

2. Problem Significance

• Fundamental Physics: NDTR phenomena reveal the basic mechanisms of non-equilibrium heat transport
• Application Value: Can be utilized to design novel thermal control devices, such as thermal transistors and thermal logic gates
• Theoretical Development: Extends NDTR theory from lattice systems to fluid systems

3. Limitations of Existing Methods

• Thermal bath-induced NDTR mechanisms are only effective in weakly interacting systems
• Lack of systematic research on gravity's effects on NDTR
• Insufficient understanding of NDTR mechanisms in binary mixtures

4. Research Motivation

Earth's gravitational field can significantly affect the thermophysical properties of various fluid systems; therefore, understanding gravity's influence on heat transport and NDTR characteristics in these systems is essential.

Core Contributions

1. Theoretical Derivation: Derives analytical heat flux formulas exhibiting NDTR effects in the presence of gravity in the integrable limit
2. Critical Temperature Difference Reduction: Demonstrates that gravity along the thermodynamic force direction can significantly lower the critical temperature difference for NDTR occurrence
3. Extension to Stronger Interactions: Shows that gravity enables the NDTR mechanism to operate effectively in systems with stronger interactions
4. Robustness in Binary Mixtures: Validates the mechanism's robustness in binary fluid mixtures
5. Theoretical Foundation: Provides theoretical basis for designing fluid thermal devices exploiting NDTR effects under gravity

Methodology Details

Task Definition

Investigation of the modulation of negative differential thermal resistance in two-dimensional fluid systems by gravitational fields, specifically including:

• Inputs: Gravitational field strength g, temperature difference ΔT, interaction strength τ
• Outputs: Heat flux J, collision frequency f, density distribution ρ(h), temperature distribution T(h)
• Constraints: System operates under steady state, satisfying energy and momentum conservation

Model Architecture

1. Physical Model Setup

• System: N equal-mass point particles moving in a rectangular domain (L × H)
• Boundary Conditions: Upper and lower boundaries coupled to thermal baths at temperatures T_U and T_L respectively, with periodic boundary conditions in the x-direction
• Gravitational Field: g = (0, -g) along the y-direction

2. Multi-Particle Collision Dynamics (MPC)

System evolution consists of two stages:

• Streaming Stage: Particles move under gravity for time τ
• Collision Stage: Space is divided into a × a grid cells; particle velocities within each cell are rotated around the center-of-mass velocity by ±θ

Equations of motion:

ṙᵢ = vᵢ + gτ
ṗᵢ = mᵢv̇ᵢ = mᵢg

Velocity update rule:

vᵢ → V_c.m. + R̂_±θ(vᵢ - V_c.m.)

3. Heat Flux Formula Derivation

In the integrable limit (no particle interactions), the analytical expression for heat flux is:

J = N(1+d)/2 |k_B(T_U-T_L)|/|⟨t_U→L⟩+⟨t_L→U⟩|

where ⟨t_U→L⟩ and ⟨t_L→U⟩ represent the average times for particles to travel from the upper bath to the lower bath and vice versa:

⟨t_U→L⟩ = (1/g)√(πk_BT_U/2m){e^(mgH/k_BT_U)[1-Erf(√(mgH/k_BT_U))]-1} + √(2H/g)

⟨t_L→U⟩ = (1/g)√(πk_BT_L/2m){e^(mgH/k_BT_L)[1+Erf(√(mgH/k_BT_L))]-1} - √(2H/g)

Technical Innovations

1. Analytical Treatment of Gravitational Effects

• Accounts for gravity's influence on particle trajectories
• Handles the complex scenario of particles potentially colliding with the lower bath multiple times
• Establishes quantitative relationships between gravitational strength and critical temperature difference

2. Parameterization of Interaction Strength

• Uses time step τ as an adjustable parameter controlling interaction strength
• τ → ∞ corresponds to the non-interacting limit
• Smaller τ corresponds to stronger interactions

3. Extension to Binary Mixtures

• Considers binary fluids with different particle masses m and M
• Investigates effects of mass ratio and component proportion on NDTR

Experimental Setup

Numerical Simulation Parameters

• Basic Parameters: m_i = k_B = T = a = 1, θ = π/2
• System Size: L = 2, H = 10
• Particle Density: ρ = N/(LH) = 6
• Temperature Settings: T_U = T, T_L = T - ΔT

Evaluation Metrics

1. Heat Flux J: Energy exchange per unit time per unit area through the system
2. Collision Frequency f: Frequency of particle collisions with thermal baths
3. Density Distribution ρ(h): Particle number density distribution along the height direction
4. Temperature Distribution T(h): Temperature distribution along the height direction

Comparative Analysis

• Gravitational Strength: g = 0, 0.01, 0.05, 0.1, 0.2, etc.
• Interaction Strength: τ = ∞ (non-interacting), 10, 1, 0.1, etc.
• Binary Mixtures: Various mass ratios M and component proportions p

Experimental Results

Main Results

1. Gravity's Promotion of NDTR

• At g = 0: Critical temperature difference (ΔT)_cr = 0.75
• As g increases: (ΔT)_cr decreases significantly, NDTR region expands substantially
• Theory vs. Numerics: Analytical formulas perfectly match simulation results

2. Effects of Interaction Strength

• Without gravity (g = 0): NDTR disappears under strong interactions (τ = 0.1)
• With gravity (g = 0.2): NDTR persists even at τ = 0.1
• Mechanism: Gravity prevents low-velocity particles from reaching the hot end, maintaining the NDTR effect

3. Robustness in Binary Mixtures

• Mass Ratio Variation: As M varies from 1 to 12, NDTR region remains unchanged
• Component Proportion: As p varies from 0 to 1, mechanism remains robust
• Universality: Demonstrates the mechanism's universality in pure and mixed fluids

Physical Mechanism Analysis

1. Mechanism of Gravity Enhancement of NDTR

As gravity increases, particles reflected from the lower thermal bath require more collisions to reach the upper bath, resulting in:

• Increased propagation time ⟨t_L→U⟩
• Decreased collision frequency f
• NDTR appearance at smaller temperature differences

2. Density and Temperature Distributions

• Density Distribution: Follows Boltzmann distribution ρ(h) = ρ₀e^(-mgh/k_BT)
• Temperature Distribution: T(h) decreases with height due to kinetic energy conversion to gravitational potential energy

Reverse Gravity Configuration

When gravity opposes the thermodynamic force:

• Critical temperature difference (ΔT)_cr increases significantly
• NDTR region shrinks substantially
• Mechanism fails under strong interactions

Related Work

History of NDTR Research

1. 2004: Li et al. first predicted NDTR in nonlinear phononic lattices
2. Theoretical Development: Revealed fundamental NDTR characteristics and mechanisms based on lattice models
3. Quantum Systems: Extended NDTR research to quantum systems
4. Fluid Systems: Luo discovered thermal bath-induced NDTR in one-dimensional hard-point gases in 2019

Relationship to Related Work

• Extension: From one-dimensional to two-dimensional fluid systems
• Innovation: First systematic investigation of gravity's effects on NDTR
• Deepening: Reveals the physical mechanisms of gravity-modulated NDTR

Conclusions and Discussion

Main Conclusions

1. Gravity along the thermodynamic force direction significantly reduces the critical temperature difference for NDTR occurrence
2. Gravity enables the NDTR mechanism to operate effectively in strongly interacting systems
3. The mechanism maintains robustness in binary mixtures
4. Provides theoretical foundation for designing fluid thermal devices

Limitations

1. Model Simplification: MPC dynamics represents a simplified description of real fluids
2. Two-Dimensional Constraint: Primarily studies two-dimensional systems; three-dimensional effects require further investigation
3. Gravitational Strength Range: Limited range of gravitational strengths studied

Future Directions

1. Higher-Dimensional Extension: Investigate Rayleigh-Bénard convection in three-dimensional systems
2. Thermochemical Baths: Introduce coupled transport of particles and heat exchange
3. Experimental Verification: Explore possibilities for experimental validation
4. Device Design: Design fluid thermal transistors based on NDTR effects

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Provides complete analytical derivations and numerical verification
2. Physical Insights: Clearly explains the physical mechanisms of gravity-modulated NDTR
3. Systematic Approach: Comprehensive investigation from integrable cases to interacting systems to binary mixtures
4. Practical Value: Provides theoretical guidance for fluid thermal device design

Weaknesses

1. Experimental Absence: Lacks experimental verification, based solely on numerical simulations
2. Parameter Range: Study of certain parameters (e.g., gravitational strength) may be limited
3. Real-World Complexity: Actual fluid systems may exhibit complexity beyond model descriptions

Impact

1. Academic Contribution: Provides new perspectives for non-equilibrium heat transport theory
2. Application Prospects: Opens new directions for microfluidic thermal management technology
3. Interdisciplinary: Connects statistical physics, fluid mechanics, and heat transport fields

Applicable Scenarios

1. Microgravity Environments: Thermal management in space stations and other microgravity settings
2. Ground-Based Applications: Thermal device design considering gravitational effects
3. Fundamental Research: Theoretical investigations in non-equilibrium statistical physics

References

The paper cites 46 relevant references covering key works in NDTR theory development, MPC dynamics methods, gravity's effects on fluid properties, and other critical domains, providing a solid theoretical foundation for the research.

Overall Assessment: This is a high-quality theoretical physics paper that introduces gravity as a novel modulation factor in the study of NDTR, an important physical phenomenon. It not only provides rigorous theoretical analysis but also demonstrates promising application prospects. The paper features complete theoretical derivations and sufficient numerical verification, making significant contributions to the field of fluid heat transport.