Granular flows are frequently observed in nature and appear in many industrial processes as well. In this numerical work the focus is mainly directed at the understanding of how the change of different grain properties, such as shape, friction and stiffness, influences the flow out of a silo. However, the heating dynamics of a granular gas of rods is also analyzed. In all these scenarios, the simulations are paired with experiments to calibrate and validate the results. The numerical analysis indicated that the discharge of soft, low-friction grains from a container exhibits a height-dependent flow rate, which is not usual for granular media. The systematic study mapped the parameter space of particle friction and stiffness, exploring the system’s macroscopic response in detail. Moreover, the examination of the coarse-graining fields helped us to explain when and why the flow rate depends on the column height. The answers include the material response to pressure gradients, but also the way stress is transmitted in the system. The outcomes allowed us to propose simple theoretical arguments, connecting the macroscopic flow rate with the pressure gradient at the orifice. As a result, we have come up with a well-reasoned explanation for the height-dependent discharge flow rate, shown experimentally and numerically by soft low-frictional grains. Our numerical investigation of a 2D silo flow of mixtures of soft and hard grains reproduced the high impact that even 5% of hard frictional grains have on the flow of an ensemble of low-friction, soft particles. Numerical results signaled the importance of the friction between the two types of grains. When the frictional hard grains are added to the soft grains, the flow gets slower, clogging becomes more frequent, and the force measured on the bottom plate decreases. Moreover, we obtained that these effects are enhanced when the interspecies friction is increased. The introduction of a rotational shear through the rotation of the flat silo bottom leads to a surprising effect on the discharge of rods: the flow rate is reduced significantly, by up to 70%. Our simulations and the application of the coarse-graining methods reveal the underlying reasons for this observation. The exit velocity of particles is the main contributing factor to this drop, which is in correlation with the vertical orientation of the grains. Our numerical tool allowed the exploration of the dependence of flow rate 𝑄 on the orifice size 𝐷. In the limit of large orifices, the classical power-law correlation 𝑄 ∼ 𝐷5/2 was reproduced. For small apertures, however, we obtained a power-law relation but with a notably larger exponent. Furthermore, the size of the so-called free fall arch region is also found to decrease due to the rotation. Finally, granular gas made up of rods and heated by the walls has been studied numerically. Our study provides additional insight into the process for instance by describing the system behavior in the non-symmetrically heated direction, which was not accessible experimentally. The velocity distributions of the particles are examined and found to be well-fitted by stretched exponential distributions, in agreement with previous experiments. Moreover, in the experimentally not accessible direction, we obtained asymmetrical distribution tails. Additionally, the collapsing of the velocity distributions lead us to conclude that the energy scales with the square of the characteristic velocity of the wall. Los flujos granulares se observan con frecuencia en procesos naturales e industriales. El enfoque de este trabajo se dirige a la comprensión de cómo el cambio de las diferentes propiedades de los granos (forma, fricción y rigidez) influye en el caudal de los silos. Complementariamente, se analiza la dinámica del calentamiento de un gas granular de partículas alargadas. En todos estos escenarios, las simulaciones se combinan con experimentos para calibrar y validar las herramientas numéricas empleadas. En primer lugar, el análisis numérico de la descarga de granos blandos-lisos de un contenedor indicó el desarrollo de un caudal másico dependiente de la altura, lo cual no es habitual en los medios granulares. Nuestro estudio sistemático examinó el espacio de parámetros de fricción y rigidez de las partículas, explorando detalladamente la respuesta macroscópica de este sistema. Complementariamente, el análisis de los campos medios nos ayudó a explicar cuándo y por qué el caudal depende de la altura de la columna. La explicación general incluye: la respuesta del material a los gradientes de presión en el orificio de salida, pero también la forma en que se transmiten los esfuerzos en el sistema. Los resultados nos permitieron proponer argumentos teóricos simples, conectando el caudal macroscópico con el gradiente de presión en el orificio. Como resultado, hemos encontrado una explicación bien razonada de los valores de caudal dependiente de la altura, lo cual ha sido encontrado, experimental y numéricamente, para granos blandos de baja fricción. En segundo lugar, realizamos una investigación numérica del flujo en un silo bidimensional, usando una de mezclas de granos blandos-lisos y duros-rugosos. Partiendo de un sistema homogéneo de granos blandos-lisos, nuestros resultados numéricos reproducen el alto impacto que produce la inclusión de solo un 5% de los granos duros, en el flujo del conjunto de partículas. Además, resaltaron la importancia de la fricción entre granos de diferente tipo, en le desarrollo de este proceso. Cuando los granos durosrugosos son agregados al sistema, el flujo se vuelve más lento, los atascos resultan más frecuentes y la fuerza medida sobre el fondo disminuye. Además, resulta que estos efectos se potencian cuando aumenta la fricción entre especies. También estudiamos, sistemáticamente, la introducción de esfuerzos cortantes, imponiendo la rotación del fondo de un silo plano. Estas condiciones de contorno producen un efecto sorprendente en la descarga de partículas alargadas: el caudal se reduce significativamente, hasta en un 70%. Nuestras simulaciones y la aplicación de los métodos de promediación espacio-temporales revelan las razones subyacentes de esta observación. La velocidad de salida de las partículas es el principal factor que contribuye a esta caída, lo cual correlaciona con la orientación vertical de los granos. Nuestra herramienta numérica permitió explorar la dependencia del caudal 𝑄 del tamaño del orificio 𝐷. En el límite de los orificios grandes, se reprodujo la clásica correlación de ley de potencia 𝑄 ∼ 𝐷5/2. Sin embargo, para aberturas pequeñas obtuvimos una ley de potencia pero con un exponente notablemente mayor. Además, también se encontró que el tamaño de la región denominada ”arco de caída libre” disminuye debido a la rotación. Finalmente, se ha estudiado numéricamente un gas granular agitado formado por partículas alargadas. Este sistema fue estudiado experimentalmente con anterioridad, y nuestro estudio proporciona información adicional sobre el proceso. Por ejemplo, describiendo el comportamiento del sistema en la dirección de calentamiento asimétrico, al que no se accedió experimentalmente. Se estudiaron las distribuciones de velocidad de las partículas, encontrando que ajustan muy bien con distribuciones exponenciales con colas largas, lo cual está en acuerdo con experimentos previos. Además, en la dirección no accesible experimentalmente, obtuvimos colas asimétricas. El colapso de las distribuciones de velocidad nos lleva a concluir que la energía media del sistema escala con el cuadrado de la velocidad característica de la pared.