Classifications

• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
  • G06Q10/00—Administration; Management
  • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Definitions

• the invention relates to computer predictions of atmospheric contaminant transport, with much faster than real-time access to results for emergency response and enables accurate new capabilities to locate and unknown source from isolated sensor readings of the contaminant.
• CBR chemical, biological, or radiological
• the existing plume prediction technology in use throughout the nation (and the world) is based on Gaussian similarity solutions and/or Lagrangian particle models.
• the Gaussian similarity solutions (“puffs” or “plumes") are extended Lagrangian approximations for a single particle or "puff.”
• These Gaussian puff models only really apply for large regions and flat terrain where large-scale vortex shedding from buildings, cliffs, or mountains is absent.
• the Lagrangian particle models treat many more moving points, and thus take much longer to compute, but they do not allow each particle or "puff to expand.
• the CBR defense of a fixed site or region has a number of important features that make it different from the predictive simulation of a contaminant plume from a l ⁇ iown set of initial conditions. The biggest difference is that very little may be known about the source, perhaps not even its location. Therefore any analysis methods for real-time, emergency response cannot require this information. It is crucial to be able to use anecdotal information, qualitative data, and any quantitative sensor data that may be available and then to instantly build a situation assessment suitable for immediate action based on this fragmented information.
• a software emergency assessment tool should be effectively instantaneous and easy to use to allow immediate assimilation of new data, instantaneous computation of exposed and soon-to-be exposed regions, and zero-delay evaluations of options for future actions.
• the software should also be capable of projecting optimal evacuation paths based on the current and evolving situation assessment. To meet these crutial requirements, a new tool is required that is both much faster than current "common use” models and accurate compared to three-dimensional, physics-based flow simulations for scenarios involving complex and urban landscapes.
• the data structure and algorithms disclosed here focus on situation assessment through sensor fusion of qualitative and incomplete data using a summary of accurate flow details, rather than integrating a computer simulation.
• ALOHA Absolute Locations of Hazardous Atmospheres
• ALOHA allows the user to estimate the downwind dispersion of a chemical cloud based on the toxicological/physical characteristics of the released chemical, atmospheric conditions, and specific circumstances of the release.
• Graphical outputs include a "cloud footprint" that can be plotted on maps.
• MIDAS-ATTM Metalological Information and Dispersion Assessment System - Anti-Terrorism
• MIDAS-AT is software that models dispersion using a Gaussian puff model of releases of industrial chemicals, chemical and biological agents, and radiological isotopes caused by accidents or intentional acts.
• MIDAS-AT is designed for use during emergencies and for planning emergency response drills.
• VLSTRACK (Vapor, Liquid, and Solid Tracking) provides approximate downwind hazard predictions for a wide range of chemical and biological agents and munitions of military interest using a Gaussian puff model.
• CATS Consequences Assessment Tool Set
• SCIPUFF Hazard Prediction and Assessment Capability
• HP AC Hazard Prediction and Assessment Capability
• HP AC estimates the effects of hazardous material releases into the atmosphere.
• the HP AC system also predicts approximate downwind hazard areas resulting from a nuclear weapon strike or reactor accident.
• NARAC National Atmospheric Release Advisory Center
• emergency response central modeling system consists of a coupled suite of meteorological and dispersion models, including Gaussian, Lagrangian and computational fluid dynamics (CFD) models. Users must initiate a problem through a phone call to their operations staff or interactively via computer.
• CFD computational fluid dynamics
• NARAC will then execute a combination of Gaussian puff/plume and approximate Lagrangian particle 3-D models to generate the requested products that depict the size and location of the plume, affected population, health risks, and proposed emergency responses.
• NARAC has recently announced a multi-year research program aimed at adding an urban modeling and sensor fusion capability.
• PEAC® Palmtop Emergency Action for Chemicals
• PEAC® Palmtop Emergency Action for Chemicals
• PEAC® can return results within seconds and requires less detailed knowledge of the source, but the resulting fixed-shape plume does not take into account any effect of complex terrain or buildings.
• none of these systems allow a user to backtrack the very limited input data likely to be available in a terrorist scenario to an undisclosed source location.
• FASTD-CT FAST3D - Contaminant Transport
• the fluid dynamics is performed with a fourth-order accurate implementation of a low-dissipation algorithm that sheds vortices from obstacles as small one cell in size.
• the region of interest or domain is divided into finite parcels called computational cells. These units limit the spatial accuracy of a particular model. Particular care has been paid to the turbulence treatments since the turbulence in the urban canyons lofts ground-level contaminant up to where the faster horizontal airflow can transport it downwind.
• FAST3D-CT also has a number of physical processes specific to contaminant transport in urban areas such as solar heating and buoyancy, solar chemical degradation, evaporation of airborne droplets, re-lofting of particles and ground evaporation of liquids.
• a data structure with associated algorithms for representing and reconstructing contaminant transport in realistic complex geometries, e.g. cities and buildings, is described.
• the data structure has a pair of two-dimensional matrices compressed and extracted from three-dimensional data sets describing the transport and dispersion of contaminants in complex geometries using full-resolution, time-dependent computational fluid dynamics or compete and detailed experimental data.
• a first matrix has numerical values, in which continuous, monotone, gauge-invariant contours of the same value represent the edge of a family of contaminant clouds transported from one edge of the domain.
• a second matrix contains numerical values, in which continuous, monotone, gauge-invariant contours of the same value represent the opposite edge of a family of contaminant clouds transported from the other edge of the domain.
• This dual-matrix data structure is used by overlapping the two matrices within the domain, and reading the areas between the two edge contours at particular locations of interest.
• the invention provides a system of storing and manipulating data within the two matrices to give much faster than real time data for first responders to CBR attacks.
• the invention provides a method of manipulating data using the two matrices to enable quick response to imminent threat of CBR attacks.
• Figure 1 Overview of how to develop and use the data structure.
• Figure 2. A depiction of the data structure (grey scale substituted for full colour representation) showing the dual-matrix and the continuous nature of the data stored within the structure.
• Figure 3 (a) A left edge matrix (half of the data structure), (b) A right edge matrix (other half of the data structure).
• This figure uses contour lines (as on a topographic map) rather than continuous shading to show the values in the two matrix components of the data structure.
• the building/tree map not actually part of the data structure, shows the relationship of the dual-matrix values and contours to the detailed geometry on which the dual matrix is based.
• Figure 4 An overlay of the two edge matrices for the complex geometry shown, showing how the dual-matrix is used.
• Figure 6(a) Plume envelope and footprint of a contaminant source and upwind danger zone of a site of interest for the full geometry of buildings, trees, and terrain.
• Figure 6(b) Same as Figure 6(a) using the flat earth dual-matrix of Figure 5. Note the simple symmetry of the plume, footprint, and danger zone about the wind direction and the qualitative differences of the results when the full geometry is considered.
• Figure 7 Plots of the Figure of Merit using the data structure versus time for six different sets of full 3D simulations.
• Figure 8. Comparison of the accuracy and costs of Experiment, Computational Fluid Dynamics, Dual-Matrix Reconstruction, Puff/Plume Models, and Simple Phenomenologies .
• Figure 9 Comparison of the amount of data processed, time to run the models, and time required to learn to use the technologies for, CFD, Dual-Matrix Reconstruction, and Gaussian puff models.
• Figure 10 A simple block diagram of the routine controlling access to and use of the data structure and associated algorithms.
• Figure 11 A logical flow diagram of the procedure that processes sites of interest using the data structure.
• Figure 12 A logic flow diagram of the procedure that processes contaminant observations using the data structure.
• Figure 13 A logic flow diagram of the procedure that processes contaminant sources using the data structure.
• Figure 1 is an overview of the development process 160 and method of using 150 the dual matrix data structure.
• the data structure comprises a pair of two-dimensional matrices compressed and extracted from three-dimensional matrices describing the transport and dispersion of contaminants in complex geometries using full-resolution, time-dependent computational fluid dynamics or compete and detailed experimental data.
• One embodiment of the dual matrix data structure is a NomographTM. Algorithms using the data structure are both faster and more accurate than prior methods of forecasting transport and dispersion in complex geometries.
• the region of interest is defined and an extensive database 130 of CFD results is assembled.
• a number of different wind directions around the compass must be simulated and results for each are tabulated in the three dimensional database 130. While the figures and description refer to CFD results generated on a supercomputer, if sufficiently complete and accurate experimental results were available they could also be used to develop the new data structure.
• the CFD result database for the region of interest is next compressed, 140.
• compression 140 by a factor of about 10,000, extraneous and redundant information from the three-dimensional, time-dependant CFD results is stripped.
• this database compression 140 detailed records the contaminant flow paths and turbulent dispersion for the urban area in question (region of interest) are reduced to the continuously nested and monotonic set of left and right cloud edge contours shown in Figures 2 through 5.
• the continuously nested set of left and right cloud edge contour levels, for all wind directions comprise the dual-matrix data structure described herein. In addition to the contour levels being continuous, the transition from one contour to the next is continuous.
• the dual- matrix data structure and associated algorithms allow users to view possible source areas based on backtracking contaminant observations to determine the unknown source location; view recommended evacuation routes; and view upwind danger zones, instantly.
• This data structure, together with the set of algorithms to access and process the data structure, are assembled in the dual-matrix library 100.
• the dual-matrix library 100 includes a pair of two-dimensional matrices holding continuous contour values indicating the location of both edges of the dispersed cloud of contaminant within a region of interest. Users access the contour flow map dual-matrix library 100, with its embedded emergency-response information through a graphical user interface (GUI) 110, but other applications can call the control procedure 120, directly for autonomous operation.
• GUI graphical user interface
• the control procedure 120 is described more fully below with respect to Figure 8.
• the control procedure 120 is a procedure that interprets scenarios defined in their entirety by two sets of numbers chosen by the GUI or by the user application program.
• the control procedure 120 is a software routine.
• One set of numbers, the environmental state vector gives all the necessary information about the wind, all the settings for the various displays required, and specifies the instantaneous time after release of the scenario.
• the second set of numbers, the node state vector lists the location and status of all the nodes such as what kind of node it is, whether it is turned on, and the like.
• the term node is used here to refer to a particular computational cell (location) where additional information is available or desired.
• a node is processed by using the information stored in the data structure, combined with the algorithms to use the data structure to give the user the requested information about the node.
• These two vectors, or arrays allow the control procedure to oversee preparing the requested displays using the composite flow map data structure and algorithms in the library 100.
• control procedure 120 will, depending on the environmental state vector, invoke an algorithm to process the complete set of sensor data 170, the information that contaminant is or is not observed at a particular location within the region of interest.
• Sensor data may be anything from an anecdotal observation or news report, to readings from specialized sensors showing quantitative contaminant level.
• the control procedure 120 also can invoke an algorithm to process the site nodes 180, based on data in the environmental state vector.
• a site is a location of critical interest and may be permanent or transitory. For example a site may be a particular building or outdoor location, a parade route, or location of a sporting event.
• the algorithm to process sites 180 displays data regarding the danger zone for a particular site or location within the region of interest.
• the process sites algorithm 180 can also display the contaminated region resulting from a leak within a particular building.
• the danger zone is the region upwind from which a contammant could reach the site in question.
• the control procedure 120 depending on the environmental state vector, can invoke an algorithm to process sources 190, specific locations where contaminant has been introduced into the region of interest.
• use of the library 100, control procedure 120 and user interface 110 comprise an application of the dual-matrix for emergency response 150.
• Figure 2 is a graphical grey scale depiction of the dual matrix data structure showing the two component matrices dual -matrix and the continuous nature of the data stored within this structure for one particular region of interest, or domain 200, and wind direction indicated by the arrow 270.
• the various gray shadings at each location within the domain 200 correspond to different values within each edge matrix. While the domain 200 is shown as a square area, it need not be square or even quadrilinear. Any contiguous shape, including irregular shapes, can be chosen for the region of interest.
• the data structure comprises two matrices of numbers, shown graphically as gray shadings, representing left edge values in the domain 200 and right edge values in the same domain 200.
• the continuous variation of the shading represents the continuity of the numbers stored in the matrices. It should be noted that the domain 200 must be the same shape and location for both the left edge values and right edge values.
• the left edge matrix of values in Figure 2 is also illustrated with contours 240, at evenly spaced intervals throughout the domain 200.
• the contours begin at the edge of the domain 200 evenly spaced both physically and numerically, but the shape and physical distance between the contours vary in response to the complex geometry of the domain 200, leaving only the contour values evenly spaced numerically.
• the right edge matrix of values is marked with contours 230 throughout the domain. These contours correspond to the left and the right edges of a family of contaminant clouds.
• the contours, as well as the buildings and trees, are shown only for ease of interpreting the figures; they are not an explicit part of the composite flow map data structure.
• the processing algorithms use only the continuous array of numbers in the left and right edge matrices, allowing the processing algorithms to find a contaminant cloud edge at any point within the domain.
• the lowest left edge matrix value is located in the upper right near 210 of the domain 200, while the highest left edge matrix value is located in the lower left region near 220 of the domain 200.
• the lowest right edge matrix value is located in the right region 260 of the domain 200, while the highest right edge matrix value is located in the upper left region 250 of the domain 200.
• the contours 230 and 240 along with the shading, within each matrix the contours never cross and are defined by monotone values.
• the resulting dual-matrix data structure is significantly smaller than it would otherwise be.
• the monotonicity requirement is an approximation of the original data that allows for a greater safety margin between the contaminant cloud edges modeled using the dual-matrix data structure and the actual cloud edges from full-resolution CFD results.
• Figure 3 a is a black and white enlargement of the depiction of the left edge matrix or half of the data structure in Figure 2, showing only contours for ease of reading.
• the contour levels are equally spaced in value; their changing direction and separation indicates graphically how the new data structure captures the underlying urban geometry.
• a family of contaminant clouds, all sharing the same left edge is defined by each contour.
• the overall domain 200 of the data structure is shown, including the location of selected buildings 310 and trees 320.
• the domain in Figure 3a corresponds entirely to domain 200 in Figure 2.
• the trees 320 and buildings 310 are only shown for reference, to indicate the complexity of the region of interest through which contaminant transport can be accurately simulated.
• the left edge contour lines 330 are also shown, along with areas where the flow concentrates 340 and other areas where the flow expands 350.
• Figure 3b is a black and white enlargement of the depiction of the corresponding right edge matrix, the second half of the data structure in Figure 2, showing only the contours for ease of reading. Another family of contaminant clouds is defined by each of these right edge contours.
• the domain 200 of the data structure is shown, including the location of buildings 310 and trees 320. This corresponds to domain of Figure 2.
• the right edge contour lines 330 are also shown, along with areas where the flow concentrates 340 and areas where the flow expands 350.
• Figure 4 is an overlay of the left and right edge matrices shown in Figures 3 a and 3b, for the complex geometry of the domain 200.
• the wind is shown by a downward arrow.
• Overlaying the left and right edge matrices, here displayed as contours illustrates geometrically how the upwind and downwind sectors, originating at the intersection of one right edge contour and one left edge contour, can be defined for any particular node location, in the domain of the data structure. While the contours shown here only intersect at a few points within the domain, the continuous nature of the values within the dual-matrix data structure ensures that a pair contours can be constructed at any point within the domain 200.
• the domain 200 must be identical for both the right and the left edge matrices in order to overlay them.
• a site 410 is depicted as square icon representing a building, facility or other location of special interest needing to be analyzed or protected.
• a danger zone 440 upwind of the site of interest is defined. Since the expansion of a contaminant cloud is contained within the cloud edge lines, any contaminant outside of this danger zone 440, cannot reach the site 410 in the current wind condition. While the site 410 is shown in one particular place, users of the data structure may place any number of sites at any location or locations within the region of interest 200.
• the location of a source of contamination 450 is marked with a star shaped icon.
• the source of contamination 450 may be an accident, sprayer, broken container of hazardous chemicals or the like that needs to be tracked and analyzed in order to coordinate emergency response.
• the source processing algorithm identifies the eventual contamination footprint 480.
• a contamination footprint is the possible extent of contamination from a source given enough time to disperse downwind throughout the region of interest.
• the downwind contamination footprint 480 of the source 450 is defined by the intersection of the right 470 and left 460 edge contours passing through the source location 450.
• the actual contaminant plume envelope starts at the source 450, in the center of this footprint 480, and expands away from the source 450 and towards the left 460 and right 470 edges with time.
• the plume envelope display is not shown in Figure 4 or Figure 5 below but is illustrated in Figure 6.
• Figure 5 is an overlay of the left and right edge matrices of the composite data structure when the complex geometry shown here, and in previous figures, is in fact neglected in assembling the data in the dual-matrices.
• Such contours represent an accurate solution to the overall contaminant transport problem in large, flat regions of interest, such as a desert or empty plain.
• Figure 5 represents an approximation, using a dual-matrix, to the solutions obtained by using Gaussian puff/plume methods to compute transport and dispersion.
• a square icon represents a site 510, such as a building, facility or other location of interest for protection or analysis.
• a danger zone 540 upwind of the site of interest is defined.
• the danger zone 540 is the upwind area between the left edge contour 520 and the right edge contour 530 that intersect at the site of interest 510. Since the cloud-edge lines bound the contaminant cloud motion, any contaminant outside of this danger zone 540, is predicted to be unable reach the site 510 in the current wind conditions. However, in the complex geometry of the domain 200, the simulation that does not account for the domain geometry significantly under-predicts the danger zone. Comparing the danger zone 440 with the danger zone 540 shows this.
• the location of the source of contamination 550 is marked with a star shaped icon.
• the algorithm identifies the predicted contamination footprint 580.
• the downwind contamination footprint 580 of the source 550 is defined by the intersection of the right 570 and left 560 edge contours passing through the source location 550.
• Figure 6a shows the reconstructed plume envelope 650 at 10 minutes after contaminant release, the entire footprint 660 of a source node, shown as star icon 600, and the upwind danger zone 640 of a site node, shown as a square icon 610.
• the assessment of the extent of these critical regions illustrates the method of use of the dual- matrix right edge and left edge matrices of Figure 4 capturing the contaminant flow using the full geometry of buildings, trees, and terrain.
• Figure 6b shows the reconstructed plume envelope 655 after 10 minutes and the entire footprint 665 for the same site 610 and source 600 in the same domain - but captured using flat-earth dual-matrices with no building or tree geometry, as illustrated in Figure 5. Note the simple symmetry of the plume, footprint, and danger zone about the wind direction in this latter case and the qualitative differences of the results when the full geometry is considered.
• Sensors are indicated in these figures with a triangular icon.
• Sensors labeled 620 have not been triggered by contaminant at the time depicted in Figures 6a and 6b as they lie outside of the instantaneous plume envelope 650 and 655.
• Sensor 625 is located outside the footprint using no geometry but inside the footprint when full geometry is used.
• a simulation such as a puff or plume model, that does not account for complex geometry gives a false negative result for the site 610, possibly resulting in high strategic losses.
• a measurable contaminant level in the plume envelope has triggered sensors 630 in Figure 6a but sensors 635 at the same location in Figure 6b have not yet been contaminated.
• an emergency response tool that does not account for complex geometry results in false negatives.
• sensor responses clearly differ depending on the configuration and situation of buildings and trees using this new data structure.
• the shaded regions 640 in Figure 6a and 645 in Figure 6b estimate the upwind danger zone for the site of interest 610 with and without capturing the effects of the geometry, respectively.
• the danger zone for a location, or site is the set of all possible positions upwind where a source of contamination could reach the indicated site. This important information is an entirely new capability made possible by the dual-matrix and algorithms for its use.
• the plume envelopes 650 and 655 show the geographical region that the contaminant plume may reach during its continuing expansion at the indicated time after release.
• the contamination footprints 660 and 665 in Figures 6a and 6b respectively, represent the full extent of the growing contamination region after the plume envelope has spread to its maximum toxic extent according to dual-matrix simulations with and without underlying urban geometry.
• Figure 7 plots the figure of merit for the dual-matrix versus time for six different sets of full 3D CFD simulations representing different contaminant release scenarios.
• the figure of merit is a quantitative comparison between the two sets of results. A perfect match would be indicated by a figure of merit of 1.00, or 100% match.
• the results from the dual-matrix match those obtained using full 3D CFD calculations to about 70-90% accuracy.
• the dual-matrix results reach to between approximately 75-85% accuracy. And with time, the dual-matrix results converge to about 80% accuracy.
• Figure 8 is a table comparing accuracy and other capabilities of experiments, full 3D Computational Fluid Dynamics (CFD), the dual-matrix data structure and associated algorithms, Gaussian or Lagrangian puff/plume models, and simple phenomenologies.
• Phenomenologies are models representing dispersion in an approximate form based on physical intuitive, qualitative understanding of the underlying physics. The generally more accurate approaches are listed in the table above the generally less accurate approaches. The different approaches are evaluated with respect to the quality of approximation of the physics models and input, fluid dynamics and turbulence, boundary conditions, realistic geometry; the variability and uncertainty of the results; whether the approach allows sensor fusion and backtracking; and whether evacuation routes are provided.
• Sensor fusion is the process by which sensor data is combined with simulations based on a potential location of a contaminant source to increase the accuracy of the simulation results.
• a rating of excellent indicates that the area of concern is quantitatively accurate.
• a rating of very good means that the results can still be classified as quantitatively accurate, but there are some aspects of the model that could benefit from refinement.
• a rating of good indicates that more refinement is needed, although the results could still be classified as quantitatively correct.
• a rating of fair to good indicates that a good semi- quantitative understanding is available from the approach.
• a rating of fair indicates that good qualitative understanding may be obtained from the approach.
• a rating of poor indicates that acceptable qualitative agreement is the best that can be expected. The indication that the area is problematic means that the approach does not adequately address these areas.
• Figure 9 is a table comparing the amount of data processed, the time required to run the models, and the approximate time required to learn how to use the respective technologies for CFD, the dual-matrix and associated algorithms, and Gaussian puff models.
• the dual-matrix is approximately 1000 times faster than even the typical implementation of the Gaussian puff models, because it does not actually calculate the contaminant transport at the time of use, but compactly stores highly accurate results for effectively immediate recall when needed.
• the dual-matrix and associated algorithms are approximately 100 times faster to learn to use than current operational models. This allows more emergency response personnel to be trained to use the system without taking excessive time from other duties.
• the data requirements for the dual-matrix are very small and approximately 10 times less than the requirements of current operational models, and could be conveniently added to personal digital assistants and the like, for ease of field use.
• Figure 10 shows a logical flow diagram of the master control procedure 120 controlling access to the dual-matrices 100 and the specific processing algorithms, or procedures 170, 180, 190 of Figure 1.
• a user, client program or Graphical User interface connects to the display and analysis capabilities through this master control procedure 120.
• a number of different user programs are envisioned using this master control procedure.
• a user working through a Graphical User Interface (GUI) 110, or other user application program initiates the processing algorithms by invoking or "entering” the master control procedure 120.
• the first step in one embodiment is to enter 1000 the master control procedure 120.
• a software embodiment of the method of using the dual- matrix is a "user" or "client.”
• One embodiment of a client is CT- AnalystTM available from the U.S. Naval Research Laboratory.
• Other embodiments of a client are described in U.S. Patent Application No. 10/612582 filed July 2, 2003.
• the master control procedure begins processing with a valid environmental state vector and a node state vector constructed externally. These two state vectors define the current scenario, or system, and the displays that are to be produced and displayed for the user by the GUI or returned for further analysis to the client program 1080.
• Moveable nodes are located within the domain and defined in the node state vector.
• the master control procedure processes the moveable nodes contained in the node state vector one at a time, considering site, sensor, and source nodes. First, a check 1010 is performed to see if all the nodes have been processed.
• the master control procedure checks to see if any of the remaining nodes are site nodes 1020. If there are any site nodes remaining, the requested displays and analyses are considered, or processed, 170 one at a time.
• the processing of each site node 1030 is illustrated in detail in Figure 9. If there are any sensor nodes to be processed 1040, the requested displays and analyses are considered, or processed, 180 one at a time.
• the processing of each sensor node 1050 is illustrated in Figure 10. If there are any source nodes to be processed 1060, the requested displays and analyses are considered, or processed, 190 one at a time.
• the processing of each source node 1070 is illustrated in Figure 11.
• FIG. 11 shows a logical flow diagram of the algorithm that processes site nodes using the disclosed dual-matrices to construct the upwind danger zone and downwind leakage zone of the sites.
• the master control procedure 1000 invokes "Enters" 1100 the site processing procedure, it provides the dual-matrix corresponding to the current wind speed and direction and passes the current environmental and node state vectors defining the scenario to the site processing procedure.
• the left and right edge matrices are invariant under addition or multiplication, that is they are gauge invariant, the matrices for different wind directions or other environmental conditions may be blended or "morphed" to describe a continuum of different conditions.
• the dual-matrix for the current wind speed and direction, or even the current foliage coverage, atmospheric stability may be a blend of two or more dual-matrices constructed from the original CFD or experimental data.
• the site processing procedure 1030 first selects 1105 the next unprocessed, site node from the list of nodes or node vector. For example the selecting step 1105 could select node 410 in Figure 4, or 510 in Figure 5.
• the procedure next checks 1110 whether the danger zone display is to be constructed for site nodes, in other words, if the danger zone is "turned on" in the environmental state vector. If it is, the danger zone computation is initialized 1115.
• the data items defining the danger zone are initialized for the selected site node by finding the left and right edge values at the location of the site node (e.g. node 410 or node 510). To be specific, the contour level values from the right edge
• Each of the cells in the overall domain is "read” from the library one at a time 1120.
• the local R(X,Y) and L(X,Y) contour level values are read from the composite data structure and evaluated using Formula 1.
• Step 1125 uses Formula 1, below to determine whether the cell X,Y lies in the danger zone.
• step 1125 checks to see if the right and left edge matrix values at the site are equal to or less than the right and left edge values at all locations within the domain.
• Formula 1 in the computer programming language FORTRAN is:
• inDangerZone r(xs,ys) .le. r(x,y) .and. l(x,y) .le. l(xs,ys)
• step 1135 goes back to step 1120 and repeats the process on the next cell in the domain of the data structure until the entire danger zone has been constructed.
• the site processing procedure 1030 checks the environmental state vector in step 1140 to determine whether the leakage display is to be constructed for site nodes.
• the leakage zone is the area around and a down wind of a site that may become contaminated if a contaminant leak occurs within or immediately adjacent to the site. If it is, the leakage zone computation is initialized in step 1145.
• the data items defining the leakage zone are initialized for the selected site node by finding the left and right edge values at the location, or computational cell, of the site node, e.g. node 410.
• Each of the cells in the overall domain is again “read” from the library one at a time 1150.
• the R(X,Y) and L(X,Y) contour level values are read from the composite data structure and evaluated using Formula 2.
• Step 1155 uses Formula 2, below, to determine whether the cell X,Y lies in the leakage zone.
• step 1155 checks to see if the right and left edge matrix values at the site are equal to or greater than the right and left edge values at all locations within the domain.
• Formula 2 is one of two formulas, the downwind formula, used to extract information from the dual-matrix. If all the cells in the dual-matrix have not been scanned, step 1165 goes back to step 1150 and repeats the process on the next unprocessed cell. Once all the sites have been processed from the node list, return to the master control procedure (1170 and Figure 8).
• Figure 12 shows a logical flow diagram of the algorithm that processes sensor nodes using the data structure to construct the upwind backtrack analysis and downwind consequence zone of the sensor observations. Sensors may be either "hot,” when something dangerous is observed, or "cold,” i.e. sensing only clear air.
• the master control procedure 1000 invokes or "enters” the sensor processing procedure 1050 in step 1200, it provides the data structure corresponding to the current wind speed and direction and passes the cu ⁇ ent environmental and node state vectors defining the scenario into the overall procedure 1050.
• the procedure 1050 first selects 1205 the next unprocessed, sensor node from the list of nodes. For example the selecting step 1205 could select one of nodes 620, 525, 630, and 635 in Figures 6a and 6b.
• the procedure next checks whether the backtrack analysis display is to be constructed 1210 for the set of sensor nodes. If the backtrack analysis is turned on in the environmental state vector, the backtrack analysis computation is initialized 1215 by calculating the wind direction and establishing a radius around each sensor location within which any detection of contaminant will trigger the simulated sensor. This initialization also sets the back count a ⁇ ay to zero at every cell in the domain.
• the data that defines the backtrack zone are initialized for each sensor node by finding the left and right edge values at the location of the sensor node.
• a cell is a possible source location if the cell lies within the upwind backtrack zone of two or more "hot” sensors but lies outside the backtrack zone of all "cold” sensor.
• the backtrack analysis is based on the full set of sensor nodes, or observations, acting in concert, rather than on each node separately but additively as in the contamination footprint, danger zone, leakage zone, and consequence zone displays described elsewhere.
• Each of the cells in the overall domain is "read” from the library one at a time in step 1220.
• the local R(X,Y) and L(X,Y) contour level values are evaluated from the composite data structure using Formula 1, below.
• Step 1225 uses Fo ⁇ nula 1, below to determine whether the cell X,Y lies in the backtrack zone.
• step 1230 checks if sensor n is "cold” in step 1230. If the cell is within the backtrack zone, or in the upwind region of a sensor or with a safety radius of any "cold" sensor, that cell cannot be the location of a source so the back count a ⁇ ay value is set to -1 in step 1235. In other words, step 1225 checks to see if the right and left edge matrix values at the sensor are equal to or less than the right and left edge values at all locations within the domain.
• Formula 1 in the computer programming language FORTRAN is:
• inBacktrackZone r(xs,ys) .le. r(x,y) .and. l(x,y) .le. l(xs,ys)
• Formula 1 is used to extract information from the dual-matrix and is applied for all sensor nodes for all cells in the entire data structure.
• step 1240 transfers the back count a ⁇ ay cell value to step 1245, where it is incremented by one.
• back count is an a ⁇ ay of numbers counting how many different "hot" sensors indicate each location, or cell, might be the source of contamination as long as there is no sensor showing, according to a back count value of-1, that the given cell cannot, in fact, be the source location.
• step 1250 checks, in step 1255, whether the consequence zone display is to be constructed for "hot” sensor nodes. If the consequence zone is "turned on” in the environmental state vector, the consequence zone computation is initialized in step 1260 by setting the consequence count a ⁇ ay everywhere to zero. The consequence zone is the total region downwind of all sensors that are reading a dangerous level of contaminant at their location, i.e. are "hot.”
• the data defining the consequence zone are initialized for the selected "hot" sensor node or nodes by finding the left and right edge values at the location of each "hot” sensor node.
• the contour level values from the right edge matrix are initialized for the selected "hot" sensor node or nodes by finding the left and right edge values at the location of each "hot” sensor node.
• Each of the cells in the overall domain is again "read” from the library one at a time 1265.
• the local R(X,Y) and L(X,Y) contour level values are read from the composite data structure and evaluated using Formula 2, below.
• Formula 2 is used to determine whether X,Y lies in the consequence zone.
• step 1275 sets the consequence count a ⁇ ay value to 1 at that location. In other words, step 1275 checks to see if the right and left edge matrix values at the sensor are equal to or greater than the right and left edge values at all locations within the domain.
• inConsequenceZone r(xs,ys) .ge. r(x,y) .and. l(x,y) .ge. l(xs,ys)
• step 1280 goes back to step 1265 and repeats the process on the next unprocessed cell.
• step 1285 returns to master control procedure of Figure 8.
• Figure 13 shows a logical flow diagram of the algorithm that processes source nodes using the disclosed composite data structure to compute the contamination footprint.
• the composite flow map data structure shown in Figures 2 through 6, and formula (1) is used to determine all points within the contaminant footprint (the gray regions 660 and 665 in Figure 6).
• the process sources procedure 1070 When the overall control procedure 1000, invokes or "Enters," in step 1300, the process sources procedure 1070, it provides the composite data structure co ⁇ esponding to the cu ⁇ ent wind speed and direction and passes through the cu ⁇ ent environmental and node state vectors defining the scenario.
• the data controlling this processing is initialized 1310 and the plume cell count a ⁇ ay is set to zero at all of the cells in the data structure domain.
• the initialization of step 1310 includes determining whether the footprint has to be recalculated since the last time the control procedure was invoked.
• Each of the NX by NY cells in the composite data structure is read from the library one at a time 1320.
• the local R(X,Y) and L(X,Y) contour level values are read from the composite data structure and evaluated using Formula 1, below.
• Step 1330 uses Formula 1 determine whether X,Y lies in the contamination footprint.
• InFootprintZone r(xs,ys) .le. r(x,y) .and. l(x,y) .le. l(xs,ys)
• inFootprintZone is a logical variable. If the escape display is active in the environmental state vector according to step 1350, the escape route lines are drawn perpendicular to the locally prevailing wind on top of the footprint display 1360. These lines show the optimal direction to walk, which is away from the local plume centerlme, to avoid exposure.

Landscapes

• Engineering & Computer Science (AREA)
• Business, Economics & Management (AREA)
• Human Resources & Organizations (AREA)
• Economics (AREA)
• Strategic Management (AREA)
• Marketing (AREA)
• Game Theory and Decision Science (AREA)
• Entrepreneurship & Innovation (AREA)
• Development Economics (AREA)
• Operations Research (AREA)
• Quality & Reliability (AREA)
• Tourism & Hospitality (AREA)
• Physics & Mathematics (AREA)
• General Business, Economics & Management (AREA)
• General Physics & Mathematics (AREA)
• Theoretical Computer Science (AREA)
• Management, Administration, Business Operations System, And Electronic Commerce (AREA)
• Image Processing (AREA)

Applications Claiming Priority (3)

Publications (2)

Family

ID=32850786

Family Applications (1)

Country Status (3)

Family Cites Families (1)

• 2004
  • 2004-01-29 EP EP04706413A patent/EP1588239A4/de not_active Withdrawn
  • 2004-01-29 WO PCT/US2004/000215 patent/WO2004070532A2/en not_active Ceased
  • 2004-01-29 CA CA002512946A patent/CA2512946A1/en not_active Abandoned

Non-Patent Citations (2)

Also Published As

Similar Documents

Legal Events