MASSIVE STARS PLAY AN IMPORTANT ROLE IN MANY ASTROPHYSICAL ENVIRONMENTS. ONE IMPORTANT PROPERTY OF MASSIVE STARS IS THAT THEY USUALLY SHOW SIGNIFICANT MASS LOSS WHICH CANNOT BE EXPLAINED BY THE NORMALLY ASSUMED RADIATION PRESSURE ON THE ULTRAVIOLET LINES. THEORETICALLY STRUCTURES OF MASSIVE STARS AND ESTIMATES OF MASS LOSS BASED ON 1D STELLAR EVOLUTION CALCULATIONS ARE VERY UNCERTAIN DUE TO THE LARGE RADIATION PRESSURE THE OPACITY AND VARIOUS RADIATION (MAGNETO-)HYDRODYNAMIC INSTABILITIES EXPECTED IN THIS REGIME. WE WILL CARRY OUT FIRST EVER FIRST PRINCIPLE THREE DIMENSIONAL GLOBAL RADIATION MHD SIMULATIONS OF MASSIVE STAR ENVELOPE BY SOLVING THE TIME DEPENDENT RADIATION TRANSFER EQUATION ALONG DISCRETE RAYS DIRECTLY FOR A WIDE RANGE OF MASSES AT DIFFERENT EVOLUTION STAGES. THESE SIMULATIONS WILL RESOLVE THE CONVECTION DEVELOPED AROUND THE IRON OPACITY BUMP REGION DUE TO THE SUPER-EDDINGTON RADIATION FLUX AND STILL CAPTURE THE GLOBAL STRUCTURE AND DIFFERENT ENERGY TRANSPORT MECHANISMS IN THE ENVELOPE. THE RESULTS OF OUR MULTI-DIMENSIONAL SIMULATIONS WILL DRAMATICALLY IMPROVE THE UNDERSTANDING OF THE SURFACE LAYERS OF MASSIVE STARS INCLUDING CONVECTIVE ENERGY TRANSPORT IN THE RADIATION PRESSURE DOMINATED REGIME AND MASS LOSS VIA RADIATION-PRESSURE DRIVEN WINDS. WE WILL ALSO USE A 3D MONTE CARLO RADIATIVE TRANSFER CODE TO CALCULATE THE LINE FORCE IN THESE STARS BASED ON THE REALISTIC TURBULENT ENVELOPE STRUCTURES WE GET FROM THE SIMULATIONS. THESE CALCULATIONS CAN THEN BE USED TO IMPROVE THE CLASSICAL CAK THEORY ON LINE DRIVEN WINDS IN THIS REGIME. OUR RESULTS WILL BE INCORPORATED INTO ONE-DIMENSIONAL STELLAR EVOLUTION MODELS TO CREATE MORE REALISTIC MASSIVE STAR MODELS AND SUPERNOVAE PROGENITORS WHICH WILL SIGNIFICANTLY IMPROVE OUR UNDERSTANDING OF STRUCTURES AND EVOLUTIONS OF MASSIVE STARS. OVER THE PAST FEW YEARS OUR TEAM HAS SHOWN OUR CAPABILITY TO ACHIEVE THESE GOALS AND WE HAVE DEVELOPED ENOUGH TOOLS TO PERFORM THESE SIMULATIONS. WE DEVELOPED THE 1D MESA STELLAR EVOLUTION CODE TO CREATE STELLAR STRUCTURES AT DIFFERENT EVOLUTION STAGES ACCORDING TO THE CLASSICAL 1D ASSUMPTIONS WHICH PROVIDES US WITH REALISTIC PARAMETERS FOR 3D SIMULATIONS. WE HAVE ALSO DEVELOPED ACCURATE RADIATIVE TRANSFER ALGORITHMS IN THE NEWLY DEVELOPED 3D MHD CODE ATHENA++ WHICH HAS ALREADY BEEN SUCCESSFULLY USED TO STUDY THE STRUCTURES OF MASSIVE STAR ENVELOPES BASED ON A SERIES OF LOCAL AND GLOBAL SIMULATIONS. FINALLY WE HAVE ALSO DEVELOPED A MATURE MONTE CARLO CODE TO POST PROCESS THE SIMULATION DATA TO CALCULATE THE RADIATION FORCE DUE TO THE LINES IN THE TURBULENT ENVELOPE DIRECTLY. THESE FIRST-PRINCIPLE CALCULATIONS WILL ADDRESS FUNDAMENTAL QUESTIONS ON STRUCTURES AND EVOLUTIONS OF MASSIVE STARS (SUCH AS CONVECTION AND MASS LOSS IN THE RADIATION PRESSURE DOMINATED REGIME) THAT HAS BEEN IN THIS FIELD FOR DECADES. THE RESULTS OF OUR CALCULATIONS WILL DELIVER IMPROVED PRESCRIPTIONS FOR THE WHOLE COMMUNITY TO USE IN 1D STELLAR EVOLUTION MODELS WHICH WILL REMOVE SEVERAL KEY UNCERTAINTIES IN THE 1D CALCULATIONS.