Search NASA⌕ Search

SEARCH · Search NASA

Results for “Site closure requirements”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

98 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Middlesex South FUSRAP Site: Collaboration Towards Beneficial Reuse - 20338

The US Department of Energy (DOE) Office of Legacy Management (LM) manages DoE's post-closure responsibilities and ensures the future protection of human health and the environment with respect to sites that have no continuing DOE mission after undergoing remediation. LM's beneficial reuse program promotes the LM strategic goal to sustainably manage and optimize use of public lands. The beneficial reuse program aims to repurpose former contaminated sites to restore the environment, protect the public health, revitalize communities, and spur economic growth. Benefits from the reuse of a site may increase the local tax base, facilitate job growth, utilize existing infrastructure, and enhance or protect natural resources. Beneficial reuse also promotes protectiveness by ensuring activities are compatible with long-term maintenance and protection of public health and the environment, as well as by retaining good stewardship of natural resources. LM actively participates in promoting the economic development vision of the surrounding communities by collaborating with local communities to promote regional or municipal initiatives. The Borough of Middlesex, New Jersey, has prepared the Lincoln Boulevard Redevelopment Plan, which will reinvent a once-viable downtown and enhance the quality of life for the community. The Middlesex South, New Jersey, Formerly Utilized Sites Remedial Action Program (FUSRAP) site lies within the boundaries of this redevelopment area. In support of this development plan, the Borough is interested in acquiring the site for its Department of Public Works and to increase street access in the future commercial district to be built near adjacent properties. The opportunity to put the FUSRAP site into productive use aligns with LM's mission and strategic goals. In support of beneficial reuse activities, DOE can dispose of excess real property using several mechanisms. The two disposal options applicable to the Middlesex south site include (1) utilizing the US General Services Administration (GSA) or (2) using the Title 10 Code of Federal Regulations Section 770, 'Transfer of Real Property at Defense Nuclear Facilities for Economic Development' (known as a '770 transfer') process. Using the 770 transfer requires economic development as the primary driver, and a specific redevelopment proposal must be submitted. If there is no interest in specific economic development, the GSA option becomes the default mechanism for disposal of the property. In the case of the Middlesex south site, the 770 transfer process is being utilized and the Borough has submitted its economic development proposal to LM. Although LM is encouraged by the opportunity to put the site into productive use, the Middlesex south site is currently on the US Environmental Protection Agency (EPA) National Priorities List and remediation is being performed by the US Army Corp of Engineers (USACE) under the Comprehensive Environmental Response, Compensation, and Liability Act (CERCLA). Consequently, any property transfer must meet the requirements under CERCLA Section 120(h), 'Property Transferred by Federal Agencies.' As the groundwater remedy will not be in place for several years, approval would be required from the US EPA administrator and the governor of the state of New Jersey to perform an early transfer of the site (before all response actions have been completed). The time frame for these approvals to occur under the normal process conflicts with the proposed redevelopment schedule; therefore, various options for expedited sale, transfer, or lease of portions of the site to the Borough were evaluated. The subsequent decision on a mutually beneficial path forward was the result of a collaborative effort between LM, USACE, the Borough and its redevelopment team, EPA, and the New Jersey Department of Environmental Protection. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Vadose Zone Model for SALDS for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort: carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP1) simulator, officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP) is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the State-Approved Land Disposal Site (SALDS) model. The scope of this ECF is to document the development and results of the SALDS vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Vadose Zone Model for B-3A/B Pond for Composite Analysis

The objectives of the vadose modeling for the updated Hanford Site composite analysis (CA) are to simulate the flow and transport of water and radionuclide releases from the surface to the water table and to provide radionuclide transfer rates for the plateau to river (P2R) model, version 8.3 (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). Water additions include natural recharge and water discharged to the ground as a result of industrial processes associated with Hanford Site operations. Contaminant sources include radionuclides in water discharged to the ground during operations and radionuclides disposed “dry” in solid waste burial grounds or other means. The following 16 radionuclides were selected for this modeling effort; carbon-14 (C-14), chlorine-36 (Cl-36), tritium (H-3), iodine-129 (I-129), neptunium-237 (Np-237), rhenium-187 (Re-187), strontium-90 (Sr-90), technetium-99 (Tc-99), uranium-232 (U-232), uranium-233 (U-233), uranium-234 (U-234), uranium-235 (U-235), uranium-236 (U-236), uranium-238 (U-238), radium-226 (Ra-226), and thorium-230 (Th-230). The simulation time starts in 1943 and ends at 12070, which is 10,000 years after assumed Hanford Site closure in 2070. The parallel version of the Subsurface Transport Over Multiple Phases (STOMP) simulator officially named the exascale Subsurface Transport Over Multiple Phases (eSTOMP), is used to simulate flow and transport for the vadose models. The documentation for the STOMP code is comprehensive. The theoretical and numerical approaches applied in the STOMP code are documented in a published theory guide (PNNL-12030, STOMP Subsurface Transport Over Multiple Phases Version 2.0 Theory Guide). The code has undergone a rigorous verification procedure against analytical solutions, laboratory-scale experiments, and field-scale demonstrations. The application guide (PNNL-11216, STOMP Subsurface Transport Over Multiple Phases Application Guide) provides instructive examples in the application of the code to classical groundwater problems. The user’s guide (PNNL-15782, STOMP: Subsurface Transport Over Multiple Phases Version 4.0: User’s Guide) describes the general use, input file formatting, compilation, and execution of the code. The primary output of the vadose zone modeling is radionuclide transfer rates to the groundwater for input into the P2R model. The rates will be summed over the 100 by 100 m P2R grid cells that fall within the vadose zone model source domain. The Hanford Site Central Plateau was subdivided into 26 individual vadose zone models, with 13 in the 200 East Area and 13 in the 200 West Area. Waste sites that have a completed performance assessment (PA) or past-leak analysis were not included as sources of radionuclides. Instead the vadose zone to groundwater transfer rates of the Environmental Restoration Disposal Facility, Integrated Disposal Facility, US Ecology, and Waste Management Area C (WMA C) PAs and the past-leak analysis for WMA C were used as direct input to the P2R model. Each of the vadose zone models is documented in separate environmental calculation files (ECFs). This ECF describes the B-3A/B Ponds model. The scope of this ECF is to document the development and results of the B-3 A/B Ponds vadose zone model. CP-63515, Model Package Report: Central Plateau Vadose Zone Models, describes the approach, assumptions, process of determining the number of models required and domain of each model, input data, and processing common to all the models.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Syracuse University Industrial Assessment Center (Final Report)

Syracuse University Industrial Assessment Center (SU-IAC) project is focused on the two primary objectives: (i) educate undergraduate and graduate students in engineering and associated disciplines on the concepts of energy sustainability and to provide hands-on training by performing energy assessment at industrial facilities to become future energy efficiency experts, and (ii) offer energy assessments to small and mid-sized manufacturing companies in SU-IAC 2 the greater New York State region, perform analytical engineering work in support of assessment recommendations, and to report the results of those analyses to the client companies, to DOE, and to Field office personnel Through this project period, SU-IAC successfully trained 77 engineering students who were undergraduate or graduate students at Syracuse University College of Engineering, or at SU- IAC’s satellite center at Clarkson University. Each of these students received comprehensive ‘in-class’ technical training, safety training, and hands-on field training. The students performed various tasks related to energy assessment at the manufacturing facilities, including billing analysis, energy use reduction recommendations, and report preparation. Through this project period, a total of 34 students have completed all requirements of the SU-IAC program and have received the US DoE Certificate on Energy Efficiency, with the remaining students continuing at various stages of field work or training. During this project period, due to Covid-related measures in New York State, Syracuse University campus and most of the manufacturing facilities were completely closed or were under severe restrictions for access, for an extended period from March 2020 through April 2021. The closures affected our project execution, with the project Sponsor offering structured relief to carry-out a modified project during the affected budget periods. During this project period, SU-IAC student teams, led by its Director, completed a total of 77 energy assessments at small- and medium-scale manufacturing facilities located in New York State. For these facilities, SU-IAC prepared and presented details for 569 assessment recommendations (ARs). The facilities implemented a total of 258 of these ARs, realizing annual average energy savings in the range of $\$ $3,300 to $\$ $55,000 per facility. Cumulatively, these facilities have reduced their “first year after implementation” energy use by a total of 4.37 million kwh, and fuel use by a total of 27.7 billion BTU. SU-IAC project has successfully met the two primary objectives (and the associated sub objectives) and helped the manufacturing facilities in New York State achieve quantitative energy use reductions as a result of this project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A Perspective on the Successes of the NNSS Underground Test Area (UGTA) Activity - 20221

Between 1951 and 1992, 828 underground nuclear tests were performed at the Nevada National Security Site (NNSS). Underground testing was done in five major testing areas, which included (1) Frenchman Flat, (2) Rainier Mesa/Shoshone Mountain, (3) Yucca Flat/Climax Mine (4) Central Pahute Mesa, and (5) Western Pahute Mesa. About one third of the underground tests were detonated near, at or below the water table, and thus radioactive contamination was introduced to the groundwater system. The U.S. Department of Energy's Underground Test Area (UGTA) Activity was established in the late 1990's to address the fate and movement of residual radioactivity in groundwater, and characterize the risk that it may pose to NNSS workers and the offsite public. It has accomplished this goal through a process of iterative groundwater sampling and numerical groundwater flow and transport modeling, and long-term monitoring. The UGTA Activity draws on the expertise of scientists in the fields of geology, hydrology, radiochemistry, and risk assessment from the U.S. DOE staff, the lead contractor (currently Navarro) and many other organizations, in cooperation with the governing regulatory body, the State of Nevada's Division of Environmental Protection (NDEP). Los Alamos National Laboratory (LANL) is one of the participating research organizations involved with the UGTA Activity. Under the direction of DOE and the lead contractor, LANL's role in UGTA has evolved over time from an initial focus on geologic characterization, to conducting field and laboratory experiments, and finally to its current role of providing modeling and geochemistry expertise to characterize the rates and directions of groundwater and radionuclide movement. In its current role, LANL has either developed or provided the numerical tools for developing flow and transport models in each of the four major testing areas. In addition, LANL has provided an independent assessment of future groundwater flow and transport behavior through the analysis of naturally-occurring geochemical and isotopic tracers in groundwater. The five major testing areas are now in different stages of investigation or closure: (1) Frenchman Flat is in its fifth year of post closure modeling; (2) The Rainier Mesa/Shoshone Mountain Closure Report has been submitted to NDEP and is awaiting approval; (3) The Yucca Flat/Climax Mine Closure Report is under development; and (4) Central and (5) Western Pahute Mesa is still undergoing investigation and flow and transport model development. The success of the UGTA Activity in reaching or approaching closure in three of the five major testing areas is primarily due to the DoE's focus on regulatory closure requirements, but also is in part attributable to a rigorous internal review process and the involvement of NDEP and Nye County water professionals as participants and observers in the reviews. County and state involvement in the internal review process has promoted trust that the U.S. DOE is deeply committed to ensuring the present and future safely of NNSS workers and the surrounding communities. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Nuclear Criticality Safety Assessment of Criticality Control Containers without Moderation Control at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) provides for safe, permanent disposal of government-owned transuranic (TRU) and TRU mixed wastes. Receipt and disposal of waste at the WIPP site began in March 1999. The Sandia report, Consideration of Nuclear Criticality When Disposing of Transuranic Waste at the Waste Isolation Pilot Plant, addressed potential nuclear criticality safety issues based on the projected inventory characteristics known at the time [1]. As designs for inventory, waste forms, and disposal packages have changed, new analyses have been performed, and updates have been made to address any potential effects to the WIPP safety basis. New analyses performed include Saylor 2017 [2] and Brickner 2019 [3], which address certain waste containers with specified loadings under post-closure conditions. Both examined several hypothetical scenarios and included analyses to bound (from a criticality potential standpoint) credible configurations that could occur at WIPP during the repository regulatory post-closure disposal time period for feature, event, and process (FEP) considerations—10,000 years. During this post-closure period at WIPP, the screening of FEPs is governed by the risk-based standards and implementing regulations of the US Environmental Protection Agency (EPA) (i.e., 40 CFR 191 and 40 CFR 194, respectively) [4,5]. An FEP screening can be based on either a low-consequence or low-probability rationale. A low-probability rationale includes either (a) a qualitative rationale that the FEP is not credible or (b) a quantitative demonstration that the probability is less than 10-4 in 104 years. In this evaluation, a qualitative lowprobability rationale of not credible is used by demonstrating that bounding configurations of the waste are not critical. The demonstration of subcriticality is through quantitative calculations, but a probability of criticality is not evaluated. Rather, the rationale for this evaluation is that bounding configurations with an effective neutron multiplication factor (keff) well below the upper subcriticality limit (USL) make criticality incredible. Reference [2] documented a nuclear criticality assessment of the WIPP repository for disposal of dilute surplus plutonium materials using the Dilute and Dispose Approach and packaging in criticality control overpacks (CCOs). The CCO is the waste disposal container recently designed to allow for up to 380 fissile gram equivalent (FGE) 239 Pu per drum, which is a higher fissile loading than typical waste containers. The CCO consists of a criticality control container (CCC) positioned by upper and lower plywood spacers within a standard 55 gal drum. The CCC is used to establish a geometry control for fissile materials during transportation and WIPP emplacement operations. The current WIPP waste acceptance criteria for CCO payloads limit beryllium to less than or equal to 1% by weight of the waste contents and require the waste form to be non-machine compacted. Reference [2] considered two scenario progressions—room closure from salt creep, hereafter referred to as the reconfigured dry scenario, and flooding with brine, hereafter referred to as the reconfigured wet scenario. The subsequent drying out of the reconfigured wet scenarios was also considered. For all scenarios, subcriticality was maintained when 50 g of B 4 C (acting as a neutron absorber) per CCC was intermixed within the plutonium disposition waste form. The analysis used a waste form description that limits the amount of moderation that could be present within the waste form (i.e., it limits the amount of water and polyethylene that could be present based on planned processing conditions). This analysis to evaluate increased limits on the amount of moderation that could be present was performed as a companion to Reference [2] to address concerns associated with verifying moisture and/or plastic contents of waste materials following packaging of dilute surplus plutonium in the CCO. To that end, this analysis used the models and methods from Reference [2] to evaluate a more generic base waste form consisting of water and polyethylene that is more similar (and nearly identical) to the generic waste forms utilized in other models/analyses supporting the TRU Package Transporter Model II (TRUPACTII) safety analysis [6] (all are without moderation controls). The waste form in this analysis uses a base mixture of 75% water and 25% polyethylene, the total amount of which is varied to determine the optimum moderation to fissile material (H/Pu) ratio. The fissile loading is maintained at up to 380 FGE 239 Pu (modeled as PuO 2 ) per CCO with an additional 545 g of beryllium (to bound the 1% by weight contents restriction) and 50 g of B 4 C intermixed per CCO. The beryllium content (1% by weight) is based on the total allowed waste weight (this does not include packaging and container weights). Figures ES-1 and ES-2 display summary results, showing that with this model including 50 g of B 4 C per CCO, the system keff remains under 0.85 for all moderator amounts and provides a significant margin against post-closure criticality under postulated bounding conditions for compaction. Figure ES-1 compares an infinite model with a room model at the initial emplacement spacing and under full radial compaction. Full radial compaction places each CCC in direct contact and does not credit any anticipated spacing associated with current post-closure geomechanical modeling of the repository [7]. The effects of variations in the H/Pu ratio were evaluated by varying the amount of the water/polyethylene component of the waste model, with fissile loading maintained at 380 239 Pu FGE. Similarly, Figure ES-2 illustrates how various amounts of B 4 C per CCO influence k eff at different radial compactions, all at the H/Pu ratio of 200 (in the room array model). Therefore, while the results from Saylor 2017 [2] modeled more realistic process limits associated with packaging of dilute surplus plutonium, this analysis demonstrates that limits on moderation (plastic and water content) are not necessary to ensure subcriticality in the WIPP repository, provided the requisite B 4 C absorber is present.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

SRS Strategy for Tank 3 Salt Dissolution - Performance and Lessons Learned - 20441

The Savannah River Site (SRS) Tank Farms have 51 underground waste tanks used to store and process liquid nuclear waste materials. There are 4 different tank types, ranging in capacity from 2,840,000 to 4,920,000 L (750,000 to 1,300,000 Gal). Twenty-four of the tanks are older style and do not meet full secondary containment standards. The older style tanks are the initial focus of waste removal efforts for tank closure at SRS. Eight of these twenty-four tanks have completed waste removal and are filled with grout. Prior to salt dissolution, Tank 3 was a dry salt tank that contained 5.06 m (199.3 in) of salt and sludge waste. Additional salt waste was present on cooling coils above the salt layer up to approximately 5.59 m (220 in). Three mixing eductors were installed in Tank 3 to aid in dissolving salt waste in three tank riser access ports. A transfer pump was installed, and the transfer pump suction was located 25.4 cm (10 in) from the tank bottom. Well water was added to the tank through a downcomer until the dry bulk salt was covered with liquid. During the initial fill of Tank 3, approximately 242,000 L (64,000 Gal) of rain water were added from periodic F-Tank Farm (FTF) Catch Tank additions and approximately 17,000 L (4,500 Gal) of well water were added. Following liquid additions to cover the dry bulk salt, well water was added through the three mixing eductors in batches during each stage. Additionally, during the salt dissolution campaigns, the FTF Catch Tank was utilized to add rain water through a downcomer in the center tank riser access port as needed for volume relief in the FTF Catch Tank. Following liquid additions, the water was recirculated (internal to the tank) using the transfer pump, and a sample was pulled to confirm the target specific gravity (SpG) of the dissolved salt solution was achieved. The dissolved salt solution was then transferred to the receipt tank (Tank 7), and the mixing eductors were lowered as close to the new bulk salt layer as possible, to support subsequent dissolution campaigns. While the mixing jets were able to dissolve salt successfully in Tank 3, they did not do so in a completely uniform manner. Throughout dissolution, mounds were discovered under Riser 1, Riser 2, and the Center Riser. The mixing jets dissolved the salt around the edges of the tank well but were not as effective toward the center of the tank. FTF Catch Tank additions through the center tank riser access port downcomer were required to impact the mound under the center tank riser. FTF Catch Tank additions were sporadic as they were dependent on rainfall in the area. Additionally, all mixing jets were rarely able to be lowered to the same elevation, indicating some mounding in the bulk salt layer. Indexing of the mixing jets was also utilized to impact the salt mounds. Overall, Tank 3 salt dissolution was successful as approximately 1,476,000 L (390,000 Gal) of dissolved salt solution was transferred to Tank 7 over six stages. Slightly less than 852,000 L (225,000 Gal) of water were added to Tank 3. The original material balance prediction estimated that the bulk salt level in Tank 3 would be 3.88 m (152.8 in) after six salt dissolution stages. After six stages the actual bulk salt level was 3.38 m (133 in). This paper will discuss the salt removal strategy, each salt dissolution stage, and lessons learned for future salt dissolution. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗