• 2-Methylfuran: A bio-derived octane booster for spark-ignition engines

      Sarathy, Mani; Shankar, Vijai; Tripathi, Rupali; Pitsch, Heinz; Sarathy, Mani (Elsevier BV, 2018-04-02)
      The efficiency of spark-ignition engines is limited by the phenomenon of knock, which is caused by auto-ignition of the fuel-air mixture ahead of the spark-initiated flame front. The resistance of a fuel to knock is quantified by its octane index; therefore, increasing the octane index of a spark-ignition engine fuel increases the efficiency of the respective engine. However, raising the octane index of gasoline increases the refining costs, as well as the energy consumption during production. The use of alternative fuels with synergistic blending effects presents an attractive option for improving octane index. In this work, the octane enhancing potential of 2-methylfuran (2-MF), a next-generation biofuel, has been examined and compared to other high-octane components (i.e., ethanol and toluene). A primary reference fuel with an octane index of 60 (PRF60) was chosen as the base fuel since it closely represents refinery naphtha streams, which are used as gasoline blend stocks. Initial screening of the fuels was done in an ignition quality tester (IQT). The PRF60/2-MF (80/20 v/v%) blend exhibited longer ignition delay times compared to PRF60/ethanol (80/20 v/v%) blend and PRF60/toluene (80/20 v/v%) blend, even though pure 2-MF is more reactive than both ethanol and toluene. The mixtures were also tested in a cooperative fuels research (CFR) engine under research octane number and motor octane number like conditions. The PRF60/2-MF blend again possesses a higher octane index than other blending components. A detailed chemical kinetic analysis was performed to understand the synergetic blending effect of 2-MF, using a well-validated PRF/2-MF kinetic model. Kinetic analysis revealed superior suppression of low-temperature chemistry with the addition of 2-MF. The results from simulations were further confirmed by homogeneous charge compression ignition engine experiments, which established its superior low-temperature heat release (LTHR) suppression compared to ethanol, resulting in better blending octane numbers. This work explores and provides a chemically sound explanation for the potential of 2-MF as an octane enhancer.
    • A 3D domain decomposition approach for the identification of spatially varying elastic material parameters

      Moussawi, Ali; Lubineau, Gilles; Xu, Jiangping; Pan, Bing (Wiley-Blackwell, 2015-02-24)
      Summary: The post-treatment of (3D) displacement fields for the identification of spatially varying elastic material parameters is a large inverse problem that remains out of reach for massive 3D structures. We explore here the potential of the constitutive compatibility method for tackling such an inverse problem, provided an appropriate domain decomposition technique is introduced. In the method described here, the statically admissible stress field that can be related through the known constitutive symmetry to the kinematic observations is sought through minimization of an objective function, which measures the violation of constitutive compatibility. After this stress reconstruction, the local material parameters are identified with the given kinematic observations using the constitutive equation. Here, we first adapt this method to solve 3D identification problems and then implement it within a domain decomposition framework which allows for reduced computational load when handling larger problems.
    • 3D imaging using X-Ray tomography and SEM combined FIB to study non isothermal creep damage of (111) oriented samples of γ / γ ′ nickel base single crystal superalloy MC2

      Jouiad, Mustapha; Ghighi, J.; Cormier, Jonathan; Ostoja-Kuczynski, Elisabeth; Lubineau, Gilles; Méndez, José (Trans Tech Publications, 2012-01)
      An unprecedented investigation consisting of the association of X-Ray tomography and Scanning Electron Microscopy combined with Focus Ion Beam (SEM-FIB) is conducted to perform a 3D reconstruction imaging. These techniques are applied to study the non-isothermal creep behavior of close (111) oriented samples of MC2 nickel base superalloys single crystal. The issue here is to develop a strategy to come out with the 3D rafting of γ' particles and its interaction whether with dislocation structures or/and with the preexisting voids. This characterization is uncommonly performed away from the conventional studied orientation [001] in order to feed the viscoplastic modeling leading to its improvement by taking into account the crystal anisotropy. The creep tests were performed at two different conditions: classical isothermal tests at 1050°C under 140 MPa and a non isothermal creep test consisting of one overheating at 1200°C and 30 seconds dwell time during the isothermal creep life. The X-Ray tomography shows a great deformation heterogeneity that is pronounced for the non-isothermal tested samples. This deformation localization seems to be linked to the preexisting voids. Nevertheless, for both tested samples, the voids coalescence is the precursor of the observed damage leading to failure. SEM-FIB investigation by means of slice and view technique gives 3D views of the rafted γ' particles and shows that γ corridors evolution seems to be the main creep rate controlling parameter. © 2012 Trans Tech Publications, Switzerland.
    • Absolute spectroscopy near 7.8 μm with a comb-locked extended-cavity quantum-cascade-laser

      Lamperti, Marco; Alsaif, Bidoor; Gatti, Davide; Fermann, Martin; Laporta, Paolo; Farooq, Aamir; Marangoni, Marco (Springer Nature, 2018-01-16)
      We report for the first time the frequency locking of an extended-cavity quantum-cascade-laser (EC-QCL) to a near-infrared frequency comb. The locked laser source is exploited to carry out molecular spectroscopy around 7.8 μm with a line-centre frequency combined uncertainty of ~63 kHz. The strength of the approach, in view of an accurate retrieval of line centre frequencies over a spectral range as large as 100 cm-1, is demonstrated on the P(40), P(18) and R(31) lines of the fundamental rovibrational band of N2O covering the centre and edges of the P and R branches. The spectrometer has the potential to be straightforwardly extended to other spectral ranges, till 12 μm, which is the current wavelength limit for commercial cw EC-QCLs.
    • Adaptive coupling between damage mechanics and peridynamics: a route for objective simulation of material degradation up to complete failure

      Han, Fei; Lubineau, Gilles; Azdoud, Yan (Elsevier BV, 2016-05-17)
      The objective (mesh-independent) simulation of evolving discontinuities, such as cracks, remains a challenge. Current techniques are highly complex or involve intractable computational costs, making simulations up to complete failure difficult. We propose a framework as a new route toward solving this problem that adaptively couples local-continuum damage mechanics with peridynamics to objectively simulate all the steps that lead to material failure: damage nucleation, crack formation and propagation. Local-continuum damage mechanics successfully describes the degradation related to dispersed microdefects before the formation of a macrocrack. However, when damage localizes, it suffers spurious mesh dependency, making the simulation of macrocracks challenging. On the other hand, the peridynamic theory is promising for the simulation of fractures, as it naturally allows discontinuities in the displacement field. Here, we present a hybrid local-continuum damage/peridynamic model. Local-continuum damage mechanics is used to describe “volume” damage before localization. Once localization is detected at a point, the remaining part of the energy is dissipated through an adaptive peridynamic model capable of the transition to a “surface” degradation, typically a crack. We believe that this framework, which actually mimics the real physical process of crack formation, is the first bridge between continuum damage theories and peridynamics. Two-dimensional numerical examples are used to illustrate that an objective simulation of material failure can be achieved by this method.
    • Adjustable static and dynamic actuation of clamped-guided beams using electrothermal axial loads

      Alcheikh, Nouha; Tella, Sherif Adekunle; Younis, Mohammad I. (Elsevier BV, 2018-02-14)
      The paper presents adjustable static and dynamic actuations of in-plane clamped-guided beams. The structures, of variable stiffness, can be used as highly tunable resonators and actuators. Axial loads are applied through electrothermal U-shaped and flexure beams actuators stacked near the edges of curved (arch) beams. The electrothermal actuators can be configurred in various ways to adjust as desired the mechanical stiffness of the structures; thereby controlling their deformation stroke as actuators and their operating resonance frequency as resonators. The experimental and finite element results demonstrate the flexibility of the designs in terms of static displacements and resonance frequencies of the first and second symmetric modes of the arches. The results show considerable increase in the resonance frequency and deflection of the microbeam upon changing end actuation conditions, which can be promising for low voltage actuation and tunable resonators applications, such as filters and memory devices. As case studies of potential device configurations of the proposed design, we demonstrate eight possibilities of achieving new static and dynamic behaviors, which produce various resonance frequencies and static displacement curves. The ability to actively shift the entire frequency response curve of a device is desirable for several applications to compensate for in-use anchor degradations and deformations. As an example, we experimentally demonstrate using the device as a resonant logic gate, with active resonance tuning, showing fundamental 2-bit logic functions, such as AND,XOR, and NOR.
    • Adsorption Characteristics of Water and Silica Gel System for Desalination Cycle

      Cevallos, Oscar R. (2012-07)
      An adsorbent suitable for adsorption desalination cycles is essentially characterized by a hydrophilic and porous structure with high surface area where water molecules are adsorbed via hydrogen bonding mechanism. Silica gel type A++ possesses the highest surface area and exhibits the highest equilibrium uptake from all the silica gels available in the market, therefore being suitable for water desalination cycles; where adsorbent’s adsorption characteristics and water vapor uptake capacity are key parameters in the compactness of the system; translated as feasibility of water desalination through adsorption technologies. The adsorption characteristics of water vapor onto silica gel type A++ over a temperature range of 30 oC to 60 oC are investigated in this research. This is done using water vapor adsorption analyzer utilizing a constant volume and variable pressure method, namely the Hydrosorb-1000 instrument by Quantachrome. The experimental uptake data is studied using numerous isotherm models, i. e. the Langmuir, Tóth, generalized Dubinin-Astakhov (D-A), Dubinin-Astakhov based on pore size distribution (PSD) and Dubinin-Serpinski (D-Se) isotherm for the whole pressure range, and for a pressure range below 10 kPa, proper for desalination cycles; isotherms type V of the International Union of Pure and Applied Chemistry (IUPAC) classification were exhibited. It is observed that the D-A based on PSD and the D-Se isotherm models describe the best fitting of the experimental uptake data for desalination cycles within a regression error of 2% and 6% respectively. All isotherm models, except the D-A based on PSD, have failed to describe the obtained experimental uptake data; an empirical isotherm model is proposed by observing the behavior of Tóth and D-A isotherm models. The new empirical model describes the water adsorption onto silica gel type A++ within a regression error of 3%. This will aid to describe the advantages of silica gel type A++ for the design of adsorption desalination processes where reducing capital cost and footprint area are highly important parameters to take into account.
    • Advancing predictive models for particulate formation in turbulent flames via massively parallel direct numerical simulations

      Bisetti, Fabrizio; Attili, Antonio; Pitsch, Heinz G. (The Royal Society, 2014-07-14)
      Combustion of fossil fuels is likely to continue for the near future due to the growing trends in energy consumption worldwide. The increase in efficiency and the reduction of pollutant emissions from combustion devices are pivotal to achieving meaningful levels of carbon abatement as part of the ongoing climate change efforts. Computational fluid dynamics featuring adequate combustion models will play an increasingly important role in the design of more efficient and cleaner industrial burners, internal combustion engines, and combustors for stationary power generation and aircraft propulsion. Today, turbulent combustion modelling is hindered severely by the lack of data that are accurate and sufficiently complete to assess and remedy model deficiencies effectively. In particular, the formation of pollutants is a complex, nonlinear and multi-scale process characterized by the interaction of molecular and turbulent mixing with a multitude of chemical reactions with disparate time scales. The use of direct numerical simulation (DNS) featuring a state of the art description of the underlying chemistry and physical processes has contributed greatly to combustion model development in recent years. In this paper, the analysis of the intricate evolution of soot formation in turbulent flames demonstrates how DNS databases are used to illuminate relevant physico-chemical mechanisms and to identify modelling needs. © 2014 The Author(s) Published by the Royal Society.
    • Aerodynamics characteristic of axisymmetric surface protuberance in supersonic regime

      Qamar, Adnan; Sanghi, Sanjeev (Informa UK Limited, 2012-01)
      The present work deals with the problem of an axi-symmetric surface protuberance mounted on a spherical nosed body of revolution. The numerical computations are carried out for laminar supersonic viscous flow for trapezoidal shape axi-symmetric protuberances. A free stream Mach number ranging from 3 to 8 in steps of 1 at a fixed free stream Reynolds number of 1.8x10(4) has been used in the present study. The steady solutions are obtained using a time marching approach. A newly developed Particle Velocity Upwinding (PVU) scheme has been used for the computation. The spatial flow pattern exhibits a strong bow shock in front of the hemispherical nose, which engulfs the entire base body. Near the protuberance, the fluid particle decelerates due to the adverse pressure created by the protuberance and thus the flow separates in front of the protuberance. This point of separation is found to be a function of Mach number and the protuberance shape. A low-pressure expansion region dominates the base region of the obstacle. The reattachment point for the base separation is also a function of Mach number. As the Mach number is increased the reattachment point shifts toward the protuberances base. A weak recompression shock is also seen in the base, which affects the separated zone behind the protuberance. The important design parameters such as skin friction, heat transfer, drag, and surface pressure coefficients are reported extensively.
    • All Organic Polymers Based Morphing Skin with Controllable Surface Texture

      Favero Bolson, Natanael (2018-05)
      Smart skins are integrating an increasing number of functionalities in order to improve the interaction between the systems they equip and their ambient environment. Here we have developed an electromechanical soft actuator with controlled surface texture due to applied thermal gradient via electrical voltage. The device was fabricated and integrated with optimized process parameters for a prepared heater element [doped PEDOT: PSS (poly-(3, 4 ethylenedioxythiophene): poly (styrene sulfonic acid))], a soft actuator (Ecoflex 00-50/ethanol) and overall packaging case [PDMS (polydimethylsiloxane)]. To study a potential application of the proposed smart skin, we analyze the fluid drag reduction in a texture controlled water flow unit. As a result, we obtained a reduction of approximately 14% in the skin drag friction coefficient during the actuation. We conclude that the proposed soft actuator device is a preferred option for a texture-controlled skin that reduces the skin drag friction coefficient.
    • Analysis of bifurcation behavior of a piecewise linear vibrator with electromagnetic coupling for energy harvesting applications

      El Aroudi, Abdelali; Ouakad, Hassen M.; Benadero, Luis; Younis, Mohammad I. (World Scientific Pub Co Pte Lt, 2014-05)
      Recently, nonlinearities have been shown to play an important role in increasing the extracted energy of vibration-based energy harvesting systems. In this paper, we study the dynamical behavior of a piecewise linear (PWL) spring-mass-damper system for vibration-based energy harvesting applications. First, we present a continuous time single degree of freedom PWL dynamical model of the system. Different configurations of the PWL model and their corresponding state-space regions are derived. Then, from this PWL model, extensive numerical simulations are carried out by computing time-domain waveforms, state-space trajectories and frequency responses under a deterministic harmonic excitation for different sets of system parameter values. Stability analysis is performed using Floquet theory combined with Filippov method, Poincaré map modeling and finite difference method (FDM). The Floquet multipliers are calculated using these three approaches and a good concordance is obtained among them. The performance of the system in terms of the harvested energy is studied by considering both purely harmonic excitation and a noisy vibrational source. A frequency-domain analysis shows that the harvested energy could be larger at low frequencies as compared to an equivalent linear system, in particular, for relatively low excitation intensities. This could be an advantage for potential use of this system in low frequency ambient vibrational-based energy harvesting applications. © 2014 World Scientific Publishing Company.
    • Analysis of Hydrogen/Air Turbulent Premixed Flames at Different Karlovitz Numbers Using Computational Singular Perturbation

      Manias, Dimitrios; Tingas, Alexandros-Efstathios; Hernandez Perez, Francisco E.; Im, Hong G.; Galassi, Riccardo Malpica; Ciottoli, Pietro Paolo; Valorani, Mauro (American Institute of Aeronautics and Astronautics, 2018-01-08)
      The dynamics and structure of two turbulent H2/air premixed flames, representative of the corrugated flamelet (Case 1) and thin reaction zone (Case 2) regimes, are analyzed and compared, using the computational singular perturbation (CSP) tools, by incorporating the tangential stretch rate (TSR) approach. First, the analysis is applied to a laminar premixed H2/air flame for reference. Then, a two-dimensional (2D) slice of Case 1 is studied at three time steps, followed by the comparison between two representative 2D slices of Case 1 and Case 2, respectively. Last, statistical analysis is performed on the full three-dimensional domain for the two cases. The dominant reaction and transport processes are identified for each case and the overall role of kinetics/transport is determined.
    • Analysis of interlaminar fracture toughness and damage mechanisms in composite laminates reinforced with sprayed multi-walled carbon nanotubes

      Almuhammadi, Khaled; Alfano, Marco; Yang, Yang; Lubineau, Gilles (Elsevier BV, 2014-01)
      The present work is focused on the nanoreinforcement of prepreg based carbon fiber composite laminates to improve delamination resistance. Functionalized multi-walled carbon nanotubes (MWCNTs) were dispersed over the interface between prepreg layers through solvent spraying and the resulting mode I interlaminar fracture toughness was determined. For comparison, baseline samples with neat prepregs were also prepared. Results indicate that the introduction of functionalized MWCNTs can favorably affect the interlaminar fracture toughness, and the associated mechanisms of failure have been investigated. The manufacturing procedures and the interfacial reinforcing mechanism were explored by analyzing (i) the wettability between CNTs-solvent solution and prepreg surface, (ii) CNTs dispersion and (iii) the fractured surfaces through high resolution scanning electron microscopy and Raman mapping. © 2013 Elsevier Ltd.
    • Analysis of the current–voltage curves and saturation currents in burner-stabilised premixed flames with detailed ion chemistry and transport models

      Belhi, Memdouh; Han, Jie; Casey, Tiernan A.; Chen, Jyh-Yuan; Im, Hong G.; Sarathy, S. Mani; Bisetti, Fabrizio (Informa UK Limited, 2018-05-22)
      Current-voltage, or i–V, curves are used in combustion to characterise the ionic structure of flames. The objective of this paper is to develop a detailed modelling framework for the quantitative prediction of the i–V curves in methane/air flames. Ion and electron transport coefficients were described using methods appropriate for charged species interactions. An ionic reaction mechanism involving cations, anions and free electrons was used, together with up-to-date rate coefficients and thermodynamic data. Because of the important role of neutral species in the ion production process, its prediction by the detailed AramcoMech 1.4 mechanism was optimised by using available experimental measurements. Model predictions were evaluated by comparing to i–V curves measured in atmospheric-pressure, premixed, burner-stabilised flames. A detailed evaluation of the reliability of ion kinetic and transport parameters adopted was performed. The model provides good quantitative agreement with experimental data for various conditions.
    • Analysis of the Nonlinear Static and Dynamic Behavior of Offshore Structures

      Alfosail, Feras (2015-07)
      Understanding static and dynamic nonlinear behavior of pipes and risers is crucial for the design aspects in offshore engineering fields. In this work, we examine two nonlinear problems in offshore engineering field: vortex Induced vibration of straight horizontal pipes, and boundary layer static solution of inclined risers. In the first study, we analyze the effect of the internal velocity of straight horizontal pipe and obtain the vortex induced vibration forces via coupling the pipe equation of motion with the recently modified Van Der Pol oscillator governing the lift coefficient. Our numerical results are obtained for two different pipe configurations: hinged-hinged, and clamped- clamped. The results show that the internal velocity reduces the vibration and the oscillation amplitudes. Also, it is shown that the clamped-clamped pipe configuration offers a wider range of internal velocities before buckling instability occurs. The results also demonstrate the effect of the end condition on the amplitudes of vibration. In the second study, we develop a boundary layer perturbation static solution to govern and simulate the static behavior of inclined risers. In the boundary layer analysis, we take in consideration the effects of the axial stretch, applied tension, and internal velocity. Our numerical simulation results show good agreement with the exact solutions for special cases. In addition, our developed method overcomes the mathematical and numerical limitations of the previous methods used before.
    • Analysis of Transition from HCCI to CI via PPC with Low Octane Gasoline Fuels Using Optical Diagnostics and Soot Particle Analysis

      An, Yanzhao; Vallinayagam, R; Vedharaj, S; Masurier, Jean-Baptiste; Dawood, Alaaeldin; Izadi Najafabadi, Mohammad; Somers, Bart; Johansson, Bengt (SAE International, 2017-10-10)
      In-cylinder visualization, combustion stratification, and engine-out particulate matter (PM) emissions were investigated in an optical engine fueled with Haltermann straight-run naphtha fuel and corresponding surrogate fuel. The combustion mode was transited from homogeneous charge compression ignition (HCCI) to conventional compression ignition (CI) via partially premixed combustion (PPC). Single injection strategy with the change of start of injection (SOI) from early to late injections was employed. The high-speed color camera was used to capture the in-cylinder combustion images. The combustion stratification was analyzed based on the natural luminosity of the combustion images. The regulated emission of unburned hydrocarbon (UHC), carbon monoxide (CO) and nitrogen oxides (NO) were measured to evaluate the combustion efficiency together with the in-cylinder rate of heat release. Soot mass concentration was measured and linked with the combustion stratification and the integrated red channel intensity of the high-speed images for the soot emissions. The nucleation nanoscale particle number and the particle size distribution were sampled to understand the effect of combustion mode switch.
    • Analytical solutions of the electrostatically actuated curled beam problem

      Younis, Mohammad I. (Springer Science + Business Media, 2014-07-24)
      This works presents analytical expressions of the electrostatically actuated initially deformed cantilever beam problem. The formulation is based on the continuous Euler-Bernoulli beam model combined with a single-mode Galerkin approximation. We derive simple analytical expressions for two commonly observed deformed beams configurations: the curled and tilted configurations. The derived analytical formulas are validated by comparing their results to experimental data and numerical results of a multi-mode reduced order model. The derived expressions do not involve any complicated integrals or complex terms and can be conveniently used by designers for quick, yet accurate, estimations. The formulas are found to yield accurate results for most commonly encountered microbeams of initial tip deflections of few microns. For largely deformed beams, we found that these formulas yield less accurate results due to the limitations of the single-mode approximation. In such cases, multi-mode reduced order models are shown to yield accurate results. © 2014 Springer-Verlag Berlin Heidelberg.
    • Apparatus and method for lagrangian precipitation sensing

      Claudel, Christian G.; Abdelkader, Mohamed; Shaqura, Mohammad; Hoteit, Ibrahim (2015-06-18)
      A weather measurement system can include an aerial vehicle. An array of accelerometers can be disposed across the surface of the aerial vehicle. The accelerometers can be configured to measure raindrops impacting the surface. The measurements can include number per unit time and intensity. The system can include a processing unit configured to receive measurement data from the array and to process the measurement data into preprocessed data.
    • Appearance of deterministic mixing behavior from ensembles of fluctuating hydrodynamics simulations of the Richtmyer-Meshkov instability

      Narayanan, Kiran; Samtaney, Ravi (American Physical Society (APS), 2018-04-19)
      We obtain numerical solutions of the two-fluid fluctuating compressible Navier-Stokes (FCNS) equations, which consistently account for thermal fluctuations from meso- to macroscales, in order to study the effect of such fluctuations on the mixing behavior in the Richtmyer-Meshkov instability (RMI). The numerical method used was successfully verified in two stages: for the deterministic fluxes by comparison against air-SF6 RMI experiment, and for the stochastic terms by comparison against the direct simulation Monte Carlo results for He-Ar RMI. We present results from fluctuating hydrodynamic RMI simulations for three He-Ar systems having length scales with decreasing order of magnitude that span from macroscopic to mesoscopic, with different levels of thermal fluctuations characterized by a nondimensional Boltzmann number (Bo). For a multidimensional FCNS system on a regular Cartesian grid, when using a discretization of a space-time stochastic flux Z(x,t) of the form Z(x,t)→1/-tN(ih,nΔt) for spatial interval h, time interval Δt, h, and Gaussian noise N should be greater than h0, with h0 corresponding to a cell volume that contains a sufficient number of molecules of the fluid such that the fluctuations are physically meaningful and produce the right equilibrium spectrum. For the mesoscale RMI systems simulated, it was desirable to use a cell size smaller than this limit in order to resolve the viscous shock. This was achieved by using a modified regularization of the noise term via Zx,t→1/-tmaxh3,h03Nih,nΔt, with h0=ξh<h0,ξZ+. Our simulations show that for systems with Bo1 deterministic mixing behavior emerges as the ensemble-averaged behavior of several fluctuating instances, whereas when Bo≈1, a deviation from deterministic behavior is observed. For all cases, the FCNS solution provides bounds on the growth rate of the amplitude of the mixing layer.
    • Application of a robust and efficient Lagrangian particle scheme to soot transport in turbulent flames

      Attili, Antonio; Bisetti, Fabrizio (Elsevier BV, 2013-09)
      A Lagrangian particle scheme is applied to the solution of soot dynamics in turbulent nonpremixed flames. Soot particulate is described using a method of moments and the resulting set of continuum advection-reaction equations is solved using the Lagrangian particle scheme. The key property of the approach is the independence between advection, described by the movement of Lagrangian notional particles along pathlines, and internal aerosol processes, evolving on each notional particle via source terms. Consequently, the method overcomes the issues in Eulerian grid-based schemes for the advection of moments: errors in the advective fluxes pollute the moments compromising their realizability and the stiffness of source terms weakens the stability of the method. The proposed scheme exhibits superior properties with respect to conventional Eulerian schemes in terms of stability, accuracy, and grid convergence. Taking into account the quality of the solution, the Lagrangian approach can be computationally more economical than commonly used Eulerian schemes as it allows the resolution requirements dictated by the different physical phenomena to be independently optimized. Finally, the scheme posseses excellent scalability on massively parallel computers. © 2013 Elsevier Ltd.