ML20204G817

From kanterella
Revision as of 03:59, 7 December 2021 by StriderTol (talk | contribs) (StriderTol Bot change)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
Technical Evaluation of Dresden Nuclear Power Station Units 2 & 3,Plant-Unique Analysis Rept
ML20204G817
Person / Time
Site: Dresden  Constellation icon.png
Issue date: 09/30/1984
From: Bienkowski G
BROOKHAVEN NATIONAL LABORATORY
To:
Office of Nuclear Reactor Regulation
Shared Package
ML17195A951 List:
References
CON-FIN-A-3713 BNL-04243, BNL-4243, NUDOCS 8410040404
Download: ML20204G817 (26)


Text

_ _ _ _ _ .

i

\ >

TECHNICAL EVALUATION OF THE DRESDEN NUCLEAR GENERATING STATION UNITS 2 AND 3 PLANT-UNIQUE ANALYSIS REPORT George Bienkowski John R. Lehner Reactor Safety Licensing Assistance Division -

Departtunt of Nuclear Energy Brookhaven National Laboratory Upton, New York 11973 September 1984 FIN A-3713 ---

BNL-04243

.- (f 40 4t000 g

l ABSTRACT w

This Technical Evaluation Report (TER) presents the results of the Post-implementation audit of the Plant Unique Analysis Report (PUAR) for the Dresden Eclear Generating Station Units 2 and 3. The contents of the PUAR were compared against the hydrodynamic load Acceptance Criteria (AC) contained in NUREG-0661. The TER summarizes the audit findings (Table 1), and discusses the nature and status of any exceptions to the AC, identified during the audit (Table 2).

W W*

I O

m oo S

e

-- m ---- _ _ _ __ _ _ _ -- - -- - - . - - - .. .

ACKNOWLEDGEMENTS The cognizant NRC Technical Monitor for this program was Dr. fa'rouk Eltawila of the Containment Systems Branch (DSI) and the NRC Project Manager was Mr. Jack N. Donohiw of the Technical Assistance Program Management Group of the Division of Licensing. Mr. Byron Siegel of Operating Reactors Branch No. 2 (DL) was Lead Project Manager.

O O

mph e

o O

O l

i , _.

l ~- - __

l List of Acronyms AC Acceptance Criteria BNL Brookhaven National Laboratory BWR Boiling Water Reactor

~

CO Condensation Oscillation DBA Design Basis Accident DL Division of Licensing DSI Division of Systems Integration FSI Fluid Structure Interaction FSTF Full Scale Test Facility LDR Load Definition Report LOCA Loss-of-Coolant Accident LTP Long Term Program NRC Nuclear Regulatory Comission PUA Plant-Unique Analysis ,

PUAR Plant-Unique Analysis Report QSTF Quarter Scale Test Facility RFI Request For Information _

SER Safety Evaluation Report SPTMS Suppression Pool Temperature Monitoring System S/RV Safety / Relief Valve S/RVDL Safety / Relief Valve Discharge Line STP Short Term Program TAP Torus Attached Piping TER Technical Evaluation Report

-- 111

i l

Table of Contents

'Page No .

ABSTRACT f 11 ACKNOWLEDGEMENTS LIST OF ACRONYMS iii 1

1. INTRODUCTION 3
2. POST-IMPLEMENTATION AUDIT

SUMMARY

11

3. EXCEPTIONS TO GENERIC ACCEPTANCE CRITERIA
4. CONCLUSIONS 18
5. REFERENCES 19 h
  • a a

O er e-o 9

-iv-

~

I

=

1. INTRODUCTION The suppression pool hydrodynamic loads associated with a postulated loss-of-coolant accident (LOCA) were first identified during large-scale l testing of an advanced design pressure-suppression containment (Mark III). These additional loads, sich had not explicitly been included in the original Mark I containment design, result from the dynamic effects of drywell air and steam being rapidly forced into the suppression pool (torus). Because these hydrody-namic loads had not been considered in the original design of the Mark I con-tainment, a detailed reevaluation of the Mark I containment system was required.

A historical development of the bases for the original Mark I design as well as a summary of the two-part overall program (i.e., Short Term and Long Term Programs) used to resolve these issues can be found in Section 1 of Refer-ence 1. Reference 2 describes the staff's evaluation of the Short Term Program (STP) used to verify that licensed Mark I facilities could continue to operate safely while the Long Term Program (LTP) was being conducted.

The objectives of the LTP were to establish design-basis (conservative)

' loads that are appropriate for the anticipated life of each Mark I BWR facility-  :

(40 years), and to restore the originally intended des +gn-safety margins for

- each Mark I containment system. The principal thrust of sne LTP has been the development of generic methods for the definition of suppression pool hydrody-namic loadings and the associated structural assessment techniques for the Mark l I configuration. The generic aspects of the Mark I Owners Group LTP were com-pleted with the submitt'al of the " Mark I Containment Program Load Definition Re-port" (Ref. 3) and the " Mark I Containment Program Structural Acceptance Guide" (Ref. 4), as' well as supporting reports on the LTP experimental and-4nalytical tasks. The Mark I containment LTP Safety Evaluation Report (NUREG-0661)

< r presented the NRC staff's review of the generic suppression pool hydrodynamic load definition and structural assessment techniques proposed in the reports cited above. It was concluded that the load definition procedures utilized by the Mark I Duners Group, as modified by NRC requirements, provide cbnservative estimates of these loading conditions and that the structural acceptance crite-ria are consistent with the requirements of the applicable codes and standards.

The generic analysis techniques are intended to be used to perform a plant-unique analysis (PUA) for each Mark I facility to verify compliance with the acceptance criteria (AC) of Appendix A to NUREG-0661. The o' o jective of this study is to perform a post-implementation audit of the Dresden plant-unique analysis (Reference 5) against the hydrodynamic load criteria in NUREG-0661.

W O

Gm

  • e O

9

==e O

O _.

i

l. - - - - - - - . - . - . .

, 2. POST-IMPLEMENTATION AUDIT

SUMMARY

The purpose of the post-implementation audit was to evaluate the hydrodyna-mic loading methodologies which were used as the basis for modifying the pres-sure suppression system of the Dresden Nuclear Generating Station Units 2 and

3. The Dresden PUAR methodologies (Reference 5) were compared with those of the LDR (Reference 3)'as approved in the AC of NUREG-0661 (Reference 1). The audit l

procedure consisted of a moderately detailed review of the plant unique analysis report (PUAR) to verify both its completeness and its compliance with the accep-tance criteria. A list of requests for further information was submitted (Ref-erence 6), and answers were obtained at a meeting with the licensee (Reference 7).

Table 1 summarizes the audit results. It lists the various load categories specified in the AC, and indicates plant-unique information through the refer-ences, in the right-hand column, to the notes which follow in the text.

m 1

e d

9 m..

9

! 0 I ..

~

CRITERIA -

g z' j wz l' w9 os g4 m -Qu 2 4 m

w .

no a J 4 .

  • 4 #

LOADS CONTAINMENT PRESSURE a TEMPERATURE 2.1 / ,

VENT SYSTEM THRUST LOADS 2.2 /

i POOL SWELL TORUS NET VERTICAL LOADS 2.5 / '

. 1' TORUS SHELL PRESSURE HISTORIES 2.4 /

VENT SYSTEM IMPACT AND DRAG 2.6 / /

IMPACT AND DRAG ON OTHER STRUCTURES 2.7 /

FROTH IMPINGEMENT 2.8 / g t 1 POOL FALLBACK 2.9 /

i LOCA JET 2.14.1 / 3 LOCA BUBBLE DRAG 2.I4.2 / 3 VENT HEADER DEFLECTOR LOADS 2.10 /

, r.-

t TABLE 1. LOAD CHECKLIST FOR POST-lMPLEMENTATION AUDIT

' . . f l

- g f, mw $y ,

M4 I 3 3' '

f',

Iu4 WQ2 J 4 /y //

3m4 a#

s A

I R

E T

I R )

C // / / D E

U N

I T

1 5 32 1

6. 3 2 N 29s 0# 1 4 1 1 242 2 O 1 1 1 1 1 1 I.

1 C

5wo

( 2Z 2 2 2 2 222 2 (

1 E

L

- B A

T S S E E R R S

, N U U R T T E O

I C C M T U S U O A R D R C L T A T N S S O S W D L L A

l D D O C S E S C SES D O S L O D AR G ADM 1,

DGD ARA N

_ A OE O O N LMLOOE N LM LSD O Y I

T L B U MD L BM A L U E O E 'S EA L

A TE G_ E S TL S HNS M N HNS _

N S OYS O I G

SO Y S A L  !

E S S C SS R

- D UD TN G UDTE N RANW U RAN T O OOEO H OO EA _

C TLVD C TLVL _

.t

- a' i . i ' , J .

l .

t .-

I CRITERIA gzo j WI m m No m os 94 z< w bw $

MET NOT @ $O gm MET a F& z 20 a 4 4 ,

24 4 -

LOADS T-QUENCHER LOADS DISCHARbE LINE CLEARING 2.13.2 /

TORUS SHELL PRESSURES 2.13.3 / 6 i

& JET LOADS 04 SUBMERGED STRUCTURES 2.14.3 /

AIR BUBBLE DRAG - 2.14.4 / S7 j THRUST LOADS ON T/0 ARMS 2.13.5 /

S/RVDL ENVIRONMENTAL TEMPERATURES 2.13.6 /

l

! , tn .

I I

i 9

l TABLE 1. (CONTINUED) f

1 .

I, -

a' g' . ,

CRITERIA WI 52 3 to 9 OH . p4 m 'Do z< m w

k a J<

ou 4 #

Z4 DESCRIPTION SUPRESSION POOL TEMPERATURE 2.13.8 /

LIMIT I

SUPRESSION POOL TEMPERATURE 2.13.9 / 6
MONITORING SYSTEM t h DIFFERENTIAL PRESSURE CONTROL SYSTEM FOR THOSE PLANTS USING A 2.16 DRYWELL-TO-WETWELL PRESSURE DIFFERENCE AS A POOL SWELL

/ 9

MITIGATOR i SRV LOAD ASSESSMENT BY IN-PLANT TEST 2.13.9 / fo l ,

i .-

4 . TABLE 1. (CONTINUED) i

~

Notes to Table 1 Number 1 The Acceptance Criteria do not provide a separate procedure for calculating pool swell impact on spherical structures such as the main vent-to-vent header junction in Dresden. In the PUAR, the spherical junction was modeled as a series of cylinders udth axes along the main vent centerline. Acceleration drag, buoyancy and velocity drag were calculated using AC methodology for cylinders. This procedure was found acceptable.

2 For some structures, Region I froth loads were calculated using the i high-speed QSTF movies. This alternative is outlined in Appendix A of the AC.

3 Instead of the equivalent cylinder procedure specified in the AC to calculate acceleration drag volumes on sharp cornered submerged structures, the PUAR selected alternate modeling of the structures and ,

used published acceleratf on volumes. The discussion in Section 3.1 explains why this procedure was found acceptable.

~

- - - 4 To calculate CO and post-chug loads on the torus shell as well as on submerged structures, the 50 individual load harmonics were combined using a random phasing technique instead of the absolute summation specified in the AC. The discussion of Section 3.2 describes why this alternate met' hod was found acceptable.

m we h

o

-~ i 1

}

>bmber 5 To account for FSI effects during CO and chugging submerged structure loads, the AC suggested adding torus boundary acceleratiohs directly to local fluid accelerations. Instead,theapplicantuseI'amethod which calculated FSI acceleration fields anywhere in the torus based on knowing the boundary accelerations. This method, which has been accepted during previous PUAR reviews, is discussed in Section 3.3.

6 The analytical model to calculate SRV torus shell loads approved in the AC was modified slightly before being applied to Dresden. The purpose of the modifications was to more closely bound the pressure traces observed in the Monticello tests on which the model is based.

These changes have been found acceptable. SRV tests conducted in the Dresden plant further confirm that the analytically obtained loadings are conservative.

7 For SRV air bubble drag loads, the applicant reduced the AC bubble pressure bounding factor of 2.5 to 1.75. This still bounded peak positive bubble pressure and maximum bubble pressure differential from the Monticello test data. Dynamic load factors were derived from Dresden's'in-plant SRV test data. These modifications have been found

~

acceptable and are discussed in more detail in Section 3.4.

8 The new SPTMS is acceptable. As stated in Section 1-5.2 of the p0AR, the applicant has committed to perform a separate analysis demon-strating that' delayed operator action based on SPTMS readings will not cause the suppression pool temperature to exceed the limit specified in NUREG-0783. . _ _

Number 9 In order to reduce torus shell pressures caused by DBA pool swell, a minimum positive pressure difference of 1.0 psi is saintat'ned between a-W the Dresden drymell including the vent system, and the torcs air space. According to Technical Specifications, the plant is required to come'to shutdown if the main ap system fails.

10 SRV tests performed in the Dresden plant, were used to confirm that the analytically derived SRV shell loads are conservative and to deduce dynamic load factors for submerged structures.

e A

O 6

4 e-.

4 l

4

3. EXCEPTIONS TO GENERIC ACCEPTANCE CRITERIA Dresden Units 2 and 3 are two of'several plants analyzed by NUTECH Engineers, Inc. based on an essentially common hydrodyr.amic loading methodology (Fermi, Duane Arnold, Monticello and Quad Cities are other plants in this group). The methodology differs from the generic acceptance criteria of NUREG-0661 in four major areas which are listed in Table 2.

In what follows, each of these areas is discussed in detail, and the bases for the resolutions of the differences indicated.

i O

Op*

e 9

O O

m ee

Table 2: Issues Identified During Audit as Exceptions to the Generic Acceptance Criteria  !

c b-Issue No. Description Status Resolved Open

1. Use of acceleration drag volumes which X differ from those approved in the AC to determine drag on sharp cornered struc-tures.
2. Phasing of load harmonics used to analyze X structures affected by C0 and post-chug loads.
3. FSI methodology used for CO and chugging X submerged structure loads.

--- 4. Use of calibration factors developed from X Dresden in-plant tests for use in defining SRV submerged structure drag loads.

8 e

0

- -- n._-____-_-__. _ _ _ _ _ _ . _ _ _ _ _ _ _ _ _ _ _ _ _

3.1 Acceleration Drag Volumes for Sharp Cornered Structures The Acceptance Criteria 2.14.2 Section 2b in NUREG-0661 states that drag forces on structures with sharp corners (e.g. rectangles and "I" b ams) must be computed by considering forces on an equivalent cylinder of diameter 1

Deq=2 /2 Qax where Leax is the maximum transverse dimension. The intent of this criterion is to provide a conservative bound (based on very lim-ited data) that includes non-potential flow effects such as vortex shedding on both the acceleration drag due to hydrodynamic mass and the " standard" drag pro-portional to velocity squared. Since the dominant load for the Ring Beam (the primary non-cylindrical structure) is acceleration drag, the issue concerns only the hydrodynamic mass or acceleration volume and not the drag coefficient in the Dresden plant-specific case.

The PUAR states that " published" acceleration drag volumes listed in Table 1-4.1-1 are used for sharp edged structures rather than the equivalent cylinder i specified in the acceptance criteria. The detailed response to a Request for 4

Information (Item 1) explains that modeling of the actual structures is neces-sary, and in particular, forces on the web of the ring beam are obtained by modelling the beam by a circumscribed rectangle. In order to evaluate the implications of this modelling, sample calculations were performed on the ring- .

beam for the post-chug loading condition.

A direct application of the Dresden PUAR methodology leads to an accel-eration volume of 17.3 ft3 for in-plane forces on the maximally loaded ring beam sepent. A model that more accurately represents the interference effect but uses acceleration volumes from Table '1-4.1-1 gives a volume of 9.27 ft3, thus providing a substantial margin for possible non-potential flow 1tffects.

For the out-of-plane forces, the PUAR model given an acceleration volume of 59.9 ft 3 thile a similar modelling of the structure but using more realistic l

l

3 interference corrections yields a transverse acceleration volume of 53.3 ft ,

While this leaves very little margin for non-potential flow correct 1ons, in the _

parameterrangeofCOandpost-chugaccelerationspectrumwherethe{majorenergy In addition, is concentrated, the flow is expected to be very nearly potential.

the use of single. mode dynamic load factors, as explained in response to RFI The item 2, provides $dditional substantial conservatism up to a factor of 2.

conservative application of the AC equivalent cylinder model and interference correction to the hydrodynamic volume alone, mhile retaining the real volume for the " effective buoyancy" effect, gives an effective acceleration volume of 87.4 ft ,3 which yields out-of-plane loads 46% higher than those predicted by the PUAR.

Because of the parameter range in the Dresden plants and the conservative application of these loads, the potential non-conservatism on the acceleration volumes is adequately balanced by the conservatisms in the interference correc-tions and the load application.

On the basis of these comparisons we conclude that while the direct use of

" published" acceleration volumes for sharp edge structures may not in general

~

lead to conservative loads, the PDAR methodology for the application of these loads to the relevant structures, has sufficient conservatism to bound any hy-drodynamically produced stresses that could arise in these structures. .

3.2 CO and Post-Chug Harmonic Phasing The DBA condensation oscillation and the post-chug load definitions on the torus shell and on subme rged structures, accepted in the NUREG-0661, were based on data from a series of blowdomms in the FSTF facility (NEDE-24539), subject to l

i additional confirmatory tests reported in' the General Electric Letter Report J * - - .

M1-LR-81-01 of April 1981.

i i

I

l I

The condensation oscillation load definition as described in NEDD-21888 is l based on taking the absolute sum of 1* Hertz components of a spectrum from 0 to 50 Hz. Three alternative spectra are to be calculated with the one[ producing w

maximum response used for load definition. The procedure was found acceptable in the supplement to the SER (NUREG-0661), because the demonstrated high degree i

! of conservatism associated with the direct summation of the Fourier components t of the spectrum was ' sufficient to compensate for any uncertainties concomitant with the data available. The post-chug load definition is based on bounding i

i FSTF chugging data but otherwise follows similar procedures to those used in the CO load definition.

The PUAR uses a factor of .65 to multiply the C0 and post-chug loads com-puted on the basis of the absolute sum of the harmonic components. The justifi-1 j cation is based on comparisons of measured and predicted stresses in the FSTF facility using statistical studies of different phasing models (References 8, 9, 10,11). The factor .65 is chosen to give 84% non-exceedance probability with a confidence level of 90%. The PUAR does use an additional spectrum, Alternate 4 for the CO loading, based on test M12 from the supplementary tests reported in -

the letter report M1-LR-81-01. The information in Table 1-4.1-4 of the PUAR

! -- provides additional ' justification to show that the computed loads (using the .65.

i factor and Alternates 1 through 3) bound the measured stresses at critical points in the FSTF facility by 11% for axial shell stress to 59% for column i

force. The use of Alternate 4 in the Dresden plants provides an additional conservatism of about 20% to the shell response. The use of random phasing in the time domain for TAP (Volumes 6 and 7 of the PDAR) coupled with a factor of ,

1.3 for alte'rnates 1, 2, and 3 and a factor 1.15 for alternate 4 is-consistent with the results of Reference 10 and conservatively bounds all FSTF data.

i l

l

. i e

  • The procedures are a conservative application of the phasing design rules evaluated in Reference 12 and are therefore found acceptable.

3.3 FSI Methodology for C0 and Chugging Drag Loads [-

u-A detailed discussion of the method used to account for FSI effects on con-densation oscillation and chugging submerged structure loads is provided in Ref-erence 13. The mIthodology described in this note is used to compute accelera-tion fields across a submerged structure anywhere in the torus resulting from The method is presented FSI, based on knowing the torus boundary acceleration.

as an alternative to the NRC Acceptance Criteria suggestion of adding the bound-ary accelerations directly to the local fluid acceleration to account for FSI effects since the latter is deemed too conservative.

The review of the method outlined in Reference 13 has shown it to be rea-sonable and acceptable. The equations derived for fluid accelerations and pres-sure fields are plausible approximations fcr the conditions prevailing in the suppression pool. Assumed boundary conditions including the driving one at the torus wall are suitable. Overall trends as well as the acceleration fields de-Therefore, the alternate pro-

~ picted in the selected results appear reasonable.

cedure used to account for FSI effects on submerged structures is considered ac-

-- ceptable in this ap' plication. _

3.4 Calibration of SRV Drag Loads Based on In-Plant Tests The staff requested clarification of the detailed procedures used to derive the calibration factors from in-plant tests for SRV submerged-structure loads.

On the basis of this response, as well as those provided in other PUAR reviews of NUTECH plants, the staff considers the procedures as an acceptable modifica-

~~

tion of the AC.

l l

_ , _ ~ - -

The SRV bubble pressure data from Monticello tests is shown to be bounded using a bounding factor of 1.75 instead of the 2.5 specified in the AC. In the Dresden plants, dynamic load factors are derived on the basis of injplant tests. A bounding DLF value of 2.5 is then used for all submerged kructures.

The staff considers these procedures to be a reasonable application of the in-plant test results, and considers any potential uncertainties associated with the limited data base to be bounded by other conservatisms associated with the design load calculation procedures.

e WWD o

l

= = .

_17 l

4. CONCLUSIONS i

A post-implementation pool dynamic load audit of the Dresden PUAR has been completed to verify compliance with the generic acceptance criterta[of NUREG-0661. Four major differences between the PUAR and the AC merbidentified along with some other minor issues needing additional clarification. Based on additional information supplied by the applicant, as detailed in the previous section, all of these issues were resolved. The review of the Dresden PUAR has been completed with no issues or concerns outstanding.

4 e

9 j hm

  • D O

9 0

, _. 1 i

e MV &

- . - , - . . . , , , . -n,-,.

5. REFERENCES References cited in this report are available as follows: .

Those items marked with one asterisk (*) are available in the .WRC Public ~

Document Room for inspection; they may be copied for a fee.

Material marked with two asterisks (**) is not publicly available because it contains proprietary information; however, a nonproprietary version is avail-able in the NRC Public Document Room for inspection and may be copied for a fee.

Those reference items marked with three asterisks (***) are available for purchase from the NRC/GPO Sales Program, U. S. Nuclear Regulatory Comission, Washington, D. C. 20555, and/or the National Technical Information Service, Springfield, Virginia 22161.

All other material referenced is in the open literature and is availeble through public technical libraries.

(1) " Safety Evaluation Report, Mark I Long Tenn Program, Resolution of Generic Technical Activity A-7", NUREG-0661 July 1980.***

(2) " Mark I Containment Short-Term Program Safety Evaluation Report".

NUREG-0408, December 1977.***

(3) General Electric Company, " Mark I Containment Program Load Definition Re-port", General Electric Topical Report NEDO-21888, Revision 2, November 1981.*

(4) Mark I Owners Group " Mark I Containment Program Structural Acceptance Cri-teria Plant-Unique Analysis Applications Guide, Task Number 3.1.3", General Electric Topical Report NEDO-24583, Revision 1. July 1979.* .

i (5) "Dresden Nuclear Power Station Units 2 and 3 Plant-Unique Analysis Report",

Vols.1-7, Prepared for Connonwealth Edison Company by NUTECH Engineers, Inc., May 1983'. **

(6) Attachment to Letter from J. R. Lehner, BNL to F. Eltawila, NRC,

Subject:

Dresden Units 2 & 3 and Quad Cities Units 1 & 2 Request For Information, 5

June 27, 1984 (7) R. Rybak to H. Denton letter dated August 24,1984, " Response to Questions concerning Mark I Containment Plant Unique Analysis, NRC Docket Nos. 50-237/249 +

and 50-254/265". -

(8) General Electric Company " Mark I Containment Program, Evaluation of Har-monic Phasing for Mark I Torus Shell Condensation Oscillation Loads",

NEDE-24840, prepared for GE by Structural Mechanics Associates, October 1980. _ _ .

l (9) " Evaluation of FSTF Tests M12 and M118 Condensation Loads and Responses",

!MA12101.04-R0010, prepared by Structural Mechanics Associates for General Eiectric Company,1982.

M E

, , . , , - - - . , , , . . , , - ,_--,n. , ,., - - .__.. ,. , - , .-

)

(10) R. P. Kennedy, 'RIsp:nsa Factors Apprcpriate fer Use with C0 Harmonic Re-sponse Combination Design Rules," SMA12101.04-R002D, prepared by Structural Mechanics Associates for General Electric Company, March 1982.

(11) R. P. Kennedy, "A Statistical Basis for Load Factors Appropriate for Use with C0 Harmonic Response Combination Design Rules " SMA 12101404-R003D, prepared by Structural Mechanics Associates for General Electrft Company, March 1982.

(12) G. Bienkowski, " Review of the Validity of Random Phasing Rules as Applied to C0 Torus Loads". Internal BNL Memo August 1983.

(13) A. J. Bilanin, " Mark 1 Methodology for FSI Induced Subme:rged Structure Fluid Acceleration Drag Loads", Continuum Dynamics Tech. Note No. 82-15 June 1982.*

J l

l 4

O e

O m ee

~y ,

-. . . ,--,