Amongst the most important constitutive relations in Mechanics, when characterizing the behavior of complex materials, one can identify the stress tensor T, the heat flux vector q (related to heat conduction) and the radiant heating (related to the radiation term in the energy equation). Of course, the expression 'complex materials' is not new. In fact, at least since the publication of the paper by Rivlin & Ericksen (1955), who discussed fluids of complexity (Truesdell & Noll, 1992), to the recently published books (Deshpande et al., 2010), the term complex fluids refers in general to fluid-like materials whose response, namely the stress tensor, is 'non-linear' in some fashion. This non-linearity can manifest itself in variety of forms such as memory effects, yield stress, creep or relaxation, normal-stress differences, etc. The emphasis in this chapter, while focusing on the constitutive modeling of complex fluids, is on granular materials (such as coal) and non-linear fluids (such as coal-slurries). One of the main areas of interest in energy related processes, such as power plants, atomization, alternative fuels, etc., is the use of slurries, specifically coal-water or coal-oil slurries, as the primary fuel. Some studies indicate that the viscosity of coal-water mixtures depends not only on the volume fraction of solids, and the mean size and the size distribution of the coal, but also on the shear rate, since the slurry behaves as shear-rate dependent fluid. There are also studies which indicate that preheating the fuel results in better performance, and as a result of such heating, the viscosity changes. Constitutive modeling of these non-linear fluids, commonly referred to as non-Newtonian fluids, has received much attention. Most of the naturally occurring and synthetic fluids are non-linear fluids, for example, polymer melts, suspensions, blood, coal-water slurries, drilling fluids, mud, etc. It should be noted that sometimes these fluids show Newtonian (linear) behavior for a given range of parameters or geometries; there are many empirical or semi-empirical constitutive equations suggested for these fluids. There have also been many non-linear constitutive relations which have been derived based on the techniques of continuum mechanics. The non-linearities oftentimes appear due to higher gradient terms or time derivatives. When thermal and or chemical effects are also important, the (coupled) momentum and energy equations can give rise to a variety of interesting problems, such as instability, for example the phenomenon of double-diffusive convection in a fluid layer. In Conclusion, we have studied the flow of a compressible (density gradient type) non-linear fluid down an inclined plane, subject to radiation boundary condition. The heat transfer is also considered where a source term, similar to the Arrhenius type reaction, is included. The non-dimensional forms of the equations are solved numerically and the competing effects of conduction, dissipation, heat generation and radiation are discussed. It is observed that the velocity increases rapidly in the region near the inclined surface and is slower in the region near the free surface. Since R{sub 7} is a measure of the heat generation due to chemical reaction, when the reaction is frozen (R{sub 7}=0.0) the temperature distributions would depend only on R{sub 1}, and R{sub 2}, representing the effects of the pressure force developed in the material due to the distribution, R{sub 3} and R{sub 4} viscous dissipation, R{sub 5} the normal stress coefficient, R{sub 6} the measure of the emissivity of the particles to the thermal conductivity, etc. When the flow is not frozen (RP{sub 7} > 0) the temperature inside the flow domain is much higher than those at the inclined and free surfaces. As a result, heat is transferred away from the flow toward both the inclined surface and the free surface with a rate that increases as R{sub 7} increases. For a given temperature, an increase in {zeta} implies that the activation energy is smaller and thus, the reaction rate is increased leading to an increase in the heat of the reaction. As a result the flow is chemically heated and its temperature increase. The results shown here indicate that for all values of {zeta} used the chemical effects are significant and the temperature is always higher than both the surface temperature and the free surface temperature. The heat transfer is always from the flow toward both the inclined surface and the free stream. It is also noticed that for all values of m chosen in this study, the temperature is higher than the surface and the free stream temperature. The heat transfer at the inclined surface and at the free stream increase slowly for negative values of m to about m=0.5, but it begins to significantly increase for m greater than 0.5.