The main objective of this work is to develop computational models to analyze electromagnetic logging-while-drilling tool response in arbitrary geophysical formations. This modeling requires the determination of electromagnetic fields in three- dimensional (3-D) complex regions and consequently, the solution of large scale non-hermitian systems. The numerical modeling is done by using Finite Volume Methods (FVM) in the frequency domain. Both isotropic and anisotropic models are developed. Maxwell's equations are solved by using both the field formulation and the coupled vector-scalar potentials formulation. The proposed FVM technique utilizes an edge-based staggered-grid scheme in cylindrical coordinates to avoid staircasing errors on the tool geometry. Four numerical techniques are incorporated in the models in order to increase the computational efficiency and the accuracy of the method. The field formulation and the coupled vector-scalar potentials formulation are compared in terms of their accuracy, convergence rate, and CPU time for three-dimensional environments. The models were validated and tested in 3-D complex environments, such as:(i) horizontal and directional boreholes; (ii) multilayered geophysical formations including mud-filtrate invasions; (iii) anisotropic formations and (iv)eccentric boreholes. The methodology is extended to analyze LWD tools that are constructed with the transmitters and/or receivers tilted with respect to the axis of the drill collar. Such tools can provide improved anisotropy measurements and azimuthal sensitivity to benefit geosteering.