Steady airflow pattern during a full breathing cycle in human upper and central respiratory tract was simulated by solving the Navier-Stokes and continuity equations. For ellipsoidal fiber trajectory analysis under cyclic breathing condition, several user defined functions (UDFs) were developed and coupled to the ANSYS-Fluent discrete phase model (DPM). The developed model accounted for solving the coupled translational and rotational equations of motion of ellipsoidal fibers. The airway passage model was extended from the vestibule to the fifth generation of the bronchial bifurcations obtained mostly from computed tomography (CT) scan. A constant flow rate of 15 L/min was used to simulate the normal breathing condition. The velocity and pressure fields for different regions of the respiratory track were evaluated and used for Lagrangian particle trajectory analysis. Total and regional depositions of each region for a range of ellipsoidal particle diameter and aspect ratios were evaluated and the results compared with the experimental data.