APPENDIX 5.5.D CHARACTERIZATION OF WIND LOADING OF TELESCOPES

Size: px
Start display at page:

Download "APPENDIX 5.5.D CHARACTERIZATION OF WIND LOADING OF TELESCOPES"

Transcription

1 APPENDIX 5.5.D CHARACTERIZATION OF WIND LOADING OF TELESCOPES Published in SPIE Proceedings Vol. 4757, Integrated Modeling of Telescoes, Lund, Sweden, February, 7-8.

2 Characterization of Wind Loading of Telescoes George Z. Angeli *, Myung K. Cho, Michael Sheehan, Larry M. Ste New Initiatives Office, AURA Inc., 95 North Cherry Ave., Tucson AZ 8579 Gemini Observatories, 95 North Cherry Ave., Tucson AZ 8579 ABSTRACT Ground-based telescoes oerate in a turbulent atmoshere that affects the otical ath across the aerture by changing both the mirror ositions (wind induced vibrations) and the air refraction index. Although the characteristics of the atmoshere are well understood in the inertial range, the validity of the homogeneous, isotroic field assumtion is questionable inside the enclosure and in the close vicinity of the structure. To understand the effect of wind on an actual telescoe, we conducted extensive wind measurements at the Gemini South Telescoe. Simultaneous measurements were made of ressures at multile oints on the mirror surface, as well as wind velocity and direction at several locations inside and outside the dome. During the test we varied the dome osition relative to the wind, the telescoe elevation angle, the osition of windscreens in the observing slit, and the size of the oenings in the ventilation gates. The data sets have been rocessed to rovide the temoral and satial characteristics of the ressure variations on the rimary mirror in comarison to the theory of atmosheric turbulence. Our investigation is art of an effort leading to the develoment of a scalable wind model for large telescoe simulations, which describes the forces due to air turbulence on the rimary mirror and telescoe structure reasonably well even inside an enclosure. Keywords: wind buffeting, wind model, large telescoes. INTRODUCTION Considering the effects of wind in the design of the next generation giant ground based telescoes is going to be vital to ensure the success of these telescoes. As the sizes of the telescoes are getting larger and larger, the concerns of wind are getting more and more serious. The direct effect of wind is buffeting, i.e. wind inflicted deformations of the mechanical structure as well as the otical elements, like the rimary mirror. The rimary mirror deformation is esecially roblematic for large segmented mirrors, where the quality of the otical surface deends only on the stiffness of the back structure and actuators. The severe effect of wind on giant telescoes is aggravated by the fact that the larger the telescoe, the lower its resonant frequencies are. Consequently, a larger telescoe can absorb a significantly higher ercentage of wind energy than a smaller one could. The same wind also affects local seeing at the telescoe through the thermal characteristics of the local turbulence. Extensive studies have roved, that maintaining thermal equilibrium inside the dome can drastically reduce this local seeing. The straightforward way to achieve isothermal surfaces is to allow efficient ventilation of the telescoe enclosure. However, the need for ventilation must be balanced against the buffeting effects of the wind. The behavior of the turbulent atmosheric boundary layer reresenting the wind close to the Earth s surface - is reasonably well understood. The velocity fluctuations at a given oint of this turbulent layer may be considered as a suerosition of theoretical eddies transorted by the mean wind, like a standing wave can be concetualized as a suerosition of traveling waves. Each eddy is assumed to cause a velocity fluctuation at the given oint with frequency of f. By analogy to the traveling wave, a wavelength and wave number can be associated with each eddy. U πf ω λ = ; κ = = f U U Here U is the mean velocity of the wind. The hyothesis that the frozen eddies are traveling with the wind, i.e. the satial frequency (wavenumber) of turbulence κ is roortional to its temoral frequency f, was first suggested by Taylor (in []).

3 The internal structure of the turbulence layer can be understood as an energy cascading rocess. When the wind seed u or more recisely the Reynolds number (Lu/ν) characterizing the airflow - reaches a critical value, the flow becomes turbulent and large eddies aear. However, inside each eddy the seed can still be high enough to generate internal turbulence, i.e. smaller eddies. This cascading goes on until the eddy size L becomes small enough to revent further turbulence. In the laminar flow then the remaining energy is dissiated as heat. Kolmogorov s two hyotheses ([] in []) state that there is a range in size of eddies where the fluid motion is locally isotroic and deends only on the rate of energy dissiation ε in the smallest eddies. The boundaries of this inertial subrange are the inner and outer scales (l and L, resectively). π π < κ < L l In order to characterize the local structure of turbulence in the atmoshere, Kolmogorov [] in [] introduced the structure tensor. After some mathematical and hysical considerations ([] in []) [], [4], [5], the structure tensor yields the onedimensional structure function, which deends on the scalar distance r and energy transfer rate ε only, with a dimensionless scaling factor of C u. () r = C ε r D u u The structure function has a close relationshi to the auto-correlation B u () and cross-correlation B u (r) functions. D () r = [ B () B () r ] () u Bernoulli s Law imlies, that - in absence of viscous stresses, - the stochastic roerties of the ressure in a given fluid are analogous to those of the square of the velocity. According to Kolmogorov s second hyothesis, in the inertial subrange we can neglect the effect of viscosity. Invoking the same dimensional argument that resulted in the velocity structure function [4], we must consider the dimensions of the velocity square structure function, dislacement and energy transfer rate as (length) 4 /(time) 4, (length) and (length) /(time), resectively. The resulting structure function has a different scaling factor C and steeer distance deendence than the velocity structure function. D u 4 4 () r = C ε r A characteristic length of correlation L can be defined for the turbulence by means of the auto-correlation B () and crosscorrelation B (r) functions, as an integral scale [6]. L = B ( ) Actually the integral can be limited to a bounded region R as long as the cross-correlation is negligible outside of this region. By using the well-known relationshi between the structure function and the correlation functions, one can exress the integral scale with the structure function only. L B ( R) u ()r r d R = R D ()r r d () D According to the Wiener-Khinchine theorem, the ower sectral density (PSD) Φ (κ) of ressure fluctuation is the Fourier transforms of the cross-correlation function B (ρ). Following Tatarski s calculation for velocity ower sectral density as an examle [4], one can determine the ressure ower sectral density in the inertial subrange. Φ 7 ( κ) = C ε κ However, the ower content of eddies certainly does not aroach infinity with growing eddy sizes, as the Kolmogorov sectrum imlies. Based on exerimental results, von Karman suggested ower saturation outside the inertial subrange ([7] in []). In his formula, the characteristic wave number κ is related to the outer scale of the turbulence, as κ =π/l.

4 vk C vk Φ ( κ ) = 7 6 κ + κ The temoral ower sectral density of the turbulence can be calculated by using Taylor s hyothesis. Here f, the characteristic frequency equals to U/L. vk C vk t Φ t ( f ) = () 7 6 f + f There are several other fitting functions for wind ower sectral densities, like Davenort [8], Antoniou [9], Harris [6], etc. but in this work we use the von Karman sectrum for its convenience. Although it is a relatively simle function, it agrees reasonably well with the test results. The von Karman sectrum is fully defined by two arameters, the magnitude C vk t and the bandwidth f. Although it is customary to assume a Kolmogorov (or von Karman) sectrum during the design and simulation of telescoes, it is widely considered as a crude estimate only. The turbulence inside the dome is obviously not homogeneous and isotroic enough to suort Kolmogorov s hyotheses, certainly not on the scale of the mountainto environment. To imrove the aroximation, a frequency deendent multilying factor, the aerodynamic attenuation was introduced [], []. It is worth to note that eole usually choose f and f a that close to each other. A = (4) + f f a The integration of the Gemini South telescoe rovided an excellent oortunity to conduct extensive wind velocity and ressure measurements in an existing large telescoe dome. Besides actual engineering interests like determining the resonse of Gemini structure on wind buffeting, or develoing a rocedure to set the wind gates on that articular enclosure, - our objective was to verify and ossibly imrove our current understanding of wind loading of telescoes, as well as to work toward a scalable telescoe wind model. Some of the results of these measurements were reorted earlier [], [], [4]. Our articular interest in our investigation resulted in this aer was to find out whether the wind buffeting of the rimary mirror was mainly influenced by the outside environment, the enclosure, or local features on and around the mirror. By identifying the source we may get a ste closer to actually redicting at least a worst-case wind scenario for telescoe designs and simulations. 4. MEASUREMENT SETUP AND DATA REDUCTION Although during the measurement camaign at Gemini South Telescoe a large amount of wind velocity and ressure data was collected, in this work we mainly focus on the ressure distribution measurements on the rimary mirror. The details of the data collection were described by Cho et al. []. Here we reeat only the most relevant information necessary to understand our conclusions. The measurements were taken inside the Gemini dome, in an environment where the wind was controlled by large wind gates (Figure.). At the time of the measurements the glass rimary mirror was not installed yet, but there was a steel dummy instead in lace. It was covered by lywood to rovide a continuous surface. The ressure sensors were installed on this lywood layer according to the configuration shown in Figure. The x-axis of the Cartesian coordinate system is arallel to the elevation axis of the telescoe with y-axis ointing uward when the telescoe is ointed toward the horizon. The actual 4

5 environment of the rimary mirror is shown in Figure. The stray light baffle is clearly noticeable in the icture, together with the mirror cover. Figure. Gemini enclosure with large ventilation gates to allow natural ventilation with ambient air Data was taken with 6 bit resolution at Hz samling rate. The total length of each measurement was 5 minutes that resulted in samling oints. The detection limit (S/N ratio ) of the ressure sensors is 8 mpa with full range of Pa. Figure. The Gemini rimary mirror and its environment 5

6 Velocity sensors (three-dimensional ultrasonic anemometers) were also laced at different locations: (i) on the edges of the rimary mirror at +X, -X, +Y and Y ositions; (ii) at the secondary mirror; and (iii) on to of the dome. The dome sensor detected the outside wind velocity for reference uroses. In the cases considered in this aer, the ambient wind velocity was ranging from 8m/s to m/s. The data sets collected have a unique label, for examle c96oo. The first letter shows the date it was collected, the first three numbers the wind direction and the last two numbers the telescoe zenith angle. The tailing two letters indicates the osition of ventilation gates: in the examle both sides are oen. In all cases described in this aer the windshield and shutter are both oen leaving the observing slit unobstructed. The temoral behavior of the turbulence on the rimary mirror is characterized by the ower sectral densities of the time series collected from each sensor. To estimate the ower sectral densities we used Welch s method [5]. This technique divides the entire time series in sections with length of 56 samles and 5 % overla, each section is windowed with a Hamming window, and indeendent eriodograms are calculated for each. The final estimate is the average of these eriodograms X Y 4 Figure. Pressure sensor locations on the dummy rimary mirror To characterize the satial correlations in the samles, we have used structure functions. Each structure function was directly calculated following the definition by grouing together the sensor airs with a given distance. Although from a strictly mathematical oint of view the structure function is defined on homogenous and isotroic fields, we alied its definition on the turbulent field inside the dome to visualize the analogy and also the differences between inside and outside ressure fields. It should be understood that the structure function averaged over the whole mirror surface is only a global measure of correlation and conceals highly turbulent local effects. The mean ressure value on the rimary mirror is certainly not constant as it would be exected outside of the enclosure and is rather deterministic with significant cross-correlation over the whole area. Since our objective was to investigate the random attributes of turbulence, we excluded the mean values from our calculations. Earlier [], when we were mainly interested in the deformation of the mirror, we reorted the structure functions incororating the mean values.. RESULTS Looking at the ressure ower sectral densities, one can recognize evident differences between high and low wind conditions (Figures 4. and 8.). The characteristic differences are clearly recurring, so it s reasonable to discuss the features of the ower sectral densities individually for oen and closed ventilation gates... Power sectral densities for high winds When both the windshield and the ventilation gates are oen, there is significant wind inside the enclosure. While the mean wind velocity outside the dome is fairly stable, on the edges of the mirror surface it varies widely, deending on the wind direction relative to the dome as well as the zenith angle of the telescoe. In the general case, when the wind is blowing on the mirror from an arbitrary direction, the flow field is quite comlex and it is rather difficult to identify common atterns. In other words, the local effects of the mirror and mirror suort system are seemingly covering u the background turbulence. 6

7 Since the detailed velocity distribution on the mirror surface is not available, these local erturbations cannot be extracted from the data to better understand the global (background) effects. We calculated the ressure ower sectral densities for all of the sensors in numerous measurement series detected under various conditions like wind direction and mean velocity, telescoe elevation and ventilation gate osition. For all tests carried out under high wind conditions, i.e. with oen gates, the PSD fits surrisingly well to the theoretical von Karman curve (Equation.). The roll-off sloe corresonds to the (-7/) Kolmogorov law with varying bandwidth (f ) at different areas and under different test conditions (Figure 4.). To recognize ossible atterns, let us consider a relatively simle case, when the telescoe is ointing to zenith, i.e. the rimary mirror is horizontal, and the wind blows from (-Y) direction (Figures 5. and 6.). In this configuration we can assume a wind velocity and direction that is fairly constant over the mirror surface. Indeed, the anemometers at (+X) and (-X) osition detected similar values (. m/s and.6 m/s) with similar (close to +Y) directions. Even for this setu, the mean ressure on the mirror is far from being constant. However, the RMS ressure lot indicates a noticeable searation between local effects and a more evenly distributed background. It is remarkable that the RMS ressure calculated without considering the mean values is significantly higher than the mean ressure itself, which imlies high turbulence. PSD [Pa /Hz] - PSD [Pa /Hz] measurement von Karman fit - - Frequency [Hz] - measurement von Karman fit - - Frequency [Hz] Pressure sensor # Pressure sensor #9 Figure 4. Pressure ower sectral densities on the rimary mirror with von Karman fit (Data set coo). The turbulence on the mirror can be characterized with the arameters of the curve fit corresonding to Equation., i.e. with the magnitude C t vk and bandwidth f. Indeed, the arameter mas show high energy content without significant turbulence (bandwidth) in the negative ressure area (sensors#8, #4 and #9), which robably indicates a local maximum in tangential wind velocity. On the other hand, there is a highly turbulent area right behind the baffle column, as it aears on the bandwidth lot. Our assumtion is that a characteristic length L - virtually an outer scale - can be derived from the bandwidth f of the background turbulence. This length calculated from the bandwidth average of the background area with the assumtion of.5 m/s wind seed is about meter. The characteristic length corresonding to the highly turbulent sot behind the baffle is 7. meter. In another case we investigated the telescoe zenith angle was 6 o with wind blowing sideways, (-X) direction. The mean wind flow was arallel to the mirror surface and reasonably constant, as the anemometers indicated. The gain and bandwidth attern on the mirror is fairly similar to the horizontal case, excet the bandwidth values are significantly higher (Figure 7.) and the attern is rotated due to the different wind direction. The rotation is less than the exected 9 o because the wind is not comletely horizontal; it has a slight elevation angle. The bandwidth values are just artly higher due to larger mean velocities (8.5 m/s and 8. m/s at +X and X ositions, consecutively). The characteristic length of the turbulence is significantly smaller: the background turbulence corresonds to about meter. The high turbulent area behind the baffle yields a characteristic length of.8 meter. 7

8 Mean ressure (ranging from to Pa) RMS ressure (ranging from to Pa) Figure 5. Pressure distribution on the rimary mirror ointing to zenith with wind direction of 9 o (+Y) (Data set coo) Y Y X X Magnitude (ranging from to Pa /Hz) Bandwidth (ranging from. to.5 Hz) Figure 6. Parameters of the best fit for ressure ower sectral densities on the rimary mirror ointing to zenith with wind direction of 9 o (+Y) (Data set coo) Y Y X X Magnitude (ranging from to 9 Pa /Hz) Bandwidth (ranging from.6 to Hz) Figure 7. Parameters of the best fit for ressure ower sectral densities on the rimary mirror with zenith angle of 6 o and wind direction of aroximately 8 o (-X) (Data set c96oo)

9 Table. Characteristic length of background turbulence on rimary mirror L U = f mean Data set (coo) Data set (c96oo) Data set (too) Data set 4 (d96oo) ZENITH ANGLE º WIND DIRECTION (+Y) ZENITH ANGLE 6º WIND DIRECTION (-X) L U mean L U mean. meter ~.5 m/s.7 meter ~ 5. m/s.4 meter ~ 8. m/s. meter ~. m/s Other data sets collected on different days with different mean wind velocities but under similar telescoe conditions yield similar numbers (Table.)... Power sectral densities for low winds The data collected with closed ventilation gates show some anomaly that needs further investigation. Many sensors detected data with PSD fitting on the Kolmogorov von Karman curves (Figure 8., sensor #). However, a significantly high number of sensors ractically the whole right side of the mirror collected data with much steeer roll-off and a tail above Hz (Figure 8., sensor #4). The figure shows a curve fit with von Karman sectrum corrected with the aerodynamic attenuation (Equation 4.). As it is clearly visible in the figure, the fit is not erfect; it estimates fairly well the increased sloe, but it does not reflect two major characteristics of the measured curve: (i) the hum around. Hz and (ii) the tail above Hz. - - PSD [Pa /Hz] - - PSD [Pa /Hz] measurement von Karman fit - - Frequency [Hz] -5-6 curve fit (von Karman & aerodynamic attenuation) measurement - - Frequency [Hz] Pressure sensor # Pressure sensor #4 Figure 8. Pressure ower sectral densities on the zenith ointing rimary mirror with week wind (Data set ccc) It is worth to note that the wind seed inside the enclosure was extremely low for this measurement. With the vent gates closed on both sides, the mean wind velocity above the mirror was.5 - m/s and the flow field was rather inhomogeneous and comlicated... Structure function The satial characteristics of the turbulence can be described with the structure function. To minimize the random effects, the temoral average of the instantaneous structure functions was calculated for the whole 5 minutes collection time. We chose to resent the square root of the structure function instead of the structure function itself (Figures 9,, and ), because the 9

10 former has a saturation level with the dimension of satial RMS ressure. Actually, the saturation level is than the temoral average of the satial RMS ressures on the mirror. times higher It is obvious from Equation. that the structure function saturates as the cross-correlation diminishes with increasing the distance between sensors. It is also obvious from the following definition of the structure function that with zero searation it goes to zero. D () r = [ ( r + r) ( )] r The correlation length defined in Equation. is clearly visible on the lots as rising distance of the structure function. The correlation lengths for the different data sets are collected in Table. Quite interestingly, the correlation length does not deend on the elevation angle of the telescoe, but it shows strong deendence on the azimuth angle. Since the mirror even with the surrounding telescoe structure seems to be fairly symmetric, it suggests that the correlation length is mainly defined by the environment of the mirror. satial.5.5 Structure function, D [Pa].5.5 Structure function, D [Pa] Sensor sacing, d [m] 7 8 Sensor sacing, d [m] Zenith angle o (data set coo) Zenith angle 6 o (data set c6oo) Figure 9. Pressure structure function on the mirror surface with wind direction of 9 o (+Y) Structure function, D [Pa] Structure function, D [Pa] Sensor sacing, d [m] 7 8 Sensor sacing, d [m] Zenith angle o (data set c45oo) Zenith angle 6 o (data set c456oo) Figure.. Pressure structure function on the mirror surface with wind direction of 5 o

11 .5 Structure function, D [Pa].5.5 Structure function, D [Pa] Sensor sacing, d [m] 7 8 Sensor sacing, d [m] Zenith angle o (data set c9oo) Zenith angle 6 o (data set c96oo) Figure. Pressure structure function on the mirror surface with wind direction of 8 o (-X) According to Simiu [6], the ratio of the characteristic and correlation lengths in turbulence is constant. L = 7.8 L Where it was ossible we derived the correlation lengths of data sets from their characteristic lengths. The results are retty close to, but slightly higher than the measured values. It is understandable though, considering that the measured structure functions contain the local turbulences that are obviously reducing the correlation, while the bandwidths were calculated only for the background. Table. Correlation length L of turbulence on the rimary mirror WIND DIRECTION 9º WIND DIRECTION 5º WIND DIRECTION 8º L data set L data set L data set o.95 meter coo. meter c45oo.8 meter c9oo zenith angle.8 meter too.7 meter t45oo o.78 meter c6oo.8 meter c456oo.9 meter c96oo zenith angle 6.88 meter d6oo.9 meter d96oo 4. CONCLUSIONS It is reasonable to conclude that the background turbulence detected on the mirror is defined by the environment of the mirror, while the local obstacles are causing turbulent effects rather limited in sace. Both the ower sectral densities and the structure functions indicate that as the wind is moving sideways i.e. its azimuth angle increasing, the airflow above the mirror is getting more turbulent; the bandwidth is increasing while the correlation length is decreasing. Although, as it was reorted earlier [], the overall RMS ressure on the mirror is higher if the mirror is ointing into the wind, the temoral bandwidth of the wind load is rather small. On the other hand, the side wind can cause considerable roblems not because of its intensity but because of its bandwidth. For a telescoe of the size of Gemini this wind buffeting frequency is still way below the resonances of the structure, but as telescoes are getting larger, the interaction is more likely. For single iece rimary mirrors like the Gemini mirror - the correlation length of ressure variations is not a significant feature, since any asymmetric load excites mainly the lowest order structural mode of the mirror, which is astigmatism [].

12 However, segmented mirrors are rone to loose their continuity (hasing) if the load is highly uncorrelated, as is the case for side winds in our exeriment. From a control oint of view a high bandwidth highly uncorrelated wind load on a segmented rimary mirror can ose a substantial challenge. It is worth noting that the correlation length we measured at Gemini is about the size of the mirrors most of the lanned extremely large segmented mirror telescoes consider. For these telescoes, either the rimary mirror control system, or the enclosure, or rather both should be designed with this fact in mind. Since any structural element or baffle inside the rimary mirror area is likely to generate high bandwidth, local turbulence, the need for them should be carefully balanced. The increased frequency of wind buffeting even on relatively small areas may cause unwanted resonances in extremely large telescoes. In this aer we considered the cases where we could estimate the wind velocity distribution on the rimary mirror from the measurements taken at the edges of the mirror. To searate local and global (background) effects in the general case when the wind blows on the mirror in an arbitrary angle, further investigation is necessary. Comutational fluid dynamics (CFD) simulations can redict the velocity distribution above the mirror, which in turn yields the ressure distribution. CFD simulations could also estimate the local turbulence due to obstacles. ACKNOWLEDGEMENTS The authors would like to acknowledge the suort of the Gemini staff on Cerro Pachon in setting u and erforming the wind tests, articularly John Roberts, Gabriel Perez, Pedro Gigoux, Manual Lazo, and Pablo Prado. We would like to thank Dr. David Smith of MERLAB for organizing the dynamic measurements and collaborating in the data reduction. Secial thanks are also extended to Dr. Oleg Likhatchev for his helful comments and suggestions. The New Initiatives Office is a artnershi between two divisions of the Association of Universities for Research in Astronomy (AURA), Inc.: the National Otical Astronomy Observatory (NOAO) and the Gemini Observatory. NOAO is oerated by AURA under cooerative agreement with the National Science Foundation (NSF). The Gemini Observatory is oerated by AURA under a cooerative agreement with the NSF on behalf of the Gemini artnershi: the National Science Foundation (United States), the Particle Physics and Astronomy Research Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), CNPq (Brazil) and CONICET (Argentina). REFERENCES [] S. K. Friedlander and L. Toer, Turbulence - Classic Paers on Statistical Theory, (Interscience Publishers, Inc., New York, 96). [] A. N. Kolmogorov, "The Local Structure of Turbulence in Incomressible Viscous Fluid for Very Large Reynolds Numbers," Comtes rendus de l'academie des sciences de l'u. R. S. S., -5 (94). [] S. F. Clifford, Laser Beam Proagation in the Atmoshere, J. W. Strohbehn, ed., (Sringer-Verlag, Berlin, 978). [4] V. I. Tatarski, Wave Proagation in a Turbulent Medium, (McGraw-Hill, New York, 96). [5] G. K. Batchelor, The Theory of Homogeneous Turbulence, (Cambridge University Press, Cambridge, 959). [6] E. Simiu and R. H. Scanlan, Wind Effects on Structures - An Introduction to Wind Engineering, (John Wiley & Sons, Inc., New York, 986). [7] T. von Karman, "Progress in the Statistical Theory of Turbulence," Journal of Marine Research 7, 5-64 (948). [8] A. G. Davenort, "The Sectrum of Horizontal Gustiness Near the Ground in High Winds," Journal of the Royal Meteorological Society 87, 94- (96). [9] D. A. Antoniou, "Turbulence Measurements on To of a Stee Hill," Journal of Wind Engineering and Industrial Aerodynamics 9, 4-55 (99). [] M. Ravensbergen, "Main Axes Servo Systems of the VLT," Proceedings of SPIE 99, (994).

13 [] I. Linares, J. Carlos, and C. Hude, "Overview of the Gran Telescoio Canarias Global Model Simulation," Proceedings of SPIE 44, (). [] Cho, Myung, Ste, Larry and Kim, Seongho., "Wind Buffeting Effects on the Gemini 8m Primary Mirrors," Proceedings of SPIE 4444, -4 (). [] Smith, David, Weech, Keith, Teutsch, Johann, Avitabile, Peter, Gwaltney, Geoff, Sheehan, Mike, and Cho, Myung. Comarison of Pressure Measurements and Oerating Data for Wind Excitation of Telescoe Structures.. Proceedings of the 9th International Modal Analysis Conference. Kissimmee, FL [4] Avitabile, Peter, Weech, Keith, Smith, David, Gwaltney, Geoff, and Sheehan, Mike. Modal and Oerating Characterization of an Otical Telescoe.. Proceedings of the 9th International Modal Analysis Conference. Kissimmee, FL [5] Signal Processing Toolbox User's Guide, (The MathWorks, Inc., ).

The directivity of the forced radiation of sound from panels and openings including the shadow zone

The directivity of the forced radiation of sound from panels and openings including the shadow zone The directivity of the forced radiation of sound from anels and oenings including the shadow zone J. Davy RMIT University, Alied Physics, GPO Box 476V, 3001 Melbourne, Victoria, Australia john.davy@rmit.edu.au

More information

Integrating Lidar and Atmospheric Boundary Layer Measurements to Determine Fluxes and Dynamics of Particulate Emissions from an Agriculture Facility

Integrating Lidar and Atmospheric Boundary Layer Measurements to Determine Fluxes and Dynamics of Particulate Emissions from an Agriculture Facility Integrating Lidar and Atmosheric Boundary Layer Measurements to Determine Fluxes and Dynamics of Particulate Emissions from an Agriculture Facility Lawrence His 1, John H. Prueger 2, Jerry Hatfield 2,

More information

Hotelling s Two- Sample T 2

Hotelling s Two- Sample T 2 Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test

More information

Active optics challenges of a thirty meter segmented mirror telescope

Active optics challenges of a thirty meter segmented mirror telescope Active optics challenges of a thirty meter segmented mirror telescope George Z. Angeli 1, Robert Upton 1, Anna Segurson 1, Brent Ellerbroek 1 1 New Initiatives Office, AURA Inc. ABSTRACT Ground-based telescopes

More information

APPENDIX 5.5.C WIND BUFFETING EFFECTS ON THE GEMINI 8M PRIMARY MIRRORS

APPENDIX 5.5.C WIND BUFFETING EFFECTS ON THE GEMINI 8M PRIMARY MIRRORS APPENDIX 5.5.C WIND BUFFETING EFFECTS ON THE GEMINI 8M PRIMARY MIRRORS Paper originally appeared in SPIE Proceedings Volume 4444, "Optomechanical Design and Engineering 2001". Wind buffeting effects on

More information

Chapter 1 Fundamentals

Chapter 1 Fundamentals Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors

More information

arxiv:cond-mat/ v2 25 Sep 2002

arxiv:cond-mat/ v2 25 Sep 2002 Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,

More information

An Improved Calibration Method for a Chopped Pyrgeometer

An Improved Calibration Method for a Chopped Pyrgeometer 96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum

More information

Estimation of the large covariance matrix with two-step monotone missing data

Estimation of the large covariance matrix with two-step monotone missing data Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo

More information

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)

More information

Numerical Simulation of Particle Concentration in a Gas Cyclone Separator *

Numerical Simulation of Particle Concentration in a Gas Cyclone Separator * 2007 Petroleum Science Vol.4 No.3 Numerical Simulation of Particle Concentration in a Gas Cyclone Searator * Xue Xiaohu, Sun Guogang **, Wan Gujun and Shi Mingxian (School of Chemical Science and Engineering,

More information

Seafloor Reflectivity A Test of an Inversion Technique

Seafloor Reflectivity A Test of an Inversion Technique Seafloor Reflectivity A Test of an Inversion Technique Adrian D. Jones 1, Justin Hoffman and Paul A. Clarke 1 1 Defence Science and Technology Organisation, Australia, Student at Centre for Marine Science

More information

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO) Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment

More information

On Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada

On Using FASTEM2 for the Special Sensor Microwave Imager (SSM/I) March 15, Godelieve Deblonde Meteorological Service of Canada On Using FASTEM2 for the Secial Sensor Microwave Imager (SSM/I) March 15, 2001 Godelieve Deblonde Meteorological Service of Canada 1 1. Introduction Fastem2 is a fast model (multile-linear regression model)

More information

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites

More information

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment

More information

VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES

VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Journal of Sound and Vibration (998) 22(5), 78 85 VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Acoustics and Dynamics Laboratory, Deartment of Mechanical Engineering, The

More information

Use of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek

Use of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.

More information

Statics and dynamics: some elementary concepts

Statics and dynamics: some elementary concepts 1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and

More information

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system

More information

4. Score normalization technical details We now discuss the technical details of the score normalization method.

4. Score normalization technical details We now discuss the technical details of the score normalization method. SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules

More information

5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:

5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell: 5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum

More information

Homogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers

Homogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas

More information

Optimal Recognition Algorithm for Cameras of Lasers Evanescent

Optimal Recognition Algorithm for Cameras of Lasers Evanescent Otimal Recognition Algorithm for Cameras of Lasers Evanescent T. Gaudo * Abstract An algorithm based on the Bayesian aroach to detect and recognise off-axis ulse laser beams roagating in the atmoshere

More information

Maximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming

Maximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August

More information

State Estimation with ARMarkov Models

State Estimation with ARMarkov Models Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,

More information

CFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS

CFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS CFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS M. BORGHI - M. MILANI Diartimento di Scienze dell Ingegneria Università degli Studi di Modena Via Cami, 213/b 41100 Modena E-mail: borghi@omero.dsi.unimo.it

More information

Damage Identification from Power Spectrum Density Transmissibility

Damage Identification from Power Spectrum Density Transmissibility 6th Euroean Worksho on Structural Health Monitoring - h.3.d.3 More info about this article: htt://www.ndt.net/?id=14083 Damage Identification from Power Sectrum Density ransmissibility Y. ZHOU, R. PERERA

More information

Churilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics

Churilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated

More information

Principles of Computed Tomography (CT)

Principles of Computed Tomography (CT) Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel

More information

Wind loading of large telescopes

Wind loading of large telescopes Wind loading of large telescopes Myung K. Cho a, Larry M. Stepp a, George Z. Angeli a, David R. Smith b a AURA New Initiatives Office; b Merlab, P.C. ABSTRACT Wind loading is a critical issue for large

More information

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,

More information

PHYSICAL REVIEW LETTERS

PHYSICAL REVIEW LETTERS PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.

More information

Chapter 2 Introductory Concepts of Wave Propagation Analysis in Structures

Chapter 2 Introductory Concepts of Wave Propagation Analysis in Structures Chater 2 Introductory Concets of Wave Proagation Analysis in Structures Wave roagation is a transient dynamic henomenon resulting from short duration loading. Such transient loadings have high frequency

More information

MATH 2710: NOTES FOR ANALYSIS

MATH 2710: NOTES FOR ANALYSIS MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite

More information

Participation Factors. However, it does not give the influence of each state on the mode.

Participation Factors. However, it does not give the influence of each state on the mode. Particiation Factors he mode shae, as indicated by the right eigenvector, gives the relative hase of each state in a articular mode. However, it does not give the influence of each state on the mode. We

More information

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment

More information

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS 4 th International Conference on Earthquake Geotechnical Engineering June 2-28, 27 KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS Misko CUBRINOVSKI 1, Hayden BOWEN 1 ABSTRACT Two methods for analysis

More information

Determination of Pressure Losses in Hydraulic Pipeline Systems by Considering Temperature and Pressure

Determination of Pressure Losses in Hydraulic Pipeline Systems by Considering Temperature and Pressure Paer received: 7.10.008 UDC 61.64 Paer acceted: 0.04.009 Determination of Pressure Losses in Hydraulic Pieline Systems by Considering Temerature and Pressure Vladimir Savi 1,* - Darko Kneževi - Darko Lovrec

More information

%(*)= E A i* eiujt > (!) 3=~N/2

%(*)= E A i* eiujt > (!) 3=~N/2 CHAPTER 58 Estimating Incident and Reflected Wave Fields Using an Arbitrary Number of Wave Gauges J.A. Zelt* A.M. ASCE and James E. Skjelbreia t A.M. ASCE 1 Abstract A method based on linear wave theory

More information

Deriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.

Deriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &

More information

Tests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test)

Tests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test) Chater 225 Tests for Two Proortions in a Stratified Design (Cochran/Mantel-Haenszel Test) Introduction In a stratified design, the subects are selected from two or more strata which are formed from imortant

More information

Time Domain Calculation of Vortex Induced Vibration of Long-Span Bridges by Using a Reduced-order Modeling Technique

Time Domain Calculation of Vortex Induced Vibration of Long-Span Bridges by Using a Reduced-order Modeling Technique 2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Time Domain Calculation of Vortex Induced Vibration of Long-San Bridges by Using a Reduced-order Modeling

More information

On the relationship between sound intensity and wave impedance

On the relationship between sound intensity and wave impedance Buenos Aires 5 to 9 Setember, 16 Acoustics for the 1 st Century PROCEEDINGS of the nd International Congress on Acoustics Sound Intensity and Inverse Methods in Acoustics: Paer ICA16-198 On the relationshi

More information

Optical Fibres - Dispersion Part 1

Optical Fibres - Dispersion Part 1 ECE 455 Lecture 05 1 Otical Fibres - Disersion Part 1 Stavros Iezekiel Deartment of Electrical and Comuter Engineering University of Cyrus HMY 445 Lecture 05 Fall Semester 016 ECE 455 Lecture 05 Otical

More information

On split sample and randomized confidence intervals for binomial proportions

On split sample and randomized confidence intervals for binomial proportions On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have

More information

Research of power plant parameter based on the Principal Component Analysis method

Research of power plant parameter based on the Principal Component Analysis method Research of ower lant arameter based on the Princial Comonent Analysis method Yang Yang *a, Di Zhang b a b School of Engineering, Bohai University, Liaoning Jinzhou, 3; Liaoning Datang international Jinzhou

More information

FE FORMULATIONS FOR PLASTICITY

FE FORMULATIONS FOR PLASTICITY G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND

More information

Effect of geometry on flow structure and pressure drop in pneumatic conveying of solids along horizontal ducts

Effect of geometry on flow structure and pressure drop in pneumatic conveying of solids along horizontal ducts Journal of Scientific LAÍN & Industrial SOMMERFELD Research: PNEUMATIC CONVEYING OF SOLIDS ALONG HORIZONTAL DUCTS Vol. 70, February 011,. 19-134 19 Effect of geometry on flow structure and ressure dro

More information

Simplifications to Conservation Equations

Simplifications to Conservation Equations Chater 5 Simlifications to Conservation Equations 5.1 Steady Flow If fluid roerties at a oint in a field do not change with time, then they are a function of sace only. They are reresented by: ϕ = ϕq 1,

More information

System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests

System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests 009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract

More information

8.7 Associated and Non-associated Flow Rules

8.7 Associated and Non-associated Flow Rules 8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more

More information

A Note on Massless Quantum Free Scalar Fields. with Negative Energy Density

A Note on Massless Quantum Free Scalar Fields. with Negative Energy Density Adv. Studies Theor. Phys., Vol. 7, 13, no. 1, 549 554 HIKARI Ltd, www.m-hikari.com A Note on Massless Quantum Free Scalar Fields with Negative Energy Density M. A. Grado-Caffaro and M. Grado-Caffaro Scientific

More information

A compression line for soils with evolving particle and pore size distributions due to particle crushing

A compression line for soils with evolving particle and pore size distributions due to particle crushing Russell, A. R. (2011) Géotechnique Letters 1, 5 9, htt://dx.doi.org/10.1680/geolett.10.00003 A comression line for soils with evolving article and ore size distributions due to article crushing A. R. RUSSELL*

More information

Uniformly best wavenumber approximations by spatial central difference operators: An initial investigation

Uniformly best wavenumber approximations by spatial central difference operators: An initial investigation Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations

More information

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College

More information

Effective conductivity in a lattice model for binary disordered media with complex distributions of grain sizes

Effective conductivity in a lattice model for binary disordered media with complex distributions of grain sizes hys. stat. sol. b 36, 65-633 003 Effective conductivity in a lattice model for binary disordered media with comlex distributions of grain sizes R. PIASECKI Institute of Chemistry, University of Oole, Oleska

More information

Linear diophantine equations for discrete tomography

Linear diophantine equations for discrete tomography Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,

More information

Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting

Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting CLAS-NOTE 4-17 Determining Momentum and Energy Corrections for g1c Using Kinematic Fitting Mike Williams, Doug Alegate and Curtis A. Meyer Carnegie Mellon University June 7, 24 Abstract We have used the

More information

Uniform Law on the Unit Sphere of a Banach Space

Uniform Law on the Unit Sphere of a Banach Space Uniform Law on the Unit Shere of a Banach Sace by Bernard Beauzamy Société de Calcul Mathématique SA Faubourg Saint Honoré 75008 Paris France Setember 008 Abstract We investigate the construction of a

More information

RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES

RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete

More information

Shadow Computing: An Energy-Aware Fault Tolerant Computing Model

Shadow Computing: An Energy-Aware Fault Tolerant Computing Model Shadow Comuting: An Energy-Aware Fault Tolerant Comuting Model Bryan Mills, Taieb Znati, Rami Melhem Deartment of Comuter Science University of Pittsburgh (bmills, znati, melhem)@cs.itt.edu Index Terms

More information

On the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined

On the Fluid Dependence of Rock Compressibility: Biot-Gassmann Refined Downloaded 0/9/3 to 99.86.4.8. Redistribution subject to SEG license or coyright; see Terms of Use at htt://library.seg.org/ On the luid Deendence of Rock Comressibility: Biot-Gassmann Refined Leon Thomsen,

More information

Elementary Analysis in Q p

Elementary Analysis in Q p Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some

More information

Sound speed profile structure and variability measured over flat terrain

Sound speed profile structure and variability measured over flat terrain Sound seed rofile structure and variability measured over flat terrain Bradley, SG, Waddington, DC and Von Hunerbein, S Title Authors Tye URL Published Date 6 Sound seed rofile structure and variability

More information

Montgomery self-imaging effect using computer-generated diffractive optical elements

Montgomery self-imaging effect using computer-generated diffractive optical elements Otics Communications 225 (2003) 13 17 www.elsevier.com/locate/otcom Montgomery self-imaging effect using comuter-generated diffractive otical elements J urgen Jahns a, *, Hans Knuertz a, Adolf W. Lohmann

More information

arxiv: v1 [physics.data-an] 26 Oct 2012

arxiv: v1 [physics.data-an] 26 Oct 2012 Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch

More information

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional

More information

Brownian Motion and Random Prime Factorization

Brownian Motion and Random Prime Factorization Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........

More information

Approximating min-max k-clustering

Approximating min-max k-clustering Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost

More information

VIBRATIONS OF SHALLOW SPHERICAL SHELLS AND GONGS: A COMPARATIVE STUDY

VIBRATIONS OF SHALLOW SPHERICAL SHELLS AND GONGS: A COMPARATIVE STUDY VIBRATIONS OF SHALLOW SPHERICAL SHELLS AND GONGS: A COMPARATIVE STUDY PACS REFERENCE: 43.75.Kk Antoine CHAIGNE ; Mathieu FONTAINE ; Olivier THOMAS ; Michel FERRE ; Cyril TOUZE UER de Mécanique, ENSTA Chemin

More information

3.4 Design Methods for Fractional Delay Allpass Filters

3.4 Design Methods for Fractional Delay Allpass Filters Chater 3. Fractional Delay Filters 15 3.4 Design Methods for Fractional Delay Allass Filters Above we have studied the design of FIR filters for fractional delay aroximation. ow we show how recursive or

More information

Probability Estimates for Multi-class Classification by Pairwise Coupling

Probability Estimates for Multi-class Classification by Pairwise Coupling Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics

More information

Statistics II Logistic Regression. So far... Two-way repeated measures ANOVA: an example. RM-ANOVA example: the data after log transform

Statistics II Logistic Regression. So far... Two-way repeated measures ANOVA: an example. RM-ANOVA example: the data after log transform Statistics II Logistic Regression Çağrı Çöltekin Exam date & time: June 21, 10:00 13:00 (The same day/time lanned at the beginning of the semester) University of Groningen, Det of Information Science May

More information

ONE. The Earth-atmosphere system CHAPTER

ONE. The Earth-atmosphere system CHAPTER CHAPTER ONE The Earth-atmoshere system 1.1 INTRODUCTION The Earth s atmoshere is the gaseous enveloe surrounding the lanet. Like other lanetary atmosheres, it figures centrally in transfers of energy between

More information

ε(ω,k) =1 ω = ω'+kv (5) ω'= e2 n 2 < 0, where f is the particle distribution function and v p f v p = 0 then f v = 0. For a real f (v) v ω (kv T

ε(ω,k) =1 ω = ω'+kv (5) ω'= e2 n 2 < 0, where f is the particle distribution function and v p f v p = 0 then f v = 0. For a real f (v) v ω (kv T High High Power Power Laser Laser Programme Programme Theory Theory and Comutation and Asects of electron acoustic wave hysics in laser backscatter N J Sircombe, T D Arber Deartment of Physics, University

More information

ON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS

ON THE LEAST SIGNIFICANT p ADIC DIGITS OF CERTAIN LUCAS NUMBERS #A13 INTEGERS 14 (014) ON THE LEAST SIGNIFICANT ADIC DIGITS OF CERTAIN LUCAS NUMBERS Tamás Lengyel Deartment of Mathematics, Occidental College, Los Angeles, California lengyel@oxy.edu Received: 6/13/13,

More information

Plotting the Wilson distribution

Plotting the Wilson distribution , Survey of English Usage, University College London Setember 018 1 1. Introduction We have discussed the Wilson score interval at length elsewhere (Wallis 013a, b). Given an observed Binomial roortion

More information

Positive decomposition of transfer functions with multiple poles

Positive decomposition of transfer functions with multiple poles Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.

More information

Feedback-error control

Feedback-error control Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller

More information

Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models

Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu

More information

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi

More information

Ultrasound Beam Focusing Considering the Cutaneous Fat Layer Effects

Ultrasound Beam Focusing Considering the Cutaneous Fat Layer Effects Ultrasound Beam Focusing Considering the Cutaneous Fat Layer Effects A. B. M. Aowlad Hossain 1*, Laehoon H. Kang 1 Deartment of Electronics and Communication Engineering Khulna University of Engineering

More information

Notes on pressure coordinates Robert Lindsay Korty October 1, 2002

Notes on pressure coordinates Robert Lindsay Korty October 1, 2002 Notes on ressure coordinates Robert Lindsay Korty October 1, 2002 Obviously, it makes no difference whether the quasi-geostrohic equations are hrased in height coordinates (where x, y,, t are the indeendent

More information

Meshless Methods for Scientific Computing Final Project

Meshless Methods for Scientific Computing Final Project Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking

More information

Chapter 6. Thermodynamics and the Equations of Motion

Chapter 6. Thermodynamics and the Equations of Motion Chater 6 hermodynamics and the Equations of Motion 6.1 he first law of thermodynamics for a fluid and the equation of state. We noted in chater 4 that the full formulation of the equations of motion required

More information

Numerical and experimental investigation on shot-peening induced deformation. Application to sheet metal forming.

Numerical and experimental investigation on shot-peening induced deformation. Application to sheet metal forming. Coyright JCPDS-International Centre for Diffraction Data 29 ISSN 197-2 511 Numerical and exerimental investigation on shot-eening induced deformation. Alication to sheet metal forming. Florent Cochennec

More information

Submicrometer Position Control of Single Trapped Neutral Atoms

Submicrometer Position Control of Single Trapped Neutral Atoms Dotsenko, I and Alt, W and Khudaverdyan, M and Kuhr, S and Meschede, D and Miroshnychenko, Y and Schrader, D and Rauschenbeutel, A (25) Submicrometer osition control of single traed neutral atoms. Physical

More information

Wolfgang POESSNECKER and Ulrich GROSS*

Wolfgang POESSNECKER and Ulrich GROSS* Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND

More information

Performance of a First-Level Muon Trigger with High Momentum Resolution Based on the ATLAS MDT Chambers for HL-LHC

Performance of a First-Level Muon Trigger with High Momentum Resolution Based on the ATLAS MDT Chambers for HL-LHC Performance of a First-Level Muon rigger with High Momentum Resolution Based on the ALAS MD Chambers for HL-LHC P. Gadow, O. Kortner, S. Kortner, H. Kroha, F. Müller, R. Richter Max-Planck-Institut für

More information

All-fiber Optical Parametric Oscillator

All-fiber Optical Parametric Oscillator All-fiber Otical Parametric Oscillator Chengao Wang Otical Science and Engineering, Deartment of Physics & Astronomy, University of New Mexico Albuquerque, NM 87131-0001, USA Abstract All-fiber otical

More information

Chapter 7 Rational and Irrational Numbers

Chapter 7 Rational and Irrational Numbers Chater 7 Rational and Irrational Numbers In this chater we first review the real line model for numbers, as discussed in Chater 2 of seventh grade, by recalling how the integers and then the rational numbers

More information

Multiparameter entanglement in quantum interferometry

Multiparameter entanglement in quantum interferometry PHYSICAL REVIEW A, 66, 023822 200 Multiarameter entanglement in quantum interferometry Mete Atatüre, 1 Giovanni Di Giusee, 2 Matthew D. Shaw, 2 Alexander V. Sergienko, 1,2 Bahaa E. A. Saleh, 2 and Malvin

More information

8 STOCHASTIC PROCESSES

8 STOCHASTIC PROCESSES 8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular

More information

Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis

Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty

More information

Analysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia

Analysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.

More information

Preliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure. Balance by FEM Calculations

Preliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure. Balance by FEM Calculations Preliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure Balance by FEM Calculations G. Molinar*, M. Bergoglio*, G. Mosso*,G. Buonanno**, M. Dell Isola** * Istituto di Metrologia

More information

Characteristics of Beam-Based Flexure Modules

Characteristics of Beam-Based Flexure Modules Shorya Awtar e-mail: shorya@mit.edu Alexander H. Slocum e-mail: slocum@mit.edu Precision Engineering Research Grou, Massachusetts Institute of Technology, Cambridge, MA 039 Edi Sevincer Omega Advanced

More information

John Weatherwax. Analysis of Parallel Depth First Search Algorithms

John Weatherwax. Analysis of Parallel Depth First Search Algorithms Sulementary Discussions and Solutions to Selected Problems in: Introduction to Parallel Comuting by Viin Kumar, Ananth Grama, Anshul Guta, & George Karyis John Weatherwax Chater 8 Analysis of Parallel

More information

Shear Layers and Aperture Effects for Aero-Optics

Shear Layers and Aperture Effects for Aero-Optics 36th AIAA Plasmadynamics and Lasers onference 6-9 June 5, Toronto, anada AIAA-5-77 Shear Layers and Aerture Effects for Aero-Otics John P Siegenthaler *, Stanislav Gordeyev, and Eric Jumer Hessert Laboratory,

More information