Recent experimental studies of the multiplicity dependence of heavy quark (HQ) production in proton-proton collisions at 7 TeV showed a strong non-linear increase of the HQ multiplicity as a function of the charged particle multiplicity. We try to understand this behavior using the EPOS3 approach. Two issues play an important role: multiple scattering, in particular its impact on multiplicity fluctuations, and the collective hydrodynamic expansion. The results are very robust with respect to many details of the modeling, which means that these data contain valuable information about very basic features of the reaction mechanism in proton-proton collisions.

We describe the details of 3+1 dimensional relativistic hydrodynamic code for the simulations of quark-gluon/hadron matter expansion in ultra-relativistic heavy ion collisions. The code solves the equations of relativistic viscous hydrodynamics in the Israel-Stewart framework. With the help of ideal-viscous splitting, we keep the ability to solve the equations of ideal hydrodynamics in the limit of zero viscosities using a Godunov-type algorithm. Milne coordinates are used to treat the predominant expansion in longitudinal (beam) direction effectively. The results are successfully tested against known analytical relativistic inviscid and viscous solutions, as well as against existing 2+1D relativistic viscous code.

Experimental transverse momentum spectra of identified particles in p-Pb collisions at 5.02 TeV show many similarities to the corresponding Pb-Pb results, the latter ones usually being interpreted in term of hydrodynamic flow. We analyse these data using EPOS3, an event generator based on a 3D+1 viscous hydrodynamical evolution starting from flux tube initial conditions, which are generated in the Gribov-Regge multiple scattering framework. An individual scattering is referred to as Pomeron, identified with a parton ladder, eventually showing up as flux tubes (or strings). Each parton ladder is composed of a pQCD hard process, plus initial and final state linear parton emission. Nonlinear effects are considered by using saturation scales $Q_{s}$, depending on the energy and the number of participants connected to the Pomeron in question. We compute transverse momentum ($p_{t}$) spectra of pions, kaons, protons, lambdas, and $\Xi$ baryons in p-Pb and p-p scattering, compared to experimental data and many other models. In this way we show in a quantitative fashion that p-Pb data (and even p-p ones) show the typical ``flow effect'' of enhanced particle production at intermediate $p_{t}$ values, more and more visible with increasing hadron mass.

We discuss a theoretical scheme that accounts for bulk matter, jets, and the interaction between the two. The aim is a complete description of particle production at all transverse momentum ($p_{t}$) scales. In this picture, the hard initial scatterings result in mainly longitudinal flux tubes, with transversely moving pieces carrying the $p_{t}$ of the partons from hard scatterings. These flux tubes constitute eventually both bulk matter (which thermalizes and flows) and jets. We introduce a criterion based on parton energy loss to decide whether a given string segment contributes to the bulk or leaves the matter to end up as a jet of hadrons. Essentially low $p_{t}$ segments from inside the volume will constitute the bulk, high $p_{t}$ segments (or segments very close to the surface) contribute to the jets. The latter ones appear after the usual flux tube breaking via q-qbar production (Schwinger mechanism). Interesting is the transition region: Intermediate $p_{t}$ segments produced inside the matter close to the surface but having enough energy to escape, are supposed to pick up q-qbar pairs from the thermal matter rather than creating them via the Schwinger mechanism. This represents a communication between jets and the flowing bulk matter (fluid-jet interaction). Also very important is the interaction between jet hadrons and the soft hadrons from the fluid freeze-out. We employ the new picture to investigate Pb-Pb collisions at 2.76 TeV. We discuss the centrality and $p_{t}$ dependence of particle production and long range dihadron correlations at small and large $p_{t}$.

Since 2006, EPOS hadronic interaction model is being used for very high energy cosmic ray analysis. Designed for minimum bias particle physics and used for having a precise description of SPS and RHIC heavy ion collisions, EPOS brought more detailed description of hadronic interactions in air shower development. Thanks to this model it was possible to understand why there were less muons in air shower simulations than observed in real data. With the start of the LHC era, a better description of hard processes and collective effects is needed to deeply understand the incoming data. We will describe the basic physics in EPOS and the new developments and constraints which are taken into account in EPOS 2.

We present a second-order upwind numerical scheme for equations of relativistic hydrodynamics with a source term. A new non-linear Riemann solver is constructed. Solution of a Riemann problem on a cells boundary is based on exact relations in case of zero tangential velocities. In this sense our solver is "exact". In case of non-zero tangential velocities a reasonable approximation is made that allows to avoid solution of very complicated exact equations. The scheme is tested on several one- and three-dimensional solutions demonstrating a good performance for both strong shocks and strong rarefactions.