Terms Of Dissipation Research Articles

NUMERICAL SIMULATION OF WAVE DAMPING INDUCED BY VEGETATION ON SOLITARY WAVE RUN-UP WITH MOMENTUM CONSERVATIVE SCHEME

In this paper, the damping effects of vegetation are investigated in shallow water model. Wave dissipation due to the inertia and drag force induced by vegetation is implemented based on Morison’s formulation. The dissipation term is added in momentum equation. The non-linear model is then solved numerically using momentum conservative scheme on staggered grid domain. To validate the numerical scheme, simulations of solitary wave propagation without vegetation are conducted and compared to the exact solution and laboratory experiment. Besides, several simulations of solitary wave with vegetation on sloping beach are also conducted with varying plant stem density. The numerical results reveal that vegetation significantly reduces the velocity of wave propagation. Additionally, an increase in vegetation stem density correlates with a decrease in solitary wave run-up.

Enhanced Retrospective Forecasting in Dissipative Dynamical Systems Using Transformer and Multi-Scale ESRGAN Models

Accurate characterization of complex dynamical systems is crucial for understanding their intrinsic behavior, and retrospective prediction provides a promising solution. However, traditional methods often fail to effectively predict dissipative terms, which are key in dissipative dynamical systems. This study introduces a deep learning method (DLM) that combines a Transformer model with a multiscale enhanced super-resolution generative adversarial network (MS-ESRGAN) to improve retrospective predictions by extracting implicit information from temporal evolution data. The Transformer excels at capturing past dynamics, while MS-ESRGAN refines the predicted fields, achieving resolutions on par with ground truth data. The effectiveness of the DLM is demonstrated using two canonical flow cases: forced isotropic turbulence and a transitional boundary layer. The model closely matches the velocity fields of the ground truth, with only minor deviations attributed to the nonlinearity of the governing equations and the inherent difficulty in resolving small-scale structures. In addition, the DLM has been applied to National Oceanic and Atmospheric Administration (NOAA) sea surface temperature (SST) data, demonstrating its practical utility for climate science. Despite the challenges associated with capturing small-scale structures, the data-driven DLM outperforms traditional numerical methods that either neglect the dissipative term or utilize negative dissipative coefficients, representing a significant advancement in retrospective predictions.

Direct numerical simulation of boundary layer transition induced by roughness element in a low Reynolds number compressor blade

The surface roughness of blades in low Reynolds number compressors significantly impacts the transition of the boundary layer and the profile loss. This paper employs direct numerical simulation (DNS) to investigate the boundary layer flow within a low Reynolds number compressor cascade, where an isolated roughness element is situated on the suction surface. A cylinder and a ramp are crafted and simulated by the embedded boundary method (EBM). The research centers on assessing the influence of varied roughness element shapes on the downstream velocity field, laminar separation bubble (LSB), and transition. DNS findings indicate that the cylinder's downstream shear and disturbances are relatively weaker, prompting separation-induced transition within its boundary layer, while the ramp's intense disturbances lead to bypass transition. The lift-up mechanism and the self-induced amplification of streamwise gradient of spanwise velocity are intrinsic to the flow instability during the bypass transition. Both cases can diminish total pressure loss, with the ramp yielding a relatively higher loss due to the contribution of a large-scale turbulent boundary layer to the total pressure loss. In terms of profile loss, the terms of mean kinetic energy dissipation and production associated with the suction-surface-normal gradient of mean streamwise velocity are predominant.

Hybrid approach to perfect and dissipative spin hydrodynamics

A hybrid framework of spin hydrodynamics is proposed that combines the results of kinetic theory for particles with spin 1/2 with the Israel-Stewart method of introducing nonequilibrium dynamics. The framework of kinetic theory is used to define the perfect-fluid description that conserves baryon number, energy, linear momentum and spin part of angular momentum. This leads to the entropy conservation although, in the presence of spin degrees of freedom, the perfect-fluid formalism includes extra terms whose structure is usually attributed to dissipation. The genuine dissipative terms appear from the condition of positive entropy production in nonequilibrium processes. They are responsible for the transfer between the spin and orbital parts of angular momentum, with the total angular momentum being conserved. Published by the American Physical Society 2024

SPH simulations of non-isothermal viscoplastic free-surface flows incorporating Herschel-Bulkley-Papanastasiou model

In this paper, an improved smoothed particle hydrodynamics (SPH) method is employed to accurately simulate non-isothermal viscoplastic free surface flows, wherein the viscoplastic behavior of the fluid is precisely captured through the incorporation of the Herschel-Bulkley-Papanastasiou constitutive model. To suppress the non-physical oscillation arising from the weakly compressible hypothesis within the pressure field, the density dissipation term is incorporated into the mass conservation equation. To address the tensile instability arising from the uneven distribution of particles, the particle shifting technique is incorporated as a solution. To enhance the precision and ensure numerical stability of the gradient operator, a kernel gradient correction algorithm is implemented. The improved SPH method is employed for numerically simulating the non-isothermal viscoplastic mixed convection, dam-break flow and droplet impacting the solid wall. The effectiveness of the improved SPH method in tackling the complexities of non-isothermal viscoplastic fluid is validated through a rigorous comparison of its outcomes with those derived from alternative numerical methodologies. The assessment of the numerical convergence of the improved SPH method is undertaken through the utilization of varying initial particle spacings. The numerical outcomes demonstrate that the improved SPH method adeptly and precisely delineates the heat transfer mechanisms, intricate rheological properties, as well as the dynamic variation characteristics of the free surface in non-isothermal viscoplastic free surface flows.

Decay character and the semilinear structurally damped σ-evolution equations with two dissipative terms

In this paper, we study the decay rates of the solution and its derivatives to the Cauchy problem for the linear damped σ-evolution equation with two dissipative dampingsutt+(−Δ)σu+(−Δ)δ1ut+(−Δ)δ2ut=0 that are expressed in terms of decay characters of the initial data. Furthermore, by using the decay rates of solutions to the linear part and the fixed-point method, we investigate the existence and decay rate of the global solution with small data of the corresponding semi-linear problems with the nonlinear term of power type |ut|p, for the supercritical case p>p⁎. For some special data spaces, the optimality of p⁎ is guaranteed by nonexistence results in the subcritical case p<p⁎ and the critical case p=p⁎.

Linear response theory for transport in non-Hermitian [formula omitted]-symmetric models

We analyzed the influence of dissipative non-Hermitian terms on electrical transport in one-dimensional non-Hermitian PT-symmetric models, which presents chiral symmetry. The combination of periodic modulation and non-Hermitian effects opens up exciting avenues for exploring new phases of matter and understanding of the role of non-Hermitian physics in topological systems. The effective Hermitian Hamiltonian has always a higher dimension than the corresponding non-Hermitian one. We analyzed the effect of periodic hopping modulation on Su-Schrieffer-Heeger chain, that exhibits the non-Hermiticity in the presence of an complex staggered potential, where this dissipative non-Hermitian extension modifies the features of the topological trivial phase and the topological nontrivial phase.

Enhancing improved advection upstream splitting method on triangular grids: A hybrid approach for improved stability and accuracy in compressible flow simulations

This paper introduces NAUSM+M+AUFS (New Improved Advection Upstream Splitting Method Plus Artificially Upstream Flux Vector Splitting), a novel hybrid computational scheme for simulating compressible flows on triangular grids. The AUSM+M (Improved Advection Upstream Splitting Method) method is enhanced through two key modifications to boost numerical stability and robustness in high Mach number and hypersonic flows. The first modification redefines the interfacial numerical sound velocity, reducing shock anomalies and improving shock-capturing by integrating velocity and characteristic sound speed parameters. The second modification addresses the insufficiency of the pressure flux dissipation term at supersonic speeds by introducing a formulation that increases dissipation proportionally to the Mach number, thereby enhancing performance in high-speed flows. These enhancements constitute the NAUSM+M method. The NAUSM+M+AUFS scheme combines the strengths of NAUSM+M and AUFS (Artificially Upstream Flux Vector Splitting) methods, particularly in overcoming the limitations of NAUSM+M in handling shock instabilities and the carbuncle phenomenon on structure triangular grids. A dynamic switching function adjusts the weighting between NAUSM+M and AUFS, optimizing accuracy and stability based on local flow conditions. Numerical tests demonstrate that NAUSM+M+AUFS significantly outperforms AUSM+M, NAUSM+M, and AUFS, effectively eliminating the carbuncle phenomenon and providing smooth shock wave contours. In steady flow analysis, the new hybrid method achieves convergence speeds comparable to AUFS and shows 15% to 45% superior convergence accelerating than AUSM+M, depending on the convergence rate. In addition, in steady flow analysis, the accuracy of NAUSM+M+AUFS is 46% better than that of AUFS. This approach represents a significant advancement, offering a robust, accurate, and efficient solution for high-speed aerodynamic simulations, with broad applicability across various compressible flow challenges.

Simulations of nerve signal propagation in axons

PurposeIn the literature, soliton solutions of the Heimburg–Jackson model have been proposed by Drab et al. (2022), but for the considered models, i.e. Eq.(1) and Eq.(2), the existence of solitons, the dispersion analysis and the pseudospectral method have been studied (Engelbrecht et al., 2006, 2018, 2020; Tamm et al., 2017, 2022; Peets et al., 2013). Therefore, the gap should be filled by this work.Design/methodology/approachWhen nonlinear terms, dissipative terms and forcing terms are ignored, the system (Eq.(2)) reduces to a single, sixth-order partial differential equation (Tamm et al., 2022). In this work, our aim is to propose analytical solutions in the explicit form via ansatz-based method. Therefore, the parameter effects in wave profile will be proposed clearly in figures.FindingsWhile progress has been made in signal propagation in nerves, thanks to many experimental studies and theoretical predictions over the last two centuries, the results obtained in this study may answer new questions that arise.Originality/valueIn the literature, the existence of solitons, dispersion analysis and pseudospectral method have been investigated for the Heimburg-Jackson model (Engelbrecht et al., 2006, 2018, 2020; Tamm et al., 2017, 2022; Peets et al., 2013), and this study fills the gap in soliton solutions. Additionally, when nonlinear terms, dissipative term and forcing terms are ignored, the system (Eq.(2)) reduced to a single equation that is sixth-order partial differential equation (Tamm et al., 2022). In this work, our aim is to propose analytical solutions in the explicit form via ansatz-based method. Therefore, the parameter effects in wave profile will be proposed clearly in figures.

Data-Guided Low-Reynolds-Number Corrections for Two-Equation Models

Abstract The baseline Launder–Spalding k−ε model cannot be integrated to the wall. This paper seeks to incorporate the entire law of the wall into the model while preserving the original k−ε framework structure. Our approach involves modifying the unclosed dissipation terms in the k and ε equations specifically within the wall layer according to direct numerical simulation (DNS) data. The resulting model effectively captures the mean flow characteristics in both the buffer layer and the logarithmic layer, resulting in robust predictions of skin friction for zero-pressure-gradient (ZPG) flat-plate boundary layers and plane channels. To further validate our formulation, we apply our model to boundary layers under varying pressure gradients, channels experiencing sudden deceleration, and flow over periodic hills, with highly favorable results. Although not the focus of this study, the methodology here applies equally to the k–ω formulation and yields improved predictions of the mean flow in the viscous sublayer and buffer layer.

Comparative study on predicting turbulent kinetic energy budget using high-order upwind scheme and non-dissipative central scheme

The accurate computation of different turbulent statistics poses different requirements on numerical methods. In this paper, we investigate the capabilities of two representative numerical schemes in predicting mean velocities, Reynolds stress and budget of turbulent kinetic energy (TKE) in low Mach number flows. With concerns on numerical order of accuracy, dissipation and dispersion properties, a high-order upwind scheme with relatively good dispersion and dissipation and a second-order non-dissipative central scheme with perfect dissipation but poor dispersion are adopted for this comparative study. By carrying out a series of numerical simulations including Taylor-Green vortex, turbulent channel flow at Reτ=180\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re_{\ au }=180$$\\end{document} and turbulent flow over periodic hill at Reb=10595\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re_{b}=10595$$\\end{document}, it can be obtained that although the high-order upwind scheme lacks perfection of dissipation in high wave number range, it still demonstrates superior predictive capability compared with the second-order non-dissipative central scheme, especially with relatively coarse grids. Finally, by taking the high-order upwind scheme as a suitable selection for turbulence simulation, the turbulent flow over a 30P30N multi-element airfoil is investigated as an application study. After briefly comparing the simulated profiles and spectrum with reference experimental results as validation, the budget of TKE is analyzed to locate the dominant flow structures and regions. It is found that the production and dissipation terms behave in a “monopole” pattern in the locations with strong shears and wakes. Whereas the advection and diffusion terms show an “inward” pattern and an “outward” pattern, which indicate the spatial transport of TKE between the center of the shear layer and nearby locations.

Investigating thermal conduction dynamics in ternary Carreau–Yasuda nanofluids under convective boundary conditions and exponential heat source/sink effects

The analysis of the Carreau–Yasuda nanofluid (CYNF) across a stretching surface has practical applications in several fields, such as heat exchange, engineering, and material science. Scholars may discover novel perspectives for increasing heat transfer performance, establishing advanced materials, and enhancing the efficiency of multiple engineering systems by studying the behavior of CYNF in this particular instance. The energy transfer through trihybrid CYNF flow with the effect of magnetic dipole across a stretching sheet is examined in the present study. The ternary hybrid nanofluid (THNF) has been prepared by the addition of ternary nanoparticles (NPs) in the water (50%) and ethylene glycol (50%). Titanium dioxide (TiO2), silicon dioxide (SiO2), and aluminum oxide (Al2O3) are used in base fluid. The fluid velocity and heat transfer are examined under the impact of Darcy–Forchheimer, chemical reaction, convective condition, activation energy, and exponential heat source. The fluid flow has been stated in form of velocity, energy, and concentration equations. The set of modeled equations is simplified to non-dimensional form of ordinary differential equations (ODEs) by using similarity substitutions. Numerically, the system of lowest order ODEs is calculated through the parametric continuation method. It has been noticed that the temperature field augments against the variation of viscous dissipation and heat source term. Furthermore, the magnetic dipole has significant impact to enhance the thermal curve of THNF whereas falloffs the velocity profile. It can be noticed that the energy propagation rate enhances with the rising numbers of ternary NPs, from 1.22% to 5.98% (nanofluid), 2.30% to 9.61% (hybrid nanofluid), and 3.23% to 11.12% (trihybrid nanofluid).

A high-efficiency Q-compensated pure-viscoacoustic reverse time migration for tilted transversely isotropic media

The attenuation and anisotropy characteristics of real earth media give rise to amplitude loss and phase dispersion during seismic wave propagation. To address these effects on seismic imaging, viscoacoustic anisotropic wave equations expressed by the fractional Laplacian have been derived. However, the huge computational expense associated with multiple Fast Fourier transforms for solving these wave equations makes them unsuitable for industrial applications, especially in three dimensions. Therefore, we first derived a cost-effective pure-viscoacoustic wave equation expressed by the memory-variable in tilted transversely isotropic (TTI) media, based on the standard linear solid model. The newly derived wave equation featuring decoupled amplitude dissipation and phase dispersion terms, can be easily solved using the finite-difference method (FDM). Computational efficiency analyses demonstrate that wavefields simulated by our newly derived wave equation are more efficient compared to the previous pure-viscoacoustic TTI wave equations. The decoupling characteristics of the phase dispersion and amplitude dissipation of the proposed wave equation are illustrated in numerical tests. Additionally, we extend the newly derived wave equation to implement Q-compensated reverse time migration (RTM) in attenuating TTI media. Synthetic examples and field data test demonstrate that the proposed Q-compensated TTI RTM effectively migrate the effects of anisotropy and attenuation, providing high-quality imaging results.

Viscous-Layer Compressibility Correction for Two-Equation Reynolds-Averaged Navier–Stokes Models

The baseline two-equation Reynolds-averaged Navier–Stokes (RANS) models include fluid density but lack calibration for compressible flows, making them inadequate for high Mach numbers. Various compressibility corrections have been proposed to integrate the van Driest transformation or the semilocal transformation, i.e., a given compressible law of the wall, into the RANS formulation. Prior work has focused mainly on the logarithmic layer, but upon evaluating these log-layer compressibility corrections, we find that they do not significantly improve skin friction estimates. To overcome this challenge, we develop viscous-layer compressibility corrections. We do that by altering the dissipation terms. These corrections conform the RANS model to the semilocal scaling, resulting in more accurate predictions of mean velocity and temperature in a posteriori tests. Although not the primary focus of this study, we also find that the baseline one-equation Spalart–Allmaras model, which was proposed before the concept of semilocal scaling, produces results consistent with the semilocal scaling in a posteriori tests without requiring any compressibility correction.

Stress and power as a response to harmonic excitation of a fractional anti‐Zener and Zener type viscoelastic body

AbstractThe stress as a response to strain prescribed as a harmonic excitation is examined in both transient and steady state regime for the viscoelastic body modeled by thermodynamically consistent fractional anti‐Zener and Zener models by the use of the Laplace transform method. Assuming strain as a sine function, the time evolution of power per unit volume, previously derived as a sum of time derivative of a conserved term, which represents the rate of change of stored energy, and a dissipative term, which represents dissipated power, is investigated when expressed through the relaxation modulus and creep compliance. Further, two forms of energy and two forms of dissipated power per unit volume are examined in order to see whether they coincide.

Investigation of irreversible losses in subsonic compressor cascade flow field: A study on viscous dissipation and temperature gradient term in entropy production rate

The entropy generation rate is a widely employed loss characterization approach in the design and optimization of fluid machinery. This approach is utilized to quantitatively determine various types of losses in the flow field, and it is considered that the temperature gradient term in the entropy generation rate can characterize irreversible losses. This research paper aims to address the question of whether the temperature gradient term in entropy production rate can represent the irreversible losses in a flow field. By employing the Reynolds-Averaged Navier-Stokes (RANS) numerical simulation method, the flow field of a subsonic compressor cascade, both with and without wall heating, was investigated. Statistical analysis was conducted to examine the contributions of viscous dissipation, the temperature gradient term in entropy production rate, wall shear stress work, and global power loss. The study reveals the distribution differences of viscous dissipation and the temperature gradient term in entropy production rate within the separated flow region. The research results demonstrate that, under wall heating conditions, viscous dissipation increased by approximately 45%, and the ratio of the temperature gradient term to global power loss increased from 4% to 125%, exceeding the mechanical energy loss in the flow field. The distribution of the temperature gradient term within the separated flow region cannot be explained by irreversible losses, and it exhibits significant discrepancies with the distribution of viscous dissipation losses. Therefore, the temperature gradient term in entropy production rate cannot represent the local irreversible losses, whereas viscous dissipation serves as an accurate measure of local irreversible losses. The research results provide theoretical support for accurately determining the irreversible losses in the flow field.

The effect of inhomogeneous terms for asymptotic profiles of solutions to second-order abstract evolution equations with time-dependent dissipation

The inhomogeneous problem of second-order evolution equations with a time-dependent dissipation term of the form u″+Au+b(t)u′=g(t) in a real Hilbert space H is considered, where A is a general nonnegative selfadjoint operator in H. The result is the explicit description of asymptotic profile of solutions of such a problem when the dissipation b provides a diffusive structure with a suitable restriction. The proof is a basic energy method with a scope of well-behaved quantity for the diffusive structure.

Heat transfer analysis of nanofluid flow with entropy generation in a corrugated heat exchanger channel partially filled with porous medium

AbstractHeat exchanger research is mainly exploited to develop and optimize new engineering systems with high thermal efficiency. Passive methods based on nanofluids, fins, wavy walls, and the porous medium are the most attractive ways to achieve this goal. This investigation focuses on heat transfer and entropy production in a nanofluid laminar flow inside a plate corrugated channel (PCC). The channel geometry comprises three sections, partially filled with a porous layer located at the intermediate corrugate channel section. The physical modeling is based on the laminar, two‐dimensional Darcy–Brinkman–Forchheimer formulation for nanofluid flow and the local thermal equilibrium model for the heat equation, including the viscous dissipation term. Numerical solutions were obtained using ANSYS Fluent software based on the finite volume technique and the appropriate meshed geometries. The numerical results are validated with theoretical, numerical, and experimental studies. The simulations are performed for CuO–water nanofluid and AISI 304 porous medium. The coupled effects of porous layer thickness (δ), Reynolds number (Re), and nanoparticle fraction (φ) on velocity, streamlines, isotherm contours, Nusselt number (Nu), and entropy generation (S) are analyzed and illustrated. The simulation results demonstrate that heat transfer enhancement in clear PCC can be achieved using a porous layer insert. For the porous thickness range of [0.1–0.6], the corresponding range of average Nusselt number increase is [35.7%–176.9%], and the average entropy generation is [105.4%–771.9%]. The effect of the Reynolds number is more important in a porous duct than in a clear one. For δ = 0.4 and φ = 5%, the increase of Re in the range of [200–500] induces an increase in average Nusselt number in the range of [80.9%–108.4%] and average entropy in [222.9%–309.1%] comparatively to clear PCC. The effect of φ is practically the same for porous and clear channels. For φ = 5%, the increase on average Nu is about 9%, and entropy generation is 5%. Accordingly, important improvements in heat transfer in PCC can be achieved through the combined effect of flow Reynolds number and porous layer thickness.

Coupled Quintessence Inspired by Warm Inflation

We investigate a coupled quintessence cosmological model in which a dark-energy scalar field with an exponential potential interacts directly with a dark-matter fluid through a dissipative term inspired by warm inflation. The evolution equations of this model give rise to a three-dimensional dynamical system for which a thorough qualitative analysis is performed for all values of the relevant parameters. We find that the model is able to replicate the observed sequence of late-time cosmological eras, namely, a long enough matter-dominated era followed by a present era of accelerated expansion. In situations where there is a significant transfer of energy from dark energy to dark matter, temporary scaling-type solutions may arise, but, asymptotically, all solutions are dominated by dark energy.

High-order accurate well-balanced energy stable finite difference schemes for multi-layer shallow water equations on fixed and adaptive moving meshes

This paper develops high-order accurate well-balanced (WB) energy stable (ES) finite difference schemes for multi-layer (the number of layers M⩾2) shallow water equations (SWEs) with non-flat bottom topography on both fixed and adaptive moving meshes, extending our previous work on the single-layer shallow water magnetohydrodynamics [25] and single-layer SWEs on adaptive moving meshes [58]. To obtain an energy inequality, the convexity of an energy function for an arbitrary M is proved by finding recurrence relations of the leading principal minors or the quadratic forms of the Hessian matrix of the energy function with respect to the conservative variables, which is more involved than the single-layer case due to the coupling between the layers in the energy function. An important ingredient in developing high-order semi-discrete ES schemes is the construction of a two-point energy conservative (EC) numerical flux. In pursuit of the WB property, a sufficient condition for such EC fluxes is given with compatible discretizations of the source terms similar to the single-layer case. It can be decoupled into M identities individually for each layer, making it convenient to construct a two-point EC flux for the multi-layer system. To suppress possible oscillations near discontinuities, WENO-based dissipation terms are added to the high-order WB EC fluxes, which gives semi-discrete high-order WB ES schemes. Fully-discrete schemes are obtained by employing high-order explicit strong stability preserving Runge-Kutta methods and proved to preserve the lake at rest. The schemes are further extended to moving meshes based on a modified energy function for a reformulated system, relying on the techniques proposed in [58]. Numerical experiments are conducted for some two- and three-layer cases to validate the high-order accuracy, WB and ES properties, and high efficiency of the schemes, with a suitable amount of dissipation chosen by estimating the maximal wave speed due to the lack of an analytical expression for the eigenstructure of the multi-layer system.