DTIC ADA261671: Investigation of Liquid Sloshing in Spin-Stabilized Satellites.

214 Pages·8.1 MB·English

by  Defense Technical Information Center

"
AD-A261 671
Final Report
Grant #AFOSR-89-0403
INVESTIGATION OF LIQUID SLOSHING
IN SPIN-STABILIZED SATELLITES
Joseph R. Baumgarten
Donald R. Flugrad
Richard H. Pletcher
Engineering
DTIC Department of Mechanical
(cid:127)
W-LRECT- E January 31,1993
. - This report has been prepared for the
Air Force Office of Scientific Research
93-04639
ISU-ERI-Ames-93113
I''~ll~ lI IIl iiilllP°' '
1(cid:127)-
engi.neenng
researTc institute
iowa state umve
OvR op? blic 8elesel 34 015
Dis ion iited
ECURITY CLASSlFICA(cid:127)ION OF THIS PAGE
REPORT DOCUMENTATION PAGE
la. REPORT SECURITY CLASSIFICATION lb. RESTRICTIVE MARKINGS
Unclassified
2a. SECURITY CLASSIFICATION AUTHORITY 3. DISTRIBUTION IAVAILABIUTY OF REPORT
2b. DECLASSIFICATION / DOWNGRADING SCHEDULE Approved for public release;
distribution unlimited
4. PERFORMING ORGANIZATION REPORT NuMBER(S) S. MONITORING ORGANIZATION REPORT NUMBER(S)
ISU-ERI-Ames-93113 Project 3166 --
6;. NAME OF PERFORMING ORGANIZATION 6b. OFFICE SYMBOL 72. NAME OF MONITORING ORGANIZATION
Iowa State University (/f app/icab/e) Air Force Office of Scientific Research
Dept. of Mechanical En(cid:127)ineerin(cid:127) NA (cid:127)AFOSR(cid:127)
Ec. ADORESS (Oty, State, end ZIP Code) 7b. ADORESS (City, State, end ZIP Code)
Air Force Office of Scientific Research
Room 2025 Black Engineering Building Direct'orate of Aerospace Science (AFOSR/NA)
Ames, IA 50011
Bolling AFB, DC 20332
:(cid:127)IL NAME OF FUNDING/SPONSORING 8b. OFFICE SYMBOL 9. PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER i
S
"ADDRESS (Dry, State;BnO'ZlPCocie) . .... 10..SOURCE.OF FUNDING NUMBERS
Bldg. 410 PROGRAM PROJECT TASKNo. (cid:127) i WORK UNIT
Bolling AFB, De 20332-6448 NO/(cid:127) NO.
F i1. TITLE (IlXluO(cid:127) Je(urn(cid:127)/ Cl,e sslficatJon)
Investigation of Liquid Sloshing in Spin-Stabilized Satellites (Unclassified)
;12. PERSONAL AUTI'IOR{.€)
Baumgarten, Joseph R.; Flugrad, Donald R.; Pletcher, Richard H.
1.t8. TYPE OF REPORT 13b. TIME COVERED j14. DATE OF REPORT (Y,,r.Mont(cid:127).oay) ts. PAGE COUNT
Final Report FROM 7/1/E9 TO 11./(cid:127)n/c.(cid:127) 1(cid:127)3 January 31 ?.g9
16. SUPPLEMENTARY NOTATION
17. COSATI CODES 18. SUBJECT TERMS (Continue on reverie if necessery end i(cid:127)ent(cid:127)fy by INock number)
FIELD GROUP SUB-43ROUP Coning, nutation, and structural deflection of satellites
due to sloshing fluid, fluid-structure interaction
19. ABSTRACT (Continue on reverie if necessary 4rid #o(cid:127)ntify by block number)
Several spin-stabilized communication satellites with liquid stores on board have exhibited
unstable nutational motion. Liquid sloshing is suspected as the cause of this undesirable
behavior. During an initial three year grant period, a test rig was built and instrumented,
a rigid body computer model was developed using a pendulum analogy to simulate the sloshing
liquid, and computational fluid dynamic (CFD) methods were used to develop a primitive
variable numerical algorithm to describe two and three dimensional liquid sloshing.
During the current three year grant period, the test rig and rigid body computer model have
been used to study the effect of various physical parameter values on the motion, and
stability conditions have been determined. The structural mathematical model has been
improved by introduction of finite element techniques to account for elastic deformation of
the test rig, and an impliclt-expllcit numerical approach has been implemented to solve the
coupled, nonlinear equations of motion. (Continued on reverse side)
20. DISTRIBUTION IAVAILABILITY OF ABSTRACT 21. ABSTRACT SECURITY CLASSIFICATION
[(cid:127)UNCLASSIFIED(cid:127)/NLIMITED I(cid:127) SAME AS RPT (cid:127)OTIC USERS Unclassified
i
222. NAME OF RESPONSIBLE INOI/VIOUA.L 221). TELEPHONE (IrwJl,(cid:127)le Area CoDe) 22c. OFFICE SYMBOL
DD FORM 1473, 84 MAR B3 APR edmon may be used untI exhiusted.J SEOaRITY CLASSIFICaTION OF THIS PAGE
ssified
All Other editions are obSsoUlentec.
iii
19. Two CFD numerical strategies for computing liquid sloshing have been developed.
One is based on surface fitting, and the other is based on surface capturing.
Both provide solutions to the full three-dimensional unsteady incompressible
Navier-Stokes equations which govern the liquid motion.
The structural and CFD software packages have been merged to study the fluid-
structure interaction. Several spin-up cases have been computed, and results have
been compared with experimental data.
iv
"Accesiorl
For
NTIS CRA&I
DTIC TAB
Un)announced
JUStitication
--------------------- -..----------------
DBi-S--tr-i-b-t-i-ti-o-n
Availability Codes
Dist Ava cindIo
Special
TABLE OF CONTENTS
ABSTR ACT ................................... iii
LIST OF TABLES ........................................ ix
LIST OF FIGURES ....................................... xi
1. INTRODUCTION ............................. 1
1.1 Research Objectives ............................ 3
1.2 Status of Research ............................ 4
1.3 Publications . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4 List of Research Personnel, Thesis Titles, and Degrees Awarded ... 10
1.5 Seminars, Presentations, and Laboratory Visits ................. 11
2. RESEARCH IN PROGRESS ............................. 15
2.1 Recent Progress in the Surface Fitting Approach ................ 15
2.1.1 Free-surface motion ............................... 15
2.1.2 Contact line boundary conditions ...................... 18
2.1.3 Code refinement efforts ............................ 19
2.2 Recent Progress in the Surface Capturing Approach ............. 21
2.2.1 Formulation of the governing equations ................. 22
2.2.2 Numerical methods ............................... 23
V
2.2.3 Test case results ................................. 26
2.2.4 Concluding remarks ............................... 28
2.3 Progress in Fluid-Structure Interaction ...................... 28
2.3.1 Equipment .................................... 28
2.3.2 Calibration ........ ............................ 42
2.3.3 Experiments ................................... 44
2.3.4 Results ....................................... 45
2.3.5 Numerical simulation ............................. 46
2.3.6 Discussion of results ...... ....................... 47
2.3.7 Concluding remarks ............................... 48
BIBLIOGRAPHY ........................................ 59
APPENDIX A. Computer Simulation of a Test Rig to Model Slosh-
ing in Spin-Stabilized Satellites ............................ 61
APPENDIX B. Motion Study of a Spin-Stabilized Satellite Test Rig. 73
APPENDIX C. A Finite Element/Lagrangian Formulation Dynamic
Motion Prediction for a Flexible Satellite Simulator with Both
Rigid and Elastic Bodies ................................. 113
APPENDIX D. Modeling of Multibody Flexible Articulated Struc-
tures with Mutually Coupled Motions: Part I - General Theory. 123
APPENDIX E. Modeling of Multibody Flexible Articulated Struc-
tures with Mutually Coupled Motions: Part II - Application and
Results ............................................. 135
vi
APPENDIX F. A Sequential Implicit-Explicit Integration Method
in Solving Nonlinear Differential Equations from Flexible System
M odeling . ................ .... ............... 147
APPENDIX G. A Primitive Variable, Strongly Implicit Calculation
Procedure for Viscous Flows at All Speeds ............... 155
APPENDIX H. Simulation of Three-Dimensional Liquid Sloshing
Flows Using a Strongly Implicit Calculation Procedure ........ 167
APPENDIX I. A Numerical and Experimental Study of Three-
Dimensional Liquid Sloshing in a Rotating Spherical Container.. 205
Vii
viii
LIST OF TABLES
Table 2.1: Description of Equipment .................... 39
ix
x
LIST OF FIGURES
Figure 2.1: Control volume within a two-fluid system ............... 24
Figure 2.2: Schematic of the broken dam problem ................. 29
Figure 2.3: Density interface for the broken dam problem selected times. 30
Figure 2.4: Velocity field for the broken dam problem selected times. . . 31
Figure 2.5: Surge front position versus dimensionless time for the broken
dam problem .................................. 32
Figure 2.6: Back wall free surface position versus dimensionless time for
the broken dam problem .......................... 33
Figure 2.7: Schematic of oscillating tank problem .................. 34
Figure 2.8: Density interface for selected times: oscillating tank test case. 35
Figure 2.9: Free surface position versus time for front and back tank walls. 36
Figure 2.10: The satellite test rig ............................. 37
Figure 2.11: Strain Gage Conditioning Circuit ..... ............... 40
Figure 2.12: Photopotentiometer Circuit ........................ 41
Figure 2.13: Data Acquisition Terminal Block Connections ........... 42
Figure 2.14: Photopotentiometer 1 Calibration .................... 49
Figure 2.15: Photopotentiometer 2 Calibration .................... 50
Figure 2.16: Tachometer Calibration ........................... 51
xi