Jung, Hyun Kim

Parth, Nagarkar

Poccia, SILVESTRO ROBERTO

Sapino, Maria Luisa

Selcuk Candan, K.

Shengyu, Huang

Sicong, Liu

Xilun, Chen

Xinsheng, Li

Yash, Garg

English

Institutional Research Information System University of Turin

07 January 2023AperTO - Archivio Istituzionale Open Access dell'Università di TorinoOriginal Citation:SIMDMS: Data Management and Analysis to Support Decision Making through Large SimulationEnsemblesPublisher:Published version:DOI:10.5441/002/edbt.2017.75Terms of use:Open Access(Article begins on next page)Anyone can freely access the full text of works made available as "Open Access". Works made availableunder a Creative Commons license can be used according to the terms and conditions of said license. Useof all other works requires consent of the right holder (author or publisher) if not exempted from copyrightprotection by the applicable law.Availability:OpenProceedings.orgThis is the author's manuscriptThis version is available http://hdl.handle.net/2318/1647179 since 2021-04-29T18:32:50ZSIMDMS: Data Management and Analysis to SupportDecision Making through Large Simulation Ensembles ∗Silvestro Poccia Maria Luisa SapinoUniversity of Torino, Italy{poccia,mlsapino}@di.unito.itSicong Liu Xilun Chen Yash Garg Shengyu HuangJung Hyun Kim Xinsheng Li Parth Nagarkar K. Selçuk CandanArizona State University, USA{ s.liu,xilun.chen,ygarg,shuang54,jkim294,lxinshen,pnagarka,candan}@asu.eduABSTRACTData- and model-driven computer simulations are increas-ingly critical in many application domains. These simula-tions may track 100s or 1000s of inter-dependent parame-ters, spanning multiple layers and spatial-temporal frames,affected by complex dynamic processes operating at differ-ent resolutions. Because of the size and complexity of thedata and the varying spatial and temporal scales at whichthe key processes operate, experts often lack the means toanalyze results of large simulation ensembles, understandrelevant processes, and assess the robustness of conclusionsdriven from the resulting simulations. Moreover, data andmodels dynamically evolve over time requiring continuousadaptation of simulation ensembles. The simDMS platformaims to address the key challenges underlying the creationand use of large simulation ensembles and enables (a) execu-tion, storage, and indexing of large ensemble simulation datasets and the corresponding models; and (b) search, analysis,and exploration of ensemble simulation data sets to enableensemble-based decision support.KeywordsSimulation ensembles, multivariate time series1. INTRODUCTIONData- and model-driven computer simulations are increas-ingly critical in many application domains.∗Research is supported by NSF#1318788 “Data Managementfor Real-Time Data Driven Epidemic Spread Simulations”,NSF#1339835 “E-SDMS: Energy Simulation Data Management Sys-tem Software”, NSF#1610282 “DataStorm: A Data Enabled Systemfor End-to-End Disaster Planning and Response”, NSF#1633381“BIGDATA: Discovering Context-Sensitive Impact in ComplexSystems”, and “FourCmodeling”: EU-H2020 Marie Sklodowska-Curiegrant agreement No 690817.c©2017, Copyright is with the authors. Published in Proc. 20th Inter-national Conference on Extending Database Technology (EDBT), March21-24, 2017 - Venice, Italy: ISBN 978-3-89318-073-8, on OpenProceed-ings.org. Distribution of this paper is permitted under the terms of the Cre-ative Commons license CC-by-nc-nd 4.0Figure 1: Simulation ensembles are (a) multi-variate, (b) multi-modal (temporal, spatial, hi-erarchical, graphical), (c) multi-layer, (d) multi-resolution, and (e) inter-dependent (i.e., observa-tions of interest depend on and impact each other)Epidemic Simulation Ensembles: For example, for pre-dicting geo-temporal evolution of epidemics and assessingthe impact of interventions, experts often rely on epidemicspread simulation software such as (e.g., GLEaM [2] andSTEM [3]). The GLEaM simulation engine, for example,consists of three layers: (a) a population layer, (b) a mo-bility layer which includes both long-range air travel andshort-range commuting patterns between adjacent subpop-ulations, and (c) an epidemic layer which allows the user tospecify parameters (such as reproductive number and sea-sonality) for the infectious disease, initial outbreak condi-tions (e.g. seeding of the epidemic and the immunity profileof the subpopulation), and intervention measures.Building Energy Simulation Ensembles: Similarly, ef-fective building energy management, leading to more sus-tainable building systems and architectural designs withmonitoring, prioritization, and adaptation of building com-ponents and subsystems, requires large data-driven simula-tions involving (a) location and climate information for thecity in which the building is located, (b) building construc-tion information, such as building geometry and surface con-structions (including exterior walls, interior walls, partitions,floors, ceilings, roofs, windows and doors), (c) building useinformation, including the lighting and other equipment (e.g.electric, gas, etc.) and the number of people in each area ofthe building, (d) building thermostatic control information,including the temperature control strategy for each area, (e)heating, ventilation, and air conditioning (HVAC) operationDemonstration  Series ISSN: 2367-2005 582 10.5441/002/edbt.2017.75Execution Support Analysis Support epiViz External Simulation Engines (STEM) epiStore      Parameter and Model Analysis Feature Extraction and Indexing epiRun Query and Search Clustering and Classification Figure 2: simDMS system overview (instantiated withepidemic simulation ensembles)and scheduling information, and (f) central plant informa-tion for specification and scheduling of boilers, chillers, andother equipment. EnergyPlus software, for example, relieson the description of the building’s physical make-up andassociated mechanical and other systems and includes time-step based simulation for many energy-related building pa-rameters [1].1.1 Challenge: Ensemble based DecisionsWhile, in most cases, very powerful simulation softwareexist, using these simulation software for decision makingfaces several significant challenges: (a) Creating correct sim-ulation models is a costly operation, and it is often the casethat the designed simulation models are incomplete or im-precise. (b) Also, the execution of a simulation can be verycostly, given the fact that complex, inter-dependent param-eters affected by complex dynamic processes at varying spa-tial and temporal scales have to be taken into account. (c) Athird major source of cost is the simulation ensemble anal-ysis: because of the size and complexity of the data andthe varying spatial and temporal scales at which the keyprocesses operate, experts often lack the means of analyzingresults of large simulation ensembles, understanding relevantprocesses, and assessing the robustness of conclusions drivenfrom the resulting simulations.As visualized in Figure 1, the key characteristics of thesimulation data sets include the following: (a) multi-variate,(b) multi-modal (temporal, spatial, hierarchical, graphical),(c) multi-layer, (d) multi-resolution, and (e) inter-dependent(i.e., observations of interest depend on and impact eachother). In particular, simulations may track 100s or 1000sof inter-dependent parameters, spanning multiple layers andspatial-temporal frames, affected by complex dynamic pro-cesses operating at different resolutions. Moreover, generat-ing an appropriate ensemble of stochastic realizations mayrequire multiple simulations, each with different parametersettings corresponding to slightly different, but plausible,scenarios. As a consequence, running simulations and inter-preting simulation results (along with the real-world obser-vations) to generate timely actionable results are difficult.We argue that these challenges can be significantly alle-viated using a data-driven approach that addresses the fol-lowing fundamental questions:• Given a large parameter space and fixed budget of sim-ulations, can we decide which simulations to executein the ensemble? Can we revise the ensemble as wereceive a stream of real world observations?• Can we compare a large number of simulation en-sembles and observations (under different parameter(a) Query and exploration interface(b) Simulation visualization interfaceFigure 3: simVIZ simulation query, visualization,and analysis interfaces (instantiated with epidemicsimulation ensembles): visualizing an epidemic sim-ulation as a multi-variate time series and the keyrobust multivariate (RMT) events [9] identified ona given simulationsettings) to identify their similarities and differences?Can we analyze one or more simulation ensemblesto discover patterns and relationships between inputparameters, key events/interventions, and simulationoutcomes? Can we discover key events and summarizea large simulation ensemble to highlight these events?Can we classify these key events?• Can we search and explore simulation ensembles basedon the underlying key events or the overall simulationsimilarities? Can we keep track of the most relevantand most outlier simulations in an ensemble as we re-ceive a stream of real world observations?1.2 simDMS OverviewThe simDMS system (Figure 2) and its visualization enginesimVIZ (Figure 3) aim at assisting users to explore largesimulation ensembles while limiting the impact of aforemen-tioned challenges [4, 5, 6, 7, 8]. In particular, simDMS sup-ports• analysis and indexing of simulation data sets, includ-ing extraction of salient multi-variate temporal fea-tures from inter-dependent parameters (spanning mul-tiple layers and spatial-temporal frames, driven bycomplex dynamic processes operating at different reso-lutions) and indexing of these features for efficient andaccurate search and alignment;583FOR $p in fn:collection(’EpidemicSimulationEnsemble’) ^LET $diseaseModel := $p/project/scenario/model/disease ^LET $triggerModel := $p/project/scenario/trigger ^LET $epidemicScenario := $p/project/scenario ^WHERE$diseaseModel/transmissionRate <= 0.6 and$diseaseModel/transmissionRate >= 0.3 and$diseaseModel/recoveryRate = 0.5 and$triggerModel/@type="Vaccination" and($epidemicScenario/infector/@targetISOKey="US-CA" or$epidemicScenario/infector/@targetISOKey="US-NY" ) and($epidemicScenario/graph = "mobility_graph_7.xml" or$epidemicScenario/graph = "mobility_graph_8.xml") ^RETURN$diseaseModel/transmissionRate,$diseaseModel/recoveryRate,$epidemicScenario/graph ^STATE={AZ,CA,NM};MODEL={SEIR,SIR};PROPERTIES={Infected,Incidence,Deaths};FROM ={01/01/2012 12:00:00}; TO={08/31/2012 12:00:00};BY={1-D}; FUNCTION ={avg};Figure 4: A metadata query over an epidemic sim-ulation ensemble• parameter and feature analysis, including identificationof unknown dependencies across the input parametersand output variables spanning the different layers ofthe observation and simulation data. These, and theprocesses they imply, can be used for understandingand refining the parameter dependencies and models.Query and visualization interfaces for the epi-demic (epiDMS) and building energy (eDMS)instantiations of the simDMS platform can befound at http://aria.asu.edu/epidms andhttp://aria.asu.edu/edms, respectively. You can watch atutorial at https://youtu.be/9w-4nDhXv3k .2. DEMONSTRATION SCENARIOSWe will demonstrate the system on (a) epidemic simula-tion data sets created using the Spatiotemporal Epidemio-logical Modeler (STEM) [3] and (b) building energy simula-tion data sets, created using the EnergyPlus building energysimulation program [1]. The simulations will be stored insimDMS and will be visually analyzed during the demonstra-tion using simVIZ.2.1 Simulation Ensemble PlanningA simulation ensemble (consisting of a set of simulationinstances sampled from an input parameter space) can beseen as defining an outcome-surface for each of the out-put variables (such as the number of deaths that will re-sult from an epidemic): each outcome-surface describes theprobability distribution of the potential outcomes for thecorresponding variable. These simulation ensembles, con-sisting of potentially tens of thousands of simulations, areexpensive to obtain: therefore we need sampling strategiesfor the input parameter spaces that eliminate irrelevant sce-narios in such a way that more accurate simulation resultsare obtained where they are more relevant. Moreover, thesesimulation ensembles need to be continuously revised andrefined as the situation on the ground changes: (a) revisionsinvolve incorporating real-world observations into existingsimulations to alter their outcomes; (b) refinements involveidentifying new simulations to run based on the changingsituation on the ground. Therefore, we will demonstratedata-driven sampling strategies to decide (given a budgetof simulations) which simulations to run and incrementalnon-uniform sample-based data construction techniques torevise outcome-surfaces. We will specifically highlight howto assign utility- and cost-functions for each potential sample(based on how well the observed data are fitting the previ-ous simulations, how likely a new simulation at the givensample improves the accuracy of fit, and how costly the cor-responding simulation would be) and use these functions todecide the optimal re-sampling strategy.2.2 Scenario- and Similarity-based QueryingA basic function of the simDMS system is to retrieve simula-tions based on a user-specified scenario description. Figure 4presents a sample query:• The“FOR”statement allows the user to select the sim-ulation dataset to query. In this example, the userfocuses on the stored simulation set “EpidemicSimula-tionEnsemble”.• The “LET” statement allows the user to associatevariables representing disease and intervention triggermodels with epidemic scenarios.• The “WHERE” clause allows the user to specify con-ditions on the simulation models to filter those simu-lations that are relevant for the current analysis. Inthis example, the user specifies that for the returnedsimulations, the transmission rate parameter shouldbe between 0.3 and 0.6, the recovery rate parametershould be set to 0.5, and that a“vaccination” type trig-ger should be included in the simulation model. Theuser also specifies that epidemic should have startedin California (CA) or New York (NY) and the “mo-bility graph 7.xml” or “mobility graph 8.xml” shouldhave been used to generate the simulations.• The “RETURN” clause lists the simulation parame-ters to be returned in the result. In this example, theuser is interested in the transmission rate, recoveryrate, the mobility graph for each returned simulation.In addition, the query asks the system to return thetime series corresponding to the“infected”, “incidence”,and“deaths”simulation output parameters for Arizona(AZ), California (CA), and New Mexico (NM).• The user further specifies that s/he is interested in onlythe first 8 months of the simulation.• Finally, the user specifies that the system returns daily(1-D) averages of the simulation parameters for thespecified duration.Note that, in order to process this single query, simDMS com-bines data of different forms (structured, semi-structured,and temporal), stored in different back-end storage engines.In addition to scenario-based filtering and search, the plat-form also enables searching and/or triggering based on par-ticular temporal patterns on the ensembles. This feature584(a) Exploration hierarchy(b) Comparing selected simulations(c) Metadata comparison (d) Feature-based comparisonFigure 5: Sample interfaces for exploring ensemblesallows the expert to identify relevant subsets of stored simu-lations that match actual real-world observations or specifictargets for intervention measures.2.3 Analysis and Exploration of EnsemblesOnce the query is executed and the relevant simulationsare identified, the system then organizes the results into anavigable hierarchy, based on the temporal dynamics of thesimulation results (Figure 5): Since simulation data sets canbe viewed as multi-variate time series, simVIZ focuses on vi-sual analysis (e.g. event detection, similarity and differenceanalysis) of single and multiple multi-variate simulation datasets. Scenarios that result in similar patterns are groupedunder the same branch, while simulations that show majordifferences in disease development are placed under differ-ent branches of the navigation hierarchy. The user can thennavigate this hierarchy using “drill-down” and “roll-up” op-erations and pick sets of simulations to study and comparethe corresponding scenarios in further detail.The interface presents both conventional series plots aswell as heatmap visualizations, where each series is shownas a row of pixels. It is important to note that, whilethe temporal (i.e., horizontal) axis is ordered, the verticalaxis corresponding to the different states is not ordered, inthat two nearby states according to user mobility may notbe neighboring rows on the interface due to the complex-ity of the mobility graph. The interface also highlights,on the heatmap, the major robust multi-variate time se-ries (RMT) features (optimized for supporting alignmentsof multi-variate time series, leveraging known correlationsand dependencies among the variates [9]) identified on theheatmap. An RMT feature is a part of the time series thatis different in structure from its immediate context in timeand/or variate relationships. A key property of these RMTfeatures is that they are robust against noise and commontransformations, such as temporal shifts or missing variates.This is illustrated in Figure 5(d), which shows two differ-ent epidemic simulations, with the same starting state, butdifferent disease parameters and interventions. While the re-sulting disease evolutions are visibly different in shape, thesame multi-variate feature (corresponding to the onset ofthe disease on the same nearby states) is identified on bothsimulations. This robustness property of RMT features en-ables various simVIZ functions, such as search, clustering,classification, and summarization of simulations and largesimulation data sets [4, 5, 6, 7, 8].3. CONCLUSIONSThe simDMS platform provides metadata and event-drivenanalysis and visualization of simulation ensembles to assistdecision makers to query and explore ensemble simulationsand decide which additional simulations to execute.4. ACKNOWLEDGMENTSWe also thank Reece Bailey, Hans Behrens, Sarah Fallah-Adl, Ashish Gadarki, Anisha Gupta, Divyanshu Khare,Mao-Lin Li, Nancy Li, Samuel Morton, Shivam Sadakar,Xiaolan Wang, Fiona Zhang, Luke Zhang, and Jerry Zhufor their contributions to the simDMS platform.5. REFERENCES[1] EnergyPlus Energy Simulation Software.http://apps1.eere.energy.gov/buildings/energyplus/[2] The Global Epidemic and Mobility Model, GLEAM.http://www.gleamviz.org[3] STEM. The Spatiotemporal Epidemiological ModelerProject. Available at http://www.eclipse.org/stem.[4] S. Huang, K.S. Candan, M.L.Sapino. BICP:Block-Incremental CP Decomposition with UpdateSensitive Refinement. CIKM 2016: 1221-1230[5] X. Li, S. Huang, K.S. Candan, M.L. Sapino. 2PCP:Two-phase CP Decomposition for Billion-Scale DenseTensors. ICDE 2016: pp. 835-846.[6] S. Liu, Y. Garg, K.S. Candan, M.L. Sapino, G.Chowell. NOTES2: Networks-Of-Traces for EpidemicSpread Simulations. AAAI Workshop on ComputationalSustainability, 2015.[7] S. Liu, S. Poccia, K.S. Candan, G. Chowell, and M.L.Sapino. epiDMS: Data Management and Analytics forDecision-Making From Epidemic Spread SimulationEnsembles. Journal of Infectious Diseases. 2016.[8] P. Nagarkar, K.S. Candan, A. Bhat. CompressedSpatial Hierarchical Bitmap (cSHB) Indexes forEfficiently Processing Spatial Range Query Workloads.PVLDB 8(12): 1382-1393 (2015)[9] X.Wang, K.S. Candan, and ML Sapino. LeveragingMetadata for Identifying Local, Robust Multi-VariateTemporal (RMT) Features. ICDE 2014, 388-399.585

SIMDMS: Data Management and Analysis to Support Decision Making through Large Simulation Ensembles

https://iris.unito.it/bitstream/2318/1647179/1/paper-438.pdf

SIMDMS: Data Management and Analysis to Support Decision Making through Large Simulation Ensembles

Abstract

Similar works

Full text

Available Versions

Institutional Research Information System University of Turin