A Mathematical Model for the Fire-extinguishing Rocket Flight in a Turbulent Atmosphere
|
|
- Elfrieda Hubbard
- 5 years ago
- Views:
Transcription
1 A Mathematical Model for the Fire-extinuishin Rocket Fliht in a Turbulent Atmosphere CRISTINA MIHAILESCU Electromecanica Ploiesti SA Soseaua Ploiesti-Tiroviste, Km 8 ROMANIA crismihailescu@yahoo.com Abstract: The paper analyses the evolution of rockets used to extinuish fires developed in open space in conditions of turbulent atmosphere. After the eneral definitions of the turbulent atmosphere and workin hypotheses, the model of and turbulent atmosphere obtained from analytical relations is presented in detail. Startin from this model, the components of velocity of as functions of time are deduced. Introducin the components of velocity into the model of simulation of the rocket fliht, a fascicle of the trajectories is obtained. The work novelty consists in the new technique used to model the turbulent atmosphere influence on unuided rocket fliht. Key-Words: turbulent atmosphere, rocket trajectory, statistical function, dynamics of fliht, models Introduction Fire-extinuishin rocket trajectory must be evaluated in the early project desin stae, when the extreme atmospheric s, the vertical air flows and the no wind must be considered as possible perturbations. The effect of lare scale perturbation is important for establishin the averae trajectory. A supplementary task would be to determine the trajectory dispersion due to perturbation radients produced by the atmospheric usts. The importance of the study lies in the results it provides: the probable trajectory, its dispersion and the rocket response to dynamic loadin. This paper tries to provide an answer to the problem of turbulent atmosphere influence on fire-extinuishin rocket fliht. Definition of Atmospheric Turbulence. Hypotheses Since is a random process that cannot be described by explicit functions of time, only a statistical, probabilistic approach can be taken. The velocity field, variable in time and space, may be considered as result of an averae value and a perturbation around it. In the case of lon time flihts, the averae value may be modelled as referencin the fliht with respect to a mobile rectanular frame which moves with the wind s averae velocity. In a point located by r x, x, the air averae velocity is: ( x3 [ u u u ] T u( r, t 3, ( It may be described properly by a correlation matrix composed of elements like this: Ri, j ( ξ, τ < ui ( r, t u j ( r+ξ, t+τ >, ( which represents the temporal and spatial averae of the product between a component ui located at a point defined by r at the time t, and the component u from the point r + ξ at the time t + τ. To this matrix we can associate a spectral function matrix whose elements are defined by the Fourier transform function: Θ i( Ωξ+ωτ j ( Ω, ω ( ξ, τ dξ dξ dξ 4 j e 3 (π R dτ(3 The inverse of these elements belon to the correlation matrix: i( Ωξ+ωτ j ( ξ, τ Θ j ( Ω, ω e dω dω dω3 R dω (4 Generally, the atmosphere s are not Gaussian. However, for practical purposes, the normal repartition may be used, producin sinificant simplifications. On the other hand, the statistical parameters R i, j and j j Θ of the are defined in each point in space, r and variable with time t. A particular aspect of the fire-extinuishin rocket is that it flies a lon time, with a sliht chane in velocity manitude, due to slow burnin of rocket propellant. This is a sinificant theoretical advantae because it allows admittin the hypothesis that the is a stationary process. Another hypothesis is that the could be considered homoenous on small intervals, which means that R i, j and Θ j do not depend on r, alonside of at least a trajectory sement. At hih altitude, the is identical in all points belonin to the same layer, which means it miht be considered homoenous. Close to the round often chanes with altitude occur. However, in the case of small launch anles, homoenous miht be considered alon the trajectories close to round. ISSN: ISBN:
2 R i, j and j Θ usually depend on wind referenced frame axes orientation. If the statistical functions do not depend on space, the is considered isotropic, which makes the velocity components square means to be equal: u u u3, (5 where the standard deviation will be called intensity. Due to the fact that rocket velocity is much hiher than their variation velocity, the may be rearded to as a frozen model. This allows us to nelect time in the function u ( r, t, which corresponds to Taylor s hypothesis. Consequently, the correlation and spectral functions become: R j ( ξ, τ R j ( ξ ; Θi, j ( Ω, ω Θ j ( Ω, (6 Consequently, the simplest model is the one of homoenous and isotropic s, Gaussian and frozen, a model which is used for the fliht outside the round proximity layer. Inside of this layer is necessary to consider anisotropic. For the isotropic, usin Batchelor s relation: Rij ( ξ ξξ [ f ( ξ ( ξ] ξ i j + ( ξ δ where ξ ξ ; δ j - Kroneker s symbol; j, (7 The lonitudinal correlation function f (ξ represents the correlation between velocity components alon an axis, iven in two points of that axis (Fiure.: R ξ f ξ ( lim u( x u( x+ξ x d ( (8 Because of the hypothesis reardin isotropy: f ( ξ f ( ξ f ( ξ3 f ( ξ, (9 and: R( ξ lim f ξ u( x u( x+ξ d x ( ( x ξ u x ( x + ξ u ( x Fi. Lonitudinal correlation of velocity components The other function (ξ, is the lateral correlation function and represents the correlation between the transversal velocity components with respect to an axis (Fiure, defined in two points: x 3 u 3 u x R33( ξ lim ξ u3( x u3( x+ξd x ( ( Based on the hypothesis reardin isotropy, it results: < u3 ( x u3 ( x+ξ >< u ( x u ( x+ξ >, ( < u3 ( x u3 ( x +ξ >... and: R33( ξ lim ξ u3( x u3( x+ξ d x ( (3 u 3 ( x + ξ u ( x 3 + ξ 3 u 3 ( x u (x + ξ x 3 u ( x u ( x u ( x 3 + ξ 3 u 3 ( x + ξ x Fi. Transversal correlation of velocity components 3 The Uniform Turbulence Model For the characteristic correlation functions, the followin relations from work [] may be used: / 3 3 f ( ξ ζ / K/ 3( ζ ; Γ(/ 3 /3 4/3 ( ξ f ( ξ ζ K/ 3( ζ (4 Γ(/ 3 ξ where: ζ, a, 339 and L representin the al characteristic lenth. Γ is the second deree Euler function and K υ is the second deree Bessel modified function. The Fourier transforms of the functions f (ξ and (ξ are spectral one-dimensional functions: Φ iωξ ( Ω ( ξ e dξ π iωξ Φ33( Ω ( ξ e d π f ; ξ (5 Fiures 3 and 4 illustrate raphical representation of functions defined by (4 and (5. For these spectral functions, throuh experimental measurements, several models have been established, amon which, the most commonly used are the von Karman model: L Φ( Ω ; + (8 3( alω Φ33( Ω π [+ ( alω π [+ ( alω ] 5 6 ] and the Dryden model: L Φ( Ω π + ( LΩ ; Φ 33 x L (6 L + 3( LΩ ( Ω π [+ ( LΩ ] 6 (7 ISSN: ISBN:
3 f,.5.5 f H> 65 m L 76 m a ξ/l Fi. 3 Characteristic correlation functions Φ L. Φ Φ33 4 The Linear Turbulence Method A hiher order approximation is the one of a field, the velocities bein considered as functions dependin on position. Authors of [], [] and [3] developed in case of airplane models in four points. These models also allow calculation of rotation velocity components p, q, r due to. Startin from the method proposed in [], we developed for rockets a model of in two points, which allow deducin of anular velocity of pitch and yaw motion. In the case, due to the slender axisymmetrical shape of the rocket, some radients may be nelected. We shall consider a model of the rocket in two points as shown in Fiure 5. (C.m O l v (F. H> 65 m L 76 m Fi. 5 Model of the rocket in two points 5 Ω L 5 Fi. 4 Characteristic spectral functions The established experimental relations [] are in ood accordance with relations (5 and may be successfully used in the case of the study models with. However, in the case of method, due to necessity to define cross correlation functions which show the mutual dependence between velocity components, the characteristic functions defined by (4 will be used. In case of altitudes over 65m the boundary layer influence is diminished and the atmosphere may be considered, a unique value for the characteristic lenth L bein recommended, both to lonitudinal and lateral spectral functions. For the model considered, accordin to [] and [] L76m. The standard deviation is the same alon all the three directions, dependin only on the intensity:,53 [ m / s] - low ; 3,5 [ m / s] - mean ; 6,4 [ m / s] - heiht (storm. The model obtained in this way is the simplest possible, bein known as the field model, in which only the components of velocity alon the three axes, reduced in the ravity centre of the rocket, are considered: u, v, w. The lenth l v, called aerodynamic lenth, represents the distance from the mass centre to the focal point. Based on this model, the pitch and yaw velocities, equivalent to the velocity radients, are: q ( w w ; r ( v v, (8 lv lv to which there are added the velocities due translation for the : w w ; u u ; v v (9 To calculate the velocities produced by, we expressed their correlation functions throuh the characteristic functions f and, proceedin as in the case of distribution, we calculated the spectral functions correspondin to the correlation functions and, from these, we found some possible forms of the translation and rotation velocities enerated by. For the case of the pitch velocity the followin initial relation was used: Rqq ( τ < q ( t q ( t+τ >, ( obtainin: ( lv R ( τ ( ξ ( ξ ( ξ3 qq, ( where: ξ Vτ ξ + l ξ ξ l ; ξ v 3 v ξξ+l v ξ3ξ-l v O I F I O F ξvτ Fi. 6 Correlation functions calculus schema ISSN: ISBN:
4 ADVANCES in MATHEMATICAL and COMPUTATIONAL METHODS Similarly the anular yaw velocity: Rrr( τ < r ( t r ( t+τ > < ( v v ( v' v' > ( lv ( [ < vv' >+< vv ' > < vv ' > < vv' > ] ( lv from which we obtained: Rrr ( τ [R ( Vτ,, R ( Vτ+ l, (,] v R Vτ lv ( l v ( lv R ( τ ( ξ ( ξ ( ξ3 rr, (3 From the previous relations and (3, the correlation functions for the rotation velocities were determined. The two rotation velocity components are equal, R, due to axial symmetry (Fiure 7. The qq R rr spectral function for the rotation velocity is represented in Fiure 8. We compare the velocity components calculated with the field theory to those obtained usin the field theory. These proved to be identical. Rqq lv/ H> 65 m L 76 m a ξ/l* Fi. 7 Correlation function for rotation velocities Φ qq lv /( L * H> 65 m L 76 m a Ω L/ Fi. 8 Spectral function for the rotation velocity Thus, for the lonitudinal component we obtained: Ruu ( τ < u ( t u ( t+τ >< uu' > R( Vτ,, f ( ξ (4 For the lateral component we obtained: Rvv ( τ < v ( t v ( t+τ >< vv' > R( Vτ,, ( ξ (5 Rww ( τ < w ( t w ( t+τ >< ww' > R33( Vτ,, ( ξ (6 Relations (4 (6 show that for model we found the same characteristic f and from Fiure 3 and the spectral functions shown in Fiure 4. This means that the translation velocities for the liner model are the same with those from the model and for these we can use theoretical relation (4 or the experimental relations (6 and (7. Relations (3 for the anular velocities may be used only in the case of the theoretical model defined by the characteristic functions (4. In the case of axissymmetric rockets, the only supplementary term specific to relations (3 is produced by the anular pitch/yaw velocities. If the aerodynamic lenth is comparable with respect to the characteristic lenth, the influence of this term is small and can be nelected. In the case of short rockets, like fire-extinuishin rocket, the term is important as we will see from the numerical comparison. 5 Simulation of Atmospheric Turbulence The simulation of an atmospheric is based on the method of sinusoidal sum, usin the spectral function modulus of the velocity to obtain the oscillation amplitude specific to a frequency bandwidth. So, for the spectral functions values presented in Fiure 4, we determine by interation the value of amplitude correspondin to a frequency bandwidth, which can be approximated with the frequency correspondin to the centre of the interated domain. To complete the model, to each frequency we associated an initial arbitrary phase. Usin the notation S for the portion of the area φ obtained by interatin the modulus of a spectral function from Fiure 4 and considerin that the function has been represented only for the positive values of the variable Ω, the amplitude has the value: A S (7, φ Two of the translation velocities enerated by in a point ξ are iven by: i i ( ; i u ξ A sin( Ω ξ+ φ i i w ( ξ A sin( Ωξ+ φ, (8 i where the index was iven to amplitude, pulsation and initial phase obtained from the modulus of the lonitudinal spectral function, and the index iven to the elements obtained from the modulus of the lateral spectral function. The third component of the velocity is the same as the second, exceptin the initial phase: i i v ( ξ A sin( Ωξ+ φ3 (9 i Makin the assumption that the rocket flies on a trajectory with a constant velocity V: ISSN: ISBN:
5 ξ V τ, (3 Consequently, the velocities oriinated by may be expressed as functions of time (Fiure H> 65 m L76 m V 4 m/s u/ w/ t[s] Fi. 9 The axial u and normal w velocities produced by The modulus of the spectral functions presented in Fiure 8 was used to find the anular velocities. By applyin similar relations to (7 and (8 we obtained a form of the anular velocities oriinated by with respect to time. Because the spectral functions are equal (Fi.8, the yaw and pitch velocities will be approximately equal r q, differin only by initial phases. In this case, in Fiure, we presented only the pitch velocity function q. of trajectories was obtained. The results are presented in the next tables. In order to analyse trajectory parameters, we define four check points: -Start point, -End of rocket motor burn, 3-Trajectory apex, 4-Impact point Table. Nominal trajectory parameter, Stae of fliht 3 4 T [s] V [m/s] γ [de] [m] Y [m] Table. Averae and standard deviation trajectory parameters, low. 53 Stae of Uniform model Linear model fliht T [s] V [m/s] γ [de] [m] Y [m] q*lv/ H> 65 m L76 m V 4 m/s In table for each parameter are two rows: first for averae and the second for standard deviation values. For a better analyse the main results are presented in a comparative form t[s] Fi. The pitch velocity produced by V [m/s] nominal 6 Fire-extinuishin Rocket Trajectory Startin from papers [3] and [4], some details reardin the procedure of buildin the simulation model for the fire-extinuishin rocket are presented in this section. For the turbulent atmosphere we used in order to make a comparison, the and the model. The initial values of the sinusoidal sums phases were obtained by eneratin random numbers, and the amplitude and frequency were determined startin from the spectral density functions usin the procedure described in the previous section. Because of the sequential initial random input of the phases, a fascicle Fi. Velocity diaram averae values, low ISSN: ISBN:
6 [m/s] Fi. Velocity diaram standard deviation values, low [de] γ nominal Fi. 3 Climb anle diaram, averae values, low [de].5 7 Conclusions Two models were presented in the paper. In the model, only the translation velocities of the mass centre are considered. In the model, besides the velocities, the anular velocities of pitch and yaw with respect to mass centre are added. In the case of axis-symmetric rockets, the term dependin on anular velocities is the only which couples the equations. A further simplification could be made in the case of short rockets to which the aerodynamic lenth is much smaller than the characteristic lenth. In this situation, the couplin term due to rotations is small and may be nelected, so that the model is an acceptable approximation. Considerin the model developed in our study, we have analysed in a simulation example the fliht dynamics of a rocket, in which the influence of was evaluated in four characteristic staes of the fliht. The work novelty consists in technique used in modellin the turbulent atmosphere influence on unuided rocket fliht: two point model used for the the technique used to obtain the velocity the comparison between and model Fi. 4 Climb anle diaram, standard deviation values, low Y[m] nominal 5 [m] Fi. 5 Trajectory diaram, averae values, low The averae parameters from nominal trajectory, model and model are relatively closed. Different from these, the parameters dispersion obtained with model and model are different, therefore the use of the model are rihtful. References: [] Etkin, Bernard, Dynamics of Atmospheric Fliht, John Wiley & Sons, Inc., New York, 97. [] Etkin, Bernard, Lloyd Duff Reid, Dynamics of Fliht: Stability and Control -3 rd ed, John Wiley & Sons, Inc., New York, 995. [3] Chelaru T.V., Barbu C. Mathematical model and technical solutions for multi-stae soundin rockets, usin the rotation anles, Proceedins of the 3 rd International Conference on Applied Mathematics, Simulation, Modelin (ASM 9, pp.94-99, ISBN: , NAUN, Vouliamen Athena, Greece, December 9-3, 9. [4] Chelaru T.V., Barbu C., Chelaru A., Mathematical model and technical solution for small multistae suborbital launchers, Proceedins of the 4 th International Conference on Recent Advances in Space Technoloies-RAST 9, pp. 86-9, ISBN: , Istanbul, Turkey, -3 Jun. 9. ISSN: ISBN:
2.3. PBL Equations for Mean Flow and Their Applications
.3. PBL Equations for Mean Flow and Their Applications Read Holton Section 5.3!.3.1. The PBL Momentum Equations We have derived the Reynolds averaed equations in the previous section, and they describe
More informationModeling for control of a three degrees-of-freedom Magnetic. Levitation System
Modelin for control of a three derees-of-freedom Manetic evitation System Rafael Becerril-Arreola Dept. of Electrical and Computer En. University of Toronto Manfredi Maiore Dept. of Electrical and Computer
More informationThe Analysis of Dispersion for Trajectories of Fire-extinguishing Rocket
The Analysis of Dispersion for Trajectories of Fire-extinguishing Rocket CRISTINA MIHAILESCU Electromecanica Ploiesti SA Soseaua Ploiesti-Tirgoviste, Km 8 ROMANIA crismihailescu@yahoo.com http://www.elmec.ro
More information(a) 1m s -2 (b) 2 m s -2 (c) zero (d) -1 m s -2
11 th Physics - Unit 2 Kinematics Solutions for the Textbook Problems One Marks 1. Which one of the followin Cartesian coordinate system is not followed in physics? 5. If a particle has neative velocity
More informationAltitude measurement for model rocketry
Altitude measurement for model rocketry David A. Cauhey Sibley School of Mechanical Aerospace Enineerin, Cornell University, Ithaca, New York 14853 I. INTRODUCTION In his book, Rocket Boys, 1 Homer Hickam
More informationMotion in Two Dimensions Sections Covered in the Text: Chapters 6 & 7, except 7.5 & 7.6
Motion in Two Dimensions Sections Covered in the Tet: Chapters 6 & 7, ecept 7.5 & 7.6 It is time to etend the definitions we developed in Note 03 to describe motion in 2D space. In doin so we shall find
More informationGround Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions
PC11 Fundamentals of Physics I Lectures 7 and 8 Motion in Two Dimensions Dr Tay Sen Chuan 1 Ground Rules Switch off your handphone and paer Switch off your laptop computer and keep it No talkin while lecture
More informationONLINE: MATHEMATICS EXTENSION 2 Topic 6 MECHANICS 6.3 HARMONIC MOTION
ONINE: MATHEMATICS EXTENSION Topic 6 MECHANICS 6.3 HARMONIC MOTION Vibrations or oscillations are motions that repeated more or less reularly in time. The topic is very broad and diverse and covers phenomena
More informationXI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com
XI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K affan_414@live.com https://promotephysics.wordpress.com [MOTION IN TWO DIMENSIONS] CHAPTER NO. 4 In this chapter we are oin to discuss motion in projectile
More informationEnergizing Math with Engineering Applications
Enerizin Math with Enineerin Applications Understandin the Math behind Launchin a Straw-Rocket throuh the use of Simulations. Activity created by Ira Rosenthal (rosenthi@palmbeachstate.edu) as part of
More information2.2 Differentiation and Integration of Vector-Valued Functions
.. DIFFERENTIATION AND INTEGRATION OF VECTOR-VALUED FUNCTIONS133. Differentiation and Interation of Vector-Valued Functions Simply put, we differentiate and interate vector functions by differentiatin
More informationv( t) g 2 v 0 sin θ ( ) ( ) g t ( ) = 0
PROJECTILE MOTION Velocity We seek to explore the velocity of the projectile, includin its final value as it hits the round, or a taret above the round. The anle made by the velocity vector with the local
More informationWire antenna model of the vertical grounding electrode
Boundary Elements and Other Mesh Reduction Methods XXXV 13 Wire antenna model of the vertical roundin electrode D. Poljak & S. Sesnic University of Split, FESB, Split, Croatia Abstract A straiht wire antenna
More informationPHY 133 Lab 1 - The Pendulum
3/20/2017 PHY 133 Lab 1 The Pendulum [Stony Brook Physics Laboratory Manuals] Stony Brook Physics Laboratory Manuals PHY 133 Lab 1 - The Pendulum The purpose of this lab is to measure the period of a simple
More informationFOOTING FIXITY EFFECT ON PIER DEFLECTION
Southern Illinois University Carbondale OpenSIUC Publications Department of Civil and Environmental Enineerin Fall 9--04 FOOTING FIXITY EFFECT ON PIER DEFECTION Jen-kan Kent Hsiao Yunyi Jian Southern Illinois
More informationParametric Equations
Parametric Equations Suppose a cricket jumps off of the round with an initial velocity v 0 at an anle θ. If we take his initial position as the oriin, his horizontal and vertical positions follow the equations:
More informationAAPT UNITED STATES PHYSICS TEAM AIP 2009
2009 F = ma Exam 1 AAPT UNITED STATES PHYSICS TEAM AIP 2009 2009 F = ma Contest 25 QUESTIONS - 75 MINUTES INSTRUCTIONS DO NOT OPEN THIS TEST UNTI YOU ARE TOD TO BEGIN Use = 10 N/k throuhout this contest.
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2016 17 ENGINEERING MATHEMATICS AND MECHANICS ENG-4004Y Time allowed: 2 Hours Attempt QUESTIONS 1 and 2, and ONE other question.
More informationAdvanced Methods Development for Equilibrium Cycle Calculations of the RBWR. Andrew Hall 11/7/2013
Advanced Methods Development for Equilibrium Cycle Calculations of the RBWR Andrew Hall 11/7/2013 Outline RBWR Motivation and Desin Why use Serpent Cross Sections? Modelin the RBWR Axial Discontinuity
More informationDesign of Chevron Gusset Plates
017 SEAOC CONENTION PROCEEDINGS Desin of Chevron Gusset Plates Rafael Sali, Director of Seismic Desin Walter P Moore San Francisco, California Leih Arber, Senior Enineer American Institute of Steel Construction
More informationHEAT TRANSFER IN EXHAUST SYSTEM OF A COLD START ENGINE AT LOW ENVIRONMENTAL TEMPERATURE
THERMAL SCIENCE: Vol. 14 (2010) Suppl. pp. S219 S232 219 HEAT TRANSFER IN EXHAUST SYSTEM OF A COLD START ENGINE AT LOW ENVIRONMENTAL TEMPERATURE by Snežana D. PETKOVIĆ a Radivoje B. PEŠIĆ b b and Jovanka
More informationGeodesics as gravity
Geodesics as ravity February 8, 05 It is not obvious that curvature can account for ravity. The orbitin path of a planet, for example, does not immediately seem to be the shortest path between points.
More informationExperiment 3 The Simple Pendulum
PHY191 Fall003 Experiment 3: The Simple Pendulum 10/7/004 Pae 1 Suested Readin for this lab Experiment 3 The Simple Pendulum Read Taylor chapter 5. (You can skip section 5.6.IV if you aren't comfortable
More informationEfficient method for obtaining parameters of stable pulse in grating compensated dispersion-managed communication systems
3 Conference on Information Sciences and Systems, The Johns Hopkins University, March 12 14, 3 Efficient method for obtainin parameters of stable pulse in ratin compensated dispersion-manaed communication
More informationHomework # 2. SOLUTION - We start writing Newton s second law for x and y components: F x = 0, (1) F y = mg (2) x (t) = 0 v x (t) = v 0x (3)
Physics 411 Homework # Due:..18 Mechanics I 1. A projectile is fired from the oriin of a coordinate system, in the x-y plane (x is the horizontal displacement; y, the vertical with initial velocity v =
More informationAnalysis of visibility level in road lighting using image processing techniques
Scientific Research and Essays Vol. 5 (18), pp. 2779-2785, 18 September, 2010 Available online at http://www.academicjournals.or/sre ISSN 1992-2248 2010 Academic Journals Full enth Research Paper Analysis
More information1 CHAPTER 7 PROJECTILES. 7.1 No Air Resistance
CHAPTER 7 PROJECTILES 7 No Air Resistance We suppose that a particle is projected from a point O at the oriin of a coordinate system, the y-axis bein vertical and the x-axis directed alon the round The
More informationGeneration of random waves in time-dependent extended mild-slope. equations using a source function method
Generation of random waves in time-dependent etended mild-slope equations usin a source function method Gunwoo Kim a, Chanhoon Lee b*, Kyun-Duck Suh c a School of Civil, Urban, and Geosystem Enineerin,
More informationViscoData & ViscoShift
Introductory Theory Manual ViscoData & ViscoShift (Preliminary Version) Contents ) Preface 2) Introduction 3) Dynamic mechanical properties 4) sin ViscoShift and ViscoData 5) Concludin remarks 6) Refereces
More informationInvestigation of ternary systems
Investiation of ternary systems Introduction The three component or ternary systems raise not only interestin theoretical issues, but also have reat practical sinificance, such as metallury, plastic industry
More informationExperiment 1: Simple Pendulum
COMSATS Institute of Information Technoloy, Islamabad Campus PHY-108 : Physics Lab 1 (Mechanics of Particles) Experiment 1: Simple Pendulum A simple pendulum consists of a small object (known as the bob)
More informationDynamical Diffraction
Dynamical versus Kinematical Di raction kinematical theory valid only for very thin crystals Dynamical Diffraction excitation error s 0 primary beam intensity I 0 intensity of other beams consider di racted
More informationIdentification of muscle forces in human lower limbs during sagittal plane movements Part I: Human body modelling
Acta of ioenineerin and iomechanics Vol. No. 00 Identification of muscle forces in human lower limbs durin saittal plane movements Part I: uman bod modellin WOJCIEC LAJER KRZYSZOF DZIEWIECKI ZENON MAZUR
More informationChapter K. Oscillatory Motion. Blinn College - Physics Terry Honan. Interactive Figure
K. - Simple Harmonic Motion Chapter K Oscillatory Motion Blinn Collee - Physics 2425 - Terry Honan The Mass-Sprin System Interactive Fiure Consider a mass slidin without friction on a horizontal surface.
More informationParameterization and Vector Fields
Parameterization and Vector Fields 17.1 Parameterized Curves Curves in 2 and 3-space can be represented by parametric equations. Parametric equations have the form x x(t), y y(t) in the plane and x x(t),
More information11 Free vibrations: one degree of freedom
11 Free vibrations: one deree of freedom 11.1 A uniform riid disk of radius r and mass m rolls without slippin inside a circular track of radius R, as shown in the fiure. The centroidal moment of inertia
More informationMechanics Cycle 3 Chapter 12++ Chapter 12++ Revisit Circular Motion
Chapter 12++ Revisit Circular Motion Revisit: Anular variables Second laws for radial and tanential acceleration Circular motion CM 2 nd aw with F net To-Do: Vertical circular motion in ravity Complete
More informationAorangi (Mt. Cook) we typically view the state of the troposphere on "isobaric surfaces" (surfaces of constant pressure) rather
Ch3 (ctd.) EAS270_Ch3_BehaviourAtmos_B.odp JDW, EAS UAlberta, last mod. 19 Sept. 2016 Atmospheric "behaviour" Aorani (Mt. Cook) "ride" we typically view the state of the troposphere on "isobaric surfaces"
More informationDynamic Analysis of the Load Lifting Mechanisms
Strojniški vestnik - Journal of Mechanical Enineerin 54(2008)10, 655-661 UDC 621.87 Paper received: 08.04.2008 Paper accepted: 07.07.2008 Dynamic Analysis of the Load Liftin Mechanisms Jovan Vladić* -
More informationREVIEW: Going from ONE to TWO Dimensions with Kinematics. Review of one dimension, constant acceleration kinematics. v x (t) = v x0 + a x t
Lecture 5: Projectile motion, uniform circular motion 1 REVIEW: Goin from ONE to TWO Dimensions with Kinematics In Lecture 2, we studied the motion of a particle in just one dimension. The concepts of
More informationFiring an Ideal Projectile
92 Chapter 13: Vector-Valued Functions and Motion in Space 13.2 Modelin Projectile Motion 921 r at time t v v cos i a j (a) v sin j Newton s second law of motion sas that the force actin on the projectile
More informationINFLUENCE OF TUBE BUNDLE GEOMETRY ON HEAT TRANSFER TO FOAM FLOW
HEFAT7 5 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics Sun City, South Africa Paper number: GJ1 INFLUENCE OF TUBE BUNDLE GEOMETRY ON HEAT TRANSFER TO FOAM FLOW Gylys
More informationModified Reductive Perturbation Method as Applied to Long Water-Waves:The Korteweg-de Vries Hierarchy
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.6(28) No.,pp.-2 Modified Reductive Perturbation Method as Applied to Lon Water-Waves:The Kortewe-de Vries Hierarchy
More informationExam 2A Solution. 1. A baseball is thrown vertically upward and feels no air resistance. As it is rising
Exam 2A Solution 1. A baseball is thrown vertically upward and feels no air resistance. As it is risin Solution: Possible answers: A) both its momentum and its mechanical enery are conserved - incorrect.
More informationarxiv: v1 [quant-ph] 30 Jan 2015
Thomas Fermi approximation and lare-n quantum mechanics Sukla Pal a,, Jayanta K. Bhattacharjee b a Department of Theoretical Physics, S.N.Bose National Centre For Basic Sciences, JD-Block, Sector-III,
More informationFirst- and Second Order Phase Transitions in the Holstein- Hubbard Model
Europhysics Letters PREPRINT First- and Second Order Phase Transitions in the Holstein- Hubbard Model W. Koller 1, D. Meyer 1,Y.Ōno 2 and A. C. Hewson 1 1 Department of Mathematics, Imperial Collee, London
More informationDirection of Hurricane Beta Drift in Horizontally Sheared Flows
146 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 54 Direction of Hurricane Beta Drift in Horizontally Sheared Flows BIN WANG, XIAOFAN LI,* AND LIGUANG WU Department of Meteoroloy, School of Ocean and Earth
More informationEquivalent rocking systems: Fundamental rocking parameters
Equivalent rockin systems: Fundamental rockin parameters M.J. DeJon University of Cambride, United Kindom E.G. Dimitrakopoulos The Hon Kon University of Science and Technoloy SUMMARY Early analytical investiations
More informationSlip-Flow and Heat Transfer in Isoflux Rectangular Microchannels with Thermal Creep Effects
Journal of Applied Fluid Mechanics, Vol. 3, No. 2, pp. 33-4, 200. Available online at www.jafmonline.net, ISSN 735-3645. Slip-Flow and Heat Transfer in Isoflux Rectanular Microchannels with Thermal Creep
More informationConical Pendulum: Part 2 A Detailed Theoretical and Computational Analysis of the Period, Tension and Centripetal Forces
European J of Physics Education Volume 8 Issue 1 1309-70 Dean onical Pendulum: Part A Detailed heoretical and omputational Analysis of the Period, ension and entripetal orces Kevin Dean Physics Department,
More informationMidterm Feb. 17, 2009 Physics 110B Secret No.=
Midterm Feb. 17, 29 Physics 11B Secret No.= PROBLEM (1) (4 points) The radient operator = x i ê i transforms like a vector. Use ɛ ijk to prove that if B( r) = A( r), then B( r) =. B i = x i x i = x j =
More informationA-LEVEL Mathematics. MM05 Mechanics 5 Mark scheme June Version 1.0: Final
A-LEVEL Mathematics MM05 Mechanics 5 Mark scheme 6360 June 016 Version 1.0: Final Mark schemes are prepared by the Lead Assessment Writer and considered, toether with the relevant questions, by a panel
More information5 Shallow water Q-G theory.
5 Shallow water Q-G theory. So far we have discussed the fact that lare scale motions in the extra-tropical atmosphere are close to eostrophic balance i.e. the Rossby number is small. We have examined
More informationDynamics 4600:203 Homework 03 Due: February 08, 2008 Name:
Dynamics 4600:03 Homework 03 Due: ebruary 08, 008 Name: Please denote your answers clearly, i.e., bo in, star, etc., and write neatly. There are no points for small, messy, unreadable work... please use
More informationVector Valued Functions
SUGGESTED REFERENCE MATERIAL: Vector Valued Functions As you work throuh the problems listed below, you should reference Chapters. &. of the recommended textbook (or the equivalent chapter in your alternative
More informationAutonomous Orbit Coordination for Two Unmanned Aerial Vehicles
Autonomous Orbit Coordination for Two Unmanned Aerial Vehicles Rolf Rysdyk, Christopher Lum and Juris Vaners Autonomous Fliht Systems Laboratory University of Washinton, Seattle, WA 98195 ABSTRACT This
More informationNumerical and Experimental Investigations of Lateral Cantilever Shaft Vibration of Passive-Pitch Vertical-Axis Ocean Current
R. Hantoro, et al. / International Enery Journal 1 (011) 191-00 191 Numerical and Experimental Investiations of Lateral Cantilever Shaft Vibration of Passive-Pitch Vertical-Axis Ocean Current R. Hantoro
More informationAccording to the Ratio of Height to Least Lateral Dimension: (ACI Code and ACI Code ) < 22, when member is braced against sideways.
CHAPTER 5 COLUMN 5.1 INTRODUCTION Columns are structural elements used primarily to support compressive loads. All practical columns are members subjected not only to axial load but also to bendin moments,
More information2.5 Velocity and Acceleration
82 CHAPTER 2. VECTOR FUNCTIONS 2.5 Velocity and Acceleration In this section, we study the motion of an object alon a space curve. In other words, as the object moves with time, its trajectory follows
More information(a) Find the function that describes the fraction of light bulbs failing by time t e (0.1)x dx = [ e (0.1)x ] t 0 = 1 e (0.1)t.
1 M 13-Lecture March 8, 216 Contents: 1) Differential Equations 2) Unlimited Population Growth 3) Terminal velocity and stea states Voluntary Quiz: The probability density function of a liht bulb failin
More informationUniversity of Alabama Department of Physics and Astronomy. PH 125 / LeClair Fall Exam III Solution
University of Alabama Department of Physics and Astronomy PH 5 / LeClair Fall 07 Exam III Solution. A child throws a ball with an initial speed of 8.00 m/s at an anle of 40.0 above the horizontal. The
More informationCh 14: Feedback Control systems
Ch 4: Feedback Control systems Part IV A is concerned with sinle loop control The followin topics are covered in chapter 4: The concept of feedback control Block diaram development Classical feedback controllers
More informationProblem Set 2 Solutions
UNIVERSITY OF ALABAMA Department of Physics and Astronomy PH 125 / LeClair Sprin 2009 Problem Set 2 Solutions The followin three problems are due 20 January 2009 at the beinnin of class. 1. (H,R,&W 4.39)
More informationV DD. M 1 M 2 V i2. V o2 R 1 R 2 C C
UNVERSTY OF CALFORNA Collee of Enineerin Department of Electrical Enineerin and Computer Sciences E. Alon Homework #3 Solutions EECS 40 P. Nuzzo Use the EECS40 90nm CMOS process in all home works and projects
More informationThe Hashin-Shtrikman-Walpole lower and upper bounds for d-dimensional isotropic materials with n phases are the following (1): i=1
Supplemental Document Mechanics of Partially Mineralized Collaen Yanxin Liu, Stavros Thomopoulos, Chanqin Chen, Victor Birman, Markus J. Buehler, & Guy M. Genin 1 Hashin-Shtrikman-Walpole bounds The Hashin-Shtrikman-Walpole
More informationFoundation on Compressible Fluid Flow
Chapter Foundation on Compressible Fluid Flow. COMRESSIBLE FLUIDS In everyday life one reconizes three states of matter, namely, solid, liquid and aseous. Solids, liquids and ases are all comprised of
More informationDielectric characteristics of glass fibre reinforced plastics and their components
Plasticheskie Massy, No. 11, 2004, pp. 20 22 Dielectric characteristics of lass fibre reinforced plastics and their components V. I. Sokolov, S. I. Shalunov, I. G. Gurtovnik, L. G. Mikheeva, and I. D.
More information4.3. Solving Friction Problems. Static Friction Problems. Tutorial 1 Static Friction Acting on Several Objects. Sample Problem 1.
Solvin Friction Problems Sometimes friction is desirable and we want to increase the coefficient of friction to help keep objects at rest. For example, a runnin shoe is typically desined to have a lare
More informationSTOCHASTICALLY GENERATED MULTIGROUP DIFFUSION COEFFICIENTS
STOCHASTICALLY GENERATED MULTIGROUP DIFFUSION COEFFICIENTS A Thesis Presented to The Academic Faculty by Justin M. Pounders In Partial Fulfillment of the Requirements for the Deree Master of Science in
More informationDynamic Response of Cantilever Retaining Walls Considering Soil Non-Linearity
Missouri University of Science and Technoloy Scholars' Mine International Conferences on Recent Advances in Geotechnical Earthquake Enineerin and Soil Dynamics 2 - Fifth International Conference on Recent
More informationANALYZE In all three cases (a) (c), the reading on the scale is. w = mg = (11.0 kg) (9.8 m/s 2 ) = 108 N.
Chapter 5 1. We are only concerned with horizontal forces in this problem (ravity plays no direct role). We take East as the +x direction and North as +y. This calculation is efficiently implemented on
More informationPrince Sultan University Physics Department First Semester 2012 /2013. PHY 105 First Major Exam Allowed Time: 60 min
Prince Sultan University Physics Department First Semester 01 /01 PHY 105 First Major Exam Allowed Time: 60 min Student Name: 1. Write your name in the specified space NOW.. Any paper without name will
More informationFig.1: Non-stationary temperature distribution in a circular plate.
LECTURE Problem 1 13. 5 : Fi.1: Non-stationary temperature distribution in a circular plate. The non-stationary, radial symmetric temperature distribution ur, t in a circular plate of radius c is determined
More informationEmotional Optimized Design of Electro-hydraulic Actuators
Sensors & Transducers, Vol. 77, Issue 8, Auust, pp. 9-9 Sensors & Transducers by IFSA Publishin, S. L. http://www.sensorsportal.com Emotional Optimized Desin of Electro-hydraulic Actuators Shi Boqian,
More informationDigital Electronics Paper-EE-204-F SECTION-A
B.Tech 4 th Semester (AEIE) F Scheme, May 24 Diital Electronics Paper-EE-24-F Note : Attempt five questions. Question is compulsory and one question from each of the four sections.. What is a loic ate?
More informationf 1. (8.1.1) This means that SI unit for frequency is going to be s 1 also known as Hertz d1hz
ecture 8-1 Oscillations 1. Oscillations Simple Harmonic Motion So far we have considered two basic types of motion: translational motion and rotational motion. But these are not the only types of motion
More informationStoichiometry of the reaction of sodium carbonate with hydrochloric acid
Stoichiometry of the reaction of sodium carbonate with hydrochloric acid Purpose: To calculate the theoretical (expected) yield of product in a reaction. To weih the actual (experimental) mass of product
More informationPhysics 20 Homework 1 SIMS 2016
Physics 20 Homework 1 SIMS 2016 Due: Wednesday, Auust 17 th Problem 1 The idea of this problem is to et some practice in approachin a situation where you miht not initially know how to proceed, and need
More informationA-Level Mathematics. MM05 Mechanics 5 Final Mark scheme June Version/Stage: v1.0
A-Level Mathematics MM0 Mechanics Final Mark scheme 6360 June 07 Version/Stae: v.0 Mark schemes are prepared by the Lead Assessment Writer and considered, toether with the relevant questions, by a panel
More informationPhysics 111 P 2 A = P 1. A + mg = P 1. A + ρ( AΔh)g. Wednesday, 8-9 pm in NSC 118/119 Sunday, 6:30-8 pm in CCLIR 468.
ics Announcements day, ember 11, 011 C5: Fluids Pascal s Principle Archimede s Principle Fluid Flows Continuity Equation Bernoulli s Equation Toricelli s Theorem Wednesday, 8-9 pm in NSC 118/119 Sunday,
More information1. According to the U.S. DOE, each year humans are. of CO through the combustion of fossil fuels (oil,
1. Accordin to the U.S. DOE, each year humans are producin 8 billion metric tons (1 metric ton = 00 lbs) of CO throuh the combustion of fossil fuels (oil, coal and natural as). . Since the late 1950's
More informationSHIP BUOYANCY AND STABILITY. Lecture 02 Ship equilibrium and introduction to ship hydrostatics
SHIP BUOYANCY AND STABILITY Lecture 02 Ship equilibrium and introduction to ship hydrostatics 1 Literature J. Matusiak: Laivan kelluvuus ja vakavuus Biran A. B., Ship Hydrostatics and Stability, 2003 J.
More informationPhysicsAndMathsTutor.com
. A particle of mass m is projected vertically upwards, at time t =, with speed. The particle is mv subject to air resistance of manitude, where v is the speed of the particle at time t and is a positive
More informationGet Solution of These Packages & Learn by Video Tutorials on PROJECTILE MOTION
FREE Download Study Packae from website: www.tekoclasses.com & www.mathsbysuha.com Get Solution of These Packaes & Learn by Video Tutorials on www.mathsbysuha.com. BASIC CONCEPT :. PROJECTILE PROJECTILE
More informationA Study of Computer Modeling Techniques to Predict the Response of Floor Systems Due to Walking
A Study of Computer Modelin Techniques to Predict the Response of Floor Systems Due to Walkin by Jason Daniel Perry Thesis submitted to the faculty of the Virinia Polytechnic Institute and State University
More informationMODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION
= = = MODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION Mauro Cettolo and Marcello Federico ITC-irst - Centro per la Ricerca Scientifica e Tecnoloica I-385 Povo, Trento, Italy ABSTRACT Robust acoustic
More informationCourse , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Energy Cycle
Course.8, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Enery Cycle Enery Forms: As we saw in our discussion of the heat budet, the enery content of the atmosphere per
More informationUSE OF FILTERED SMITH PREDICTOR IN DMC
Proceedins of the th Mediterranean Conference on Control and Automation - MED22 Lisbon, Portual, July 9-2, 22. USE OF FILTERED SMITH PREDICTOR IN DMC C. Ramos, M. Martínez, X. Blasco, J.M. Herrero Predictive
More informationRESISTANCE STRAIN GAGES FILLAMENTS EFFECT
RESISTANCE STRAIN GAGES FILLAMENTS EFFECT Nashwan T. Younis, Younis@enr.ipfw.edu Department of Mechanical Enineerin, Indiana University-Purdue University Fort Wayne, USA Bonsu Kan, kan@enr.ipfw.edu Department
More informationCode_Aster. Element CABLE_GAINE
Titre : Élément CABLE_GAINE Date : 28/07/2015 Pae : 1/12 Element CABLE_GAINE Summary: The element CABLE_GAINE presented in this document is to model cables of prestressin bein able to rub or slip into
More informationCOUPLED OSCILLATORS. Two identical pendulums
COUPED OSCIATORS A real physical object can be rearded as a lare nuber of siple oscillators coupled toether (atos and olecules in solids. The question is: how does the couplin affect the behavior of each
More informationThree Spreadsheet Models Of A Simple Pendulum
Spreadsheets in Education (ejsie) Volume 3 Issue 1 Article 5 10-31-008 Three Spreadsheet Models Of A Simple Pendulum Jan Benacka Constantine the Philosopher University, Nitra, Slovakia, jbenacka@ukf.sk
More informationThermo-Mechanical Damage Modeling of Polymer Matrix Composite Structures in Fire
Thermo-Mechanical Damae Modelin of Polymer Matrix Composite Structures in Fire CHANGSONG LUO, and JIM LUA Global Enineerin and Materials, Inc. One Airport Place, Suite One Princeton, NJ 08540 USA ABSTRACT
More informationthe equations for the motion of the particle are written as
Dynamics 4600:203 Homework 02 Due: ebruary 01, 2008 Name: Please denote your answers clearly, ie, box in, star, etc, and write neatly There are no points for small, messy, unreadable work please use lots
More information7.2 Maximization of the Range of a Rocket
138 CHAPTER 7. SOME APPLICATIONS The counterintuitive answer that a supersonic aircraft must dive first in order to climb to a iven altitude in minimum time was first discovered by Walter Denham and Art
More informationLightning Transients on Branched Distribution Lines Considering Frequency-Dependent Ground Parameters
214 International Conference on Lihtnin Protection (ICLP), Shanhai, China Lihtnin Transients on Branched Distribution Lines Considerin Frequency-Dependent Ground Parameters Alberto De Conti LRC Lihtnin
More informationAnother possibility is a rotation or reflection, represented by a matrix M.
1 Chapter 25: Planar defects Planar defects: orientation and types Crystalline films often contain internal, 2-D interfaces separatin two reions transformed with respect to one another, but with, otherwise,
More informationDESIGN OF PI CONTROLLERS FOR MIMO SYSTEM WITH DECOUPLERS
Int. J. Chem. Sci.: 4(3), 6, 598-6 ISSN 97-768X www.sadurupublications.com DESIGN OF PI CONTROLLERS FOR MIMO SYSTEM WITH DECOUPLERS A. ILAKKIYA, S. DIVYA and M. MANNIMOZHI * School of Electrical Enineerin,
More informationConical Pendulum Linearization Analyses
European J of Physics Education Volume 7 Issue 3 309-70 Dean et al. Conical Pendulum inearization Analyses Kevin Dean Jyothi Mathew Physics Department he Petroleum Institute Abu Dhabi, PO Box 533 United
More informationAmplitude Adaptive ASDM without Envelope Encoding
Amplitude Adaptive ASDM without Envelope Encodin Kaspars Ozols and Rolands Shavelis Institute of Electronics and Computer Science 14 Dzerbenes Str., LV-16, Ria, Latvia Email: kaspars.ozols@edi.lv, shavelis@edi.lv
More information