Test Method of Trap Performance for Induced Siphonage

Similar documents
GEA FOR ADVANCED STRUCTURAL DYNAMIC ANALYSIS

The Estimation Of Compressor Performance Using A Theoretical Analysis Of The Gas Flow Through the Muffler Combined With Valve Motion

Critical Gust Pressures on Tall Building Frames-Review of Codal Provisions

DIAGNOSTICS OF IMPULSE LINE BLOCKAGE WITH A MULTI-SENSING DIFFERENTIAL PRESSURE TRANSMITTER AT THE AIR LINE

CFD Simulation and Experimental Validation of a Diaphragm Pressure Wave Generator

Sizing Pulsation Dampeners Is Critical to Effectiveness

Acoustical Modeling of Reciprocating Compressors With Stepless Valve Unloaders

Experimental determination of deflagration explosion characteristics of methane air mixture and their verification by advanced numerical simulation

Quantification of the Effects of Turbulence in Wind on the Flutter Stability of Suspension Bridges

Analysis of Pressure Rise During Internal Arc Faults in Switchgear

The University of Hong Kong Department of Physics Experimental Physics Laboratory

Noise Characteristics Of A Check Valve Installed In R22 And R410A Scroll Compressors

Special edition paper

CHARACTERISTICS OF LEAK DETECTION BASED ON DIFERENTIAL PRESSURE MEASUREMENT

WATER HYDRAULIC HIGH SPEED SOLENOID VALVE AND ITS APPLICATION

Citation Journal of Thermal Science, 18(4),

The Study of Bubbly Gas-Water Flow Using

The Use of a Process Simulator to Model Aeration Control Valve Position and System Pressure

Modal Analysis of Propulsion Shafting of a 48,000 tons Bulk Carrier

Gerald D. Anderson. Education Technical Specialist

sensors ISSN

AN ANALYSIS ON HIGH PRESSURE DYNAMIC CALIBRATORS USED IN THE DEFENSE AREAS

Vibration of floors and footfall analysis

SEMATECH Provisional Test Method for Pressure Cycle Testing Filter Cartridges Used in UPW Distribution Systems

Waves & Interference

Introduction to Pumps

2. Determine how the mass transfer rate is affected by gas flow rate and liquid flow rate.

Analysis of pressure losses in the diffuser of a control valve

Currents measurements in the coast of Montevideo, Uruguay

Conditions for occurrence of vortex shedding on a large cable stayed bridge. Full scale data from monitoring system

PROPAGATION OF LONG-PERIOD WAVES INTO AN ESTUARY THROUGH A NARROW INLET

AIRMOUNT VIBRATION ISOLATION

WHITEPAPER. Integrated process monitoring of dosing pressures

Structural Design and Analysis of the New Mobile Refuge Chamber

Yasuyuki Hirose 1. Abstract

Conventional Ship Testing

A review of best practices for Selection, Installation, Operation and Maintenance of Gas meters for Flare Applications used for Managing facility

Special edition paper

Development of a Non-Destructive Test Machine for Testing Tunnel Lining Concrete

Aeroelasticity of Turbomachines

Application of Computational Fluid Dynamics to Compressor Efficiency Improvement

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS AP PHYSICS

Physics Extended Essay. A research about the effect of sound waves on standing waves by using Ruben s Tube.

Components for air preparation and pressure adjustment. OUT port position ( ) connected Rear side. of IN port. Air tank. directly.

V10K GAS FEED SYSTEM WALLACE & TIERNAN PRODUCTS

ADVANCES in NATURAL and APPLIED SCIENCES

Exercise 4-2. Centrifugal Pumps EXERCISE OBJECTIVE DISCUSSION OUTLINE DISCUSSION. Pumps

Friction properties of the face of a hand-held tennis racket

TESTING AND ANALYZING ON P-V DIAGRAM OF CO 2 ROLLING PISTON EXPANDER

Offshore platforms survivability to underwater explosions: part I

Dynamic Component of Ship s Heeling Moment due to Sloshing vs. IMO IS-Code Recommendations

AC : MEASUREMENT OF HYDROGEN IN HELIUM FLOW

Inlet Influence on the Pressure and Temperature Distortion Entering the Compressor of an Air Vehicle

VACUVENT Valve functions - benefits - flowpath - VG series - Sewage

Causes and Remedies to 7 Instability Patterns in Blown Film Extrusion

Research of Variable Volume and Gas Injection DC Inverter Air Conditioning Compressor

Gas Pressure Regulator HON 200

Effect of Nozzle Twisted Lance on Jet Behavior and Spitting Rate in Top Blown Process

Pilot HON 625. Entwurf. Product information. serving the gas industry worldwide

INSTRUMENT INSTRUMENTAL ERROR (of full scale) INSTRUMENTAL RESOLUTION. Tutorial simulation. Tutorial simulation

3. REQUIRED EQUIPMENT

Influence of rounding corners on unsteady flow and heat transfer around a square cylinder

The Discussion of this exercise covers the following points: Pumps Basic operation of a liquid pump Types of liquid pumps The centrifugal pump.

FORECASTING OF ROLLING MOTION OF SMALL FISHING VESSELS UNDER FISHING OPERATION APPLYING A NON-DETERMINISTIC METHOD

SOME TECHNICAL AND PRACTICAL RECOMMENDATIONS ABOUT PULSATION DAMPENERS IN CIRCUITS WITH DOSING OR VOLUMETRIC PUMPS

Experimental Investigation on Transient Response of Two Stage Pressure Relief Valve

Exercise 2-2. Second-Order Interacting Processes EXERCISE OBJECTIVE DISCUSSION OUTLINE. The actual setup DISCUSSION

DEVELOPMENTS IN PIPELINE FILLING AND EMPTYING EXPERIMENTATION IN A LABORATORY PIPELINE APPARATUS

New Findings for the Application of Systems for Explosion- Isolation with Explosion Venting

TEK-SUB 4800B. Submersible Level Transmitter. Flow Level Temperature Pressure Valves Analyzers Accessories TekValSys LEVEL

VORTEX SHEDDING AND VORTEX FORMATION FROM A PAIR OF IN-LINE FORCED OSCILLATING TANDEM ARRANGED CIRCULAR CYINDERS IN A UNIFORM FLOW

This educational seminar discusses creating, measuring, and troubleshooting Rough Vacuum.

I.CHEM.E. SYMPOSIUM SERIES NO. 97 BUOYANCY-DRIVEN NATURAL VENTILATION OP ENCLOSED SPACES

Study on Fire Plume in Large Spaces Using Ground Heating

CHEMICAL ENGINEERING LABORATORY CHEG 239W. Control of a Steam-Heated Mixing Tank with a Pneumatic Process Controller

FLOW CONSIDERATIONS IN INDUSTRIAL SILENCER DESIGN

RESEARCH OF BLOCKAGE SEGMENT DETECTION IN WATER SUPPLY PIPELINE BASED ON FLUID TRANSIENT ANALYSIS ABSTRACT

CORESTA RECOMMENDED METHOD N 6

Ratio Regulator Zero Governor. FRG/6 Series

An underwater explosion is an explosion where the point of detonation is below the surface of the water.

Prof AH Basson, Pr Eng March Departement Meganiese en Megatroniese Ingenieurswese. Department of Mechanical and Mechatronic Engineering

Available online at Procedia Engineering 00 2 (2010) (2009)

International Journal of Technical Research and Applications e-issn: , Volume 4, Issue 3 (May-June, 2016), PP.

A MODEL FOR ANALYSIS OF THE IMPACT BETWEEN A TENNIS RACKET AND A BALL

Acoustical approach to analysis of energy conversions in an oscillating bubble

WHITEPAPER. The flowplus 16 Pressure Sensor. Physical Principles and Function

Cover Page for Lab Report Group Portion. Pump Performance

Constant pressure air charging cascade control system based on fuzzy PID controller. Jun Zhao & Zheng Zhang

CFD Analysis and Experiment Study of the Rotary Two-Stage Inverter Compressor with Vapor Injection

Vibration Analysis and Test of Backup Roll in Temper Mill

Pressure relief valve DHV 712 DN 65-80: 0,5-10 bar, DN : 0,3-4 bar, DN 100: 0,5-6 bar

13. Evaporative Emission Control System

Experimental investigation of the check valve behaviour when the flow is reversing

Applied Fluid Mechanics

Application of pushover analysis in estimating seismic demands for large-span spatial structure

Application of Simulation Technology to Mitsubishi Air Lubrication System

CanNeed-DSST-100 Digital Secure Seal Tester --- To evaluate the secure seal quality of PET bottles, glass bottles and cans

LABORATORY EXERCISE 1 CONTROL VALVE CHARACTERISTICS

Standing Waves in a String

A Research on the Airflow Efficiency Analysis according to the Variation of the Geometry Tolerance of the Sirocco Fan Cut-off for Air Purifier

Transcription:

VII.4 Test Method of Trap Performance for Induced Siphonage K. Sakaue (1), H. Kuriyama (2), H. Iizuka (3), M. Kamata (4) (1) sakaue@ isc.meiji.ac.jp School of Science and Technology, Meiji University, Japan (2)kuriyama.hana@nikken.co.jp Nikken Sekkei Ltd., Japan (3) iizukah@nikken.co.jp Nikken Sekkei Ltd., Japan (4) ft11743@kanagawa-u.ac.jp Faculty of Engineering, Kanagawa University, Japan Abstract The drainage system is designed in such a way that drain traps are not prone to seal break. Seal break is caused by seal loss phenomena, and induced siphonage plays an important role. A reliable drainage system cannot be designed without first defining trap performance for induced siphonage. However, no effective methods have been established to test the performance of traps. Based on the premise that seal loss phenomena are produced as response of seal water to pneumatic pressure in drain, we conducted free vibration test, single sine wave response test, and 2 times natural frequency response test using a pressure generating device with two types each of fixture traps, floor drain traps, and WC to elucidate their characteristics and strengths. The results clearly indicated that methods defining minimum pressure required to cause seal break was most appropriate for testing traps performance to withstand pneumatic pressure. It was also found that natural 47

frequencies obtained by the free vibration test, and max. response magnifications by the single sine wave response test were relevant and applicable. Keywords Drainage system, trap, induced siphonage, test method, single sine wave 1. Introduction The first consideration in designing any drainage system including vent system is to make it in such a way that it does not cause seal break. So-called seal break phenomena include such conditions as induced siphonage, self-siphonage, and evaporation. Of these, however, one that is pertinent to the entire system would be induced siphonage, which is characterized by seal loss caused by pneumatic pressure (referred to as pressure below) affecting the trap seal. In other words, induced siphonage is a dynamic vibrational response phenomenon caused by the flow of water in the drainage pipe. To prevent seal break caused by this induced siphonage, it is necessary to establish a way to predict pressure resulting from discharge load and entire system configuration such as piping, pipe diameter, type of vent system, as well as to clarify anti-pressure performance of trap (referred to as trap performance below). Theoretically it can be achieved by establishing permissible pressure based on trap performance, and by defining permissible flow rate against it for each system type. The more clearly we define the performance of each trap, the more flexibly we can design the section of piping in which the trap is placed. In the current designing scheme, trap performance is only defined by minimum seal depth with no pressure prediction made, and at the very best rough permissible flow rate for each system type is defined. Of all plumbing codes in the world, only the former DIN1986 in Germany and SHASE-S-218 in Japan stipulate permissible seal loss. As for permissible residual seal depth, only SHASE-S-218 has stipulations. In the ever-evolving world of plumbing system, designing of drainage system, extremely ill-defined compared to that of water supply system and hot water supply system, has been unable to respond to the trend of rationalization. Though the selfsealing traps without requiring seal water have been developed in Great Britain, water seal trap is still the mainstream in the world. It seems imperative that we clearly define trap performance to improve designing methods of drainage system. Based on the above considerations, the authors have tested various methods of defining trap performance [1-4]. In this study, as continuation of the ongoing research, the authors have examined testing methods using a simple device generating single sine waves with two types each of fixture traps, floor drain traps, and WC. 48

2. Experimental apparatus and measuring equipment 2.1 Experimental apparatus The experimental apparatus used in this study consists of a pressure generator (pistons, a frequency variable device, various chambers), drainage pipes, and an analyzer (Figure 1). It has a triple-piston structure with cranks adjustable at 8 steps between 15 ~ 5mm in increments of 5mm. Frequencies are variable in increments of.33hz within the range of.166 ~ 4.5Hz. It is also equipped with a small blower as a bias pressure device, which reproduces the steady pressure component. Amplitudes were adjusted by changing the water contents of vertical cylindrical chamber and pressure adjustment chambers, which in turn changed air volume in the chambers. Only one of the three pistons on the apparatus was operated this experiment. 2.2 Test traps Basic configuration parameters and cross sectional shapes of test traps are shown in Table 1 and Figure 2. The and for fixtures, the bell trap and contrary bell trap for floor drains are made of transparent resin, and seal conditions inside can be observed from outside (Figure 2). Actual models of two types of WC, supper water saving WC (WC -6L ) and water saving WC (WC -8L ) were used. 49

2.3 Measuring equipment We used a diffuse semiconductor type pressure sensor for measuring pressure, and a capacitive water level sensor for seal water level. The water level sensor had been modified so that it could be inserted into the inlet leg of a trap. Both of them had adequate response frequency of more than about 1kHz. 8 9 12 13 1 7 1 11 2 6 3 3 3 4 4 4 5 1; Changing device of frequency 2; Blower for bias pressure device 3; Chamber with piston 4; Pressure control valve 5; Mixing chamber 6; Pressure control chamber 7; Test trap 8; Water level sensor 9; Amplifier 1; Pressure sensor 11; Amplifier 12; Data roger 13; PC Figure 1 - Test apparatus 41

Table 1 Basic parameters of test raps Test traps Volume of seal Seal depth Ratio of leg s sectional area F 1) [-] For Fixtures 15 15 6 6 1. 1. For floor drains Contrary trap bell 55 41 53 5 1.41 1.26 WC WC -6L WC -8L 2,47 1,37 55 67-2) - 2) Note 1) F = (mean sectional area of inlet leg) / (mean sectional area of outlet leg) of a trap 2) These weren't confirmed because an internal form was complex. 88 45 31 14 6 6 53 5 3 4 3 3 4 3 92 112 1 33 43 Contrary bell trap Figure 2 - Test traps 3. Free vibration test 3.1 Purpose Natural frequency and damping ratio are the two main basic response characteristics in forced vibration phenomena. Having damping ratio of.4 ~.6, seal fluctuation is known to cause weak damping oscillation [5, 6]. Considering the effects on trap performance, we used only natural frequencies as a parameter. Though there are several 411

methods of measuring natural frequencies, we adopted a method based on power spectrum of response of water level in this study. 3.2 Experimental method We first caused half-full seal to free-vibrate, and measured free-vibration wave patterns of seal. Then we obtained the density distribution of power spectrum of response of water level by giving FFT treatment to the wave patterns. The frequency at which the power spectrum density became maximum value was determined to be the natural frequency f. 3.3 Results and discussion As shown in Figure 3, power spectrum density distributions of each test trap were as follows: : 1.95Hz, : 1.93Hz, : 2.43Hz, Contrary bell trap: 2.57Hz. These values corresponded with the values obtained from the calculate equation of natural frequency. Although the density distribution of WC generally forms two peaks as that of WC -6L [7], the second predominant frequency of WC -8L was extremely low forming only one peak. If the most predominant frequency is regarded as natural frequency, the peaks of density distribution for WC -6L and WC -8L were 1.61Hz and 1.46Hz respectively. From this, the range of natural frequencies when test traps were half-full was found to be to be 1.4 ~ 2.6Hz. Power spectrum Power spectrum 12 9 6 2 6 3. 2. 4. 6. 12 9 6 2 6 F : 195Hz F 193H 3 7. 2. 4. 6. F : 243Hz. 2. 4. 6. Contrary F : 243Hz. 2. 4. 6. WC 6L F : 161. 2. 4. 6. WC 8L F : 146. 2. 4. 6. Figure 3 - Power spectrum of water level and natural frequency f 412

4. Single sine wave response test 4.1 Purpose This experiment was intended to define trap performance by obtaining response characteristics (maximum response magnification) when frequencies of sine wave pneumatic pressure were fluctuated with its amplitudes kept constant. 4.2 Experimental method Single sine waves mainly consisting of atmospheric pressure with various frequencies and constant amplitude were applied to half-full seal. Then the response magnification curve as a function of frequency was obtained by based on calculation of the static seal level (X) of trap seal and the ratios of response fluctuations (X t ) against single sine wave pressure (response magnification M=X t /X) at each frequency were obtained. 4.3 Results and discussion Response magnification curves for each test trap are shown in Figure 4. Response magnification curves corresponded with density distribution curves in the free vibration test, and frequencies at which response magnification M reached maximum roughly corresponded with natural frequency f. However, with WC -8L, response magnification comparable to natural frequency f was noted at the second predominant frequency (2.4Hz) Vibration characteristics of test trap seal seemed to manifest more evidently in the single sine wave response experiment than in the free vibration experiment. Experimental method is simpler with only single sine waves of constant amplitudes used. As a result an experimental apparatus is not required to generate a wide range of pressure load, and the experimental apparatus itself can be reduced to small size. Response Response 6 5 4 3 M = 4.8 (1.9 Hz) 2 1. 2. 4. 6. 6 5 4 3 M = 5. (2. Hz) 2 1. 2. 4. 6. ifi ti M ifi ti M. 2. 4. 6. Contrary bell trap. 2. 4. 6. WC 6L. 2. 4. 6. WC 8L. 2. 4. 6. Figure 4 - Response magnification M of test raps 413

Vibration characteristics of test trap seal seemed to manifest more evidently in the single sine wave response experiment than in the free vibration experiment. Experimental method is simpler with only single sine waves of constant amplitudes used. As a result an experimental apparatus is not required to generate a wide range of pressure load, and the experimental apparatus itself can be reduced to small size. 5. Seal break test 5.1 Purpose This experiment is intended to define trap performance by measuring amplitudes at the time of instantaneous seal break when frequencies and amplitudes of single sine waves are changed. Instantaneous seal break refers to the conditions when air bubbles pass through top dip of trap. 5.2 Experimental method Single sine waves consisting of amplitudes of atmospheric pressure with a constant frequency were applied to half-full seal, and the minimum amplitude leading to instantaneous seal break within 1 seconds (referred to as seal break pressure) at the given frequency was measured. Then the seal break characteristics curve was obtained as a function of frequency and seal break pressure. 5.3 Results and discussion Seal break characteristics curves (seal break pressure curves) for each test trap are shown in Figure 5. The seal break characteristics curves roughly corresponded with density distribution curves obtained from the free vibration test and response magnification curves from single sine wave response test. The frequency where seal break pressure was minimum also corresponded with the natural frequency f. However, the minimum value of seal break pressure was extremely small in WC. The reason for this seems to be that WC was made of opaque ceramic, and that instantaneous seal break was not detected accurately as the conditions of seal could not be observed directly. Based on dynamic oscillation, seal break test is geared to focus on seal break, the most important element in this study, and seal break pressure provides convenient data to evaluate trap performance. However, as mentioned above, seal break test has a drawback in that instantaneous seal break cannot be detected accurately if traps are made of opaque materials. 414

Seal break pressure [Pa] 6. 2 times natural frequency response test 6.1 Purpose In the field of oscillation technology, frequency 2 times the natural frequency f of experimental subject has been used as eigenvalue of force. This experiment is -2-4 -6-8 -2-4. 2. 4. 6.. 2. 4. 6.. 2. 4. 6. P s min.= 196 Pa. 2. 4. 6.. 2. 4. 6.. 2. 4. 6. P s min.= P s min.= 216 Pa Contrary bell trap P s min.= WC 6L P s min.= 183 Pa WC 8L P s min.= -6 186 Pa -8 P s min : Minimum seal break pressure 24 Pa 152 Pa Figure 5 Seal break curve of test raps intended to define trap performance by studying response levels of seal to single sine wave with 2 times natural frequency f. 6.2 Experimental method Single sine waves mainly consisting of atmospheric pressure with various amplitudes and constant frequency were applied to half-full seal for 3 seconds, and seal fluctuations and seal loss H loss were measured. Pressure waves with natural frequency f and 2 times f were applied. 6.3 Results and discussion The relationship between applied pressure and maximum seal fluctuation for each trap is shown in Figure 6. The input pressures (pressures waves inputted into the pressure generator) were retained at a constant level until the amplitude of applied pressure reached approximately ±15Pa, but the larger the amplitude, the more differences in output pressures (pressures actually measured in the pipe of test apparatus) were actually measured in traps (Table 2, Table 3). This phenomenon seems to have more to do with seal volume than with the size of cross-sectional areas of trap legs. It can be assumed that input pressure was lessened by resonance phenomenon when frequency f is applied, but on the other hand, air inside the apparatus became compressed before 415

applied pressure had any influence on seal fluctuation as seal volume increased when 2 times f frequency was applied. As far as seal loss is concerned, no evaluation of trap performance was possible as seal loss remained zero regardless of amplitude of applied pressure at 2 times f frequency (Figure 7). The experiment using 2 times natural frequency f proved to be unproductive as it required cumbersome procedures and considerable time, and it was difficult to adjust amplitude of applied pressure because of compression of air inside the apparatus. It was also found difficult to obtain clear indices by which performance of each trap can be evaluated and ranked. 12 9 WC 6L Maximum seal water level [mm] 6 3 12 9 6 3-2 -4-6 7-2 -4-6 Pressure [Pa] -2-4 -6 Contrary bell trap -2-4 -6 Pressure [Pa] -5-5 Pressure [Pa] WC 8L -1-1 Figure 6 Relations ship between Pressure and max. water level of test raps 416

Table 2 Input pressure and output pressure (f ) Input pressure [Pa] Test traps 15 2 25 35 4 45 55 65 75 85 Output pressure [Pa] 225 229 287 352 425 184 237 252 359 424 164 26 22 337 371 Contrary tarp bell 169 235 25 336 378 WC -6L 246 262 318 4 48 45 469 545 564 637 WC -8L 15 12 174 262 31 283 355 418 442 456 Table 3 Input pressure and output pressure ( 2 f ) Input pressure [Pa] Test traps 15 2 25 35 4 45 55 65 75 85 Output pressure [Pa] 167 226 288 323 453 214 21 313 367 445 222 255 333 431 497 Contrary tarp bell 179 223 38 48 49 WC -6L 213 238 314 383 479 511 631 765 941 1.57 WC -8L 187 241 319 422 52 523 6 695 816 9 417

Seal loss [mm] 4 35 3 25 2 15 1 5-2 Pressure [Pa] -4 f 2 f -6 Figure 7 An example of output pressure against input pressure (f, 2 f ) Seal loss [mm] 4 35 3 25 2 15 1 5-2 Pressure [Pa] -4 f 2 f -6 7. Conclusion 2 times natural frequency response test was found inappropriate for testing trap performance. Considering the ease of experimental method and the size of experimental apparatus, minimum seal break pressures obtained by the seal break test seems to provide most appropriate criteria for evaluating trap performance. In addition, it is desirable to include natural frequencies and maximum response magnifications obtained by the free vibration test and the single sine wave response test respectively. An issue that needs to be addressed in the future is how to apply the seal break test to opaque traps. One way to achieve this would be not to judge instantaneous seal break by visually, but to use a water level indicator, and regard seal break as having occurred when water level reaches near seal depth mm when pressure is applied. In this study performance evaluation was only made in terms of negative pressure, where seal break is most likely to occur, as we used amplitudes of pressure waves based on atmospheric pressure. Therefore it is necessary to develop a method of performance evaluation in terms of positive pressure by applying bias inside the apparatus. 418

8. References 1. Sakaue, K., Kamata,M., Kuroda, K., Wang Y. (21). Studies on Dynamic Characteristics of Trap Seal, Proceedings of CIB W62 International Symposium (pp. E4-1-4-5) 2. Sakaue, K., Kamata, M., Zyang, Y. (27). A study on the Test method of Trap performance, Proceedings of CIB W62 International Symposium (pp.321-332) 3. Kuriyama H., Sakaue, K., Yanagisawa, Y., Kamata, M., Sudo, H. (26), Studies on the Test method of Trap performance (Part 11), Technical papers of Annual meeting, SHASE (pp.793-796). 4. Kuriyama, H., Sakaue, K., Yanagisawa Y., Kamata, M., Iizuka, H. (28). A study on the Test method of Trap performance using Simple Test Apparatus, Proceedings of CIB W62 International Symposium (pp.252-263) 5. Sakaue, K., Shinohara, T., Kaizuka, M. (1977). A Study on the Dynamic Characteristics of the Trap Seal, Proceedings of CIB W62 International Symposium. 6. Sakaue, K., Shinohara, T., Kaizuka, M. (1982). A Study on the Dynamic Characteristics of the Seal in Deformed Traps, Proceedings of CIB W62 International Symposium. 7. Tomonari, H., Wang, Y., Kamata, M., Sakaue, K. (22). Numerical Studies of Seal Movement in Traps, Transactions of SHASE (pp.87-96) 419