Science.gov

Sample records for general event-driven simulator

  1. Event-driven simulation in SELMON: An overview of EDSE

    NASA Technical Reports Server (NTRS)

    Rouquette, Nicolas F.; Chien, Steve A.; Charest, Leonard, Jr.

    1992-01-01

    EDSE (event-driven simulation engine), a model-based event-driven simulator implemented for SELMON, a tool for sensor selection and anomaly detection in real-time monitoring is described. The simulator is used in conjunction with a causal model to predict future behavior of the model from observed data. The behavior of the causal model is interpreted as equivalent to the behavior of the physical system being modeled. An overview of the functionality of the simulator and the model-based event-driven simulation paradigm on which it is based is provided. Included are high-level descriptions of the following key properties: event consumption and event creation, iterative simulation, synchronization and filtering of monitoring data from the physical system. Finally, how EDSE stands with respect to the relevant open issues of discrete-event and model-based simulation is discussed.

  2. Cellular Dynamic Simulator: An Event Driven Molecular Simulation Environment for Cellular Physiology

    PubMed Central

    Byrne, Michael J.; Waxham, M. Neal; Kubota, Yoshihisa

    2010-01-01

    In this paper, we present the Cellular Dynamic Simulator (CDS) for simulating diffusion and chemical reactions within crowded molecular environments. CDS is based on a novel event driven algorithm specifically designed for precise calculation of the timing of collisions, reactions and other events for each individual molecule in the environment. Generic mesh based compartments allow the creation / importation of very simple or detailed cellular structures that exist in a 3D environment. Multiple levels of compartments and static obstacles can be used to create a dense environment to mimic cellular boundaries and the intracellular space. The CDS algorithm takes into account volume exclusion and molecular crowding that may impact signaling cascades in small sub-cellular compartments such as dendritic spines. With the CDS, we can simulate simple enzyme reactions; aggregation, channel transport, as well as highly complicated chemical reaction networks of both freely diffusing and membrane bound multi-protein complexes. Components of the CDS are generally defined such that the simulator can be applied to a wide range of environments in terms of scale and level of detail. Through an initialization GUI, a simple simulation environment can be created and populated within minutes yet is powerful enough to design complex 3D cellular architecture. The initialization tool allows visual confirmation of the environment construction prior to execution by the simulator. This paper describes the CDS algorithm, design implementation, and provides an overview of the types of features available and the utility of those features are highlighted in demonstrations. PMID:20361275

  3. Simulating large-scale pedestrian movement using CA and event driven model: Methodology and case study

    NASA Astrophysics Data System (ADS)

    Li, Jun; Fu, Siyao; He, Haibo; Jia, Hongfei; Li, Yanzhong; Guo, Yi

    2015-11-01

    Large-scale regional evacuation is an important part of national security emergency response plan. Large commercial shopping area, as the typical service system, its emergency evacuation is one of the hot research topics. A systematic methodology based on Cellular Automata with the Dynamic Floor Field and event driven model has been proposed, and the methodology has been examined within context of a case study involving the evacuation within a commercial shopping mall. Pedestrians walking is based on Cellular Automata and event driven model. In this paper, the event driven model is adopted to simulate the pedestrian movement patterns, the simulation process is divided into normal situation and emergency evacuation. The model is composed of four layers: environment layer, customer layer, clerk layer and trajectory layer. For the simulation of movement route of pedestrians, the model takes into account purchase intention of customers and density of pedestrians. Based on evacuation model of Cellular Automata with Dynamic Floor Field and event driven model, we can reflect behavior characteristics of customers and clerks at the situations of normal and emergency evacuation. The distribution of individual evacuation time as a function of initial positions and the dynamics of the evacuation process is studied. Our results indicate that the evacuation model using the combination of Cellular Automata with Dynamic Floor Field and event driven scheduling can be used to simulate the evacuation of pedestrian flows in indoor areas with complicated surroundings and to investigate the layout of shopping mall.

  4. An Event-Driven Hybrid Molecular Dynamics and Direct Simulation Monte Carlo Algorithm

    SciTech Connect

    Donev, A; Garcia, A L; Alder, B J

    2007-07-30

    A novel algorithm is developed for the simulation of polymer chains suspended in a solvent. The polymers are represented as chains of hard spheres tethered by square wells and interact with the solvent particles with hard core potentials. The algorithm uses event-driven molecular dynamics (MD) for the simulation of the polymer chain and the interactions between the chain beads and the surrounding solvent particles. The interactions between the solvent particles themselves are not treated deterministically as in event-driven algorithms, rather, the momentum and energy exchange in the solvent is determined stochastically using the Direct Simulation Monte Carlo (DSMC) method. The coupling between the solvent and the solute is consistently represented at the particle level, however, unlike full MD simulations of both the solvent and the solute, the spatial structure of the solvent is ignored. The algorithm is described in detail and applied to the study of the dynamics of a polymer chain tethered to a hard wall subjected to uniform shear. The algorithm closely reproduces full MD simulations with two orders of magnitude greater efficiency. Results do not confirm the existence of periodic (cycling) motion of the polymer chain.

  5. NEVESIM: event-driven neural simulation framework with a Python interface

    PubMed Central

    Pecevski, Dejan; Kappel, David; Jonke, Zeno

    2014-01-01

    NEVESIM is a software package for event-driven simulation of networks of spiking neurons with a fast simulation core in C++, and a scripting user interface in the Python programming language. It supports simulation of heterogeneous networks with different types of neurons and synapses, and can be easily extended by the user with new neuron and synapse types. To enable heterogeneous networks and extensibility, NEVESIM is designed to decouple the simulation logic of communicating events (spikes) between the neurons at a network level from the implementation of the internal dynamics of individual neurons. In this paper we will present the simulation framework of NEVESIM, its concepts and features, as well as some aspects of the object-oriented design approaches and simulation strategies that were utilized to efficiently implement the concepts and functionalities of the framework. We will also give an overview of the Python user interface, its basic commands and constructs, and also discuss the benefits of integrating NEVESIM with Python. One of the valuable capabilities of the simulator is to simulate exactly and efficiently networks of stochastic spiking neurons from the recently developed theoretical framework of neural sampling. This functionality was implemented as an extension on top of the basic NEVESIM framework. Altogether, the intended purpose of the NEVESIM framework is to provide a basis for further extensions that support simulation of various neural network models incorporating different neuron and synapse types that can potentially also use different simulation strategies. PMID:25177291

  6. NEVESIM: event-driven neural simulation framework with a Python interface.

    PubMed

    Pecevski, Dejan; Kappel, David; Jonke, Zeno

    2014-01-01

    NEVESIM is a software package for event-driven simulation of networks of spiking neurons with a fast simulation core in C++, and a scripting user interface in the Python programming language. It supports simulation of heterogeneous networks with different types of neurons and synapses, and can be easily extended by the user with new neuron and synapse types. To enable heterogeneous networks and extensibility, NEVESIM is designed to decouple the simulation logic of communicating events (spikes) between the neurons at a network level from the implementation of the internal dynamics of individual neurons. In this paper we will present the simulation framework of NEVESIM, its concepts and features, as well as some aspects of the object-oriented design approaches and simulation strategies that were utilized to efficiently implement the concepts and functionalities of the framework. We will also give an overview of the Python user interface, its basic commands and constructs, and also discuss the benefits of integrating NEVESIM with Python. One of the valuable capabilities of the simulator is to simulate exactly and efficiently networks of stochastic spiking neurons from the recently developed theoretical framework of neural sampling. This functionality was implemented as an extension on top of the basic NEVESIM framework. Altogether, the intended purpose of the NEVESIM framework is to provide a basis for further extensions that support simulation of various neural network models incorporating different neuron and synapse types that can potentially also use different simulation strategies. PMID:25177291

  7. Self-Adaptive Event-Driven Simulation of Multi-Scale Plasma Systems

    NASA Astrophysics Data System (ADS)

    Omelchenko, Yuri; Karimabadi, Homayoun

    2005-10-01

    Multi-scale plasmas pose a formidable computational challenge. The explicit time-stepping models suffer from the global CFL restriction. Efficient application of adaptive mesh refinement (AMR) to systems with irregular dynamics (e.g. turbulence, diffusion-convection-reaction, particle acceleration etc.) may be problematic. To address these issues, we developed an alternative approach to time stepping: self-adaptive discrete-event simulation (DES). DES has origin in operations research, war games and telecommunications. We combine finite-difference and particle-in-cell techniques with this methodology by assuming two caveats: (1) a local time increment, dt for a discrete quantity f can be expressed in terms of a physically meaningful quantum value, df; (2) f is considered to be modified only when its change exceeds df. Event-driven time integration is self-adaptive as it makes use of causality rules rather than parametric time dependencies. This technique enables asynchronous flux-conservative update of solution in accordance with local temporal scales, removes the curse of the global CFL condition, eliminates unnecessary computation in inactive spatial regions and results in robust and fast parallelizable codes. It can be naturally combined with various mesh refinement techniques. We discuss applications of this novel technology to diffusion-convection-reaction systems and hybrid simulations of magnetosonic shocks.

  8. Efficient event-driven simulations shed new light on microtubule organization in the plant cortical array

    NASA Astrophysics Data System (ADS)

    Tindemans, Simon H.; Deinum, Eva E.; Lindeboom, Jelmer J.; Mulder, Bela M.

    2014-04-01

    The dynamics of the plant microtubule cytoskeleton is a paradigmatic example of the complex spatiotemporal processes characterising life at the cellular scale. This system is composed of large numbers of spatially extended particles, each endowed with its own intrinsic stochastic dynamics, and is capable of non-equilibrium self-organisation through collisional interactions of these particles. To elucidate the behaviour of such a complex system requires not only conceptual advances, but also the development of appropriate computational tools to simulate it. As the number of parameters involved is large and the behaviour is stochastic, it is essential that these simulations be fast enough to allow for an exploration of the phase space and the gathering of sufficient statistics to accurately pin down the average behaviour as well as the magnitude of fluctuations around it. Here we describe a simulation approach that meets this requirement by adopting an event-driven methodology that encompasses both the spontaneous stochastic changes in microtubule state as well as the deterministic collisions. In contrast with finite time step simulations this technique is intrinsically exact, as well as several orders of magnitude faster, which enables ordinary PC hardware to simulate systems of ˜ 10^3 microtubules on a time scale ˜ 10^{3} faster than real time. In addition we present new tools for the analysis of microtubule trajectories on curved surfaces. We illustrate the use of these methods by addressing a number of outstanding issues regarding the importance of various parameters on the transition from an isotropic to an aligned and oriented state.

  9. A new concept for simulation of vegetated land surface dynamics - Part 1: The event driven phenology model

    NASA Astrophysics Data System (ADS)

    Kovalskyy, V.; Henebry, G. M.

    2012-01-01

    Phenologies of the vegetated land surface are being used increasingly for diagnosis and prognosis of climate change consequences. Current prospective and retrospective phenological models stand far apart in their approaches to the subject. We report on an exploratory attempt to implement a phenological model based on a new event driven concept which has both diagnostic and prognostic capabilities in the same modeling framework. This Event Driven Phenological Model (EDPM) is shown to simulate land surface phenologies and phenophase transition dates in agricultural landscapes based on assimilation of weather data and land surface observations from spaceborne sensors. The model enables growing season phenologies to develop in response to changing environmental conditions and disturbance events. It also has the ability to ingest remotely sensed data to adjust its output to improve representation of the modeled variable. We describe the model and report results of initial testing of the EDPM using Level 2 flux tower records from the Ameriflux sites at Mead, Nebraska, USA, and at Bondville, Illinois, USA. Simulating the dynamics of normalized difference vegetation index based on flux tower data, the predictions by the EDPM show good agreement (RMSE < 0.08; r2 > 0.8) for maize and soybean during several growing seasons at different locations. This study presents the EDPM used in the companion paper (Kovalskyy and Henebry, 2011) in a coupling scheme to estimate daily actual evapotranspiration over multiple growing seasons.

  10. A new concept for simulation of vegetated land surface dynamics - Part 1: The event driven phenology model

    NASA Astrophysics Data System (ADS)

    Kovalskyy, V.; Henebry, G. M.

    2011-05-01

    Phenologies of the vegetated land surface are being used increasingly for diagnosis and prognosis of climate change consequences. Current prospective and retrospective phenological models stand far apart in their approaches to the subject. We report on an exploratory attempt to implement a phenological model based on a new event driven concept which has both diagnostic and prognostic capabilities in the same modeling framework. This Event Driven Phenological Model (EDPM) is shown to simulate land surface phenologies and phenophase transition dates in agricultural landscapes based on assimilation of weather data and land surface observations from spaceborne sensors. The model enables growing season phenologies to develop in response to changing environmental conditions and disturbance events. It also has the ability to ingest remotely sensed data to adjust its output to improve representation of the modeled variable. We describe the model and report results of initial testing of the EDPM using Level 2 flux tower records from the Ameriflux sites at Mead, Nebraska, USA, and at Bondville, Illinois, USA. Simulating the dynamics of normalized difference vegetation index based on flux tower data, the predictions by the EDPM show good agreement (RMSE < 0.08; r2>0.8) for maize and soybean during several growing seasons at different locations. This study presents the EDPM used in the companion paper (Kovalskyy and Henebry, 2011) in a coupling scheme to estimate daily actual evapotranspiration over multiple growing seasons.

  11. A combined Event-Driven/Time-Driven molecular dynamics algorithm for the simulation of shock waves in rarefied gases

    SciTech Connect

    Valentini, Paolo Schwartzentruber, Thomas E.

    2009-12-10

    A novel combined Event-Driven/Time-Driven (ED/TD) algorithm to speed-up the Molecular Dynamics simulation of rarefied gases using realistic spherically symmetric soft potentials is presented. Due to the low density regime, the proposed method correctly identifies the time that must elapse before the next interaction occurs, similarly to Event-Driven Molecular Dynamics. However, each interaction is treated using Time-Driven Molecular Dynamics, thereby integrating Newton's Second Law using the sufficiently small time step needed to correctly resolve the atomic motion. Although infrequent, many-body interactions are also accounted for with a small approximation. The combined ED/TD method is shown to correctly reproduce translational relaxation in argon, described using the Lennard-Jones potential. For densities between {rho}=10{sup -4}kg/m{sup 3} and {rho}=10{sup -1}kg/m{sup 3}, comparisons with kinetic theory, Direct Simulation Monte Carlo, and pure Time-Driven Molecular Dynamics demonstrate that the ED/TD algorithm correctly reproduces the proper collision rates and the evolution toward thermal equilibrium. Finally, the combined ED/TD algorithm is applied to the simulation of a Mach 9 shock wave in rarefied argon. Density and temperature profiles as well as molecular velocity distributions accurately match DSMC results, and the shock thickness is within the experimental uncertainty. For the problems considered, the ED/TD algorithm ranged from several hundred to several thousand times faster than conventional Time-Driven MD. Moreover, the force calculation to integrate the molecular trajectories is found to contribute a negligible amount to the overall ED/TD simulation time. Therefore, this method could pave the way for the application of much more refined and expensive interatomic potentials, either classical or first-principles, to Molecular Dynamics simulations of shock waves in rarefied gases, involving vibrational nonequilibrium and chemical reactivity.

  12. Asynchronous Event-Driven Particle Algorithms

    SciTech Connect

    Donev, A

    2007-02-28

    We present in a unifying way the main components of three examples of asynchronous event-driven algorithms for simulating physical systems of interacting particles. The first example, hard-particle molecular dynamics (MD), is well-known. We also present a recently-developed diffusion kinetic Monte Carlo (DKMC) algorithm, as well as a novel event-driven algorithm for Direct Simulation Monte Carlo (DSMC). Finally, we describe how to combine MD with DSMC in an event-driven framework, and discuss some promises and challenges for event-driven simulation of realistic physical systems.

  13. Anticipating the Chaotic Behaviour of Industrial Systems Based on Stochastic, Event-Driven Simulations

    NASA Astrophysics Data System (ADS)

    Bruzzone, Agostino G.; Revetria, Roberto; Simeoni, Simone; Viazzo, Simone; Orsoni, Alessandra

    2004-08-01

    In logistics and industrial production managers must deal with the impact of stochastic events to improve performances and reduce costs. In fact, production and logistics systems are generally designed considering some parameters as deterministically distributed. While this assumption is mostly used for preliminary prototyping, it is sometimes also retained during the final design stage, and especially for estimated parameters (i.e. Market Request). The proposed methodology can determine the impact of stochastic events in the system by evaluating the chaotic threshold level. Such an approach, based on the application of a new and innovative methodology, can be implemented to find the condition under which chaos makes the system become uncontrollable. Starting from problem identification and risk assessment, several classification techniques are used to carry out an effect analysis and contingency plan estimation. In this paper the authors illustrate the methodology with respect to a real industrial case: a production problem related to the logistics of distributed chemical processing.

  14. Asynchronous Event-Driven Particle Algorithms

    SciTech Connect

    Donev, A

    2007-08-30

    We present, in a unifying way, the main components of three asynchronous event-driven algorithms for simulating physical systems of interacting particles. The first example, hard-particle molecular dynamics (MD), is well-known. We also present a recently-developed diffusion kinetic Monte Carlo (DKMC) algorithm, as well as a novel stochastic molecular-dynamics algorithm that builds on the Direct Simulation Monte Carlo (DSMC). We explain how to effectively combine event-driven and classical time-driven handling, and discuss some promises and challenges for event-driven simulation of realistic physical systems.

  15. Stochastic Event-Driven Molecular Dynamics

    SciTech Connect

    Donev, Aleksandar Garcia, Alejandro L.; Alder, Berni J.

    2008-02-01

    A novel Stochastic Event-Driven Molecular Dynamics (SEDMD) algorithm is developed for the simulation of polymer chains suspended in a solvent. SEDMD combines event-driven molecular dynamics (EDMD) with the Direct Simulation Monte Carlo (DSMC) method. The polymers are represented as chains of hard-spheres tethered by square wells and interact with the solvent particles with hard-core potentials. The algorithm uses EDMD for the simulation of the polymer chain and the interactions between the chain beads and the surrounding solvent particles. The interactions between the solvent particles themselves are not treated deterministically as in EDMD, rather, the momentum and energy exchange in the solvent is determined stochastically using DSMC. The coupling between the solvent and the solute is consistently represented at the particle level retaining hydrodynamic interactions and thermodynamic fluctuations. However, unlike full MD simulations of both the solvent and the solute, in SEDMD the spatial structure of the solvent is ignored. The SEDMD algorithm is described in detail and applied to the study of the dynamics of a polymer chain tethered to a hard-wall subjected to uniform shear. SEDMD closely reproduces results obtained using traditional EDMD simulations with two orders of magnitude greater efficiency. Results question the existence of periodic (cycling) motion of the polymer chain.

  16. A General and Efficient Method for Incorporating Precise Spike Times in Globally Time-Driven Simulations

    PubMed Central

    Hanuschkin, Alexander; Kunkel, Susanne; Helias, Moritz; Morrison, Abigail; Diesmann, Markus

    2010-01-01

    Traditionally, event-driven simulations have been limited to the very restricted class of neuronal models for which the timing of future spikes can be expressed in closed form. Recently, the class of models that is amenable to event-driven simulation has been extended by the development of techniques to accurately calculate firing times for some integrate-and-fire neuron models that do not enable the prediction of future spikes in closed form. The motivation of this development is the general perception that time-driven simulations are imprecise. Here, we demonstrate that a globally time-driven scheme can calculate firing times that cannot be discriminated from those calculated by an event-driven implementation of the same model; moreover, the time-driven scheme incurs lower computational costs. The key insight is that time-driven methods are based on identifying a threshold crossing in the recent past, which can be implemented by a much simpler algorithm than the techniques for predicting future threshold crossings that are necessary for event-driven approaches. As run time is dominated by the cost of the operations performed at each incoming spike, which includes spike prediction in the case of event-driven simulation and retrospective detection in the case of time-driven simulation, the simple time-driven algorithm outperforms the event-driven approaches. Additionally, our method is generally applicable to all commonly used integrate-and-fire neuronal models; we show that a non-linear model employing a standard adaptive solver can reproduce a reference spike train with a high degree of precision. PMID:21031031

  17. Event-Driven Process Chains (EPC)

    NASA Astrophysics Data System (ADS)

    Mendling, Jan

    This chapter provides a comprehensive overview of Event-driven Process Chains (EPCs) and introduces a novel definition of EPC semantics. EPCs became popular in the 1990s as a conceptual business process modeling language in the context of reference modeling. Reference modeling refers to the documentation of generic business operations in a model such as service processes in the telecommunications sector, for example. It is claimed that reference models can be reused and adapted as best-practice recommendations in individual companies (see [230, 168, 229, 131, 400, 401, 446, 127, 362, 126]). The roots of reference modeling can be traced back to the Kölner Integrationsmodell (KIM) [146, 147] that was developed in the 1960s and 1970s. In the 1990s, the Institute of Information Systems (IWi) in Saarbrücken worked on a project with SAP to define a suitable business process modeling language to document the processes of the SAP R/3 enterprise resource planning system. There were two results from this joint effort: the definition of EPCs [210] and the documentation of the SAP system in the SAP Reference Model (see [92, 211]). The extensive database of this reference model contains almost 10,000 sub-models: 604 of them non-trivial EPC business process models. The SAP Reference model had a huge impact with several researchers referring to it in their publications (see [473, 235, 127, 362, 281, 427, 415]) as well as motivating the creation of EPC reference models in further domains including computer integrated manufacturing [377, 379], logistics [229] or retail [52]. The wide-spread application of EPCs in business process modeling theory and practice is supported by their coverage in seminal text books for business process management and information systems in general (see [378, 380, 49, 384, 167, 240]). EPCs are frequently used in practice due to a high user acceptance [376] and extensive tool support. Some examples of tools that support EPCs are ARIS Toolset by IDS

  18. The three-dimensional Event-Driven Graphics Environment (3D-EDGE)

    NASA Technical Reports Server (NTRS)

    Freedman, Jeffrey; Hahn, Roger; Schwartz, David M.

    1993-01-01

    Stanford Telecom developed the Three-Dimensional Event-Driven Graphics Environment (3D-EDGE) for NASA GSFC's (GSFC) Communications Link Analysis and Simulation System (CLASS). 3D-EDGE consists of a library of object-oriented subroutines which allow engineers with little or no computer graphics experience to programmatically manipulate, render, animate, and access complex three-dimensional objects.

  19. Feasibility study for a generalized gate logic software simulator

    NASA Technical Reports Server (NTRS)

    Mcgough, J. G.

    1983-01-01

    Unit-delay simulation, event driven simulation, zero-delay simulation, simulation techniques, 2-valued versus multivalued logic, network initialization, gate operations and alternate network representations, parallel versus serial mode simulation fault modelling, extension of multiprocessor systems, and simulation timing are discussed. Functional level networks, gate equivalent circuits, the prototype BDX-930 network model, fault models, identifying detected faults for BGLOSS are discussed. Preprocessor tasks, postprocessor tasks, executive tasks, and a library of bliss coded macros for GGLOSS are also discussed.

  20. An Event Driven Phenology Model: Results from initial testing

    NASA Astrophysics Data System (ADS)

    Kovalskyy, V.; Henebry, G.

    2008-12-01

    A new model was developed to simulate land surface phenology and seasonality dynamics based on assimilation of weather data and operational land surface observations. The model was built to establish a much needed interface between boundary layer dynamics represented in weather models and surface observations and the temporal variability captured in series of remotely sensed images. The Event Driven Phenology Model (EDPM) takes advantage of both ground-based weather and climate records and spaceborne sensors to provide retrospective and predictive power. The new model enables phenologies to unfold in response to changing environmental conditions and disturbance events. It also has the ability to ingest contemporaneous discrete external records of land surface dynamics to adjust its output to achieve a better representation of the observed process. The EDPM presents an alternative to contemporary methods such as retrospective curve-fitting (either to time or to a temporal proxy, such as thermal time), long term averages (or climatologies), and phenologies from look-up tables based on land use/land cover or plant functional types. We describe the model and report results of initial testing of the EDPM using level 1 flux tower records from the Mead, Nebraska Ameriflux sites in conjunction with associated MODIS subsets from the Carbon DAAC at ORNL. We assess the EDPM predictions by comparing and contrasting the results with reserved ground records and with outcomes of other phenology models. Finally, we point out prospects for future use of descriptive and prescriptive EDPM capabilities in the work of climate models, production of continuous remote sensing records, and other scientific applications.

  1. Two-ball problem revisited: Limitations of event-driven modeling

    NASA Astrophysics Data System (ADS)

    Müller, Patric; Pöschel, Thorsten

    2011-04-01

    The main precondition of simulating systems of hard particles by means of event-driven modeling is the assumption of instantaneous collisions. The aim of this paper is to quantify the deviation of event-driven modeling from the solution of Newton’s equation of motion using a paradigmatic example: If a tennis ball is held above a basketball with their centers vertically aligned, and the balls are released to collide with the floor, the tennis ball may rebound at a surprisingly high speed. We show in this article that the simple textbook explanation of this effect is an oversimplification, even for the limit of perfectly elastic particles. Instead, there may occur a rather complex scenario including multiple collisions which may lead to a very different final velocity as compared with the velocity resulting from the oversimplified model.

  2. General Data Simulation Program.

    ERIC Educational Resources Information Center

    Burns, Edward

    Described is a computer program written in FORTRAN IV which offers considerable flexibility in generating simulated data pertinent to education and educational psychology. The user is allowed to specify the number of samples, data sets, and variables, together with the population means, standard deviations and intercorrelations. In addition the…

  3. Event management for large scale event-driven digital hardware spiking neural networks.

    PubMed

    Caron, Louis-Charles; D'Haene, Michiel; Mailhot, Frédéric; Schrauwen, Benjamin; Rouat, Jean

    2013-09-01

    The interest in brain-like computation has led to the design of a plethora of innovative neuromorphic systems. Individually, spiking neural networks (SNNs), event-driven simulation and digital hardware neuromorphic systems get a lot of attention. Despite the popularity of event-driven SNNs in software, very few digital hardware architectures are found. This is because existing hardware solutions for event management scale badly with the number of events. This paper introduces the structured heap queue, a pipelined digital hardware data structure, and demonstrates its suitability for event management. The structured heap queue scales gracefully with the number of events, allowing the efficient implementation of large scale digital hardware event-driven SNNs. The scaling is linear for memory, logarithmic for logic resources and constant for processing time. The use of the structured heap queue is demonstrated on a field-programmable gate array (FPGA) with an image segmentation experiment and a SNN of 65,536 neurons and 513,184 synapses. Events can be processed at the rate of 1 every 7 clock cycles and a 406×158 pixel image is segmented in 200 ms. PMID:23522624

  4. A tutorial on creating logfiles for event-driven applications.

    PubMed

    Breinholt, G; Krueger, H

    1999-08-01

    This paper describes the practical steps necessary to write logfiles for recording user actions in event-driven applications. Data logging has long been used as a reliable method to record all user actions, whether assessing new software or running a behavioral experiment. With the widespread introduction of event-driven software, the logfile must enable accurate recording of all the user's actions, whether with the keyboard or another input device. Logging is only an effective tool when it can accurately and consistently record all actions in a format that aids the extraction of useful information from the mass of data collected. Logfiles are often presented as one of many methods that could be used, and here a technique is proposed for the construction of logfiles for the quantitative assessment of software from the user's point of view. PMID:10502862

  5. Event-Driven Cyberinfrastructure Technologies Supporting the Disaster Life Cycle

    NASA Astrophysics Data System (ADS)

    Graves, S. J.; Maskey, M.; Keiser, K.

    2014-12-01

    The cyberinfrastructure components to be discussed include Event-Driven Data Delivery (ED3) and Data Albums. These are complementary technologies that combine to provide comprehensive access to timely and relevant data for disaster events. ED3 provides a cyber framework that allows situational awareness and decision systems to prepare data plans consisting of data access, generation, workflows, etc., that meet the users' data needs in the event of a future disaster event. Data Albums provides a resulting container of relevant data and functionality for an overall information package for a specific event. The combination of these technologies provides useful data capabilities as part of the disaster life cycle.

  6. Intelligent fuzzy controller for event-driven real time systems

    NASA Technical Reports Server (NTRS)

    Grantner, Janos; Patyra, Marek; Stachowicz, Marian S.

    1992-01-01

    Most of the known linguistic models are essentially static, that is, time is not a parameter in describing the behavior of the object's model. In this paper we show a model for synchronous finite state machines based on fuzzy logic. Such finite state machines can be used to build both event-driven, time-varying, rule-based systems and the control unit section of a fuzzy logic computer. The architecture of a pipelined intelligent fuzzy controller is presented, and the linguistic model is represented by an overall fuzzy relation stored in a single rule memory. A VLSI integrated circuit implementation of the fuzzy controller is suggested. At a clock rate of 30 MHz, the controller can perform 3 MFLIPS on multi-dimensional fuzzy data.

  7. Multirate and event-driven Kalman filters for helicopter flight

    NASA Technical Reports Server (NTRS)

    Sridhar, Banavar; Smith, Phillip; Suorsa, Raymond E.; Hussien, Bassam

    1993-01-01

    A vision-based obstacle detection system that provides information about objects as a function of azimuth and elevation is discussed. The range map is computed using a sequence of images from a passive sensor, and an extended Kalman filter is used to estimate range to obstacles. The magnitude of the optical flow that provides measurements for each Kalman filter varies significantly over the image depending on the helicopter motion and object location. In a standard Kalman filter, the measurement update takes place at fixed intervals. It may be necessary to use a different measurement update rate in different parts of the image in order to maintain the same signal to noise ratio in the optical flow calculations. A range estimation scheme that accepts the measurement only under certain conditions is presented. The estimation results from the standard Kalman filter are compared with results from a multirate Kalman filter and an event-driven Kalman filter for a sequence of helicopter flight images.

  8. A Full Parallel Event Driven Readout Technique for Area Array SPAD FLIM Image Sensors

    PubMed Central

    Nie, Kaiming; Wang, Xinlei; Qiao, Jun; Xu, Jiangtao

    2016-01-01

    This paper presents a full parallel event driven readout method which is implemented in an area array single-photon avalanche diode (SPAD) image sensor for high-speed fluorescence lifetime imaging microscopy (FLIM). The sensor only records and reads out effective time and position information by adopting full parallel event driven readout method, aiming at reducing the amount of data. The image sensor includes four 8 × 8 pixel arrays. In each array, four time-to-digital converters (TDCs) are used to quantize the time of photons’ arrival, and two address record modules are used to record the column and row information. In this work, Monte Carlo simulations were performed in Matlab in terms of the pile-up effect induced by the readout method. The sensor’s resolution is 16 × 16. The time resolution of TDCs is 97.6 ps and the quantization range is 100 ns. The readout frame rate is 10 Mfps, and the maximum imaging frame rate is 100 fps. The chip’s output bandwidth is 720 MHz with an average power of 15 mW. The lifetime resolvability range is 5–20 ns, and the average error of estimated fluorescence lifetimes is below 1% by employing CMM to estimate lifetimes. PMID:26828490

  9. A Full Parallel Event Driven Readout Technique for Area Array SPAD FLIM Image Sensors.

    PubMed

    Nie, Kaiming; Wang, Xinlei; Qiao, Jun; Xu, Jiangtao

    2016-01-01

    This paper presents a full parallel event driven readout method which is implemented in an area array single-photon avalanche diode (SPAD) image sensor for high-speed fluorescence lifetime imaging microscopy (FLIM). The sensor only records and reads out effective time and position information by adopting full parallel event driven readout method, aiming at reducing the amount of data. The image sensor includes four 8 × 8 pixel arrays. In each array, four time-to-digital converters (TDCs) are used to quantize the time of photons' arrival, and two address record modules are used to record the column and row information. In this work, Monte Carlo simulations were performed in Matlab in terms of the pile-up effect induced by the readout method. The sensor's resolution is 16 × 16. The time resolution of TDCs is 97.6 ps and the quantization range is 100 ns. The readout frame rate is 10 Mfps, and the maximum imaging frame rate is 100 fps. The chip's output bandwidth is 720 MHz with an average power of 15 mW. The lifetime resolvability range is 5-20 ns, and the average error of estimated fluorescence lifetimes is below 1% by employing CMM to estimate lifetimes. PMID:26828490

  10. Event-Driven Random-Access-Windowing CCD Imaging System

    NASA Technical Reports Server (NTRS)

    Monacos, Steve; Portillo, Angel; Ortiz, Gerardo; Alexander, James; Lam, Raymond; Liu, William

    2004-01-01

    A charge-coupled-device (CCD) based high-speed imaging system, called a realtime, event-driven (RARE) camera, is undergoing development. This camera is capable of readout from multiple subwindows [also known as regions of interest (ROIs)] within the CCD field of view. Both the sizes and the locations of the ROIs can be controlled in real time and can be changed at the camera frame rate. The predecessor of this camera was described in High-Frame-Rate CCD Camera Having Subwindow Capability (NPO- 30564) NASA Tech Briefs, Vol. 26, No. 12 (December 2002), page 26. The architecture of the prior camera requires tight coupling between camera control logic and an external host computer that provides commands for camera operation and processes pixels from the camera. This tight coupling limits the attainable frame rate and functionality of the camera. The design of the present camera loosens this coupling to increase the achievable frame rate and functionality. From a host computer perspective, the readout operation in the prior camera was defined on a per-line basis; in this camera, it is defined on a per-ROI basis. In addition, the camera includes internal timing circuitry. This combination of features enables real-time, event-driven operation for adaptive control of the camera. Hence, this camera is well suited for applications requiring autonomous control of multiple ROIs to track multiple targets moving throughout the CCD field of view. Additionally, by eliminating the need for control intervention by the host computer during the pixel readout, the present design reduces ROI-readout times to attain higher frame rates. This camera (see figure) includes an imager card consisting of a commercial CCD imager and two signal-processor chips. The imager card converts transistor/ transistor-logic (TTL)-level signals from a field programmable gate array (FPGA) controller card. These signals are transmitted to the imager card via a low-voltage differential signaling (LVDS) cable

  11. General Purpose Heat Source Simulator

    NASA Technical Reports Server (NTRS)

    Emrich, William J., Jr.

    2008-01-01

    The General Purpose Heat Source (GPHS) project seeks to combine the development of an electrically heated, single GPHS module simulator with the evaluation of potential nuclear surface power systems. The simulator is designed to match the form, fit, and function of actual GPHS modules which normally generate heat through the radioactive decay of Pu238. The use of electrically heated modules rather than modules containing Pu238 facilitates the testing of the subsystems and systems without sacrificing the quantity and quality of the test data gathered. Current GPHS activities are centered on developing robust heater designs with sizes and weights which closely match those of actual Pu238 fueled GPHS blocks. Designs are being pursued which will allow operation up to 1100 C.

  12. Event Driven Messaging with Role-Based Subscriptions

    NASA Technical Reports Server (NTRS)

    Bui, Tung; Bui, Bach; Malhotra, Shantanu; Chen, Fannie; Kim, rachel; Allen, Christopher; Luong, Ivy; Chang, George; Zendejas, Silvino; Sadaqathulla, Syed

    2009-01-01

    Event Driven Messaging with Role-Based Subscriptions (EDM-RBS) is a framework integrated into the Service Management Database (SMDB) to allow for role-based and subscription-based delivery of synchronous and asynchronous messages over JMS (Java Messaging Service), SMTP (Simple Mail Transfer Protocol), or SMS (Short Messaging Service). This allows for 24/7 operation with users in all parts of the world. The software classifies messages by triggering data type, application source, owner of data triggering event (mission), classification, sub-classification and various other secondary classifying tags. Messages are routed to applications or users based on subscription rules using a combination of the above message attributes. This program provides a framework for identifying connected users and their applications for targeted delivery of messages over JMS to the client applications the user is logged into. EDMRBS provides the ability to send notifications over e-mail or pager rather than having to rely on a live human to do it. It is implemented as an Oracle application that uses Oracle relational database management system intrinsic functions. It is configurable to use Oracle AQ JMS API or an external JMS provider for messaging. It fully integrates into the event-logging framework of SMDB (Subnet Management Database).

  13. Event-driven contrastive divergence for spiking neuromorphic systems

    PubMed Central

    Neftci, Emre; Das, Srinjoy; Pedroni, Bruno; Kreutz-Delgado, Kenneth; Cauwenberghs, Gert

    2014-01-01

    Restricted Boltzmann Machines (RBMs) and Deep Belief Networks have been demonstrated to perform efficiently in a variety of applications, such as dimensionality reduction, feature learning, and classification. Their implementation on neuromorphic hardware platforms emulating large-scale networks of spiking neurons can have significant advantages from the perspectives of scalability, power dissipation and real-time interfacing with the environment. However, the traditional RBM architecture and the commonly used training algorithm known as Contrastive Divergence (CD) are based on discrete updates and exact arithmetics which do not directly map onto a dynamical neural substrate. Here, we present an event-driven variation of CD to train a RBM constructed with Integrate & Fire (I&F) neurons, that is constrained by the limitations of existing and near future neuromorphic hardware platforms. Our strategy is based on neural sampling, which allows us to synthesize a spiking neural network that samples from a target Boltzmann distribution. The recurrent activity of the network replaces the discrete steps of the CD algorithm, while Spike Time Dependent Plasticity (STDP) carries out the weight updates in an online, asynchronous fashion. We demonstrate our approach by training an RBM composed of leaky I&F neurons with STDP synapses to learn a generative model of the MNIST hand-written digit dataset, and by testing it in recognition, generation and cue integration tasks. Our results contribute to a machine learning-driven approach for synthesizing networks of spiking neurons capable of carrying out practical, high-level functionality. PMID:24574952

  14. Multiagent Attitude Control System for Satellites Based in Momentum Wheels and Event-Driven Synchronization

    NASA Astrophysics Data System (ADS)

    Garcia, Juan L.; Moreno, Jose Sanchez

    2012-12-01

    Attitude control is a requirement always present in spacecraft design. Several kinds of actuators exist to accomplish this control, being momentum wheels one of the most employed. Usually satellites carry redundant momentum wheels to handle any possible single failure, but the controller remains as a single centralized element, posing problems in case of failures. In this work a decentralized agent-based event-driven algorithm for attitude control is presented as a possible solution. Several agents based in momentum wheels will interact among them to accomplish the satellite control. A simulation environment has been developed to analyze the behavior of this architecture. This environment has been made available through the web page http://www.dia.uned.es.

  15. A hybrid adaptive routing algorithm for event-driven wireless sensor networks.

    PubMed

    Figueiredo, Carlos M S; Nakamura, Eduardo F; Loureiro, Antonio A F

    2009-01-01

    Routing is a basic function in wireless sensor networks (WSNs). For these networks, routing algorithms depend on the characteristics of the applications and, consequently, there is no self-contained algorithm suitable for every case. In some scenarios, the network behavior (traffic load) may vary a lot, such as an event-driven application, favoring different algorithms at different instants. This work presents a hybrid and adaptive algorithm for routing in WSNs, called Multi-MAF, that adapts its behavior autonomously in response to the variation of network conditions. In particular, the proposed algorithm applies both reactive and proactive strategies for routing infrastructure creation, and uses an event-detection estimation model to change between the strategies and save energy. To show the advantages of the proposed approach, it is evaluated through simulations. Comparisons with independent reactive and proactive algorithms show improvements on energy consumption. PMID:22423207

  16. General Purpose Heat Source Simulator

    NASA Technical Reports Server (NTRS)

    Emrich, Bill

    2008-01-01

    The General Purpose Heat Source (GPHS) simulator project is designed to replicate through the use of electrical heaters, the form, fit, and function of actual GPHS modules which generate heat through the radioactive decay of Pu238. The use of electrically heated modules rather than modules containing Pu238 facilitates the testing of spacecraft subsystems and systems without sacrificing the quantity and quality of the test data gathered. Previous GPHS activities are centered around developing robust heater designs with sizes and weights that closely matched those of actual Pu238 fueled GPHS blocks. These efforts were successful, although their maximum temperature capabilities were limited to around 850 C. New designs are being pursued which also replicate the sizes and weights of actual Pu238 fueled GPHS blocks but will allow operation up to 1100 C.

  17. Generalized Fluid System Simulation Program

    NASA Technical Reports Server (NTRS)

    Majumdar, Alok Kumar (Inventor); Bailey, John W. (Inventor); Schallhorn, Paul Alan (Inventor); Steadman, Todd E. (Inventor)

    2004-01-01

    A general purpose program implemented on a computer analyzes steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal force. A preprocessor provides for the inter- active development of a fluid network simulation having nodes and branches. Mass, energy, and specie conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. Contained herein are subroutines for computing "real fluid" thermodynamic and thermophysical properties for 12 fluids, and a number of different source options are provided for model- ing momentum sources or sinks in the branches. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. Application and verification of this invention are provided through an example problem, which demonstrates that the predictions of the present invention compare most reasonably with test data.

  18. General Reactive Atomistic Simulation Program

    Energy Science and Technology Software Center (ESTSC)

    2004-09-22

    GRASP (General Reactive Atomistic Simulation Program) is primarily intended as a molecular dynamics package for complex force fields, The code is designed to provide good performance for large systems, either in parallel or serial execution mode, The primary purpose of the code is to realistically represent the structural and dynamic properties of large number of atoms on timescales ranging from picoseconds up to a microsecond. Typically the atoms form a representative sample of some material,more » such as an interface between polycrystalline silicon and amorphous silica. GRASP differs from other parallel molecular dynamics codes primarily due to it’s ability to handle relatively complicated interaction potentials and it’s ability to use more than one interaction potential in a single simulation. Most of the computational effort goes into the calculation of interatomic forces, which depend in a complicated way on the positions of all the atoms. The forces are used to integrate the equations of motion forward in time using the so-called velocity Verlet integration scheme. Alternatively, the forces can be used to find a minimum energy configuration, in which case a modified steepest descent algorithm is used.« less

  19. Mapping from frame-driven to frame-free event-driven vision systems by low-rate rate coding and coincidence processing--application to feedforward ConvNets.

    PubMed

    Pérez-Carrasco, José Antonio; Zhao, Bo; Serrano, Carmen; Acha, Begoña; Serrano-Gotarredona, Teresa; Chen, Shouchun; Linares-Barranco, Bernabé

    2013-11-01

    Event-driven visual sensors have attracted interest from a number of different research communities. They provide visual information in quite a different way from conventional video systems consisting of sequences of still images rendered at a given "frame rate." Event-driven vision sensors take inspiration from biology. Each pixel sends out an event (spike) when it senses something meaningful is happening, without any notion of a frame. A special type of event-driven sensor is the so-called dynamic vision sensor (DVS) where each pixel computes relative changes of light or "temporal contrast." The sensor output consists of a continuous flow of pixel events that represent the moving objects in the scene. Pixel events become available with microsecond delays with respect to "reality." These events can be processed "as they flow" by a cascade of event (convolution) processors. As a result, input and output event flows are practically coincident in time, and objects can be recognized as soon as the sensor provides enough meaningful events. In this paper, we present a methodology for mapping from a properly trained neural network in a conventional frame-driven representation to an event-driven representation. The method is illustrated by studying event-driven convolutional neural networks (ConvNet) trained to recognize rotating human silhouettes or high speed poker card symbols. The event-driven ConvNet is fed with recordings obtained from a real DVS camera. The event-driven ConvNet is simulated with a dedicated event-driven simulator and consists of a number of event-driven processing modules, the characteristics of which are obtained from individually manufactured hardware modules. PMID:24051730

  20. Mapping from Frame-Driven to Frame-Free Event-Driven Vision Systems by Low-Rate Rate-Coding and Coincidence Processing. Application to Feed Forward ConvNets.

    PubMed

    Perez-Carrasco, J A; Zhao, B; Serrano, C; Acha, B; Serrano-Gotarredona, T; Chen, S; Linares-Barranco, B

    2013-04-10

    Event-driven visual sensors have attracted interest from a number of different research communities. They provide visual information in quite a different way from conventional video systems consisting of sequences of still images rendered at “frame rate”. Event-driven vision sensors take inspiration from biology. A special type of Event-driven sensor is the so called Dynamic-Vision-Sensor (DVS) where each pixel computes relative changes of light, or “temporal contrast”. Pixel events become available with micro second delays with respect to “reality”. These events can be processed “as they flow” by a cascade of event (convolution) processors. As a result, input and output event flows are practically coincident, and objects can be recognized as soon as the sensor provides enough meaningful events. In this paper we present a methodology for mapping from a properly trained neural network in a conventional Frame-driven representation, to an Event-driven representation. The method is illustrated by studying Event-driven Convolutional Neural Networks (ConvNet) trained to recognize rotating human silhouettes or high speed poker card symbols. The Event-driven ConvNet is fed with recordings obtained from a real DVS camera. The Event-driven ConvNet is simulated with a dedicated Event-driven simulator, and consists of a number of Event-driven processing modules the characteristics of which are obtained from individually manufactured hardware modules. PMID:23589589

  1. Flow simulation on generalized grids

    SciTech Connect

    Koomullil, R.P.; Soni, B.K.; Huang, Chi Ti

    1996-12-31

    A hybrid grid generation methodology and flow simulation on grids having an arbitrary number of sided polygons is presented. A hyperbolic type marching scheme is used for generating structured grids near the solid boundaries. A local elliptic solver is utilized for smoothing the grid lines and for avoiding grid line crossing. A new method for trimming the overlaid structured grid is presented. Delaunay triangulation is employed to generate an unstructured grid in the regions away from the body. The structured and unstructured grid regions are integrated together to form a single grid for the flow solver. An edge based data structure is used to store the grid information to ease the handling of general polygons. Integral form of the Navier-Stokes equations makes up the governing equations. A Roe averaged Riemann solver is utilized to evaluate the numerical flux at cell faces. Higher order accuracy is achieved by applying Taylor`s series expansion to the conserved variables, and the gradient is calculated by using Green`s theorem. For the implicit scheme, the sparse matrix resulting from the linearization is solved using GMRES method. The flux Jacobians are calculated numerically or by an approximate analytic method. Results are presented to validate the current methodology.

  2. Gaming Simulation: A General Classification.

    ERIC Educational Resources Information Center

    Cecchini, Arnaldo; Frisenna, Adriana

    1987-01-01

    Reviews the problems of classifying gaming techniques and suggests a heuristic approach as one solution. Definitions of simulation, models, role, and game and play are discussed to help develop a classification based on a technique called gaming simulation. (Author/LRW)

  3. Event-Driven X-Ray CCD Detectors for High Energy Astrophysics

    NASA Technical Reports Server (NTRS)

    Ricker, George R.

    2004-01-01

    A viewgraph presentation describing the Event-Driven X- Ray CCD (EDCCD) detector system for high energy astrophysics is presented. The topics include: 1) EDCCD: Description and Advantages; 2) Summary of Grant Activity Carried Out; and 3) EDCCD Test System.

  4. Event-driven model predictive control of sewage pumping stations for sulfide mitigation in sewer networks.

    PubMed

    Liu, Yiqi; Ganigué, Ramon; Sharma, Keshab; Yuan, Zhiguo

    2016-07-01

    Chemicals such as Mg(OH)2 and iron salts are widely dosed to sewage for mitigating sulfide-induced corrosion and odour problems in sewer networks. The chemical dosing rate is usually not automatically controlled but profiled based on experience of operators, often resulting in over- or under-dosing. Even though on-line control algorithms for chemical dosing in single pipes have been developed recently, network-wide control algorithms are currently not available. The key challenge is that a sewer network is typically wide-spread comprising many interconnected sewer pipes and pumping stations, making network-wide sulfide mitigation with a relatively limited number of dosing points challenging. In this paper, we propose and demonstrate an Event-driven Model Predictive Control (EMPC) methodology, which controls the flows of sewage streams containing the dosed chemical to ensure desirable distribution of the dosed chemical throughout the pipe sections of interests. First of all, a network-state model is proposed to predict the chemical concentration in a network. An EMPC algorithm is then designed to coordinate sewage pumping station operations to ensure desirable chemical distribution in the network. The performance of the proposed control methodology is demonstrated by applying the designed algorithm to a real sewer network simulated with the well-established SeweX model using real sewage flow and characteristics data. The EMPC strategy significantly improved the sulfide mitigation performance with the same chemical consumption, compared to the current practice. PMID:27124127

  5. Field Evaluation of a General Purpose Simulator.

    ERIC Educational Resources Information Center

    Spangenberg, Ronald W.

    The use of a general purpose simulator (GPS) to teach Air Force technicians diagnostic and repair procedures for specialized aircraft radar systems is described. An EC II simulator manufactured by Educational Computer Corporation was adapted to resemble the actual configuration technicians would encounter in the field. Data acquired in the…

  6. Simulations in generalized ensembles through noninstantaneous switches

    NASA Astrophysics Data System (ADS)

    Giovannelli, Edoardo; Cardini, Gianni; Chelli, Riccardo

    2015-10-01

    Generalized-ensemble simulations, such as replica exchange and serial generalized-ensemble methods, are powerful simulation tools to enhance sampling of free energy landscapes in systems with high energy barriers. In these methods, sampling is enhanced through instantaneous transitions of replicas, i.e., copies of the system, between different ensembles characterized by some control parameter associated with thermodynamical variables (e.g., temperature or pressure) or collective mechanical variables (e.g., interatomic distances or torsional angles). An interesting evolution of these methodologies has been proposed by replacing the conventional instantaneous (trial) switches of replicas with noninstantaneous switches, realized by varying the control parameter in a finite time and accepting the final replica configuration with a Metropolis-like criterion based on the Crooks nonequilibrium work (CNW) theorem. Here we revise these techniques focusing on their correlation with the CNW theorem in the framework of Markovian processes. An outcome of this report is the derivation of the acceptance probability for noninstantaneous switches in serial generalized-ensemble simulations, where we show that explicit knowledge of the time dependence of the weight factors entering such simulations is not necessary. A generalized relationship of the CNW theorem is also provided in terms of the underlying equilibrium probability distribution at a fixed control parameter. Illustrative calculations on a toy model are performed with serial generalized-ensemble simulations, especially focusing on the different behavior of instantaneous and noninstantaneous replica transition schemes.

  7. Simulations in generalized ensembles through noninstantaneous switches.

    PubMed

    Giovannelli, Edoardo; Cardini, Gianni; Chelli, Riccardo

    2015-10-01

    Generalized-ensemble simulations, such as replica exchange and serial generalized-ensemble methods, are powerful simulation tools to enhance sampling of free energy landscapes in systems with high energy barriers. In these methods, sampling is enhanced through instantaneous transitions of replicas, i.e., copies of the system, between different ensembles characterized by some control parameter associated with thermodynamical variables (e.g., temperature or pressure) or collective mechanical variables (e.g., interatomic distances or torsional angles). An interesting evolution of these methodologies has been proposed by replacing the conventional instantaneous (trial) switches of replicas with noninstantaneous switches, realized by varying the control parameter in a finite time and accepting the final replica configuration with a Metropolis-like criterion based on the Crooks nonequilibrium work (CNW) theorem. Here we revise these techniques focusing on their correlation with the CNW theorem in the framework of Markovian processes. An outcome of this report is the derivation of the acceptance probability for noninstantaneous switches in serial generalized-ensemble simulations, where we show that explicit knowledge of the time dependence of the weight factors entering such simulations is not necessary. A generalized relationship of the CNW theorem is also provided in terms of the underlying equilibrium probability distribution at a fixed control parameter. Illustrative calculations on a toy model are performed with serial generalized-ensemble simulations, especially focusing on the different behavior of instantaneous and noninstantaneous replica transition schemes. PMID:26565367

  8. Connection between Newtonian simulations and general relativity

    SciTech Connect

    Chisari, Nora Elisa; Zaldarriaga, Matias

    2011-06-15

    On large scales, comparable to the horizon, the observable clustering properties of galaxies are affected by various general relativistic effects. To calculate these effects one needs to consistently solve for the metric, densities, and velocities in a specific coordinate system or gauge. The method of choice for simulating large-scale structure is numerical N-body simulations which are performed in the Newtonian limit. Even though one might worry that the use of the Newtonian approximation would make it impossible to use these simulations to compute properties on very large scales, we show that the simulations are still solving the dynamics correctly even for long modes and we give formulas to obtain the position of particles in the conformal Newtonian gauge given the positions computed in the simulation. We also give formulas to convert from the output coordinates of N-body simulations to the observable coordinates of the particles.

  9. Notification Event Architecture for Traveler Screening: Predictive Traveler Screening Using Event Driven Business Process Management

    ERIC Educational Resources Information Center

    Lynch, John Kenneth

    2013-01-01

    Using an exploratory model of the 9/11 terrorists, this research investigates the linkages between Event Driven Business Process Management (edBPM) and decision making. Although the literature on the role of technology in efficient and effective decision making is extensive, research has yet to quantify the benefit of using edBPM to aid the…

  10. A general software reliability process simulation technique

    NASA Technical Reports Server (NTRS)

    Tausworthe, Robert C.

    1991-01-01

    The structure and rationale of the generalized software reliability process, together with the design and implementation of a computer program that simulates this process are described. Given assumed parameters of a particular project, the users of this program are able to generate simulated status timelines of work products, numbers of injected anomalies, and the progress of testing, fault isolation, repair, validation, and retest. Such timelines are useful in comparison with actual timeline data, for validating the project input parameters, and for providing data for researchers in reliability prediction modeling.

  11. Spectral Methods in General Relativistic MHD Simulations

    NASA Astrophysics Data System (ADS)

    Garrison, David

    2012-03-01

    In this talk I discuss the use of spectral methods in improving the accuracy of a General Relativistic Magnetohydrodynamic (GRMHD) computer code. I introduce SpecCosmo, a GRMHD code developed as a Cactus arrangement at UHCL, and show simulation results using both Fourier spectral methods and finite differencing. This work demonstrates the use of spectral methods with the FFTW 3.3 Fast Fourier Transform package integrated with the Cactus Framework to perform spectral differencing using MPI.

  12. Simulation of General Physics laboratory exercise

    NASA Astrophysics Data System (ADS)

    Aceituno, P.; Hernández-Aceituno, J.; Hernández-Cabrera, A.

    2015-01-01

    Laboratory exercises are an important part of general Physics teaching, both during the last years of high school and the first year of college education. Due to the need to acquire enough laboratory equipment for all the students, and the widespread access to computers rooms in teaching, we propose the development of computer simulated laboratory exercises. A representative exercise in general Physics is the calculation of the gravity acceleration value, through the free fall motion of a metal ball. Using a model of the real exercise, we have developed an interactive system which allows students to alter the starting height of the ball to obtain different fall times. The simulation was programmed in ActionScript 3, so that it can be freely executed in any operative system; to ensure the accuracy of the calculations, all the input parameters of the simulations were modelled using digital measurement units, and to allow a statistical management of the resulting data, measurement errors are simulated through limited randomization.

  13. General Relativistic Magnetohydrodynamic Simulations of Collapsars

    NASA Technical Reports Server (NTRS)

    Mizuno, Yosuke; Yamada, S.; Koider, S.; Shipata, K.

    2005-01-01

    We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. Initially, we assume that the core collapse has failed in this star. A rotating black hole of a few solar masses is inserted by hand into the calculation. The simulation results show the formation of a disklike structure and the generation of a jetlike outflow near the central black hole. The jetlike outflow propagates and accelerated mainly by the magnetic field. The total jet velocity is approximately 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is directly converted to kinetic energy of the jet rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster.

  14. A general formulation for compositional reservoir simulation

    SciTech Connect

    Rodriguez, F.; Guzman, J.; Galindo-Nava, A.

    1994-12-31

    In this paper the authors present a general formulation to solve the non-linear difference equations that arise in compositional reservoir simulation. The general approach here presented is based on newton`s method and provides a systematic approach to generate several formulations to solve the compositional problem, each possessing a different degree of implicitness and stability characteristics. The Fully-Implicit method is at the higher end of the implicitness spectrum while the IMPECS method, implicit in pressure-explicit in composition and saturation, is at the lower end. They show that all methods may be obtained as particular cases of the fully-implicit method. Regarding the matrix problem, all methods have a similar matrix structure; the composition of the Jacobian matrix is however unique in each case, being in some instances amenable to reductions for optimal solution of the matrix problem. Based on this, a different approach to derive IMPECS type methods is proposed; in this case, the whole set of 2nc + 6 equations, that apply in each gridblock, is reduced to a single pressure equation through matrix reduction operations; this provides a more stable numerical scheme, compared to other published IMPCS methods, in which the subset of thermodynamic equilibrium equations is arbitrarily decoupled form the set of gridblock equations to perform such reduction. The authors discuss how the general formulation here presented can be used to formulate and construct an adaptive-implicit compositional simulators. They also present results on the numerical performance of FI, IMPSEC and IMPECS methods on some test problems.

  15. Modeling the Energy Performance of Event-Driven Wireless Sensor Network by Using Static Sink and Mobile Sink

    PubMed Central

    Chen, Jiehui; Salim, Mariam B.; Matsumoto, Mitsuji

    2010-01-01

    Wireless Sensor Networks (WSNs) designed for mission-critical applications suffer from limited sensing capacities, particularly fast energy depletion. Regarding this, mobile sinks can be used to balance the energy consumption in WSNs, but the frequent location updates of the mobile sinks can lead to data collisions and rapid energy consumption for some specific sensors. This paper explores an optimal barrier coverage based sensor deployment for event driven WSNs where a dual-sink model was designed to evaluate the energy performance of not only static sensors, but Static Sink (SS) and Mobile Sinks (MSs) simultaneously, based on parameters such as sensor transmission range r and the velocity of the mobile sink v, etc. Moreover, a MS mobility model was developed to enable SS and MSs to effectively collaborate, while achieving spatiotemporal energy performance efficiency by using the knowledge of the cumulative density function (cdf), Poisson process and M/G/1 queue. The simulation results verified that the improved energy performance of the whole network was demonstrated clearly and our eDSA algorithm is more efficient than the static-sink model, reducing energy consumption approximately in half. Moreover, we demonstrate that our results are robust to realistic sensing models and also validate the correctness of our results through extensive simulations. PMID:22163503

  16. Event-driven charge-coupled device design and applications therefor

    NASA Technical Reports Server (NTRS)

    Doty, John P. (Inventor); Ricker, Jr., George R. (Inventor); Burke, Barry E. (Inventor); Prigozhin, Gregory Y. (Inventor)

    2005-01-01

    An event-driven X-ray CCD imager device uses a floating-gate amplifier or other non-destructive readout device to non-destructively sense a charge level in a charge packet associated with a pixel. The output of the floating-gate amplifier is used to identify each pixel that has a charge level above a predetermined threshold. If the charge level is above a predetermined threshold the charge in the triggering charge packet and in the charge packets from neighboring pixels need to be measured accurately. A charge delay register is included in the event-driven X-ray CCD imager device to enable recovery of the charge packets from neighboring pixels for accurate measurement. When a charge packet reaches the end of the charge delay register, control logic either dumps the charge packet, or steers the charge packet to a charge FIFO to preserve it if the charge packet is determined to be a packet that needs accurate measurement. A floating-diffusion amplifier or other low-noise output stage device, which converts charge level to a voltage level with high precision, provides final measurement of the charge packets. The voltage level is eventually digitized by a high linearity ADC.

  17. Event-Driven Observations and Comprehensive Evaluation for Natural Disaster Assessment in China

    NASA Astrophysics Data System (ADS)

    Wang, Z.; Li, Z.; Shen, Y.; Wu, L.; Li, H.

    2012-08-01

    The Chinese event-driven observations and disaster assessment system is established so as to make information related to environmental risk and vulnerability easily accessible to decision-makers through this centralized platform. At 7:49 AM on April 14, 2010, an earthquake of 7.1 magnitude collapsed buildings in Yushu County, Yushu Tibetan Autonomous Prefecture in Qinghai Province. For emergence response, we presented a method for generating seismic intensity map based on seismological mechanism solutions. The disaster assessment system automatically drew the distribution map of affected population 1 hour after the Yushu earthquake. In the case of distribution map of affected population, we made the judgment that the Gyêgu town maybe the worst hit town in the Yushu earthquake because it is not only near the epicenter, but also the capital of Yushu Tibetan Autonomous Prefecture. Then event-driven observations are taken on Gyêgu town. Referring to usable data, the chains of rapid assessment on earthquake disaster were analyzed, and different models were established for assessing affected population, damaged houses and lifelines and comprehensive earthquake loss evaluation.

  18. EVENT DRIVEN AUTOMATIC STATE MODIFICATION OF BNL'S BOOSTER FOR NASA SPACE RADIATION LABORATORY SOLAR PARTICLE SIMULATOR.

    SciTech Connect

    BROWN, D.; BINELLO, S.; HARVEY, M.; MORRIS, J.; RUSEK, A.; TSOUPAS, N.

    2005-05-16

    The NASA Space Radiation Laboratory (NSRL) was constructed in collaboration with NASA for the purpose of performing radiation effect studies for the NASA space program. The NSRL makes use of heavy ions in the range of 0.05 to 3 GeV/n slow extracted from BNL's AGS Booster. NASA is interested in reproducing the energy spectrum from a solar flare in the space environment for a single ion species. To do this we have built and tested a set of software tools which allow the state of the Booster and the NSRL beam line to be changed automatically. In this report we will describe the system and present results of beam tests.

  19. Event-driven visual attention for the humanoid robot iCub.

    PubMed

    Rea, Francesco; Metta, Giorgio; Bartolozzi, Chiara

    2013-01-01

    Fast reaction to sudden and potentially interesting stimuli is a crucial feature for safe and reliable interaction with the environment. Here we present a biologically inspired attention system developed for the humanoid robot iCub. It is based on input from unconventional event-driven vision sensors and an efficient computational method. The resulting system shows low-latency and fast determination of the location of the focus of attention. The performance is benchmarked against an instance of the state of the art in robotics artificial attention system used in robotics. Results show that the proposed system is two orders of magnitude faster that the benchmark in selecting a new stimulus to attend. PMID:24379753

  20. Event-driven approach of layered multicast to network adaptation in RED-based IP networks

    NASA Astrophysics Data System (ADS)

    Nahm, Kitae; Li, Qing; Kuo, C.-C. J.

    2003-11-01

    In this work, we investigate the congestion control problem for layered video multicast in IP networks of active queue management (AQM) using a simple random early detection (RED) queue model. AQM support from networks improves the visual quality of video streaming but makes network adaptation more di+/-cult for existing layered video multicast proticols that use the event-driven timer-based approach. We perform a simplified analysis on the response of the RED algorithm to burst traffic. The analysis shows that the primary problem lies in the weak correlation between the network feedback and the actual network congestion status when the RED queue is driven by burst traffic. Finally, a design guideline of the layered multicast protocol is proposed to overcome this problem.

  1. Event-driven visual attention for the humanoid robot iCub

    PubMed Central

    Rea, Francesco; Metta, Giorgio; Bartolozzi, Chiara

    2013-01-01

    Fast reaction to sudden and potentially interesting stimuli is a crucial feature for safe and reliable interaction with the environment. Here we present a biologically inspired attention system developed for the humanoid robot iCub. It is based on input from unconventional event-driven vision sensors and an efficient computational method. The resulting system shows low-latency and fast determination of the location of the focus of attention. The performance is benchmarked against an instance of the state of the art in robotics artificial attention system used in robotics. Results show that the proposed system is two orders of magnitude faster that the benchmark in selecting a new stimulus to attend. PMID:24379753

  2. Design of an Event-Driven Random-Access-Windowing CCD-Based Camera

    NASA Technical Reports Server (NTRS)

    Monacos, Steve P.; Lam, Raymond K.; Portillo, Angel A.; Ortiz, Gerardo G.

    2003-01-01

    Commercially available cameras are not design for the combination of single frame and high-speed streaming digital video with real-time control of size and location of multiple regions-of-interest (ROI). A new control paradigm is defined to eliminate the tight coupling between the camera logic and the host controller. This functionality is achieved by defining the indivisible pixel read out operation on a per ROI basis with in-camera time keeping capability. This methodology provides a Random Access, Real-Time, Event-driven (RARE) camera for adaptive camera control and is will suited for target tracking applications requiring autonomous control of multiple ROI's. This methodology additionally provides for reduced ROI read out time and higher frame rates compared to the original architecture by avoiding external control intervention during the ROI read out process.

  3. FusionAnalyser: a new graphical, event-driven tool for fusion rearrangements discovery.

    PubMed

    Piazza, Rocco; Pirola, Alessandra; Spinelli, Roberta; Valletta, Simona; Redaelli, Sara; Magistroni, Vera; Gambacorti-Passerini, Carlo

    2012-09-01

    Gene fusions are common driver events in leukaemias and solid tumours; here we present FusionAnalyser, a tool dedicated to the identification of driver fusion rearrangements in human cancer through the analysis of paired-end high-throughput transcriptome sequencing data. We initially tested FusionAnalyser by using a set of in silico randomly generated sequencing data from 20 known human translocations occurring in cancer and subsequently using transcriptome data from three chronic and three acute myeloid leukaemia samples. in all the cases our tool was invariably able to detect the presence of the correct driver fusion event(s) with high specificity. In one of the acute myeloid leukaemia samples, FusionAnalyser identified a novel, cryptic, in-frame ETS2-ERG fusion. A fully event-driven graphical interface and a flexible filtering system allow complex analyses to be run in the absence of any a priori programming or scripting knowledge. Therefore, we propose FusionAnalyser as an efficient and robust graphical tool for the identification of functional rearrangements in the context of high-throughput transcriptome sequencing data. PMID:22570408

  4. GTOSS: Generalized Tethered Object Simulation System

    NASA Technical Reports Server (NTRS)

    Lang, David D.

    1987-01-01

    GTOSS represents a tether analysis complex which is described by addressing its family of modules. TOSS is a portable software subsystem specifically designed to be introduced into the environment of any existing vehicle dynamics simulation to add the capability of simulating multiple interacting objects (via multiple tethers). These objects may interact with each other as well as with the vehicle into whose environment TOSS is introduced. GTOSS is a stand alone tethered system analysis program, representing an example of TOSS having been married to a host simulation. RTOSS is the Results Data Base (RDB) subsystem designed to archive TOSS simulation results for future display processing. DTOSS is a display post processors designed to utilize the RDB. DTOSS extracts data from the RDB for multi-page printed time history displays. CTOSS is similar to DTOSS, but is designed to create ASCII plot files. The same time history data formats provided for DTOSS (for printing) are available via CTOSS for plotting. How these and other modules interact with each other is discussed.

  5. A 300-mV 220-nW event-driven ADC with real-time QRS detection for wearable ECG sensors.

    PubMed

    Zhang, Xiaoyang; Lian, Yong

    2014-12-01

    This paper presents an ultra-low-power event-driven analog-to-digital converter (ADC) with real-time QRS detection for wearable electrocardiogram (ECG) sensors in wireless body sensor network (WBSN) applications. Two QRS detection algorithms, pulse-triggered (PUT) and time-assisted PUT (t-PUT), are proposed based on the level-crossing events generated from the ADC. The PUT detector achieves 97.63% sensitivity and 97.33% positive prediction in simulation on the MIT-BIH Arrhythmia Database. The t-PUT improves the sensitivity and positive prediction to 97.76% and 98.59% respectively. Fabricated in 0.13 μm CMOS technology, the ADC with QRS detector consumes only 220 nW measured under 300 mV power supply, making it the first nanoWatt compact analog-to-information (A2I) converter with embedded QRS detector. PMID:25608283

  6. Real-time gesture interface based on event-driven processing from stereo silicon retinas.

    PubMed

    Lee, Jun Haeng; Delbruck, Tobi; Pfeiffer, Michael; Park, Paul K J; Shin, Chang-Woo; Ryu, Hyunsurk Eric; Kang, Byung Chang

    2014-12-01

    We propose a real-time hand gesture interface based on combining a stereo pair of biologically inspired event-based dynamic vision sensor (DVS) silicon retinas with neuromorphic event-driven postprocessing. Compared with conventional vision or 3-D sensors, the use of DVSs, which output asynchronous and sparse events in response to motion, eliminates the need to extract movements from sequences of video frames, and allows significantly faster and more energy-efficient processing. In addition, the rate of input events depends on the observed movements, and thus provides an additional cue for solving the gesture spotting problem, i.e., finding the onsets and offsets of gestures. We propose a postprocessing framework based on spiking neural networks that can process the events received from the DVSs in real time, and provides an architecture for future implementation in neuromorphic hardware devices. The motion trajectories of moving hands are detected by spatiotemporally correlating the stereoscopically verged asynchronous events from the DVSs by using leaky integrate-and-fire (LIF) neurons. Adaptive thresholds of the LIF neurons achieve the segmentation of trajectories, which are then translated into discrete and finite feature vectors. The feature vectors are classified with hidden Markov models, using a separate Gaussian mixture model for spotting irrelevant transition gestures. The disparity information from stereovision is used to adapt LIF neuron parameters to achieve recognition invariant of the distance of the user to the sensor, and also helps to filter out movements in the background of the user. Exploiting the high dynamic range of DVSs, furthermore, allows gesture recognition over a 60-dB range of scene illuminance. The system achieves recognition rates well over 90% under a variety of variable conditions with static and dynamic backgrounds with naïve users. PMID:25420246

  7. General Case Simulation Instruction of Generalized Housekeeping Skills in Blind, Multihandicapped Adults.

    ERIC Educational Resources Information Center

    Lengyel, L. M.; And Others

    1990-01-01

    The study with three blind and mentally retarded adults with additional disabilities found that general case simulation instruction in housekeeping skills led to generalization to untrained settings. Degree of generalization was inversely related to the severity and complexity of participant disability. (Author/DB)

  8. General Aviation Cockpit Weather Information System Simulation Studies

    NASA Technical Reports Server (NTRS)

    McAdaragh, Ray; Novacek, Paul

    2003-01-01

    This viewgraph presentation provides information on two experiments on the effectiveness of a cockpit weather information system on a simulated general aviation flight. The presentation covers the simulation hardware configuration, the display device screen layout, a mission scenario, conclusions, and recommendations. The second experiment, with its own scenario and conclusions, is a follow-on experiment.

  9. Test and evaluation of the generalized gate logic system simulator

    NASA Technical Reports Server (NTRS)

    Miner, Paul S.

    1991-01-01

    The results of the initial testing of the Generalized Gate Level Logic Simulator (GGLOSS) are discussed. The simulator is a special purpose fault simulator designed to assist in the analysis of the effects of random hardware failures on fault tolerant digital computer systems. The testing of the simulator covers two main areas. First, the simulation results are compared with data obtained by monitoring the behavior of hardware. The circuit used for these comparisons is an incomplete microprocessor design based upon the MIL-STD-1750A Instruction Set Architecture. In the second area of testing, current simulation results are compared with experimental data obtained using precursors of the current tool. In each case, a portion of the earlier experiment is confirmed. The new results are then viewed from a different perspective in order to evaluate the usefulness of this simulation strategy.

  10. A Comparison of General Case In Vivo and General Case Simulation Plus In Vivo Training.

    ERIC Educational Resources Information Center

    McDonnell, John J.; Ferguson, Brad

    1988-01-01

    The study examined the relative effectiveness and efficiency of general case in vivo and general case simulation plus in vivo training in teaching six students with moderate and severe disabilities to purchase items in fast-food restaurants. Although both strategies led to reliable performance in nontrained settings, the in vivo instruction…

  11. Sampling of general correlators in worm-algorithm based simulations

    NASA Astrophysics Data System (ADS)

    Rindlisbacher, Tobias; Åkerlund, Oscar; de Forcrand, Philippe

    2016-08-01

    Using the complex ϕ4-model as a prototype for a system which is simulated by a worm algorithm, we show that not only the charged correlator <ϕ* (x) ϕ (y) >, but also more general correlators such as < | ϕ (x) | | ϕ (y) | > or < arg ⁡ (ϕ (x)) arg ⁡ (ϕ (y)) >, as well as condensates like < | ϕ | >, can be measured at every step of the Monte Carlo evolution of the worm instead of on closed-worm configurations only. The method generalizes straightforwardly to other systems simulated by worms, such as spin or sigma models.

  12. The development of an interim generalized gate logic software simulator

    NASA Technical Reports Server (NTRS)

    Mcgough, J. G.; Nemeroff, S.

    1985-01-01

    A proof-of-concept computer program called IGGLOSS (Interim Generalized Gate Logic Software Simulator) was developed and is discussed. The simulator engine was designed to perform stochastic estimation of self test coverage (fault-detection latency times) of digital computers or systems. A major attribute of the IGGLOSS is its high-speed simulation: 9.5 x 1,000,000 gates/cpu sec for nonfaulted circuits and 4.4 x 1,000,000 gates/cpu sec for faulted circuits on a VAX 11/780 host computer.

  13. Simulations of accretion disks in pseudo-complex General Relativity

    NASA Astrophysics Data System (ADS)

    Hess, P. O.; Algalán B., M.; Schönenbach, T.; Greiner, W.

    2015-11-01

    After a summary on pseudo-complex General Relativity (pc-GR), circular orbits and stable orbits in general are discussed, including predictions compared to observations. Using a modified version of a model for accretions disks, presented by Page and Thorne in 1974, we apply the raytracing technique in order to simulate the appearance of an accretion disk as it should be observed in a detector. In pc-GR we predict a dark ring near a very massive, rapidly rotating object.

  14. The power of event-driven analytics in Large Scale Data Processing

    ScienceCinema

    None

    2011-04-25

    FeedZai is a software company specialized in creating high-­-throughput low-­-latency data processing solutions. FeedZai develops a product called "FeedZai Pulse" for continuous event-­-driven analytics that makes application development easier for end users. It automatically calculates key performance indicators and baselines, showing how current performance differ from previous history, creating timely business intelligence updated to the second. The tool does predictive analytics and trend analysis, displaying data on real-­-time web-­-based graphics. In 2010 FeedZai won the European EBN Smart Entrepreneurship Competition, in the Digital Models category, being considered one of the "top-­-20 smart companies in Europe". The main objective of this seminar/workshop is to explore the topic for large-­-scale data processing using Complex Event Processing and, in particular, the possible uses of Pulse in the scope of the data processing needs of CERN. Pulse is available as open-­-source and can be licensed both for non-­-commercial and commercial applications. FeedZai is interested in exploring possible synergies with CERN in high-­-volume low-­-latency data processing applications. The seminar will be structured in two sessions, the first one being aimed to expose the general scope of FeedZai's activities, and the second focused on Pulse itself: 10:00-11:00 FeedZai and Large Scale Data Processing Introduction to FeedZai FeedZai Pulse and Complex Event Processing Demonstration Use-Cases and Applications Conclusion and Q&A 11:00-11:15 Coffee break 11:15-12:30 FeedZai Pulse Under the Hood A First FeedZai Pulse Application PulseQL overview Defining KPIs and Baselines Conclusion and Q&A About the speakers Nuno Sebastião is the CEO of FeedZai. Having worked for many years for the European Space Agency (ESA), he was responsible the overall design and development of Satellite Simulation Infrastructure of the agency. Having left ESA to found FeedZai, Nuno is

  15. The power of event-driven analytics in Large Scale Data Processing

    SciTech Connect

    2011-02-24

    FeedZai is a software company specialized in creating high-­-throughput low-­-latency data processing solutions. FeedZai develops a product called "FeedZai Pulse" for continuous event-­-driven analytics that makes application development easier for end users. It automatically calculates key performance indicators and baselines, showing how current performance differ from previous history, creating timely business intelligence updated to the second. The tool does predictive analytics and trend analysis, displaying data on real-­-time web-­-based graphics. In 2010 FeedZai won the European EBN Smart Entrepreneurship Competition, in the Digital Models category, being considered one of the "top-­-20 smart companies in Europe". The main objective of this seminar/workshop is to explore the topic for large-­-scale data processing using Complex Event Processing and, in particular, the possible uses of Pulse in the scope of the data processing needs of CERN. Pulse is available as open-­-source and can be licensed both for non-­-commercial and commercial applications. FeedZai is interested in exploring possible synergies with CERN in high-­-volume low-­-latency data processing applications. The seminar will be structured in two sessions, the first one being aimed to expose the general scope of FeedZai's activities, and the second focused on Pulse itself: 10:00-11:00 FeedZai and Large Scale Data Processing Introduction to FeedZai FeedZai Pulse and Complex Event Processing Demonstration Use-Cases and Applications Conclusion and Q&A 11:00-11:15 Coffee break 11:15-12:30 FeedZai Pulse Under the Hood A First FeedZai Pulse Application PulseQL overview Defining KPIs and Baselines Conclusion and Q&A About the speakers Nuno Sebastião is the CEO of FeedZai. Having worked for many years for the European Space Agency (ESA), he was responsible the overall design and development of Satellite Simulation Infrastructure of the agency. Having left ESA to found FeedZai, Nuno is

  16. The architecture of Newton, a general-purpose dynamics simulator

    NASA Technical Reports Server (NTRS)

    Cremer, James F.; Stewart, A. James

    1989-01-01

    The architecture for Newton, a general-purpose system for simulating the dynamics of complex physical objects, is described. The system automatically formulates and analyzes equations of motion, and performs automatic modification of this system equations when necessitated by changes in kinematic relationships between objects. Impact and temporary contact are handled, although only using simple models. User-directed influence of simulations is achieved using Newton's module, which can be used to experiment with the control of many-degree-of-freedom articulated objects.

  17. Applications of a general thermal/hydraulic simulation tool

    NASA Technical Reports Server (NTRS)

    Cullimore, B. A.

    1989-01-01

    The analytic techniques, sample applications, and development status of a general-purpose computer program called SINDA '85/FLUINT (for systems improved numerical differencing analyzer, 1985 version with fluid integrator), designed for simulating thermal structures and internal fluid systems, are described, with special attention given to the applications of the fluid system capabilities. The underlying assumptions, methodologies, and modeling capabilities of the system are discussed. Sample applications include component-level and system-level simulations. A system-level analysis of a cryogenic storage system is presented.

  18. Generalized Fluid System Simulation Program (GFSSP) Version 6 - General Purpose Thermo-Fluid Network Analysis Software

    NASA Technical Reports Server (NTRS)

    Majumdar, Alok; Leclair, Andre; Moore, Ric; Schallhorn, Paul

    2011-01-01

    GFSSP stands for Generalized Fluid System Simulation Program. It is a general-purpose computer program to compute pressure, temperature and flow distribution in a flow network. GFSSP calculates pressure, temperature, and concentrations at nodes and calculates flow rates through branches. It was primarily developed to analyze Internal Flow Analysis of a Turbopump Transient Flow Analysis of a Propulsion System. GFSSP development started in 1994 with an objective to provide a generalized and easy to use flow analysis tool for thermo-fluid systems.

  19. A General Relativistic Magnetohydrodynamic Simulation of Jet Formation

    NASA Astrophysics Data System (ADS)

    Nishikawa, K.-I.; Richardson, G.; Koide, S.; Shibata, K.; Kudoh, T.; Hardee, P.; Fishman, G. J.

    2005-05-01

    We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity ~0.3c) is created, as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The three-dimensional simulation ran over 100 light crossing time units (τS=rS/c, where rS≡2GM/c2), which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted owing in part to magnetic pressure from the twisting of the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r=3rS. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outward with a wider angle than the initial jet. The widening of the jet is consistent with the outward-moving torsional Alfvén waves. This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk because of the lack of streaming material from an accompanying star.

  20. A General Mission Independent Simulator (GMIS) and Simulator Control Program (SCP)

    NASA Technical Reports Server (NTRS)

    Baker, Paul L.; Moore, J. Michael; Rosenberger, John

    1994-01-01

    GMIS is a general-purpose simulator for testing ground system software. GMIS can be adapted to any mission to simulate changes in the data state maintained by the mission's computers. GMIS was developed in Code 522 NASA Goddard Space Flight Center. The acronym GMIS stands for GOTT Mission Independent Simulator, where GOTT is the Ground Operations Technology Testbed. Within GOTT, GMIS is used to provide simulated data to an installation of TPOCC - the Transportable Payload Operations Control Center. TPOCC was developed by Code 510 as a reusable control center. GOTT uses GMIS and TPOCC to test new technology and new operator procedures.

  1. Generalized directed loop method for quantum Monte Carlo simulations.

    PubMed

    Alet, Fabien; Wessel, Stefan; Troyer, Matthias

    2005-03-01

    Efficient quantum Monte Carlo update schemes called directed loops have recently been proposed, which improve the efficiency of simulations of quantum lattice models. We propose to generalize the detailed balance equations at the local level during the loop construction by accounting for the matrix elements of the operators associated with open world-line segments. Using linear programming techniques to solve the generalized equations, we look for optimal construction schemes for directed loops. This also allows for an extension of the directed loop scheme to general lattice models, such as high-spin or bosonic models. The resulting algorithms are bounce free in larger regions of parameter space than the original directed loop algorithm. The generalized directed loop method is applied to the magnetization process of spin chains in order to compare its efficiency to that of previous directed loop schemes. In contrast to general expectations, we find that minimizing bounces alone does not always lead to more efficient algorithms in terms of autocorrelations of physical observables, because of the nonuniqueness of the bounce-free solutions. We therefore propose different general strategies to further minimize autocorrelations, which can be used as supplementary requirements in any directed loop scheme. We show by calculating autocorrelation times for different observables that such strategies indeed lead to improved efficiency; however, we find that the optimal strategy depends not only on the model parameters but also on the observable of interest. PMID:15903632

  2. A High-Speed, Event-Driven, Active Pixel Sensor Readout for Photon-Counting Microchannel Plate Detectors

    NASA Technical Reports Server (NTRS)

    Kimble, Randy A.; Pain, Bedabrata; Norton, Timothy J.; Haas, J. Patrick; Oegerle, William R. (Technical Monitor)

    2002-01-01

    Silicon array readouts for microchannel plate intensifiers offer several attractive features. In this class of detector, the electron cloud output of the MCP intensifier is converted to visible light by a phosphor; that light is then fiber-optically coupled to the silicon array. In photon-counting mode, the resulting light splashes on the silicon array are recognized and centroided to fractional pixel accuracy by off-chip electronics. This process can result in very high (MCP-limited) spatial resolution while operating at a modest MCP gain (desirable for dynamic range and long term stability). The principal limitation of intensified CCD systems of this type is their severely limited local dynamic range, as accurate photon counting is achieved only if there are not overlapping event splashes within the frame time of the device. This problem can be ameliorated somewhat by processing events only in pre-selected windows of interest of by using an addressable charge injection device (CID) for the readout array. We are currently pursuing the development of an intriguing alternative readout concept based on using an event-driven CMOS Active Pixel Sensor. APS technology permits the incorporation of discriminator circuitry within each pixel. When coupled with suitable CMOS logic outside the array area, the discriminator circuitry can be used to trigger the readout of small sub-array windows only when and where an event splash has been detected, completely eliminating the local dynamic range problem, while achieving a high global count rate capability and maintaining high spatial resolution. We elaborate on this concept and present our progress toward implementing an event-driven APS readout.

  3. A High-Speed, Event-Driven, Active Pixel Sensor Readout for Photon-Counting Microchannel Plate Detectors

    NASA Technical Reports Server (NTRS)

    Kimble, Randy A.; Pain, B.; Norton, T. J.; Haas, P.; Fisher, Richard R. (Technical Monitor)

    2001-01-01

    Silicon array readouts for microchannel plate intensifiers offer several attractive features. In this class of detector, the electron cloud output of the MCP intensifier is converted to visible light by a phosphor; that light is then fiber-optically coupled to the silicon array. In photon-counting mode, the resulting light splashes on the silicon array are recognized and centroided to fractional pixel accuracy by off-chip electronics. This process can result in very high (MCP-limited) spatial resolution for the readout while operating at a modest MCP gain (desirable for dynamic range and long term stability). The principal limitation of intensified CCD systems of this type is their severely limited local dynamic range, as accurate photon counting is achieved only if there are not overlapping event splashes within the frame time of the device. This problem can be ameliorated somewhat by processing events only in pre-selected windows of interest or by using an addressable charge injection device (CID) for the readout array. We are currently pursuing the development of an intriguing alternative readout concept based on using an event-driven CMOS Active Pixel Sensor. APS technology permits the incorporation of discriminator circuitry within each pixel. When coupled with suitable CMOS logic outside the array area, the discriminator circuitry can be used to trigger the readout of small sub-array windows only when and where an event splash has been detected, completely eliminating the local dynamic range problem, while achieving a high global count rate capability and maintaining high spatial resolution. We elaborate on this concept and present our progress toward implementing an event-driven APS readout.

  4. Architecture for event-driven real-time distributed computer systems

    SciTech Connect

    McDonald, J.E.

    1983-01-01

    The author describes a proposed preliminary system design that includes hardware and software for real-time distributed computer systems. This new system is appropriate as a digital avionics architecture or as a real-time multi-computer simulation system using a mixture of computers, mainframes to micros. The hardware contains a network that employs high-speed serial data transmission concepts in emulating a multicomputer shared memory system. The distributed multicomputer system then capitalizes on the attributes of the hardware by structuring the real-time software as the data-driven input-output system. The real-time software executes only on demand and not synchronously as in conventional real-time systems. Background information concerning multi-computer systems using serial and parallel data transmission networks is given. This information supports the design rationale of the proposed hardware system which is basically a technology blend of conventional serial and parallel transmission schemes. 2 references.

  5. Comparison of Cenozoic atmospheric general circulation model simulations

    SciTech Connect

    Barron, E.J.

    1985-01-01

    Paleocene, Eocene, Miocene and present day (with polar ice) geography are specified as the lower boundary condition in a mean annual, energy balance ocean version of the Community Climate Model (CCM), a spectral General Circulation Model of the Atmosphere developed at the National Center for Atmospheric Research. This version of the CCM has a 4.5/sup 0/ latitudinal and 7.5/sup 0/ longitudinal resolution with 9 vertical levels and includes predictions for pressure, winds, temperature, evaporation, precipitation, cloud cover, snow cover and sea ice. The model simulations indicate little geographically-induced climates changes from the Paleocene to the Miocene, but substantial differences between the Miocene and the present simulations. The simulated climate differences between the Miocene and present day include: 1) cooler present temperatures (2/sup 0/C in tropics, 15-35 C in polar latitudes) with the exception of warmer subtropical desert conditions, 2) a generally weaker present hydrologic cycle, with greater subtropical aridity, 3) strengthened present day westerly jets with a slight poleward displacement, and 4) the largest regional climate changes associated with Antarctica. The results of the climate model sensitivity experiments have considerable implications for understanding how geography influences climate.

  6. Automatic CT simulation optimization for radiation therapy: A general strategy

    SciTech Connect

    Li, Hua Chen, Hsin-Chen; Tan, Jun; Gay, Hiram; Michalski, Jeff M.; Mutic, Sasa; Yu, Lifeng; Anastasio, Mark A.; Low, Daniel A.

    2014-03-15

    Purpose: In radiation therapy, x-ray computed tomography (CT) simulation protocol specifications should be driven by the treatment planning requirements in lieu of duplicating diagnostic CT screening protocols. The purpose of this study was to develop a general strategy that allows for automatically, prospectively, and objectively determining the optimal patient-specific CT simulation protocols based on radiation-therapy goals, namely, maintenance of contouring quality and integrity while minimizing patient CT simulation dose. Methods: The authors proposed a general prediction strategy that provides automatic optimal CT simulation protocol selection as a function of patient size and treatment planning task. The optimal protocol is the one that delivers the minimum dose required to provide a CT simulation scan that yields accurate contours. Accurate treatment plans depend on accurate contours in order to conform the dose to actual tumor and normal organ positions. An image quality index, defined to characterize how simulation scan quality affects contour delineation, was developed and used to benchmark the contouring accuracy and treatment plan quality within the predication strategy. A clinical workflow was developed to select the optimal CT simulation protocols incorporating patient size, target delineation, and radiation dose efficiency. An experimental study using an anthropomorphic pelvis phantom with added-bolus layers was used to demonstrate how the proposed prediction strategy could be implemented and how the optimal CT simulation protocols could be selected for prostate cancer patients based on patient size and treatment planning task. Clinical IMRT prostate treatment plans for seven CT scans with varied image quality indices were separately optimized and compared to verify the trace of target and organ dosimetry coverage. Results: Based on the phantom study, the optimal image quality index for accurate manual prostate contouring was 4.4. The optimal tube

  7. A General Relativistic Magnetohydrodynamic Simulation of Jet Formation

    NASA Technical Reports Server (NTRS)

    Nishikawa, K.-I.; Richardson, G.; Koide, S.; Shibata, K.; Kudoh, T.; Hardee, P.; Fishman, G. J.

    2005-01-01

    We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation ofjet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity approx.0.3c) is created, as shown by previous two-dimensional axi- symmetric simulations with mirror symmetry at the equator. The three-dimensional simulation ran over 100 light crossing time units (T(sub s) = r(sub s)/c, where r(sub s = 2GM/c(sup 2), which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted owing in part to magnetic pressure from the twisting of the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r(sub s). At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface ofthe thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outward with a wider angle than the initial jet. The widening of the jet is consistent with the outward-moving torsional Alfven waves. This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk because of the iack of streaming materiai from an accompanying star.

  8. Localized and generalized simulated wear of resin composites.

    PubMed

    Barkmeier, W W; Takamizawa, T; Erickson, R L; Tsujimoto, A; Latta, M; Miyazaki, M

    2015-01-01

    A laboratory study was conducted to examine the wear of resin composite materials using both a localized and generalized wear simulation model. Twenty specimens each of seven resin composites (Esthet•X HD [HD], Filtek Supreme Ultra [SU], Herculite Ultra [HU], SonicFill [SF], Tetric EvoCeram Bulk Fill [TB], Venus Diamond [VD], and Z100 Restorative [Z]) were subjected to a wear challenge of 400,000 cycles for both localized and generalized wear in a Leinfelder-Suzuki wear simulator (Alabama machine). The materials were placed in custom cylinder-shaped stainless steel fixtures. A stainless steel ball bearing (r=2.387 mm) was used as the antagonist for localized wear, and a stainless steel, cylindrical antagonist with a flat tip was used for generalized wear. A water slurry of polymethylmethacrylate (PMMA) beads was used as the abrasive media. A noncontact profilometer (Proscan 2100) with Proscan software was used to digitize the surface contours of the pretest and posttest specimens. AnSur 3D software was used for wear assessment. For localized testing, maximum facet depth (μm) and volume loss (mm(3)) were used to compare the materials. The mean depth of the facet surface (μm) and volume loss (mm(3)) were used for comparison of the generalized wear specimens. A one-way analysis of variance (ANOVA) and Tukey post hoc test were used for data analysis of volume loss for both localized and generalized wear, maximum facet depth for localized wear, and mean depth of the facet for generalized wear. The results for localized wear simulation were as follows [mean (standard deviation)]: maximum facet depth (μm)--Z, 59.5 (14.7); HU, 99.3 (16.3); SU, 102.8 (13.8); HD, 110.2 (13.3); VD, 114.0 (10.3); TB, 125.5 (12.1); SF, 195.9 (16.9); volume loss (mm(3))--Z, 0.013 (0.002); SU, 0.026 (0.006); HU, 0.043 (0.008); VD, 0.057 (0.009); HD, 0.058 (0.014); TB, 0.061 (0.010); SF, 0.135 (0.024). Generalized wear simulation results were as follows: mean depth of facet (μm)--Z, 9.3 (3

  9. The Speedster-EXD- A New Event-Driven Hybrid CMOS X-ray Detector

    NASA Astrophysics Data System (ADS)

    Griffith, Christopher V.; Falcone, Abraham D.; Prieskorn, Zachary R.; Burrows, David N.

    2016-01-01

    The Speedster-EXD is a new 64×64 pixel, 40-μm pixel pitch, 100-μm depletion depth hybrid CMOS x-ray detector with the capability of reading out only those pixels containing event charge, thus enabling fast effective frame rates. A global charge threshold can be specified, and pixels containing charge above this threshold are flagged and read out. The Speedster detector has also been designed with other advanced in-pixel features to improve performance, including a low-noise, high-gain capacitive transimpedance amplifier that eliminates interpixel capacitance crosstalk (IPC), and in-pixel correlated double sampling subtraction to reduce reset noise. We measure the best energy resolution on the Speedster-EXD detector to be 206 eV (3.5%) at 5.89 keV and 172 eV (10.0%) at 1.49 keV. The average IPC to the four adjacent pixels is measured to be 0.25%±0.2% (i.e., consistent with zero). The pixel-to-pixel gain variation is measured to be 0.80%±0.03%, and a Monte Carlo simulation is applied to better characterize the contributions to the energy resolution.

  10. Event-driven Monte Carlo: Exact dynamics at all time scales for discrete-variable models

    NASA Astrophysics Data System (ADS)

    Mendoza-Coto, Alejandro; Díaz-Méndez, Rogelio; Pupillo, Guido

    2016-06-01

    We present an algorithm for the simulation of the exact real-time dynamics of classical many-body systems with discrete energy levels. In the same spirit of kinetic Monte Carlo methods, a stochastic solution of the master equation is found, with no need to define any other phase-space construction. However, unlike existing methods, the present algorithm does not assume any particular statistical distribution to perform moves or to advance the time, and thus is a unique tool for the numerical exploration of fast and ultra-fast dynamical regimes. By decomposing the problem in a set of two-level subsystems, we find a natural variable step size, that is well defined from the normalization condition of the transition probabilities between the levels. We successfully test the algorithm with known exact solutions for non-equilibrium dynamics and equilibrium thermodynamical properties of Ising-spin models in one and two dimensions, and compare to standard implementations of kinetic Monte Carlo methods. The present algorithm is directly applicable to the study of the real-time dynamics of a large class of classical Markovian chains, and particularly to short-time situations where the exact evolution is relevant.

  11. Solute transport processes in flow-event-driven stream-aquifer interaction

    NASA Astrophysics Data System (ADS)

    Xie, Yueqing; Cook, Peter G.; Simmons, Craig T.

    2016-07-01

    The interaction between streams and groundwater controls key features of the stream hydrograph and chemograph. Since surface runoff is usually less saline than groundwater, flow events are usually accompanied by declines in stream salinity. In this paper, we use numerical modelling to show that, at any particular monitoring location: (i) the increase in stream stage associated with a flow event will precede the decrease in solute concentration (arrival time lag for solutes); and (ii) the decrease in stream stage following the flow peak will usually precede the subsequent return (increase) in solute concentration (return time lag). Both arrival time lag and return time lag increase with increasing wave duration. However, arrival time lag decreases with increasing wave amplitude, whereas return time lag increases. Furthermore, while arrival time lag is most sensitive to parameters that control river velocity (channel roughness and stream slope), return time lag is most sensitive to groundwater parameters (aquifer hydraulic conductivity, recharge rate, and dispersitivity). Additionally, the absolute magnitude of the decrease in river concentration is sensitive to both river and groundwater parameters. Our simulations also show that in-stream mixing is dominated by wave propagation and bank storage processes, and in-stream dispersion has a relatively minor effect on solute concentrations. This has important implications for spreading of contaminants released to streams. Our work also demonstrates that a high contribution of pre-event water (or groundwater) within the flow hydrograph can be caused by the combination of in-stream and bank storage exchange processes, and does not require transport of pre-event water through the catchment.

  12. An event-driven approach for studying gene block evolution in bacteria

    PubMed Central

    Ream, David C.; Bankapur, Asma R.; Friedberg, Iddo

    2015-01-01

    Motivation: Gene blocks are genes co-located on the chromosome. In many cases, gene blocks are conserved between bacterial species, sometimes as operons, when genes are co-transcribed. The conservation is rarely absolute: gene loss, gain, duplication, block splitting and block fusion are frequently observed. An open question in bacterial molecular evolution is that of the formation and breakup of gene blocks, for which several models have been proposed. These models, however, are not generally applicable to all types of gene blocks, and consequently cannot be used to broadly compare and study gene block evolution. To address this problem, we introduce an event-based method for tracking gene block evolution in bacteria. Results: We show here that the evolution of gene blocks in proteobacteria can be described by a small set of events. Those include the insertion of genes into, or the splitting of genes out of a gene block, gene loss, and gene duplication. We show how the event-based method of gene block evolution allows us to determine the evolutionary rateand may be used to trace the ancestral states of their formation. We conclude that the event-based method can be used to help us understand the formation of these important bacterial genomic structures. Availability and implementation: The software is available under GPLv3 license on http://github.com/reamdc1/gene_block_evolution.git. Supplementary online material: http://iddo-friedberg.net/operon-evolution Contact: i.friedberg@miamioh.edu Supplementary information: Supplementary data are available at Bioinformatics online. PMID:25717195

  13. Characterization and development of an event-driven hybrid CMOS x-ray detector

    NASA Astrophysics Data System (ADS)

    Griffith, Christopher

    2015-06-01

    Hybrid CMOS detectors (HCD) have provided great benefit to the infrared and optical fields of astronomy, and they are poised to do the same for X-ray astronomy. Infrared HCDs have already flown on the Hubble Space Telescope and the Wide-Field Infrared Survey Explorer (WISE) mission and are slated to fly on the James Webb Space Telescope (JWST). Hybrid CMOS X-ray detectors offer low susceptibility to radiation damage, low power consumption, and fast readout time to avoid pile-up. The fast readout time is necessary for future high throughput X-ray missions. The Speedster-EXD X-ray HCD presented in this dissertation offers new in-pixel features and reduces known noise sources seen on previous generation HCDs. The Speedster-EXD detector makes a great step forward in the development of these detectors for future space missions. This dissertation begins with an overview of future X-ray space mission concepts and their detector requirements. The background on the physics of semiconductor devices and an explanation of the detection of X-rays with these devices will be discussed followed by a discussion on CCDs and CMOS detectors. Next, hybrid CMOS X-ray detectors will be explained including their advantages and disadvantages. The Speedster-EXD detector and its new features will be outlined including its ability to only read out pixels which contain X-ray events. Test stand design and construction for the Speedster-EXD detector is outlined and the characterization of each parameter on two Speedster-EXD detectors is detailed including read noise, dark current, interpixel capacitance crosstalk (IPC), and energy resolution. Gain variation is also characterized, and a Monte Carlo simulation of its impact on energy resolution is described. This analysis shows that its effect can be successfully nullified with proper calibration, which would be important for a flight mission. Appendix B contains a study of the extreme tidal disruption event, Swift J1644+57, to search for

  14. A Distributed Laboratory for Event-Driven Coastal Prediction and Hazard Planning

    NASA Astrophysics Data System (ADS)

    Bogden, P.; Allen, G.; MacLaren, J.; Creager, G. J.; Flournoy, L.; Sheng, Y. P.; Graber, H.; Graves, S.; Conover, H.; Luettich, R.; Perrie, W.; Ramakrishnan, L.; Reed, D. A.; Wang, H. V.

    2006-12-01

    The 2005 Atlantic hurricane season was the most active in recorded history. Collectively, 2005 hurricanes caused more than 2,280 deaths and record damages of over 100 billion dollars. Of the storms that made landfall, Dennis, Emily, Katrina, Rita, and Wilma caused most of the destruction. Accurate predictions of storm-driven surge, wave height, and inundation can save lives and help keep recovery costs down, provided the information gets to emergency response managers in time. The information must be available well in advance of landfall so that responders can weigh the costs of unnecessary evacuation against the costs of inadequate preparation. The SURA Coastal Ocean Observing and Prediction (SCOOP) Program is a multi-institution collaboration implementing a modular, distributed service-oriented architecture for real time prediction and visualization of the impacts of extreme atmospheric events. The modular infrastructure enables real-time prediction of multi- scale, multi-model, dynamic, data-driven applications. SURA institutions are working together to create a virtual and distributed laboratory integrating coastal models, simulation data, and observations with computational resources and high speed networks. The loosely coupled architecture allows teams of computer and coastal scientists at multiple institutions to innovate complex system components that are interconnected with relatively stable interfaces. The operational system standardizes at the interface level to enable substantial innovation by complementary communities of coastal and computer scientists. This architectural philosophy solves a long-standing problem associated with the transition from research to operations. The SCOOP Program thereby implements a prototype laboratory consistent with the vision of a national, multi-agency initiative called the Integrated Ocean Observing System (IOOS). Several service- oriented components of the SCOOP enterprise architecture have already been designed and

  15. Data Albums: An Event Driven Search, Aggregation and Curation Tool for Earth Science

    NASA Technical Reports Server (NTRS)

    Ramachandran, Rahul; Kulkarni, Ajinkya; Maskey, Manil; Bakare, Rohan; Basyal, Sabin; Li, Xiang; Flynn, Shannon

    2014-01-01

    One of the largest continuing challenges in any Earth science investigation is the discovery and access of useful science content from the increasingly large volumes of Earth science data and related information available. Approaches used in Earth science research such as case study analysis and climatology studies involve gathering discovering and gathering diverse data sets and information to support the research goals. Research based on case studies involves a detailed description of specific weather events using data from different sources, to characterize physical processes in play for a specific event. Climatology-based research tends to focus on the representativeness of a given event, by studying the characteristics and distribution of a large number of events. This allows researchers to generalize characteristics such as spatio-temporal distribution, intensity, annual cycle, duration, etc. To gather relevant data and information for case studies and climatology analysis is both tedious and time consuming. Current Earth science data systems are designed with the assumption that researchers access data primarily by instrument or geophysical parameter. Those who know exactly the datasets of interest can obtain the specific files they need using these systems. However, in cases where researchers are interested in studying a significant event, they have to manually assemble a variety of datasets relevant to it by searching the different distributed data systems. In these cases, a search process needs to be organized around the event rather than observing instruments. In addition, the existing data systems assume users have sufficient knowledge regarding the domain vocabulary to be able to effectively utilize their catalogs. These systems do not support new or interdisciplinary researchers who may be unfamiliar with the domain terminology. This paper presents a specialized search, aggregation and curation tool for Earth science to address these existing

  16. A generalized Poisson solver for first-principles device simulations

    NASA Astrophysics Data System (ADS)

    Bani-Hashemian, Mohammad Hossein; Brück, Sascha; Luisier, Mathieu; VandeVondele, Joost

    2016-01-01

    Electronic structure calculations of atomistic systems based on density functional theory involve solving the Poisson equation. In this paper, we present a plane-wave based algorithm for solving the generalized Poisson equation subject to periodic or homogeneous Neumann conditions on the boundaries of the simulation cell and Dirichlet type conditions imposed at arbitrary subdomains. In this way, source, drain, and gate voltages can be imposed across atomistic models of electronic devices. Dirichlet conditions are enforced as constraints in a variational framework giving rise to a saddle point problem. The resulting system of equations is then solved using a stationary iterative method in which the generalized Poisson operator is preconditioned with the standard Laplace operator. The solver can make use of any sufficiently smooth function modelling the dielectric constant, including density dependent dielectric continuum models. For all the boundary conditions, consistent derivatives are available and molecular dynamics simulations can be performed. The convergence behaviour of the scheme is investigated and its capabilities are demonstrated.

  17. A generalized Poisson solver for first-principles device simulations.

    PubMed

    Bani-Hashemian, Mohammad Hossein; Brück, Sascha; Luisier, Mathieu; VandeVondele, Joost

    2016-01-28

    Electronic structure calculations of atomistic systems based on density functional theory involve solving the Poisson equation. In this paper, we present a plane-wave based algorithm for solving the generalized Poisson equation subject to periodic or homogeneous Neumann conditions on the boundaries of the simulation cell and Dirichlet type conditions imposed at arbitrary subdomains. In this way, source, drain, and gate voltages can be imposed across atomistic models of electronic devices. Dirichlet conditions are enforced as constraints in a variational framework giving rise to a saddle point problem. The resulting system of equations is then solved using a stationary iterative method in which the generalized Poisson operator is preconditioned with the standard Laplace operator. The solver can make use of any sufficiently smooth function modelling the dielectric constant, including density dependent dielectric continuum models. For all the boundary conditions, consistent derivatives are available and molecular dynamics simulations can be performed. The convergence behaviour of the scheme is investigated and its capabilities are demonstrated. PMID:26827208

  18. Better Space Construction Decisions by Instructional Program Simulation Utilizing the Generalized Academic Simulation Programs.

    ERIC Educational Resources Information Center

    Apker, Wesley

    This school district utilized the generalized academic simulation programs (GASP) to assist in making decisions regarding the kinds of facilities that should be constructed at Pilchuck Senior High School. Modular scheduling was one of the basic educational parameters used in determining the number and type of facilities needed. The objectives of…

  19. A General Simulator for Reaction-Based Biogeochemical Processes

    SciTech Connect

    Fang, Yilin; Yabusaki, Steven B.; Yeh, George

    2006-02-01

    As more complex biogeochemical situations are being investigated (e.g., evolving reactivity, passivation of reactive surfaces, dissolution of sorbates), there is a growing need for biogeochemical simulators to flexibly and facilely address new reaction forms and rate laws. This paper presents an approach that accommodates this need to efficiently simulate general biogeochemical processes, while insulating the user from additional code development. The approach allows for the automatic extraction of fundamental reaction stoichiometry and thermodynamics from a standard chemistry database, and the symbolic entry of arbitrarily complex user-specified reaction forms, rate laws, and equilibria. The user-specified equilibrium and kinetic reactions (i.e., reactions not defined in the format of the standardized database) are interpreted by the Maple symbolic mathematical software package. FORTRAN 90 code is then generated by Maple for (1) the analytical Jacobian matrix (if preferred over the numerical Jacobian matrix) used in the Newton-Raphson solution procedure, and (2) the residual functions for user-specified equilibrium expressions and rate laws. Matrix diagonalization eliminates the need to conceptualize the system of reactions as a tableau, while identifying a minimum rank set of basis species with enhanced numerical convergence properties. The newly generated code, which is designed to operate in the BIOGEOCHEM biogeochemical simulator, is then compiled and linked into the BIOGEOCHEM executable. With these features, users can avoid recoding the simulator to accept new equilibrium expressions or kinetic rate laws, while still taking full advantage of the stoichiometry and thermodynamics provided by an existing chemical database. Thus, the approach introduces efficiencies in the specification of biogeochemical reaction networks and eliminates opportunities for mistakes in preparing input files and coding errors. Test problems are used to demonstrate the features of

  20. Event-driven, pattern-based methodology for cost-effective development of standardized personal health devices.

    PubMed

    Martínez-Espronceda, Miguel; Trigo, Jesús D; Led, Santiago; Barrón-González, H Gilberto; Redondo, Javier; Baquero, Alfonso; Serrano, Luis

    2014-11-01

    Experiences applying standards in personal health devices (PHDs) show an inherent trade-off between interoperability and costs (in terms of processing load and development time). Therefore, reducing hardware and software costs as well as time-to-market is crucial for standards adoption. The ISO/IEEE11073 PHD family of standards (also referred to as X73PHD) provides interoperable communication between PHDs and aggregators. Nevertheless, the responsibility of achieving inexpensive implementations of X73PHD in limited resource microcontrollers falls directly on the developer. Hence, the authors previously presented a methodology based on patterns to implement X73-compliant PHDs into devices with low-voltage low-power constraints. That version was based on multitasking, which required additional features and resources. This paper therefore presents an event-driven evolution of the patterns-based methodology for cost-effective development of standardized PHDs. The results of comparing between the two versions showed that the mean values of decrease in memory consumption and cycles of latency are 11.59% and 45.95%, respectively. In addition, several enhancements in terms of cost-effectiveness and development time can be derived from the new version of the methodology. Therefore, the new approach could help in producing cost-effective X73-compliant PHDs, which in turn could foster the adoption of standards. PMID:25123101

  1. Event driven executive

    NASA Technical Reports Server (NTRS)

    Tulpule, Bhalchandra R. (Inventor); Collins, Robert E. (Inventor); Cheetham, John (Inventor); Cornwell, Smith (Inventor)

    1990-01-01

    Tasks may be planned for execution on a single processor or are split up by the designer for execution among a plurality of signal processors. The tasks are modeled using a design aid called a precedence graph, from which a dependency table and a prerequisite table are established for reference within each processor. During execution, at the completion of a given task, an end of task interrupt is provided from any processor which has completed a task to any and all other processors including itself in which completion of that task is a prerequisite for commencement of any dependent tasks. The relevant updated data may be transferred by the processor either before or after signalling task completion to the processors needing the updated data prior to commencing execution of the dependent tasks. Coherency may be ensured, however, by sending the data before the interrupt. When the end of task interrupt is received in a processor, its dependency table is consulted to determine those tasks dependent upon completion of the task which has just been signalled as completed, and task dependency signals indicative thereof are provided and stored in a current status list of a prerequisite table. The current status of all current prerequisites are compared to the complete prerequisites listed for all affected tasks and those tasks for which the comparison indicates that all prerequisites have been met are queued for execution in a selected order.

  2. MPSim: A Massively Parallel General Simulation Program for Materials

    NASA Astrophysics Data System (ADS)

    Iotov, Mihail; Gao, Guanghua; Vaidehi, Nagarajan; Cagin, Tahir; Goddard, William A., III

    1997-08-01

    In this talk, we describe a general purpose Massively Parallel Simulation (MPSim) program used for computational materials science and life sciences. We also will present scaling aspects of the program along with several case studies. The program incorporates highly efficient CMM method to accurately calculate the interactions. For studying bulk materials, the program uses the Reduced CMM to account for infinite range sums. The software embodies various advanced molecular dynamics algorithms, energy and structure optimization techniques with a set of analysis tools suitable for large scale structures. The applications using the program range amorphous polymers, liquid-polymer interfaces, large viruses, million atom clusters, surfaces, gas diffusion in polymers. Program is originally developed on KSR in an object oriented fashion and is ported to SGI-PC, and HP-Examplar. Message Passing version is originally implemented on Intel Paragon using NX, then MPI and later tested on Cray T3D, and IBM SP2 platforms.

  3. Magnetohydrodynamical simulations of a deep tidal disruption in general relativity

    NASA Astrophysics Data System (ADS)

    Sądowski, Aleksander; Tejeda, Emilio; Gafton, Emanuel; Rosswog, Stephan; Abarca, David

    2016-06-01

    We perform hydro- and magnetohydrodynamical general-relativistic simulations of a tidal disruption of a 0.1 M⊙ red dwarf approaching a 105 M⊙ non-rotating massive black hole on a close (impact parameter β = 10) elliptical (eccentricity e = 0.97) orbit. We track the debris self-interaction, circularization and the accompanying accretion through the black hole horizon. We find that the relativistic precession leads to the formation of a self-crossing shock. The dissipated kinetic energy heats up the incoming debris and efficiently generates a quasi-spherical outflow. The self-interaction is modulated because of the feedback exerted by the flow on itself. The debris quickly forms a thick, almost marginally bound disc that remains turbulent for many orbital periods. Initially, the accretion through the black hole horizon results from the self-interaction, while in the later stages it is dominated by the debris originally ejected in the shocked region, as it gradually falls back towards the hole. The effective viscosity in the debris disc stems from the original hydrodynamical turbulence, which dominates over the magnetic component. The radiative efficiency is very low because of low energetics of the gas crossing the horizon and large optical depth that results in photon trapping. Although the parameters of the simulated tidal disruption are probably not representative of most observed events, it is possible to extrapolate some of its properties towards more common configurations.

  4. Amyloid oligomer structure characterization from simulations: A general method

    SciTech Connect

    Nguyen, Phuong H.; Li, Mai Suan

    2014-03-07

    Amyloid oligomers and plaques are composed of multiple chemically identical proteins. Therefore, one of the first fundamental problems in the characterization of structures from simulations is the treatment of the degeneracy, i.e., the permutation of the molecules. Second, the intramolecular and intermolecular degrees of freedom of the various molecules must be taken into account. Currently, the well-known dihedral principal component analysis method only considers the intramolecular degrees of freedom, and other methods employing collective variables can only describe intermolecular degrees of freedom at the global level. With this in mind, we propose a general method that identifies all the structures accurately. The basis idea is that the intramolecular and intermolecular states are described in terms of combinations of single-molecule and double-molecule states, respectively, and the overall structures of oligomers are the product basis of the intramolecular and intermolecular states. This way, the degeneracy is automatically avoided. The method is illustrated on the conformational ensemble of the tetramer of the Alzheimer's peptide Aβ{sub 9−40}, resulting from two atomistic molecular dynamics simulations in explicit solvent, each of 200 ns, starting from two distinct structures.

  5. Generalized Fluid System Simulation Program, Version 6.0

    NASA Technical Reports Server (NTRS)

    Majumdar, A. K.; LeClair, A. C.; Moore, A.; Schallhorn, P. A.

    2013-01-01

    The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependant flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 25 demonstrated example problems.

  6. Generalized Fluid System Simulation Program (GFSSP) - Version 6

    NASA Technical Reports Server (NTRS)

    Majumdar, Alok; LeClair, Andre; Moore, Ric; Schallhorn, Paul

    2015-01-01

    The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, flow control valves and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. Users can introduce new physics, non-linear and time-dependent boundary conditions through user-subroutine.

  7. Generalized Fluid System Simulation Program, Version 5.0-Educational

    NASA Technical Reports Server (NTRS)

    Majumdar, A. K.

    2011-01-01

    The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the point, drag and click method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids and 21 different resistance/source options are provided for modeling momentum sources or sinks in the branches. This Technical Memorandum illustrates the application and verification of the code through 12 demonstrated example problems.

  8. Alternative methods to predict actual evapotranspiration illustrate the importance of accounting for phenology - Part 2: The event driven phenology model

    NASA Astrophysics Data System (ADS)

    Kovalskyy, V.; Henebry, G. M.

    2011-05-01

    Evapotranspiration (ET) flux constitutes a major component of both the water and energy balances at the land surface. Among the many factors that control evapotranspiration, phenology poses a major source of uncertainty in attempts to predict ET. Contemporary approaches to ET modeling and monitoring frequently summarize the complexity of the seasonal development of vegetation cover into static phenological trajectories (or climatologies) that lack sensitivity to changing environmental conditions. The Event Driven Phenology Model (EDPM) offers an alternative, interactive approach to representing phenology. This study presents the results of an experiment designed to illustrate the differences in ET arising from various techniques used to mimic phenology in models of land surface processes. The experiment compares and contrasts two realizations of static phenologies derived from long-term satellite observations of the Normalized Difference Vegetation Index (NDVI) against canopy trajectories produced by the interactive EDPM trained on flux tower observations. The assessment was carried out through validation of predicted ET against records collected by flux tower instruments. The VegET model (Senay, 2008) was used as a framework to estimate daily actual evapotranspiration and supplied with seasonal canopy trajectories produced by the EDPM and traditional techniques. The interactive approach presented the following advantages over phenology modeled with static climatologies: (a) lower prediction bias in crops; (b) smaller root mean square error in daily ET - 0.5 mm per day on average; (c) stable level of errors throughout the season similar among different land cover types and locations; and (d) better estimation of season duration and total seasonal ET.

  9. Alternative methods to predict actual evapotranspiration illustrate the importance of accounting for phenology - Part 2: The event driven phenology model

    NASA Astrophysics Data System (ADS)

    Kovalskyy, V.; Henebry, G. M.

    2012-01-01

    Evapotranspiration (ET) flux constitutes a major component of both the water and energy balances at the land surface. Among the many factors that control evapotranspiration, phenology poses a major source of uncertainty in attempts to predict ET. Contemporary approaches to ET modeling and monitoring frequently summarize the complexity of the seasonal development of vegetation cover into static phenological trajectories (or climatologies) that lack sensitivity to changing environmental conditions. The Event Driven Phenology Model (EDPM) offers an alternative, interactive approach to representing phenology. This study presents the results of an experiment designed to illustrate the differences in ET arising from various techniques used to mimic phenology in models of land surface processes. The experiment compares and contrasts two realizations of static phenologies derived from long-term satellite observations of the Normalized Difference Vegetation Index (NDVI) against canopy trajectories produced by the interactive EDPM trained on flux tower observations. The assessment was carried out through validation of predicted ET against records collected by flux tower instruments. The VegET model (Senay, 2008) was used as a framework to estimate daily actual evapotranspiration and supplied with seasonal canopy trajectories produced by the EDPM and traditional techniques. The interactive approach presented the following advantages over phenology modeled with static climatologies: (a) lower prediction bias in crops; (b) smaller root mean square error in daily ET - 0.5 mm per day on average; (c) stable level of errors throughout the season similar among different land cover types and locations; and (d) better estimation of season duration and total seasonal ET.

  10. Hospitable archean climates simulated by a general circulation model.

    PubMed

    Wolf, E T; Toon, O B

    2013-07-01

    Evidence from ancient sediments indicates that liquid water and primitive life were present during the Archean despite the faint young Sun. To date, studies of Archean climate typically utilize simplified one-dimensional models that ignore clouds and ice. Here, we use an atmospheric general circulation model coupled to a mixed-layer ocean model to simulate the climate circa 2.8 billion years ago when the Sun was 20% dimmer than it is today. Surface properties are assumed to be equal to those of the present day, while ocean heat transport varies as a function of sea ice extent. Present climate is duplicated with 0.06 bar of CO2 or alternatively with 0.02 bar of CO2 and 0.001 bar of CH4. Hot Archean climates, as implied by some isotopic reconstructions of ancient marine cherts, are unattainable even in our warmest simulation having 0.2 bar of CO2 and 0.001 bar of CH4. However, cooler climates with significant polar ice, but still dominated by open ocean, can be maintained with modest greenhouse gas amounts, posing no contradiction with CO2 constraints deduced from paleosols or with practical limitations on CH4 due to the formation of optically thick organic hazes. Our results indicate that a weak version of the faint young Sun paradox, requiring only that some portion of the planet's surface maintain liquid water, may be resolved with moderate greenhouse gas inventories. Thus, hospitable late Archean climates are easily obtained in our climate model. PMID:23808659

  11. Extension of Generalized Fluid System Simulation Program's Fluid Property Database

    NASA Technical Reports Server (NTRS)

    Patel, Kishan

    2011-01-01

    This internship focused on the development of additional capabilities for the General Fluid Systems Simulation Program (GFSSP). GFSSP is a thermo-fluid code used to evaluate system performance by a finite volume-based network analysis method. The program was developed primarily to analyze the complex internal flow of propulsion systems and is capable of solving many problems related to thermodynamics and fluid mechanics. GFSSP is integrated with thermodynamic programs that provide fluid properties for sub-cooled, superheated, and saturation states. For fluids that are not included in the thermodynamic property program, look-up property tables can be provided. The look-up property tables of the current release version can only handle sub-cooled and superheated states. The primary purpose of the internship was to extend the look-up tables to handle saturated states. This involves a) generation of a property table using REFPROP, a thermodynamic property program that is widely used, and b) modifications of the Fortran source code to read in an additional property table containing saturation data for both saturated liquid and saturated vapor states. Also, a method was implemented to calculate the thermodynamic properties of user-fluids within the saturation region, given values of pressure and enthalpy. These additions required new code to be written, and older code had to be adjusted to accommodate the new capabilities. Ultimately, the changes will lead to the incorporation of this new capability in future versions of GFSSP. This paper describes the development and validation of the new capability.

  12. GLoBES: General Long Baseline Experiment Simulator

    NASA Astrophysics Data System (ADS)

    Huber, Patrick; Kopp, Joachim; Lindner, Manfred; Rolinec, Mark; Winter, Walter

    2007-09-01

    GLoBES (General Long Baseline Experiment Simulator) is a flexible software package to simulate neutrino oscillation long baseline and reactor experiments. On the one hand, it contains a comprehensive abstract experiment definition language (AEDL), which allows to describe most classes of long baseline experiments at an abstract level. On the other hand, it provides a C-library to process the experiment information in order to obtain oscillation probabilities, rate vectors, and Δχ-values. Currently, GLoBES is available for GNU/Linux. Since the source code is included, the port to other operating systems is in principle possible. GLoBES is an open source code that has previously been described in Computer Physics Communications 167 (2005) 195 and in Ref. [7]). The source code and a comprehensive User Manual for GLoBES v3.0.8 is now available from the CPC Program Library as described in the Program Summary below. The home of GLobES is http://www.mpi-hd.mpg.de/~globes/. Program summaryProgram title: GLoBES version 3.0.8 Catalogue identifier: ADZI_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZI_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 145 295 No. of bytes in distributed program, including test data, etc.: 1 811 892 Distribution format: tar.gz Programming language: C Computer: GLoBES builds and installs on 32bit and 64bit Linux systems Operating system: 32bit or 64bit Linux RAM: Typically a few MBs Classification: 11.1, 11.7, 11.10 External routines: GSL—The GNU Scientific Library, www.gnu.org/software/gsl/ Nature of problem: Neutrino oscillations are now established as the leading flavor transition mechanism for neutrinos. In a long history of many experiments, see, e.g., [1], two oscillation frequencies have been identified: The fast atmospheric

  13. Sensitivity simulations of superparameterised convection in a general circulation model

    NASA Astrophysics Data System (ADS)

    Rybka, Harald; Tost, Holger

    2015-04-01

    Cloud Resolving Models (CRMs) covering a horizontal grid spacing from a few hundred meters up to a few kilometers have been used to explicitly resolve small-scale and mesoscale processes. Special attention has been paid to realistically represent cloud dynamics and cloud microphysics involving cloud droplets, ice crystals, graupel and aerosols. The entire variety of physical processes on the small-scale interacts with the larger-scale circulation and has to be parameterised on the coarse grid of a general circulation model (GCM). Since more than a decade an approach to connect these two types of models which act on different scales has been developed to resolve cloud processes and their interactions with the large-scale flow. The concept is to use an ensemble of CRM grid cells in a 2D or 3D configuration in each grid cell of the GCM to explicitly represent small-scale processes avoiding the use of convection and large-scale cloud parameterisations which are a major source for uncertainties regarding clouds. The idea is commonly known as superparameterisation or cloud-resolving convection parameterisation. This study presents different simulations of an adapted Earth System Model (ESM) connected to a CRM which acts as a superparameterisation. Simulations have been performed with the ECHAM/MESSy atmospheric chemistry (EMAC) model comparing conventional GCM runs (including convection and large-scale cloud parameterisations) with the improved superparameterised EMAC (SP-EMAC) modeling one year with prescribed sea surface temperatures and sea ice content. The sensitivity of atmospheric temperature, precipiation patterns, cloud amount and types is observed changing the embedded CRM represenation (orientation, width, no. of CRM cells, 2D vs. 3D). Additionally, we also evaluate the radiation balance with the new model configuration, and systematically analyse the impact of tunable parameters on the radiation budget and hydrological cycle. Furthermore, the subgrid

  14. Generalized Fluid System Simulation Program, Version 6.0

    NASA Technical Reports Server (NTRS)

    Majumdar, A. K.; LeClair, A. C.; Moore, R.; Schallhorn, P. A.

    2016-01-01

    The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct

  15. Evaluation of a coupled event-driven phenology and evapotranspiration model for croplands in the United States northern Great Plains

    NASA Astrophysics Data System (ADS)

    Kovalskyy, V.; Henebry, G. M.; Roy, D. P.; Adusei, B.; Hansen, M.; Senay, G.; Mocko, D. M.

    2013-06-01

    A new model coupling scheme with remote sensing data assimilation was developed for estimation of daily actual evapotranspiration (ET). The scheme consists of the VegET, a model to estimate ET from meteorological and water balance data, and an Event Driven Phenology Model (EDPM), an empirical crop specific model trained on multiple years of flux tower data transformed into six types of environmental forcings that are called "events" to emphasize their temporally discrete character, which has advantages for modeling multiple contingent influences. The EDPM in prognostic mode supplies seasonal trajectories of normalized difference vegetation index (NDVI); whereas in diagnostic mode, it can adjust the NDVI prediction with assimilated remotely sensed observations. The scheme was deployed within the croplands of the Northern Great Plains. The evaluation used 2007-2009 land surface forcing data from the North American Land Data Assimilation System and crop maps derived from remotely sensed data of NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). We compared the NDVI produced by the EDPM with NDVI data derived from the MODIS nadir bidirectional reflectance distribution function adjusted reflectance product. The EDPM performance in prognostic mode yielded a coefficient of determination (r2) of 0.8 ± 0.15and the root mean square error (RMSE) of 0.1 ± 0.035 across the entire study area. Retrospective correction of canopy attributes using assimilated MODIS NDVI values improved EDPM NDVI estimates, bringing the errors down to the average level of 0.1. The ET estimates produced by the coupled scheme were compared with the MODIS evapotranspiration product and with ET from NASA's Mosaic land surface model. The expected r2 = 0.7 ± 0.15 and RMSE = 11.2 ± 4 mm per 8 days achieved in earlier point-based validations were met in this study by the coupling scheme functioning in both prognostic and retrospective modes. Coupled model performance was diminished at the

  16. NON-SPATIAL CALIBRATIONS OF A GENERAL UNIT MODEL FOR ECOSYSTEM SIMULATIONS. (R827169)

    EPA Science Inventory

    General Unit Models simulate system interactions aggregated within one spatial unit of resolution. For unit models to be applicable to spatial computer simulations, they must be formulated generally enough to simulate all habitat elements within the landscape. We present the d...

  17. NON-SPATIAL CALIBRATIONS OF A GENERAL UNIT MODEL FOR ECOSYSTEM SIMULATIONS. (R825792)

    EPA Science Inventory

    General Unit Models simulate system interactions aggregated within one spatial unit of resolution. For unit models to be applicable to spatial computer simulations, they must be formulated generally enough to simulate all habitat elements within the landscape. We present the d...

  18. Generalized simulation technique for turbojet engine system analysis

    NASA Technical Reports Server (NTRS)

    Seldner, K.; Mihaloew, J. R.; Blaha, R. J.

    1972-01-01

    A nonlinear analog simulation of a turbojet engine was developed. The purpose of the study was to establish simulation techniques applicable to propulsion system dynamics and controls research. A schematic model was derived from a physical description of a J85-13 turbojet engine. Basic conservation equations were applied to each component along with their individual performance characteristics to derive a mathematical representation. The simulation was mechanized on an analog computer. The simulation was verified in both steady-state and dynamic modes by comparing analytical results with experimental data obtained from tests performed at the Lewis Research Center with a J85-13 engine. In addition, comparison was also made with performance data obtained from the engine manufacturer. The comparisons established the validity of the simulation technique.

  19. Simulated patients in general practice: a different look at the consultation.

    PubMed Central

    Rethans, J J; van Boven, C P

    1987-01-01

    To develop a better empirical basis for developing quality assessment in general practice three simulated patients made appointments with 48 general practitioners during actual surgery hours and collected facts about their performance. The simulated patients were indistinguishable from real patients and presented a standardised story of a symptomatic urinary tract infection. Two months later the same general practitioners received a written simulation about a patient who had the same urinary tract infection and were asked how they would handle this in real practice. Both results were scored against an existing consensus standard. The overall score for both methods did not show any substantial differences. A more differentiated analysis, however, showed that general practitioners performed significantly better with simulated patients. It also showed that general practitioners answering the written simulation performed significantly more unnecessary and superfluous actions. The results of this study show that the use of simulated patients seems to show the efficient performance of general practitioners in practice. PMID:3105753

  20. Seasonal changes in the atmospheric heat balance simulated by the GISS general circulation model

    NASA Technical Reports Server (NTRS)

    Stone, P. H.; Chow, S.; Helfand, H. M.; Quirk, W. J.; Somerville, R. C. J.

    1975-01-01

    Tests of the ability of numerical general circulation models to simulate the atmosphere have focussed so far on simulations of the January climatology. These models generally present boundary conditions such as sea surface temperature, but this does not prevent testing their ability to simulate seasonal changes in atmospheric processes that accompany presented seasonal changes in boundary conditions. Experiments to simulate changes in the zonally averaged heat balance are discussed since many simplified models of climatic processes are based solely on this balance.

  1. Predicting analysis time in events-driven clinical trials using accumulating time-to-event surrogate information.

    PubMed

    Wang, Jianming; Ke, Chunlei; Yu, Zhinuan; Fu, Lei; Dornseif, Bruce

    2016-05-01

    For clinical trials with time-to-event endpoints, predicting the accrual of the events of interest with precision is critical in determining the timing of interim and final analyses. For example, overall survival (OS) is often chosen as the primary efficacy endpoint in oncology studies, with planned interim and final analyses at a pre-specified number of deaths. Often, correlated surrogate information, such as time-to-progression (TTP) and progression-free survival, are also collected as secondary efficacy endpoints. It would be appealing to borrow strength from the surrogate information to improve the precision of the analysis time prediction. Currently available methods in the literature for predicting analysis timings do not consider utilizing the surrogate information. In this article, using OS and TTP as an example, a general parametric model for OS and TTP is proposed, with the assumption that disease progression could change the course of the overall survival. Progression-free survival, related both to OS and TTP, will be handled separately, as it can be derived from OS and TTP. The authors seek to develop a prediction procedure using a Bayesian method and provide detailed implementation strategies under certain assumptions. Simulations are performed to evaluate the performance of the proposed method. An application to a real study is also provided. Copyright © 2015 John Wiley & Sons, Ltd. PMID:26689725

  2. General approach to boat simulation in virtual reality systems

    NASA Astrophysics Data System (ADS)

    Aranov, Vladislav Y.; Belyaev, Sergey Y.

    2002-02-01

    The paper is dedicated to real time simulation of sport boats, particularly a kayak and high-speed skimming boat, for training goals. This training is issue of the day, since kayaking and riding a high-speed skimming boat are both extreme sports. Participating in such types of competitions puts sportsmen into danger, particularly due to rapids, waterfalls, different water streams, and other obstacles. In order to make the simulation realistic, it is necessary to calculate data for at least 30 frames per second. These calculations may take not more than 5% CPU time, because very time-consuming 3D rendering process takes the rest - 95% CPU time. This paper describes an approach for creating minimal boat simulator models that satisfy the mentioned requirements. Besides, this approach can be used for other watercraft models of this kind.

  3. Verifying Algorithms for Autonomous Aircraft by Simulation Generalities and Example

    NASA Technical Reports Server (NTRS)

    White, Allan L.

    2010-01-01

    An open question in Air Traffic Management is what procedures can be validated by simulation where the simulation shows that the probability of undesirable events is below the required level at some confidence level. The problem is including enough realism to be convincing while retaining enough efficiency to run the large number of trials needed for high confidence. The paper first examines the probabilistic interpretation of a typical requirement by a regulatory agency and computes the number of trials needed to establish the requirement at an equivalent confidence level. Since any simulation is likely to consider only one type of event and there are several types of events, the paper examines under what conditions this separate consideration is valid. The paper establishes a separation algorithm at the required confidence level where the aircraft operates under feedback control as is subject to perturbations. There is a discussion where it is shown that a scenario three of four orders of magnitude more complex is feasible. The question of what can be validated by simulation remains open, but there is reason to be optimistic.

  4. Projectile General Motion in a Vacuum and a Spreadsheet Simulation

    ERIC Educational Resources Information Center

    Benacka, Jan

    2015-01-01

    This paper gives the solution and analysis of projectile motion in a vacuum if the launch and impact heights are not equal. Formulas for the maximum horizontal range and the corresponding angle are derived. An Excel application that simulates the motion is also presented, and the result of an experiment in which 38 secondary school students…

  5. New techniques for meshless flow simulation generalizing moving least squares

    NASA Astrophysics Data System (ADS)

    Trask, Nathaniel; Maxey, Martin

    2015-11-01

    While the Lagrangian nature of SPH offers unique flexibility in application problems, practitioners are forced to choose between compatibility in div/grad operators or low accuracy limiting the scope of the method. In this work, two new discretization frameworks are introduced that extend concepts from finite difference methods to a meshless context: one generalizing the high-order convergence of compact finite differences and another generalizing the enhanced stability of staggered marker-and-cell schemes. The discretizations are based on a novel polynomial reconstruction process that allows arbitrary order polynomial accuracy for both the differential operators and general boundary conditions while maintaining stability and computational efficiency. We demonstrate how the method fits neatly into the ISPH framework and offers a new degree of fidelity and accuracy in Lagrangian particle methods. Supported by the Collaboratory on Mathematics for Mesoscopic Modeling of Materials (CM4), DOE Award DE-SC0009247.

  6. Large eddy simulation using the general circulation model ICON

    NASA Astrophysics Data System (ADS)

    Dipankar, Anurag; Stevens, Bjorn; Heinze, Rieke; Moseley, Christopher; Zängl, Günther; Giorgetta, Marco; Brdar, Slavko

    2015-09-01

    ICON (ICOsahedral Nonhydrostatic) is a unified modeling system for global numerical weather prediction (NWP) and climate studies. Validation of its dynamical core against a test suite for numerical weather forecasting has been recently published by Zängl et al. (2014). In the present work, an extension of ICON is presented that enables it to perform as a large eddy simulation (LES) model. The details of the implementation of the LES turbulence scheme in ICON are explained and test cases are performed to validate it against two standard LES models. Despite the limitations that ICON inherits from being a unified modeling system, it performs well in capturing the mean flow characteristics and the turbulent statistics of two simulated flow configurations—one being a dry convective boundary layer and the other a cumulus-topped planetary boundary layer.

  7. Projectile general motion in a vacuum and a spreadsheet simulation

    NASA Astrophysics Data System (ADS)

    Benacka, Jan

    2015-01-01

    This paper gives the solution and analysis of projectile motion in a vacuum if the launch and impact heights are not equal. Formulas for the maximum horizontal range and the corresponding angle are derived. An Excel application that simulates the motion is also presented, and the result of an experiment in which 38 secondary school students developed the application and investigated the system is given. A questionnaire survey was carried out to find out whether the students found the lessons interesting, learned new skills and wanted to model projectile motion in the air as an example of more realistic motion. The results are discussed.

  8. A General Simulation Method for Multiple Bodies in Proximate Flight

    NASA Technical Reports Server (NTRS)

    Meakin, Robert L.

    2003-01-01

    Methods of unsteady aerodynamic simulation for an arbitrary number of independent bodies flying in close proximity are considered. A novel method to efficiently detect collision contact points is described. A method to compute body trajectories in response to aerodynamic loads, applied loads, and inter-body collisions is also given. The physical correctness of the methods are verified by comparison to a set of analytic solutions. The methods, combined with a Navier-Stokes solver, are used to demonstrate the possibility of predicting the unsteady aerodynamics and flight trajectories of moving bodies that involve rigid-body collisions.

  9. BIRD: A general interface for sparse distributed memory simulators

    NASA Technical Reports Server (NTRS)

    Rogers, David

    1990-01-01

    Kanerva's sparse distributed memory (SDM) has now been implemented for at least six different computers, including SUN3 workstations, the Apple Macintosh, and the Connection Machine. A common interface for input of commands would both aid testing of programs on a broad range of computer architectures and assist users in transferring results from research environments to applications. A common interface also allows secondary programs to generate command sequences for a sparse distributed memory, which may then be executed on the appropriate hardware. The BIRD program is an attempt to create such an interface. Simplifying access to different simulators should assist developers in finding appropriate uses for SDM.

  10. Plasma Jet Simulations Using a Generalized Ohm's Law

    NASA Technical Reports Server (NTRS)

    Ebersohn, Frans; Shebalin, John V.; Girimaji, Sharath S.

    2012-01-01

    Plasma jets are important physical phenomena in astrophysics and plasma propulsion devices. A currently proposed dual jet plasma propulsion device to be used for ISS experiments strongly resembles a coronal loop and further draws a parallel between these physical systems [1]. To study plasma jets we use numerical methods that solve the compressible MHD equations using the generalized Ohm s law [2]. Here, we will discuss the crucial underlying physics of these systems along with the numerical procedures we utilize to study them. Recent results from our numerical experiments will be presented and discussed.

  11. Optimal generalized multistep integration formulae for real-time digital simulation

    NASA Technical Reports Server (NTRS)

    Moerder, D. D.; Halyo, N.

    1985-01-01

    The problem of discretizing a dynamical system for real-time digital simulation is considered. Treating the system and its simulation as stochastic processes leads to a statistical characterization of simulator fidelity. A plant discretization procedure based on an efficient matrix generalization of explicit linear multistep discrete integration formulae is introduced, which minimizes a weighted sum of the mean squared steady-state and transient error between the system and simulator outputs.

  12. GOOSE, a generalized object-oriented simulation environment

    SciTech Connect

    Ford, C.E.; March-Leuba, C. ); Guimaraes, L.; Ugolini, D. . Dept. of Nuclear Engineering)

    1991-01-01

    GOOSE, prototype software for a fully interactive, object-oriented simulation environment, is being developed as part of the Advanced Controls Program at Oak Ridge National Laboratory. Dynamic models may easily be constructed and tested; fully interactive capabilities allow the user to alter model parameters and complexity without recompilation. This environment provides access to powerful tools, such as numerical integration packages, graphical displays, and online help. Portability has been an important design goal; the system was written in Objective-C in order to run on a wide variety of computers and operating systems, including UNIX workstations and personal computers. A detailed library of nuclear reactor components, currently under development, will also be described. 5 refs., 4 figs.

  13. Parametrizing linear generalized Langevin dynamics from explicit molecular dynamics simulations

    SciTech Connect

    Gottwald, Fabian; Karsten, Sven; Ivanov, Sergei D. Kühn, Oliver

    2015-06-28

    Fundamental understanding of complex dynamics in many-particle systems on the atomistic level is of utmost importance. Often the systems of interest are of macroscopic size but can be partitioned into a few important degrees of freedom which are treated most accurately and others which constitute a thermal bath. Particular attention in this respect attracts the linear generalized Langevin equation, which can be rigorously derived by means of a linear projection technique. Within this framework, a complicated interaction with the bath can be reduced to a single memory kernel. This memory kernel in turn is parametrized for a particular system studied, usually by means of time-domain methods based on explicit molecular dynamics data. Here, we discuss that this task is more naturally achieved in frequency domain and develop a Fourier-based parametrization method that outperforms its time-domain analogues. Very surprisingly, the widely used rigid bond method turns out to be inappropriate in general. Importantly, we show that the rigid bond approach leads to a systematic overestimation of relaxation times, unless the system under study consists of a harmonic bath bi-linearly coupled to the relevant degrees of freedom.

  14. General relativistic magnetohydrodynamical simulations of the jet in M 87

    NASA Astrophysics Data System (ADS)

    Mościbrodzka, Monika; Falcke, Heino; Shiokawa, Hotaka

    2016-02-01

    Context. The connection between black hole, accretion disk, and radio jet can be constrained best by fitting models to observations of nearby low-luminosity galactic nuclei, in particular the well-studied sources Sgr A* and M 87. There has been considerable progress in modeling the central engine of active galactic nuclei by an accreting supermassive black hole coupled to a relativistic plasma jet. However, can a single model be applied to a range of black hole masses and accretion rates? Aims: Here we want to compare the latest three-dimensional numerical model, originally developed for Sgr A* in the center of the Milky Way, to radio observations of the much more powerful and more massive black hole in M 87. Methods: We postprocess three-dimensional GRMHD models of a jet-producing radiatively inefficient accretion flow around a spinning black hole using relativistic radiative transfer and ray-tracing to produce model spectra and images. As a key new ingredient in these models, we allow the proton-electron coupling in these simulations depend on the magnetic properties of the plasma. Results: We find that the radio emission in M 87 is described well by a combination of a two-temperature accretion flow and a hot single-temperature jet. Most of the radio emission in our simulations comes from the jet sheath. The model fits the basic observed characteristics of the M 87 radio core: it is "edge-brightened", starts subluminally, has a flat spectrum, and increases in size with wavelength. The best fit model has a mass-accretion rate of Ṁ ~ 9 × 10-3M⊙ yr-1 and a total jet power of Pj ~ 1043 erg s-1. Emission at λ = 1.3 mm is produced by the counter-jet close to the event horizon. Its characteristic crescent shape surrounding the black hole shadow could be resolved by future millimeter-wave VLBI experiments. Conclusions: The model was successfully derived from one for the supermassive black hole in the center of the Milky Way by appropriately scaling mass and

  15. A general Kirchhoff approximation for echo simulation in ultrasonic NDT

    NASA Astrophysics Data System (ADS)

    Dorval, V.; Chatillon, S.; Lu, B.; Darmon, M.; Mahaut, S.

    2012-05-01

    The Kirchhoff approximation is commonly used for the modeling of echoes in ultrasonic NDE. It consists in locally approximating the illuminated surface by an infinite plane to compute elastic fields. A model based on this approximation is used in the CIVA software, developed at CEA LIST, to compute echoes from cracks and backwalls. In its current version, it is limited to stress-free surfaces. A new model using a more general formalism has been developed. It is based on reciprocity principles and is valid for any host and flaw materials (liquids, isotropic and anisotropic solids). Experimental validations confirm that this new model can be used for a wider range of applications than the previous one. A second part of this communication deals with the improvement of the Kirchhoff approximation in the aim of predicting diffraction echoes. It is based on an approach called refined Kirchhoff, which combines the Kirchhoff and Geometrical Theory of Diffraction (GTD) models. An illustration of this method for the case of a rigid obstacle in a fluid is given.

  16. Mars' Thermal Structure From The Lower To Middle Atmosphere: NASA Ames Mars General Circulation Simulations

    NASA Astrophysics Data System (ADS)

    Brecht, A. S.; Hollingsworth, J. L.; Kahre, M. A.

    2014-07-01

    The NASA Ames Mars General Ciculation Model (MGCM) has been extended to incorporate the middle atmosphere (~80 km to ~120 km). The extended MGCM simulated thermal structure will be compared to MRO-MCS and MEx-SPICAM observations.

  17. No Vent Tank Fill and Transfer Line Chilldown Analysis by Generalized Fluid System Simulation Program (GFSSP)

    NASA Technical Reports Server (NTRS)

    Majumdar, Alok

    2013-01-01

    The purpose of the paper is to present the analytical capability developed to model no vent chill and fill of cryogenic tank to support CPST (Cryogenic Propellant Storage and Transfer) program. Generalized Fluid System Simulation Program (GFSSP) was adapted to simulate charge-holdvent method of Tank Chilldown. GFSSP models were developed to simulate chilldown of LH2 tank in K-site Test Facility and numerical predictions were compared with test data. The report also describes the modeling technique of simulating the chilldown of a cryogenic transfer line and GFSSP models were developed to simulate the chilldown of a long transfer line and compared with test data.

  18. Instructor and student pilots' subjective evaluation of a general aviation simulator with a terrain visual system

    NASA Technical Reports Server (NTRS)

    Kiteley, G. W.; Harris, R. L., Sr.

    1978-01-01

    Ten student pilots were given a 1 hour training session in the NASA Langley Research Center's General Aviation Simulator by a certified flight instructor and a follow-up flight evaluation was performed by the student's own flight instructor, who has also flown the simulator. The students and instructors generally felt that the simulator session had a positive effect on the students. They recommended that a simulator with a visual scene and a motion base would be useful in performing such maneuvers as: landing approaches, level flight, climbs, dives, turns, instrument work, and radio navigation, recommending that the simulator would be an efficient means of introducing the student to new maneuvers before doing them in flight. The students and instructors estimated that about 8 hours of simulator time could be profitably devoted to the private pilot training.

  19. A Comparative Analysis of General Case Simulation Instruction and Naturalistic Instruction.

    ERIC Educational Resources Information Center

    Domaracki, Joseph W.; Lyon, Steven R.

    1992-01-01

    This study, which involved training four young adults with moderate or severe mental retardation on housekeeping and janitorial work skills, found that simulation instruction based on general case methodology can be used to teach complex sequences, that naturalistic instruction seemed more efficient than simulation instruction, and that neither…

  20. General specifications for the development of a PC-based simulator of the NASA RECON system

    NASA Technical Reports Server (NTRS)

    Dominick, Wayne D. (Editor); Triantafyllopoulos, Spiros

    1984-01-01

    The general specifications for the design and implementation of an IBM PC/XT-based simulator of the NASA RECON system, including record designs, file structure designs, command language analysis, program design issues, error recovery considerations, and usage monitoring facilities are discussed. Once implemented, such a simulator will be utilized to evaluate the effectiveness of simulated information system access in addition to actual system usage as part of the total educational programs being developed within the NASA contract.

  1. A General Simulator Using State Estimation for a Space Tug Navigation System. [computerized simulation, orbital position estimation and flight mechanics

    NASA Technical Reports Server (NTRS)

    Boland, J. S., III

    1975-01-01

    A general simulation program is presented (GSP) involving nonlinear state estimation for space vehicle flight navigation systems. A complete explanation of the iterative guidance mode guidance law, derivation of the dynamics, coordinate frames, and state estimation routines are given so as to fully clarify the assumptions and approximations involved so that simulation results can be placed in their proper perspective. A complete set of computer acronyms and their definitions as well as explanations of the subroutines used in the GSP simulator are included. To facilitate input/output, a complete set of compatable numbers, with units, are included to aid in data development. Format specifications, output data phrase meanings and purposes, and computer card data input are clearly spelled out. A large number of simulation and analytical studies were used to determine the validity of the simulator itself as well as various data runs.

  2. General purpose simulation system of the data management system for space shuttle mission 18

    NASA Technical Reports Server (NTRS)

    Bengtson, N. M.; Mellichamp, J. M.; Crenshaw, J.

    1975-01-01

    The simulation program of the science and engineering data management system for the space shuttle is presented. The programming language used was General Purpose Simulation System V (OS). The data flow was modeled from its origin at the experiments or subsystems to transmission from the space shuttle. Mission 18 was the particular flight chosen for simulation. First, the general structure of the program is presented and the trade studies which were performed are identified. Inputs required to make runs are discussed followed by identification of the output statistics. Some areas for model modifications are pointed out. A detailed model configuration, program listing and results are included.

  3. Cloud-radiative effects on implied oceanic energy transport as simulated by atmospheric general circulation models

    NASA Technical Reports Server (NTRS)

    Gleckler, P. J.; Randall, D. A.; Boer, G.; Colman, R.; Dix, M.; Galin, V.; Helfand, M.; Kiehl, J.; Kitoh, A.; Lau, W.

    1995-01-01

    This paper summarizes the ocean surface net energy flux simulated by fifteen atmospheric general circulation models constrained by realistically-varying sea surface temperatures and sea ice as part of the Atmospheric Model Intercomparison Project. In general, the simulated energy fluxes are within the very large observational uncertainties. However, the annual mean oceanic meridional heat transport that would be required to balance the simulated surface fluxes is shown to be critically sensitive to the radiative effects of clouds, to the extent that even the sign of the Southern Hemisphere ocean heat transport can be affected by the errors in simulated cloud-radiation interactions. It is suggested that improved treatment of cloud radiative effects should help in the development of coupled atmosphere-ocean general circulation models.

  4. Generalized image charge solvation model for electrostatic interactions in molecular dynamics simulations of aqueous solutions

    NASA Astrophysics Data System (ADS)

    Deng, Shaozhong; Xue, Changfeng; Baumketner, Andriy; Jacobs, Donald; Cai, Wei

    2013-07-01

    This paper extends the image charge solvation model (ICSM) [Y. Lin, A. Baumketner, S. Deng, Z. Xu, D. Jacobs, W. Cai, An image-based reaction field method for electrostatic interactions in molecular dynamics simulations of aqueous solutions, J. Chem. Phys. 131 (2009) 154103], a hybrid explicit/implicit method to treat electrostatic interactions in computer simulations of biomolecules formulated for spherical cavities, to prolate spheroidal and triaxial ellipsoidal cavities, designed to better accommodate non-spherical solutes in molecular dynamics (MD) simulations. In addition to the utilization of a general truncated octahedron as the MD simulation box, central to the proposed extension is an image approximation method to compute the reaction field for a point charge placed inside such a non-spherical cavity by using a single image charge located outside the cavity. The resulting generalized image charge solvation model (GICSM) is tested in simulations of liquid water, and the results are analyzed in comparison with those obtained from the ICSM simulations as a reference. We find that, for improved computational efficiency due to smaller simulation cells and consequently a less number of explicit solvent molecules, the generalized model can still faithfully reproduce known static and dynamic properties of liquid water at least for systems considered in the present paper, indicating its great potential to become an accurate but more efficient alternative to the ICSM when bio-macromolecules of irregular shapes are to be simulated.

  5. A Multi-mission Event-Driven Component-Based System for Support of Flight Software Development, ATLO, and Operations first used by the Mars Science Laboratory (MSL) Project

    NASA Technical Reports Server (NTRS)

    Dehghani, Navid; Tankenson, Michael

    2006-01-01

    This paper details an architectural description of the Mission Data Processing and Control System (MPCS), an event-driven, multi-mission ground data processing components providing uplink, downlink, and data management capabilities which will support the Mars Science Laboratory (MSL) project as its first target mission. MPCS is developed based on a set of small reusable components, implemented in Java, each designed with a specific function and well-defined interfaces. An industry standard messaging bus is used to transfer information among system components. Components generate standard messages which are used to capture system information, as well as triggers to support the event-driven architecture of the system. Event-driven systems are highly desirable for processing high-rate telemetry (science and engineering) data, and for supporting automation for many mission operations processes.

  6. Tryout of a General Purpose Simulator in an Air National Guard Training Environment. Interim Report, June 1974-August 1974.

    ERIC Educational Resources Information Center

    Spangenberg, Ronald W.

    An evaluation of the usability, effectiveness, and acceptance in a job environment was performed on a general purpose simulator using a simulation of a radar system. General purpose simulators permit sharing of a programable capacity among simulations, thus providing economical hands-on training and training not usually economically available by…

  7. A GeneralizedWeight-Based Particle-In-Cell Simulation Scheme

    SciTech Connect

    W.W. Lee, T.G. Jenkins and S. Ethier

    2010-02-02

    A generalized weight-based particle simulation scheme suitable for simulating magnetized plasmas, where the zeroth-order inhomogeneity is important, is presented. The scheme is an extension of the perturbative simulation schemes developed earlier for particle-in-cell (PIC) simulations. The new scheme is designed to simulate both the perturbed distribution (δf) and the full distribution (full-F) within the same code. The development is based on the concept of multiscale expansion, which separates the scale lengths of the background inhomogeneity from those associated with the perturbed distributions. The potential advantage for such an arrangement is to minimize the particle noise by using δf in the linear stage stage of the simulation, while retaining the flexibility of a full-F capability in the fully nonlinear stage of the development when signals associated with plasma turbulence are at a much higher level than those from the intrinsic particle noise.

  8. Simulation of the great plains low-level jet and associated clouds by general circulation models

    SciTech Connect

    Ghan, S.J.; Bian, X.; Corsetti, L.

    1996-07-01

    The low-level jet frequently observed in the Great Plains of the United States forms preferentially at night and apparently influences the timing of the thunderstorms in the region. The authors have found that both the European Centre for Medium-Range Weather Forecasts general circulation model and the National Center for Atmospheric Research Community Climate Model simulate the low-level jet rather well, although the spatial distribution of the jet frequency simulated by the two GCM`s differ considerably. Sensitivity experiments have demonstrated that the simulated low-level jet is surprisingly robust, with similar simulations at much coarser horizontal and vertical resolutions. However, both GCM`s fail to simulate the observed relationship between clouds and the low-level jet. The pronounced nocturnal maximum in thunderstorm frequency associated with the low-level jet is not simulated well by either GCM, with only weak evidence of a nocturnal maximum in the Great Plains. 36 refs., 20 figs.

  9. Evaluating the GPSS simulation model for the Viking batch computer system. [General Purpose Simulation System

    NASA Technical Reports Server (NTRS)

    Lee, J.-J.

    1976-01-01

    In anticipation of extremely heavy loading requirements by the Viking mission during the post-landing periods, a GPSS model has been developed for the purpose of simulating these requirements on the Viking batch computer system. This paper presents the effort pursued in evaluating such a model and results thereby obtained. The evaluation effort consists of selecting the evaluation approach, collecting actual test run data, making comparisons and deriving conclusions.

  10. General Relativistic Radiative Transfer and GeneralRelativistic MHD Simulations of Accretion and Outflows of Black Holes

    SciTech Connect

    Fuerst, Steven V.; Mizuno, Yosuke; Nishikawa, Ken-Ichi; Wu, Kinwah; /Mullard Space Sci. Lab.

    2007-01-05

    We calculate the emission from relativistic flows in black hole systems using a fully general relativistic radiative transfer formulation, with flow structures obtained by general relativistic magneto-hydrodynamic simulations. We consider thermal free-free emission and thermal synchrotron emission. Bright filament-like features protrude (visually) from the accretion disk surface, which are enhancements of synchrotron emission where the magnetic field roughly aligns with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We propose that variations and drifts of the features produce certain X-ray quasi-periodic oscillations (QPOs) observed in black-hole X-ray binaries.

  11. Simulating the Illuminance and Efficiency of the LEDs Used in General Household Lighting

    NASA Astrophysics Data System (ADS)

    Sun, Wen-Shing; Tsuei, Chih-Hsuan; Huang, Yi-Han

    The advantage of the LEDs illumination in general household lighting was proposed. High efficiency white LED as the light source was provided the energy saving illumination of the general household lightings. Different spaces in general household with different standards of average illuminance were designed and simulated by LightTools and DIALux software. The power consumptions and efficiency of traditional illuminated light sources and LED light source in lighting the household environment were analyzed and compared with each other at the same standard of average illuminance. Finally, it provided the advantage of using white LEDs in different spaces of the general household lighting.

  12. General purpose simulation system of the data management system for Space Shuttle mission 18

    NASA Technical Reports Server (NTRS)

    Bengtson, N. M.; Mellichamp, J. M.; Smith, O. C.

    1976-01-01

    A simulation program for the flow of data through the Data Management System of Spacelab and Space Shuttle was presented. The science, engineering, command and guidance, navigation and control data were included. The programming language used was General Purpose Simulation System V (OS). The science and engineering data flow was modeled from its origin at the experiments and subsystems to transmission from Space Shuttle. Command data flow was modeled from the point of reception onboard and from the CDMS Control Panel to the experiments and subsystems. The GN&C data flow model handled data between the General Purpose Computer and the experiments and subsystems. Mission 18 was the particular flight chosen for simulation. The general structure of the program is presented, followed by a user's manual. Input data required to make runs are discussed followed by identification of the output statistics. The appendices contain a detailed model configuration, program listing and results.

  13. Using a million cell simulation of the cerebellum: network scaling and task generality

    PubMed Central

    Li, Wen-Ke; Hausknecht, Matthew J.; Stone, Peter H.; Mauk, Michael D.

    2012-01-01

    Several factors combine to make it feasible to build computer simulations of the cerebellum and to test them in biologically realistic ways. These simulations can be used to help understand the computational contributions of various cerebellar components, including the relevance of the enormous number of neurons in the granule cell layer. In previous work we have used a simulation containing 12000 granule cells to develop new predictions and to account for various aspects of eyelid conditioning, a form of motor learning mediated by the cerebellum. Here we demonstrate the feasibility of scaling up this simulation to over one million granule cells using parallel graphics processing unit (GPU) technology. We observe that this increase in number of granule cells requires only twice the execution time of the smaller simulation on the GPU. We demonstrate that this simulation, like its smaller predecessor, can emulate certain basic features of conditioned eyelid responses, with a slight improvement in performance in one measure. We also use this simulation to examine the generality of the computation properties that we have derived from studying eyelid conditioning. We demonstrate that this scaled up simulation can learn a high level of performance in a classic machine learning task, the cart-pole balancing task. These results suggest that this parallel GPU technology can be used to build very large-scale simulations whose connectivity ratios match those of the real cerebellum and that these simulations can be used guide future studies on cerebellar mediated tasks and on machine learning problems. PMID:23200194

  14. Experiments in monthly mean simulation of the atmosphere with a coarse-mesh general circulation model

    NASA Technical Reports Server (NTRS)

    Lutz, R. J.; Spar, J.

    1978-01-01

    The Hansen atmospheric model was used to compute five monthly forecasts (October 1976 through February 1977). The comparison is based on an energetics analysis, meridional and vertical profiles, error statistics, and prognostic and observed mean maps. The monthly mean model simulations suffer from several defects. There is, in general, no skill in the simulation of the monthly mean sea-level pressure field, and only marginal skill is indicated for the 850 mb temperatures and 500 mb heights. The coarse-mesh model appears to generate a less satisfactory monthly mean simulation than the finer mesh GISS model.

  15. General circulation model simulations of winter and summer sea-level pressures over North America

    USGS Publications Warehouse

    McCabe, G.J., Jr.; Legates, D.R.

    1992-01-01

    In this paper, observed sea-level pressures were used to evaluate winter and summer sea-level pressures over North America simulated by the Goddard Institute for Space Studies (GISS) and the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation models. The objective of the study is to determine how similar the spatial and temporal distributions of GCM-simulated daily sea-level pressures over North America are to observed distributions. Overall, both models are better at reproducing observed within-season variance of winter and summer sea-level pressures than they are at simulating the magnitude of mean winter and summer sea-level pressures. -from Authors

  16. Development and evaluation of a general aviation real world noise simulator

    NASA Technical Reports Server (NTRS)

    Galanter, E.; Popper, R.

    1980-01-01

    An acoustic playback system is described which realistically simulates the sounds experienced by the pilot of a general aviation aircraft during engine idle, take-off, climb, cruise, descent, and landing. The physical parameters of the signal as they appear in the simulator environment are compared to analogous parameters derived from signals recorded during actual flight operations. The acoustic parameters of the simulated and real signals during cruise conditions are within plus or minus two dB in third octave bands from 0.04 to 4 kHz. The overall A-weighted levels of the signals are within one dB of signals generated in the actual aircraft during equivalent maneuvers. Psychoacoustic evaluations of the simulator signal are compared with similar measurements based on transcriptions of actual aircraft signals. The subjective judgments made by human observers support the conclusion that the simulated sound closely approximates transcribed sounds of real aircraft.

  17. Development and evaluation of a general aviation real world noise simulator

    NASA Astrophysics Data System (ADS)

    Galanter, E.; Popper, R.

    1980-03-01

    An acoustic playback system is described which realistically simulates the sounds experienced by the pilot of a general aviation aircraft during engine idle, take-off, climb, cruise, descent, and landing. The physical parameters of the signal as they appear in the simulator environment are compared to analogous parameters derived from signals recorded during actual flight operations. The acoustic parameters of the simulated and real signals during cruise conditions are within plus or minus two dB in third octave bands from 0.04 to 4 kHz. The overall A-weighted levels of the signals are within one dB of signals generated in the actual aircraft during equivalent maneuvers. Psychoacoustic evaluations of the simulator signal are compared with similar measurements based on transcriptions of actual aircraft signals. The subjective judgments made by human observers support the conclusion that the simulated sound closely approximates transcribed sounds of real aircraft.

  18. The General-Use Nodal Network Solver (GUNNS) Modeling Package for Space Vehicle Flow System Simulation

    NASA Technical Reports Server (NTRS)

    Harvey, Jason; Moore, Michael

    2013-01-01

    The General-Use Nodal Network Solver (GUNNS) is a modeling software package that combines nodal analysis and the hydraulic-electric analogy to simulate fluid, electrical, and thermal flow systems. GUNNS is developed by L-3 Communications under the TS21 (Training Systems for the 21st Century) project for NASA Johnson Space Center (JSC), primarily for use in space vehicle training simulators at JSC. It has sufficient compactness and fidelity to model the fluid, electrical, and thermal aspects of space vehicles in real-time simulations running on commodity workstations, for vehicle crew and flight controller training. It has a reusable and flexible component and system design, and a Graphical User Interface (GUI), providing capability for rapid GUI-based simulator development, ease of maintenance, and associated cost savings. GUNNS is optimized for NASA's Trick simulation environment, but can be run independently of Trick.

  19. Estimating plant available water for general crop simulations in ALMANAC/APEX/EPIC/SWAT

    Technology Transfer Automated Retrieval System (TEKTRAN)

    Process-based simulation models ALMANAC/APEX/EPIC/SWAT contain generalized plant growth subroutines to predict biomass and crop yield. Environmental constraints typically restrict plant growth and yield. Water stress is often an important limiting factor; it is calculated as the sum of water use f...

  20. Computer considerations for real time simulation of a generalized rotor model

    NASA Technical Reports Server (NTRS)

    Howe, R. M.; Fogarty, L. E.

    1977-01-01

    Scaled equations were developed to meet requirements for real time computer simulation of the rotor system research aircraft. These equations form the basis for consideration of both digital and hybrid mechanization for real time simulation. For all digital simulation estimates of the required speed in terms of equivalent operations per second are developed based on the complexity of the equations and the required intergration frame rates. For both conventional hybrid simulation and hybrid simulation using time-shared analog elements the amount of required equipment is estimated along with a consideration of the dynamic errors. Conventional hybrid mechanization using analog simulation of those rotor equations which involve rotor-spin frequencies (this consititutes the bulk of the equations) requires too much analog equipment. Hybrid simulation using time-sharing techniques for the analog elements appears possible with a reasonable amount of analog equipment. All-digital simulation with affordable general-purpose computers is not possible because of speed limitations, but specially configured digital computers do have the required speed and consitute the recommended approach.

  1. Hybrid General Pattern Search and Simulated Annealing for Industrail Production Planning Problems

    NASA Astrophysics Data System (ADS)

    Vasant, P.; Barsoum, N.

    2010-06-01

    In this paper, the hybridization of GPS (General Pattern Search) method and SA (Simulated Annealing) incorporated in the optimization process in order to look for the global optimal solution for the fitness function and decision variables as well as minimum computational CPU time. The real strength of SA approach been tested in this case study problem of industrial production planning. This is due to the great advantage of SA for being easily escaping from trapped in local minima by accepting up-hill move through a probabilistic procedure in the final stages of optimization process. Vasant [1] in his Ph. D thesis has provided 16 different techniques of heuristic and meta-heuristic in solving industrial production problems with non-linear cubic objective functions, eight decision variables and 29 constraints. In this paper, fuzzy technological problems have been solved using hybrid techniques of general pattern search and simulated annealing. The simulated and computational results are compared to other various evolutionary techniques.

  2. Serial Generalized Ensemble Simulations of Biomolecules with Self-Consistent Determination of Weights.

    PubMed

    Chelli, Riccardo; Signorini, Giorgio F

    2012-03-13

    Serial generalized ensemble simulations, such as simulated tempering, enhance phase space sampling through non-Boltzmann weighting protocols. The most critical aspect of these methods with respect to the popular replica exchange schemes is the difficulty in determining the weight factors which enter the criterion for accepting replica transitions between different ensembles. Recently, a method, called BAR-SGE, was proposed for estimating optimal weight factors by resorting to a self-consistent procedure applied during the simulation (J. Chem. Theory Comput.2010, 6, 1935-1950). Calculations on model systems have shown that BAR-SGE outperforms other approaches proposed for determining optimal weights in serial generalized ensemble simulations. However, extensive tests on real systems and on convergence features with respect to the replica exchange method are lacking. Here, we report on a thorough analysis of BAR-SGE by performing molecular dynamics simulations of a solvated alanine dipeptide, a system often used as a benchmark to test new computational methodologies, and comparing results to the replica exchange method. To this aim, we have supplemented the ORAC program, a FORTRAN suite for molecular dynamics simulations (J. Comput. Chem.2010, 31, 1106-1116), with several variants of the BAR-SGE technique. An illustration of the specific BAR-SGE algorithms implemented in the ORAC program is also provided. PMID:26593345

  3. Gyrokinetic particle simulation of microturbulence for general magnetic geometry and experimental profiles

    SciTech Connect

    Xiao, Yong; Holod, Ihor; Wang, Zhixuan; Lin, Zhihong; Zhang, Taige

    2015-02-15

    Developments in gyrokinetic particle simulation enable the gyrokinetic toroidal code (GTC) to simulate turbulent transport in tokamaks with realistic equilibrium profiles and plasma geometry, which is a critical step in the code–experiment validation process. These new developments include numerical equilibrium representation using B-splines, a new Poisson solver based on finite difference using field-aligned mesh and magnetic flux coordinates, a new zonal flow solver for general geometry, and improvements on the conventional four-point gyroaverage with nonuniform background marker loading. The gyrokinetic Poisson equation is solved in the perpendicular plane instead of the poloidal plane. Exploiting these new features, GTC is able to simulate a typical DIII-D discharge with experimental magnetic geometry and profiles. The simulated turbulent heat diffusivity and its radial profile show good agreement with other gyrokinetic codes. The newly developed nonuniform loading method provides a modified radial transport profile to that of the conventional uniform loading method.

  4. The Tropical Subseasonal Variability Simulated in the NASA GISS General Circulation Model

    NASA Technical Reports Server (NTRS)

    Kim, Daehyun; Sobel, Adam H.; DelGenio, Anthony D.; Chen, Yonghua; Camargo, Suzana J.; Yao, Mao-Sung; Kelley, Maxwell; Nazarenko, Larissa

    2012-01-01

    The tropical subseasonal variability simulated by the Goddard Institute for Space Studies general circulation model, Model E2, is examined. Several versions of Model E2 were developed with changes to the convective parameterization in order to improve the simulation of the Madden-Julian oscillation (MJO). When the convective scheme is modified to have a greater fractional entrainment rate, Model E2 is able to simulate MJO-like disturbances with proper spatial and temporal scales. Increasing the rate of rain reevaporation has additional positive impacts on the simulated MJO. The improvement in MJO simulation comes at the cost of increased biases in the mean state, consistent in structure and amplitude with those found in other GCMs when tuned to have a stronger MJO. By reinitializing a relatively poor-MJO version with restart files from a relatively better-MJO version, a series of 30-day integrations is constructed to examine the impacts of the parameterization changes on the organization of tropical convection. The poor-MJO version with smaller entrainment rate has a tendency to allow convection to be activated over a broader area and to reduce the contrast between dry and wet regimes so that tropical convection becomes less organized. Besides the MJO, the number of tropical-cyclone-like vortices simulated by the model is also affected by changes in the convection scheme. The model simulates a smaller number of such storms globally with a larger entrainment rate, while the number increases significantly with a greater rain reevaporation rate.

  5. Robustness of a high-resolution central scheme for hydrodynamic simulations in full general relativity

    NASA Astrophysics Data System (ADS)

    Shibata, Masaru; Font, José A.

    2005-08-01

    A recent paper by Lucas-Serrano et al. [A. Lucas-Serrano, J. A. Font, J. M. Ibánez, and J. M. Martí, Astron. Astrophys. 428, 703 (2004)] indicates that a high-resolution central (HRC) scheme is robust enough to yield accurate hydrodynamical simulations of special relativistic flows in the presence of ultrarelativistic speeds and strong shock waves. In this paper we apply this scheme in full general relativity (involving dynamical spacetimes), and assess its suitability by performing test simulations for oscillations of rapidly rotating neutron stars and merger of binary neutron stars. It is demonstrated that this HRC scheme can yield results as accurate as those by the so-called high-resolution shock-capturing (HRSC) schemes based upon Riemann solvers. Furthermore, the adopted HRC scheme has increased computational efficiency as it avoids the costly solution of Riemann problems and has practical advantages in the modeling of neutron star spacetimes. Namely, it allows simulations with stiff equations of state by successfully dealing with very low-density unphysical atmospheres. These facts not only suggest that such a HRC scheme may be a desirable tool for hydrodynamical simulations in general relativity, but also open the possibility to perform accurate magnetohydrodynamical simulations in curved dynamic spacetimes.

  6. The Early Jurassic climate: General circulation model simulations and the paleoclimate record

    SciTech Connect

    Chandler, M.A.

    1992-01-01

    This thesis presents the results of several general circulation model simulations of the Early Jurassic climate. The general circulation model employed was developed at the Goddard Institute for Space Studies while most paleoclimate data were provided by the Paleographic Atlas Project of the University of Chicago. The first chapter presents an Early Jurassic base simulation, which uses detailed reconstructions of paleogeography, vegetation, and sea surface temperature as boundary condition data sets. The resulting climatology reveals an Earth 5.2[degrees]C warmer, globally, than at present and a latitudinal temperature gradient dominated by high-latitude warming (+20[degrees]C) and little tropical change (+1[degrees]C). Comparisons show a good correlation between simulated results and paleoclimate data. Sensitivity experiments are used to investigate any model-data mismatches. Chapters two and three discuss two important aspects of Early Jurassic climate, continental aridity and global warming. Chapter two focuses on the hydrological capabilities of the general circulation model. The general circulation model's hydrologic diagnostics are evaluated, using the distribution of modern deserts and Early Jurassic paleoclimate data as validating constraints. A new method, based on general circulation model diagnostics and empirical formulae, is proposed for evaluating moisture balance. Chapter three investigates the cause of past global warming, concentrating on the role of increased ocean heat transport. Early Jurassic simulations show that increased ocean heat transports may have been a major factor in past climates. Increased ocean heat transports create latitudinal temperature gradients that closely approximate paleoclimate data and solve the problem of tropical overheating that results from elevated atmospheric carbon dioxide. Increased carbon dioxide cannot duplicate the Jurassic climate without also including increased ocean heat transports.

  7. Cloud-radiative effects on implied oceanic energy transports as simulated by atmospheric general circulation models

    SciTech Connect

    Gleckler, P.J.; Randall, D.A.; Boer, G.

    1994-03-01

    This paper reports on energy fluxes across the surface of the ocean as simulated by fifteen atmospheric general circulation models in which ocean surface temperatures and sea-ice boundaries are prescribed. The oceanic meridional energy transport that would be required to balance these surface fluxes is computed, and is shown to be critically sensitive to the radiative effects of clouds, to the extent that even the sign of the Southern Hemisphere ocean energy transport can be affected by the errors in simulated cloud-radiation interactions.

  8. General-relativistic simulations of binary black hole-neutron stars: Precursor electromagnetic signals

    NASA Astrophysics Data System (ADS)

    Paschalidis, Vasileios; Etienne, Zachariah B.; Shapiro, Stuart L.

    2013-07-01

    We perform the first general relativistic force-free simulations of neutron star magnetospheres in orbit about spinning and nonspinning black holes. We find promising precursor electromagnetic emission: typical Poynting luminosities at, e.g., an orbital separation of r=6.6RNS are LEM˜6×1042(BNS,p/1013G)2(MNS/1.4M⊙)2erg/s. The Poynting flux peaks within a broad beam of ˜40° in the azimuthal direction and within ˜60° from the orbital plane, establishing a possible lighthouse effect. Our calculations, though preliminary, preview more detailed simulations of these systems that we plan to perform in the future.

  9. Generalized math model for simulation of high-altitude balloon systems

    NASA Technical Reports Server (NTRS)

    Nigro, N. J.; Elkouh, A. F.; Hinton, D. E.; Yang, J. K.

    1985-01-01

    Balloon systems have proved to be a cost-effective means for conducting research experiments (e.g., infrared astronomy) in the earth's atmosphere. The purpose of this paper is to present a generalized mathematical model that can be used to simulate the motion of these systems once they have attained float altitude. The resulting form of the model is such that the pendulation and spin motions of the system are uncoupled and can be analyzed independently. The model is evaluated by comparing the simulation results with data obtained from an actual balloon system flown by NASA.

  10. Evaluation of the Event Driven Phenology Model Coupled to the VegET Evapotranspiration Model Using Spatially Explicit Comparisons with Independent Reference Data

    NASA Astrophysics Data System (ADS)

    Kovalskyy, V.; Henebry, G. M.; Roy, D. P.; Senay, G. B.

    2011-12-01

    Vegetation growing cycles have a profound influence on regional evapotranspiration regimes. The recently developed Event Driven Phenology Model (EDPM) is an empirical crop-specific phenology model with data assimilation capabilities. Deployed in prognostic mode, the EDPM uses weather forcing data to produce daily estimates of phenology coefficients; and in diagnostic mode a one-dimensional Kalman filter is used to adjust EDPM estimates with satellite normalized difference vegetation index (NDVI) retrievals. In this study the EDPM is coupled to the VegET model that uses the Penman-Monteith equation to calculate reference ET and a water balance model for water stress coefficients to derive daily actual evapotranspiration. The coupled models were run for the croplands of the U.S. Northern Great Plains for three annual growing seasons to derive 8-day total actual evapotranspiration (ETa) estimates at 0.05° spatial resolution. The models were driven by North American Land Data Assimilation System (NLDAS) weather forcing and parameterized using annual MODIS cropland cover maps. Regional validation of the modeled NDVI and ETa were undertaken by comparison with MODIS NDVI and MODIS ETa products respectively. The modeled NDVI had a median coefficient of determination (r2) of 0.83 and a root mean square error (RMSE) of 0.15 within study area. With the EDPM deployed in both prognostic and diagnostic modes, the modeled ETa had r2 of 0.75 and RMSE of about 25% of season average ETa per observation period. With small computational effort these results yield comparable accuracy to those from computationally complex models of ETa which require more parameterization. The performance of the coupling scheme demonstrates that the modeling approach is a promising avenue for regional application studies.

  11. High frequency scattering by a smooth coated cylinder simulated with generalized impedance boundary conditions

    NASA Technical Reports Server (NTRS)

    Syed, Hasnain H.; Volakis, John L.

    1991-01-01

    Rigorous uniform geometrical theory of diffraction (UGTD) diffraction coefficients are presented for a coated convex cylinder simulated with generalized impedance boundary conditions. In particular, ray solutions are obtained which remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. These involve new transition functions in place of the usual Fock-type integrals, characteristic to the impedance cylinder. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. As usual, the diffraction coefficients for the convex cylinder are obtained via a generalization of the corresponding ones for the circular cylinder.

  12. A General Relativistic Magnetohydrodynamics Simulation of Jet Formation with a State Transition

    NASA Technical Reports Server (NTRS)

    Nishikawa, K. I.; Richardson, G.; Koide, S.; Shibata, K.; Kudoh, T.; Hardee, P.; Fushman, G. J.

    2004-01-01

    We have performed the first fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity sim 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving shock wave. This evolution of jet-disk coupling suggests that the low/hard state of the jet system may switch to the high/soft state with a wind, as the accretion rate diminishes.

  13. Efficient classical simulation of matchgate circuits with generalized inputs and measurements

    NASA Astrophysics Data System (ADS)

    Brod, Daniel J.

    2016-06-01

    Matchgates are a restricted set of two-qubit gates known to be classically simulable under particular conditions. Specifically, if a circuit consists only of nearest-neighbor matchgates, an efficient classical simulation is possible if either (i) the input is a computational-basis state and the simulation requires computing probabilities of multiqubit outcomes (including also adaptive measurements) or (ii) if the input is an arbitrary product state, but the output of the circuit consists of a single qubit. In this paper we extend these results to show that matchgates are classically simulable even in the most general combination of these settings, namely, if the inputs are arbitrary product states, if the measurements are over arbitrarily many output qubits, and if adaptive measurements are allowed. This remains true even for arbitrary single-qubit measurements, albeit only in a weaker notion of classical simulation. These results make for an interesting contrast with other restricted models of computation, such as Clifford circuits or (bosonic) linear optics, where the complexity of simulation varies greatly under similar modifications.

  14. Well-posedness and generalized plane waves simulations of a 2D mode conversion model

    NASA Astrophysics Data System (ADS)

    Imbert-Gérard, Lise-Marie

    2015-12-01

    Certain types of electro-magnetic waves propagating in a plasma can undergo a mode conversion process. In magnetic confinement fusion, this phenomenon is very useful to heat the plasma, since it permits to transfer the heat at or near the plasma center. This work focuses on a mathematical model of wave propagation around the mode conversion region, from both theoretical and numerical points of view. It aims at developing, for a well-posed equation, specific basis functions to study a wave mode conversion process. These basis functions, called generalized plane waves, are intrinsically based on variable coefficients. As such, they are particularly adapted to the mode conversion problem. The design of generalized plane waves for the proposed model is described in detail. Their implementation within a discontinuous Galerkin method then provides numerical simulations of the process. These first 2D simulations for this model agree with qualitative aspects studied in previous works.

  15. Using MASS for AO simulations: a note on the comparison between MASS and Generalized SCIDAR techniques

    NASA Astrophysics Data System (ADS)

    Lombardi, G.; Sarazin, M.

    2016-01-01

    Recent studies on the comparison between the Multi Aperture Scintillation Sensor (MASS) and Generalized Scintillation Detection and Ranging (G-SCIDAR) profiler techniques have suggested significant discrepancies between the results delivered by the two instruments. MASS has been largely used in the recent site testing campaigns for the future next generation giant telescopes [i.e. the European Extremely Large Telescope, the Thirty Meter Telescope (TMT) and the Giant Magellan Telescope (GMT)] and is still used to monitor the conditions of world-class astronomical sites, as well as to deliver free atmosphere turbulence profiles to feed Adaptive Optics performance simulations. In this paper, we explore a different approach in the comparison between MASS and Generalized SCIDAR techniques with respect to previous studies, in order to provide a method for the use of the MASS data bases accumulated at European Southern Obseratory Paranal Observatory in Adaptive Optics simulations.

  16. Simulator Evaluation of Runway Incursion Prevention Technology for General Aviation Operations

    NASA Technical Reports Server (NTRS)

    Jones, Denise R.; Prinzel, Lawrence J., III

    2011-01-01

    A Runway Incursion Prevention System (RIPS) has been designed under previous research to enhance airport surface operations situation awareness and provide cockpit alerts of potential runway conflict, during transport aircraft category operations, in order to prevent runway incidents while also improving operations capability. This study investigated an adaptation of RIPS for low-end general aviation operations using a fixed-based simulator at the National Aeronautics and Space Administration (NASA) Langley Research Center (LaRC). The purpose of the study was to evaluate modified RIPS aircraft-based incursion detection algorithms and associated alerting and airport surface display concepts for low-end general aviation operations. This paper gives an overview of the system, simulation study, and test results.

  17. SIMULATION OF GENERAL ANESTHESIA ON THE "SIMMAN 3G" AND ITS EFFICIENCY.

    PubMed

    Potapov, A F; Matveev, A S; Ignatiev, V G; Ivanova, A A; Aprosimov, L A

    2015-01-01

    In recent years in medical educational process new innovative technologies are widely used with computer simulation, providing the reality of medical intervations and procedures. Practice-training teaching with using of simulation allows to improve the efficiency of learning material at the expense of creating imaginary professional activity and leading barring material to practical activity. The arm of the investigation is evaluation of the students training efficiency at the Medical Institute on the topic "General Anesthesia with using a modern simulation "SimMan 3 G". The material of the investigation is the results, carried out on the basis of the Centre of Practical skills and medical virtual educational technologies (Simulation Centre) at the Medical Institute of NEFU by M.K. Ammosov. The Object of the investigation was made up by 55 students of the third (3) course of the Faculty of General Medicine of the Medical Institute of NEFU. The investigation was hold during practical trainings (April-May 2014) of the General Surgery Department on the topic "General Anesthesia". A simulation practical course "General Anesthesia" consisted of 12 academic hours. Practical training was carried out using instruments, equipments and facilities to install anesthesia on the SimMan 3G with shooting the process and further discussions of the results. The methods of the investigations were the appreciation of students background knowledge before and after practical training (by 5 points scale) and the analysis of the results. The results of the investigation showed that before the practical course only 23 students (41.8%) had dot positive marks: "Good"--7 students (12.7%) and "Satisfactory"--16 (29.1%) students. The rest 22 (58.2%) students had bad results. The practical trainings using real instruments, equipments and facilities with imitation of installation of preparations for introductory anesthesia, main analgesics and muscle relaxants showed a patients reaction on the

  18. Simulation and flight evaluation of a heads-up display for general aviation

    NASA Technical Reports Server (NTRS)

    Harris, R. L., Sr.

    1974-01-01

    A landing-site indicator (LASI) has been devised as a relatively simple heads-up display to show the pilot the magnitude and direction of the aircraft's velocity vector superimposed on the pilot's view of the landing area. A total of 160 landings were performed in a fixed-base simulation program by four pilots with and without the LASI display. These tests showed the display to be of beneficial use in making the approaches more consistent. Some inferences were also made that the physical workload would also be less with its use. The pilots generally agreed that the LASI, as represented in the simulation was a useful landing aid. Additional pilot comments from preliminary flight tests of a breadboard LASI display unit tend to confirm the simulator results.

  19. Improved Carbohydrate Structure Generalization Scheme for (1)H and (13)C NMR Simulations.

    PubMed

    Kapaev, Roman R; Toukach, Philip V

    2015-07-21

    The improved Carbohydrate Structure Generalization Scheme has been developed for the simulation of (13)C and (1)H NMR spectra of oligo- and polysaccharides and their derivatives, including those containing noncarbohydrate constituents found in natural glycans. Besides adding the (1)H NMR calculations, we improved the accuracy and performance of prediction and optimized the mathematical model of the precision estimation. This new approach outperformed other methods of chemical shift simulation, including database-driven, neural net-based, and purely empirical methods and quantum-mechanical calculations at high theory levels. It can process structures with rarely occurring and noncarbohydrate constituents unsupported by the other methods. The algorithm is transparent to users and allows tracking used reference NMR data to original publications. It was implemented in the Glycan-Optimized Dual Empirical Spectrum Simulation (GODESS) web service, which is freely available at the platform of the Carbohydrate Structure Database (CSDB) project ( http://csdb.glycoscience.ru). PMID:26087011

  20. GOOSE 1. 4 -- Generalized Object-Oriented Simulation Environment user's manual

    SciTech Connect

    Nypaver, D.J. ); Abdalla, M.A. ); Guimaraes, L. , Sao Jose dos Campos, SP . Inst. de Estudos Avancados)

    1992-11-01

    The Generalized Object-Oriented Simulation Environment (GOOSE) is a new and innovative simulation tool that is being developed by the Simulation Group of the Advanced Controls Program at Oak Ridge National Laboratory. GOOSE is a fully interactive prototype software package that provides users with the capability of creating sophisticated mathematical models of physical systems. GOOSE uses an object-oriented approach to modeling and combines the concept of modularity (building a complex model easily from a collection of previously written components) with the additional features of allowing precompilation, optimization, and testing and validation of individual modules. Once a library of components has been defined and compiled, models can be built and modified without recompilation. This user's manual provides detailed descriptions of the structure and component features of GOOSE, along with a comprehensive example using a simplified model of a pressurized water reactor.

  1. GOOSE 1.4 -- Generalized Object-Oriented Simulation Environment user`s manual

    SciTech Connect

    Nypaver, D.J.; Abdalla, M.A.; Guimaraes, L.

    1992-11-01

    The Generalized Object-Oriented Simulation Environment (GOOSE) is a new and innovative simulation tool that is being developed by the Simulation Group of the Advanced Controls Program at Oak Ridge National Laboratory. GOOSE is a fully interactive prototype software package that provides users with the capability of creating sophisticated mathematical models of physical systems. GOOSE uses an object-oriented approach to modeling and combines the concept of modularity (building a complex model easily from a collection of previously written components) with the additional features of allowing precompilation, optimization, and testing and validation of individual modules. Once a library of components has been defined and compiled, models can be built and modified without recompilation. This user`s manual provides detailed descriptions of the structure and component features of GOOSE, along with a comprehensive example using a simplified model of a pressurized water reactor.

  2. Simulation and flight evaluation of a head-up landing aid for general aviation

    NASA Technical Reports Server (NTRS)

    Harris, R. L., Sr.; Goode, M. W.; Yenni, K. R.

    1978-01-01

    A head-up general aviation landing aid called a landing site indicator (LASI) was tested in a fixed-base, visual simulator and in an airplane to determine the effectiveness of the LASI. The display, which had a simplified format and method of implementation, presented to the pilot in his line of sight through the windshield a graphic representation of the airplane's velocity vector. In each testing model (simulation of flight), each of 4 pilots made 20 landing approaches with the LASI and 20 approaches without it. The standard deviations of approach and touchdown parameters were considered an indication of pilot consistency. Use of the LASI improved consistency and also reduced elevator, aileron, and rudder control activity. Pilots' comments indicated that the LASI reduced work load. An appendix is included with a discussion of the simulator effectiveness for visual flight tasks.

  3. Simulation of charge breeding of rubidium using Monte Carlo charge breeding code and generalized ECRIS model

    SciTech Connect

    Zhao, L.; Cluggish, B.; Kim, J. S.; Pardo, R.; Vondrasek, R.

    2010-02-15

    A Monte Carlo charge breeding code (MCBC) is being developed by FAR-TECH, Inc. to model the capture and charge breeding of 1+ ion beam in an electron cyclotron resonance ion source (ECRIS) device. The ECRIS plasma is simulated using the generalized ECRIS model which has two choices of boundary settings, free boundary condition and Bohm condition. The charge state distribution of the extracted beam ions is calculated by solving the steady state ion continuity equations where the profiles of the captured ions are used as source terms. MCBC simulations of the charge breeding of Rb+ showed good agreement with recent charge breeding experiments at Argonne National Laboratory (ANL). MCBC correctly predicted the peak of highly charged ion state outputs under free boundary condition and similar charge state distribution width but a lower peak charge state under the Bohm condition. The comparisons between the simulation results and ANL experimental measurements are presented and discussed.

  4. Towards Observational Astronomy of Jets in Active Galaxies from General Relativistic Magnetohydrodynamic Simulations

    NASA Astrophysics Data System (ADS)

    Anantua, Richard; Roger Blandford, Jonathan McKinney and Alexander Tchekhovskoy

    2016-01-01

    We carry out the process of "observing" simulations of active galactic nuclei (AGN) with relativistic jets (hereafter called jet/accretion disk/black hole (JAB) systems) from ray tracing between image plane and source to convolving the resulting images with a point spread function. Images are generated at arbitrary observer angle relative to the black hole spin axis by implementing spatial and temporal interpolation of conserved magnetohydrodynamic flow quantities from a time series of output datablocks from fully general relativistic 3D simulations. We also describe the evolution of simulations of JAB systems' dynamical and kinematic variables, e.g., velocity shear and momentum density, respectively, and the variation of these variables with respect to observer polar and azimuthal angles. We produce, at frequencies from radio to optical, fixed observer time intensity and polarization maps using various plasma physics motivated prescriptions for the emissivity function of physical quantities from the simulation output, and analyze the corresponding light curves. Our hypothesis is that this approach reproduces observed features of JAB systems such as superluminal bulk flow projections and quasi-periodic oscillations in the light curves more closely than extant stylized analytical models, e.g., cannonball bulk flows. Moreover, our development of user-friendly, versatile C++ routines for processing images of state-of-the-art simulations of JAB systems may afford greater flexibility for observing a wide range of sources from high power BL-Lacs to low power quasars (possibly with the same simulation) without requiring years of observation using multiple telescopes. Advantages of observing simulations instead of observing astrophysical sources directly include: the absence of a diffraction limit, panoramic views of the same object and the ability to freely track features. Light travel time effects become significant for high Lorentz factor and small angles between

  5. Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters

    SciTech Connect

    Dong, S.

    2015-02-15

    We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.

  6. Evaluation of the Event Driven Phenology Model Coupled with the VegET Evapotranspiration Model Through Comparisons with Reference Datasets in a Spatially Explicit Manner

    NASA Technical Reports Server (NTRS)

    Kovalskyy, V.; Henebry, G. M.; Adusei, B.; Hansen, M.; Roy, D. P.; Senay, G.; Mocko, D. M.

    2011-01-01

    A new model coupling scheme with remote sensing data assimilation was developed for estimation of daily actual evapotranspiration (ET). The scheme represents a mix of the VegET, a physically based model to estimate ET from a water balance, and an event driven phenology model (EDPM), where the EDPM is an empirically derived crop specific model capable of producing seasonal trajectories of canopy attributes. In this experiment, the scheme was deployed in a spatially explicit manner within the croplands of the Northern Great Plains. The evaluation was carried out using 2007-2009 land surface forcing data from the North American Land Data Assimilation System (NLDAS) and crop maps derived from remotely sensed data of NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). We compared the canopy parameters produced by the phenology model with normalized difference vegetation index (NDVI) data derived from the MODIS nadir bi-directional reflectance distribution function (BRDF) adjusted reflectance (NBAR) product. The expectations of the EDPM performance in prognostic mode were met, producing determination coefficient (r2) of 0.8 +/-.0.15. Model estimates of NDVI yielded root mean square error (RMSE) of 0.1 +/-.0.035 for the entire study area. Retrospective correction of canopy dynamics with MODIS NDVI brought the errors down to just below 10% of observed data range. The ET estimates produced by the coupled scheme were compared with ones from the MODIS land product suite. The expected r2=0.7 +/-.15 and RMSE = 11.2 +/-.4 mm per 8 days were met and even exceeded by the coupling scheme0 functioning in both prognostic and retrospective modes. Minor setbacks of the EDPM and VegET performance (r2 about 0.5 and additional 30 % of RMSR) were found on the peripheries of the study area and attributed to the insufficient EDPM training and to spatially varying accuracy of crop maps. Overall the experiment provided sufficient evidence of soundness and robustness of the EDPM and

  7. Continuity-based model interfacing for plant-wide simulation: a general approach.

    PubMed

    Volcke, Eveline I P; van Loosdrecht, Mark C M; Vanrolleghem, Peter A

    2006-08-01

    In plant-wide simulation studies of wastewater treatment facilities, often existing models from different origin need to be coupled. However, as these submodels are likely to contain different state variables, their coupling is not straightforward. The continuity-based interfacing method (CBIM) provides a general framework to construct model interfaces for models of wastewater systems, taking into account conservation principles. In this contribution, the CBIM approach is applied to study the effect of sludge digestion reject water treatment with a SHARON-Anammox process on a plant-wide scale. Separate models were available for the SHARON process and for the Anammox process. The Benchmark simulation model no. 2 (BSM2) is used to simulate the behaviour of the complete WWTP including sludge digestion. The CBIM approach is followed to develop three different model interfaces. At the same time, the generally applicable CBIM approach was further refined and particular issues when coupling models in which pH is considered as a state variable, are pointed out. PMID:16846629

  8. General relativistic magnetohydrodynamic simulations of accretion on to Sgr A*: how important are radiative losses?

    NASA Astrophysics Data System (ADS)

    Dibi, S.; Drappeau, S.; Fragile, P. C.; Markoff, S.; Dexter, J.

    2012-11-01

    We present general relativistic magnetohydrodynamic numerical simulations of the accretion flow around the supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*). The simulations include for the first time radiative cooling processes (synchrotron, bremsstrahlung and inverse Compton) self-consistently in the dynamics, allowing us to test the common simplification of ignoring all cooling losses in the modelling of Sgr A*. We confirm that for Sgr A*, neglecting the cooling losses is a reasonable approximation if the Galactic Centre is accreting below ˜10-8 M⊙ yr-1, i.e. M⊙<10-7M⊙ Edd . However, above this limit, we show that radiative losses should be taken into account as significant differences appear in the dynamics and the resulting spectra when comparing simulations with and without cooling. This limit implies that most nearby low-luminosity active galactic nuclei are in the regime where cooling should be taken into account. We further make a parameter study of axisymmetric gas accretion around the supermassive black hole at the Galactic Centre. This approach allows us to investigate the physics of gas accretion in general, while confronting our results with the well-studied and observed source, Sgr A*, as a test case. We confirm that the nature of the accretion flow and outflow is strongly dependent on the initial geometry of the magnetic field. For example, we find it difficult, even with very high spins, to generate powerful outflows from discs threaded with multiple, separate poloidal field loops.

  9. Generalized nonequilibrium vertex correction method in coherent medium theory for quantum transport simulation of disordered nanoelectronics

    NASA Astrophysics Data System (ADS)

    Yan, Jiawei; Ke, Youqi

    2016-07-01

    Electron transport properties of nanoelectronics can be significantly influenced by the inevitable and randomly distributed impurities/defects. For theoretical simulation of disordered nanoscale electronics, one is interested in both the configurationally averaged transport property and its statistical fluctuation that tells device-to-device variability induced by disorder. However, due to the lack of an effective method to do disorder averaging under the nonequilibrium condition, the important effects of disorders on electron transport remain largely unexplored or poorly understood. In this work, we report a general formalism of Green's function based nonequilibrium effective medium theory to calculate the disordered nanoelectronics. In this method, based on a generalized coherent potential approximation for the Keldysh nonequilibrium Green's function, we developed a generalized nonequilibrium vertex correction method to calculate the average of a two-Keldysh-Green's-function correlator. We obtain nine nonequilibrium vertex correction terms, as a complete family, to express the average of any two-Green's-function correlator and find they can be solved by a set of linear equations. As an important result, the averaged nonequilibrium density matrix, averaged current, disorder-induced current fluctuation, and averaged shot noise, which involve different two-Green's-function correlators, can all be derived and computed in an effective and unified way. To test the general applicability of this method, we applied it to compute the transmission coefficient and its fluctuation with a square-lattice tight-binding model and compared with the exact results and other previously proposed approximations. Our results show very good agreement with the exact results for a wide range of disorder concentrations and energies. In addition, to incorporate with density functional theory to realize first-principles quantum transport simulation, we have also derived a general form of

  10. Multiple processor accelerator for logic simulation

    SciTech Connect

    Catlin, G.M.

    1989-10-17

    This patent describes a computer system coupled to a plurality of users for implementing an event driven algorithm of each of the users. It comprises: a master processor coupled to the users for providing overall control of the computer system and executing the event driven algorithm of each of the users, the master processor further including a master memory; a unidirectional ring bus coupled to the master processor; a plurality of processor modules; an interprocessor bus coupled to the plurality of processors within the module for transferring the simulation data among the processors; and an interface means.

  11. General relativistic N-body simulations in the weak field limit

    NASA Astrophysics Data System (ADS)

    Adamek, Julian; Daverio, David; Durrer, Ruth; Kunz, Martin

    2013-11-01

    We develop a formalism for general relativistic N-body simulations in the weak field regime, suitable for cosmological applications. The problem is kept tractable by retaining the metric perturbations to first order, the first derivatives to second order, and second derivatives to all orders, thus taking into account the most important nonlinear effects of Einstein gravity. It is also expected that any significant “backreaction” should appear at this order. We show that the simulation scheme is feasible in practice by implementing it for a plane-symmetric situation and running two test cases, one with only cold dark matter, and one which also includes a cosmological constant. For these plane-symmetric situations, the deviations from the usual Newtonian N-body simulations remain small and, apart from a nontrivial correction to the background, can be accurately estimated within the Newtonian framework. The correction to the background scale factor, which is a genuine backreaction effect, can be robustly obtained with our algorithm. Our numerical approach is also naturally suited for the inclusion of extra relativistic fields and thus for dark energy or modified gravity simulations.

  12. Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born

    PubMed Central

    2012-01-01

    We present an implementation of generalized Born implicit solvent all-atom classical molecular dynamics (MD) within the AMBER program package that runs entirely on CUDA enabled NVIDIA graphics processing units (GPUs). We discuss the algorithms that are used to exploit the processing power of the GPUs and show the performance that can be achieved in comparison to simulations on conventional CPU clusters. The implementation supports three different precision models in which the contributions to the forces are calculated in single precision floating point arithmetic but accumulated in double precision (SPDP), or everything is computed in single precision (SPSP) or double precision (DPDP). In addition to performance, we have focused on understanding the implications of the different precision models on the outcome of implicit solvent MD simulations. We show results for a range of tests including the accuracy of single point force evaluations and energy conservation as well as structural properties pertainining to protein dynamics. The numerical noise due to rounding errors within the SPSP precision model is sufficiently large to lead to an accumulation of errors which can result in unphysical trajectories for long time scale simulations. We recommend the use of the mixed-precision SPDP model since the numerical results obtained are comparable with those of the full double precision DPDP model and the reference double precision CPU implementation but at significantly reduced computational cost. Our implementation provides performance for GB simulations on a single desktop that is on par with, and in some cases exceeds, that of traditional supercomputers. PMID:22582031

  13. Nonparametric simulation-based statistics for detecting linkage in general pedigrees

    SciTech Connect

    Davis, S.; Schroeder, M.; Weeks, D.E.; Goldin, L.R.

    1996-04-01

    We present here four nonparametric statistics for linkage analysis that test whether pairs of affected relatives share marker alleles more often than expected. These statistics are based on simulating the null distribution of a given statistic conditional on the unaffecteds` marker genotypes. Each statistic uses a different measure of marker sharing: the SimAPM statistic uses the simulation-based affected-pedigree-member measure based on identity-by-state (IBS) sharing. The SimKIN (kinship) measure is 1.0 for identity-by-descent (IBD) sharing, 0.0 for no IBD sharing, and the kinship coefficient when the IBD status is ambiguous. The simulation-based IBD (SimIBD) statistic uses a recursive algorithm to determine the probability of two affecteds sharing a specific allele IBD. The SimISO statistic is identical to SimIBD, except that it also measures marker similarity between unaffected pairs. We evaluated our statistics on data simulated under different two-locus disease models, comparing our results to those obtained with several other nonparametric statistics. Use of IBD information produces dramatic increases in power over the SimAPM method, which uses only IBS information. The power of our best statistic in most cases meets or exceeds the power of the other nonparametric statistics. Furthermore, our statistics perform comparisons between all affected relative pairs within general pedigrees and are not restricted to sib pairs or nuclear families. 32 refs., 5 figs., 6 tabs.

  14. Dust Emissions, Transport, and Deposition Simulated with the NASA Finite-Volume General Circulation Model

    NASA Technical Reports Server (NTRS)

    Colarco, Peter; daSilva, Arlindo; Ginoux, Paul; Chin, Mian; Lin, S.-J.

    2003-01-01

    Mineral dust aerosols have radiative impacts on Earth's atmosphere, have been implicated in local and regional air quality issues, and have been identified as vectors for transporting disease pathogens and bringing mineral nutrients to terrestrial and oceanic ecosystems. We present for the first time dust simulations using online transport and meteorological analysis in the NASA Finite-Volume General Circulation Model (FVGCM). Our dust formulation follows the formulation in the offline Georgia Institute of Technology-Goddard Global Ozone Chemistry Aerosol Radiation and Transport Model (GOCART) using a topographical source for dust emissions. We compare results of the FVGCM simulations with GOCART, as well as with in situ and remotely sensed observations. Additionally, we estimate budgets of dust emission and transport into various regions.

  15. 2D simulations based on general time-dependent reciprocal relation for LFEIT.

    PubMed

    Karadas, Mursel; Gencer, Nevzat Guneri

    2015-08-01

    Lorentz field electrical impedance tomography (LFEIT) is a newly proposed technique for imaging the conductivity of the tissues by measuring the electromagnetic induction under the ultrasound pressure field. In this paper, the theory and numerical simulations of the LFEIT are reported based on the general time dependent formulation. In LFEIT, a phased array ultrasound probe is used to introduce a current distribution inside a conductive body. The velocity current occurs, due to the movement of the conductive particles under a static magnetic field. In order to sense this current, a receiver coil configuration that surrounds the volume conductor is utilized. Finite Element Method (FEM) is used to carry out the simulations of LFEIT. It is shown that, LFEIT can be used to reconstruct the conductivity even up to 50% perturbation in the initial conductivity distribution. PMID:26736569

  16. GENASIS: General Astrophysical Simulation System. I. Refinable Mesh and Nonrelativistic Hydrodynamics

    NASA Astrophysics Data System (ADS)

    Cardall, Christian Y.; Budiardja, Reuben D.; Endeve, Eirik; Mezzacappa, Anthony

    2014-02-01

    GenASiS (General Astrophysical Simulation System) is a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper—the first in a series—demonstrates a centrally refined coordinate patch suitable for gravitational collapse and documents methods for compressible nonrelativistic hydrodynamics. We benchmark the hydrodynamics capabilities of GenASiS against many standard test problems; the results illustrate the basic competence of our implementation, demonstrate the strengths and limitations of the HLLC relative to the HLL Riemann solver in a number of interesting cases, and provide preliminary indications of the code's ability to scale and to function with cell-by-cell fixed-mesh refinement.

  17. Automated procedure for developing hybrid computer simulations of turbofan engines. Part 1: General description

    NASA Technical Reports Server (NTRS)

    Szuch, J. R.; Krosel, S. M.; Bruton, W. M.

    1982-01-01

    A systematic, computer-aided, self-documenting methodology for developing hybrid computer simulations of turbofan engines is presented. The methodology that is pesented makes use of a host program that can run on a large digital computer and a machine-dependent target (hybrid) program. The host program performs all the calculations and data manipulations that are needed to transform user-supplied engine design information to a form suitable for the hybrid computer. The host program also trims the self-contained engine model to match specified design-point information. Part I contains a general discussion of the methodology, describes a test case, and presents comparisons between hybrid simulation and specified engine performance data. Part II, a companion document, contains documentation, in the form of computer printouts, for the test case.

  18. User's guide for a general purpose dam-break flood simulation model (K-634)

    USGS Publications Warehouse

    Land, Larry F.

    1981-01-01

    An existing computer program for simulating dam-break floods for forecast purposes has been modified with an emphasis on general purpose applications. The original model was formulated, developed and documented by the National Weather Service. This model is based on the complete flow equations and uses a nonlinear implicit finite-difference numerical method. The first phase of the simulation routes a flood wave through the reservoir and computes an outflow hydrograph which is the sum of the flow through the dam 's structures and the gradually developing breach. The second phase routes this outflow hydrograph through the stream which may be nonprismatic and have segments with subcritical or supercritical flow. The results are discharge and stage hydrographs at the dam as well as all of the computational nodes in the channel. From these hydrographs, peak discharge and stage profiles are tabulated. (USGS)

  19. A Generalized Fast Frequency Sweep Algorithm for Coupled Circuit-EM Simulations

    SciTech Connect

    Rockway, J D; Champagne, N J; Sharpe, R M; Fasenfest, B

    2004-01-14

    Frequency domain techniques are popular for analyzing electromagnetics (EM) and coupled circuit-EM problems. These techniques, such as the method of moments (MoM) and the finite element method (FEM), are used to determine the response of the EM portion of the problem at a single frequency. Since only one frequency is solved at a time, it may take a long time to calculate the parameters for wideband devices. In this paper, a fast frequency sweep based on the Asymptotic Wave Expansion (AWE) method is developed and applied to generalized mixed circuit-EM problems. The AWE method, which was originally developed for lumped-load circuit simulations, has recently been shown to be effective at quasi-static and low frequency full-wave simulations. Here it is applied to a full-wave MoM solver, capable of solving for metals, dielectrics, and coupled circuit-EM problems.

  20. General Relativistic Magnetohydrodynamics Simulations of Tilted Black Hole Accretion Flows and Their Radiative Properties

    NASA Astrophysics Data System (ADS)

    Shiokawa, Hotaka; Gammie, C. F.; Dolence, J.; Noble, S. C.

    2013-01-01

    We perform global General Relativistic Magnetohydrodynamics (GRMHD) simulations of non-radiative, magnetized disks that are initially tilted with respect to the black hole's spin axis. We run the simulations with different size and tilt angle of the tori for 2 different resolutions. We also perform radiative transfer using Monte Carlo based code that includes synchrotron emission, absorption and Compton scattering to obtain spectral energy distribution and light curves. Similar work was done by Fragile et al. (2007) and Dexter & Fragile (2012) to model the super massive black hole SgrA* with tilted accretion disks. We compare our results of fully conservative hydrodynamic code and spectra that include X-ray, with their results.

  1. GENERAL-RELATIVISTIC SIMULATIONS OF THREE-DIMENSIONAL CORE-COLLAPSE SUPERNOVAE

    SciTech Connect

    Ott, Christian D.; Abdikamalov, Ernazar; Moesta, Philipp; Haas, Roland; Drasco, Steve; O'Connor, Evan P.; Reisswig, Christian; Meakin, Casey A.; Schnetter, Erik

    2013-05-10

    We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27 M{sub Sun} star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a three-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27 M{sub Sun} progenitor was studied in 2D by Mueller et al., who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.

  2. SIMPSON: A General Simulation Program for Solid-State NMR Spectroscopy

    NASA Astrophysics Data System (ADS)

    Bak, Mads; Rasmussen, Jimmy T.; Nielsen, Niels Chr.

    2000-12-01

    A computer program for fast and accurate numerical simulation of solid-state NMR experiments is described. The program is designed to emulate a NMR spectrometer by letting the user specify high-level NMR concepts such as spin systems, nuclear spin interactions, RF irradiation, free precession, phase cycling, coherence-order filtering, and implicit/explicit acquisition. These elements are implemented using the Tcl scripting language to ensure a minimum of programming overhead and direct interpretation without the need for compilation, while maintaining the flexibility of a full-featured programming language. Basicly, there are no intrinsic limitations to the number of spins, types of interactions, sample conditions (static or spinning, powders, uniaxially oriented molecules, single crystals, or solutions), and the complexity or number of spectral dimensions for the pulse sequence. The applicability ranges from simple 1D experiments to advanced multiple-pulse and multiple-dimensional experiments, series of simulations, parameter scans, complex data manipulation/visualization, and iterative fitting of simulated to experimental spectra. A major effort has been devoted to optimizing the computation speed using state-of-the-art algorithms for the time-consuming parts of the calculations implemented in the core of the program using the C programming language. Modification and maintenance of the program are facilitated by releasing the program as open source software (General Public License) currently at http://nmr.imsb.au.dk. The general features of the program are demonstrated by numerical simulations of various aspects for REDOR, rotational resonance, DRAMA, DRAWS, HORROR, C7, TEDOR, POST-C7, CW decoupling, TPPM, F-SLG, SLF, SEMA-CP, PISEMA, RFDR, QCPMG-MAS, and MQ-MAS experiments.

  3. SIMPSON: A general simulation program for solid-state NMR spectroscopy

    NASA Astrophysics Data System (ADS)

    Bak, Mads; Rasmussen, Jimmy T.; Nielsen, Niels Chr.

    2011-12-01

    A computer program for fast and accurate numerical simulation of solid-state NMR experiments is described. The program is designed to emulate a NMR spectrometer by letting the user specify high-level NMR concepts such as spin systems, nuclear spin interactions, RF irradiation, free precession, phase cycling, coherence-order filtering, and implicit/explicit acquisition. These elements are implemented using the Tel scripting language to ensure a minimum of programming overhead and direct interpretation without the need for compilation, while maintaining the flexibility of a full-featured programming language. Basicly, there are no intrinsic limitations to the number of spins, types of interactions, sample conditions (static or spinning, powders, uniaxially oriented molecules, single crystals, or solutions), and the complexity or number of spectral dimensions for the pulse sequence. The applicability ranges from simple ID experiments to advanced multiple-pulse and multiple-dimensional experiments, series of simulations, parameter scans, complex data manipulation/visualization, and iterative fitting of simulated to experimental spectra. A major effort has been devoted to optimizing the computation speed using state-of-the-art algorithms for the time-consuming parts of the calculations implemented in the core of the program using the C programming language. Modification and maintenance of the program are facilitated by releasing the program as open source software (General Public License) currently at http://nmr.imsb.au.dk. The general features of the program are demonstrated by numerical simulations of various aspects for REDOR, rotational resonance, DRAMA, DRAWS, HORROR, C7, TEDOR, POST-C7, CW decoupling, TPPM, F-SLG, SLF, SEMA-CP, PISEMA, RFDR, QCPMG-MAS, and MQ-MAS experiments.

  4. General-relativistic Simulations of Three-dimensional Core-collapse Supernovae

    NASA Astrophysics Data System (ADS)

    Ott, Christian D.; Abdikamalov, Ernazar; Mösta, Philipp; Haas, Roland; Drasco, Steve; O'Connor, Evan P.; Reisswig, Christian; Meakin, Casey A.; Schnetter, Erik

    2013-05-01

    We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27 M ⊙ star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a three-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27 M ⊙ progenitor was studied in 2D by Müller et al., who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.

  5. An in-flight simulation of lateral control nonlinearities. [for general aviation aircraft

    NASA Technical Reports Server (NTRS)

    Ellis, D. R.; Tilak, N. W.

    1975-01-01

    An in-flight simulation program was conducted to explore, in a generalized way, the influence of spoiler-type roll-control nonlinearities on handling qualities. The roll responses studied typically featured a dead zone or very small effectiveness for small control inputs, a very high effectiveness for mid-range deflections, and low effectiveness again for large inputs. A linear force gradient with no detectable breakout force was provided. Given otherwise good handling characteristics, it was found that moderate nonlinearities of the types tested might yield acceptable roll control, but the best level of handling qualities is obtained with linear, aileron-like control.

  6. Terahertz spectroscopic polarimetry of generalized anisotropic media composed of Archimedean spiral arrays: Experiments and simulations

    NASA Astrophysics Data System (ADS)

    Aschaffenburg, Daniel J.; Williams, Michael R. C.; Schmuttenmaer, Charles A.

    2016-05-01

    Terahertz time-domain spectroscopic polarimetry has been used to measure the polarization state of all spectral components in a broadband THz pulse upon transmission through generalized anisotropic media consisting of two-dimensional arrays of lithographically defined Archimedean spirals. The technique allows a full determination of the frequency-dependent, complex-valued transmission matrix and eigenpolarizations of the spiral arrays. Measurements were made on a series of spiral array orientations. The frequency-dependent transmission matrix elements as well as the eigenpolarizations were determined, and the eigenpolarizations were found be to elliptically corotating, as expected from their symmetry. Numerical simulations are in quantitative agreement with measured spectra.

  7. Mars atmospheric dynamics as simulated by the NASA AMES General Circulation Model. II - Transient baroclinic eddies

    NASA Astrophysics Data System (ADS)

    Barnes, J. R.; Pollack, J. B.; Haberle, R. M.; Leovy, C. B.; Zurek, R. W.; Lee, H.; Schaeffer, J.

    1993-02-01

    A large set of experiments performed with the NASA Ames Mars General Circulation Model is analyzed to determine the properties, structure, and dynamics of the simulated transient baroclinic eddies. There is strong transient baroclinic eddy activity in the extratropics of the Northern Hemisphere during the northern autumn, winter, and spring seasons. The eddy activity remains strong for very large dust loadings, though it shifts northward. The eastward propagating eddies are characterized by zonal wavenumbers of 1-4 and periods of about 2-10 days. The properties of the GCM baroclinic eddies in the northern extratropics are compared in detail with analogous properties inferred from Viking Lander meteorology observations.

  8. Reconstruction of bremsstrahlung spectra from attenuation data using generalized simulated annealing.

    PubMed

    Menin, O H; Martinez, A S; Costa, A M

    2016-05-01

    A generalized simulated annealing algorithm, combined with a suitable smoothing regularization function is used to solve the inverse problem of X-ray spectrum reconstruction from attenuation data. The approach is to set the initial acceptance and visitation temperatures and to standardize the terms of objective function to automate the algorithm to accommodate different spectra ranges. Experiments with both numerical and measured attenuation data are presented. Results show that the algorithm reconstructs spectra shapes accurately. It should be noted that in this algorithm, the regularization function was formulated to guarantee a smooth spectrum, thus, the presented technique does not apply to X-ray spectrum where characteristic radiation are present. PMID:26943902

  9. General purpose computational tools for simulation and analysis of medium-energy backscattering spectra

    NASA Astrophysics Data System (ADS)

    Weller, Robert A.

    1999-06-01

    This paper describes a suite of computational tools for general-purpose ion-solid calculations, which has been implemented in the platform-independent computational environment Mathematica®. Although originally developed for medium energy work (beam energies < 300 keV), they are suitable for general, classical, non-relativistic calculations. Routines are available for stopping power, Rutherford and Lenz-Jensen (screened) cross sections, sputtering yields, small-angle multiple scattering, and back-scattering-spectrum simulation and analysis. Also included are a full range of supporting functions, as well as easily accessible atomic mass and other data on all the stable isotopes in the periodic table. The functions use common calling protocols, recognize elements and isotopes by symbolic names and, wherever possible, return symbolic results for symbolic inputs, thereby facilitating further computation. A new paradigm for the representation of backscattering spectra is introduced.

  10. Development of generalized mapping tools to improve implementation of data driven computer simulations (04-ERD-083)

    SciTech Connect

    Ramirez, A; Pasyanos, M; Franz, G A

    2004-09-17

    The Stochastic Engine (SE) is a data driven computer simulation tool for predicting the characteristics of complex systems. The SE integrates accurate simulators with the Monte Carlo Markov Chain (MCMC) approach (a stochastic inverse technique) to identify alternative models that are consistent with available data and ranks these alternatives according to their probabilities. Implementation of the SE is currently cumbersome owing to the need to customize the pre-processing and processing steps that are required for a specific application. This project widens the applicability of the Stochastic Engine by generalizing some aspects of the method (i.e. model-to-data transformation types, configuration, model representation). We have generalized several of the transformations that are necessary to match the observations to proposed models. These transformations are sufficiently general not to pertain to any single application. This approach provides a framework that increases the efficiency of the SE implementation. The overall goal is to reduce response time and make the approach as ''plug-and-play'' as possible, and will result in the rapid accumulation of new data types for a host of both earth science and non-earth science problems. When adapting the SE approach to a specific application, there are various pre-processing and processing steps that are typically needed to run a specific problem. Many of these steps are common to a wide variety of specific applications. Here we list and describe several data transformations that are common to a variety of subsurface inverse problems. A subset of these steps have been developed in a generalized form such that they could be used with little or no modifications in a wide variety of specific applications. This work was funded by the LDRD Program (tracking number 04-ERD-083).

  11. The impact of a realistic vertical dust distribution on the simulation of the Martian General Circulation

    NASA Astrophysics Data System (ADS)

    Guzewich, Scott D.; Toigo, Anthony D.; Richardson, Mark I.; Newman, Claire E.; Talaat, Elsayed R.; Waugh, Darryn W.; McConnochie, Timothy H.

    2013-05-01

    Limb-scanning observations with the Mars Climate Sounder and Thermal Emission Spectrometer (TES) have identified discrete layers of enhanced dust opacity well above the boundary layer and a mean vertical structure of dust opacity very different from the expectation of well-mixed dust in the lowest 1-2 scale heights. To assess the impact of this vertical dust opacity profile on atmospheric properties, we developed a TES limb-scan observation-based three-dimensional and time-evolving dust climatology for use in forcing general circulation models (GCMs). We use this to force the MarsWRF GCM and compare with simulations that use a well-mixed (Conrath-ν) vertical dust profile and Mars Climate Database version 4 (MCD) horizontal distribution dust opacity forcing function. We find that simulated temperatures using the TES-derived forcing yield a 1.18 standard deviation closer match to TES temperature retrievals than a MarsWRF simulation using MCD forcing. The climatological forcing yields significant changes to many large-scale features of the simulated atmosphere. Notably the high-latitude westerly jet speeds are 10-20 m/s higher, polar warming collar temperatures are 20-30 K warmer near northern winter solstice and tilted more strongly poleward, the middle and lower atmospheric meridional circulations are partially decoupled, the migrating diurnal tide exhibits destructive interference and is weakened by 50% outside of equinox, and the southern hemisphere wave number 1 stationary wave is strengthened by up to 4 K (45%). We find the vertical dust distribution is an important factor for Martian lower and middle atmospheric thermal structure and circulation that cannot be neglected in analysis and simulation of the Martian atmosphere.

  12. A Novel Approach for Modeling Chemical Reaction in Generalized Fluid System Simulation Program

    NASA Technical Reports Server (NTRS)

    Sozen, Mehmet; Majumdar, Alok

    2002-01-01

    The Generalized Fluid System Simulation Program (GFSSP) is a computer code developed at NASA Marshall Space Flight Center for analyzing steady state and transient flow rates, pressures, temperatures, and concentrations in a complex flow network. The code, which performs system level simulation, can handle compressible and incompressible flows as well as phase change and mixture thermodynamics. Thermodynamic and thermophysical property programs, GASP, WASP and GASPAK provide the necessary data for fluids such as helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, water, a hydrogen, isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, several refrigerants, nitrogen trifluoride and ammonia. The program which was developed out of need for an easy to use system level simulation tool for complex flow networks, has been used for the following purposes to name a few: Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump Secondary Flow Circuits, Axial Thrust Balance of the Fastrac Engine Turbopump, Pressurized Propellant Feed System for the Propulsion Test Article at Stennis Space Center, X-34 Main Propulsion System, X-33 Reaction Control System and Thermal Protection System, and International Space Station Environmental Control and Life Support System design. There has been an increasing demand for implementing a combustion simulation capability into GFSSP in order to increase its system level simulation capability of a liquid rocket propulsion system starting from the propellant tanks up to the thruster nozzle for spacecraft as well as launch vehicles. The present work was undertaken for addressing this need. The chemical equilibrium equations derived from the second law of thermodynamics and the energy conservation equation derived from the first law of thermodynamics are solved simultaneously by a Newton-Raphson method. The numerical scheme was implemented as a User

  13. Simulation of reactive nanolaminates using reduced models: III. Ingredients for a general multidimensional formulation

    SciTech Connect

    Salloum, Maher; Knio, Omar M.

    2010-06-15

    A transient multidimensional reduced model is constructed for the simulation of reaction fronts in Ni/Al multilayers. The formulation is based on the generalization of earlier methodologies developed for quasi-1D axial and normal propagation, specifically by adapting the reduced formalism for atomic mixing and heat release. This approach enables us to focus on resolving the thermal front structure, whose evolution is governed by thermal diffusion and heat release. A mixed integration scheme is used for this purpose, combining an extended-stability, Runge-Kutta-Chebychev (RKC) integration of the diffusion term with exact treatment of the chemical source term. Thus, a detailed description of atomic mixing within individual layers is avoided, which enables transient modeling of the reduced equations of motion in multiple dimensions. Two-dimensional simulations are first conducted of front propagation in composites combining two bilayer periods. Results are compared with the experimental measurements of Knepper et al., which reveal that the reaction velocity can depend significantly on layering frequency. The comparison indicates that, using a concentration-dependent conductivity model, the transient 2D computations can reasonably reproduce the experimental behavior. Additional tests are performed based on 3D computations of surface initiated reactions. Comparison of computed predictions with laser ignition measurements indicates that the computations provide reasonable estimates of ignition thresholds. A detailed discussion is finally provided of potential generalizations and associated hurdles. (author)

  14. Variable-resolution frameworks for the simulation of tropical cyclones in global atmospheric general circulation models

    NASA Astrophysics Data System (ADS)

    Zarzycki, Colin

    The ability of atmospheric General Circulation Models (GCMs) to resolve tropical cyclones in the climate system has traditionally been difficult. The challenges include adequately capturing storms which are small in size relative to model grids and the fact that key thermodynamic processes require a significant level of parameterization. At traditional GCM grid spacings of 50-300 km tropical cyclones are severely under-resolved, if not completely unresolved. This thesis explores a variable-resolution global model approach that allows for high spatial resolutions in areas of interest, such as low-latitude ocean basins where tropical cyclogenesis occurs. Such GCM designs with multi-resolution meshes serve to bridge the gap between globally-uniform grids and limited area models and have the potential to become a future tool for regional climate assessments. A statically-nested, variable-resolution option has recently been introduced into the Department of Energy/National Center for Atmospheric Research (DoE/NCAR) Community Atmosphere Model's (CAM) Spectral Element (SE) dynamical core. Using an idealized tropical cyclone test, variable-resolution meshes are shown to significantly lessen computational requirements in regional GCM studies. Furthermore, the tropical cyclone simulations are free of spurious numerical errors at the resolution interfaces. Utilizing aquaplanet simulations as an intermediate test between idealized simulations and fully-coupled climate model runs, climate statistics within refined patches are shown to be well-matched to globally-uniform simulations of the same grid spacing. Facets of the CAM version 4 (CAM4) subgrid physical parameterizations are likely too scale sensitive for variable-resolution applications, but the newer CAM5 package is vastly improved in performance at multiple grid spacings. Multi-decadal simulations following 'Atmospheric Model Intercomparison Project' protocols have been conducted with variable-resolution grids. Climate

  15. Radioscience simulations in general relativity and in alternative theories of gravity

    NASA Astrophysics Data System (ADS)

    Hees, A.; Lamine, B.; Reynaud, S.; Jaekel, M.-T.; Le Poncin-Lafitte, C.; Lainey, V.; Füzfa, A.; Courty, J.-M.; Dehant, V.; Wolf, P.

    2012-12-01

    This paper deals with tests of general relativity (GR) in the Solar System using tracking observables from planetary spacecraft. We present a new software that simulates the Range and Doppler signals resulting from a given spacetime metric. This flexible approach allows one to perform simulations in GR as well as in alternative metric theories of gravity. The outputs of this software provide templates of anomalous residuals that should show up in real data if the underlying theory of gravity is not GR. Those templates can be used to give a rough estimation of constraints on additional parameters entering alternative theory of gravity and also signatures that can be searched for in data from past or future space missions aiming at testing gravitational laws in the Solar System. As an application of the potentiality of this software, we present some simulations performed for Cassini-like mission in post-Einsteinian gravity and in the context of MOND external field effect. We derive signatures arising from these alternative theories of gravity and estimate expected amplitudes of the anomalous residuals.

  16. Numerical simulation of the general circulation of the atmosphere of Titan.

    PubMed

    Hourdin, F; Talagrand, O; Sadourny, R; Courtin, R; Gautier, D; McKay, C P

    1995-10-01

    The atmospheric circulation of Titan is investigated with a general circulation model. The representation of the large-scale dynamics is based on a grid point model developed and used at Laboratoire de Météorologie Dynamique for climate studies. The code also includes an accurate representation of radiative heating and cooling by molecular gases and haze as well as a parametrization of the vertical turbulent mixing of momentum and potential temperature. Long-term simulations of the atmospheric circulation are presented. Starting from a state of rest, the model spontaneously produces a strong superrotation with prograde equatorial winds (i.e., in the same sense as the assumed rotation of the solid body) increasing from the surface to reach 100 m sec-1 near the 1-mbar pressure level. Those equatorial winds are in very good agreement with some indirect observations, especially those of the 1989 occultation of Star 28-Sgr by Titan. On the other hand, the model simulates latitudinal temperature contrasts in the stratosphere that are significantly weaker than those observed by Voyager 1 which, we suggest, may be partly due to the nonrepresentation of the spatial and temporal variations of the abundances of molecular species and haze. We present diagnostics of the simulated atmospheric circulation underlying the importance of the seasonal cycle and a tentative explanation for the creation and maintenance of the atmospheric superrotation based on a careful angular momentum budget. PMID:11538593

  17. General relativistic simulations of slowly rotating, magnetized stars: A perturbative metric approach

    NASA Astrophysics Data System (ADS)

    Etienne, Zachariah; Liu, Y. T.; Shapiro, S.

    2007-04-01

    Understanding the role general relativistic magnetohydrodynamic (GRMHD) effects play in the evolution of nascent neutron stars is a problem at the forefront of theoretical astrophysics. To this end, we performed long-term (˜10^4 M) axisymmetric simulations of differentially rotating magnetized neutron stars in the slow-rotation, weak magnetic field limit using a dynamically updated perturbative metric evolution technique. Although the perturbative metric approach yields results comparable to those obtained via a nonperturbative (BSSN) metric evolution technique, simulations performed with the perturbative metric solver require about 1/4 the computational resources at a given resolution. This computational efficiency enabled us to observe and analyze the effects of magnetic braking and the magnetorotational instability (MRI) at very high resolution. Our GRMHD simulations demonstrate that (1) MRI is not observed unless the estimated fastest-growing mode wavelength is resolved by >˜ 10 gridpoints; (2) as resolution is improved, the MRI growth rate converges, but due to the small-scale nature of MRI-induced turbulence, the maximum growth amplitude increases, but does not exhibit convergence, even at the highest resolution; and (3) independent of resolution, magnetic braking drives the star toward uniform rotation as energy is sapped from differential rotation by winding magnetic fields.

  18. Towards a General Turbulence Model for Planetary Boundary Layers Based on Direct Statistical Simulation

    NASA Astrophysics Data System (ADS)

    Marston, Brad; Fox-Kemper, Baylor; Skitka, Joe

    Sub-grid turbulence models for planetary boundary layers are typically constructed additively, starting with local flow properties and including non-local (KPP) or higher order (Mellor-Yamada) parameters until a desired level of predictive capacity is achieved or a manageable threshold of complexity is surpassed. Such approaches are necessarily limited in general circumstances, like global circulation models, by their being optimized for particular flow phenomena. By using direct statistical simulation (DSS) that is based upon expansion in equal-time cumulants we offer the prospect of a turbulence model and an investigative tool that is equally applicable to all flow types and able to take advantage of the wealth of nonlocal information in any flow. We investigate the feasibility of a second-order closure (CE2) by performing simulations of the ocean boundary layer in a quasi-linear approximation for which CE2 is exact. As oceanographic examples, wind-driven Langmuir turbulence and thermal convection are studied by comparison of the statistics of quasi-linear and fully nonlinear simulations. We also characterize the computational advantages and physical uncertainties of CE2 defined on a reduced basis determined via proper orthogonal decomposition (POD) of the flow fields. Supported in part by NSF DMR-1306806.

  19. Martian atmospheric gravity waves simulated by a high-resolution general circulation model

    NASA Astrophysics Data System (ADS)

    Kuroda, Takeshi; Yiǧit, Erdal; Medvedev, Alexander S.; Hartogh, Paul

    2016-07-01

    Gravity waves (GWs) significantly affect temperature and wind fields in the Martian middle and upper atmosphere. They are also one of the observational targets of the MAVEN mission. We report on the first simulations with a high-resolution general circulation model (GCM) and present a global distributions of small-scale GWs in the Martian atmosphere. The simulated GW-induced temperature variances are in a good agreement with available radio occultation data in the lower atmosphere between 10 and 30 km. For the northern winter solstice, the model reveals a latitudinal asymmetry with stronger wave generation in the winter hemisphere and two distinctive sources of GWs: mountainous regions and the meandering winter polar jet. Orographic GWs are filtered upon propagating upward, and the mesosphere is primarily dominated by harmonics with faster horizontal phase velocities. Wave fluxes are directed mainly against the local wind. GW dissipation in the upper mesosphere generates a body force per unit mass of tens of m s^{-1} per Martian solar day (sol^{-1}), which tends to close the simulated jets. The results represent a realistic surrogate for missing observations, which can be used for constraining GW parameterizations and validating GCMs.

  20. A comparison between general circulation model simulations using two sea surface temperature datasets for January 1979

    NASA Technical Reports Server (NTRS)

    Ose, Tomoaki; Mechoso, Carlos; Halpern, David

    1994-01-01

    Simulations with the UCLA atmospheric general circulation model (AGCM) using two different global sea surface temperature (SST) datasets for January 1979 are compared. One of these datasets is based on Comprehensive Ocean-Atmosphere Data Set (COADS) (SSTs) at locations where there are ship reports, and climatology elsewhere; the other is derived from measurements by instruments onboard NOAA satellites. In the former dataset (COADS SST), data are concentrated along shipping routes in the Northern Hemisphere; in the latter dataset High Resolution Infrared Sounder (HIRS SST), data cover the global domain. Ensembles of five 30-day mean fields are obtained from integrations performed in the perpetual-January mode. The results are presented as anomalies, that is, departures of each ensemble mean from that produced in a control simulation with climatological SSTs. Large differences are found between the anomalies obtained using COADS and HIRS SSTs, even in the Northern Hemisphere where the datasets are most similar to each other. The internal variability of the circulation in the control simulation and the simulated atmospheric response to anomalous forcings appear to be linked in that the pattern of geopotential height anomalies obtained using COADS SSTs resembles the first empirical orthogonal function (EOF 1) in the control simulation. The corresponding pattern obtained using HIRS SSTs is substantially different and somewhat resembles EOF 2 in the sector from central North America to central Asia. To gain insight into the reasons for these results, three additional simulations are carried out with SST anomalies confined to regions where COADS SSTs are substantially warmer than HIRS SSTs. The regions correspond to warm pools in the northwest and northeast Pacific, and the northwest Atlantic. These warm pools tend to produce positive geopotential height anomalies in the northeastern part of the corresponding oceans. Both warm pools in the Pacific produce large

  1. A simulation study of control and display requirements for zero-experience general aviation pilots

    NASA Technical Reports Server (NTRS)

    Stewart, Eric C.

    1993-01-01

    The purpose of this simulation study was to define the basic human factor requirements for operating an airplane in all weather conditions. The basic human factors requirements are defined as those for an operator who is a complete novice for airplane operations but who is assumed to have automobile driving experience. These operators thus have had no piloting experience or training of any kind. The human factor requirements are developed for a practical task which includes all of the basic maneuvers required to go from one airport to another airport in limited visibility conditions. The task was quite demanding including following a precise path with climbing and descending turns while simultaneously changing airspeed. The ultimate goal of this research is to increase the utility of general aviation airplanes - that is, to make them a practical mode of transportation for a much larger segment of the general population. This can be accomplished by reducing the training and proficiency requirements of pilots while improving the level of safety. It is believed that advanced technologies such as fly-by-wire (or light), and head-up pictorial displays can be of much greater benefit to the general aviation pilot than to the full-time, professional pilot.

  2. Balancing an accurate representation of the molecular surface in generalized born formalisms with integrator stability in molecular dynamics simulations.

    PubMed

    Chocholousová, Jana; Feig, Michael

    2006-04-30

    Different integrator time steps in NVT and NVE simulations of protein and nucleic acid systems are tested with the GBMV (Generalized Born using Molecular Volume) and GBSW (Generalized Born with simple SWitching) methods. The simulation stability and energy conservation is investigated in relation to the agreement with the Poisson theory. It is found that very close agreement between generalized Born methods and the Poisson theory based on the commonly used sharp molecular surface definition results in energy drift and simulation artifacts in molecular dynamics simulation protocols with standard 2-fs time steps. New parameters are proposed for the GBMV method, which maintains very good agreement with the Poisson theory while providing energy conservation and stable simulations at time steps of 1 to 1.5 fs. PMID:16518883

  3. Relations between winter precipitation and atmospheric circulation simulated by the Geophysical Fluid Dynamics Laboratory general circulation model

    USGS Publications Warehouse

    McCabe, G.J., Jr.; Dettinger, M.D.

    1995-01-01

    General circulation model (GCM) simulations of atmospheric circulation are more reliable than GCM simulations of temperature and precipitation. In this study, temporal correlations between 700 hPa height anomalies simulated winter precipitation at eight locations in the conterminous United States are compared with corresponding correlations in observations. The objectives are to 1) characterize the relations between atmospheric circulation and winter precipitation simulated by the GFDL, GCM for selected locations in the conterminous USA, ii) determine whether these relations are similar to those found in observations of the actual climate system, and iii) determine if GFDL-simulated precipitation is forced by the same circulation patterns as in the real atmosphere. -from Authors

  4. Simulation of Venus polar vortices with the non-hydrostatic general circulation model

    NASA Astrophysics Data System (ADS)

    Rodin, Alexander V.; Mingalev, Oleg; Orlov, Konstantin

    2012-07-01

    The dynamics of Venus atmosphere in the polar regions presents a challenge for general circulation models. Numerous images and hyperspectral data from Venus Express mission shows that above 60 degrees latitude atmospheric motion is substantially different from that of the tropical and extratropical atmosphere. In particular, extended polar hoods composed presumably of fine haze particles, as well as polar vortices revealing mesoscale wave perturbations with variable zonal wavenumbers, imply the significance of vertical motion in these circulation elements. On these scales, however, hydrostatic balance commonly used in the general circulation models is no longer valid, and vertical forces have to be taken into account to obtain correct wind field. We present the first non-hydrostatic general circulation model of the Venus atmosphere based on the full set of gas dynamics equations. The model uses uniform grid with the resolution of 1.2 degrees in horizontal and 200 m in the vertical direction. Thermal forcing is simulated by means of relaxation approximation with specified thermal profile and time scale. The model takes advantage of hybrid calculations on graphical processors using CUDA technology in order to increase performance. Simulations show that vorticity is concentrated at high latitudes within planetary scale, off-axis vortices, precessing with a period of 30 to 40 days. The scale and position of these vortices coincides with polar hoods observed in the UV images. The regions characterized with high vorticity are surrounded by series of small vortices which may be caused by shear instability of the zonal flow. Vertical velocity component implies that in the central part of high vorticity areas atmospheric flow is downwelling and perturbed by mesoscale waves with zonal wavenumbers 1-4, resembling observed wave structures in the polar vortices. Simulations also show the existence of areas with strong vertical flow, concentrated in spiral branches extending

  5. Venus atmosphere simulated by a high-resolution general circulation model

    NASA Astrophysics Data System (ADS)

    Sugimoto, Norihiko

    2016-07-01

    An atmospheric general circulation model (AGCM) for Venus on the basis of AFES (AGCM For the Earth Simulator) have been developed (e.g., Sugimoto et al., 2014a) and a very high-resolution simulation is performed. The highest resolution of the model is T319L120; 960 times 480 horizontal grids (grid intervals are about 40 km) with 120 vertical layers (layer intervals are about 1 km). In the model, the atmosphere is dry and forced by the solar heating with the diurnal and semi-diurnal components. The infrared radiative process is simplified by adopting Newtonian cooling approximation. The temperature is relaxed to a prescribed horizontally uniform temperature distribution, in which a layer with almost neutral static stability observed in the Venus atmosphere presents. A fast zonal wind in a solid-body rotation is given as the initial state. Starting from this idealized superrotation, the model atmosphere reaches a quasi-equilibrium state within 1 Earth year and this state is stably maintained for more than 10 Earth years. The zonal-mean zonal flow with weak midlatitude jets has almost constant velocity of 120 m/s in latitudes between 45°S and 45°N at the cloud top levels, which agrees very well with observations. In the cloud layer, baroclinic waves develop continuously at midlatitudes and generate Rossby-type waves at the cloud top (Sugimoto et al., 2014b). At the polar region, warm polar vortex zonally surrounded by a cold latitude band (cold collar) is well reproduced (Ando et al., 2016). As for horizontal kinetic energy spectra, divergent component is broadly (k>10) larger than rotational component compared with that on Earth (Kashimura et al., in preparation). Finally, recent results for thermal tides and small-scale waves will be shown in the presentation. Sugimoto, N. et al. (2014a), Baroclinic modes in the Venus atmosphere simulated by GCM, Journal of Geophysical Research: Planets, Vol. 119, p1950-1968. Sugimoto, N. et al. (2014b), Waves in a Venus general

  6. Development and Implementation of Non-Newtonian Rheology Into the Generalized Fluid System Simulation Program (GFSSP)

    NASA Technical Reports Server (NTRS)

    DiSalvo, Roberto; Deaconu, Stelu; Majumdar, Alok

    2006-01-01

    One of the goals of this program was to develop the experimental and analytical/computational tools required to predict the flow of non-Newtonian fluids through the various system components of a propulsion system: pipes, valves, pumps etc. To achieve this goal we selected to augment the capabilities of NASA's Generalized Fluid System Simulation Program (GFSSP) software. GFSSP is a general-purpose computer program designed to calculate steady state and transient pressure and flow distributions in a complex fluid network. While the current version of the GFSSP code is able to handle various systems components the implicit assumption in the code is that the fluids in the system are Newtonian. To extend the capability of the code to non-Newtonian fluids, such as silica gelled fuels and oxidizers, modifications to the momentum equations of the code have been performed. We have successfully implemented in GFSSP flow equations for fluids with power law behavior. The implementation of the power law fluid behavior into the GFSSP code depends on knowledge of the two fluid coefficients, n and K. The determination of these parameters for the silica gels used in this program was performed experimentally. The n and K parameters for silica water gels were determined experimentally at CFDRC's Special Projects Laboratory, with a constant shear rate capillary viscometer. Batches of 8:1 (by weight) water-silica gel were mixed using CFDRC s 10-gallon gelled propellant mixer. Prior to testing the gel was allowed to rest in the rheometer tank for at least twelve hours to ensure that the delicate structure of the gel had sufficient time to reform. During the tests silica gel was pressure fed and discharged through stainless steel pipes ranging from 1", to 36", in length and three diameters; 0.0237", 0.032", and 0.047". The data collected in these tests included pressure at tube entrance and volumetric flowrate. From these data the uncorrected shear rate, shear stress, residence time

  7. Field Evaluation of the Generalized Maintenance Trainer-Simulator: I. Fleet Communications System. ; Technical Report No. 89.

    ERIC Educational Resources Information Center

    Rigney, J. W.; And Others

    This report describes the Generalized Maintenance Trainer-Simulator (GMTS), an instructional system designed to give electronics students intensive troubleshooting practice in a simulated hands-on training environment, and reports on a field evaluation of the GMTS applied to systems level troubleshooting in radio communications. The GMTS can be…

  8. A generalized model for simulating adsorption on porous media and checking for reversibility by desorption

    NASA Astrophysics Data System (ADS)

    Batzias, Fragiskos; Bountri, Athanasia; Sidiras, Dimitris

    2012-12-01

    Most adsorption kinetic models are of integer order (mainly of first and to a lesser extent of second order) with two parameters (rate constant and equilibrium parameter) and without an intercept, when used in their analytic form. In this work, we derive a four-parameter nth-order (n being not an integer, in general) model, simulating adsorption on porous media. We proved that this model implied best fitting to experimental data of dye adsorption on fir sawdust. Subsequently, a criterion of competitiveness is presented to find out which simplified form of a pre-set order is the second best, in order to obtain parameter values comparable to results already stored in corresponding Data Bases. Partial reversibility was also confirmed by desorption, from saturated-with-dye biomass to aquatic solution, using a Friendlichtype desorption isotherm.

  9. Verification of a generalized Aboav-Weaire law via experiment and large-scale simulation

    NASA Astrophysics Data System (ADS)

    Wang, H.; Liu, G. Q.; Chen, Y.; Li, W. W.

    2014-03-01

    Topological correlations in grain boundary networks are investigated on the basis of more than 14000 experimental grains and 9000 Monte Carlo-Potts model simulation grains. A generalized Aboav-Weaire law which serves as a description of the short- and long-range nearest-neighbor topological correlations, is proved to hold both in 2D grain structures and 2D cross-section structures. However, the nearest-neighbor topological correlations have no obvious influence on the rate of 2D grain growth, which is explicitly different from the case in three dimensions that was previously reported in Wang H., Liu G. Q., Song X. Y. and Luan J. H., EPL, 96 (2011) 38003.

  10. Real time simulation of nonlinear generalized predictive control for wind energy conversion system with nonlinear observer.

    PubMed

    Ouari, Kamel; Rekioua, Toufik; Ouhrouche, Mohand

    2014-01-01

    In order to make a wind power generation truly cost-effective and reliable, an advanced control techniques must be used. In this paper, we develop a new control strategy, using nonlinear generalized predictive control (NGPC) approach, for DFIG-based wind turbine. The proposed control law is based on two points: NGPC-based torque-current control loop generating the rotor reference voltage and NGPC-based speed control loop that provides the torque reference. In order to enhance the robustness of the controller, a disturbance observer is designed to estimate the aerodynamic torque which is considered as an unknown perturbation. Finally, a real-time simulation is carried out to illustrate the performance of the proposed controller. PMID:24021543

  11. A fully general relativistic numerical simulation code for spherically symmetric matter

    NASA Astrophysics Data System (ADS)

    Park, Dong-Ho; Cho, Inyong; Kang, Gungwon; Lee, Hyung Mok

    2013-02-01

    We present a fully general relativistic open-source code that can be used for simulating a system of spherically symmetric perfect fluid matter. It is based on the Arnowitt-Deser-Misner 3+1 formalism with maximal slicing and isotropic spatial coordinates. For hydrodynamic matter High Resolution Shock Capturing (HRSC) schemes with a monotonized central-difference limiter and approximated Riemann solvers are used in the Eulerian viewpoint. The accuracy and the convergence of our numerical code are verified by performing several test problems. These include a relativistic blast wave, relativistic spherical accretion of matter into a black hole, Tolman-Oppenheimer-Volkoff (TOV) stars and Oppenheimer-Snyder (OS) dust collapses. In particular, a dynamical code test is done for the OS collapse by explicitly performing numerical coordinate transformations between our coordinate 8system and the one used for the analytic solution. Finally, some TOV star solutions are presented for the Eddington-inspired Born-Infeld gravity theory.

  12. Strong scaling of general-purpose molecular dynamics simulations on GPUs

    NASA Astrophysics Data System (ADS)

    Glaser, Jens; Nguyen, Trung Dac; Anderson, Joshua A.; Lui, Pak; Spiga, Filippo; Millan, Jaime A.; Morse, David C.; Glotzer, Sharon C.

    2015-07-01

    We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, 2013). Our approach is inspired by a traditional CPU-based code, LAMMPS (Plimpton, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5 ×.

  13. Finite-difference simulation and visualization of elastodynamics in time-evolving generalized curvilinear coordinates

    NASA Technical Reports Server (NTRS)

    Kaul, Upender K. (Inventor)

    2009-01-01

    Modeling and simulation of free and forced structural vibrations is essential to an overall structural health monitoring capability. In the various embodiments, a first principles finite-difference approach is adopted in modeling a structural subsystem such as a mechanical gear by solving elastodynamic equations in generalized curvilinear coordinates. Such a capability to generate a dynamic structural response is widely applicable in a variety of structural health monitoring systems. This capability (1) will lead to an understanding of the dynamic behavior of a structural system and hence its improved design, (2) will generate a sufficiently large space of normal and damage solutions that can be used by machine learning algorithms to detect anomalous system behavior and achieve a system design optimization and (3) will lead to an optimal sensor placement strategy, based on the identification of local stress maxima all over the domain.

  14. Simulation of Indian Monsoon Variability in the Medieval Warm Period using ECHAM5 General Circulation Model

    NASA Astrophysics Data System (ADS)

    Polanski, Stefan; Fallah, Bijan; Prasad, Sushma; Cubasch, Ulrich

    2013-04-01

    Within the framework of the DFG research group HIMPAC, the general circulation model ECHAM5 has been used to simulate the Indian monsoon and its variability during the Medieval Warm Period (MWP; 900-1100 AD) and for recent climate (REC; 1800-2000 AD). The focus is on the analysis of internal and external drivers leading to extreme rainfall events over India from interannual to multidecadal time scale. An evaluation of spatio-temporal monsoon patterns with present-day observation data is in agreement with other state-of-the-art monsoon modeling studies. The simulated monsoon intensity on multidecadal time scale is weakened (enhanced) in summer (winter) due to colder (warmer) SSTs in the Indian Ocean. Variations in solar insolation are the main drivers for these SST anomalies, verified by very high temporal correlations between Total Solar Irradiance and All-India-Monsoon-Rainfall in summer monsoon months (-0.95). The external solar forcing is coupled and overlain by internal climate modes of the Ocean (ENSO and IOD) with asynchronous intensities and lengths of periods. In addition, the model simulations have been compared with a relative moisture index derived from paleoclimatic reconstructions based on various proxies and archives in India (Anoop et al., 2012 (under revision); Bhattacharya et al., 2007; Chauhan et al., 2000; Denniston et al., 2000; Ely et al., 1999; Kar et al., 2002; Ponton et al., 2012; Prasad et al., 2012 (under revision)). In this context, the reconstructed climate of the well-dated Lonar record in Central India has been highlighted and evaluated the first time (Anoop et al., 2012 (under revision); Prasad et al., 2012 (under revision)). Particularly with regard to the long continuously chronology of the last 11000 years, the Lonar site gives a unique possibility for a comparison of long-term climate time series. The simulated relative annual rainfall anomalies ("MWP" minus "REC") are in agreement with the reconstructed moisture index. The dry

  15. Large-eddy simulation of airflow and heat transfer in a general ward of hospital

    NASA Astrophysics Data System (ADS)

    Hasan, Md. Farhad; Himika, Taasnim Ahmed; Molla, Md. Mamun

    2016-07-01

    In this paper, a very popular alternative computational technique, the Lattice Boltzmann Method (LBM) has been used for Large-Eddy Simulation (LES) of airflow and heat transfer in general ward of hospital. Different Reynolds numbers have been used to study the airflow pattern. In LES, Smagorinsky turbulence model has been considered and a discussion has been conducted in brief. A code validation has been performed comparing the present results with benchmark results for lid-driven cavity problem and the results are found to agree very well. LBM is demonstrated through simulation in forced convection inside hospital ward with six beds with a partition in the middle, which acted like a wall. Changes in average rate of heat transfer in terms of average Nusselt numbers have also been recorded in tabular format and necessary comparison has been showed. It was found that partition narrowed the path for airflow and once the air overcame this barrier, it got free space and turbulence appeared. For higher turbulence, the average rate of heat transfer increased and patients near the turbulence zone released maximum heat and felt more comfortable.

  16. General Force-Field Parametrization Scheme for Molecular Dynamics Simulations of Conjugated Materials in Solution.

    PubMed

    Wildman, Jack; Repiščák, Peter; Paterson, Martin J; Galbraith, Ian

    2016-08-01

    We describe a general scheme to obtain force-field parameters for classical molecular dynamics simulations of conjugated polymers. We identify a computationally inexpensive methodology for calculation of accurate intermonomer dihedral potentials and partial charges. Our findings indicate that the use of a two-step methodology of geometry optimization and single-point energy calculations using DFT methods produces potentials which compare favorably to high level theory calculation. We also report the effects of varying the conjugated backbone length and alkyl side-chain lengths on the dihedral profiles and partial charge distributions and determine the existence of converged lengths above which convergence is achieved in the force-field parameter sets. We thus determine which calculations are required for accurate parametrization and the scope of a given parameter set for variations to a given molecule. We perform simulations of long oligomers of dioctylfluorene and hexylthiophene in explicit solvent and find peristence lengths and end-length distributions consistent with experimental values. PMID:27397762

  17. Application of the general thermal field model to simulate the behaviour of nanoscale Cu field emitters

    SciTech Connect

    Eimre, Kristjan; Aabloo, Alvo; Parviainen, Stefan Djurabekova, Flyura; Zadin, Vahur

    2015-07-21

    Strong field electron emission from a nanoscale tip can cause a temperature rise at the tip apex due to Joule heating. This becomes particularly important when the current value grows rapidly, as in the pre-breakdown (the electrostatic discharge) condition, which may occur near metal surfaces operating under high electric fields. The high temperatures introduce uncertainties in calculations of the current values when using the Fowler–Nordheim equation, since the thermionic component in such conditions cannot be neglected. In this paper, we analyze the field electron emission currents as the function of the applied electric field, given by both the conventional Fowler–Nordheim field emission and the recently developed generalized thermal field emission formalisms. We also compare the results in two limits: discrete (atomistic simulations) and continuum (finite element calculations). The discrepancies of both implementations and their effect on final results are discussed. In both approaches, the electric field, electron emission currents, and Joule heating processes are simulated concurrently and self-consistently. We show that the conventional Fowler–Nordheim equation results in significant underestimation of electron emission currents. We also show that Fowler–Nordheim plots used to estimate the field enhancement factor may lead to significant overestimation of this parameter especially in the range of relatively low electric fields.

  18. Evaluating Parameterizations in General Circulation Models: Climate Simulation Meets Weather Prediction

    SciTech Connect

    Phillips, T J; Potter, G L; Williamson, D L; Cederwall, R T; Boyle, J S; Fiorino, M; Hnilo, J J; Olson, J G; Xie, S; Yio, J J

    2004-05-06

    To significantly improve the simulation of climate by general circulation models (GCMs), systematic errors in representations of relevant processes must first be identified, and then reduced. This endeavor demands that the GCM parameterizations of unresolved processes, in particular, should be tested over a wide range of time scales, not just in climate simulations. Thus, a numerical weather prediction (NWP) methodology for evaluating model parameterizations and gaining insights into their behavior may prove useful, provided that suitable adaptations are made for implementation in climate GCMs. This method entails the generation of short-range weather forecasts by a realistically initialized climate GCM, and the application of six-hourly NWP analyses and observations of parameterized variables to evaluate these forecasts. The behavior of the parameterizations in such a weather-forecasting framework can provide insights on how these schemes might be improved, and modified parameterizations then can be tested in the same framework. In order to further this method for evaluating and analyzing parameterizations in climate GCMs, the U.S. Department of Energy is funding a joint venture of its Climate Change Prediction Program (CCPP) and Atmospheric Radiation Measurement (ARM) Program: the CCPP-ARM Parameterization Testbed (CAPT). This article elaborates the scientific rationale for CAPT, discusses technical aspects of its methodology, and presents examples of its implementation in a representative climate GCM.

  19. Monte Carlo Simulations of PAC-Spectra as a General Approach to Dynamic Interactions

    NASA Astrophysics Data System (ADS)

    Danielsen, Eva; Jørgensen, Lars Elkjær; Sestoft, Peter

    Time Dependent Perturbed Angular Correlations of γ-rays (PAC) can be used to study hyperfine interactions of a dynamic nature. However, the exact effect of the dynamic interaction on the PAC-spectrum is sometimes difficult to derive analytically. A new approach based on Monte Carlo simulations is therefore suggested, here implemented as a Fortran 90 program for simulating PAC spectra of dynamic electric field gradients of any origin. The program is designed for the most common experimental condition where the intermediate level has spin 5/2, but the approach can equally well be used for other spin states. Codes for 4 different situations have been developed: (1) Rotational diffusion by jumps; used as a test case. (2) Jumps between two states with different electric field gradients, different lifetimes and different orientations of the electric field gradient principal axes. (3) Relaxation of one state to another. (4) Molecules adhering to a surface with random rotational jumps around the axis perpendicular to the surface. To illustrate how this approach can be used to improve data-interpretation, previously published data on 111mCd-plastocyanin and 111Ag-plastocyanin are re-considered. The strength of this novel approach is its simplicity and generality so that other dynamic processes can easily be included by only adding new program units describing the random process behind the dynamics. The program is hereby made publicly available.

  20. A general method for spatially coarse-graining Metropolis Monte Carlo simulations onto a lattice

    NASA Astrophysics Data System (ADS)

    Liu, Xiao; Seider, Warren D.; Sinno, Talid

    2013-03-01

    A recently introduced method for coarse-graining standard continuous Metropolis Monte Carlo simulations of atomic or molecular fluids onto a rigid lattice of variable scale [X. Liu, W. D. Seider, and T. Sinno, Phys. Rev. E 86, 026708 (2012)], 10.1103/PhysRevE.86.026708 is further analyzed and extended. The coarse-grained Metropolis Monte Carlo technique is demonstrated to be highly consistent with the underlying full-resolution problem using a series of detailed comparisons, including vapor-liquid equilibrium phase envelopes and spatial density distributions for the Lennard-Jones argon and simple point charge water models. In addition, the principal computational bottleneck associated with computing a coarse-grained interaction function for evolving particle positions on the discretized domain is addressed by the introduction of new closure approximations. In particular, it is shown that the coarse-grained potential, which is generally a function of temperature and coarse-graining level, can be computed at multiple temperatures and scales using a single set of free energy calculations. The computational performance of the method relative to standard Monte Carlo simulation is also discussed.

  1. General continuum boundary conditions for miscible binary fluids from molecular dynamics simulations.

    PubMed

    Denniston, Colin; Robbins, Mark O

    2006-12-01

    Molecular dynamics simulations are used to explore the flow behavior and diffusion of miscible fluids near solid surfaces. The solid produces deviations from bulk fluid behavior that decay over a distance of the order of the fluid correlation length. Atomistic results are mapped onto two types of continuum model: Mesoscopic models that follow this decay and conventional sharp interface boundary conditions for the stress and velocity. The atomistic results, and mesoscopic models derived from them, are consistent with the conventional Marangoni stress boundary condition. However, there are deviations from the conventional Navier boundary condition that states that the slip velocity between wall and fluid is proportional to the strain rate. A general slip boundary condition is derived from the mesoscopic model that contains additional terms associated with the Marangoni stress and diffusion, and is shown to describe the atomistic simulations. The additional terms lead to strong flows when there is a concentration gradient. The potential for using this effect to make a nanomotor or pump is evaluated. PMID:17166010

  2. General Force-Field Parametrization Scheme for Molecular Dynamics Simulations of Conjugated Materials in Solution

    PubMed Central

    2016-01-01

    We describe a general scheme to obtain force-field parameters for classical molecular dynamics simulations of conjugated polymers. We identify a computationally inexpensive methodology for calculation of accurate intermonomer dihedral potentials and partial charges. Our findings indicate that the use of a two-step methodology of geometry optimization and single-point energy calculations using DFT methods produces potentials which compare favorably to high level theory calculation. We also report the effects of varying the conjugated backbone length and alkyl side-chain lengths on the dihedral profiles and partial charge distributions and determine the existence of converged lengths above which convergence is achieved in the force-field parameter sets. We thus determine which calculations are required for accurate parametrization and the scope of a given parameter set for variations to a given molecule. We perform simulations of long oligomers of dioctylfluorene and hexylthiophene in explicit solvent and find peristence lengths and end-length distributions consistent with experimental values. PMID:27397762

  3. Longitudinal biases in the Seychelles Dome simulated by 35 ocean-atmosphere coupled general circulation models

    NASA Astrophysics Data System (ADS)

    Nagura, Motoki; Sasaki, Wataru; Tozuka, Tomoki; Luo, Jing-Jia; Behera, Swadhin K.; Yamagata, Toshio

    2013-02-01

    Seychelles Dome refers to the shallow climatological thermocline in the southwestern Indian Ocean, where ocean wave dynamics efficiently affect sea surface temperature, allowing sea surface temperature anomalies to be predicted up to 1-2 years in advance. Accurate reproduction of the dome by ocean-atmosphere coupled general circulation models (CGCMs) is essential for successful seasonal predictions in the Indian Ocean. This study examines the Seychelles Dome as simulated by 35 CGCMs, including models used in phase five of the Coupled Model Intercomparison Project (CMIP5). Among the 35 CGCMs, 14 models erroneously produce an upwelling dome in the eastern half of the basin whereas the observed Seychelles Dome is located in the southwestern tropical Indian Ocean. The annual mean Ekman pumping velocity in these models is found to be almost zero in the southern off-equatorial region. This result is inconsistent with observations, in which Ekman upwelling acts as the main cause of the Seychelles Dome. In the models reproducing an eastward-displaced dome, easterly biases are prominent along the equator in boreal summer and fall, which result in shallow thermocline biases along the Java and Sumatra coasts via Kelvin wave dynamics and a spurious upwelling dome in the region. Compared to the CMIP3 models, the CMIP5 models are even worse in simulating the dome longitudes.

  4. General relativistic magnetohydrodynamic simulations of binary neutron star mergers with the APR4 equation of state

    NASA Astrophysics Data System (ADS)

    Endrizzi, A.; Ciolfi, R.; Giacomazzo, B.; Kastaun, W.; Kawamura, T.

    2016-08-01

    We present new results of fully general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers performed with the Whisky code. All the models use a piecewise polytropic approximation of the APR4 equation of state for cold matter, together with a ‘hybrid’ part to incorporate thermal effects during the evolution. We consider both equal and unequal-mass models, with total masses such that either a supramassive NS or a black hole is formed after merger. Each model is evolved with and without a magnetic field initially confined to the stellar interior. We present the different gravitational wave (GW) signals as well as a detailed description of the matter dynamics (magnetic field evolution, ejected mass, post-merger remnant/disk properties). Our simulations provide new insights into BNS mergers, the associated GW emission and the possible connection with the engine of short gamma-ray bursts (both in the ‘standard’ and in the ‘time-reversal’ scenarios) and other electromagnetic counterparts.

  5. TOUGH2: A general-purpose numerical simulator for multiphase nonisothermal flows

    SciTech Connect

    Pruess, K.

    1991-06-01

    Numerical simulators for multiphase fluid and heat flows in permeable media have been under development at Lawrence Berkeley Laboratory for more than 10 yr. Real geofluids contain noncondensible gases and dissolved solids in addition to water, and the desire to model such `compositional` systems led to the development of a flexible multicomponent, multiphase simulation architecture known as MULKOM. The design of MULKOM was based on the recognition that the mass-and energy-balance equations for multiphase fluid and heat flows in multicomponent systems have the same mathematical form, regardless of the number and nature of fluid components and phases present. Application of MULKOM to different fluid mixtures, such as water and air, or water, oil, and gas, is possible by means of appropriate `equation-of-state` (EOS) modules, which provide all thermophysical and transport parameters of the fluid mixture and the permeable medium as a function of a suitable set of primary thermodynamic variables. Investigations of thermal and hydrologic effects from emplacement of heat-generating nuclear wastes into partially water-saturated formations prompted the development and release of a specialized version of MULKOM for nonisothermal flow of water and air, named TOUGH. TOUGH is an acronym for `transport of unsaturated groundwater and heat` and is also an allusion to the tuff formations at Yucca Mountain, Nevada. The TOUGH2 code is intended to supersede TOUGH. It offers all the capabilities of TOUGH and includes a considerably more general subset of MULKOM modules with added capabilities. The paper briefly describes the simulation methodology and user features.

  6. Internal versus SST-forced atmospheric variability as simulated by an atmospheric general circulation model

    SciTech Connect

    Harzallah, A.; Sadourny, R.

    1995-03-01

    The variability of atmospheric flow is analyzed by separating it into an internal part due to atmospheric dynamics only and an external (or forced) part due to the variability of sea surface temperature forcing. The two modes of variability are identified by performing an ensemble of seven independent long-term simulations of the atmospheric response to observed SST (1970-1988) with the LMD atmospheric general circulation model. The forced variability is defined from the analysis of the ensemble mean and the internal variability from the analysis of deviations from the ensemble mean. Emphasis is put on interannual variability of sea level pressure and 500-hPa geopotential height for the Northern Hemisphere winter. In view of the large systematic errors related to the relatively small number of realizations, unbiased variance estimators have been developed. Although statistical significance is not reached in some extratropical regions, large significant extratropical responses are found at the North Pacific-Alaska sector for SLP and over western Canada and the Aleutians for 500-hPa geopotential height. The influence of SST variations on internal variability is also examined by using a 7-year simulation using the climatological SST seasonal cycle. It is found that interannual SST changes strongly influence the geographical distribution of internal variability; in particular, it tends to increase it over oceans. EOF decompositions, showing that the model realistically simulates the leading observed variability modes. The geographical structure of internal variability patterns is found to be similar to that of total variability, although similar modes tend to evolve rather differently in time. The zonally symmetric seesaw dominates the internal variability for both observed and climatologically prescribed SST. 46 refs., 15 figs., 3 tabs.

  7. An Experimental Analysis of General Case Simulation Instruction and the Establishment and Maintenance of Work Performance in Severely Handicapped Students.

    ERIC Educational Resources Information Center

    Woolcock, William Woodrow

    This doctoral dissertation examines the extent to which general case simulation instruction on a janitorial task sequence and a housekeeping task sequence conducted with four secondary and postsecondary age persons with moderate mental retardation resulted in generalized performance. A multiple baseline design across subjects and behaviors was…

  8. A Multi-mission Event-Driven Component-Based System for Support of Flight Software Development, ATLO, and Operations first used by the Mars Science Laboratory (MSL) Project

    NASA Technical Reports Server (NTRS)

    Dehghani, Navid; Tankenson, Michael

    2006-01-01

    This viewgraph presentation reviews the architectural description of the Mission Data Processing and Control System (MPCS). MPCS is an event-driven, multi-mission ground data processing components providing uplink, downlink, and data management capabilities which will support the Mars Science Laboratory (MSL) project as its first target mission. MPCS is designed with these factors (1) Enabling plug and play architecture (2) MPCS has strong inheritance from GDS components that have been developed for other Flight Projects (MER, MRO, DAWN, MSAP), and are currently being used in operations and ATLO, and (3) MPCS components are Java-based, platform independent, and are designed to consume and produce XML-formatted data

  9. General Fluid System Simulation Program to Model Secondary Flows in Turbomachinery

    NASA Technical Reports Server (NTRS)

    Majumdar, Alok K.; Van Hoosier, Katherine P.

    1995-01-01

    The complexity and variety of turbomachinery flow circuits created a need for a general fluid system simulation program for test data anomaly resolution as well as design review. The objective of the paper is to present a computer program that has been developed to support Marshall Space Flight Center's turbomachinery internal flow analysis efforts. The computer program solves for the mass. energy and species conservation equation at each node and flow rate equation at each branch of the network by a novel numerical procedure which is a combination of both Newton-Ralphson and successive substitution method and uses a thermodynamic property program for computing real gas properties. A generalized, robust, modular, and 'user-friendly' computer program has been developed to model internal flow rates, pressures, temperatures, concentrations of gas mixtures and axial thrusts. The program can be used for any network for compressible and incompressible flows, choked flow, change of phase and gaseous mixturecs. The code has been validated by comparing the predictions with Space Shuttle Main Engine test data.

  10. Generalized SIMD algorithm for efficient EM-PIC simulations on modern CPUs

    NASA Astrophysics Data System (ADS)

    Fonseca, Ricardo; Decyk, Viktor; Mori, Warren; Silva, Luis

    2012-10-01

    There are several relevant plasma physics scenarios where highly nonlinear and kinetic processes dominate. Further understanding of these scenarios is generally explored through relativistic particle-in-cell codes such as OSIRIS [1], but this algorithm is computationally intensive, and efficient use high end parallel HPC systems, exploring all levels of parallelism available, is required. In particular, most modern CPUs include a single-instruction-multiple-data (SIMD) vector unit that can significantly speed up the calculations. In this work we present a generalized PIC-SIMD algorithm that is shown to work efficiently with different CPU (AMD, Intel, IBM) and vector unit types (2-8 way, single/double). Details on the algorithm will be given, including the vectorization strategy and memory access. We will also present performance results for the various hardware variants analyzed, focusing on floating point efficiency. Finally, we will discuss the applicability of this type of algorithm for EM-PIC simulations on GPGPU architectures [2]. [4pt] [1] R. A. Fonseca et al., LNCS 2331, 342, (2002)[0pt] [2] V. K. Decyk, T. V. Singh; Comput. Phys. Commun. 182, 641-648 (2011)