Asymptotically entropy-conservative and kinetic-energy preserving numerical fluxes for compressible Euler equations

This paper proposes a hierarchy of numerical fluxes for the compressible flow equations which are kinetic-energy and pressure equilibrium preserving and asymptotically entropy conservative, i.e., they are able to arbitrarily reduce the numerical error on entropy production due to the spatial discretization. The fluxes are based on the use of the harmonic mean for internal energy and only use algebraic operations, making them less computationally expensive than the entropy-conserving fluxes based on the logarithmic mean. The use of the geometric mean is also explored and identified to be well-suited to reduce errors on entropy evolution. Results of numerical tests confirmed the theoretical predictions and the entropy-conserving capabilities of a selection of schemes have been compared.Comment: 9 pages, 4 figure

Asymptotically entropy conservative discretization of convective terms in compressible Euler equations

A new class of Asymptotically Entropy Conservative schemes is proposed for the numerical simulation of compressible (shock-free) turbulent flows. These schemes consist of a suitable spatial discretization of the convective terms in the Euler equations, which retains at the discrete level many important properties of the continuous formulation, resulting in enhanced reliability and robustness of the overall numerical method. In addition to the Kinetic Energy Preserving property, the formulation guarantees the preservation of pressure equilibrium in the case of uniform pressure and velocity distributions, and arbitrarily reduces the spurious production of entropy. The main feature of the proposed schemes is that, in contrast to existing Entropy Conservative schemes, which are based on the evaluation of costly transcendental functions, they are based on the specification of numerical fluxes involving only algebraic operations, resulting in an efficient and economical procedure. Numerical tests on a highly controlled one-dimensional problem, as well as on more realistic turbulent three-dimensional cases, are shown, together with a cost-efficiency study

Numerical treatment of the energy equation in compressible flows simulations

We analyze the conservation properties of various discretizations of the system of compressible Euler equations for shock-free flows, with special focus on the treatment of the energy equation and on the induced discrete equations for other thermodynamic quantities. The analysis is conducted both theoretically and numerically and considers two important factors characterizing the various formulations, namely the choice of the energy equation and the splitting used in the discretization of the convective terms. The energy equations analyzed are total and internal energy, total enthalpy, pressure, speed of sound and entropy. In all the cases examined the discretization of the convective terms is made with locally conservative and kinetic-energy preserving schemes. Some important relations between the various formulations are highlighted and the performances of the various schemes are assessed by considering two widely used test cases. Together with some popular formulations from the literature, also new and potentially useful ones are analyzed

Fast-projection methods for the incompressible navier–stokes equations

An analysis of existing and newly derived fast-projection methods for the numerical integration of incompressible Navier–Stokes equations is proposed. Fast-projection methods are based on the explicit time integration of the semi-discretized Navier–Stokes equations with a Runge–Kutta (RK) method, in which only one Pressure Poisson Equation is solved at each time step. The methods are based on a class of interpolation formulas for the pseudo-pressure computed inside the stages of the RK procedure to enforce the divergence-free constraint on the velocity field. The procedure is independent of the particular multi-stage method, and numerical tests are performed on some of the most commonly employed RK schemes. The proposed methodology includes, as special cases, some fast-projection schemes already presented in the literature. An order-of-accuracy analysis of the family of interpolations here presented reveals that the method generally has second-order accuracy, though it is able to attain third-order accuracy only for specific interpolation schemes. Applications to wall-bounded 2D (driven cavity) and 3D (turbulent channel flow) cases are presented to assess the performances of the schemes in more realistic configurations.Peer ReviewedPostprint (published version

Performance and error analysis of structure-preserving time-integration procedures for incompressible-flow simulations

The effects of kinetic-energy preservation errors due to Runge-Kutta (RK) temporal integrators have been analyzed for the case of large-eddy simulations of incompressible turbulent channel flow. Simulations have been run using the open-source solver Xcompact3D with an implicit spectral vanishing viscosity model and a variety of temporal Runge-Kutta integrators. Explicit pseudo-symplectic schemes, with improved energy preservation properties, have been compared to standard RK methods. The results show a marked decrease in the temporal error for higher-order pseudo-symplectic methods, and suggest that these schemes could be used to attain results comparable to traditional methods at a reduced computational cost

Performance and error analysis of structure-preserving time-integration procedures for incompressible-flow simulations

The effects of kinetic-energy preservation errors due to Runge-Kutta (RK) temporal integrators have been analyzed for the case of large-eddy simulations of incompressible turbulent channel flow. Simulations have been run using the open-source solver Xcompact3D with an implicit spectral vanishing viscosity model and a variety of temporal Runge-Kutta integrators. Explicit pseudo-symplectic schemes, with improved energy preservation properties, have been compared to standard RK methods. The results show a marked decrease in the temporal error for higher-order pseudo-symplectic methods, and suggest that these schemes could be used to attain results comparable to traditional methods at a reduced computational cost.Postprint (published version

Preventing postsurgical venous thromboembolism: pharmacological approaches

The use of antithrombotic drugs for the prevention of venous thromboembolism (VTE) in patients undergoing surgery is presently based on solid principles and high-level scientific evidence. This article reviews current strategies of pharmacological thromboprophylaxis. The level of VTE risk following surgery depends on a variety of factors that the surgeon should take into account, including the type of surgery and the presence of additional risk factors, such as elderly age and cancer. In patients undergoing minor general surgery, early mobilization is sufficient as prophylaxis, whereas in those undergoing major general surgery, thromboprophylaxis with low molecular weight heparin (LMWH), low-dose unfractionated heparin, or the pentasaccharide fondaparinux is recommended. Patients undergoing major orthopedic surgery have a particularly high risk of VTE, and routine thromboprophylaxis with LMWH, fondaparinux, or a vitamin K antagonist (international normalized ratio target: 2.0 to 3.0) is the standard of care in this group of patients. Recently, two new oral anticoagulants, rivaroxaban (a factor Xa inhibitor) and dabigatran etexilate (a direct thrombin inhibitor) have been licensed to be used for thromboprophylaxis after orthopedic surgery in Europe. Mechanical methods of thromboprophylaxis (compression stockings, intermittent pneumatic compression, vena cava filters), not discussed in detail in this review, should always be considered in patients at high thrombotic risk, in association with the pharmacolocical strategies, or in cases of contraindications to antocoagulants, as in patients or procedures at high risk of bleeding

Percutaneous Mechanical Circulatory Support Devices: Systems and Clinical Options

Cardiogenic shock (CS) still remains a leading cause of hospital death. The adoption of percutaneous ventricular assist devices (pVADs) as treatment of CS is an option which continues to rise. Several types of pVADs have been developed by time to provide full cardiac support with few related complications and easy implantation settings. pVADs are used to support the failing heart as a bridge to recovery, decision, durable device or heart transplantation. None of these devices adopted in the clinical practice is ideal for all patients. Disadvantages may be related to the risk of limb/arm ischaemia or cerebral stroke or haemolysis. The most important choice is to identify the best device for each patient depending on haemodynamics, clinical scenario and patient anatomical/pathological issues. This chapter discusses the current pVAD options to treat CS patients

A long-term chronology of Pinus pinea L. from Parco della Versiliana (Pietrasanta, Italy) derived from treefall induced by a windstorm on March 4th-5th, 2015

Abstract Pinewoods are distinctive environmental elements in the Mediterranean coastal area and have both natural and historical significance. From the evening of March 4th to the morning of March 5th, 2015, a severe and unusual windstorm occurred in the Tuscany region of central Italy with wind gusts over 120 km/h. The windstorm caused vast damage to the anthropic and natural environment and wounded numerous trees in the renowned pinewoods of Parco della Versiliana in the Tyrrhenian coastal area. The meteorological calamity provided the opportunity to i) date the onset of the artificial plantation of the present Italian stone pine (Pinus pinea L.) forest to the 1820s, ii) build a long-term tree-ring chronology of the Italian stone pines in the area and iii) analyze the climate-growth relationship of the Italian stone pine in the study area. The resulting Versiliana chronology was derived from 60 trees and spanned from 1828 to 2014 (187 years), representing one of the longest living Italian stone pine forests on the Italian Peninsula. Finally, the climate-growth analysis highlighted that at this site the latewood width is positively influenced by summer temperature, a peculiarity worthy of further investigations

Dataset on gait patterns in degenerative neurological diseases

We collected the gait parameters and lower limb joint kinematics of patients with three different types of primary degenerative neurological diseases: (i) cerebellar ataxia (19 patients), (ii) hereditary spastic paraparesis (26 patients), and (iii) Parkinson's disease (32 patients). Sixty-five gender-age matched healthy subjects were enrolled as control group. An optoelectronic motion analysis system was used to measure time-distance parameters and lower limb joint kinematics during gait in both patients and healthy controls