The approval of the current dissertation by the Department of Mechanical Engineering of the University of Thessaly does not imply acceptance of the author's opinions (Law 5343/32 number 202 paragraph 2). Also, the views and opinions expressed herein do not necessarily reflect those of the European Commission. There are several other people who I feel that have contributed in the writing of this dissertation.The most important contributions, along with several helpful conversations, were made by Prof. F.Sharipov, who provided the basic concept of the channel end effect study, and Dr. Chr. Day along with the KIT personnel, who briefly introduced me in the world of experimental rarefied gas dynamics. It was also a priviledge to discuss and work for a brief period of time with Prof. A. Grecos, who tirelessly discussed with me on various subjects of statistical mechanics. Financial support for this Ph.D. has been provided by the Euratom Association -Hellenic Republic, to which I am deeply indepted. Also, partial support for conference mobility expenses has been received from the GASMEMS Marie Curie programme. In this work, non-equilibrium transport phenomena have been examined via the kinetic theory of gases. The behaviour of the gas in low pressure conditions or in low-dimensionality systems can not be captured by the usual Navier-Stokes formulation, in conjunction with the constitutive laws of Newton-Fourier-Fick. The limited number of intermolecular collisions results in departure from equilibrium and the particulate nature of the gas must be taken into account, greatly increasing the computational effort.The problem is described by the integro-differential Boltzmann equation, which is used to determine the distribution of particles in physical and molecular velocity space, as well as in time.There are difficulties associated with the seven-dimensional nature of the distribution function, as well as with the complexity of the collision term, which is usually substituted by an appropriate model.The most widely used and successful numerical methodologies, the Discrete Velocity Method (DVM) and the Direct Simulation Monte Carlo (DSMC), are applied here in the whole range of the Knudsen number. It is important to employ approaches based on kinetic theory, since these are the only ones which are valid for any rarefaction level.In this work, several problems including non-equilibrium phenomena are considered. The interaction of gases with solid surfaces is considered for several problems according to the CercignaniLampis boundary conditions. The non-linear form of this scattering kernel, which had not been pre- viously used in the literature, is applied on the problem of heat transfer between parallel plates and coaxial cylinders. Its linearized form has also been employed for both flow and heat transfer problems and a comparison with relevant experiments has lead to surface characterization with respect to argon and helium flows.Non-linear heat conduction phenomena are also studied for the parallel plates and coaxial c...