THE PRINCIPAL OBJECTIVE OF THIS WORK IS TO INVESTIGATE DEVELOP AND DEMONSTRATE ADVANCED SOLUTION STRATEGIES FOR CURRENT AND EMERGING HIGH-ORDER CFD DISCRETIZATIONS WHICH ARE BOTH ROBUST AND EFFICIENT FOR THE SOLUTION OF LARGE-SCALE TRANSIENT PROBLEMS ON MASSIVELY PARALLEL EXASCALE ARCHITECTURES. ALTHOUGH LOCAL NONLINEAR SOLVERS ENJOYED EARLY SUCCESS IN CFD SOLVERS PARTICULARLY WHEN USED AS SMOOTHERS WITHIN A NON-LINEAR MULTIGRID APPROACH THESE HAVE INCREASINGLY BEEN REPLACED BY LINEAR SOLVER TECHNIQUES USED AS DRIVERS FOR NEWTON OR QUASI-NEWTON METHODS DUE TO ROBUSTNESS ISSUES ENCOUNTERED WITH LOCAL NONLINEAR METHODS FOR MORE DEMANDING CFD PROBLEMS INCLUDING HIGHER RESOLUTION CASES MORE CHALLENGING FLOW PHYSICS AND STIFFER NONLINEAR SYSTEMS RESULTING FROM EMERGING HIGH-ORDER DISCRETIZATION SCHEMES. WHILE NEWTON-KRYLOV METHODS ARE WELL SUITED FOR DELIVERING TIGHT CONVERGENCE TOLERANCES PARTICULARLY IN THE FINAL STAGES OF CONVERGENCE THESE METHODS CAN BE EXTREMELY COSTLY OVERALL DUE TO THE REQUIREMENT OF OVERCOMING STRONG NONLINEAR TRANSIENTS IN THE INITIAL STAGES OF CONVERGENCE. ON THE OTHER HAND LOCAL NONLINEAR SOLVERS ARE GENERALLY VERY EFFECTIVE AT OVERCOMING THESE INITIAL STRONG NONLINEARITIES BUT LACK THE RIGOROUS PROPERTIES OF CONTRACTIVITY THAT CAN BE DEVELOPED FOR LINEAR METHODS SUCH AS KRYLOV TECHNIQUES WHICH IN TURN CAN BE USED TO DRIVE A NEWTON METHOD AT LOW RESIDUAL LEVELS WHERE NONLINEAR EFFECTS BECOME VANISHINGLY SMALL. FURTHERMORE THE EMERGENCE OF HETEROGENEOUS ARCHITECTURES WITH LOW MEMORY BANDWIDTH AND MASSIVE PARALLELISM WOULD SEEM TO PUT GLOBAL LINEARIZATION METHODS AT A DISADVANTAGE DUE TO THE LARGE MEMORY REQUIREMENTS OF (APPROXIMATE) MATRIX FACTORIZATION AND POOR SCALABILITY POTENTIAL AT THE HIGH END. HOWEVER THE USE OF LINEAR AND NONLINEAR METHODS NEED NOT BE MUTUALLY EXCLUSIVE AND TECHNIQUES FOR COMBINING THE ADVANTAGES OF BOTH APPROACHES OFFER THE POSSIBILITY OF DEVELOPING SUPERIOR SOLVERS THAT ARE PROVABLY ROBUST NUMERICALLY EFFICIENT AND HIGHLY SCALABLE. THE PRINCIPAL IDEA OF THIS PROPOSAL IS TO DEVISE STRATEGIES THAT MAKE USE OF A SPECTRUM OF BOTH LINEAR AND NONLINEAR SOLUTION TECHNIQUES AT DIFFERENT SCALES PROVIDING AN AVENUE FOR OPTIMIZING ROBUSTNESS EFFICIENCY AND SCALABILITY FOR VARIOUS CFD DISCRETIZATIONS ON VARIOUS TYPES OF HARDWARE. EXAMPLES INCLUDE THE USE OF NONLINEAR RESIDUAL SMOOTHING METHODS USED IN CONJUNCTION WITH NEWTON-KRYLOV METHODS THE APPLICATION OF NONLINEAR (FAS) MULTIGRID METHODS FOR NEWTON METHOD CONTINUATION AND THE TAILORING OF LOCAL LINEAR SOLVERS AND/OR NONLINEAR SMOOTHERS TO BLOCK SIZES THAT ENABLE OPTIMAL USE OF EMERGING HARDWARE CONFIGURATIONS. ADDITIONALLY THESE METHODS ARE DEVISED FOR USE IN MULTILEVEL OR MULTIGRID METHODS IN ORDER TO ENSURE AN OVERALL NUMERICALLY OPTIMAL SOLUTION SCHEME. FOR HIGH-ORDER TEMPORAL DISCRETIZATIONS AN ADDITIONAL DEGREE OF TIME PARALLELISM CAN BE LEVERAGED TO ENHANCE OVERALL SCALABILITY. THE DEVELOPED SOLVERS ARE TO BE DEMONSTRATED WITH DIFFERENT DISCRETIZATIONS INCLUDING CONTINUOUS AND DISCONTINUOUS HIGH-ORDER FINITE-ELEMENT DISCRETIZATIONS AS WELL AS HIGHORDER KINETIC ENERGY OR ENTROPY PRESERVING SCHEMES FOR AERODYNAMIC PROBLEMS OF RELEVANCE TO THE AEROSPACE ENGINEERING COMMUNITY. RATHER THAN ATTEMPTING TO DEVELOP A MONOLITHIC NEW SOLVER APPROACH FOR SPECIFIC CFD PROBLEMS WE SEEK TO EXTEND KNOWLEDGE ABOUT THE PERFORMANCE STRENGTH AND WEAKNESSES OF VARIOUS LINEAR AND NONLINEAR SOLVER TECHNIQUES WHILE ATTEMPTING TO COMBINE THESE COMPONENT TECHNOLOGIES IN OPTIMAL WAYS FOR SPECIFIC CFD DISCRETIZATIONS AND APPLICATIONS. OUR GOAL IS TO REVOLUTIONIZE AEROSPACE SIMULATION CAPABILITIES THROUGH THE DEVELOPMENT AND ADOPTION OF SOLVER TECHNOLOGY THAT IS OPTIMALLY SUITED FOR EXISTING AND NOVEL DISCRETIZATIONS AS WELL AS EMERGING HPC HARDWARE.