Software Engineering for Scientific Computer Simulations

NASA Astrophysics Data System (ADS)

Post, Douglass E.; Henderson, Dale B.; Kendall, Richard P.; Whitney, Earl M.

2004-11-01

Computer simulation is becoming a very powerful tool for analyzing and predicting the performance of fusion experiments. Simulation efforts are evolving from including only a few effects to many effects, from small teams with a few people to large teams, and from workstations and small processor count parallel computers to massively parallel platforms. Successfully making this transition requires attention to software engineering issues. We report on the conclusions drawn from a number of case studies of large scale scientific computing projects within DOE, academia and the DoD. The major lessons learned include attention to sound project management including setting reasonable and achievable requirements, building a good code team, enforcing customer focus, carrying out verification and validation and selecting the optimum computational mathematics approaches.

The International Negotiation Seminars Project. Project ICONS.

ERIC Educational Resources Information Center

Wilkenfeld, Jonathan; Kaufman, Joyce; Starkey, Brigid

This report of a study at the University of Maryland describes an international, interactive, and interdisciplinary project for first- and second-year students, which combines a large lecture format with small-group, seminar-type sessions organized around a computer-assisted simulation model, the International Communication and Negotiation…

Projected changes in extreme precipitation over Scotland and Northern England using a high-resolution regional climate model

NASA Astrophysics Data System (ADS)

Chan, Steven C.; Kahana, Ron; Kendon, Elizabeth J.; Fowler, Hayley J.

2018-03-01

The UK Met Office has previously conducted convection-permitting climate simulations over the southern UK (Kendon et al. in Nat Clim Change 4:570-576, 2014). The southern UK simulations have been followed up by a new set of northern UK simulations using the same model configuration. Here we present the mean and extreme precipitation projections from these new simulations. Relative to the southern UK, the northern UK projections show a greater summertime increase of return levels and extreme precipitation intensity in both 1.5 km convection-permitting and 12 km convection-parameterised simulations, but this increase is against a backdrop of large decreases in summertime mean precipitation and precipitation frequency. Similar to the southern UK, projected change is model resolution dependent and the convection-permitting simulation projects a larger intensification. For winter, return level increases are somewhat lower than for the southern UK. Analysis of model biases highlight challenges in simulating the diurnal cycle over high terrain, sensitivity to domain size and driving-GCM biases, and quality issues of radar precipitation observations, which are relevant to the wider regional climate modelling community.

Dynamic system simulation of small satellite projects

NASA Astrophysics Data System (ADS)

Raif, Matthias; Walter, Ulrich; Bouwmeester, Jasper

2010-11-01

A prerequisite to accomplish a system simulation is to have a system model holding all necessary project information in a centralized repository that can be accessed and edited by all parties involved. At the Institute of Astronautics of the Technische Universitaet Muenchen a modular approach for modeling and dynamic simulation of satellite systems has been developed called dynamic system simulation (DySyS). DySyS is based on the platform independent description language SysML to model a small satellite project with respect to the system composition and dynamic behavior. A library of specific building blocks and possible relations between these blocks have been developed. From this library a system model of the satellite of interest can be created. A mapping into a C++ simulation allows the creation of an executable system model with which simulations are performed to observe the dynamic behavior of the satellite. In this paper DySyS is used to model and simulate the dynamic behavior of small satellites, because small satellite projects can act as a precursor to demonstrate the feasibility of a system model since they are less complex compared to a large scale satellite project.

Large scale rigidity-based flexibility analysis of biomolecules

PubMed Central

Streinu, Ileana

2016-01-01

KINematics And RIgidity (KINARI) is an on-going project for in silico flexibility analysis of proteins. The new version of the software, Kinari-2, extends the functionality of our free web server KinariWeb, incorporates advanced web technologies, emphasizes the reproducibility of its experiments, and makes substantially improved tools available to the user. It is designed specifically for large scale experiments, in particular, for (a) very large molecules, including bioassemblies with high degree of symmetry such as viruses and crystals, (b) large collections of related biomolecules, such as those obtained through simulated dilutions, mutations, or conformational changes from various types of dynamics simulations, and (c) is intended to work as seemlessly as possible on the large, idiosyncratic, publicly available repository of biomolecules, the Protein Data Bank. We describe the system design, along with the main data processing, computational, mathematical, and validation challenges underlying this phase of the KINARI project. PMID:26958583

A Center for Extraterrestrial Engineering and Construction (CETEC)

NASA Technical Reports Server (NTRS)

Leigh, Gerald G.

1992-01-01

A group of knowledgeable scientists and engineers in New Mexico has recognized the need for such a testing capability and has proposed a project to evelop an extraterrestrial surface simulation facility. A group of universities, national laboratories, and private industrial firms is proposing to establish a Center for Extraterrestrial Engineering and Construction (CETEC) and to develop large extraterrestrial surface simulation facilities in which this needed testing can be realistically performed. The CETEC is envisioned to be both a center of knowledge and data regarding engineering, construction, mining, and material process operations on extraterrestrial bodies and a set of extraterrestrial surface simulation facilities. The primary CETEC facility is proposed to be a large domed building made of steel reinforced concrete with more than one acre of test floor area covered with several feet of simulated lunar soil and dust. Various aspects of the project are presented in viewgraph form.

Cloud radiative effects and changes simulated by the Coupled Model Intercomparison Project Phase 5 models

NASA Astrophysics Data System (ADS)

Shin, Sun-Hee; Kim, Ok-Yeon; Kim, Dongmin; Lee, Myong-In

2017-07-01

Using 32 CMIP5 (Coupled Model Intercomparison Project Phase 5) models, this study examines the veracity in the simulation of cloud amount and their radiative effects (CREs) in the historical run driven by observed external radiative forcing for 1850-2005, and their future changes in the RCP (Representative Concentration Pathway) 4.5 scenario runs for 2006-2100. Validation metrics for the historical run are designed to examine the accuracy in the representation of spatial patterns for climatological mean, and annual and interannual variations of clouds and CREs. The models show large spread in the simulation of cloud amounts, specifically in the low cloud amount. The observed relationship between cloud amount and the controlling large-scale environment are also reproduced diversely by various models. Based on the validation metrics, four models—ACCESS1.0, ACCESS1.3, HadGEM2-CC, and HadGEM2-ES—are selected as best models, and the average of the four models performs more skillfully than the multimodel ensemble average. All models project global-mean SST warming at the increase of the greenhouse gases, but the magnitude varies across the simulations between 1 and 2 K, which is largely attributable to the difference in the change of cloud amount and distribution. The models that simulate more SST warming show a greater increase in the net CRE due to reduced low cloud and increased incoming shortwave radiation, particularly over the regions of marine boundary layer in the subtropics. Selected best-performing models project a significant reduction in global-mean cloud amount of about -0.99% K-1 and net radiative warming of 0.46 W m-2 K-1, suggesting a role of positive feedback to global warming.

SIERRA Low Mach Module: Fuego User Manual Version 4.46.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sierra Thermal/Fluid Team

2017-09-01

The SIERRA Low Mach Module: Fuego along with the SIERRA Participating Media Radiation Module: Syrinx, henceforth referred to as Fuego and Syrinx, respectively, are the key elements of the ASCI fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Syrinx represents the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the coremore » architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.« less

SIERRA Low Mach Module: Fuego Theory Manual Version 4.44

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sierra Thermal /Fluid Team

2017-04-01

The SIERRA Low Mach Module: Fuego along with the SIERRA Participating Media Radiation Module: Syrinx, henceforth referred to as Fuego and Syrinx, respectively, are the key elements of the ASCI fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Syrinx represents the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the coremore » architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.« less

SIERRA Low Mach Module: Fuego Theory Manual Version 4.46.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sierra Thermal/Fluid Team

The SIERRA Low Mach Module: Fuego along with the SIERRA Participating Media Radiation Module: Syrinx, henceforth referred to as Fuego and Syrinx, respectively, are the key elements of the ASCI fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Syrinx represents the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the coremore » architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.« less

Efficient coarse simulation of a growing avascular tumor

PubMed Central

Kavousanakis, Michail E.; Liu, Ping; Boudouvis, Andreas G.; Lowengrub, John; Kevrekidis, Ioannis G.

2013-01-01

The subject of this work is the development and implementation of algorithms which accelerate the simulation of early stage tumor growth models. Among the different computational approaches used for the simulation of tumor progression, discrete stochastic models (e.g., cellular automata) have been widely used to describe processes occurring at the cell and subcell scales (e.g., cell-cell interactions and signaling processes). To describe macroscopic characteristics (e.g., morphology) of growing tumors, large numbers of interacting cells must be simulated. However, the high computational demands of stochastic models make the simulation of large-scale systems impractical. Alternatively, continuum models, which can describe behavior at the tumor scale, often rely on phenomenological assumptions in place of rigorous upscaling of microscopic models. This limits their predictive power. In this work, we circumvent the derivation of closed macroscopic equations for the growing cancer cell populations; instead, we construct, based on the so-called “equation-free” framework, a computational superstructure, which wraps around the individual-based cell-level simulator and accelerates the computations required for the study of the long-time behavior of systems involving many interacting cells. The microscopic model, e.g., a cellular automaton, which simulates the evolution of cancer cell populations, is executed for relatively short time intervals, at the end of which coarse-scale information is obtained. These coarse variables evolve on slower time scales than each individual cell in the population, enabling the application of forward projection schemes, which extrapolate their values at later times. This technique is referred to as coarse projective integration. Increasing the ratio of projection times to microscopic simulator execution times enhances the computational savings. Crucial accuracy issues arising for growing tumors with radial symmetry are addressed by applying the coarse projective integration scheme in a cotraveling (cogrowing) frame. As a proof of principle, we demonstrate that the application of this scheme yields highly accurate solutions, while preserving the computational savings of coarse projective integration. PMID:22587128

The western Pacific monsoon in CMIP5 models: Model evaluation and projections

NASA Astrophysics Data System (ADS)

Brown, Josephine R.; Colman, Robert A.; Moise, Aurel F.; Smith, Ian N.

2013-11-01

ability of 35 models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) to simulate the western Pacific (WP) monsoon is evaluated over four representative regions around Timor, New Guinea, the Solomon Islands and Palau. Coupled model simulations are compared with atmosphere-only model simulations (with observed sea surface temperatures, SSTs) to determine the impact of SST biases on model performance. Overall, the CMIP5 models simulate the WP monsoon better than previous-generation Coupled Model Intercomparison Project Phase 3 (CMIP3) models, but some systematic biases remain. The atmosphere-only models are better able to simulate the seasonal cycle of zonal winds than the coupled models, but display comparable biases in the rainfall. The CMIP5 models are able to capture features of interannual variability in response to the El Niño-Southern Oscillation. In climate projections under the RCP8.5 scenario, monsoon rainfall is increased over most of the WP monsoon domain, while wind changes are small. Widespread rainfall increases at low latitudes in the summer hemisphere appear robust as a large majority of models agree on the sign of the change. There is less agreement on rainfall changes in winter. Interannual variability of monsoon wet season rainfall is increased in a warmer climate, particularly over Palau, Timor and the Solomon Islands. A subset of the models showing greatest skill in the current climate confirms the overall projections, although showing markedly smaller rainfall increases in the western equatorial Pacific. The changes found here may have large impacts on Pacific island countries influenced by the WP monsoon.

Large-Eddy Simulation of the Impact of Great Garuda Project on Wind and Thermal Environment over Built-Up Area in Jakarta

NASA Astrophysics Data System (ADS)

Yucel, M.; Sueishi, T.; Inagaki, A.; Kanda, M.

2017-12-01

`Great Garuda' project is an eagle-shaped offshore structure with 17 artificial islands. This project has been designed for the coastal protection and land reclamation of Jakarta due to catastrophic flooding in the city. It offers an urban generation for 300.000 inhabitants and 600.000 workers in addition to its water safety goal. A broad coalition of Indonesian scientists has criticized the project for being negative impacts on the surrounding environment. Despite the vast research by Indonesian scientist on maritime environment, studies on wind and thermal environment over built-up area are still lacking. However, the construction of the various islands off the coast may result changes in wind patterns and thermal environment due to the alteration of the coastline and urbanization in the Jakarta Bay. Therefore, it is important to understand the airflow within the urban canopy in case of unpredictable gust events. These gust events may occur through the closely-packed high-rise buildings and pedestrians may be harmed from such gusts. Accordingly, we used numerical simulations to investigate the impact of the sea wall and the artificial islands over built-up area and, the intensity of wind gusts at the pedestrian level. Considering the fact that the size of turbulence organized structure sufficiently large computational domain is required. Therefore, a 19.2km×4.8km×1.0 km simulation domain with 2-m resolution in all directions was created to explicitly resolve the detailed shapes of buildings and the flow at the pedestrian level. This complex computation was accomplished by implementing a large-eddy simulation (LES) model. Two case studies were conducted considering the effect of realistic surface roughness and upward heat flux. Case_1 was conducted based on the current built environment and Case_2 for investigating the effect of the project on the chosen coastal region of the city. Fig.1 illustrates the schematic of the large-eddy simulation domains of two cases with and without Great Garuda Sea Wall and 17 artificial islands. 3D model of Great Garuda is shown in Fig.2. In addition to the cases mentioned above, the simulation will be generated assigning more realistic heat flux outputs from energy balance model and, inflow boundary conditions coupling with mesoscale model (Weather Research and Forecast model).

Simulator Network Project Report: A tool for improvement of teaching materials and targeted resource usage in Skills Labs

PubMed Central

Damanakis, Alexander; Blaum, Wolf E.; Stosch, Christoph; Lauener, Hansjörg; Richter, Sabine; Schnabel, Kai P.

2013-01-01

During the last decade, medical education in the German-speaking world has been striving to become more practice-oriented. This is currently being achieved in many schools through the implementation of simulation-based instruction in Skills Labs. Simulators are thus an essential part of this type of medical training, and their acquisition and operation by a Skills Lab require a large outlay of resources. Therefore, the Practical Skills Committee of the Medical Education Society (GMA) introduced a new project, which aims to improve the flow of information between the Skills Labs and enable a transparent assessment of the simulators via an online database (the Simulator Network). PMID:23467581

Development of carbon response trajectories using FIA plot data and FVS growth simulator: challenges of a large scale simulation project

Treesearch

James B. McCarter; Sean Healey

2015-01-01

The Forest Carbon Management Framework (ForCaMF) integrates Forest Inventory and Analysis (FIA) plot inventory data, disturbance histories, and carbon response trajectories to develop estimates of disturbance and management effects on carbon pools for the National Forest System. All appropriate FIA inventory plots are simulated using the Forest Vegetation Simulator (...

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Gang

Mid-latitude extreme weather events are responsible for a large part of climate-related damage. Yet large uncertainties remain in climate model projections of heat waves, droughts, and heavy rain/snow events on regional scales, limiting our ability to effectively use these projections for climate adaptation and mitigation. These uncertainties can be attributed to both the lack of spatial resolution in the models, and to the lack of a dynamical understanding of these extremes. The approach of this project is to relate the fine-scale features to the large scales in current climate simulations, seasonal re-forecasts, and climate change projections in a very widemore » range of models, including the atmospheric and coupled models of ECMWF over a range of horizontal resolutions (125 to 10 km), aqua-planet configuration of the Model for Prediction Across Scales and High Order Method Modeling Environments (resolutions ranging from 240 km – 7.5 km) with various physics suites, and selected CMIP5 model simulations. The large scale circulation will be quantified both on the basis of the well tested preferred circulation regime approach, and very recently developed measures, the finite amplitude Wave Activity (FAWA) and its spectrum. The fine scale structures related to extremes will be diagnosed following the latest approaches in the literature. The goal is to use the large scale measures as indicators of the probability of occurrence of the finer scale structures, and hence extreme events. These indicators will then be applied to the CMIP5 models and time-slice projections of a future climate.« less

Signal to noise quantification of regional climate projections

NASA Astrophysics Data System (ADS)

Li, S.; Rupp, D. E.; Mote, P.

2016-12-01

One of the biggest challenges in interpreting climate model outputs for impacts studies and adaptation planning is understanding the sources of disagreement among models (which is often used imperfectly as a stand-in for system uncertainty). Internal variability is a primary source of uncertainty in climate projections, especially for precipitation, for which models disagree about even the sign of changes in large areas like the continental US. Taking advantage of a large initial-condition ensemble of regional climate simulations, this study quantifies the magnitude of changes forced by increasing greenhouse gas concentrations relative to internal variability. Results come from a large initial-condition ensemble of regional climate model simulations generated by weather@home, a citizen science computing platform, where the western United States climate was simulated for the recent past (1985-2014) and future (2030-2059) using a 25-km horizontal resolution regional climate model (HadRM3P) nested in global atmospheric model (HadAM3P). We quantify grid point level signal-to-noise not just in temperature and precipitation responses, but also the energy and moisture flux terms that are related to temperature and precipitation responses, to provide important insights regarding uncertainty in climate change projections at local and regional scales. These results will aid modelers in determining appropriate ensemble sizes for different climate variables and help users of climate model output with interpreting climate model projections.

Using a Computerized Classroom Simulation to Prepare Pre-Service Teachers

ERIC Educational Resources Information Center

McPherson, Rebekah; Tyler-Wood, Tandra; McEnturff Ellison, Amber; Peak, Pamela

2011-01-01

This study at a large midwestern university evaluated the use of a web-based simulated classroom, simSchool, with pre-service and in-service special education students, to determine if use of the simulated classroom influences students' perceptions of inclusion and teacher preparation. The project used a nonequivalent comparison group,…

Field scale test of multi-dimensional flow and morphodynamic simulations used for restoration design analysis

USGS Publications Warehouse

McDonald, Richard R.; Nelson, Jonathan M.; Fosness, Ryan L.; Nelson, Peter O.; Constantinescu, George; Garcia, Marcelo H.; Hanes, Dan

2016-01-01

Two- and three-dimensional morphodynamic simulations are becoming common in studies of channel form and process. The performance of these simulations are often validated against measurements from laboratory studies. Collecting channel change information in natural settings for model validation is difficult because it can be expensive and under most channel forming flows the resulting channel change is generally small. Several channel restoration projects designed in part to armor large meanders with several large spurs constructed of wooden piles on the Kootenai River, ID, have resulted in rapid bed elevation change following construction. Monitoring of these restoration projects includes post- restoration (as-built) Digital Elevation Models (DEMs) as well as additional channel surveys following high channel forming flows post-construction. The resulting sequence of measured bathymetry provides excellent validation data for morphodynamic simulations at the reach scale of a real river. In this paper we test the performance a quasi-three-dimensional morphodynamic simulation against the measured elevation change. The resulting simulations predict the pattern of channel change reasonably well but many of the details such as the maximum scour are under predicted.

The APOSTLE project: Local Group kinematic mass constraints and simulation candidate selection

NASA Astrophysics Data System (ADS)

Fattahi, Azadeh; Navarro, Julio F.; Sawala, Till; Frenk, Carlos S.; Oman, Kyle A.; Crain, Robert A.; Furlong, Michelle; Schaller, Matthieu; Schaye, Joop; Theuns, Tom; Jenkins, Adrian

2016-03-01

We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group (LG) of galaxies. We find, in agreement with the `timing argument' and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of ˜5 × 1012 M⊙. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of ˜2 times smaller than suggested by the timing argument, but with large dispersion. Guided by these results, we have selected 12 halo pairs with total mass in the range 1.6-3.6 × 1012 M⊙ for the APOSTLE project (A Project Of Simulating The Local Environment), a suite of hydrodynamical resimulations at various numerical resolution levels (reaching up to ˜104 M⊙ per gas particle) that use the subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the Λ cold dark matter cosmology with observations of our local Universe.

Projection Effects of Large-scale Structures on Weak-lensing Peak Abundances

NASA Astrophysics Data System (ADS)

Yuan, Shuo; Liu, Xiangkun; Pan, Chuzhong; Wang, Qiao; Fan, Zuhui

2018-04-01

High peaks in weak lensing (WL) maps originate dominantly from the lensing effects of single massive halos. Their abundance is therefore closely related to the halo mass function and thus a powerful cosmological probe. However, besides individual massive halos, large-scale structures (LSS) along lines of sight also contribute to the peak signals. In this paper, with ray-tracing simulations, we investigate the LSS projection effects. We show that for current surveys with a large shape noise, the stochastic LSS effects are subdominant. For future WL surveys with source galaxies having a median redshift z med ∼ 1 or higher, however, they are significant. For the cosmological constraints derived from observed WL high-peak counts, severe biases can occur if the LSS effects are not taken into account properly. We extend the model of Fan et al. by incorporating the LSS projection effects into the theoretical considerations. By comparing with simulation results, we demonstrate the good performance of the improved model and its applicability in cosmological studies.

Why build a virtual brain? Large-scale neural simulations as jump start for cognitive computing

NASA Astrophysics Data System (ADS)

Colombo, Matteo

2017-03-01

Despite the impressive amount of financial resources recently invested in carrying out large-scale brain simulations, it is controversial what the pay-offs are of pursuing this project. One idea is that from designing, building, and running a large-scale neural simulation, scientists acquire knowledge about the computational performance of the simulating system, rather than about the neurobiological system represented in the simulation. It has been claimed that this knowledge may usher in a new era of neuromorphic, cognitive computing systems. This study elucidates this claim and argues that the main challenge this era is facing is not the lack of biological realism. The challenge lies in identifying general neurocomputational principles for the design of artificial systems, which could display the robust flexibility characteristic of biological intelligence.

Practical recipes for the model order reduction, dynamical simulation and compressive sampling of large-scale open quantum systems

NASA Astrophysics Data System (ADS)

Sidles, John A.; Garbini, Joseph L.; Harrell, Lee E.; Hero, Alfred O.; Jacky, Jonathan P.; Malcomb, Joseph R.; Norman, Anthony G.; Williamson, Austin M.

2009-06-01

Practical recipes are presented for simulating high-temperature and nonequilibrium quantum spin systems that are continuously measured and controlled. The notion of a spin system is broadly conceived, in order to encompass macroscopic test masses as the limiting case of large-j spins. The simulation technique has three stages: first the deliberate introduction of noise into the simulation, then the conversion of that noise into an equivalent continuous measurement and control process, and finally, projection of the trajectory onto state-space manifolds having reduced dimensionality and possessing a Kähler potential of multilinear algebraic form. These state-spaces can be regarded as ruled algebraic varieties upon which a projective quantum model order reduction (MOR) is performed. The Riemannian sectional curvature of ruled Kählerian varieties is analyzed, and proved to be non-positive upon all sections that contain a rule. These manifolds are shown to contain Slater determinants as a special case and their identity with Grassmannian varieties is demonstrated. The resulting simulation formalism is used to construct a positive P-representation for the thermal density matrix. Single-spin detection by magnetic resonance force microscopy (MRFM) is simulated, and the data statistics are shown to be those of a random telegraph signal with additive white noise. Larger-scale spin-dust models are simulated, having no spatial symmetry and no spatial ordering; the high-fidelity projection of numerically computed quantum trajectories onto low dimensionality Kähler state-space manifolds is demonstrated. The reconstruction of quantum trajectories from sparse random projections is demonstrated, the onset of Donoho-Stodden breakdown at the Candès-Tao sparsity limit is observed, a deterministic construction for sampling matrices is given and methods for quantum state optimization by Dantzig selection are given.

Innovative Instruction in Higher Education: Thirty Exemplary Projects Conducted in Selected Institutions of Post-Secondary Education--State of Oregon.

ERIC Educational Resources Information Center

Harper, Ronald; And Others

This manual reviews thirty projects selected by the Oregon Educational Coordinating Council (ECC) as exemplary in method, operation, and development. The projects are organized into 9 broad classifications: large group-small group alternatives, autotutorial programmed instruction, process centered, computer and simulation, on-site/field study,…

Indian LSSC (Large Space Simulation Chamber) facility

NASA Technical Reports Server (NTRS)

Brar, A. S.; Prasadarao, V. S.; Gambhir, R. D.; Chandramouli, M.

1988-01-01

The Indian Space Agency has undertaken a major project to acquire in-house capability for thermal and vacuum testing of large satellites. This Large Space Simulation Chamber (LSSC) facility will be located in Bangalore and is to be operational in 1989. The facility is capable of providing 4 meter diameter solar simulation with provision to expand to 4.5 meter diameter at a later date. With such provisions as controlled variations of shroud temperatures and availability of infrared equipment as alternative sources of thermal radiation, this facility will be amongst the finest anywhere. The major design concept and major aspects of the LSSC facility are presented here.

Projected Future Vegetation Changes for the Northwest United States and Southwest Canada at a Fine Spatial Resolution Using a Dynamic Global Vegetation Model.

PubMed

Shafer, Sarah L; Bartlein, Patrick J; Gray, Elizabeth M; Pelltier, Richard T

2015-01-01

Future climate change may significantly alter the distributions of many plant taxa. The effects of climate change may be particularly large in mountainous regions where climate can vary significantly with elevation. Understanding potential future vegetation changes in these regions requires methods that can resolve vegetation responses to climate change at fine spatial resolutions. We used LPJ, a dynamic global vegetation model, to assess potential future vegetation changes for a large topographically complex area of the northwest United States and southwest Canada (38.0-58.0°N latitude by 136.6-103.0°W longitude). LPJ is a process-based vegetation model that mechanistically simulates the effect of changing climate and atmospheric CO2 concentrations on vegetation. It was developed and has been mostly applied at spatial resolutions of 10-minutes or coarser. In this study, we used LPJ at a 30-second (~1-km) spatial resolution to simulate potential vegetation changes for 2070-2099. LPJ was run using downscaled future climate simulations from five coupled atmosphere-ocean general circulation models (CCSM3, CGCM3.1(T47), GISS-ER, MIROC3.2(medres), UKMO-HadCM3) produced using the A2 greenhouse gases emissions scenario. Under projected future climate and atmospheric CO2 concentrations, the simulated vegetation changes result in the contraction of alpine, shrub-steppe, and xeric shrub vegetation across the study area and the expansion of woodland and forest vegetation. Large areas of maritime cool forest and cold forest are simulated to persist under projected future conditions. The fine spatial-scale vegetation simulations resolve patterns of vegetation change that are not visible at coarser resolutions and these fine-scale patterns are particularly important for understanding potential future vegetation changes in topographically complex areas.

Australian snowpack in the NARCliM ensemble: evaluation, bias correction and future projections

NASA Astrophysics Data System (ADS)

Luca, Alejandro Di; Evans, Jason P.; Ji, Fei

2017-10-01

In this study we evaluate the ability of an ensemble of high-resolution Regional Climate Model simulations to represent snow cover characteristics over the Australian Alps and go on to asses future projections of snowpack characteristics. Our results show that the ensemble presents a cold temperature bias and overestimates total precipitation leading to a general overestimation of the snow cover as compared with MODIS satellite data. We then produce a new set of snowpack characteristics by running a temperature based snow melt/accumulation model forced by bias corrected temperature and precipitation fields. While some positive snow cover biases remain, the bias corrected (BC) dataset show large improvements regarding the simulation of total amounts, seasonality and spatial distribution of the snow cover compared with MODIS products. Both the raw and BC datasets are then used to assess future changes in the snowpack characteristics. Both datasets show robust increases in near-surface temperatures and decreases in snowfall that lead to a substantial reduction of the snowpack over the Australian Alps. The snowpack decreases by about 15 and 60% by 2030 and 2070 respectively. While the BC data introduce large differences in the simulation of the present climate snowpack, in relative terms future changes appear to be similar to those obtained using the raw data. Future temperature projections show a clear dependence with elevation through the snow-albedo feedback effect that affects snowpack projections. Uncertainties in future projections of the snowpack are large in both datasets and are mainly dominated by the choice of the lateral boundary conditions.

Expert systems and simulation models; Proceedings of the Seminar, Tucson, AZ, November 18, 19, 1985

NASA Technical Reports Server (NTRS)

1986-01-01

The seminar presents papers on modeling and simulation methodology, artificial intelligence and expert systems, environments for simulation/expert system development, and methodology for simulation/expert system development. Particular attention is given to simulation modeling concepts and their representation, modular hierarchical model specification, knowledge representation, and rule-based diagnostic expert system development. Other topics include the combination of symbolic and discrete event simulation, real time inferencing, and the management of large knowledge-based simulation projects.

Future Climate Change in the Baltic Sea Area

NASA Astrophysics Data System (ADS)

Bøssing Christensen, Ole; Kjellström, Erik; Zorita, Eduardo; Sonnenborg, Torben; Meier, Markus; Grinsted, Aslak

2015-04-01

Regional climate models have been used extensively since the first assessment of climate change in the Baltic Sea region published in 2008, not the least for studies of Europe (and including the Baltic Sea catchment area). Therefore, conclusions regarding climate model results have a better foundation than was the case for the first BACC report of 2008. This presentation will report model results regarding future climate. What is the state of understanding about future human-driven climate change? We will cover regional models, statistical downscaling, hydrological modelling, ocean modelling and sea-level change as it is projected for the Baltic Sea region. Collections of regional model simulations from the ENSEMBLES project for example, financed through the European 5th Framework Programme and the World Climate Research Programme Coordinated Regional Climate Downscaling Experiment, have made it possible to obtain an increasingly robust estimation of model uncertainty. While the first Baltic Sea assessment mainly used four simulations from the European 5th Framework Programme PRUDENCE project, an ensemble of 13 transient regional simulations with twice the horizontal resolution reaching the end of the 21st century has been available from the ENSEMBLES project; therefore it has been possible to obtain more quantitative assessments of model uncertainty. The literature about future climate change in the Baltic Sea region is largely built upon the ENSEMBLES project. Also within statistical downscaling, a considerable number of papers have been published, encompassing now the application of non-linear statistical models, projected changes in extremes and correction of climate model biases. The uncertainty of hydrological change has received increasing attention since the previous Baltic Sea assessment. Several studies on the propagation of uncertainties originating in GCMs, RCMs, and emission scenarios are presented. The number of studies on uncertainties related to downscaling and impact models is relatively small, but more are emerging. A large number of coupled climate-environmental scenario simulations for the Baltic Sea have been performed within the BONUS+ projects (ECOSUPPORT, INFLOW, AMBER and Baltic-C (2009-2011)), using various combinations of output from GCMs, RCMs, hydrological models and scenarios for load and emission of nutrients as forcing for Baltic Sea models. Such a large ensemble of scenario simulations for the Baltic Sea has never before been produced and enables for the first time an estimation of uncertainties.

Large Scale Simulation Platform for NODES Validation Study

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sotorrio, P.; Qin, Y.; Min, L.

2017-04-27

This report summarizes the Large Scale (LS) simulation platform created for the Eaton NODES project. The simulation environment consists of both wholesale market simulator and distribution simulator and includes the CAISO wholesale market model and a PG&E footprint of 25-75 feeders to validate the scalability under a scenario of 33% RPS in California with additional 17% of DERS coming from distribution and customers. The simulator can generate hourly unit commitment, 5-minute economic dispatch, and 4-second AGC regulation signals. The simulator is also capable of simulating greater than 10k individual controllable devices. Simulated DERs include water heaters, EVs, residential and lightmore » commercial HVAC/buildings, and residential-level battery storage. Feeder-level voltage regulators and capacitor banks are also simulated for feeder-level real and reactive power management and Vol/Var control.« less

System dynamic simulation: A new method in social impact assessment (SIA)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Karami, Shobeir, E-mail: shobeirkarami@gmail.com; Karami, Ezatollah, E-mail: ekarami@shirazu.ac.ir; Buys, Laurie, E-mail: l.buys@qut.edu.au

Many complex social questions are difficult to address adequately with conventional methods and techniques, due to the complicated dynamics, and hard to quantify social processes. Despite these difficulties researchers and practitioners have attempted to use conventional methods not only in evaluative modes but also in predictive modes to inform decision making. The effectiveness of SIAs would be increased if they were used to support the project design processes. This requires deliberate use of lessons from retrospective assessments to inform predictive assessments. Social simulations may be a useful tool for developing a predictive SIA method. There have been limited attempts tomore » develop computer simulations that allow social impacts to be explored and understood before implementing development projects. In light of this argument, this paper aims to introduce system dynamic (SD) simulation as a new predictive SIA method in large development projects. We propose the potential value of the SD approach to simulate social impacts of development projects. We use data from the SIA of Gareh-Bygone floodwater spreading project to illustrate the potential of SD simulation in SIA. It was concluded that in comparison to traditional SIA methods SD simulation can integrate quantitative and qualitative inputs from different sources and methods and provides a more effective and dynamic assessment of social impacts for development projects. We recommend future research to investigate the full potential of SD in SIA in comparing different situations and scenarios.« less

Enabling Predictive Simulation and UQ of Complex Multiphysics PDE Systems by the Development of Goal-Oriented Variational Sensitivity Analysis and a-Posteriori Error Estimation Methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Estep, Donald

2015-11-30

This project addressed the challenge of predictive computational analysis of strongly coupled, highly nonlinear multiphysics systems characterized by multiple physical phenomena that span a large range of length- and time-scales. Specifically, the project was focused on computational estimation of numerical error and sensitivity analysis of computational solutions with respect to variations in parameters and data. In addition, the project investigated the use of accurate computational estimates to guide efficient adaptive discretization. The project developed, analyzed and evaluated new variational adjoint-based techniques for integration, model, and data error estimation/control and sensitivity analysis, in evolutionary multiphysics multiscale simulations.

Commnity Petascale Project for Accelerator Science And Simulation: Advancing Computational Science for Future Accelerators And Accelerator Technologies

DOE Office of Scientific and Technical Information (OSTI.GOV)

Spentzouris, Panagiotis; /Fermilab; Cary, John

The design and performance optimization of particle accelerators are essential for the success of the DOE scientific program in the next decade. Particle accelerators are very complex systems whose accurate description involves a large number of degrees of freedom and requires the inclusion of many physics processes. Building on the success of the SciDAC-1 Accelerator Science and Technology project, the SciDAC-2 Community Petascale Project for Accelerator Science and Simulation (ComPASS) is developing a comprehensive set of interoperable components for beam dynamics, electromagnetics, electron cooling, and laser/plasma acceleration modelling. ComPASS is providing accelerator scientists the tools required to enable the necessarymore » accelerator simulation paradigm shift from high-fidelity single physics process modeling (covered under SciDAC1) to high-fidelity multiphysics modeling. Our computational frameworks have been used to model the behavior of a large number of accelerators and accelerator R&D experiments, assisting both their design and performance optimization. As parallel computational applications, the ComPASS codes have been shown to make effective use of thousands of processors.« less

Discrete event simulation of NASA's Remote Exploration and Experimentation Project (REE)

NASA Technical Reports Server (NTRS)

Dunphy, J.; Rogstad, S.

2001-01-01

The Remote Exploration and Experimentation Project (REE) is a new initiative at JPL to be able to place a supercomputer on board a spacecraft and allow large amounts of data reduction and compression to be done before science results are returned to Earth.

A Full-Color, High-Resolution Laser Projector for a Flight Simulator Visual Display

DTIC Science & Technology

1993-08-01

current projection technology. Laser projection promises an increased color gamut , higher luminance, zero persistence, and increased line rate. The...1 7 List of Figures Fig. I Color Gamut of Laser Projection Display...primary wavelengths plot on the boundary of the ICI (CIE) chromaticity diagram; a characteristic that allows a very large gamut of colors. Unlike

Simulation based evaluation of the designs of the Advanced Gamma-ray Imageing System (AGIS)

NASA Astrophysics Data System (ADS)

Bugaev, Slava; Buckley, James; Digel, Seth; Funk, Stephen; Konopelko, Alex; Krawczynski, Henric; Lebohec, Steohan; Maier, Gernot; Vassiliev, Vladimir

2009-05-01

The AGIS project under design study, is a large array of imaging atmospheric Cherenkov telescopes for gamma-rays astronomy between 40GeV and 100 TeV. In this paper we present the ongoing simulation effort to model the considered design approaches as a function of the main parameters such as array geometry, telescope optics and camera design in such a way the gamma ray observation capabilities can be optimized against the overall project cost.

Development of a database system for near-future climate change projections under the Japanese National Project SI-CAT

NASA Astrophysics Data System (ADS)

Nakagawa, Y.; Kawahara, S.; Araki, F.; Matsuoka, D.; Ishikawa, Y.; Fujita, M.; Sugimoto, S.; Okada, Y.; Kawazoe, S.; Watanabe, S.; Ishii, M.; Mizuta, R.; Murata, A.; Kawase, H.

2017-12-01

Analyses of large ensemble data are quite useful in order to produce probabilistic effect projection of climate change. Ensemble data of "+2K future climate simulations" are currently produced by Japanese national project "Social Implementation Program on Climate Change Adaptation Technology (SI-CAT)" as a part of a database for Policy Decision making for Future climate change (d4PDF; Mizuta et al. 2016) produced by Program for Risk Information on Climate Change. Those data consist of global warming simulations and regional downscaling simulations. Considering that those data volumes are too large (a few petabyte) to download to a local computer of users, a user-friendly system is required to search and download data which satisfy requests of the users. We develop "a database system for near-future climate change projections" for providing functions to find necessary data for the users under SI-CAT. The database system for near-future climate change projections mainly consists of a relational database, a data download function and user interface. The relational database using PostgreSQL is a key function among them. Temporally and spatially compressed data are registered on the relational database. As a first step, we develop the relational database for precipitation, temperature and track data of typhoon according to requests by SI-CAT members. The data download function using Open-source Project for a Network Data Access Protocol (OPeNDAP) provides a function to download temporally and spatially extracted data based on search results obtained by the relational database. We also develop the web-based user interface for using the relational database and the data download function. A prototype of the database system for near-future climate change projections are currently in operational test on our local server. The database system for near-future climate change projections will be released on Data Integration and Analysis System Program (DIAS) in fiscal year 2017. Techniques of the database system for near-future climate change projections might be quite useful for simulation and observational data in other research fields. We report current status of development and some case studies of the database system for near-future climate change projections.

Virtual reality to simulate large lighting with high efficiency LEDs

NASA Astrophysics Data System (ADS)

Blandet, Thierry; Coutelier, Gilles; Meyrueis, Patrick

2011-05-01

When a city or a local authority wishes to emphasize its historical heritage, for the lighting of its streets, setting up lights during the festive season, they call upon the skills of a lighting designer. The lighting designer proposes concepts, ideas, lighting, and to be able to present them, he makes use of simulation. On the other hand lighting technologies are evolving very rapidly and new lighting systems offer features that lighting designers are now integrating their projects. The street lights consume lot of energy; light projects are now taking into account the energy saving aspect. Lighting systems based on LEDs today provide good lighting needs, taking into account sustainable development issues while enabling new creative dimension. The lighting simulation can handle these parameters. Images or video simulation are no longer sufficient: stereoscopy and virtual reality techniques allow better communication and better understanding of projects. Virtual reality offers new possibilities of interaction, the freedom of movement in a scene, the presentation of variants or interactive simulations.

Future changes in the climatology of the Great Plains low-level jet derived from fine resolution multi-model simulations.

PubMed

Tang, Ying; Winkler, Julie; Zhong, Shiyuan; Bian, Xindi; Doubler, Dana; Yu, Lejiang; Walters, Claudia

2017-07-10

The southerly Great Plains low-level jet (GPLLJ) is one of the most significant circulation features of the central U.S. linking large-scale atmospheric circulation with the regional climate. GPLLJs transport heat and moisture, contribute to thunderstorm and severe weather formation, provide a corridor for the springtime migration of birds and insects, enhance wind energy availability, and disperse air pollution. We assess future changes in GPLLJ frequency using an eight member ensemble of dynamically-downscaled climate simulations for the mid-21st century. Nocturnal GPLLJ frequency is projected to increase in the southern plains in spring and in the central plains in summer, whereas current climatological patterns persist into the future for daytime and cool season GPLLJs. The relationship between future GPLLJ frequency and the extent and strength of anticyclonic airflow over eastern North America varies with season. Most simulations project a westward shift of anticyclonic airflow in summer, but uncertainty is larger for spring with only half of the simulations suggesting a westward expansion. The choice of regional climate model and the driving lateral boundary conditions have a large influence on the projected future changes in GPLLJ frequency and highlight the importance of multi-model ensembles to estimate the uncertainty surrounding the future GPLLJ climatology.

Computer-aided software development process design

NASA Technical Reports Server (NTRS)

Lin, Chi Y.; Levary, Reuven R.

1989-01-01

The authors describe an intelligent tool designed to aid managers of software development projects in planning, managing, and controlling the development process of medium- to large-scale software projects. Its purpose is to reduce uncertainties in the budget, personnel, and schedule planning of software development projects. It is based on dynamic model for the software development and maintenance life-cycle process. This dynamic process is composed of a number of time-varying, interacting developmental phases, each characterized by its intended functions and requirements. System dynamics is used as a modeling methodology. The resulting Software LIfe-Cycle Simulator (SLICS) and the hybrid expert simulation system of which it is a subsystem are described.

Modern (1992–2011) and projected (2012–99) peak snowpack and May–July runoff for the Fort Peck Lake and Lake Sakakawea watersheds in the Upper Missouri River Basin

USGS Publications Warehouse

Stamm, John F.; Todey, Dennis; Mayes Bousted, Barbara; Rossi, Shawn; Norton, Parker A.; Carter, Janet M.

2016-02-09

Annual peak snowpack was projected to have a downward trend for the Fort Peck Lake watershed and an upward trend for the lower Lake Sakakawea watershed. Projections of May–July runoff had a significant downward trend for the Fort Peck Lake, lower Lake Sakakawea, and Lake Sakakawea (combination of Fort Peck Lake and lower Lake Sakakawea) watersheds. Downward trends in projected May–July runoff indicated that power production at Fort Peck Dam might be affected particularly in the later part of the simulation (2061–99); however, confidence in projected May–July runoff for the later part of the simulation was less certain because bias-corrected air temperatures from CCSM3 and CCSM4 commonly fell outside of the observed range used for calibration. Projected May–July runoff combined for the Fort Peck Lake and lower Lake Sakakawea watersheds were on the order of magnitude of the 2011 flood for 1 simulation year for each of the CCSM-based simulations. High peak snowpack and precipitation in April, May, and June in the plains was associated with large May–July runoff events; therefore, high precipitation at lower elevations in the Fort Peck Lake and lower Lake Sakakawea watersheds was a factor in the simulation of extreme runoff events at the magnitude of the 2011 flood.

Projection rule for complex-valued associative memory with large constant terms

NASA Astrophysics Data System (ADS)

Kitahara, Michimasa; Kobayashi, Masaki

Complex-valued Associative Memory (CAM) has an inherent property of rotation invariance. Rotation invariance produces many undesirable stable states and reduces the noise robustness of CAM. Constant terms may remove rotation invariance, but if the constant terms are too small, rotation invariance does not vanish. In this paper, we eliminate rotation invariance by introducing large constant terms to complex-valued neurons. We have to make constant terms sufficiently large to improve the noise robustness. We introduce a parameter to control the amplitudes of constant terms into projection rule. The large constant terms are proved to be effective by our computer simulations.

How do the methodological choices of your climate change study affect your results? A hydrologic case study across the Pacific Northwest

NASA Astrophysics Data System (ADS)

Chegwidden, O.; Nijssen, B.; Rupp, D. E.; Kao, S. C.; Clark, M. P.

2017-12-01

We describe results from a large hydrologic climate change dataset developed across the Pacific Northwestern United States and discuss how the analysis of those results can be seen as a framework for other large hydrologic ensemble investigations. This investigation will better inform future modeling efforts and large ensemble analyses across domains within and beyond the Pacific Northwest. Using outputs from the Coupled Model Intercomparison Project Phase 5 (CMIP5), we provide projections of hydrologic change for the domain through the end of the 21st century. The dataset is based upon permutations of four methodological choices: (1) ten global climate models (2) two representative concentration pathways (3) three meteorological downscaling methods and (4) four unique hydrologic model set-ups (three of which entail the same hydrologic model using independently calibrated parameter sets). All simulations were conducted across the Columbia River Basin and Pacific coastal drainages at a 1/16th ( 6 km) resolution and at a daily timestep. In total, the 172 distinct simulations offer an updated, comprehensive view of climate change projections through the end of the 21st century. The results consist of routed streamflow at 400 sites throughout the domain as well as distributed spatial fields of relevant hydrologic variables like snow water equivalent and soil moisture. In this presentation, we discuss the level of agreement with previous hydrologic projections for the study area and how these projections differ with specific methodological choices. By controlling for some methodological choices we can show how each choice affects key climatic change metrics. We discuss how the spread in results varies across hydroclimatic regimes. We will use this large dataset as a case study for distilling a wide range of hydroclimatological projections into useful climate change assessments.

Concurrent increases in wet and dry extremes projected in Texas and combined effects on groundwater

NASA Astrophysics Data System (ADS)

Yoon, Jin-Ho; Wang, S.-Y. Simon; Lo, Min-Hui; Wu, Wen-Ying

2018-05-01

The US state of Texas has experienced consecutive flooding events since spring 2015 with devastating consequences, yet these happened only a few years after the record drought of 2011. Identifying the effect of climate variability on regional water cycle extremes, such as the predicted occurrence of La Niña in winter 2017–2018 and its association with drought in Texas, remains a challenge. The present analyses use large-ensemble simulations to project the future of water cycle extremes in Texas and assess their connection with the changing El Niño–Southern Oscillation (ENSO) teleconnection under global warming. Large-ensemble simulations indicate that both intense drought and excessive precipitation are projected to increase towards the middle of the 21st century, associated with a strengthened effect from ENSO. Despite the precipitation increase projected for the southern Great Plains, groundwater storage is likely to decrease in the long run with diminishing groundwater recharge; this is due to the concurrent increases and strengthening in drought offsetting the effect of added rains. This projection provides implications to short-term climate anomaly in the face of the La Niña and to long-term water resources planning.

Acceleration techniques for dependability simulation. M.S. Thesis

NASA Technical Reports Server (NTRS)

Barnette, James David

1995-01-01

As computer systems increase in complexity, the need to project system performance from the earliest design and development stages increases. We have to employ simulation for detailed dependability studies of large systems. However, as the complexity of the simulation model increases, the time required to obtain statistically significant results also increases. This paper discusses an approach that is application independent and can be readily applied to any process-based simulation model. Topics include background on classical discrete event simulation and techniques for random variate generation and statistics gathering to support simulation.

Long-term Ozone Changes and Associated Climate Impacts in CMIP5 Simulations

NASA Technical Reports Server (NTRS)

Eyring, V.; Arblaster, J. M.; Cionni, I.; Sedlacek, J.; Perlwitz, J.; Young, P. J.; Bekki, S.; Bergmann, D.; Cameron-Smith, P.; Collins, W. J.;

Compactified cosmological simulations of the infinite universe

NASA Astrophysics Data System (ADS)

Rácz, Gábor; Szapudi, István; Csabai, István; Dobos, László

2018-06-01

We present a novel N-body simulation method that compactifies the infinite spatial extent of the Universe into a finite sphere with isotropic boundary conditions to follow the evolution of the large-scale structure. Our approach eliminates the need for periodic boundary conditions, a mere numerical convenience which is not supported by observation and which modifies the law of force on large scales in an unrealistic fashion. We demonstrate that our method outclasses standard simulations executed on workstation-scale hardware in dynamic range, it is balanced in following a comparable number of high and low k modes and, its fundamental geometry and topology match observations. Our approach is also capable of simulating an expanding, infinite universe in static coordinates with Newtonian dynamics. The price of these achievements is that most of the simulated volume has smoothly varying mass and spatial resolution, an approximation that carries different systematics than periodic simulations. Our initial implementation of the method is called StePS which stands for Stereographically projected cosmological simulations. It uses stereographic projection for space compactification and naive O(N^2) force calculation which is nevertheless faster to arrive at a correlation function of the same quality than any standard (tree or P3M) algorithm with similar spatial and mass resolution. The N2 force calculation is easy to adapt to modern graphics cards, hence our code can function as a high-speed prediction tool for modern large-scale surveys. To learn about the limits of the respective methods, we compare StePS with GADGET-2 running matching initial conditions.

A Novel Multi-scale Simulation Strategy for Turbulent Reacting Flows

DOE Office of Scientific and Technical Information (OSTI.GOV)

James, Sutherland C.

In this project, a new methodology was proposed to bridge the gap between Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES). This novel methodology, titled Lattice-Based Multiscale Simulation (LBMS), creates a lattice structure of One-Dimensional Turbulence (ODT) models. This model has been shown to capture turbulent combustion with high fidelity by fully resolving interactions between turbulence and diffusion. By creating a lattice of ODT models, which are then coupled, LBMS overcomes the shortcomings of ODT, which are its inability to capture large scale three dimensional flow structures. However, by spacing these lattices significantly apart, LBMS can avoid the cursemore » of dimensionality that creates untenable computational costs associated with DNS. This project has shown that LBMS is capable of reproducing statistics of isotropic turbulent flows while coarsening the spacing between lines significantly. It also investigates and resolves issues that arise when coupling ODT lines, such as flux reconstruction perpendicular to a given ODT line, preservation of conserved quantities when eddies cross a course cell volume and boundary condition application. Robust parallelization is also investigated.« less

GeNN: a code generation framework for accelerated brain simulations

NASA Astrophysics Data System (ADS)

Yavuz, Esin; Turner, James; Nowotny, Thomas

2016-01-01

Large-scale numerical simulations of detailed brain circuit models are important for identifying hypotheses on brain functions and testing their consistency and plausibility. An ongoing challenge for simulating realistic models is, however, computational speed. In this paper, we present the GeNN (GPU-enhanced Neuronal Networks) framework, which aims to facilitate the use of graphics accelerators for computational models of large-scale neuronal networks to address this challenge. GeNN is an open source library that generates code to accelerate the execution of network simulations on NVIDIA GPUs, through a flexible and extensible interface, which does not require in-depth technical knowledge from the users. We present performance benchmarks showing that 200-fold speedup compared to a single core of a CPU can be achieved for a network of one million conductance based Hodgkin-Huxley neurons but that for other models the speedup can differ. GeNN is available for Linux, Mac OS X and Windows platforms. The source code, user manual, tutorials, Wiki, in-depth example projects and all other related information can be found on the project website http://genn-team.github.io/genn/.

GeNN: a code generation framework for accelerated brain simulations.

PubMed

Yavuz, Esin; Turner, James; Nowotny, Thomas

2016-01-07

Large-scale numerical simulations of detailed brain circuit models are important for identifying hypotheses on brain functions and testing their consistency and plausibility. An ongoing challenge for simulating realistic models is, however, computational speed. In this paper, we present the GeNN (GPU-enhanced Neuronal Networks) framework, which aims to facilitate the use of graphics accelerators for computational models of large-scale neuronal networks to address this challenge. GeNN is an open source library that generates code to accelerate the execution of network simulations on NVIDIA GPUs, through a flexible and extensible interface, which does not require in-depth technical knowledge from the users. We present performance benchmarks showing that 200-fold speedup compared to a single core of a CPU can be achieved for a network of one million conductance based Hodgkin-Huxley neurons but that for other models the speedup can differ. GeNN is available for Linux, Mac OS X and Windows platforms. The source code, user manual, tutorials, Wiki, in-depth example projects and all other related information can be found on the project website http://genn-team.github.io/genn/.

GeNN: a code generation framework for accelerated brain simulations

PubMed Central

Yavuz, Esin; Turner, James; Nowotny, Thomas

2016-01-01

Large-scale numerical simulations of detailed brain circuit models are important for identifying hypotheses on brain functions and testing their consistency and plausibility. An ongoing challenge for simulating realistic models is, however, computational speed. In this paper, we present the GeNN (GPU-enhanced Neuronal Networks) framework, which aims to facilitate the use of graphics accelerators for computational models of large-scale neuronal networks to address this challenge. GeNN is an open source library that generates code to accelerate the execution of network simulations on NVIDIA GPUs, through a flexible and extensible interface, which does not require in-depth technical knowledge from the users. We present performance benchmarks showing that 200-fold speedup compared to a single core of a CPU can be achieved for a network of one million conductance based Hodgkin-Huxley neurons but that for other models the speedup can differ. GeNN is available for Linux, Mac OS X and Windows platforms. The source code, user manual, tutorials, Wiki, in-depth example projects and all other related information can be found on the project website http://genn-team.github.io/genn/. PMID:26740369

Development of a Spot-Application Tool for Rapid, High-Resolution Simulation of Wave-Driven Nearshore Hydrodynamics

DTIC Science & Technology

2013-09-30

flow models, such as Delft3D, with our developed Boussinesq -type model. The vision of this project is to develop an operational tool for the...situ measurements or large-scale wave models. This information will be used to drive the offshore wave boundary condition. • Execute the Boussinesq ...model to match with the Boussinesq -type theory would be one which can simulate sheared and stratified currents due to large-scale (non-wave) forcings

Convergence in France facing Big Data era and Exascale challenges for Climate Sciences

NASA Astrophysics Data System (ADS)

Denvil, Sébastien; Dufresne, Jean-Louis; Salas, David; Meurdesoif, Yann; Valcke, Sophie; Caubel, Arnaud; Foujols, Marie-Alice; Servonnat, Jérôme; Sénési, Stéphane; Derouillat, Julien; Voury, Pascal

2014-05-01

The presentation will introduce a french national project : CONVERGENCE that has been funded for four years. This project will tackle big data and computational challenges faced by climate modeling community in HPC context. Model simulations are central to the study of complex mechanisms and feedbacks in the climate system and to provide estimates of future and past climate changes. Recent trends in climate modelling are to add more physical components in the modelled system, increasing the resolution of each individual component and the more systematic use of large suites of simulations to address many scientific questions. Climate simulations may therefore differ in their initial state, parameter values, representation of physical processes, spatial resolution, model complexity, and degree of realism or degree of idealisation. In addition, there is a strong need for evaluating, improving and monitoring the performance of climate models using a large ensemble of diagnostics and better integration of model outputs and observational data. High performance computing is currently reaching the exascale and has the potential to produce this exponential increase of size and numbers of simulations. However, post-processing, analysis, and exploration of the generated data have stalled and there is a strong need for new tools to cope with the growing size and complexity of the underlying simulations and datasets. Exascale simulations require new scalable software tools to generate, manage and mine those simulations ,and data to extract the relevant information and to take the correct decision. The primary purpose of this project is to develop a platform capable of running large ensembles of simulations with a suite of models, to handle the complex and voluminous datasets generated, to facilitate the evaluation and validation of the models and the use of higher resolution models. We propose to gather interdisciplinary skills to design, using a component-based approach, a specific programming environment for scalable scientific simulations and analytics, integrating new and efficient ways of deploying and analysing the applications on High Performance Computing (HPC) system. CONVERGENCE, gathering HPC and informatics expertise that cuts across the individual partners and the broader HPC community, will allow the national climate community to leverage information technology (IT) innovations to address its specific needs. Our methodology consists in developing an ensemble of generic elements needed to run the French climate models with different grids and different resolution, ensuring efficient and reliable execution of these models, managing large volume and number of data and allowing analysis of the results and precise evaluation of the models. These elements include data structure definition and input-output (IO), code coupling and interpolation, as well as runtime and pre/post-processing environments. A common data and metadata structure will allow transferring consistent information between the various elements. All these generic elements will be open source and publicly available. The IPSL-CM and CNRM-CM climate models will make use of these elements that will constitute a national platform for climate modelling. This platform will be used, in its entirety, to optimise and tune the next version of the IPSL-CM model and to develop a global coupled climate model with a regional grid refinement. It will also be used, at least partially, to run ensembles of the CNRM-CM model at relatively high resolution and to run a very-high resolution prototype of this model. The climate models we developed are already involved in many international projects. For instance we participate to the CMIP (Coupled Model Intercomparison Project) project that is very demanding but has a high visibility: its results are widely used and are in particular synthesised in the IPCC (Intergovernmental Panel on Climate Change) assessment reports. The CONVERGENCE project will constitute an invaluable step for the French climate community to prepare and better contribute to the next phase of the CMIP project.

Projected future vegetation changes for the northwest United States and southwest Canada at a fine spatial resolution using a dynamic global vegetation model.

USGS Publications Warehouse

Shafer, Sarah; Bartlein, Patrick J.; Gray, Elizabeth M.; Pelltier, Richard T.

2015-01-01

Future climate change may significantly alter the distributions of many plant taxa. The effects of climate change may be particularly large in mountainous regions where climate can vary significantly with elevation. Understanding potential future vegetation changes in these regions requires methods that can resolve vegetation responses to climate change at fine spatial resolutions. We used LPJ, a dynamic global vegetation model, to assess potential future vegetation changes for a large topographically complex area of the northwest United States and southwest Canada (38.0–58.0°N latitude by 136.6–103.0°W longitude). LPJ is a process-based vegetation model that mechanistically simulates the effect of changing climate and atmospheric CO2 concentrations on vegetation. It was developed and has been mostly applied at spatial resolutions of 10-minutes or coarser. In this study, we used LPJ at a 30-second (~1-km) spatial resolution to simulate potential vegetation changes for 2070–2099. LPJ was run using downscaled future climate simulations from five coupled atmosphere-ocean general circulation models (CCSM3, CGCM3.1(T47), GISS-ER, MIROC3.2(medres), UKMO-HadCM3) produced using the A2 greenhouse gases emissions scenario. Under projected future climate and atmospheric CO2 concentrations, the simulated vegetation changes result in the contraction of alpine, shrub-steppe, and xeric shrub vegetation across the study area and the expansion of woodland and forest vegetation. Large areas of maritime cool forest and cold forest are simulated to persist under projected future conditions. The fine spatial-scale vegetation simulations resolve patterns of vegetation change that are not visible at coarser resolutions and these fine-scale patterns are particularly important for understanding potential future vegetation changes in topographically complex areas.

Projected Future Vegetation Changes for the Northwest United States and Southwest Canada at a Fine Spatial Resolution Using a Dynamic Global Vegetation Model

PubMed Central

Shafer, Sarah L.; Bartlein, Patrick J.; Gray, Elizabeth M.; Pelltier, Richard T.

2015-01-01

Future climate change may significantly alter the distributions of many plant taxa. The effects of climate change may be particularly large in mountainous regions where climate can vary significantly with elevation. Understanding potential future vegetation changes in these regions requires methods that can resolve vegetation responses to climate change at fine spatial resolutions. We used LPJ, a dynamic global vegetation model, to assess potential future vegetation changes for a large topographically complex area of the northwest United States and southwest Canada (38.0–58.0°N latitude by 136.6–103.0°W longitude). LPJ is a process-based vegetation model that mechanistically simulates the effect of changing climate and atmospheric CO2 concentrations on vegetation. It was developed and has been mostly applied at spatial resolutions of 10-minutes or coarser. In this study, we used LPJ at a 30-second (~1-km) spatial resolution to simulate potential vegetation changes for 2070–2099. LPJ was run using downscaled future climate simulations from five coupled atmosphere-ocean general circulation models (CCSM3, CGCM3.1(T47), GISS-ER, MIROC3.2(medres), UKMO-HadCM3) produced using the A2 greenhouse gases emissions scenario. Under projected future climate and atmospheric CO2 concentrations, the simulated vegetation changes result in the contraction of alpine, shrub-steppe, and xeric shrub vegetation across the study area and the expansion of woodland and forest vegetation. Large areas of maritime cool forest and cold forest are simulated to persist under projected future conditions. The fine spatial-scale vegetation simulations resolve patterns of vegetation change that are not visible at coarser resolutions and these fine-scale patterns are particularly important for understanding potential future vegetation changes in topographically complex areas. PMID:26488750

Large-Scale Optimization for Bayesian Inference in Complex Systems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Willcox, Karen; Marzouk, Youssef

2013-11-12

The SAGUARO (Scalable Algorithms for Groundwater Uncertainty Analysis and Robust Optimization) Project focused on the development of scalable numerical algorithms for large-scale Bayesian inversion in complex systems that capitalize on advances in large-scale simulation-based optimization and inversion methods. The project was a collaborative effort among MIT, the University of Texas at Austin, Georgia Institute of Technology, and Sandia National Laboratories. The research was directed in three complementary areas: efficient approximations of the Hessian operator, reductions in complexity of forward simulations via stochastic spectral approximations and model reduction, and employing large-scale optimization concepts to accelerate sampling. The MIT--Sandia component of themore » SAGUARO Project addressed the intractability of conventional sampling methods for large-scale statistical inverse problems by devising reduced-order models that are faithful to the full-order model over a wide range of parameter values; sampling then employs the reduced model rather than the full model, resulting in very large computational savings. Results indicate little effect on the computed posterior distribution. On the other hand, in the Texas--Georgia Tech component of the project, we retain the full-order model, but exploit inverse problem structure (adjoint-based gradients and partial Hessian information of the parameter-to-observation map) to implicitly extract lower dimensional information on the posterior distribution; this greatly speeds up sampling methods, so that fewer sampling points are needed. We can think of these two approaches as ``reduce then sample'' and ``sample then reduce.'' In fact, these two approaches are complementary, and can be used in conjunction with each other. Moreover, they both exploit deterministic inverse problem structure, in the form of adjoint-based gradient and Hessian information of the underlying parameter-to-observation map, to achieve their speedups.« less

Large near-term projected snowpack loss over the western United States

NASA Astrophysics Data System (ADS)

Fyfe, John C.; Derksen, Chris; Mudryk, Lawrence; Flato, Gregory M.; Santer, Benjamin D.; Swart, Neil C.; Molotch, Noah P.; Zhang, Xuebin; Wan, Hui; Arora, Vivek K.; Scinocca, John; Jiao, Yanjun

2017-04-01

Peak runoff in streams and rivers of the western United States is strongly influenced by melting of accumulated mountain snowpack. A significant decline in this resource has a direct connection to streamflow, with substantial economic and societal impacts. Observations and reanalyses indicate that between the 1980s and 2000s, there was a 10-20% loss in the annual maximum amount of water contained in the region's snowpack. Here we show that this loss is consistent with results from a large ensemble of climate simulations forced with natural and anthropogenic changes, but is inconsistent with simulations forced by natural changes alone. A further loss of up to 60% is projected within the next 30 years. Uncertainties in loss estimates depend on the size and the rate of response to continued anthropogenic forcing and the magnitude and phasing of internal decadal variability. The projected losses have serious implications for the hydropower, municipal and agricultural sectors in the region.

The uncertainty of crop yield projections is reduced by improved temperature response functions.

PubMed

Wang, Enli; Martre, Pierre; Zhao, Zhigan; Ewert, Frank; Maiorano, Andrea; Rötter, Reimund P; Kimball, Bruce A; Ottman, Michael J; Wall, Gerard W; White, Jeffrey W; Reynolds, Matthew P; Alderman, Phillip D; Aggarwal, Pramod K; Anothai, Jakarat; Basso, Bruno; Biernath, Christian; Cammarano, Davide; Challinor, Andrew J; De Sanctis, Giacomo; Doltra, Jordi; Fereres, Elias; Garcia-Vila, Margarita; Gayler, Sebastian; Hoogenboom, Gerrit; Hunt, Leslie A; Izaurralde, Roberto C; Jabloun, Mohamed; Jones, Curtis D; Kersebaum, Kurt C; Koehler, Ann-Kristin; Liu, Leilei; Müller, Christoph; Naresh Kumar, Soora; Nendel, Claas; O'Leary, Garry; Olesen, Jørgen E; Palosuo, Taru; Priesack, Eckart; Eyshi Rezaei, Ehsan; Ripoche, Dominique; Ruane, Alex C; Semenov, Mikhail A; Shcherbak, Iurii; Stöckle, Claudio; Stratonovitch, Pierre; Streck, Thilo; Supit, Iwan; Tao, Fulu; Thorburn, Peter; Waha, Katharina; Wallach, Daniel; Wang, Zhimin; Wolf, Joost; Zhu, Yan; Asseng, Senthold

2017-07-17

Increasing the accuracy of crop productivity estimates is a key element in planning adaptation strategies to ensure global food security under climate change. Process-based crop models are effective means to project climate impact on crop yield, but have large uncertainty in yield simulations. Here, we show that variations in the mathematical functions currently used to simulate temperature responses of physiological processes in 29 wheat models account for >50% of uncertainty in simulated grain yields for mean growing season temperatures from 14 °C to 33 °C. We derived a set of new temperature response functions that when substituted in four wheat models reduced the error in grain yield simulations across seven global sites with different temperature regimes by 19% to 50% (42% average). We anticipate the improved temperature responses to be a key step to improve modelling of crops under rising temperature and climate change, leading to higher skill of crop yield projections.

The Uncertainty of Crop Yield Projections Is Reduced by Improved Temperature Response Functions

NASA Technical Reports Server (NTRS)

Wang, Enli; Martre, Pierre; Zhao, Zhigan; Ewert, Frank; Maiorano, Andrea; Rotter, Reimund P.; Kimball, Bruce A.; Ottman, Michael J.; White, Jeffrey W.; Reynolds, Matthew P.;

Numerical simulation of exploding pusher targets

NASA Astrophysics Data System (ADS)

Atzeni, S.; Rosenberg, M. J.; Gatu Johnson, M.; Petrasso, R. D.

2017-10-01

Exploding pusher targets, i.e. gas-filled large aspect-ratio glass or plastic shells, driven by a strong laser-generated shock, are widely used as pulsed sources of neutrons and fast charged particles. Recent experiments on exploding pushers provided evidence for the transition from a purely fluid behavior to a kinetic one. Indeed, fluid models largely overpredict yield and temperature as the Knudsen number Kn (ratio of ion mean-free path to compressed gas radius) is comparable or larger than one. At Kn = 0.3 - 1, fluid codes reasonably estimate integral quantities as yield and neutron-averaged temperatures, but do not reproduce burn radii, burn profiles and DD/DHe3 yield ratio. This motivated a detailed simulation study of intermediate-Kn exploding pushers. We will show how simulation results depend on models for laser-interaction, electron conductivity (flux-limited local vs nonlocal), viscosity (physical vs artificial), and ion mixing. Work partially supported by Sapienza Project C26A15YTMA, Sapienza 2016 (n. 257584), and Eurofusion Project AWP17-ENR-IFE-CEA-01.

Future changes in the climatology of the Great Plains low-level jet derived from fine resolution multi-model simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tang, Ying; Winkler, Julie; Zhong, Shiyuan

The southerly Great Plains low-level jet (GPLLJ) is one of the most significant circulation features of the central U.S. linking large-scale atmospheric circulation with the regional climate. GPLLJs transport heat and moisture, contribute to thunderstorm and severe weather formation, provide a corridor for the springtime migration of birds and insects, enhance wind energy availability, and disperse air pollution. We also assess future changes in GPLLJ frequency using an eight member ensemble of dynamically-downscaled climate simulations for the mid-21st century. Nocturnal GPLLJ frequency is projected to increase in the southern plains in spring and in the central plains in summer, whereasmore » current climatological patterns persist into the future for daytime and cool season GPLLJs. The relationship between future GPLLJ frequency and the extent and strength of anticyclonic airflow over eastern North America varies with season. Most simulations project a westward shift of anticyclonic airflow in summer, but uncertainty is larger for spring with only half of the simulations suggesting a westward expansion. Furthemore, the choice of regional climate model and the driving lateral boundary conditions have a large influence on the projected future changes in GPLLJ frequency and highlight the importance of multi-model ensembles to estimate the uncertainty surrounding the future GPLLJ climatology.« less

Future changes in the climatology of the Great Plains low-level jet derived from fine resolution multi-model simulations

DOE PAGES

Tang, Ying; Winkler, Julie; Zhong, Shiyuan; ...

2017-07-10

The southerly Great Plains low-level jet (GPLLJ) is one of the most significant circulation features of the central U.S. linking large-scale atmospheric circulation with the regional climate. GPLLJs transport heat and moisture, contribute to thunderstorm and severe weather formation, provide a corridor for the springtime migration of birds and insects, enhance wind energy availability, and disperse air pollution. We also assess future changes in GPLLJ frequency using an eight member ensemble of dynamically-downscaled climate simulations for the mid-21st century. Nocturnal GPLLJ frequency is projected to increase in the southern plains in spring and in the central plains in summer, whereasmore » current climatological patterns persist into the future for daytime and cool season GPLLJs. The relationship between future GPLLJ frequency and the extent and strength of anticyclonic airflow over eastern North America varies with season. Most simulations project a westward shift of anticyclonic airflow in summer, but uncertainty is larger for spring with only half of the simulations suggesting a westward expansion. Furthemore, the choice of regional climate model and the driving lateral boundary conditions have a large influence on the projected future changes in GPLLJ frequency and highlight the importance of multi-model ensembles to estimate the uncertainty surrounding the future GPLLJ climatology.« less

Simulating ungulate herbivory across forest landscapes: A browsing extension for LANDIS-II

Treesearch

Nathan R. De Jager; Patrick J. Drohan; Brian M. Miranda; Brian R. Sturtevant; Susan L. Stout; Alejandro A. Royo; Eric J. Gustafson; Mark C. Romanski

2017-01-01

Browsing ungulates alter forest productivity and vegetation succession through selective foraging onspecies that often dominate early succession. However, the long-term and large-scale effects of browsing on forest succession are not possible to project without the use of simulation models. To explore the effects of ungulates on succession in a spatially explicit...

Improving longleaf pine mortality predictions in the Southern Variant of the Forest Vegetation Simulator

Treesearch

R. Justin DeRose; John D. Shaw; Giorgio Vacchiano; James N. Long

2008-01-01

The Southern Variant of the Forest Vegetation Simulator (FVS-SN) is made up of individual submodels that predict tree growth, recruitment and mortality. Forest managers on Ft. Bragg, North Carolina, discovered biologically unrealistic longleaf pine (Pinus palustris) size-density predictions at large diameters when using FVS-SN to project red-cockaded...

Incorporating variability in simulations of seasonally forced phenology using integral projection models

DOE PAGES

Goodsman, Devin W.; Aukema, Brian H.; McDowell, Nate G.; ...

2017-11-26

Phenology models are becoming increasingly important tools to accurately predict how climate change will impact the life histories of organisms. We propose a class of integral projection phenology models derived from stochastic individual-based models of insect development and demography. Our derivation, which is based on the rate summation concept, produces integral projection models that capture the effect of phenotypic rate variability on insect phenology, but which are typically more computationally frugal than equivalent individual-based phenology models. We demonstrate our approach using a temperature-dependent model of the demography of the mountain pine beetle (Dendroctonus ponderosae Hopkins), an insect that kills maturemore » pine trees. This work illustrates how a wide range of stochastic phenology models can be reformulated as integral projection models. Due to their computational efficiency, these integral projection models are suitable for deployment in large-scale simulations, such as studies of altered pest distributions under climate change.« less

Incorporating variability in simulations of seasonally forced phenology using integral projection models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Goodsman, Devin W.; Aukema, Brian H.; McDowell, Nate G.

Phenology models are becoming increasingly important tools to accurately predict how climate change will impact the life histories of organisms. We propose a class of integral projection phenology models derived from stochastic individual-based models of insect development and demography. Our derivation, which is based on the rate summation concept, produces integral projection models that capture the effect of phenotypic rate variability on insect phenology, but which are typically more computationally frugal than equivalent individual-based phenology models. We demonstrate our approach using a temperature-dependent model of the demography of the mountain pine beetle (Dendroctonus ponderosae Hopkins), an insect that kills maturemore » pine trees. This work illustrates how a wide range of stochastic phenology models can be reformulated as integral projection models. Due to their computational efficiency, these integral projection models are suitable for deployment in large-scale simulations, such as studies of altered pest distributions under climate change.« less

Navigating Earthquake Physics with High-Resolution Array Back-Projection

NASA Astrophysics Data System (ADS)

Meng, Lingsen

Understanding earthquake source dynamics is a fundamental goal of geophysics. Progress toward this goal has been slow due to the gap between state-of-art earthquake simulations and the limited source imaging techniques based on conventional low-frequency finite fault inversions. Seismic array processing is an alternative source imaging technique that employs the higher frequency content of the earthquakes and provides finer detail of the source process with few prior assumptions. While the back-projection provides key observations of previous large earthquakes, the standard beamforming back-projection suffers from low resolution and severe artifacts. This thesis introduces the MUSIC technique, a high-resolution array processing method that aims to narrow the gap between the seismic observations and earthquake simulations. The MUSIC is a high-resolution method taking advantage of the higher order signal statistics. The method has not been widely used in seismology yet because of the nonstationary and incoherent nature of the seismic signal. We adapt MUSIC to transient seismic signal by incorporating the Multitaper cross-spectrum estimates. We also adopt a "reference window" strategy that mitigates the "swimming artifact," a systematic drift effect in back projection. The improved MUSIC back projections allow the imaging of recent large earthquakes in finer details which give rise to new perspectives on dynamic simulations. In the 2011 Tohoku-Oki earthquake, we observe frequency-dependent rupture behaviors which relate to the material variation along the dip of the subduction interface. In the 2012 off-Sumatra earthquake, we image the complicated ruptures involving orthogonal fault system and an usual branching direction. This result along with our complementary dynamic simulations probes the pressure-insensitive strength of the deep oceanic lithosphere. In another example, back projection is applied to the 2010 M7 Haiti earthquake recorded at regional distance. The high-frequency subevents are located at the edges of geodetic slip regions, which are correlated to the stopping phases associated with rupture speed reduction when the earthquake arrests.

Climate SPHINX: evaluating the impact of resolution and stochastic physics parameterisations in the EC-Earth global climate model

NASA Astrophysics Data System (ADS)

Davini, Paolo; von Hardenberg, Jost; Corti, Susanna; Christensen, Hannah M.; Juricke, Stephan; Subramanian, Aneesh; Watson, Peter A. G.; Weisheimer, Antje; Palmer, Tim N.

2017-03-01

The Climate SPHINX (Stochastic Physics HIgh resolutioN eXperiments) project is a comprehensive set of ensemble simulations aimed at evaluating the sensitivity of present and future climate to model resolution and stochastic parameterisation. The EC-Earth Earth system model is used to explore the impact of stochastic physics in a large ensemble of 30-year climate integrations at five different atmospheric horizontal resolutions (from 125 up to 16 km). The project includes more than 120 simulations in both a historical scenario (1979-2008) and a climate change projection (2039-2068), together with coupled transient runs (1850-2100). A total of 20.4 million core hours have been used, made available from a single year grant from PRACE (the Partnership for Advanced Computing in Europe), and close to 1.5 PB of output data have been produced on SuperMUC IBM Petascale System at the Leibniz Supercomputing Centre (LRZ) in Garching, Germany. About 140 TB of post-processed data are stored on the CINECA supercomputing centre archives and are freely accessible to the community thanks to an EUDAT data pilot project. This paper presents the technical and scientific set-up of the experiments, including the details on the forcing used for the simulations performed, defining the SPHINX v1.0 protocol. In addition, an overview of preliminary results is given. An improvement in the simulation of Euro-Atlantic atmospheric blocking following resolution increase is observed. It is also shown that including stochastic parameterisation in the low-resolution runs helps to improve some aspects of the tropical climate - specifically the Madden-Julian Oscillation and the tropical rainfall variability. These findings show the importance of representing the impact of small-scale processes on the large-scale climate variability either explicitly (with high-resolution simulations) or stochastically (in low-resolution simulations).

Dynamical downscaling with the fifth-generation Canadian regional climate model (CRCM5) over the CORDEX Arctic domain: effect of large-scale spectral nudging and of empirical correction of sea-surface temperature

NASA Astrophysics Data System (ADS)

Takhsha, Maryam; Nikiéma, Oumarou; Lucas-Picher, Philippe; Laprise, René; Hernández-Díaz, Leticia; Winger, Katja

2017-10-01

As part of the CORDEX project, the fifth-generation Canadian Regional Climate Model (CRCM5) is used over the Arctic for climate simulations driven by reanalyses and by the MPI-ESM-MR coupled global climate model (CGCM) under the RCP8.5 scenario. The CRCM5 shows adequate skills capturing general features of mean sea level pressure (MSLP) for all seasons. Evaluating 2-m temperature (T2m) and precipitation is more problematic, because of inconsistencies between observational reference datasets over the Arctic that suffer of a sparse distribution of weather stations. In our study, we additionally investigated the effect of large-scale spectral nudging (SN) on the hindcast simulation driven by reanalyses. The analysis shows that SN is effective in reducing the spring MSLP bias, but otherwise it has little impact. We have also conducted another experiment in which the CGCM-simulated sea-surface temperature (SST) is empirically corrected and used as lower boundary conditions over the ocean for an atmosphere-only global simulation (AGCM), which in turn provides the atmospheric lateral boundary conditions to drive the CRCM5 simulation. This approach, so-called 3-step approach of dynamical downscaling (CGCM-AGCM-RCM), which had considerably improved the CRCM5 historical simulations over Africa, exhibits reduced impact over the Arctic domain. The most notable positive effect over the Arctic is a reduction of the T2m bias over the North Pacific Ocean and the North Atlantic Ocean in all seasons. Future projections using this method are compared with the results obtained with the traditional 2-step dynamical downscaling (CGCM-RCM) to assess the impact of correcting systematic biases of SST upon future-climate projections. The future projections are mostly similar for the two methods, except for precipitation.

Realism of the Indian Ocean Dipole in CMIP5 models, and the Implication for Climate Projections

NASA Astrophysics Data System (ADS)

Weller, E.; Cai, W.; Cowan, T.

2012-12-01

An assessment of how well climate models simulate the Indian Ocean Dipole (IOD) is undertaken using coupled models that have partaken in the Coupled Model Intercomparison Project Phase 5 (CMIP5). Compared to CMIP3 models, no substantial improvement is evident in the simulation of the IOD pattern and/or amplitude during its peak season in austral spring (September-October-November, or SON). The majority of CMIP5 models generate a larger variance of sea surface temperature (SST) in the Sumatra-Java upwelling region and an IOD amplitude that is far greater than what is observed. Although the relationship between precipitation and the tropical Indian Ocean SST is well simulated, future projections of SON rainfall changes over IOD-influenced regions are intrinsically linked to the IOD-rainfall teleconnection and IOD amplitude in the model present-day climate. The diversity of the simulated IOD amplitudes in CMIP5 (and CMIP3) models which tend to be overly large, results in a wide range of future modelled SON rainfall trends over IOD-influenced regions. Our results highlight the importance of realistically simulating the present-day IOD properties and the caveat that needs to be exercised in interpreting climate projections in the IOD-affected regions.

The CMIP6 Sea-Ice Model Intercomparison Project (SIMIP): Understanding sea ice through climate-model simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Notz, Dirk; Jahn, Alexandra; Holland, Marika

A better understanding of the role of sea ice for the changing climate of our planet is the central aim of the diagnostic Coupled Model Intercomparison Project 6 (CMIP6)-endorsed Sea-Ice Model Intercomparison Project (SIMIP). To reach this aim, SIMIP requests sea-ice-related variables from climate-model simulations that allow for a better understanding and, ultimately, improvement of biases and errors in sea-ice simulations with large-scale climate models. This then allows us to better understand to what degree CMIP6 model simulations relate to reality, thus improving our confidence in answering sea-ice-related questions based on these simulations. Furthermore, the SIMIP protocol provides a standardmore » for sea-ice model output that will streamline and hence simplify the analysis of the simulated sea-ice evolution in research projects independent of CMIP. To reach its aims, SIMIP provides a structured list of model output that allows for an examination of the three main budgets that govern the evolution of sea ice, namely the heat budget, the momentum budget, and the mass budget. Furthermore, we explain the aims of SIMIP in more detail and outline how its design allows us to answer some of the most pressing questions that sea ice still poses to the international climate-research community.« less

The CMIP6 Sea-Ice Model Intercomparison Project (SIMIP): Understanding sea ice through climate-model simulations

DOE PAGES

Notz, Dirk; Jahn, Alexandra; Holland, Marika; ...

2016-09-23

A better understanding of the role of sea ice for the changing climate of our planet is the central aim of the diagnostic Coupled Model Intercomparison Project 6 (CMIP6)-endorsed Sea-Ice Model Intercomparison Project (SIMIP). To reach this aim, SIMIP requests sea-ice-related variables from climate-model simulations that allow for a better understanding and, ultimately, improvement of biases and errors in sea-ice simulations with large-scale climate models. This then allows us to better understand to what degree CMIP6 model simulations relate to reality, thus improving our confidence in answering sea-ice-related questions based on these simulations. Furthermore, the SIMIP protocol provides a standardmore » for sea-ice model output that will streamline and hence simplify the analysis of the simulated sea-ice evolution in research projects independent of CMIP. To reach its aims, SIMIP provides a structured list of model output that allows for an examination of the three main budgets that govern the evolution of sea ice, namely the heat budget, the momentum budget, and the mass budget. Furthermore, we explain the aims of SIMIP in more detail and outline how its design allows us to answer some of the most pressing questions that sea ice still poses to the international climate-research community.« less

Damping characterization in large structures

NASA Technical Reports Server (NTRS)

Eke, Fidelis O.; Eke, Estelle M.

1991-01-01

This research project has as its main goal the development of methods for selecting the damping characteristics of components of a large structure or multibody system, in such a way as to produce some desired system damping characteristics. The main need for such an analytical device is in the simulation of the dynamics of multibody systems consisting, at least partially, of flexible components. The reason for this need is that all existing simulation codes for multibody systems require component-by-component characterization of complex systems, whereas requirements (including damping) often appear at the overall system level. The main goal was met in large part by the development of a method that will in fact synthesize component damping matrices from a given system damping matrix. The restrictions to the method are that the desired system damping matrix must be diagonal (which is almost always the case) and that interbody connections must be by simple hinges. In addition to the technical outcome, this project contributed positively to the educational and research infrastructure of Tuskegee University - a Historically Black Institution.

eScience for molecular-scale simulations and the eMinerals project.

PubMed

Salje, E K H; Artacho, E; Austen, K F; Bruin, R P; Calleja, M; Chappell, H F; Chiang, G-T; Dove, M T; Frame, I; Goodwin, A L; Kleese van Dam, K; Marmier, A; Parker, S C; Pruneda, J M; Todorov, I T; Trachenko, K; Tyer, R P; Walker, A M; White, T O H

2009-03-13

We review the work carried out within the eMinerals project to develop eScience solutions that facilitate a new generation of molecular-scale simulation work. Technological developments include integration of compute and data systems, developing of collaborative frameworks and new researcher-friendly tools for grid job submission, XML data representation, information delivery, metadata harvesting and metadata management. A number of diverse science applications will illustrate how these tools are being used for large parameter-sweep studies, an emerging type of study for which the integration of computing, data and collaboration is essential.

Freak waves in random oceanic sea states.

PubMed

Onorato, M; Osborne, A R; Serio, M; Bertone, S

2001-06-18

Freak waves are very large, rare events in a random ocean wave train. Here we study their generation in a random sea state characterized by the Joint North Sea Wave Project spectrum. We assume, to cubic order in nonlinearity, that the wave dynamics are governed by the nonlinear Schrödinger (NLS) equation. We show from extensive numerical simulations of the NLS equation how freak waves in a random sea state are more likely to occur for large values of the Phillips parameter alpha and the enhancement coefficient gamma. Comparison with linear simulations is also reported.

The Roland Maze Project school-based extensive air shower network

NASA Astrophysics Data System (ADS)

Feder, J.; Jȩdrzejczak, K.; Karczmarczyk, J.; Lewandowski, R.; Swarzyński, J.; Szabelska, B.; Szabelski, J.; Wibig, T.

2006-01-01

We plan to construct the large area network of extensive air shower detectors placed on the roofs of high school buildings in the city of Łódź. Detection points will be connected by INTERNET to the central server and their work will be synchronized by GPS. The main scientific goal of the project are studies of ultra high energy cosmic rays. Using existing town infrastructure (INTERNET, power supply, etc.) will significantly reduce the cost of the experiment. Engaging high school students in the research program should significantly increase their knowledge of science and modern technologies, and can be a very efficient way of science popularisation. We performed simulations of the projected network capabilities of registering Extensive Air Showers and reconstructing energies of primary particles. Results of the simulations and the current status of project realisation will be presented.

Using historical simulations of vegetation to assess departure of current vegetation conditions across large landscapes[Chapter 11

Treesearch

Lisa Holsinger; Robert E. Keane; Brian Steele; Matthew C. Reeves; Sarah Pratt

2006-01-01

The Landscape Fire and Resource Management Planning Tools Prototype Project, or LANDFIRE Prototype Project, was conceived, in part, to identify areas across the nation where existing landscape conditions are markedly different from historical conditions (Keane and Rollins, Ch. 3). This objective arose from the recognition that over 100 years of land use and wildland...

North Atlantic winter eddy-driven jet and atmospheric blocking variability in the Community Earth System Model version 1 Large Ensemble simulations

NASA Astrophysics Data System (ADS)

Kwon, Young-Oh; Camacho, Alicia; Martinez, Carlos; Seo, Hyodae

2018-01-01

The atmospheric jet and blocking distributions, especially in the North Atlantic sector, have been challenging features for a climate model to realistically reproduce. This study examines climatological distributions of winter (December-February) daily jet latitude and blocking in the North Atlantic from the 40-member Community Earth System Model version 1 Large Ensemble (CESM1LE) simulations. This analysis aims at examining whether a broad range of internal climate variability encompassed by a large ensemble of simulations results in an improved representation of the jet latitude distributions and blocking days in CESM1LE. In the historical runs (1951-2005), the daily zonal wind at 850 hPa exhibits three distinct preferred latitudes for the eddy-driven jet position as seen in the reanalysis datasets, which represents a significant improvement from the previous version of the same model. However, the meridional separations between the three jet latitudes are much smaller than those in the reanalyses. In particular, the jet rarely migrates to the observed southernmost position around 37°N. This leads to the bias in blocking frequency that is too low over Greenland and too high over the Azores. These features are shown to be remarkably stable across the 40 ensemble members with negligible member-to-member spread. This result implies the range of internal variability of winter jet latitude and blocking frequency within the 55-year segment from each ensemble member is comparable to that represented by the full large ensemble. Comparison with 2046-2100 from the RCP8.5 future projection runs suggests that the daily jet position is projected to maintain the same three preferred latitudes, with a slightly higher frequency of occurrence over the central latitude around 50°N, instead of shifting poleward in the future as documented in some previous studies. In addition, the daily jet speed is projected not to change significantly between 1951-2005 and 2046-2100. On the other hand, the climatological mean jet is projected to become slightly more elongated and stronger on its southern flank, and the blocking frequency over the Azores is projected to decrease.

Subgrid-scale models for large-eddy simulation of rotating turbulent channel flows

NASA Astrophysics Data System (ADS)

Silvis, Maurits H.; Bae, Hyunji Jane; Trias, F. Xavier; Abkar, Mahdi; Moin, Parviz; Verstappen, Roel

2017-11-01

We aim to design subgrid-scale models for large-eddy simulation of rotating turbulent flows. Rotating turbulent flows form a challenging test case for large-eddy simulation due to the presence of the Coriolis force. The Coriolis force conserves the total kinetic energy while transporting it from small to large scales of motion, leading to the formation of large-scale anisotropic flow structures. The Coriolis force may also cause partial flow laminarization and the occurrence of turbulent bursts. Many subgrid-scale models for large-eddy simulation are, however, primarily designed to parametrize the dissipative nature of turbulent flows, ignoring the specific characteristics of transport processes. We, therefore, propose a new subgrid-scale model that, in addition to the usual dissipative eddy viscosity term, contains a nondissipative nonlinear model term designed to capture transport processes, such as those due to rotation. We show that the addition of this nonlinear model term leads to improved predictions of the energy spectra of rotating homogeneous isotropic turbulence as well as of the Reynolds stress anisotropy in spanwise-rotating plane-channel flows. This work is financed by the Netherlands Organisation for Scientific Research (NWO) under Project Number 613.001.212.

Dakota Uncertainty Quantification Methods Applied to the CFD code Nek5000

DOE Office of Scientific and Technical Information (OSTI.GOV)

Delchini, Marc-Olivier; Popov, Emilian L.; Pointer, William David

This report presents the state of advancement of a Nuclear Energy Advanced Modeling and Simulation (NEAMS) project to characterize the uncertainty of the computational fluid dynamics (CFD) code Nek5000 using the Dakota package for flows encountered in the nuclear engineering industry. Nek5000 is a high-order spectral element CFD code developed at Argonne National Laboratory for high-resolution spectral-filtered large eddy simulations (LESs) and unsteady Reynolds-averaged Navier-Stokes (URANS) simulations.

Effects of different regional climate model resolution and forcing scales on projected hydrologic changes

NASA Astrophysics Data System (ADS)

Mendoza, Pablo A.; Mizukami, Naoki; Ikeda, Kyoko; Clark, Martyn P.; Gutmann, Ethan D.; Arnold, Jeffrey R.; Brekke, Levi D.; Rajagopalan, Balaji

2016-10-01

We examine the effects of regional climate model (RCM) horizontal resolution and forcing scaling (i.e., spatial aggregation of meteorological datasets) on the portrayal of climate change impacts. Specifically, we assess how the above decisions affect: (i) historical simulation of signature measures of hydrologic behavior, and (ii) projected changes in terms of annual water balance and hydrologic signature measures. To this end, we conduct our study in three catchments located in the headwaters of the Colorado River basin. Meteorological forcings for current and a future climate projection are obtained at three spatial resolutions (4-, 12- and 36-km) from dynamical downscaling with the Weather Research and Forecasting (WRF) regional climate model, and hydrologic changes are computed using four different hydrologic model structures. These projected changes are compared to those obtained from running hydrologic simulations with current and future 4-km WRF climate outputs re-scaled to 12- and 36-km. The results show that the horizontal resolution of WRF simulations heavily affects basin-averaged precipitation amounts, propagating into large differences in simulated signature measures across model structures. The implications of re-scaled forcing datasets on historical performance were primarily observed on simulated runoff seasonality. We also found that the effects of WRF grid resolution on projected changes in mean annual runoff and evapotranspiration may be larger than the effects of hydrologic model choice, which surpasses the effects from re-scaled forcings. Scaling effects on projected variations in hydrologic signature measures were found to be generally smaller than those coming from WRF resolution; however, forcing aggregation in many cases reversed the direction of projected changes in hydrologic behavior.

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

NASA Astrophysics Data System (ADS)

Haarsma, Reindert J.; Roberts, Malcolm J.; Vidale, Pier Luigi; Senior, Catherine A.; Bellucci, Alessio; Bao, Qing; Chang, Ping; Corti, Susanna; Fučkar, Neven S.; Guemas, Virginie; von Hardenberg, Jost; Hazeleger, Wilco; Kodama, Chihiro; Koenigk, Torben; Leung, L. Ruby; Lu, Jian; Luo, Jing-Jia; Mao, Jiafu; Mizielinski, Matthew S.; Mizuta, Ryo; Nobre, Paulo; Satoh, Masaki; Scoccimarro, Enrico; Semmler, Tido; Small, Justin; von Storch, Jin-Song

2016-11-01

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950-2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. HighResMIP thereby focuses on one of the CMIP6 broad questions, "what are the origins and consequences of systematic model biases?", but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

Climate change impact on streamflow in large-scale river basins: projections and their uncertainties sourced from GCMs and RCP scenarios

NASA Astrophysics Data System (ADS)

Nasonova, Olga N.; Gusev, Yeugeniy M.; Kovalev, Evgeny E.; Ayzel, Georgy V.

2018-06-01

Climate change impact on river runoff was investigated within the framework of the second phase of the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP2) using a physically-based land surface model Soil Water - Atmosphere - Plants (SWAP) (developed in the Institute of Water Problems of the Russian Academy of Sciences) and meteorological projections (for 2006-2099) simulated by five General Circulation Models (GCMs) (including GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, and NorESM1-M) for each of four Representative Concentration Pathway (RCP) scenarios (RCP2.6, RCP4.5, RCP6.0, and RCP8.5). Eleven large-scale river basins were used in this study. First of all, SWAP was calibrated and validated against monthly values of measured river runoff with making use of forcing data from the WATCH data set and all GCMs' projections were bias-corrected to the WATCH. Then, for each basin, 20 projections of possible changes in river runoff during the 21st century were simulated by SWAP. Analysis of the obtained hydrological projections allowed us to estimate their uncertainties resulted from application of different GCMs and RCP scenarios. On the average, the contribution of different GCMs to the uncertainty of the projected river runoff is nearly twice larger than the contribution of RCP scenarios. At the same time the contribution of GCMs slightly decreases with time.

Python-based geometry preparation and simulation visualization toolkits for STEPS

PubMed Central

Chen, Weiliang; De Schutter, Erik

2014-01-01

STEPS is a stochastic reaction-diffusion simulation engine that implements a spatial extension of Gillespie's Stochastic Simulation Algorithm (SSA) in complex tetrahedral geometries. An extensive Python-based interface is provided to STEPS so that it can interact with the large number of scientific packages in Python. However, a gap existed between the interfaces of these packages and the STEPS user interface, where supporting toolkits could reduce the amount of scripting required for research projects. This paper introduces two new supporting toolkits that support geometry preparation and visualization for STEPS simulations. PMID:24782754

The CESM Large Ensemble Project: Inspiring New Ideas and Understanding

NASA Astrophysics Data System (ADS)

Kay, J. E.; Deser, C.

2016-12-01

While internal climate variability is known to affect climate projections, its influence is often underappreciated and confused with model error. Why? In general, modeling centers contribute a small number of realizations to international climate model assessments [e.g., phase 5 of the Coupled Model Intercomparison Project (CMIP5)]. As a result, model error and internal climate variability are difficult, and at times impossible, to disentangle. In response, the Community Earth System Model (CESM) community designed the CESM Large Ensemble (CESM-LE) with the explicit goal of enabling assessment of climate change in the presence of internal climate variability. All CESM-LE simulations use a single CMIP5 model (CESM with the Community Atmosphere Model, version 5). The core simulations replay the twenty to twenty-first century (1920-2100) 40+ times under historical and representative concentration pathway 8.5 external forcing with small initial condition differences. Two companion 2000+-yr-long preindustrial control simulations (fully coupled, prognostic atmosphere and land only) allow assessment of internal climate variability in the absence of climate change. Comprehensive outputs, including many daily fields, are available as single-variable time series on the Earth System Grid for anyone to use. Examples of scientists and stakeholders that are using the CESM-LE outputs to help interpret the observational record, to understand projection spread and to plan for a range of possible futures influenced by both internal climate variability and forced climate change will be highlighted the presentation.

Ground motion simulations in Marmara (Turkey) region from 3D finite difference method

NASA Astrophysics Data System (ADS)

Aochi, Hideo; Ulrich, Thomas; Douglas, John

2016-04-01

In the framework of the European project MARSite (2012-2016), one of the main contributions from our research team was to provide ground-motion simulations for the Marmara region from various earthquake source scenarios. We adopted a 3D finite difference code, taking into account the 3D structure around the Sea of Marmara (including the bathymetry) and the sea layer. We simulated two moderate earthquakes (about Mw4.5) and found that the 3D structure improves significantly the waveforms compared to the 1D layer model. Simulations were carried out for different earthquakes (moderate point sources and large finite sources) in order to provide shake maps (Aochi and Ulrich, BSSA, 2015), to study the variability of ground-motion parameters (Douglas & Aochi, BSSA, 2016) as well as to provide synthetic seismograms for the blind inversion tests (Diao et al., GJI, 2016). The results are also planned to be integrated in broadband ground-motion simulations, tsunamis generation and simulations of triggered landslides (in progress by different partners). The simulations are freely shared among the partners via the internet and the visualization of the results is diffused on the project's homepage. All these simulations should be seen as a reference for this region, as they are based on the latest knowledge that obtained during the MARSite project, although their refinement and validation of the model parameters and the simulations are a continuing research task relying on continuing observations. The numerical code used, the models and the simulations are available on demand.

Neuroscience thinks big (and collaboratively).

PubMed

Kandel, Eric R; Markram, Henry; Matthews, Paul M; Yuste, Rafael; Koch, Christof

2013-09-01

Despite cash-strapped times for research, several ambitious collaborative neuroscience projects have attracted large amounts of funding and media attention. In Europe, the Human Brain Project aims to develop a large-scale computer simulation of the brain, whereas in the United States, the Brain Activity Map is working towards establishing a functional connectome of the entire brain, and the Allen Institute for Brain Science has embarked upon a 10-year project to understand the mouse visual cortex (the MindScope project). US President Barack Obama's announcement of the BRAIN Initiative (Brain Research through Advancing Innovative Neurotechnologies Initiative) in April 2013 highlights the political commitment to neuroscience and is expected to further foster interdisciplinary collaborations, accelerate the development of new technologies and thus fuel much needed medical advances. In this Viewpoint article, five prominent neuroscientists explain the aims of the projects and how they are addressing some of the questions (and criticisms) that have arisen.

Final Report on DOE Project entitled Dynamic Optimized Advanced Scheduling of Bandwidth Demands for Large-Scale Science Applications

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ramamurthy, Byravamurthy

2014-05-05

In this project, developed scheduling frameworks for dynamic bandwidth demands for large-scale science applications. In particular, we developed scheduling algorithms for dynamic bandwidth demands in this project. Apart from theoretical approaches such as Integer Linear Programming, Tabu Search and Genetic Algorithm heuristics, we have utilized practical data from ESnet OSCARS project (from our DOE lab partners) to conduct realistic simulations of our approaches. We have disseminated our work through conference paper presentations and journal papers and a book chapter. In this project we addressed the problem of scheduling of lightpaths over optical wavelength division multiplexed (WDM) networks. We published severalmore » conference papers and journal papers on this topic. We also addressed the problems of joint allocation of computing, storage and networking resources in Grid/Cloud networks and proposed energy-efficient mechanisms for operatin optical WDM networks.« less

A Computational Study of the Flow Physics of Acoustic Liners

NASA Technical Reports Server (NTRS)

Tam, Christopher

2006-01-01

The present investigation is a continuation of a previous joint project between the Florida State University and the NASA Langley Research Center Liner Physics Team. In the previous project, a study of acoustic liners, in two dimensions, inside a normal incidence impedance tube was carried out. The study consisted of two parts. The NASA team was responsible for the experimental part of the project. This involved performing measurements in an impedance tube with a large aspect ratio slit resonator. The FSU team was responsible for the computation part of the project. This involved performing direct numerical simulation (DNS) of the NASA experiment in two dimensions using CAA methodology. It was agreed that upon completion of numerical simulation, the computed values of the liner impedance were to be sent to NASA for validation with experimental results. On following this procedure good agreements were found between numerical results and experimental measurements over a wide range of frequencies and sound-pressure-level. Broadband incident sound waves were also simulated numerically and measured experimentally. Overall, good agreements were also found.

Search for hidden high-Z materials inside containers with the Muon Portal Project

NASA Astrophysics Data System (ADS)

La Rocca, P.; Antonuccio, V.; Bandieramonte, M.; Becciani, U.; Belluomo, F.; Belluso, M.; Billotta, S.; Blancato, A. A.; Bonanno, D.; Bonanno, G.; Costa, A.; Fallica, G.; Garozzo, S.; Indelicato, V.; Leonora, E.; Longhitano, F.; Longo, S.; Lo Presti, D.; Massimino, P.; Petta, C.; Pistagna, C.; Pugliatti, C.; Puglisi, M.; Randazzo, N.; Riggi, F.; Riggi, S.; Romeo, G.; Russo, G. V.; Santagati, G.; Valvo, G.; Vitello, F.; Zaia, A.; Zappalà, G.

2014-01-01

The Muon Portal is a recently born project that plans to build a large area muon detector for a noninvasive inspection of shipping containers in the ports, searching for the presence of potential fissile (U, Pu) threats. The technique employed by the project is the well-known muon tomography, based on cosmic muon scattering from high-Z materials. The design and operational parameters of the muon portal under construction will be described in this paper, together with preliminary simulation and test results.

Realtime monitoring of bridge scour using remote monitoring technology

DOT National Transportation Integrated Search

2011-02-01

The research performed in this project focuses on the application of instruments including accelerometers : and tiltmeters to monitor bridge scour. First, two large scale laboratory experiments were performed. One : experiment is the simulation of a ...

A Highly Resolved Large-Eddy Simulation of a Wind Turbine using an Actuator Line Model with Optimal Body Force Projection

DOE PAGES

Martinez-Tossas, Luis A.; Churchfield, Matthew J.; Meneveau, Charles

2016-10-03

When representing the blade aerodynamics with rotating actuator lines, the computed forces have to be projected back to the CFD flow field as a volumetric body force. That has been done in the past with a geometrically simple uniform three-dimensional Gaussian at each point along the blade. Here, we argue that the body force can be shaped in a way that better predicts the blade local flow field, the blade load distribution, and the formation of the tip/root vortices. In previous work, we have determined the optimal scales of circular and elliptical Gaussian kernels that best reproduce the local flowmore » field in two-dimensions. Lastly, in this work we extend the analysis and applications by considering the full three-dimensional blade to test our hypothesis in a highly resolved Large Eddy Simulation.« less

A Highly Resolved Large-Eddy Simulation of a Wind Turbine using an Actuator Line Model with Optimal Body Force Projection

DOE Office of Scientific and Technical Information (OSTI.GOV)

Martinez-Tossas, Luis A.; Churchfield, Matthew J.; Meneveau, Charles

When representing the blade aerodynamics with rotating actuator lines, the computed forces have to be projected back to the CFD flow field as a volumetric body force. That has been done in the past with a geometrically simple uniform three-dimensional Gaussian at each point along the blade. Here, we argue that the body force can be shaped in a way that better predicts the blade local flow field, the blade load distribution, and the formation of the tip/root vortices. In previous work, we have determined the optimal scales of circular and elliptical Gaussian kernels that best reproduce the local flowmore » field in two-dimensions. Lastly, in this work we extend the analysis and applications by considering the full three-dimensional blade to test our hypothesis in a highly resolved Large Eddy Simulation.« less

World, We Have Problems: Simulation for Large Complex, Risky Projects, and Events

NASA Technical Reports Server (NTRS)

Elfrey, Priscilla

2010-01-01

Prior to a spacewalk during the NASA STS/129 mission in November 2009, Columbia Broadcasting System (CBS) correspondent William Harwood reported astronauts, "were awakened again", as they had been the day previously. Fearing something not properly connected was causing a leak, the crew, both on the ground and in space, stopped and checked everything. The alarm proved false. The crew did complete its work ahead of schedule, but the incident reminds us that correctly connecting hundreds and thousands of entities, subsystems and systems, finding leaks, loosening stuck valves, and adding replacements to very large complex systems over time does not occur magically. Everywhere major projects present similar pressures. Lives are at - risk. Responsibility is heavy. Large natural and human-created disasters introduce parallel difficulties as people work across boundaries their countries, disciplines, languages, and cultures with known immediate dangers as well as the unexpected. NASA has long accepted that when humans have to go where humans cannot go that simulation is the sole solution. The Agency uses simulation to achieve consensus, reduce ambiguity and uncertainty, understand problems, make decisions, support design, do planning and troubleshooting, as well as for operations, training, testing, and evaluation. Simulation is at the heart of all such complex systems, products, projects, programs, and events. Difficult, hazardous short and, especially, long-term activities have a persistent need for simulation from the first insight into a possibly workable idea or answer until the final report perhaps beyond our lifetime is put in the archive. With simulation we create a common mental model, try-out breakdowns of machinery or teamwork, and find opportunity for improvement. Lifecycle simulation proves to be increasingly important as risks and consequences intensify. Across the world, disasters are increasing. We anticipate more of them, as the results of global warming prove more and more ominous-glaciers melting in Bolivia, floods in Saudi Arabia, the Maldives sinking and salt rising along the Nile. Fear grows about potential asteroid crashes and nightly television images raise awareness of victims of floods, hurricanes, cyclones and typhoons, fire, tornado, tsunami, bombings, landslides, and cross-boundary criminality. The Red Cross says that disasters impact 250 million people each year. That means that 700,000 people are having a very bad day today. Modeling and simulation is and must be part of the solution.

Training in Methods in Computational Neuroscience

DTIC Science & Technology

1992-08-29

in Tritonia. Roger Traub Models with realistic neurons , with an emphasis on large-scale modeling of epileptic phenomena in hippocampus. Rodolpho...Cell Model Plan: 1) Convert some of my simulations from NEURON to GENESIS (and thus learn GENESIS). 2) Develop a realistic inhibtory model . 3) Further...General Hospital, MA Course Project: Membrane Properties of a Neostriatal Neuron and Dopamine Modulation The purpose of my project was to model the

Drought in the Horn of Africa: attribution of a damaging and repeating extreme event

NASA Astrophysics Data System (ADS)

Marthews, Toby; Otto, Friederike; Mitchell, Daniel; Dadson, Simon; Jones, Richard

2015-04-01

We have applied detection and attribution techniques to the severe drought that hit the Horn of Africa in 2014. The short rains failed in late 2013 in Kenya, South Sudan, Somalia and southern Ethiopia, leading to a very dry growing season January to March 2014, and subsequently to the current drought in many agricultural areas of the sub-region. We have made use of the weather@home project, which uses publicly-volunteered distributed computing to provide a large ensemble of simulations sufficient to sample regional climate uncertainty. Based on this, we have estimated the occurrence rates of the kinds of the rare and extreme events implicated in this large-scale drought. From land surface model runs based on these ensemble simulations, we have estimated the impacts of climate anomalies during this period and therefore we can reliably identify some factors of the ongoing drought as attributable to human-induced climate change. The UNFCCC's Adaptation Fund is attempting to support projects that bring about an adaptation to "the adverse effects of climate change", but in order to formulate such projects we need a much clearer way to assess how much climate change is human-induced and how much is a consequence of climate anomalies and large-scale teleconnections, which can only be provided by robust attribution techniques.

Large-scale tropospheric transport in the Chemistry-Climate Model Initiative (CCMI) simulations

NASA Astrophysics Data System (ADS)

Orbe, Clara; Yang, Huang; Waugh, Darryn W.; Zeng, Guang; Morgenstern, Olaf; Kinnison, Douglas E.; Lamarque, Jean-Francois; Tilmes, Simone; Plummer, David A.; Scinocca, John F.; Josse, Beatrice; Marecal, Virginie; Jöckel, Patrick; Oman, Luke D.; Strahan, Susan E.; Deushi, Makoto; Tanaka, Taichu Y.; Yoshida, Kohei; Akiyoshi, Hideharu; Yamashita, Yousuke; Stenke, Andreas; Revell, Laura; Sukhodolov, Timofei; Rozanov, Eugene; Pitari, Giovanni; Visioni, Daniele; Stone, Kane A.; Schofield, Robyn; Banerjee, Antara

2018-05-01

Understanding and modeling the large-scale transport of trace gases and aerosols is important for interpreting past (and projecting future) changes in atmospheric composition. Here we show that there are large differences in the global-scale atmospheric transport properties among the models participating in the IGAC SPARC Chemistry-Climate Model Initiative (CCMI). Specifically, we find up to 40 % differences in the transport timescales connecting the Northern Hemisphere (NH) midlatitude surface to the Arctic and to Southern Hemisphere high latitudes, where the mean age ranges between 1.7 and 2.6 years. We show that these differences are related to large differences in vertical transport among the simulations, in particular to differences in parameterized convection over the oceans. While stronger convection over NH midlatitudes is associated with slower transport to the Arctic, stronger convection in the tropics and subtropics is associated with faster interhemispheric transport. We also show that the differences among simulations constrained with fields derived from the same reanalysis products are as large as (and in some cases larger than) the differences among free-running simulations, most likely due to larger differences in parameterized convection. Our results indicate that care must be taken when using simulations constrained with analyzed winds to interpret the influence of meteorology on tropospheric composition.

The Space-Time Conservative Schemes for Large-Scale, Time-Accurate Flow Simulations with Tetrahedral Meshes

NASA Technical Reports Server (NTRS)

Venkatachari, Balaji Shankar; Streett, Craig L.; Chang, Chau-Lyan; Friedlander, David J.; Wang, Xiao-Yen; Chang, Sin-Chung

2016-01-01

Despite decades of development of unstructured mesh methods, high-fidelity time-accurate simulations are still predominantly carried out on structured, or unstructured hexahedral meshes by using high-order finite-difference, weighted essentially non-oscillatory (WENO), or hybrid schemes formed by their combinations. In this work, the space-time conservation element solution element (CESE) method is used to simulate several flow problems including supersonic jet/shock interaction and its impact on launch vehicle acoustics, and direct numerical simulations of turbulent flows using tetrahedral meshes. This paper provides a status report for the continuing development of the space-time conservation element solution element (CESE) numerical and software framework under the Revolutionary Computational Aerosciences (RCA) project. Solution accuracy and large-scale parallel performance of the numerical framework is assessed with the goal of providing a viable paradigm for future high-fidelity flow physics simulations.

Using microstructure observations to quantify fracture properties and improve reservoir simulations. Final report, September 1998

DOE Office of Scientific and Technical Information (OSTI.GOV)

Laubach, S.E.; Marrett, R.; Rossen, W.

The research for this project provides new technology to understand and successfully characterize, predict, and simulate reservoir-scale fractures. Such fractures have worldwide importance because of their influence on successful extraction of resources. The scope of this project includes creation and testing of new methods to measure, interpret, and simulate reservoir fractures that overcome the challenge of inadequate sampling. The key to these methods is the use of microstructures as guides to the attributes of the large fractures that control reservoir behavior. One accomplishment of the project research is a demonstration that these microstructures can be reliably and inexpensively sampled. Specificmore » goals of this project were to: create and test new methods of measuring attributes of reservoir-scale fractures, particularly as fluid conduits, and test the methods on samples from reservoirs; extrapolate structural attributes to the reservoir scale through rigorous mathematical techniques and help build accurate and useful 3-D models of the interwell region; and design new ways to incorporate geological and geophysical information into reservoir simulation and verify the accuracy by comparison with production data. New analytical methods developed in the project are leading to a more realistic characterization of fractured reservoir rocks. Testing diagnostic and predictive approaches was an integral part of the research, and several tests were successfully completed.« less

Simulation of the effects of ground-water withdrawals and recharge on ground-water flow in Cape Cod, Martha's Vineyard, and Nantucket Island basins, Massachusetts

USGS Publications Warehouse

Masterson, John P.; Barlow, Paul M.

1994-01-01

The effects of changing patterns of ground-water pumping and aquifer recharge on the surface-water and ground-water hydrologic systems were determined for the Cape Cod, Martha's Vineyard, and Nantucket Island Basins. Three-dimensional, transient, ground-water-flow modelS that simulate both freshwater and saltwater flow were developed for the f1ow cells of Cape Cod which currently have large-capacity public-supply wells. Only the freshwater-flow system was simulated for the Cape Cod flow cells where public-water supply demands are satisfied by small-capacity domestic wells. Two- dimensional, finite-difference, change models were developed for Martha's Vineyard and Nantucket Island to determine the projected drawdowns in response to projected in-season pumping rates for 180 days of no aquifer recharge. Results of the simulations indicate very little change in the position of the freshwater-saltwater interface from predevelopment flow conditions to projected ground-water pumping and recharge rates for Cape Cod in the year 2020. Results of change model simulations for Martha's Vineyard and Nantucket Island indicate that the greatest impact in response to projected in-season ground-water pumping occurs at the pumping centers and the magnitude of the drawdowns are minimal with respect to the total thickness of the aquifers.

More robust regional precipitation projection from selected CMIP5 models based on multiple-dimensional metrics

NASA Astrophysics Data System (ADS)

Qian, Y.; Wang, L.; Leung, L. R.; Lin, G.; Lu, J.; Gao, Y.; Zhang, Y.

2017-12-01

Projecting precipitation changes is challenging because of incomplete understanding of the climate system and biases and uncertainty in climate models. In East Asia where summer precipitation is dominantly influenced by the monsoon circulation and the global models from Coupled Model Intercomparison Project Phase 5 (CMIP5), however, give various projection of precipitation change for 21th century. It is critical for community to know which models' projection are more reliable in response to natural and anthropogenic forcings. In this study we defined multiple-dimensional metrics, measuring the model performance in simulating the present-day of large-scale circulation, regional precipitation and relationship between them. The large-scale circulation features examined in this study include the lower tropospheric southwesterly winds, the western North Pacific subtropical high, the South China Sea Subtropical High, and the East Asian westerly jet in the upper troposphere. Each of these circulation features transport moisture to East Asia, enhancing the moist static energy and strengthening the Meiyu moisture front that is the primary mechanism for precipitation generation in eastern China. Based on these metrics, 30 models in CMIP5 ensemble are classified into three groups. Models in the top performing group projected regional precipitation patterns that are more similar to each other than the bottom or middle performing group and consistently projected statistically significant increasing trends in two of the large-scale circulation indices and precipitation. In contrast, models in the bottom or middle performing group projected small drying or no trends in precipitation. We also find the models that only reasonably reproduce the observed precipitation climatology does not guarantee more reliable projection of future precipitation because good simulation skill could be achieved through compensating errors from multiple sources. Herein the potential for more robust projections of precipitation changes at regional scale is demonstrated through the use of discriminating metric to subsample the multi-model ensemble. The results from this study provides insights for how to select models from CMIP ensemble to project regional climate and hydrological cycle changes.

Western Pacific emergent constraint lowers projected increase in Indian summer monsoon rainfall

NASA Astrophysics Data System (ADS)

Li, Gen; Xie, Shang-Ping; He, Chao; Chen, Zesheng

2017-10-01

The agrarian-based socioeconomic livelihood of densely populated South Asian countries is vulnerable to modest changes in Indian summer monsoon (ISM) rainfall. How the ISM rainfall will evolve is a question of broad scientific and socioeconomic importance. In response to increased greenhouse gas (GHG) forcing, climate models commonly project an increase in ISM rainfall. This wetter ISM projection, however, does not consider large model errors in both the mean state and ocean warming pattern. Here we identify a relationship between biases in simulated present climate and future ISM projections in a multi-model ensemble: models with excessive present-day precipitation over the tropical western Pacific tend to project a larger increase in ISM rainfall under GHG forcing because of too strong a negative cloud-radiation feedback on sea surface temperature. The excessive negative feedback suppresses the local ocean surface warming, strengthening ISM rainfall projections via atmospheric circulation. We calibrate the ISM rainfall projections using this `present-future relationship’ and observed western Pacific precipitation. The correction reduces by about 50% of the projected rainfall increase over the broad ISM region. Our study identifies an improved simulation of western Pacific convection as a priority for reliable ISM projections.

How Robust Is Your Project? From Local Failures to Global Catastrophes: A Complex Networks Approach to Project Systemic Risk.

PubMed

Ellinas, Christos; Allan, Neil; Durugbo, Christopher; Johansson, Anders

2015-01-01

Current societal requirements necessitate the effective delivery of complex projects that can do more while using less. Yet, recent large-scale project failures suggest that our ability to successfully deliver them is still at its infancy. Such failures can be seen to arise through various failure mechanisms; this work focuses on one such mechanism. Specifically, it examines the likelihood of a project sustaining a large-scale catastrophe, as triggered by single task failure and delivered via a cascading process. To do so, an analytical model was developed and tested on an empirical dataset by the means of numerical simulation. This paper makes three main contributions. First, it provides a methodology to identify the tasks most capable of impacting a project. In doing so, it is noted that a significant number of tasks induce no cascades, while a handful are capable of triggering surprisingly large ones. Secondly, it illustrates that crude task characteristics cannot aid in identifying them, highlighting the complexity of the underlying process and the utility of this approach. Thirdly, it draws parallels with systems encountered within the natural sciences by noting the emergence of self-organised criticality, commonly found within natural systems. These findings strengthen the need to account for structural intricacies of a project's underlying task precedence structure as they can provide the conditions upon which large-scale catastrophes materialise.

Diurnal Cycle of Convection and Interaction with the Large-Scale Circulation

NASA Technical Reports Server (NTRS)

Salby, Murry L.

2002-01-01

The science in this effort was scheduled in the project's third and fourth years, after a long record of high-resolution Global Cloud Imagery (GCI) had been produced. Unfortunately, political disruptions that interfered with this project led to its funding being terminated after only two years of support. Nevertheless, the availability of intermediate data opened the door to a number of important scientific studies. Beyond considerations of the diurnal cycle addressed in this grant, the GCI makes possible a wide range of studies surrounding convection, cloud, and precipitation. Several are already underway with colleagues in the US and abroad, including global cloud simulations, a global precipitation product, global precipitation simulations, upper tropospheric humidity, asynoptic sampling studies, convective organization studies, equatorial wave simulations, and the tropical tropopause.

MEchatronic REspiratory System SImulator for Neonatal Applications (MERESSINA) project: a novel bioengineering goal

PubMed Central

Scaramuzzo, Rosa T; Ciantelli, Massimiliano; Baldoli, Ilaria; Bellanti, Lisa; Gentile, Marzia; Cecchi, Francesca; Sigali, Emilio; Tognarelli, Selene; Ghirri, Paolo; Mazzoleni, Stefano; Menciassi, Arianna; Cuttano, Armando; Boldrini, Antonio; Laschi, Cecilia; Dario, Paolo

2013-01-01

Respiratory function is mandatory for extrauterine life, but is sometimes impaired in newborns due to prematurity, congenital malformations, or acquired pathologies. Mechanical ventilation is standard care, but long-term complications, such as bronchopulmonary dysplasia, are still largely reported. Therefore, continuous medical education is mandatory to correctly manage devices for assistance. Commercially available breathing function simulators are rarely suitable for the anatomical and physiological realities. The aim of this study is to develop a high-fidelity mechatronic simulator of neonatal airways and lungs for staff training and mechanical ventilator testing. The project is divided into three different phases: (1) a review study on respiratory physiology and pathophysiology and on already available single and multi-compartment models; (2) the prototyping phase; and (3) the on-field system validation. PMID:23966804

WE-AB-207A-12: HLCC Based Quantitative Evaluation Method of Image Artifact in Dental CBCT

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Y; Wu, S; Qi, H

Purpose: Image artifacts are usually evaluated qualitatively via visual observation of the reconstructed images, which is susceptible to subjective factors due to the lack of an objective evaluation criterion. In this work, we propose a Helgason-Ludwig consistency condition (HLCC) based evaluation method to quantify the severity level of different image artifacts in dental CBCT. Methods: Our evaluation method consists of four step: 1) Acquire Cone beam CT(CBCT) projection; 2) Convert 3D CBCT projection to fan-beam projection by extracting its central plane projection; 3) Convert fan-beam projection to parallel-beam projection utilizing sinogram-based rebinning algorithm or detail-based rebinning algorithm; 4) Obtain HLCCmore » profile by integrating parallel-beam projection per view and calculate wave percentage and variance of the HLCC profile, which can be used to describe the severity level of image artifacts. Results: Several sets of dental CBCT projections containing only one type of artifact (i.e. geometry, scatter, beam hardening, lag and noise artifact), were simulated using gDRR, a GPU tool developed for efficient, accurate, and realistic simulation of CBCT Projections. These simulated CBCT projections were used to test our proposed method. HLCC profile wave percentage and variance induced by geometry distortion are about 3∼21 times and 16∼393 times as large as that of the artifact-free projection, respectively. The increase factor of wave percentage and variance are 6 and133 times for beam hardening, 19 and 1184 times for scatter, and 4 and16 times for lag artifacts, respectively. In contrast, for noisy projection the wave percentage, variance and inconsistency level are almost the same with those of the noise-free one. Conclusion: We have proposed a quantitative evaluation method of image artifact based on HLCC theory. According to our simulation results, the severity of different artifact types is found to be in a following order: Scatter>Geometry>Beam hardening>Lag>Noise>Artifact-free in dental CBCT.« less

Incorporating variability in simulations of seasonally forced phenology using integral projection models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Goodsman, Devin W.; Aukema, Brian H.; McDowell, Nate G.

Phenology models are becoming increasingly important tools to accurately predict how climate change will impact the life histories of organisms. We propose a class of integral projection phenology models derived from stochastic individual-based models of insect development and demography.Our derivation, which is based on the rate-summation concept, produces integral projection models that capture the effect of phenotypic rate variability on insect phenology, but which are typically more computationally frugal than equivalent individual-based phenology models. We demonstrate our approach using a temperature-dependent model of the demography of the mountain pine beetle (Dendroctonus ponderosae Hopkins), an insect that kills mature pine trees.more » This work illustrates how a wide range of stochastic phenology models can be reformulated as integral projection models. Due to their computational efficiency, these integral projection models are suitable for deployment in large-scale simulations, such as studies of altered pest distributions under climate change.« less

Multi-Scale Simulations of Past and Future Projections of Hydrology in Lake Tahoe Basin, California-Nevada (Invited)

NASA Astrophysics Data System (ADS)

Niswonger, R. G.; Huntington, J. L.; Dettinger, M. D.; Rajagopal, S.; Gardner, M.; Morton, C. G.; Reeves, D. M.; Pohll, G. M.

2013-12-01

Water resources in the Tahoe basin are susceptible to long-term climate change and extreme events because it is a middle-altitude, snow-dominated basin that experiences large inter-annual climate variations. Lake Tahoe provides critical water supply for its basin and downstream populations, but changes in water supply are obscured by complex climatic and hydrologic gradients across the high relief, geologically complex basin. An integrated surface and groundwater model of the Lake Tahoe basin has been developed using GSFLOW to assess the effects of climate change and extreme events on surface and groundwater resources. Key hydrologic mechanisms are identified with this model that explains recent changes in water resources of the region. Critical vulnerabilities of regional water-supplies and hazards also were explored. Maintaining a balance between (a) accurate representation of spatial features (e.g., geology, streams, and topography) and hydrologic response (i.e., groundwater, stream, lake, and wetland flows and storages), and (b) computational efficiency, is a necessity for the desired model applications. Potential climatic influences on water resources are analyzed here in simulations of long-term water-availability and flood responses to selected 100-year climate-model projections. GSFLOW is also used to simulate a scenario depicting an especially extreme storm event that was constructed from a combination of two historical atmospheric-river storm events as part of the USGS MultiHazards Demonstration Project. Historical simulated groundwater levels, streamflow, wetlands, and lake levels compare well with measured values for a 30-year historical simulation period. Results are consistent for both small and large model grid cell sizes, due to the model's ability to represent water table altitude, streams, and other hydrologic features at the sub-grid scale. Simulated hydrologic responses are affected by climate change, where less groundwater resources will be available during more frequent droughts. Simulated floods for the region indicate issues related to drainage in the developed areas around Lake Tahoe, and necessary dam releases that create downstream flood risks.

The Importance of Simulation Workflow and Data Management in the Accelerated Climate Modeling for Energy Project

NASA Astrophysics Data System (ADS)

Bader, D. C.

2015-12-01

The Accelerated Climate Modeling for Energy (ACME) Project is concluding its first year. Supported by the Office of Science in the U.S. Department of Energy (DOE), its vision is to be "an ongoing, state-of-the-science Earth system modeling, modeling simulation and prediction project that optimizes the use of DOE laboratory resources to meet the science needs of the nation and the mission needs of DOE." Included in the "laboratory resources," is a large investment in computational, network and information technologies that will be utilized to both build better and more accurate climate models and broadly disseminate the data they generate. Current model diagnostic analysis and data dissemination technologies will not scale to the size of the simulations and the complexity of the models envisioned by ACME and other top tier international modeling centers. In this talk, the ACME Workflow component plans to meet these future needs will be described and early implementation examples will be highlighted.

Impact of large-scale circulation changes in the North Atlantic sector on the current and future Mediterranean winter hydroclimate

NASA Astrophysics Data System (ADS)

Barcikowska, Monika J.; Kapnick, Sarah B.; Feser, Frauke

2018-03-01

The Mediterranean region, located in the transition zone between the dry subtropical and wet European mid-latitude climate, is very sensitive to changes in the global mean climate state. Projecting future changes of the Mediterranean hydroclimate under global warming therefore requires dynamic climate models to reproduce the main mechanisms controlling regional hydroclimate with sufficiently high resolution to realistically simulate climate extremes. To assess future winter precipitation changes in the Mediterranean region we use the Geophysical Fluid Dynamics Laboratory high-resolution general circulation model for control simulations with pre-industrial greenhouse gas and aerosol concentrations which are compared to future scenario simulations. Here we show that the coupled model is able to reliably simulate the large-scale winter circulation, including the North Atlantic Oscillation and Eastern Atlantic patterns of variability, and its associated impacts on the mean Mediterranean hydroclimate. The model also realistically reproduces the regional features of daily heavy rainfall, which are absent in lower-resolution simulations. A five-member future projection ensemble, which assumes comparatively high greenhouse gas emissions (RCP8.5) until 2100, indicates a strong winter decline in Mediterranean precipitation for the coming decades. Consistent with dynamical and thermodynamical consequences of a warming atmosphere, derived changes feature a distinct bipolar behavior, i.e. wetting in the north—and drying in the south. Changes are most pronounced over the northwest African coast, where the projected winter precipitation decline reaches 40% of present values. Despite a decrease in mean precipitation, heavy rainfall indices show drastic increases across most of the Mediterranean, except the North African coast, which is under the strong influence of the cold Canary Current.

Impact of internal variability on projections of Sahel precipitation change

NASA Astrophysics Data System (ADS)

Monerie, Paul-Arthur; Sanchez-Gomez, Emilia; Pohl, Benjamin; Robson, Jon; Dong, Buwen

2017-11-01

The impact of the increase of greenhouse gases on Sahelian precipitation is very uncertain in both its spatial pattern and magnitude. In particular, the relative importance of internal variability versus external forcings depends on the time horizon considered in the climate projection. In this study we address the respective roles of the internal climate variability versus external forcings on Sahelian precipitation by using the data from the CESM Large Ensemble Project, which consists of a 40 member ensemble performed with the CESM1-CAM5 coupled model for the period 1920-2100. We show that CESM1-CAM5 is able to simulate the mean and interannual variability of Sahel precipitation, and is representative of a CMIP5 ensemble of simulations (i.e. it simulates the same pattern of precipitation change along with equivalent magnitude and seasonal cycle changes as the CMIP5 ensemble mean). However, CESM1-CAM5 underestimates the long-term decadal variability in Sahel precipitation. For short-term (2010-2049) and mid-term (2030-2069) projections the simulated internal variability component is able to obscure the projected impact of the external forcing. For long-term (2060-2099) projections external forcing induced change becomes stronger than simulated internal variability. Precipitation changes are found to be more robust over the central Sahel than over the western Sahel, where climate change effects struggle to emerge. Ten (thirty) members are needed to separate the 10 year averaged forced response from climate internal variability response in the western Sahel for a long-term (short-term) horizon. Over the central Sahel two members (ten members) are needed for a long-term (short-term) horizon.

Soapy: an adaptive optics simulation written purely in Python for rapid concept development

NASA Astrophysics Data System (ADS)

Reeves, Andrew

2016-07-01

Soapy is a newly developed Adaptive Optics (AO) simulation which aims be a flexible and fast to use tool-kit for many applications in the field of AO. It is written purely in the Python language, adding to and taking advantage of the already rich ecosystem of scientific libraries and programs. The simulation has been designed to be extremely modular, such that each component can be used stand-alone for projects which do not require a full end-to-end simulation. Ease of use, modularity and code clarity have been prioritised at the expense of computational performance. Though this means the code is not yet suitable for large studies of Extremely Large Telescope AO systems, it is well suited to education, exploration of new AO concepts and investigations of current generation telescopes.

15 MW HArdware-in-the-loop Grid Simulation Project

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rigas, Nikolaos; Fox, John Curtiss; Collins, Randy

2014-10-31

The 15MW Hardware-in-the-loop (HIL) Grid Simulator project was to (1) design, (2) construct and (3) commission a state-of-the-art grid integration testing facility for testing of multi-megawatt devices through a ‘shared facility’ model open to all innovators to promote the rapid introduction of new technology in the energy market to lower the cost of energy delivered. The 15 MW HIL Grid Simulator project now serves as the cornerstone of the Duke Energy Electric Grid Research, Innovation and Development (eGRID) Center. This project leveraged the 24 kV utility interconnection and electrical infrastructure of the US DOE EERE funded WTDTF project at themore » Clemson University Restoration Institute in North Charleston, SC. Additionally, the project has spurred interest from other technology sectors, including large PV inverter and energy storage testing and several leading edge research proposals dealing with smart grid technologies, grid modernization and grid cyber security. The key components of the project are the power amplifier units capable of providing up to 20MW of defined power to the research grid. The project has also developed a one of a kind solution to performing fault ride-through testing by combining a reactive divider network and a large power converter into a hybrid method. This unique hybrid method of performing fault ride-through analysis will allow for the research team at the eGRID Center to investigate the complex differences between the alternative methods of performing fault ride-through evaluations and will ultimately further the science behind this testing. With the final goal of being able to perform HIL experiments and demonstration projects, the eGRID team undertook a significant challenge with respect to developing a control system that is capable of communicating with several different pieces of equipment with different communication protocols in real-time. The eGRID team developed a custom fiber optical network that is based upon FPGA hardware that allows for communication between the key real-time interfaces and reduces the latency between these interfaces to acceptable levels for HIL experiments.« less

Real-time inextensible surgical thread simulation.

PubMed

Xu, Lang; Liu, Qian

2018-03-27

This paper discusses a real-time simulation method of inextensible surgical thread based on the Cosserat rod theory using position-based dynamics (PBD). The method realizes stable twining and knotting of surgical thread while including inextensibility, bending, twisting and coupling effects. The Cosserat rod theory is used to model the nonlinear elastic behavior of surgical thread. The surgical thread model is solved with PBD to achieve a real-time, extremely stable simulation. Due to the one-dimensional linear structure of surgical thread, the direct solution of the distance constraint based on tridiagonal matrix algorithm is used to enhance stretching resistance in every constraint projection iteration. In addition, continuous collision detection and collision response guarantee a large time step and high performance. Furthermore, friction is integrated into the constraint projection process to stabilize the twining of multiple threads and complex contact situations. Through comparisons with existing methods, the surgical thread maintains constant length under large deformation after applying the direct distance constraint in our method. The twining and knotting of multiple threads correspond to stable solutions to contact and friction forces. A surgical suture scene is also modeled to demonstrate the practicality and simplicity of our method. Our method achieves stable and fast simulation of inextensible surgical thread. Benefiting from the unified particle framework, the rigid body, elastic rod, and soft body can be simultaneously simulated. The method is appropriate for applications in virtual surgery that require multiple dynamic bodies.

Simulating the potential effects of climate change in two Colorado basins and at two Colorado ski areas

USGS Publications Warehouse

Battaglin, William; Hay, Lauren E.; Markstrom, Steve

2011-01-01

The mountainous areas of Colorado are used for tourism and recreation, and they provide water storage and supply for municipalities, industries, and agriculture. Recent studies suggest that water supply and tourist industries such as skiing are at risk from climate change. In this study, a distributed-parameter watershed model, the Precipitation-Runoff Modeling System (PRMS), is used to identify the potential effects of future climate on hydrologic conditions for two Colorado basins, the East River at Almont and the Yampa River at Steamboat Springs, and at the subbasin scale for two ski areas within those basins.Climate-change input files for PRMS were generated by modifying daily PRMS precipitation and temperature inputs with mean monthly climate-change fields of precipitation and temperature derived from five general circulation model (GCM) simulations using one current and three future carbon emission scenarios. All GCM simulations of mean daily minimum and maximum air temperature for the East and Yampa River basins indicate a relatively steady increase of up to several degrees Celsius from baseline conditions by 2094. GCM simulations of precipitation in the two basins indicate little change or trend in precipitation, but there is a large range associated with these projections. PRMS projections of basin mean daily streamflow vary by scenario but indicate a central tendency toward slight decreases, with a large range associated with these projections.Decreases in water content or changes in the spatial extent of snowpack in the East and Yampa River basins are important because of potential adverse effects on water supply and recreational activities. PRMS projections of each future scenario indicate a central tendency for decreases in basin mean snow-covered area and snowpack water equivalent, with the range in the projected decreases increasing with time. However, when examined on a monthly basis, the projected decreases are most dramatic during fall and spring. Presumably, ski area locations are picked because of a tendency to receive snow and keep snowpack relative to the surrounding area. This effect of ski area location within the basin was examined by comparing projections of March snow-covered area and snowpack water equivalent for the entire basin with more local projections for the portion of the basin that represents the ski area in the PRMS models. These projections indicate a steady decrease in March snow-covered area for the basins but only small changes in March snow-covered area at both ski areas for the three future scenarios until around 2050. After 2050, larger decreases are possible, but there is a large range in the projections of future scenarios. The rates of decrease for snowpack water equivalent and precipitation that falls as snow are similar at the basin and subbasin scale in both basins. Results from this modeling effort show that there is a wide range of possible outcomes for future snowpack conditions in Colorado. The results also highlight the differences between projections for entire basins and projections for local areas or subbasins within those basins.

The formation of disc galaxies in high-resolution moving-mesh cosmological simulations

NASA Astrophysics Data System (ADS)

Marinacci, Federico; Pakmor, Rüdiger; Springel, Volker

2014-01-01

We present cosmological hydrodynamical simulations of eight Milky Way-sized haloes that have been previously studied with dark matter only in the Aquarius project. For the first time, we employ the moving-mesh code AREPO in zoom simulations combined with a comprehensive model for galaxy formation physics designed for large cosmological simulations. Our simulations form in most of the eight haloes strongly disc-dominated systems with realistic rotation curves, close to exponential surface density profiles, a stellar mass to halo mass ratio that matches expectations from abundance matching techniques, and galaxy sizes and ages consistent with expectations from large galaxy surveys in the local Universe. There is no evidence for any dark matter core formation in our simulations, even so they include repeated baryonic outflows by supernova-driven winds and black hole quasar feedback. For one of our haloes, the object studied in the recent `Aquila' code comparison project, we carried out a resolution study with our techniques, covering a dynamic range of 64 in mass resolution. Without any change in our feedback parameters, the final galaxy properties are reassuringly similar, in contrast to other modelling techniques used in the field that are inherently resolution dependent. This success in producing realistic disc galaxies is reached, in the context of our interstellar medium treatment, without resorting to a high density threshold for star formation, a low star formation efficiency, or early stellar feedback, factors deemed crucial for disc formation by other recent numerical studies.

Validating the simulation of large-scale parallel applications using statistical characteristics

DOE PAGES

Zhang, Deli; Wilke, Jeremiah; Hendry, Gilbert; ...

2016-03-01

Simulation is a widely adopted method to analyze and predict the performance of large-scale parallel applications. Validating the hardware model is highly important for complex simulations with a large number of parameters. Common practice involves calculating the percent error between the projected and the real execution time of a benchmark program. However, in a high-dimensional parameter space, this coarse-grained approach often suffers from parameter insensitivity, which may not be known a priori. Moreover, the traditional approach cannot be applied to the validation of software models, such as application skeletons used in online simulations. In this work, we present a methodologymore » and a toolset for validating both hardware and software models by quantitatively comparing fine-grained statistical characteristics obtained from execution traces. Although statistical information has been used in tasks like performance optimization, this is the first attempt to apply it to simulation validation. Lastly, our experimental results show that the proposed evaluation approach offers significant improvement in fidelity when compared to evaluation using total execution time, and the proposed metrics serve as reliable criteria that progress toward automating the simulation tuning process.« less

A TREETOPS simulation of the Hubble Space Telescope-High Gain Antenna interaction

NASA Technical Reports Server (NTRS)

Sharkey, John P.

1987-01-01

Virtually any project dealing with the control of a Large Space Structure (LSS) will involve some level of verification by digital computer simulation. While the Hubble Space Telescope might not normally be included in a discussion of LSS, it is presented to highlight a recently developed simulation and analysis program named TREETOPS. TREETOPS provides digital simulation, linearization, and control system interaction of flexible, multibody spacecraft which admit to a point-connected tree topology. The HST application of TREETOPS is intended to familiarize the LSS community with TREETOPS by presenting a user perspective of its key features.

Final Report: Large-Scale Optimization for Bayesian Inference in Complex Systems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ghattas, Omar

2013-10-15

The SAGUARO (Scalable Algorithms for Groundwater Uncertainty Analysis and Robust Optimiza- tion) Project focuses on the development of scalable numerical algorithms for large-scale Bayesian inversion in complex systems that capitalize on advances in large-scale simulation-based optimiza- tion and inversion methods. Our research is directed in three complementary areas: efficient approximations of the Hessian operator, reductions in complexity of forward simulations via stochastic spectral approximations and model reduction, and employing large-scale optimization concepts to accelerate sampling. Our efforts are integrated in the context of a challenging testbed problem that considers subsurface reacting flow and transport. The MIT component of the SAGUAROmore » Project addresses the intractability of conventional sampling methods for large-scale statistical inverse problems by devising reduced-order models that are faithful to the full-order model over a wide range of parameter values; sampling then employs the reduced model rather than the full model, resulting in very large computational savings. Results indicate little effect on the computed posterior distribution. On the other hand, in the Texas-Georgia Tech component of the project, we retain the full-order model, but exploit inverse problem structure (adjoint-based gradients and partial Hessian information of the parameter-to- observation map) to implicitly extract lower dimensional information on the posterior distribution; this greatly speeds up sampling methods, so that fewer sampling points are needed. We can think of these two approaches as "reduce then sample" and "sample then reduce." In fact, these two approaches are complementary, and can be used in conjunction with each other. Moreover, they both exploit deterministic inverse problem structure, in the form of adjoint-based gradient and Hessian information of the underlying parameter-to-observation map, to achieve their speedups.« less

Reliability ensemble averaging of 21st century projections of terrestrial net primary productivity reduces global and regional uncertainties

NASA Astrophysics Data System (ADS)

Exbrayat, Jean-François; Bloom, A. Anthony; Falloon, Pete; Ito, Akihiko; Smallman, T. Luke; Williams, Mathew

2018-02-01

Multi-model averaging techniques provide opportunities to extract additional information from large ensembles of simulations. In particular, present-day model skill can be used to evaluate their potential performance in future climate simulations. Multi-model averaging methods have been used extensively in climate and hydrological sciences, but they have not been used to constrain projected plant productivity responses to climate change, which is a major uncertainty in Earth system modelling. Here, we use three global observationally orientated estimates of current net primary productivity (NPP) to perform a reliability ensemble averaging (REA) method using 30 global simulations of the 21st century change in NPP based on the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP) business as usual emissions scenario. We find that the three REA methods support an increase in global NPP by the end of the 21st century (2095-2099) compared to 2001-2005, which is 2-3 % stronger than the ensemble ISIMIP mean value of 24.2 Pg C y-1. Using REA also leads to a 45-68 % reduction in the global uncertainty of 21st century NPP projection, which strengthens confidence in the resilience of the CO2 fertilization effect to climate change. This reduction in uncertainty is especially clear for boreal ecosystems although it may be an artefact due to the lack of representation of nutrient limitations on NPP in most models. Conversely, the large uncertainty that remains on the sign of the response of NPP in semi-arid regions points to the need for better observations and model development in these regions.

Evolution of Precipitation Extremes in Three Large Ensembles of Climate Simulations - Impact of Spatial and Temporal Resolutions

NASA Astrophysics Data System (ADS)

Martel, J. L.; Brissette, F.; Mailhot, A.; Wood, R. R.; Ludwig, R.; Frigon, A.; Leduc, M.; Turcotte, R.

2017-12-01

Recent studies indicate that the frequency and intensity of extreme precipitation will increase in future climate due to global warming. In this study, we compare annual maxima precipitation series from three large ensembles of climate simulations at various spatial and temporal resolutions. The first two are at the global scale: the Canadian Earth System Model (CanESM2) 50-member large ensemble (CanESM2-LE) at a 2.8° resolution and the Community Earth System Model (CESM1) 40-member large ensemble (CESM1-LE) at a 1° resolution. The third ensemble is at the regional scale over both Eastern North America and Europe: the Canadian Regional Climate Model (CRCM5) 50-member large ensemble (CRCM5-LE) at a 0.11° resolution, driven at its boundaries by the CanESM-LE. The CRCM5-LE is a new ensemble issued from the ClimEx project (http://www.climex-project.org), a Québec-Bavaria collaboration. Using these three large ensembles, change in extreme precipitations over the historical (1980-2010) and future (2070-2100) periods are investigated. This results in 1 500 (30 years x 50 members for CanESM2-LE and CRCM5-LE) and 1200 (30 years x 40 members for CESM1-LE) simulated years over both the historical and future periods. Using these large datasets, the empirical daily (and sub-daily for CRCM5-LE) extreme precipitation quantiles for large return periods ranging from 2 to 100 years are computed. Results indicate that daily extreme precipitations generally will increase over most land grid points of both domains according to the three large ensembles. Regarding the CRCM5-LE, the increase in sub-daily extreme precipitations will be even more important than the one observed for daily extreme precipitations. Considering that many public infrastructures have lifespans exceeding 75 years, the increase in extremes has important implications on service levels of water infrastructures and public safety.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Panebianco, S.; Dore, D.; Giomataris, I.

Time Projection Chambers are widely used since many years for tracking and identification of charged particles in high energy physics. We present a new R and D project to investigate the feasibility of a Micromegas TPC for low energy heavy ions detection. Two physics cases are relevant for this project. The first is the study of the nuclear fission of actinides by measuring the fission fragments properties (mass, nuclear charge, kinetic energy) that will be performed at different installations and in particular at the NFS facility to be built in the framework of the SPIRAL2 project in GANIL. The secondmore » physics case is the study of heavy ion reactions, like ({alpha},{gamma}), ({alpha},p), ({alpha},n) and all the inverse reactions in the energy range between 1.5 and 3 AMeV using both stable and radioactive beams. These reactions have a key role in p process in nuclear astrophysics to explain the synthesis of heavy proton-rich nuclei. Within the project, a large effort is devoted to Monte-Carlo simulations and a detailed benchmark of different simulation codes on the energy loss and range in gas of heavy ions at low energy has been performed. A new approach for simulating the ion charge state evolution in GEANT4 is also presented. Finally, preliminary results of an experimental test campaign on prototype are discussed.« less

Acceleration of Monte Carlo SPECT simulation using convolution-based forced detection

NASA Astrophysics Data System (ADS)

de Jong, H. W. A. M.; Slijpen, E. T. P.; Beekman, F. J.

2001-02-01

Monte Carlo (MC) simulation is an established tool to calculate photon transport through tissue in Emission Computed Tomography (ECT). Since the first appearance of MC a large variety of variance reduction techniques (VRT) have been introduced to speed up these notoriously slow simulations. One example of a very effective and established VRT is known as forced detection (FD). In standard FD the path from the photon's scatter position to the camera is chosen stochastically from the appropriate probability density function (PDF), modeling the distance-dependent detector response. In order to speed up MC the authors propose a convolution-based FD (CFD) which involves replacing the sampling of the PDF by a convolution with a kernel which depends on the position of the scatter event. The authors validated CFD for parallel-hole Single Photon Emission Computed Tomography (SPECT) using a digital thorax phantom. Comparison of projections estimated with CFD and standard FD shows that both estimates converge to practically identical projections (maximum bias 0.9% of peak projection value), despite the slightly different photon paths used in CFD and standard FD. Projections generated with CFD converge, however, to a noise-free projection up to one or two orders of magnitude faster, which is extremely useful in many applications such as model-based image reconstruction.

Projected Response of Low-Level Convergence and Associated Precipitation to Greenhouse Warming

NASA Astrophysics Data System (ADS)

Weller, Evan; Jakob, Christian; Reeder, Michael J.

2017-10-01

The parameterization of convection in climate models is a large source of uncertainty in projecting future precipitation changes. Here an objective method to identify organized low-level convergence lines has been used to better understand how atmospheric convection is organized and projected to change, as low-level convergence plays an important role in the processes leading to precipitation. The frequency and strength of convergence lines over both ocean and land in current climate simulations is too low compared to reanalysis data. Projections show a further reduction in the frequency and strength of convergence lines over the midlatitudes. In the tropics, the largest changes in frequency are generally associated with shifts in major low-latitude convergence zones, consistent with changes in the precipitation. Further, examining convergence lines when in the presence or absence of precipitation results in large spatial contrasts, providing a better understanding of regional changes in terms of thermodynamic and dynamic effects.

The global gridded crop model intercomparison: Data and modeling protocols for Phase 1 (v1.0)

DOE PAGES

Elliott, J.; Müller, C.; Deryng, D.; ...

2015-02-11

We present protocols and input data for Phase 1 of the Global Gridded Crop Model Intercomparison, a project of the Agricultural Model Intercomparison and Improvement Project (AgMIP). The project consist of global simulations of yields, phenologies, and many land-surface fluxes using 12–15 modeling groups for many crops, climate forcing data sets, and scenarios over the historical period from 1948 to 2012. The primary outcomes of the project include (1) a detailed comparison of the major differences and similarities among global models commonly used for large-scale climate impact assessment, (2) an evaluation of model and ensemble hindcasting skill, (3) quantification ofmore » key uncertainties from climate input data, model choice, and other sources, and (4) a multi-model analysis of the agricultural impacts of large-scale climate extremes from the historical record.« less

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

DOE PAGES

Haarsma, Reindert J.; Roberts, Malcolm J.; Vidale, Pier Luigi; ...

2016-11-22

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relativelymore » few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. Lastly, HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.« less

Quantitative analysis of artifacts in 4D DSA: the relative contributions of beam hardening and scatter to vessel dropout behind highly attenuating structures

NASA Astrophysics Data System (ADS)

Hermus, James; Szczykutowicz, Timothy P.; Strother, Charles M.; Mistretta, Charles

2014-03-01

When performing Computed Tomographic (CT) image reconstruction on digital subtraction angiography (DSA) projections, loss of vessel contrast has been observed behind highly attenuating anatomy, such as dental implants and large contrast filled aneurysms. Because this typically occurs only in a limited range of projection angles, the observed contrast time course can potentially be altered. In this work, we have developed a model for acquiring DSA projections that models both the polychromatic nature of the x-ray spectrum and the x-ray scattering interactions to investigate this problem. In our simulation framework, scatter and beam hardening contributions to vessel dropout can be analyzed separately. We constructed digital phantoms with large clearly defined regions containing iodine contrast, bone, soft issue, titanium (dental implants) or combinations of these materials. As the regions containing the materials were large and rectangular, when the phantoms were forward projected, the projections contained uniform regions of interest (ROI) and enabled accurate vessel dropout analysis. Two phantom models were used, one to model the case of a vessel behind a large contrast filled aneurysm and the other to model a vessel behind a dental implant. Cases in which both beam hardening and scatter were turned off, only scatter was turned on, only beam hardening was turned on, and both scatter and beam hardening were turned on, were simulated for both phantom models. The analysis of this data showed that the contrast degradation is primarily due to scatter. When analyzing the aneurysm case, 90.25% of the vessel contrast was lost in the polychromatic scatter image, however only 50.5% of the vessel contrast was lost in the beam hardening only image. When analyzing the teeth case, 44.2% of the vessel contrast was lost in the polychromatic scatter image and only 26.2% of the vessel contrast was lost in the beam hardening only image.

Numerical Simulations of Airflows and Tracer Transport in the Southwestern United States.

NASA Astrophysics Data System (ADS)

Yamada, Tetsuji

2000-03-01

Project MOHAVE (Measurement of Haze and Visual Effects) produced a unique set of tracer data over the southwestern United States. During the summer of 1992, a perfluorocarbon tracer gas was released from the Mohave Power Project (MPP), a large coal-fired facility in southern Nevada. Three-dimensional atmospheric models, the Higher-Order Turbulence Model for Atmospheric Circulation-Random Puff Transport and Diffusion (HOTMAC-RAPTAD), were used to simulate the concentrations of tracer gas that were observed during a portion of the summer intensive period of Project MOHAVE. The study area extended from northwestern Arizona to southern Nevada and included Lake Mead, the Colorado River Valley, the Grand Canyon National Park, and MPP. The computational domain was 368 km in the east-west direction by 252 km in the north-south direction. Rawinsonde and radar wind profiler data were used to provide initial and boundary conditions to HOTMAC simulations. HOTMAC with a horizontal grid spacing of 4 km was able to simulate the diurnal variations of drainage and upslope flows along the Grand Canyon and Colorado River Valley. HOTMAC also captured the diurnal variations of turbulence, which played important roles for the transport and diffusion simulations by RAPTAD. The modeled tracer gas concentrations were compared with observations. The model's performance was evaluated statistically.

Observations give us CLUES to Cosmic Flows' origins

NASA Astrophysics Data System (ADS)

Sorce, Jenny; Courtois, H.; Gottloeber, S.; Hoffman, Y.; Pomarede, D.; Tully, R. B.; Flows, Cosmic; CLUES

2014-01-01

In an era where the wealth of telescope-data and the development of computer superclusters keep increasing, the knowledge of Large Scale Structures' formation and evolution constitutes a tremendous challenge. Within this context the project Cosmic Flows has recently produced a catalog of peculiar velocities up to 150 Mpc. These velocities, obtained from direct distance measurements, are ideal markers of the underlying gravitational potential. They form a fantastic input to perform constrained simulations of the Local Universe within the CLUES project. A new method has recently been elaborated to achieve these simulations which prove to be excellent replicas of our neighborhood. The Wiener-Filter, the Reverse Zel'dovich Approximation and the Constrained Realization techniques are combined to build Initial Conditions. The resulting second generation of constrained simulations presents us the formidable history of the Great Attractor's and nearby supercluster's formation.

High resolution global climate modelling; the UPSCALE project, a large simulation campaign

NASA Astrophysics Data System (ADS)

Mizielinski, M. S.; Roberts, M. J.; Vidale, P. L.; Schiemann, R.; Demory, M.-E.; Strachan, J.; Edwards, T.; Stephens, A.; Lawrence, B. N.; Pritchard, M.; Chiu, P.; Iwi, A.; Churchill, J.; del Cano Novales, C.; Kettleborough, J.; Roseblade, W.; Selwood, P.; Foster, M.; Glover, M.; Malcolm, A.

2014-01-01

The UPSCALE (UK on PRACE: weather-resolving Simulations of Climate for globAL Environmental risk) project constructed and ran an ensemble of HadGEM3 (Hadley centre Global Environment Model 3) atmosphere-only global climate simulations over the period 1985-2011, at resolutions of N512 (25 km), N216 (60 km) and N96 (130 km) as used in current global weather forecasting, seasonal prediction and climate modelling respectively. Alongside these present climate simulations a parallel ensemble looking at extremes of future climate was run, using a time-slice methodology to consider conditions at the end of this century. These simulations were primarily performed using a 144 million core hour, single year grant of computing time from PRACE (the Partnership for Advanced Computing in Europe) in 2012, with additional resources supplied by the Natural Environmental Research Council (NERC) and the Met Office. Almost 400 terabytes of simulation data were generated on the HERMIT supercomputer at the high performance computing center Stuttgart (HLRS), and transferred to the JASMIN super-data cluster provided by the Science and Technology Facilities Council Centre for Data Archival (STFC CEDA) for analysis and storage. In this paper we describe the implementation of the project, present the technical challenges in terms of optimisation, data output, transfer and storage that such a project involves and include details of the model configuration and the composition of the UPSCALE dataset. This dataset is available for scientific analysis to allow assessment of the value of model resolution in both present and potential future climate conditions.

High-resolution global climate modelling: the UPSCALE project, a large-simulation campaign

NASA Astrophysics Data System (ADS)

Mizielinski, M. S.; Roberts, M. J.; Vidale, P. L.; Schiemann, R.; Demory, M.-E.; Strachan, J.; Edwards, T.; Stephens, A.; Lawrence, B. N.; Pritchard, M.; Chiu, P.; Iwi, A.; Churchill, J.; del Cano Novales, C.; Kettleborough, J.; Roseblade, W.; Selwood, P.; Foster, M.; Glover, M.; Malcolm, A.

2014-08-01

The UPSCALE (UK on PRACE: weather-resolving Simulations of Climate for globAL Environmental risk) project constructed and ran an ensemble of HadGEM3 (Hadley Centre Global Environment Model 3) atmosphere-only global climate simulations over the period 1985-2011, at resolutions of N512 (25 km), N216 (60 km) and N96 (130 km) as used in current global weather forecasting, seasonal prediction and climate modelling respectively. Alongside these present climate simulations a parallel ensemble looking at extremes of future climate was run, using a time-slice methodology to consider conditions at the end of this century. These simulations were primarily performed using a 144 million core hour, single year grant of computing time from PRACE (the Partnership for Advanced Computing in Europe) in 2012, with additional resources supplied by the Natural Environment Research Council (NERC) and the Met Office. Almost 400 terabytes of simulation data were generated on the HERMIT supercomputer at the High Performance Computing Center Stuttgart (HLRS), and transferred to the JASMIN super-data cluster provided by the Science and Technology Facilities Council Centre for Data Archival (STFC CEDA) for analysis and storage. In this paper we describe the implementation of the project, present the technical challenges in terms of optimisation, data output, transfer and storage that such a project involves and include details of the model configuration and the composition of the UPSCALE data set. This data set is available for scientific analysis to allow assessment of the value of model resolution in both present and potential future climate conditions.

Spatially distributed potential evapotranspiration modeling and climate projections.

PubMed

Gharbia, Salem S; Smullen, Trevor; Gill, Laurence; Johnston, Paul; Pilla, Francesco

2018-08-15

Evapotranspiration integrates energy and mass transfer between the Earth's surface and atmosphere and is the most active mechanism linking the atmosphere, hydrosphsophere, lithosphere and biosphere. This study focuses on the fine resolution modeling and projection of spatially distributed potential evapotranspiration on the large catchment scale as response to climate change. Six potential evapotranspiration designed algorithms, systematically selected based on a structured criteria and data availability, have been applied and then validated to long-term mean monthly data for the Shannon River catchment with a 50m 2 cell size. The best validated algorithm was therefore applied to evaluate the possible effect of future climate change on potential evapotranspiration rates. Spatially distributed potential evapotranspiration projections have been modeled based on climate change projections from multi-GCM ensembles for three future time intervals (2020, 2050 and 2080) using a range of different Representative Concentration Pathways producing four scenarios for each time interval. Finally, seasonal results have been compared to baseline results to evaluate the impact of climate change on the potential evapotranspiration and therefor on the catchment dynamical water balance. The results present evidence that the modeled climate change scenarios would have a significant impact on the future potential evapotranspiration rates. All the simulated scenarios predicted an increase in potential evapotranspiration for each modeled future time interval, which would significantly affect the dynamical catchment water balance. This study addresses the gap in the literature of using GIS-based algorithms to model fine-scale spatially distributed potential evapotranspiration on the large catchment systems based on climatological observations and simulations in different climatological zones. Providing fine-scale potential evapotranspiration data is very crucial to assess the dynamical catchment water balance to setup management scenarios for the water abstractions. This study illustrates a transferable systematic method to design GIS-based algorithms to simulate spatially distributed potential evapotranspiration on the large catchment systems. Copyright © 2018 Elsevier B.V. All rights reserved.

The Building Blocks for JWST I and T (Integrations and Test) to Operations - From Simulator to Flight Units

NASA Technical Reports Server (NTRS)

Fatig, Curtis; Ochs, William; Johns, Alan; Seaton, Bonita; Adams, Cynthia; Wasiak, Francis; Jones, Ronald; Jackson, Wallace

2012-01-01

The James Webb Space Telescope (JWST) Project has an extended integration and test (I&T) phase due to long procurement and development times of various components as well as recent launch delays. The JWST Ground Segment and Operations group has developed a roadmap of the various ground and flight elements and their use in the various JWST I&T test programs. The JWST Project s building block approach to the eventual operational systems, while not new, is complex and challenging; a large-scale mission like JWST involves international partners, many vendors across the United States, and competing needs for the same systems. One of the challenges is resource balancing so simulators and flight products for various elements congeal into integrated systems used for I&T and flight operations activities. This building block approach to an incremental buildup provides for early problem identification with simulators and exercises the flight operations systems, products, and interfaces during the JWST I&T test programs. The JWST Project has completed some early I&T with the simulators, engineering models and some components of the operational ground system. The JWST Project is testing the various flight units as they are delivered and will continue to do so for the entire flight and operational system. The JWST Project has already and will continue to reap the value of the building block approach on the road to launch and flight operations.

Simulating the impact of the large-scale circulation on the 2-m temperature and precipitation climatology

NASA Astrophysics Data System (ADS)

Bowden, Jared H.; Nolte, Christopher G.; Otte, Tanya L.

2013-04-01

The impact of the simulated large-scale atmospheric circulation on the regional climate is examined using the Weather Research and Forecasting (WRF) model as a regional climate model. The purpose is to understand the potential need for interior grid nudging for dynamical downscaling of global climate model (GCM) output for air quality applications under a changing climate. In this study we downscale the NCEP-Department of Energy Atmospheric Model Intercomparison Project (AMIP-II) Reanalysis using three continuous 20-year WRF simulations: one simulation without interior grid nudging and two using different interior grid nudging methods. The biases in 2-m temperature and precipitation for the simulation without interior grid nudging are unreasonably large with respect to the North American Regional Reanalysis (NARR) over the eastern half of the contiguous United States (CONUS) during the summer when air quality concerns are most relevant. This study examines how these differences arise from errors in predicting the large-scale atmospheric circulation. It is demonstrated that the Bermuda high, which strongly influences the regional climate for much of the eastern half of the CONUS during the summer, is poorly simulated without interior grid nudging. In particular, two summers when the Bermuda high was west (1993) and east (2003) of its climatological position are chosen to illustrate problems in the large-scale atmospheric circulation anomalies. For both summers, WRF without interior grid nudging fails to simulate the placement of the upper-level anticyclonic (1993) and cyclonic (2003) circulation anomalies. The displacement of the large-scale circulation impacts the lower atmosphere moisture transport and precipitable water, affecting the convective environment and precipitation. Using interior grid nudging improves the large-scale circulation aloft and moisture transport/precipitable water anomalies, thereby improving the simulated 2-m temperature and precipitation. The results demonstrate that constraining the RCM to the large-scale features in the driving fields improves the overall accuracy of the simulated regional climate, and suggest that in the absence of such a constraint, the RCM will likely misrepresent important large-scale shifts in the atmospheric circulation under a future climate.

Evaluating synoptic systems in the CMIP5 climate models over the Australian region

NASA Astrophysics Data System (ADS)

Gibson, Peter B.; Uotila, Petteri; Perkins-Kirkpatrick, Sarah E.; Alexander, Lisa V.; Pitman, Andrew J.

2016-10-01

Climate models are our principal tool for generating the projections used to inform climate change policy. Our confidence in projections depends, in part, on how realistically they simulate present day climate and associated variability over a range of time scales. Traditionally, climate models are less commonly assessed at time scales relevant to daily weather systems. Here we explore the utility of a self-organizing maps (SOMs) procedure for evaluating the frequency, persistence and transitions of daily synoptic systems in the Australian region simulated by state-of-the-art global climate models. In terms of skill in simulating the climatological frequency of synoptic systems, large spread was observed between models. A positive association between all metrics was found, implying that relative skill in simulating the persistence and transitions of systems is related to skill in simulating the climatological frequency. Considering all models and metrics collectively, model performance was found to be related to model horizontal resolution but unrelated to vertical resolution or representation of the stratosphere. In terms of the SOM procedure, the timespan over which evaluation was performed had some influence on model performance skill measures, as did the number of circulation types examined. These findings have implications for selecting models most useful for future projections over the Australian region, particularly for projections related to synoptic scale processes and phenomena. More broadly, this study has demonstrated the utility of the SOMs procedure in providing a process-based evaluation of climate models.

Climate SPHINX: High-resolution present-day and future climate simulations with an improved representation of small-scale variability

NASA Astrophysics Data System (ADS)

Davini, Paolo; von Hardenberg, Jost; Corti, Susanna; Subramanian, Aneesh; Weisheimer, Antje; Christensen, Hannah; Juricke, Stephan; Palmer, Tim

2016-04-01

The PRACE Climate SPHINX project investigates the sensitivity of climate simulations to model resolution and stochastic parameterization. The EC-Earth Earth-System Model is used to explore the impact of stochastic physics in 30-years climate integrations as a function of model resolution (from 80km up to 16km for the atmosphere). The experiments include more than 70 simulations in both a historical scenario (1979-2008) and a climate change projection (2039-2068), using RCP8.5 CMIP5 forcing. A total amount of 20 million core hours will be used at end of the project (March 2016) and about 150 TBytes of post-processed data will be available to the climate community. Preliminary results show a clear improvement in the representation of climate variability over the Euro-Atlantic following resolution increase. More specifically, the well-known atmospheric blocking negative bias over Europe is definitely resolved. High resolution runs also show improved fidelity in representation of tropical variability - such as the MJO and its propagation - over the low resolution simulations. It is shown that including stochastic parameterization in the low resolution runs help to improve some of the aspects of the MJO propagation further. These findings show the importance of representing the impact of small scale processes on the large scale climate variability either explicitly (with high resolution simulations) or stochastically (in low resolution simulations).

Cretin Memory Flow on Sierra

DOE Office of Scientific and Technical Information (OSTI.GOV)

Langer, S. H.; Scott, H. A.

2016-08-05

The Cretin iCOE project has a goal of enabling the efficient generation of Non-LTE opacities for use in radiation-hydrodynamic simulation codes using the Nvidia boards on LLNL’s upcoming Sierra system. Achieving the desired level of accuracy for some simulations require the use of a vary large number of atomic configurations (a configuration includes the atomic level for all electrons and how they are coupled together). The NLTE rate matrix needs to be solved separately in each zone. Calculating NLTE opacities can consume more time than all other physics packages used in a simulation.

Convolutional neural networks applied to neutrino events in a liquid argon time projection chamber

NASA Astrophysics Data System (ADS)

Acciarri, R.; Adams, C.; An, R.; Asaadi, J.; Auger, M.; Bagby, L.; Baller, B.; Barr, G.; Bass, M.; Bay, F.; Bishai, M.; Blake, A.; Bolton, T.; Bugel, L.; Camilleri, L.; Caratelli, D.; Carls, B.; Castillo Fernandez, R.; Cavanna, F.; Chen, H.; Church, E.; Cianci, D.; Collin, G. H.; Conrad, J. M.; Convery, M.; Crespo-Anadón, J. I.; Del Tutto, M.; Devitt, D.; Dytman, S.; Eberly, B.; Ereditato, A.; Escudero Sanchez, L.; Esquivel, J.; Fleming, B. T.; Foreman, W.; Furmanski, A. P.; Garvey, G. T.; Genty, V.; Goeldi, D.; Gollapinni, S.; Graf, N.; Gramellini, E.; Greenlee, H.; Grosso, R.; Guenette, R.; Hackenburg, A.; Hamilton, P.; Hen, O.; Hewes, J.; Hill, C.; Ho, J.; Horton-Smith, G.; James, C.; de Vries, J. Jan; Jen, C.-M.; Jiang, L.; Johnson, R. A.; Jones, B. J. P.; Joshi, J.; Jostlein, H.; Kaleko, D.; Karagiorgi, G.; Ketchum, W.; Kirby, B.; Kirby, M.; Kobilarcik, T.; Kreslo, I.; Laube, A.; Li, Y.; Lister, A.; Littlejohn, B. R.; Lockwitz, S.; Lorca, D.; Louis, W. C.; Luethi, M.; Lundberg, B.; Luo, X.; Marchionni, A.; Mariani, C.; Marshall, J.; Martinez Caicedo, D. A.; Meddage, V.; Miceli, T.; Mills, G. B.; Moon, J.; Mooney, M.; Moore, C. D.; Mousseau, J.; Murrells, R.; Naples, D.; Nienaber, P.; Nowak, J.; Palamara, O.; Paolone, V.; Papavassiliou, V.; Pate, S. F.; Pavlovic, Z.; Porzio, D.; Pulliam, G.; Qian, X.; Raaf, J. L.; Rafique, A.; Rochester, L.; von Rohr, C. Rudolf; Russell, B.; Schmitz, D. W.; Schukraft, A.; Seligman, W.; Shaevitz, M. H.; Sinclair, J.; Snider, E. L.; Soderberg, M.; Söldner-Rembold, S.; Soleti, S. R.; Spentzouris, P.; Spitz, J.; St. John, J.; Strauss, T.; Szelc, A. M.; Tagg, N.; Terao, K.; Thomson, M.; Toups, M.; Tsai, Y.-T.; Tufanli, S.; Usher, T.; Van de Water, R. G.; Viren, B.; Weber, M.; Weston, J.; Wickremasinghe, D. A.; Wolbers, S.; Wongjirad, T.; Woodruff, K.; Yang, T.; Zeller, G. P.; Zennamo, J.; Zhang, C.

2017-03-01

We present several studies of convolutional neural networks applied to data coming from the MicroBooNE detector, a liquid argon time projection chamber (LArTPC). The algorithms studied include the classification of single particle images, the localization of single particle and neutrino interactions in an image, and the detection of a simulated neutrino event overlaid with cosmic ray backgrounds taken from real detector data. These studies demonstrate the potential of convolutional neural networks for particle identification or event detection on simulated neutrino interactions. We also address technical issues that arise when applying this technique to data from a large LArTPC at or near ground level.

The UK Earth System Model project

NASA Astrophysics Data System (ADS)

Tang, Yongming

2016-04-01

In this talk we will describe the development and current status of the UK Earth System Model (UKESM). This project is a NERC/Met Office collaboration and has two objectives; to develop and apply a world-leading Earth System Model, and to grow a community of UK Earth System Model scientists. We are building numerical models that include all the key components of the global climate system, and contain the important process interactions between global biogeochemistry, atmospheric chemistry and the physical climate system. UKESM will be used to make key CMIP6 simulations as well as long-time (e.g. millennium) simulations, large ensemble experiments and investigating a range of future carbon emission scenarios.

Center for Extended Magnetohydrodynamics Modeling

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ramos, Jesus

This researcher participated in the DOE-funded Center for Extended Magnetohydrodynamics Modeling (CEMM), a multi-institutional collaboration led by the Princeton Plasma Physics Laboratory with Dr. Stephen Jardin as the overall Principal Investigator. This project developed advanced simulation tools to study the non-linear macroscopic dynamics of magnetically confined plasmas. The collaborative effort focused on the development of two large numerical simulation codes, M3D-C1 and NIMROD, and their application to a wide variety of problems. Dr. Ramos was responsible for theoretical aspects of the project, deriving consistent sets of model equations applicable to weakly collisional plasmas and devising test problems for verification ofmore » the numerical codes. This activity was funded for twelve years.« less

Panoramic cone beam computed tomography

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chang Jenghwa; Zhou Lili; Wang Song

2012-05-15

Purpose: Cone-beam computed tomography (CBCT) is the main imaging tool for image-guided radiotherapy but its functionality is limited by a small imaging volume and restricted image position (imaged at the central instead of the treatment position for peripheral lesions to avoid collisions). In this paper, the authors present the concept of ''panoramic CBCT,'' which can image patients at the treatment position with an imaging volume as large as practically needed. Methods: In this novel panoramic CBCT technique, the target is scanned sequentially from multiple view angles. For each view angle, a half scan (180 deg. + {theta}{sub cone} where {theta}{submore » cone} is the cone angle) is performed with the imaging panel positioned in any location along the beam path. The panoramic projection images of all views for the same gantry angle are then stitched together with the direct image stitching method (i.e., according to the reported imaging position) and full-fan, half-scan CBCT reconstruction is performed using the stitched projection images. To validate this imaging technique, the authors simulated cone-beam projection images of the Mathematical Cardiac Torso (MCAT) thorax phantom for three panoramic views. Gaps, repeated/missing columns, and different exposure levels were introduced between adjacent views to simulate imperfect image stitching due to uncertainties in imaging position or output fluctuation. A modified simultaneous algebraic reconstruction technique (modified SART) was developed to reconstruct CBCT images directly from the stitched projection images. As a gold standard, full-fan, full-scan (360 deg. gantry rotation) CBCT reconstructions were also performed using projection images of one imaging panel large enough to encompass the target. Contrast-to-noise ratio (CNR) and geometric distortion were evaluated to quantify the quality of reconstructed images. Monte Carlo simulations were performed to evaluate the effect of scattering on the image quality and imaging dose for both standard and panoramic CBCT. Results: Truncated images with artifacts were observed for the CBCT reconstruction using projection images of the central view only. When the image stitching was perfect, complete reconstruction was obtained for the panoramic CBCT using the modified SART with the image quality similar to the gold standard (full-scan, full-fan CBCT using one large imaging panel). Imperfect image stitching, on the other hand, lead to (streak, line, or ring) reconstruction artifacts, reduced CNR, and/or distorted geometry. Results from Monte Carlo simulations showed that, for identical imaging quality, the imaging dose was lower for the panoramic CBCT than that acquired with one large imaging panel. For the same imaging dose, the CNR of the three-view panoramic CBCT was 50% higher than that of the regular CBCT using one big panel. Conclusions: The authors have developed a panoramic CBCT technique and demonstrated with simulation data that it can image tumors of any location for patients of any size at the treatment position with comparable or less imaging dose and time. However, the image quality of this CBCT technique is sensitive to the reconstruction artifacts caused by imperfect image stitching. Better algorithms are therefore needed to improve the accuracy of image stitching for panoramic CBCT.« less

Regional projections of North Indian climate for adaptation studies.

PubMed

Mathison, Camilla; Wiltshire, Andrew; Dimri, A P; Falloon, Pete; Jacob, Daniela; Kumar, Pankaj; Moors, Eddy; Ridley, Jeff; Siderius, Christian; Stoffel, Markus; Yasunari, T

2013-12-01

Adaptation is increasingly important for regions around the world where large changes in climate could have an impact on populations and industry. The Brahmaputra-Ganges catchments have a large population, a main industry of agriculture and a growing hydro-power industry, making the region susceptible to changes in the Indian Summer Monsoon, annually the main water source. The HighNoon project has completed four regional climate model simulations for India and the Himalaya at high resolution (25km) from 1960 to 2100 to provide an ensemble of simulations for the region. In this paper we have assessed the ensemble for these catchments, comparing the simulations with observations, to give credence that the simulations provide a realistic representation of atmospheric processes and therefore future climate. We have illustrated how these simulations could be used to provide information on potential future climate impacts and therefore aid decision-making using climatology and threshold analysis. The ensemble analysis shows an increase in temperature between the baseline (1970-2000) and the 2050s (2040-2070) of between 2 and 4°C and an increase in the number of days with maximum temperatures above 28°C and 35°C. There is less certainty for precipitation and runoff which show considerable variability, even in this relatively small ensemble, spanning zero. The HighNoon ensemble is the most complete data for the region providing useful information on a wide range of variables for the regional climate of the Brahmaputra-Ganges region, however there are processes not yet included in the models that could have an impact on the simulations of future climate. We have discussed these processes and show that the range from the HighNoon ensemble is similar in magnitude to potential changes in projections where these processes are included. Therefore strategies for adaptation must be robust and flexible allowing for advances in the science and natural environmental changes. Copyright © 2012 Elsevier B.V. All rights reserved.

Changes in Seasonal and Extreme Hydrologic Conditions of the Georgia Basin/Puget Sound in an Ensemble Regional Climate Simulation for the Mid-Century

DOE Office of Scientific and Technical Information (OSTI.GOV)

Leung, Lai R.; Qian, Yun

This study examines an ensemble of climate change projections simulated by a global climate model (GCM) and downscaled with a region climate model (RCM) to 40 km spatial resolution for the western North America. One control and three ensemble future climate simulations were produced by the GCM following a business as usual scenario for greenhouse gases and aerosols emissions from 1995 to 2100. The RCM was used to downscale the GCM control simulation (1995-2015) and each ensemble future GCM climate (2040-2060) simulation. Analyses of the regional climate simulations for the Georgia Basin/Puget Sound showed a warming of 1.5-2oC and statisticallymore » insignificant changes in precipitation by the mid-century. Climate change has large impacts on snowpack (about 50% reduction) but relatively smaller impacts on the total runoff for the basin as a whole. However, climate change can strongly affect small watersheds such as those located in the transient snow zone, causing a higher likelihood of winter flooding as a higher percentage of precipitation falls in the form of rain rather than snow, and reduced streamflow in early summer. In addition, there are large changes in the monthly total runoff above the upper 1% threshold (or flood volume) from October through May, and the December flood volume of the future climate is 60% above the maximum monthly flood volume of the control climate. Uncertainty of the climate change projections, as characterized by the spread among the ensemble future climate simulations, is relatively small for the basin mean snowpack and runoff, but increases in smaller watersheds, especially in the transient snow zone, and associated with extreme events. This emphasizes the importance of characterizing uncertainty through ensemble simulations.« less

High Fidelity Simulations of Large-Scale Wireless Networks

DOE Office of Scientific and Technical Information (OSTI.GOV)

Onunkwo, Uzoma; Benz, Zachary

The worldwide proliferation of wireless connected devices continues to accelerate. There are 10s of billions of wireless links across the planet with an additional explosion of new wireless usage anticipated as the Internet of Things develops. Wireless technologies do not only provide convenience for mobile applications, but are also extremely cost-effective to deploy. Thus, this trend towards wireless connectivity will only continue and Sandia must develop the necessary simulation technology to proactively analyze the associated emerging vulnerabilities. Wireless networks are marked by mobility and proximity-based connectivity. The de facto standard for exploratory studies of wireless networks is discrete event simulationsmore » (DES). However, the simulation of large-scale wireless networks is extremely difficult due to prohibitively large turnaround time. A path forward is to expedite simulations with parallel discrete event simulation (PDES) techniques. The mobility and distance-based connectivity associated with wireless simulations, however, typically doom PDES and fail to scale (e.g., OPNET and ns-3 simulators). We propose a PDES-based tool aimed at reducing the communication overhead between processors. The proposed solution will use light-weight processes to dynamically distribute computation workload while mitigating communication overhead associated with synchronizations. This work is vital to the analytics and validation capabilities of simulation and emulation at Sandia. We have years of experience in Sandia’s simulation and emulation projects (e.g., MINIMEGA and FIREWHEEL). Sandia’s current highly-regarded capabilities in large-scale emulations have focused on wired networks, where two assumptions prevent scalable wireless studies: (a) the connections between objects are mostly static and (b) the nodes have fixed locations.« less

Interactive Correlation Analysis and Visualization of Climate Data

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ma, Kwan-Liu

The relationship between our ability to analyze and extract insights from visualization of climate model output and the capability of the available resources to make those visualizations has reached a crisis point. The large volume of data currently produced by climate models is overwhelming the current, decades-old visualization workflow. The traditional methods for visualizing climate output also have not kept pace with changes in the types of grids used, the number of variables involved, and the number of different simulations performed with a climate model or the feature-richness of high-resolution simulations. This project has developed new and faster methods formore » visualization in order to get the most knowledge out of the new generation of high-resolution climate models. While traditional climate images will continue to be useful, there is need for new approaches to visualization and analysis of climate data if we are to gain all the insights available in ultra-large data sets produced by high-resolution model output and ensemble integrations of climate models such as those produced for the Coupled Model Intercomparison Project. Towards that end, we have developed new visualization techniques for performing correlation analysis. We have also introduced highly scalable, parallel rendering methods for visualizing large-scale 3D data. This project was done jointly with climate scientists and visualization researchers at Argonne National Laboratory and NCAR.« less

Evaluation of regional climate simulations for air quality modelling purposes

NASA Astrophysics Data System (ADS)

Menut, Laurent; Tripathi, Om P.; Colette, Augustin; Vautard, Robert; Flaounas, Emmanouil; Bessagnet, Bertrand

2013-05-01

In order to evaluate the future potential benefits of emission regulation on regional air quality, while taking into account the effects of climate change, off-line air quality projection simulations are driven using weather forcing taken from regional climate models. These regional models are themselves driven by simulations carried out using global climate models (GCM) and economical scenarios. Uncertainties and biases in climate models introduce an additional "climate modeling" source of uncertainty that is to be added to all other types of uncertainties in air quality modeling for policy evaluation. In this article we evaluate the changes in air quality-related weather variables induced by replacing reanalyses-forced by GCM-forced regional climate simulations. As an example we use GCM simulations carried out in the framework of the ERA-interim programme and of the CMIP5 project using the Institut Pierre-Simon Laplace climate model (IPSLcm), driving regional simulations performed in the framework of the EURO-CORDEX programme. In summer, we found compensating deficiencies acting on photochemistry: an overestimation by GCM-driven weather due to a positive bias in short-wave radiation, a negative bias in wind speed, too many stagnant episodes, and a negative temperature bias. In winter, air quality is mostly driven by dispersion, and we could not identify significant differences in either wind or planetary boundary layer height statistics between GCM-driven and reanalyses-driven regional simulations. However, precipitation appears largely overestimated in GCM-driven simulations, which could significantly affect the simulation of aerosol concentrations. The identification of these biases will help interpreting results of future air quality simulations using these data. Despite these, we conclude that the identified differences should not lead to major difficulties in using GCM-driven regional climate simulations for air quality projections.

Large projected increases in rain-on-snow flood potential over western North America

NASA Astrophysics Data System (ADS)

Musselman, K. N.; Ikeda, K.; Barlage, M. J.; Lehner, F.; Liu, C.; Newman, A. J.; Prein, A. F.; Mizukami, N.; Gutmann, E. D.; Clark, M. P.; Rasmussen, R.

2017-12-01

In the western US and Canada, some of the largest annual flood events occur when warm storm systems drop substantial rainfall on extensive snow-cover. For example, last winter's Oroville dam crisis in California was exacerbated by rapid snowmelt during a rain-on-snow (ROS) event. We present an analysis of ROS events with flood-generating potential over western North America simulated at high-resolution by the Weather Research and Forecasting (WRF) model run for both a 13-year control time period and re-run with a `business-as-usual' future (2071-2100) climate scenario. Daily ROS with flood-generating potential is defined as rainfall of at least 10 mm per day falling on snowpack of at least 10 mm water equivalent, where the sum of rainfall and snowmelt contains at least 20% snowmelt. In a warmer climate, ROS is less frequent in regions where it is historically common, and more frequent elsewhere. This is evidenced by large simulated reductions in snow-cover and ROS frequency at lower elevations, particularly in warmer, coastal regions, and greater ROS frequency at middle elevations and in inland regions. The same trend is reflected in the annual-average ROS runoff volume (rainfall + snowmelt) aggregated to major watersheds; large reductions of 25-75% are projected for much of the U.S. Pacific Northwest, while large increases are simulated for the Colorado River basin, western Canada, and the higher elevations of the Sierra Nevada. In the warmer climate, snowmelt contributes substantially less to ROS runoff per unit rainfall, particularly in inland regions. The reduction in snowmelt contribution is due to a shift in ROS timing from warm spring events to cooler winter conditions and/or from warm, lower elevations to cool, higher elevations. However, the slower snowmelt is offset by an increase in rainfall intensity, maintaining the flood potential of ROS at or above historical levels. In fact, we report large projected increases in the intensity of extreme ROS events. The projected increases in ROS flood potential are highest in historically flood-prone mountain basins and the Canadian Prairies. Increases in extreme ROS event intensity, together with a greater proportion of precipitation falling as rain, have critical implications on the climate resilience of regional flood control systems.

Estimation of the mid-century Etesians wind pattern from EURO-CORDEX models

NASA Astrophysics Data System (ADS)

Dafka, Stella; Toreti, Andrea; Luterbacher, Juerg; Zanis, Prodromos; Tyrlis, Evangelos; Xoplaki, Elena

2017-04-01

The Etesians are one of the major and most prominent wind system, prevailing over the Aegean Sea during summer and early autumn. Here, projections of changes in 30-year (2021-2050) wind speeds relative to 1971-2000, under the 8.5 and 4.5 Representative Concentration Pathways, have been produced for Etesians. Future changes in the number of Etesian days and the associated large scale dynamics are also considered. We analyze seven simulations from three EURO-CORDEX regional climate models at a 12 km grid resolution. Both scenarios indicate that in most RCMs daily wind speeds are projected to increase by 1-1.5m/s over the Aegean Sea, suggesting that the current estimate of wind power potential for Aegean Sea will be increased with the greenhouse gas forcing in the coming decades (2021-2050). Wind direction at 10-m as well as the number of Etesian days have shown to undergo minor changes. The projected changes in sea level pressure and geopotential height anomalies at 500 hPa have a large spread among the seven simulations with a disperse tendency of strengthening of the ridge over the Balkans.

Uncertain soil moisture feedbacks in model projections of Sahel precipitation

NASA Astrophysics Data System (ADS)

Berg, Alexis; Lintner, Benjamin R.; Findell, Kirsten; Giannini, Alessandra

2017-06-01

Given the uncertainties in climate model projections of Sahel precipitation, at the northern edge of the West African Monsoon, understanding the factors governing projected precipitation changes in this semiarid region is crucial. This study investigates how long-term soil moisture changes projected under climate change may feedback on projected changes of Sahel rainfall, using simulations with and without soil moisture change from five climate models participating in the Global Land Atmosphere Coupling Experiment-Coupled Model Intercomparison Project phase 5 experiment. In four out of five models analyzed, soil moisture feedbacks significantly influence the projected West African precipitation response to warming; however, the sign of these feedbacks differs across the models. These results demonstrate that reducing uncertainties across model projections of the West African Monsoon requires, among other factors, improved mechanistic understanding and constraint of simulated land-atmosphere feedbacks, even at the large spatial scales considered here.