An integrated approach, combining the continuum theory of sintering with a kinetic Monte-Carlo (KMC) model-based mesostructure evolution simulation is reviewed. The effective sintering stress and the normalized bulk viscosity are derived from mesoscale simulations. A KMC model is presented to simulate microstructural evolution during sintering of complex microstructures taking into consideration grain growth, pore migration, and densification. The results of these simulations are used to generate sintering stress and normalized bulk viscosity for use in continuum level simulation of sintering. The advantage of these simulations is that they can be employed to generate more accurate constitutive parameters based on most general assumptions regarding mesostructure geometry and transport mechanisms of sintering. These constitutive parameters are used as input data for the continuum simulation of the sintering of powder bilayers. Two types of bilayered structures are considered: layers of the same particle material but with different initial porosity, and layers of two different materials. The simulation results are verified by comparing them with shrinkage and warping during the sintering of bilayer ZnO powder compacts.