Simulation Framework (simulation + framework)

Distribution by Scientific Domains


Selected Abstracts


Enabling interactive and collaborative oil reservoir simulations on the Grid

CONCURRENCY AND COMPUTATION: PRACTICE & EXPERIENCE, Issue 11 2005
Manish Parashar
Abstract Grid-enabled infrastructures and problem-solving environments can significantly increase the scale, cost-effectiveness and utility of scientific simulations, enabling highly accurate simulations that provide in-depth insight into complex phenomena. This paper presents a prototype of such an environment, i.e. an interactive and collaborative problem-solving environment for the formulation, development, deployment and management of oil reservoir and environmental flow simulations in computational Grid environments. The project builds on three independent research efforts: (1) the IPARS oil reservoir and environmental flow simulation framework; (2) the NetSolve Grid engine; and (3) the Discover Grid-based computational collaboratory. Its primary objective is to demonstrate the advantages of an integrated simulation infrastructure towards effectively supporting scientific investigation on the Grid, and to investigate the components and capabilities of such an infrastructure. Copyright © 2005 John Wiley & Sons, Ltd. [source]


Efficient sampling and data reduction techniques for probabilistic seismic lifeline risk assessment

EARTHQUAKE ENGINEERING AND STRUCTURAL DYNAMICS, Issue 10 2010
Nirmal Jayaram
Abstract Probabilistic seismic risk assessment for spatially distributed lifelines is less straightforward than for individual structures. While procedures such as the ,PEER framework' have been developed for risk assessment of individual structures, these are not easily applicable to distributed lifeline systems, due to difficulties in describing ground-motion intensity (e.g. spectral acceleration) over a region (in contrast to ground-motion intensity at a single site, which is easily quantified using Probabilistic Seismic Hazard Analysis), and since the link between the ground-motion intensities and lifeline performance is usually not available in closed form. As a result, Monte Carlo simulation (MCS) and its variants are well suited for characterizing ground motions and computing resulting losses to lifelines. This paper proposes a simulation-based framework for developing a small but stochastically representative catalog of earthquake ground-motion intensity maps that can be used for lifeline risk assessment. In this framework, Importance Sampling is used to preferentially sample ,important' ground-motion intensity maps, and K -Means Clustering is used to identify and combine redundant maps in order to obtain a small catalog. The effects of sampling and clustering are accounted for through a weighting on each remaining map, so that the resulting catalog is still a probabilistically correct representation. The feasibility of the proposed simulation framework is illustrated by using it to assess the seismic risk of a simplified model of the San Francisco Bay Area transportation network. A catalog of just 150 intensity maps is generated to represent hazard at 1038 sites from 10 regional fault segments causing earthquakes with magnitudes between five and eight. The risk estimates obtained using these maps are consistent with those obtained using conventional MCS utilizing many orders of magnitudes more ground-motion intensity maps. Therefore, the proposed technique can be used to drastically reduce the computational expense of a simulation-based risk assessment, without compromising the accuracy of the risk estimates. This will facilitate computationally intensive risk analysis of systems such as transportation networks. Finally, the study shows that the uncertainties in the ground-motion intensities and the spatial correlations between ground-motion intensities at various sites must be modeled in order to obtain unbiased estimates of lifeline risk. Copyright © 2010 John Wiley & Sons, Ltd. [source]


Coupled Navier,Stokes,Molecular dynamics simulations using a multi-physics flow simulation framework

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2010
R. Steijl
Abstract Simulation of nano-scale channel flows using a coupled Navier,Stokes/Molecular Dynamics (MD) method is presented. The flow cases serve as examples of the application of a multi-physics computational framework put forward in this work. The framework employs a set of (partially) overlapping sub-domains in which different levels of physical modelling are used to describe the flow. This way, numerical simulations based on the Navier,Stokes equations can be extended to flows in which the continuum and/or Newtonian flow assumptions break down in regions of the domain, by locally increasing the level of detail in the model. Then, the use of multiple levels of physical modelling can reduce the overall computational cost for a given level of fidelity. The present work describes the structure of a parallel computational framework for such simulations, including details of a Navier,Stokes/MD coupling, the convergence behaviour of coupled simulations as well as the parallel implementation. For the cases considered here, micro-scale MD problems are constructed to provide viscous stresses for the Navier,Stokes equations. The first problem is the planar Poiseuille flow, for which the viscous fluxes on each cell face in the finite-volume discretization are evaluated using MD. The second example deals with fully developed three-dimensional channel flow, with molecular level modelling of the shear stresses in a group of cells in the domain corners. An important aspect in using shear stresses evaluated with MD in Navier,Stokes simulations is the scatter in the data due to the sampling of a finite ensemble over a limited interval. In the coupled simulations, this prevents the convergence of the system in terms of the reduction of the norm of the residual vector of the finite-volume discretization of the macro-domain. Solutions to this problem are discussed in the present work, along with an analysis of the effect of number of realizations and sample duration. The averaging of the apparent viscosity for each cell face, i.e. the ratio of the shear stress predicted from MD and the imposed velocity gradient, over a number of macro-scale time steps is shown to be a simple but effective method to reach a good level of convergence of the coupled system. Finally, the parallel efficiency of the developed method is demonstrated. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Risk management for a global supply chain planning under uncertainty: Models and algorithms

AICHE JOURNAL, Issue 4 2009
Fengqi You
Abstract In this article, we consider the risk management for mid-term planning of a global multi-product chemical supply chain under demand and freight rate uncertainty. A two-stage stochastic linear programming approach is proposed within a multi-period planning model that takes into account the production and inventory levels, transportation modes, times of shipments, and customer service levels. To investigate the potential improvement by using stochastic programming, we describe a simulation framework that relies on a rolling horizon approach. The studies suggest that at least 5% savings in the total real cost can be achieved compared with the deterministic case. In addition, an algorithm based on the multi-cut L-shaped method is proposed to effectively solve the resulting large scale industrial size problems. We also introduce risk management models by incorporating risk measures into the stochastic programming model, and multi-objective optimization schemes are implemented to establish the tradeoffs between cost and risk. To demonstrate the effectiveness of the proposed stochastic models and decomposition algorithms, a case study of a realistic global chemical supply chain problem is presented. © 2009 American Institute of Chemical Engineers AIChE J, 2009 [source]


The role of pre- and post- alighting detection mechanisms in the responses to patch size by specialist herbivores

OIKOS, Issue 3 2005
Tibor Bukovinszky
Experimental data on the relationship between plant patch size and population density of herbivores within fields often deviates from predictions of the theory of island biogeography and the resource concentration hypothesis. Here we argue that basic features of foraging behaviour can explain different responses of specialist herbivores to habitat heterogeneity. In a combination of field and simulation studies, we applied basic knowledge on the foraging strategies of three specialist herbivores: the cabbage aphid (Brevicoryne brassicae), the cabbage butterfly (Pieris rapae L.) and the diamondback moth (Plutella xylostella L.), to explain differences in their responses to small scale fragmentation of their habitat. In our field study, populations of the three species responded to different sizes of host plant patches (9 plants and 100 plants) in different ways. Densities of winged cabbage aphids were independent of patch size. Egg-densities of the cabbage butterfly were higher in small than in large patches. Densities of diamondback moth adults were higher in large patches than in small patches. When patches in a background of barley were compared with those in grass, densities of the cabbage aphid and the diamondback moth were reduced, but not cabbage butterfly densities. To explore the role of foraging behaviour of herbivores on their response to patch size, a spatially explicit individual-based simulation framework was used. The sensory abilities of the insects to detect and respond to contact, olfactory or visual cues were varied. Species with a post-alighting host recognition behaviour (cabbage aphid) could only use contact cues from host plants encountered after landing. In contrast, species capable with a pre-alighting recognition behaviour, based on visual (cabbage butterfly) or olfactory (diamondback moth) cues, were able to recognise a preferred host plant whilst in flight. These three searching modalities were studied by varying the in flight detection abilities, the displacement speed and the arrestment response to host plants by individuals. Simulated patch size , density relationships were similar to those observed in the field. The importance of pre- and post- alighting detection in the responses of herbivores to spatial heterogeneity of the habitat is discussed. [source]