6B.7 EXPERIMENTAL VALIDATION OF THE WEBB CORRECTION FOR CO2 FLUX WITH AN OPEN-PATH GAS ANALYZER
|
|
- Joan Wells
- 5 years ago
- Views:
Transcription
1 6B.7 EXPERIMENAL VALIDAION OF HE WEBB CORRECION FOR CO2 FLUX WIH AN OPEN-PAH GAS ANALYZER Osamu sukamoto* an Fumiyoshi Kono Okayama University, Okayama, Japan 1. INRODUCION urbulent flux by the ey covariance technique is the only irect measurement over earth surfaces. CO2 an water vapor are the only trace gases that can measure turbulent flux with the open-path gas analyzer. he open-path CO2/H2O gas analyzer oes not measure non-imensional CO2 an water vapor concentrations as mixing ratios but it measures CO2 an water vapor ensities. For this reason, the trace gas flux using this analyzer nees to be correcte for the mean vertical flow ue to air ensity fluctuation. Webb et al. (1980) suggeste that the flux ue to the mean vertical flow cannot be neglecte for trace gases such as water vapor an CO2. o evaluate the magnitue of the influence by the mean vertical flow, Webb et al. (1980) assume that the vertical flux of ry air shoul be zero for the mass conservation. Practically, sensible an latent heat fluxes evaluate by the ey covariance technique are use to calculate water vapor an CO2 fluxes by the mean vertical flow (Webb or WPL correction). he Webb correction is important over sea surface CO2 flux as it is much smaller than the vegetate surface. Ohtaki et al. (1989) suggeste that the Webb correction term of CO2 flux is larger than the raw CO2 flux an changes the transport irection of the CO2 flux over the coastal sea surface. An recently, Kono an sukamoto (2007) also confirme the importance of the Webb correction over the open ocean surface. he basic Webb correction theory has been iscusse well incluing the recent iscussions (e.g. Leuning,2007), but this correction is still poorly unerstoo from actual observe ata in various conitions. herefore, it is important to evaluate the original Webb correction scheme using the experimental ata.. For this purpose, we measure the turbulent flux over an asphalt surface (large parking lot), where CO2 an water vapor fluxes were consiere almost zero an sensible heat flux was significant *Corresponing author ; Osamu sukamoto, Okayama University, Okayama, Japan tsuka@cc.okayama-u.ac.jp 2. EXPERIMENAL MEHOD he observation was mae in the large parking lot on August 24, his site is locate at a coast in Okayama, Japan. his parking lot extens 110 m x 180 m an provie a uniform an ry asphalt surface. he fine weather continue an the surface was very ry. Sea breeze ominate uring this observation perio. Fluctuations of three-imensional win velocities an soun virtual temperature were measure using a sonic anemometer thermometer (KAIJO, DA-600-3V). he SA probe was irecte against the ominant sea breeze. CO2 an water vapor ensity fluctuations were measure using an infrare open-path CO2/H2O gas analyzer (LI- COR, LI-7500). Both instruments were mounte on a tripo at a 1.62 m height above the asphalt surface an separate horizontally by 0.25 m. he observation point was locate near the northwest ege of the parking lot to ensure the sufficient fetch from the sea breeze. he fetch to the southeast coast was more than 180 m. he analog output signals from these instruments were sample an igitize at the rate of 10 Hz by a ata acquisition system (LabVIEW, National Instruments) on a PC. 3. EDDY FLUX EVALUAION INCLUDING HE WEBB CORRECION H = he sensible heat flux (H) is given by C ρ w ' ' p p = C ρ ( w ' ' sv sv w' q') (1) = sv ( q) (2) where, C p : specific heat of air at constant pressure (J kg -1 K -1 ) ρ : ensity of ry air (kg m -3 ) w : vertical win velocity (m s -1 ) : air temperature (K) sv : soun virtual temperature (K) Q : specific humiity (kg kg -1 ).
2 Win velocities measure by the sonic anemometer were rotate into a coorinate system aligne with the mean win (McMillen, 1988). he soun virtual temperature by the sonic thermometer inclue the crosswin correction (Hignett, 1992). Humiity correction was applie to the soun virtual temperature to evaluate the real sensible heat flux (Eq. 1) an air temperature (Eq. 2). In this paper, the specific humiity fluctuation in Eq. 1 was calculate as the ratio of water vapor ensity fluctuation measure by the infrare open-path CO 2 /H 2 O gas analyzer to air ensity. As shown in Fig. 6, the correcte sensible heat flux was almost ientical to the raw sensible heat flux, as the humiity correction was negligibly small in the present experiment. hen, the total latent heat (λe) an CO 2 (F c ) fluxes incluing the Webb correction terms were calculate as follows. ρv λ E λ ( 1+ μσ )( w' ρ' v + w' ') ρc ρc F = w ρ c + μ w ρ v + 1+ μσ ) w ρ = (3) ( (4) c where, λ : latent heat of vaporization (J kg -1 ) μ : ratio of the molecular weights of ry air an water vapor σ : ratio of water vapor an ry air ensities ρ v an ρ c : H 2 O an CO 2 ensities by the open-path CO 2 /H 2 O gas analyzer (kg m -3 ). he first an secon terms of the right parenthesis on the right-han sie in Eq. 3 are the raw water vapor flux an the water vapor Webb correction term ue to sensible heat flux. he first, secon, an thir terms on the right-han sie in Eq. 4 are the raw CO 2 flux an the CO 2 Webb correction terms ue to latent an sensible heat fluxes, respectively. Detaile erivations of Eqs. 3 an 4 can be foun in Webb et al. (1980). 4. RESULS 4.1 urbulent fluctuations Figure 1 shows a typical example of turbulent fluctuations for vertical win velocity, soun virtual temperature, an water vapor an CO 2 ensities measure over the asphalt surface uring the aytime. hese fluctuations are plotte as 10 Hz raw ata as the eviation from the 5 minutes mean value. he soun virtual temperature showe positively skewe fluctuation with an amplitue of 5 K, an was positively correlate with vertical win velocity. he CO 2 ensity showe negatively skewe fluctuation with an amplitue of 10 mg m -3 (= 6 ppm), an was negatively correlate with vertical win velocity. Fluctuations of soun virtual temperature an CO 2 ensity were note to be istinctly an negatively correlate with each other using an inepenent instrument. hese results showe that the ey transport of CO 2 was apparently ownwar from the atmosphere to the asphalt surface with the large upwar sensible heat flux. he fluctuation of water vapor ensity was small an i not show correlation with the fluctuation of soun virtual temperature or CO 2 ensity. hese results showe that there is no cross-sensitivity between H 2 O an CO 2 gases by an infrare open-path CO 2 /H 2 O gas analyzer. In the present stuy, the total CO 2 flux was etermine by the raw CO 2 flux an the CO 2 Webb correction term ue to sensible heat flux, because the Webb correction term ue to latent heat flux can be negligible in the present ata. 1.2 Vertical Win Velocity (ms-1) Soun Virtual emperature (K) Water Vapor Density (g m-3) Carbon Dioxie Density (mg m-3) IME (0.1s) 10:45-10:50 Aug. 24, 2006 Figure 1 urbulent fluctuations of vertical win velocity, soun virtual temperature, an water vapor an CO2 ensities 4.2 Power spectra an co-spectra Figure 2 shows the normalize power spectra of the vertical win velocity, soun virtual temperature, an water vapor an CO 2 ensity fluctuations observe by the open-path instruments between 10:45 an 11:00. he normalize power spectral ensities S xx (n)/σ 2, were plotte as functions of nonimensional frequency f = nz/u. hese power spectra ha similar frequency structure, an followe the 5/3 power law in the inertial subrange above As shown in Figure 3, the
3 normalize power spectra of CO 2 ensity fluctuation show similar frequency structures uring this observation perio. the CO 2 Webb correction term ue to sensible heat flux showe positive for the whole frequency range, an was negatively correlate with that of raw CO 2 flux. In contrast, the cospectrum of the CO 2 Webb correction term ue to latent heat flux was negligible for the whole frequency range. he cospectrum of the resulte total CO 2 flux showe a positive value for the whole frequency range. Figure 2 Normalize power spectra of the vertical win velocity, soun virtual temperature, an water vapor an CO 2 ensity Figure 4 Cospectra of CO 2 ensity an vertical win velocity normalize by the raw CO 2 flux as functions of non-imensional frequency Figure 3 Normalize power spectra of the CO 2 ensity fluctuations as functions of non-imensional frequency Figure 4 shows the cospectra of CO 2 ensity an vertical win velocity normalize by raw CO 2 flux as functions of the nonimensional frequency. hese were analyze using the same runs as shown in Fig. 3. hese cospectra also ha similar frequency structure between each run, an the flux ominant frequency ranges were aroun Figure 5 compares the ensemble-average cospectrum of raw CO 2 flux with those of the Webb CO 2 correction terms ue to latent an sensible heat fluxes, an total CO 2 flux. he cospectrum of Figure 5 Ensemble-average cospectra of the raw CO 2 flux, CO 2 Webb correction terms ue to latent an sensible heat fluxes, an total CO 2 flux 4.3 Webb correction on CO 2 flux over the asphalt surface Figure 6 represents the time series of raw an correcte sensible heat fluxes(eq.1), an raw an total latent heat fluxes(eq.3). he raw an correcte sensible heat fluxes are almost ientical
4 as the surface was very ry. he latent heat flux represent some positive flux when inclue the Webb correction as to be iscusse later. Figure 6 ime-series ata of raw an correcte sensible heat fluxes, an raw an total latent heat fluxes. he raw an correcte sensible heat fluxes are almost ientical. between 87 an 243 Wm -2 ). he Webb correction term ue to latent heat flux was between an mg m -2 s -1 (raw latent heat flux was between 21 an 41 Wm -2 ). otal(correcte) CO 2 flux was always positive between an mg m -2 s -1 (the average was mg m -2 s -1 ), an i not show a clear iurnal variation, as i raw CO 2 an sensible heat fluxes. Ham an Heilman (2003) performe a similar experiment over an asphalt surface, an reporte total CO 2 flux by the ey covariance technique of 0.02 mg m -2 s -1 (= 1.8 g m -2 ay -1 ) an CO 2 flux by the chamber technique of 0.03 mg m -2 s -1 (= 2.8 g m -2 ay -1 ) over the asphalt surface. he average total CO 2 flux in the present stuy was about two or three times larger than their result. Although total CO 2 flux measure by the ey covariance technique as reporte by Ham an Heilman (2003) inclue negative values, total CO 2 flux in the present stuy retaine a significantly positive value. 5. DISCUSSIONS As shown in Fig. 6, total latent heat flux incluing the Webb correction term was between 6 an 63 Wm -2 (the average was 32 Wm -2 ), an larger than the result reporte by Ham an Heilman (2003), showing systematic iurnal variation. Over the ry asphalt surface latent heat flux is consiere almost zero, total latent heat flux resulte in a significantly positive value. Figure 7 ime-series ata of the raw CO 2 flux, Webb CO 2 correction terms for the sensible an latent heat fluxes, an total CO 2 fluxes (uncorrecte an correcte the horizontal separation). Figure 7 shows the time-series ata of the raw CO 2 flux, CO 2 Webb correction terms ue to sensible an latent heat fluxes, an total CO 2 flux. he raw CO 2 flux showe a iurnal variation, an was always negative (ownwar from the atmosphere to the asphalt surface) between an mg m -2 s -1 (the average was mg m -2 s -1 ). he CO 2 Webb correction term ue to sensible heat flux also showe a iurnal variation, an was always positive between an mg m -2 s -1 (the sensible heat flux was Figure 8 Ratio of cospectrum of the raw CO 2 flux to that of the CO 2 Webb correction term ue to sensible heat flux as a function of frequency. he error bars are the stanar eviation. In Fig. 5, the ensemble-average cospectrum of total CO 2 flux ha ominant frequency ranges aroun 0.1 Hz, an the range was consistent with
5 those of raw CO 2 flux an the CO 2 Webb correction term ue to sensible heat flux. Above 0.62 Hz, the cospectrum of raw CO 2 flux attenuate faster than that of the CO 2 Webb correction term ue to sensible heat flux. As a result, the cospectrum of total CO 2 flux kept a positive value in this frequency range. his can be cause by horizontal separation between the sonic anemometer an the open-path CO 2 /H 2 O gas analyzer. Kristensen et al. (1997) showe that the flux causes a loss of 13% at a height of 1 m with a sensor horizontal separation of 0.2 m. o improve the unerestimation of raw CO 2 flux by this separation, we have assume cospectral similarty. Figure 8 shows the ratio of the cospectrum of raw CO 2 flux to CO 2 Webb correction term ue to sensible heat flux as a function of frequency. he ratio is almost -1 between 0.04 an 0.62 Hz, suggesting negative cospectral similarity. While above 0.62 Hz the ratio increases from -1, suggesting the issimilarity. he issimilarity comes from the sensor separation mentione above. Assuming that the cospectral similarity is kept in the high frequency above 0.62Hz, the CO2 raw flux was ajuste to keep the ratio as -1. As the result, the raw (negative) CO2 flux increase by mg m -2 s -1 an correcte total (positive) CO 2 flux ecrease by mg m -2 s -1 on average as shown in Fig.7. his correction value was about 33% of the separation uncorrecte total CO 2 flux (0.064 mg m -2 s -1 ) or less than 8% of the raw CO 2 flux ( mg m -2 s -1 ). However, the average total CO 2 flux kept a positive value (0.043 mg m -2 s -1 ). CO 2 flux by the ey covariance technique using an infrare open-path CO 2 /H 2 O gas analyzer may inclue some uncertainties of the Webb correction, cross-sensitivity, energy balance closure, an so on. Our result suggeste that the Webb correction may cause an overestimate correction of CO 2 an H 2 O fluxes over a ry asphalt surface where there is no cross-sensitivity between H 2 O an CO 2 gases by an infrare open-path CO 2 /H 2 O gas analyzer. Alternative metho of the flux valiation can be chamber technique an/or close-path eycovariance. hese further stuies are require to evaluate the accuracy of the Webb correction by the ey covariance technique using an open-path CO 2 /H 2 O gas analyzer. REFERENCES Ham, J.M., an Heilman, J.L., 2003: Experimental test of ensity an energy-balance corrections on carbon ioxie flux as measure using open-path ey covariance, Agron. J., 95, Hignett, P., 1992: Corrections to temperature measurements with a sonic anemometer, Bounary-Layer Meteorol., 61, Kono, F., an sukamoto, O., 2007: Air-sea CO 2 flux by ey covariance technique in the equatorial Inian Ocean, J. Oceanogr., 63, Kono,F. an O.sukamoto, 2007: Evaluation of Webb correction on CO2 flux by ey covariance technique using open-path gas analyzer over asphalt surface, J.Agric. Meteorology, 64, 1-8. Kristensen, L., Mann, J., Oncley, S.P., an Wyngaar, J.C., 1997: How close is close enough when measuring scalar fluxes with isplace sensors?, J. Atmos. Oceanic echnol., 14, Leuning, R., 2007: he correct form of the Webb, Pearman an Leuning equation for ey fluxes of trace gases in steay an non-steay state, horizontally homogeneous flows, Bounary-Layer Meteorol., 123, McMillen, R.., 1988: An ey correlation technique with extene applicability non-simple terrain, Bounary-Layer Meteorol., 43, Ohtaki, E., sukamoto, O., Iwatani, Y., an Mitsuta, Y., 1989: Measurements of the carbon ioxie flux over the ocean, J. Meteorol. Soc. Jap., 67, Webb, E.K., Pearman, G.I., an Leuning, R., 1980: Correction of flux measurement for ensity effects ue to heat an water vapour transfer, Q. J. R. Meteoro. Soc., 106,
Air-Sea CO 2 Flux by Eddy Covariance Technique in the Equatorial Indian Ocean
Journal of Oceanography, Vol. 63, pp. 449 to 456, 2007 Air-Sea CO 2 Flux by Eddy Covariance Technique in the Equatorial Indian Ocean FUMIYOSHI KONDO* and OSAMU TSUKAMOTO Graduate School of Natural Science
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationProject 3 Convection
12.307 Project 3 Convection Loovica Illari an John Marshall March 2010 1 Backgroun The Earth is bathe in raiation from the Sun whose intensity peaks in the visible. In orer to maintain energy balance the
More informationSources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences
Sources an Sinks of Available Potential Energy in a Moist Atmosphere Olivier Pauluis 1 Courant Institute of Mathematical Sciences New York University Submitte to the Journal of the Atmospheric Sciences
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationATMO. Theoretical considerations on the energy balance closure. Frederik De Roo and Matthias Mauder
ATMO Theoretical considerations on the energy balance closure Frederik De Roo and Matthias Mauder Karlsruhe Institute of Technology, Campus Alpin Atmospheric Environmental Research (KIT/IMK-IFU) Garmisch-Partenkirchen
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationProblem 3.84 of Bergman. Consider one-dimensional conduction in a plane composite wall. The outer surfaces are exposed to a fluid at T
1/10 bergman3-84.xmc Problem 3.84 of Bergman. Consier one-imensional conuction in a plane composite wall. The outer surfaces are expose to a flui at T 5 C an a convection heat transfer coefficient of h1000
More informationMAE 210A FINAL EXAM SOLUTIONS
1 MAE 21A FINAL EXAM OLUTION PROBLEM 1: Dimensional analysis of the foling of paper (2 points) (a) We wish to simplify the relation between the fol length l f an the other variables: The imensional matrix
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More informationThe Standard Atmosphere. Dr Andrew French
The Stanar Atmosphere Dr Anrew French 1 The International Stanar Atmosphere (ISA) is an iealize moel of the variation of average air pressure an temperature with altitue. Assumptions: The atmosphere consists
More informationEvaporating droplets tracking by holographic high speed video in turbulent flow
Evaporating roplets tracking by holographic high spee vieo in turbulent flow Loïc Méès 1*, Thibaut Tronchin 1, Nathalie Grosjean 1, Jean-Louis Marié 1 an Corinne Fournier 1: Laboratoire e Mécanique es
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationDiscussion Reply to comment by Rannik on A simple method for estimating frequency response corrections for eddy covariance systems. W.J.
Agricultural and Forest Meteorology 07 (200) 247 25 Discussion Reply to comment by Rannik on A simple method for estimating frequency response corrections for eddy covariance systems W.J. Massman USDA/Forest
More informationFabrizio Pelliccia *, Stefania Bonafoni, Patrizia Basili University of Perugia, Perugia, Italy
2.6 COSMIC GPS RADIO OCCULAION: NEURAL NEWORKS FOR ROPOSPHERIC PROFILING OVER HE INERROPICAL OCEAN AREA Fabrizio Pelliccia *, Stefania Bonafoni, Patrizia Basili University of Perugia, Perugia, Italy Nazzareno
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationPARALLEL-PLATE CAPACITATOR
Physics Department Electric an Magnetism Laboratory PARALLEL-PLATE CAPACITATOR 1. Goal. The goal of this practice is the stuy of the electric fiel an electric potential insie a parallelplate capacitor.
More informationDamping Temperature Fluctuations in the EC155 Intake Tube
WHITE PAPER Damping Temperature Fluctuations in the EC155 Intake Tube Steve Sargent, Campbell Scientific, Inc. Introduction The eddy-covariance (EC) technique is widely used to quantify the exchange of
More informationCharacterization of lead zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy
JOURNAL OF APPLIED PHYSICS VOLUME 85, NUMBER 1 15 JUNE 1999 Characterization of lea zirconate titanate piezoceramic using high frequency ultrasonic spectroscopy Haifeng Wang, Wenhua Jiang, a) an Wenwu
More informationOctober 1991 J. Wang and Y. Mitsuta 587 NOTES AND CORRESPONDENCE. Turbulence Structure and Transfer Characteristics
October 1991 J. Wang and Y. Mitsuta 587 NOTES AND CORRESPONDENCE Turbulence Structure and Transfer Characteristics in the Surface Layer of the HEIFE Gobi Area By Jiemin Wang Lanzhou Institute of Plateau
More informationThis section outlines the methodology used to calculate the wave load and wave wind load values.
COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationInternational Conference KNOWLEDGE-BASED ORGANIZATION Vol. XXIII No
International Conference KNOWLEDGE-BAED ORGANIZATION Vol. XXIII No 3 2017 METHOD FOR DETERMINATION OF THE PARAMETER OF A MACHINE GUN UPENION MOUNTED ON AN ARMOURED VEHICLE vilen PIRDONOV, vilen TEFANOV
More informationADIT DEBRIS PROJECTION DUE TO AN EXPLOSION IN AN UNDERGROUND AMMUNITION STORAGE MAGAZINE
ADIT DEBRIS PROJECTION DUE TO AN EXPLOSION IN AN UNDERGROUND AMMUNITION STORAGE MAGAZINE Froe Opsvik, Knut Bråtveit Holm an Svein Rollvik Forsvarets forskningsinstitutt, FFI Norwegian Defence Research
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationATSC 5010 Physical Meteorology I Lab Lab 2 numerical integration; height/pressure relationships. à Begin with IDEAL GAS Law:
ASC 5010 Physical Meteology I Lab Lab 2 numerical integration; height/pressure relationships à Begin with IDEAL GAS Law: PV = mr Pα = R note : 1 = thus : α P = R his is true f a mixture of gases, f ry
More informationLQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR
5 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR Tatsunori Kaneko* an Yasuto Asano* * Department of Mechanical Engineering,
More informationTo understand how scrubbers work, we must first define some terms.
SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More information1. Filling an initially porous tube under a constant head imposed at x =0
Notes on Moving Bounary problems, Voller U o M, volle00@umn.eu. Filling an initially porous tube uner a constant hea impose at x =0 Governing equation is base on calculating the water volume lux by the
More information5-4 Electrostatic Boundary Value Problems
11/8/4 Section 54 Electrostatic Bounary Value Problems blank 1/ 5-4 Electrostatic Bounary Value Problems Reaing Assignment: pp. 149-157 Q: A: We must solve ifferential equations, an apply bounary conitions
More informationPerturbation Analysis and Optimization of Stochastic Flow Networks
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. XX, NO. Y, MMM 2004 1 Perturbation Analysis an Optimization of Stochastic Flow Networks Gang Sun, Christos G. Cassanras, Yorai Wari, Christos G. Panayiotou,
More information2.25 m. (a) Using Newton s laws of motion, explain why the student can gain an initial speed to leave the ground vertically.
NAME : F.5 ( ) MARS: /70 FORM FIVE PHYSICS TEST on MECHANICS Time Allowe: 70 minutes This test consists of two sections: Section A (structure type questions, 50 marks); Section B (multiple choice, 20 marks)
More informationA note on the Mooney-Rivlin material model
A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel
More informationChapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies
hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation
More informationA 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS
A 3D SEDIMENT TRANSPORT MODEL FOR COMBINED WAVE-CURRENT FLOWS Peifeng Ma 1 an Ole Secher Masen Accurate preiction of current velocity an bottom shear stress, which both can be significantly influence by
More informationNarrowband Oscillations in the Upper Equatorial Ocean. Part I: Interpretation as Shear Instabilities
VOLUME 41 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y MARCH 2011 Narrowban Oscillations in the Upper Equatorial Ocean. Part I: Interpretation as Shear Instabilities J. N. MOUM, J.D.NASH,
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationExperiment 2, Physics 2BL
Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor
More informationOutline. Calculus for the Life Sciences II. Introduction. Tides Introduction. Lecture Notes Differentiation of Trigonometric Functions
Calculus for the Life Sciences II c Functions Joseph M. Mahaffy, mahaffy@math.ssu.eu Department of Mathematics an Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State
More informationParameter estimation: A new approach to weighting a priori information
Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a
More informationNumerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number
Journal of Applie Flui Mechanics, Vol. 9, No. 2, pp. 643-651, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Investigation of Non-Stationary Parameters on Effective
More informationPerformance evaluation of an integrated open-path eddy covariance system in a cold desert environment
Performance evaluation of an integrated open-path eddy covariance system in a cold desert environment Wei Wang, Jiaping Xu, Yunqiu Gao, Ivan Bogoev, Jian Cui, Lichen Deng, Cheng Hu, Cheng Liu, Shoudong
More informationInfluence of Radiation on Product Yields in a Film Boiling Reactor
R&D NOTES Influence of Raiation on Prouct Yiels in a Film Boiling Reactor C. Thomas Aveisian, Wing Tsang, Terence Daviovits, an Jonah R. Allaben Sibley School of Mechanical an Aerospace Engineering, Cornell
More informationSimulation of Angle Beam Ultrasonic Testing with a Personal Computer
Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation
More informationConvective Fluxes: Sensible and Latent Heat Convective Fluxes Convective fluxes require Vertical gradient of temperature / water AND Turbulence ( mixing ) Vertical gradient, but no turbulence: only very
More informationShape Effect on Blind Frequency for Depth Inversion in Pulsed Thermography
Shape Effect on Blin Frequency for Depth Inversion in Pulse Thermography M. Genest 1, E. Grinzato 2, P. Bison 2, S. Marinetti 2 C. Ibarra-Castaneo 1, X. Malague 1 1 Electrical an Computing Eng. Dept.,
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationThe continuity equation
Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal
More informationProblem Set 2: Solutions
UNIVERSITY OF ALABAMA Department of Physics an Astronomy PH 102 / LeClair Summer II 2010 Problem Set 2: Solutions 1. The en of a charge rubber ro will attract small pellets of Styrofoam that, having mae
More informationChapter 9 Method of Weighted Residuals
Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationAnalytical approach of Brillouin amplification over threshold
Research Article Vol. 57, No. 4 / February 208 / Applie Optics 607 Analytical approach of Brillouin amplification over threshol FİKRİ SERDAR GÖKHAN,, *HASAN GÖKTAŞ, 2,3 AND VOLKER J. SORGER 2 Department
More informationResilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study
Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:
More informationObjective: To introduce the equations of motion and describe the forces that act upon the Atmosphere
Objective: To introuce the equations of motion an escribe the forces that act upon the Atmosphere Reaing: Rea pp 18 6 in Chapter 1 of Houghton & Hakim Problems: Work 1.1, 1.8, an 1.9 on pp. 6 & 7 at the
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationInteraction force in a vertical dust chain inside a glass box
Interaction force in a vertical ust chain insie a glass box Jie Kong, Ke Qiao, Lorin S. Matthews an Truell W. Hye Center for Astrophysics, Space Physics, an Engineering Research (CASPER) Baylor University
More informationfv = ikφ n (11.1) + fu n = y v n iσ iku n + gh n. (11.3) n
Chapter 11 Rossby waves Supplemental reaing: Pelosky 1 (1979), sections 3.1 3 11.1 Shallow water equations When consiering the general problem of linearize oscillations in a static, arbitrarily stratifie
More informationMoving Charges And Magnetism
AIND SINGH ACADEMY Moving Charges An Magnetism Solution of NCET Exercise Q -.: A circular coil of wire consisting of turns, each of raius 8. cm carries a current of. A. What is the magnitue of the magnetic
More informationAn unstructured-grid, finite-volume model of the Irish Sea and Celtic Sea: Application to tides, surges, and temperature evolution.
An unstructure-gri, finite-volume moel of the Irish Sea an Celtic Sea: Application to ties, surges, an temperature evolution. Jiuing Xing, Alan avies, Phil Hall an Eric Jones National Oceanography Centre,
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationA SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical
More informationMoist Component Potential Vorticity
166 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 60 Moist Component Potential Vorticity R. MCTAGGART-COWAN, J.R.GYAKUM, AND M. K. YAU Department of Atmospheric an Oceanic Sciences, McGill University, Montreal,
More informationACCELERATION, FORCE, MOMENTUM, ENERGY : solutions to higher level questions
ACCELERATION, FORCE, MOMENTUM, ENERGY : solutions to higher level questions 015 Question 1 (a) (i) State Newton s secon law of motion. Force is proportional to rate of change of momentum (ii) What is the
More informationEstimation amount of snow deposits on the road
Estimation amount of snow eposits on the roa Olga. Glaysheva oronezh State University of Architecture an Civil Engineering, Russia Email: glaov@ox.vsi.ru ABSTRACT The article uner consieration gives the
More informationProblem 1 (20 points)
ME 309 Fall 01 Exam 1 Name: C Problem 1 0 points Short answer questions. Each question is worth 5 points. Don t spen too long writing lengthy answers to these questions. Don t use more space than is given.
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationDiagonalization of Matrices Dr. E. Jacobs
Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is
More informationP1D.6 IMPACTS OF THE OCEAN SURFACE VELOCITY ON WIND STRESS COEFFICIENT AND WIND STRESS OVER GLOBAL OCEAN DURING
P1D.6 IMPATS OF THE OEAN SURFAE VELOITY ON WIND STRESS OEFFIIENT AND WIND STRESS OVER GLOBAL OEAN DURING 1958-001 Zengan Deng 1* Lian Xie Ting Yu 1 an Kejian Wu 1 1. Physical Oceanography Laboratory, Ocean
More informationCenter of Gravity and Center of Mass
Center of Gravity an Center of Mass 1 Introuction. Center of mass an center of gravity closely parallel each other: they both work the same way. Center of mass is the more important, but center of gravity
More informationGoldstein Chapter 1 Exercises
Golstein Chapter 1 Exercises Michael Goo July 17, 2004 1 Exercises 11. Consier a uniform thin isk that rolls without slipping on a horizontal plane. A horizontal force is applie to the center of the isk
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationExperimental Determination of Mechanical Parameters in Sensorless Vector-Controlled Induction Motor Drive
Experimental Determination of Mechanical Parameters in Sensorless Vector-Controlle Inuction Motor Drive V. S. S. Pavan Kumar Hari, Avanish Tripathi 2 an G.Narayanan 3 Department of Electrical Engineering,
More informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
More informationS10.G.1. Fluid Flow Around the Brownian Particle
Rea Reichl s introuction. Tables & proofs for vector calculus formulas can be foun in the stanar textbooks G.Arfken s Mathematical Methos for Physicists an J.D.Jackson s Classical Electroynamics. S0.G..
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (II)
Laplace s Equation in Cylinrical Coorinates an Bessel s Equation (II Qualitative properties of Bessel functions of first an secon kin In the last lecture we foun the expression for the general solution
More informationMath 1271 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly
More informationNON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS
ANZIAM J. 46(2004), 1 16 NON-ADIABATIC COMBUSTION WAVES FOR GENERAL LEWIS NUMBERS: WAVE SPEED AND EXTINCTION CONDITIONS A. C. MCINTOSH 1,R.O.WEBER 2 ang.n.mercer 2 (Receive 14 August, 2002) Abstract This
More informationAN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.
AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD Mathca Release 14 Khyruin Akbar Ansari, Ph.D., P.E. Professor of Mechanical Engineering School of Engineering an Applie Science Gonzaga University SDC
More informationPHYS1169: Tutorial 8 Solutions
PHY69: Tutorial 8 olutions Wae Motion ) Let us consier a point P on the wae with a phase φ, so y cosϕ cos( x ± ωt) At t0, this point has position x0, so ϕ x0 ± ωt0 Now, at some later time t, the position
More informationJournal of Engineering Science and Technology Review 7 (1) (2014) Research Article
Jestr Journal of Engineering Science an Technology Review 7 () (04) 8 Research Article JOURNA OF Engineering Science an Technology Review www.jestr.org Research on the Measurement Error of MWIR Average
More informationParametrization of non-convective condensation processes May 1987
Parametrization of non-convective conensation processes May 987 By M. Tietke European Centre for Meium-Range Weather Forecasts Table of contents. Thermoynamics of moist air 2. Clou physical processes 2.
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More information