TRAINING APPARATUS

Abstract: Training apparatus for use in contamination surveyance uses a mathematical model of a hypothetical contamination source (e.g. nuclear, bacteriological or chemical explosion or leak) to determine from input data defining the contamination source, the contamination level at any location within a defined exercise area. The contamination level to be displayed by the apparatus is corrected to real time from a real time clock or may be displayed in response to a time input from a keyboard. In a preferred embodiment the location is defined by entering UTM grid reference co-ordinates from a keyboard. The mathematical model used by a microprocessor of the apparatus for simulation of contamination levels in the event of a nuclear explosion is described.

Description: TRAINING APPARATUS The present invention relates to training apparatus and in particular to such apparatus for training personnel in surveying for nuclear, bacteriological or chemical contamination. There is a requirement for personnel to be trained in surveying large areas of terrain for nuclear fallout and/or contamination which m^y be airborne and/or at ground level. Deliberately contaminating a large aj?ea for training purposes is not a practical proposition. Accordingly, methods of training have been in the past purely hypo¬ thetical or have used low ' level radio active sources con¬ cealed within the field of the exercise. One more recent nuclear. radiation simulator uses a frequency modulated radio signal transmitting an encoded digital data signal representing the contamination level. The coded signal was, decoded by replica radiation meters all of which, had the same value "measured" regardless of their respective geographical locations. It is an object of the present invention to provide training apparatus for use in more realistic simulation exercises involving nuclear, bacteriological or chemical conta ination. The term "contamination" as used hereinafter means nuclear and/or bacteriological and/or chenical contamin¬ ation. According to the present invention in training aj^araus ϊov use in contamination surveillance training, a data store is arranged to store data derived from a mathematical model of an assumed source of contami- nation, said data defining contamination levels at a plurality of locations at respective elapsed times after the assumed commencement of contamination and interpolation means is responsive to input signals to said apparatus which define the present location of the apparatus and to signals which define the present time to derive from the data held in the data store a value indicative of the assumed present contamination level. Preferably the input signals which define the location of the apparatus are from a manually operable keyboard and the apparatus may include a real time clock arranged to supply the signal which define the present time. The interpolation means may be a microprocessor which may be programmed to derive the data to be stored from input parameter of an assumed contamination source. One embodiment of training apparatus in accordance with the invention will now be described by way of example only with reference to the accompanying drawings of which:- Figure 1 is a block schematic diagram of the apparatus, Figure 2 shows a keyboard layout of the apparatus, Figures 3* ^ and 6 are fallout plans used in the description of an alga ithm used by the -apparatus, and Figure is a schematic map grid pattern. Referring to Figure 1 the apparatus comprises a microprocessor 1 having an associated program store 2 and data store 3- Ωie program store 2 stores a program which enables the microprocessor 1 to derive (from data held in the data store 3) an assumed contamination level for dis¬ play on .a display device 4- which may be a meter but is preferably an alpha numeric display. i'o enable a user of the apparatus to enter a loca- tion a keyboard 5 is also provided. The keyboard 5 is also used by an instructor to enter data defining the cause of contamination and the time of its occurrence in the data store 3» The apparatus, whilst being capable of responding to a time input from the keyboard 5 also includes a real OM time clock 6 which is set before commencement of an exercise by the instructor. So that a reπuired number of simulation meters may all receive the same data, such that in use all trainees receive compatible contamination level readings from their respective meters, the apparatus includes a serial communications link 7 which enables the simula¬ tion meters to- be preprogrammed by entering data in one of them and transferring it to the rest by use of the serial communications link 7« In order that the operation of the apparatus may be readily understood an example of the actions performed during an exercise simulating the effects of a nuclear explosion will nov be given. The mathematical model of a nuclear explosion used by the microprocessor 1 as a low altitude ai_ι?'blast or as a groundlevel blast is described hereinafter. The .calcu¬ lations require the following integers namely: the yield of the assumed bomb, the location of the blast and time of occurrence, the exercise area and the wind speed and ' direction at varying intervals of altitude and time. Addi¬ tionally it is necessary for the real time clock 6 to be set to actual time. Accordingly the microprocessor 1 is programmed to guide an instructor through entry of the required data. Thus referring also to Figure 2 the user, having switched on the power and causing a built in test function (des¬ cribed hereinafter) to be carried out if required depresses the "ΞXΞB SET UP'" key of the keyboard 5. The micropro- cesεor now causes the display 4- to display YIΞD? The user now keys in the kilotonnage of the fission products of the hypothetical burst using the keyboard* 5, the value of which is displayed on the display 4- until the enter/ clear key of the keyboard 4 is depressed. (It is here noted that the mathematical model described hereinafter is capable of use for training in surveyance for contami¬ nation caused by small tactical nuclear weapons ranging up to large strategic nuclear weapons.) The microprocessor 1 now causes "G.Z.GHID" to be displayed which invites the user to enter the map grid reference of the location of the hypothetical burst (the 'ground zero" reference). The G.Z reference is entered by depressing the enter/clear key again at which time the display 4- is caused to display "3. D/TIK?". The planned date (day only e.g. 5 for the 25th) and time in the twenty four hour clock format (e.g. 0900 for nine a.m. and 2100 for nine p.m.) of the hypothetical burst and a letter reference for the time zone (e.g. 'z* for Greenwich Hea Time) are now keyed in. After checking, the operator again depresses the enter/clear key. The display 4 now shows "GHID S.jU-R S?". This is a request for the exercise area to be defined as an area of 100 kilometre (EH) squares to a maximum of ten squares in each direction. These are keyed in in apgrid character pairs commencing from the south-west corner of the exercise area and ensuring that the defined area includes the "ground zero" location previously referred to. Thus keying in "AA AB AC enter" defines an exercise area extending east of the south-w τ est corner by 300 EK. Further sets of pairs of digits are entered to define 100 EK exercise areas extending north of the south-east corner of the exercise area; (See e.g. Figure 5) Thus keying "AA AB AC enter BA BB BC enter CA CB CC enter" defines a 300 EK square sided exercise area. The microprocessor 1 having received the definition of the exercise area causes the display 4- to display "wTIvD LΞVEL H?". The mathematical model assumes that the wind speed and direction is constant throughout the exercise area. However it is arranged to handle wind speeds and directions which vary with height in 2000 metre bands up to 30,000 metres. "Wind level K" is a request for the mean wind speed and direction in the 0-2000 metres above ground level band. This is entered as a six digit number, the first three digits representing the wind direction in degrees from Eorth and the latter three digits being the wind speed in kilometres per hour. Thus keying in "180020 enter" would indicate a twenty ' kilometre per hour wind from the south in the 0-2000 metre band. When the first band wind level has been entered the microprocessor 1 requests the wind for bands I (2000 to 4000 metre) to W (28,000 to 30,000 metres) in turn. Having entered all of the data it is necessary for the real time clock 6 to be set. The microprocessor 1 causes the display 4- to show "SYNC D/TIH?" at which the instructor enters the current date and time in the same format as for entering the planned date and time of the hypothetical burst. The real time clock 6 commences running from the time at which the enter/clear key is depressed. Thus for an exercise based on a fifty kiloton nuclear burst occurring close to the south-west corner of a 300 EH x 300 K exercise area at 0900 Hrs on the 25th day of the month wit the simulator being programmed at 1800 Hrs on the previous day the actions are as follows:- DISPLAY ACTIONS Press ΞXΞR. SET button YIELD? Eey in yield YIELD 50 ENTER/CLEAR G.Z. GRID? Eey in ground zero map grid reference, G.Z. AA 100200 ENTER. OMPI DISPLAY ACTICNS B. D/TIM? Ee y in date/time of burst , B.D/TIM25 O900Z ENTER GRID SQUARES? Eey in letter pairs for first three lOOE squares, AA AB AC Eey in letter pairs for second three 0 lOOEm squares, BA BB BC ENTER Eey in letter pairs for third lOOEm squares, CA CB CC ENTER s5 WIND LEVEL H? Eey in wind direction/speed for level H, H 160020 ENTER Eey in wind direction/speed for level 0 I, Eey in wind direction/speed for level J 165025 ENTER 5 Eey in wind direction/speed for level --» E 165025 L ...... 0 up to level * u» SYNC D/TIM? Eey in ^present date/time to next minute D/TIM 24- 1800Z When operators watch reaches the entere time press ENTER 5 FXΞR. SET ENDS The "ΞXER SET ENDS" message is displayed for five seconds after which the simulator is ready for use. Pressing the "Ξ2ΞR SET UP" key again enables the operator to check the data which has been entered simply by depressing the enter/clear key in respect cf each line. The microprocessor 1 causes the previously entered data to be displayed by the display 4- which may be amended by the operator if necessary by keying in the correct data. The mathematical model of a nuclear explosion occurring at a burst time ("H") which is held in program store 2 operates in the microprocessor 1 in the following manner. A general cloud model composed of discrete elements spaced vertically in the "mushroom cloud" of the burst is assumed. Each element is assumed to fall from its starting height through the various wind levels until it reaches the ground. Each element is assumed to have a given vertical velocity such that .if the wind speed were twenty-four kilometres per hour the element is assumed to have fallen at the correct down-wind distance from ground zero. Variations of the wind speed from twenty- four kilometres/hour entered are taken into account by generating a synthetic wind speed as hereinafter described. All simulations and calculations are based on the Universal Transverse Mercator (UTM) Map Grid system (although it will be appreciated that any Map Grid system could be employed) and calculations are carried out to an accuracy of 100 metres (1/10 EM). In the map grid system two digits designate a unique grid square of 100 EM by 100 EK. The following three digits are eastings and the last three digits are northings Thus referring to Figure a grid reference of AA 000000 is the south-west corner of the exercise area previously specified, whilst the ground zero location of the burst AAloo 200 is ten kilometres east and twenty kilometres north of the south-west corner of the exercise area. The microprocessor 1 is programmed to calculate a synthetic wind speed for each 2000 metre band. As has been previously mentioned the fallout vertical velocity algarithm has been corrected for a windspeed of twenty-four EM/Er. Wind speeds of less than twenty-four EK/Hr give a synthetic wind speed of twelve plus (0.5 x entered wind speed)whilst . wind speeds in excess of twenty-four EK/Hr give a synthetic wind speed of eighteen plus (0.25 x entered wind speed). The wind direction and strength is resolved into Easterly and Northerly components, with these components having positive values. It will be appreciated that winds having a westerly component will give a negative easterly value and similarly winds having a southerly component will have a negative northerly value. The easterly and northerly components are derived by multi¬ plying synthetic wind strength by the sine ' and cosine respectively of the wind direction. Referring now to Figure 3 the fall of each of the forty elements of the cloud 300 are considered. he first of the forty elements (referred to as element C) is considered to be at two-thirds of the stem height of the cloud as indicated by the chain dash line 301 whilst the last of the elements (element 39) is at the top of the cloud arising from the burst. Each element is given a predetermined rate of fall with element 'C having the highest rate and element 39 the lowest rate. As each element falls toward ground level its direc- • tion and lateral displacement is processed by application of the wind components in respect of each of the bands through which it falls. The path 302 of the forty elements on reaching ground level is thus-determined. It is here noted that the path 302 may be curved. As the data is calculated for each of the forty elements a ground position and time of arrival is deter¬ mined together with the peak intensity of radiation at -9- that point and the'width of a line extending at right angles either side of the path 302 on which the radiation calculation (appropriately scaled from the peak intensity) for the particular element is effective. 5 Consider element 20 and element 24- of the cloud 300, each has a respective path of fall 120.124- and on arrival at ground level gives a respective lateral width slice 220, 224- having a peak intensity immediately above the line of fall which reduces in a sinusoidal manner 10 towards the ends of the slices as indicated by intensity profiles 320,324. The vertical velocity of each element of the cloud is used to calculate the time of arrival of each element at ground level and each element is allocated a duration 15 of arrival, that is a period of time over which the radiation attributable to the slice will arrive and which is calculated as ' being 0.75 of the period between the burst time H .and the time of arrival. It is here noted that the cloud is assumed to 20 stabilise at H + ten minutes at a height dependent on the kilotonnage (W) of the bomb. The peak intensity of each element is stored as the intensity at H + one hour ("H + 1" time). Thus the data calculated and stored in respect of 25 the numbered elements is given in the following table: Start Vertical Peak Max Time Dura¬ Position No, Height Velocity Intensity Width Of tion of (KM) KM/HR Rad/Hr KM Arri¬ of Arrival I val Arri¬ CO val Hrs (D) 5 Having calculated and stored data in respect of all forty elements, which calculations may be carried out by the microprocessor prior to commencement of the exercise, the microprocessor 1 is ready to calculate the intensity (rads/Hr) at any map grid position within the exercise area at any time. The exercise may be carried out in either one of two modes. The principal mode of operation relies on the real time clock 6, the user of the appa- ratus entering only a map grid position (of his loca¬ tion) and, if necessary a shielding factor by way of the keyboard 5- The shielding factor is a measure of the protection given by material surrounding the user, typical factors being ten to twenty for tanks, two to four' for armoured personnel carrier and one to two for other skinned vehicles. If no shielding factor is entered then the sheilding factor is assumed to be one. The microprocessor is programmed to cal- culate the radiation (rads/hr) and to divide the figure by the shielding factor, (if any). The calculated figure is displayed on the display 4- except when the 'CUM. OSE' key of the keyboard 5 is depressed when the accumulated (time integrals) dose received by the user of the meter since the burst time (H) is displayed. It is here noted that dosages in excess of 3,000 rads/hr are displayed solely as "over 3000" although the microprocessor 1 is programmed to handle higher doseage rates and could be arranged to display the higher (actual) figures. Although the orientation of the slices for which calculations have been made will vary on a curved fallout path, all calculations are based on the orien- tation of the central slice (that resulting from the twentieth element of the cloud). When a map grid re erence is entered from the keyboard 5 the microprocessor 1 generates the para¬ meters of an artificial slice which passes through the point specified (the survey point) from the previously calculated and stored parameters of the ones of the forty slices which encompass the survey point. For this purpose it is necessary to identify the pairs of slices which encompass the survey point, Thus, referring to Figure 6 and assuming a slice orientation ofjclcgrees a calculation is carried out in respect of each of the forty slices and the arti¬ ficial slice to project same to a point on a base line having a north axis. (However if 0 is close to zero the projection may be done to a base line having an east axis.) Since all slice orientations are the same the equation of the extended slice intercept Cx on the north axis is given by Cx = N - Cot 0 Ex. Cx is evaluated for each of the slices giving C^ to C^ Q for the orty original elements of the cloud and C for the survey point E N . A check is then carried out to determine between which (if any) pairs of slices the survey point lies. • OM It is here noted that if G is on the opposed side of C Q to C, it is necessary to consider a three di en- sional extrapolation of slice zero as shown in Figures 4A and 4-B. However as shown in Figure 5 the survey point ΣsUs lies betveen slice two and slice three for which parameters are stored in the data store 3« Having determined that the survey point lies between slices two and three it is necessary to use the parameters of these two slices to generate the parameters of the artificial slice. Thus the proportional position of the point E' E' is evaluated as follows:- The• ec j uation of the line j"oining^ Ξ2N2 to Ξ y-,E_ is given by 2 = _2 Therefore N • = The line through the survey point is given by N = C o + E cot Solving these equations simultaneously gives the co-ordinates of the crossing point E' N' , The longitudinal position ratio of E g Ε % on the line E N,-, to E-, N, is now determined from to give a ratio figure E of the distance between the slices two and three at which the artificial slice lies. The width(S) of the artificial slice may now be found from the width S of slice two and the width (S,) of slice three and is given by S = S 1 + H (S 2 - s . The microprocessor 1 determines the distance (6.) of E I< from E' » from I = J(E g - E ≤ + (K g - K g ) 2 and determines the ratio £- S Since the width of the artificial slice is assumed to be symetrical about E s' . N s 1 - it will be ap ~ preciated that if ϊ is greater than 0.5 the survey point is outside the S contour of the assumed radiation and has a zero radiation intensity. In such circumstances, further consideration by the microprocessor is not required. If the survey point is within the contamination area then the final peak intensity (Is) is calculated from the peak intensity of the two adjacent slices. Thus I- g = I 2 + E (I- - I 2 ). The time of arrival is similarly found from Ts = 2 H- S(T, - τ 2 ) and the duration of accumulation may be found either from 0.75 x Ts or from Ds = D 2 + 3 D^ - D 2 ). Tsis now compared with the time in the real time cloc 6 to determine whether irradiation has commenced. If not then again zero is displayed on the display 4. Provided that the real time is later than Ts, the microprocessor 1 is programmed to determine the final intensity (^) of the radiation at the survey point from the finalpe Thus .the following parameters of the artificial slice have.been found:- The time of interrogation (t) (from the real time clock), The final intensity (is) at the survey point, The time of arrival (Ts) of the initial fallout, The duration of arrival (Ds) of the fallout and the shielding factor (if any) from the keyboard. To determine the radiation intensity to be displayed the microprocessor examines t_ to determine whether it is after Ts + Ds. If so further calculation is based on is. If t_ lies between Ts and Ts + Ds then an accumulation rati factor t - Ts is used to determine from i g the intensity value on w ich the final calculation is based. It will be realised that the accumulation ratio factor will change from zero at Ts to 1.0 at Ts + Ds. It will be remembered that the intensity values have been normalised for radioactivity back to one hour after the burst time ("H + 1")- Since the radiation intensit decays with time, is must be multiplied by a decay factor. The microprocessor 1 therefore decays i_ according to t —1 * 2 (assuming t is after H + 1). The figure is displayed in rads/hr on the display. Periodically, say, every two minutes, the microprocessor 1 recalculates the radiation intensity taking account of time variations and any newly entered grid re erences and causes a new value to be displayed such that over a period of time the display 4 reflects the reduction (or increase) in intensity which would occur in reality. It is noted that if a shielding factor is entered it is also taken into account prior to display. The program is also arranged to integrate the doseage in real time so that the cumulative dose may be displayed on demand. Referring, fco Figure 4B if the calculation of Cs as hereinbefore described indicates that the survey point is on the opposed side of Co to CL , such as is shown by the survey point B it is necessary to determine only the distance of B from Eo No and to apply this as '_£' on slice 'C in the algorithm described. In the immediate vicinity of ground zero, neutron induced gar..sa activity (NIA) must also be t-'aken into account by the program. The effect of NIA beyond a distance of, say, five kilometres from ground zero is ignored. Otherwise the distance (S) of the survey point from ground zero is calculated and the H + 1 intensity (I) attributed to NIA is calculated from log, Λ I « - 2.2564S + 0.6021. 10 ^TCT A look up table is provided in the store 2 from which a multiplication factor for NIA dependent on the time after 'E 1 is recovered by the program and applied to the NIA "H + 1" intensity. Although the factor is dependent on soil composition a typical table is t (Hours after 0.1 0.2 0.3 0.4 0.5 0.6 0.7 l-O burst) Factor 30 2.5 1-3 1.2 1.1 1.051-0 t (Hours after 2 5 10 20 50 burst) Factor .85 -65 .42 .23 .06 The NIA intensity is added to the fallout radiation intensity before display. For realism the algorithm may be arranged to intro¬ duce pseudo-random 'hot spots' " of radiation intensity within the exercise area. If the fallout path to be assumed is exceedingly curved it is possible for the survey point to be encom¬ passed by more than one pair of slices. In this case the calculations described above will be carried out in respect of each cf the pairs and their radiation intensity effect at the survey point will be added together prior to display of the intensity figure. Referring again to Figure 1, in order to facilitate programming of a large number of simulated survey meters, the apparatus is provided with a serial communications link 7- number of meters may be linked together in serial fashion and the input data having been entered in a first one may be copied under control of the micro¬ processor 1 into each of the others. For classroom use the second mode of operation which is provided is one in which the time of survey is entered from the keyboard 5 rather than being derived from the real time clock 6. Thus an exercise which could last several days may be carried out in a few hours if desired. It is emphasised that all meters set up for a particular exercise are arranged to give consistent result It will be realised that the calculation of the para meters of the forty slices could be carried out in a centr computer and the data transferred to the apparatus thus reducing the program required by the microprocessor 1. In an alternative embodiment each of the survey apparatuses' incorporates a radio receiver and the results of the contamination calculations and updates are broad- cast in a serial data stream using a standard FM digital transmitter including radio relay stations if necessary. In this case the microprocessor selects from the digital data stream the data which is correct for display at the keyed in grid reference. It will be appreciated that additional bursts and changes in wind speed and direction may be added during the course of an exercise. In this case the effect of each burst on a particular survey point is determined and the values added together for display. It will he appreciated that the map grid positions could also be determined by use of multiple radio beacons in the exercise area the signals from which could be used in the apparatus to derive the location of the apparatus. In an alternative mode of operation position signals may be received from an aircraft navigation system. For classroom use, such as in training control centre personnel in plotting and prediction techniques, the time (t) may be entered from the keyboard 5 ε∑td the map grid position entered as before. The microprocessor 1 carries out the calculation in the manner hereinbefore described but substitutes the entered time for the real time throughout the algorithm. The following commercially available devices have been found suitable and may be used for implementation of the invention:- The microprocessor 1; Intel 80C51/31, the data store 3 (RAM); "Hitachi" 62641P, the program store 2 (Ultra Violet Eraseable PROM); MBM 27C64 and the real ime clock: National Semiconductor MM 58174. The keyboard 5 may be implemented by use of two four by four keyboards suitably mounted one such commer¬ cially available keyboard being manufactured by Sadar (of France) under part no. OCX 44MCN3. The display 4 may be of the kind known as a "Liquid Crystal Diode" display one example of which is available from PCI (Europe) Limited of England under part number PC1M200-M.

Claims: Clai s 1. Training apparatus for use in contamination surveillance training wherein a data store is arranged to store data derived from a mathematical model of an assumed source of contamination, said data defining contamination levels at a plurality of locations at respective elapsed times after the assumed commencement of contamination and interpolation means is responsive to input signals to said apparatus which define the present location of the apparatus and to signals which define the present time to derive from the data held in the data store a value indicative of the assumed present contamination level. 2. Training apparatus as claimed in Claim 1 wherein said input signals which define the location of the apparatus are from a manually operable keyboard. 3. Training apparatus as claimed in Claim 1 or Claim 2 including a real time clock arranged to supply the signals which define the present time. 4. Training apparatus as claimed in Claim 2 in which said signals which define the present time are input from the manually operable keyboard. 5- Training apparatus as claimed in Claim 1, Claim 2 or Claim 3 in which said interpolation means is a microprocessor. 6. Training apparatus as claimed in Claim 5 in which the stored data is derived from input parameters of an assumed contamination source. 7- Training apparatus as claimed in Claim 6 in which the input parameters are used to derive a respective , plurality of parameters for each of a plurality of linear slices of an exercise area and corresponding parameters for a specified survey point are derived by interpolation of the respective parameters of at least one pair of said plurality of linear slices which encompass said survey point. OMPI 8. Training apparatus as claimed in Claim 7 in which said plurality of parameters in respect of each of said linear slices define the centre of location of arrival and width of the slice, the starting time of contamin- ation at the location, the contamination intensity of the centre of the slice extrapolated to a particular time and the period of time after said starting time for the contamination to reach the contamination inten¬ sity. 9- . Training apparatus according to Claim 8 in which the contamination intensity at a survey point is inter¬ polated from the contamination intensity parameter by application of a time related decay factor. 10. Training apparatus according to Claim 8 or Claim 9 in which the contamination intensity at a survey point is interpolated from the contamination intensity para¬ meter in dependence upon the distance of the survey point from the centre of location parameter. 11. Training apparatus as claimed in any preceding claim including display means arranged to display said value indicative of the assumed present contamination level. 12. Training apparatus as claimed in any preceding claim in which said interpolation means is arranged to accumulate a value indicative of the total contamination doseage assumed to have been received by the user of the apparatus over a period of time. 13- Training apparatus as claimed in any preceding claim incorporated in a portable simulation meter. 14. Training apparatus as claimed in Claim 12 Including input/output means such that a plurality of such apparatus • may be interconnected to permit the transfer of data between them. O

Inventors: MONTEITH WILFRED DAWSON, MONTEITH, WILFRED, DAWSON
Applicants: MARCONI INSTRUMENTS LTD, MARCONI INSTRUMENTS LIMITED
Application year: 1983
Main classification: G09B9
IPC codes: G09B9