Blade-resolved numerical simulations of wind energy applications using full blade and tower models are presented. The computational methodology combines solution technologies in a multi-mesh, multi-solver paradigm through a dynamic overset framework. The coupling of a finite volume solver and a high-order, hp-adaptive finite element solver is utilized. Additional technologies including in-situ visualization and atmospheric microscale modeling are incorporated into the analysis environment. Validation of the computational framework is performed on the National Renewable Energy Laboratory (NREL) 5MW baseline wind turbine, the unsteady aerodynamics experimental NREL Phase VI turbine, and the Siemens SWT-2.3-93 wind turbine. The power and thrust results of all single turbine simulations agree well with low-fidelity model simulation results and field experiments when available. Scalability of the computational framework is demonstrated using 6, 12, 24, 48, and 96 wind turbine setups including the 48 turbine wind plant known as Lillgrund. The largest case consisting of 96 wind turbines and a total of 385 overset grids are run on 44,928 cores at a weak scaling efficiency of 86%. Demonstration of the coupling of atmospheric microscale and Computational Fluid Dynamics (CFD) solvers is presented using the National Center for Atmospheric Research (NCAR) Weather Research and Forecasting Model (WRF) solver and the NREL Simulator fOr Wind Farm Applications (SOWFA) solver.