MASSIVELY PARALLEL FULLY COUPLED IMPLICIT MODELING OF COUPLED THERMAL-HYDROLOGICAL-MECHANICAL PROCESSES FOR ENHANCED GEOTHERMAL SYSTEM RESERVOIRS

Development of enhanced geothermal systems (EGS) will require creation of a reservoir of sufficient volume to enable commercial-scale heat transfer from the reservoir rocks to the working fluid. A key assumption associated with reservoir creation/stimulation is that sufficient rock volumes can be hydraulically fractured via both tensile and shear failure, and more importantly by reactivation of naturally existing fractures (by shearing) to create the reservoir. The advancement of EGS greatly depends on our understanding of the dynamics of the intimately coupled rock-fracture-fluid system and our ability to reliably predict how reservoirs behave under stimulation and production. In order to increase our understanding of how reservoirs behave under these conditions, we have developed a physics-based rock deformation and fracture propagation simulator by coupling a discrete element model (DEM) for fracturing with a continuum multiphase flow and heat transport model. In DEM simulations, solid rock is represented by a network of discrete elements (often referred as particles) connected by various types of mechanical bonds such as springs, elastic beams or bonds that have more complex properties (such as stress-dependent elastic constants). Fracturing is represented explicitly as broken bonds (microcracks), which form and coalesce into macroscopic fractures when external load is applied. DEM models have been applied to a very wide range of fracturing processes from the molecular scale (where thermal fluctuations play an important role) to scales on the order of 1 km or greater. In this approach, the continuum flow and heat transport equations are solved on an underlying fixed finite element grid with evolving porosity and permeability for each grid cell that depends on the local structure of the discrete element network (such as DEM particle density). The fluid pressure gradient exerts forces on individual elements of the DEM network, which therefore deforms and fractures. Such deformation/fracturing in turn changes the permeability, which again changes the evolution of fluid pressure, coupling the two phenomena. The intimate coupling between fracturing and fluid flow makes the meso-scale DEM simulations necessary, as these methods have substantial advantages over conventional continuum mechanical models of elastic rock deformation. The challenges that must be overcome to simulate EGS reservoir stimulation, preliminary results, progress to date and near future research directions and opportunities will be discussed

RELAP-7 Level 2 Milestone Report: Demonstration of a Steady State Single Phase PWR Simulation with RELAP-7

The document contains the simulation results of a steady state model PWR problem with the RELAP-7 code. The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at Idaho National Laboratory (INL). The code is based on INL's modern scientific software development framework - MOOSE (Multi-Physics Object-Oriented Simulation Environment). This report summarizes the initial results of simulating a model steady-state single phase PWR problem using the current version of the RELAP-7 code. The major purpose of this demonstration simulation is to show that RELAP-7 code can be rapidly developed to simulate single-phase reactor problems. RELAP-7 is a new project started on October 1st, 2011. It will become the main reactor systems simulation toolkit for RISMC (Risk Informed Safety Margin Characterization) and the next generation tool in the RELAP reactor safety/systems analysis application series (the replacement for RELAP5). The key to the success of RELAP-7 is the simultaneous advancement of physical models, numerical methods, and software design while maintaining a solid user perspective. Physical models include both PDEs (Partial Differential Equations) and ODEs (Ordinary Differential Equations) and experimental based closure models. RELAP-7 will eventually utilize well posed governing equations for multiphase flow, which can be strictly verified. Closure models used in RELAP5 and newly developed models will be reviewed and selected to reflect the progress made during the past three decades. RELAP-7 uses modern numerical methods, which allow implicit time integration, higher order schemes in both time and space, and strongly coupled multi-physics simulations. RELAP-7 is written with object oriented programming language C++. Its development follows modern software design paradigms. The code is easy to read, develop, maintain, and couple with other codes. Most importantly, the modern software design allows the RELAP-7 code to evolve with time. RELAP-7 is a MOOSE-based application. MOOSE (Multiphysics Object-Oriented Simulation Environment) is a framework for solving computational engineering problems in a well-planned, managed, and coordinated way. By leveraging millions of lines of open source software packages, such as PETSC (a nonlinear solver developed at Argonne National Laboratory) and LibMesh (a Finite Element Analysis package developed at University of Texas), MOOSE significantly reduces the expense and time required to develop new applications. Numerical integration methods and mesh management for parallel computation are provided by MOOSE. Therefore RELAP-7 code developers only need to focus on physics and user experiences. By using the MOOSE development environment, RELAP-7 code is developed by following the same modern software design paradigms used for other MOOSE development efforts. There are currently over 20 different MOOSE based applications ranging from 3-D transient neutron transport, detailed 3-D transient fuel performance analysis, to long-term material aging. Multi-physics and multiple dimensional analyses capabilities can be obtained by coupling RELAP-7 and other MOOSE based applications and by leveraging with capabilities developed by other DOE programs. This allows restricting the focus of RELAP-7 to systems analysis-type simulations and gives priority to retain and significantly extend RELAP5's capabilities

Physics-based multiscale coupling for full core nuclear reactor simulation

Numerical simulation of nuclear reactors is a key technology in the quest for improvements in efficiency, safety, and reliability of both existing and future reactor designs. Historically, simulation of an entire reactor was accomplished by linking together multiple existing codes that each simulated a subset of the relevant multiphysics phenomena. Recent advances in the MOOSE (Multiphysics Object Oriented Simulation Environment) framework have enabled a new approach: multiple domain-specific applications, all built on the same software framework, are efficiently linked to create a cohesive application. This is accomplished with a flexible coupling capability that allows for a variety of different data exchanges to occur simultaneously on high performance parallel computational hardware. Examples based on the KAIST-3A benchmark core, as well as a simplified Westinghouse AP-1000 configuration, demonstrate the power of this new framework for tackling—in a coupled, multiscale manner—crucial reactor phenomena such as CRUD-induced power shift and fuel shuffle.Massachusetts Institute of Technology. Department of Nuclear Science and EngineeringIdaho National Laboratory (Contract DE-AC07-05ID14517

How Many Species Are There on Earth and in the Ocean?

The diversity of life is one of the most striking aspects of our planet; hence knowing how many species inhabit Earth is among the most fundamental questions in science. Yet the answer to this question remains enigmatic, as efforts to sample the world's biodiversity to date have been limited and thus have precluded direct quantification of global species richness, and because indirect estimates rely on assumptions that have proven highly controversial. Here we show that the higher taxonomic classification of species (i.e., the assignment of species to phylum, class, order, family, and genus) follows a consistent and predictable pattern from which the total number of species in a taxonomic group can be estimated. This approach was validated against well-known taxa, and when applied to all domains of life, it predicts ∼8.7 million (±1.3 million SE) eukaryotic species globally, of which ∼2.2 million (±0.18 million SE) are marine. In spite of 250 years of taxonomic classification and over 1.2 million species already catalogued in a central database, our results suggest that some 86% of existing species on Earth and 91% of species in the ocean still await description. Renewed interest in further exploration and taxonomy is required if this significant gap in our knowledge of life on Earth is to be closed

Current and Future Patterns of Global Marine Mammal Biodiversity

Quantifying the spatial distribution of taxa is an important prerequisite for the preservation of biodiversity, and can provide a baseline against which to measure the impacts of climate change. Here we analyse patterns of marine mammal species richness based on predictions of global distributional ranges for 115 species, including all extant pinnipeds and cetaceans. We used an environmental suitability model specifically designed to address the paucity of distributional data for many marine mammal species. We generated richness patterns by overlaying predicted distributions for all species; these were then validated against sightings data from dedicated long-term surveys in the Eastern Tropical Pacific, the Northeast Atlantic and the Southern Ocean. Model outputs correlated well with empirically observed patterns of biodiversity in all three survey regions. Marine mammal richness was predicted to be highest in temperate waters of both hemispheres with distinct hotspots around New Zealand, Japan, Baja California, the Galapagos Islands, the Southeast Pacific, and the Southern Ocean. We then applied our model to explore potential changes in biodiversity under future perturbations of environmental conditions. Forward projections of biodiversity using an intermediate Intergovernmental Panel for Climate Change (IPCC) temperature scenario predicted that projected ocean warming and changes in sea ice cover until 2050 may have moderate effects on the spatial patterns of marine mammal richness. Increases in cetacean richness were predicted above 40° latitude in both hemispheres, while decreases in both pinniped and cetacean richness were expected at lower latitudes. Our results show how species distribution models can be applied to explore broad patterns of marine biodiversity worldwide for taxa for which limited distributional data are available