Reduced-Order Models of Zero-Net Mass-Flux Jets for Large-Scale Flow Control Simulations

Size: px
Start display at page:

Download "Reduced-Order Models of Zero-Net Mass-Flux Jets for Large-Scale Flow Control Simulations"

Transcription

1 Reuce-Orer Moels of Zero-Net Mass-Flux ets for Large-Scale Flow Control Simulations Reni Raju, Ehsan Aram an Rajat Mittal Department of Mechanical an Aerospace Engineering, The George Washington University, Washington, DC an Louis Cattafesta, Interisciplinary Microsystems Group Department of Mechanical an Aerospace Engineering, University of Floria, Gainesville, FL - Although computational tools are well suite for moeling the ynamics of zero-net mass-flux actuators, the computational costs involve in large-scale flow control simulations necessitate the use of reuce-orer moels for these evices. A new reuce-orer moel base on moeling only the slot of ZNMF jets in grazing flows is propose. A parametric stuy of imensionless parameters governing the characteristics of a ZNMF jet in grazing flow was conucte for the reuce-orer moel along with the commonly use sinusoial plug flow bounary conition an full cavity simulations. Comparisons of the current moel are presente in terms of the vortex ynamics an mean integral quantities, an the performance of the new reuce-orer moel is encouraging. In aition, the fielity of the moel has also been explore for a canonical separate flow. On the global scale of the flow, the reuce-orer moel is able to provie a more accurate representation of the full cavity simulation as compare to the conventional bounary conition. Nomenclature C () t = n moment of the jet, (m /s ) C () t = r moment of the jet, (m /s ) δ = Bounary layer thickness, (m) = Slot with/orifice iameter, (m) f = Separation bubble frequency, (Hz) f SL = Shear layer frequency, (Hz) h = Slot/Orifice height, (m) H = Cavity height, (m) L = Length of the prescribe velocity profile, (m) Ω () t = Instantaneous vorticity flux, (m /s ) = Phase angle, ( ) p = pressure, (N) Re δ = Bounary layer Reynols number, = U δ υ Re = et Reynols number, = V υ S = Stokes number, = πf υ Grauate Stuent, Department of Mechanical & Aerospace Engineering, Stuent Member AIAA. Grauate Stuent, Department of Mechanical & Aerospace Engineering, Stuent Member AIAA. Professor, Department of Mechanical & Aerospace Engineering, Associate Fellow AIAA, mittal@gwu.eu Professor, Department of Mechanical & Aerospace Engineering, Associate Fellow AIAA, cattafes@ufl.eu. American Institute of Aeronautics an Astronautics

2 St = Strouhal number, = π f u = Streamwise velocity (in x-irection), (m/s) U V U v V, top = Velocity amplitues, (m/s) = Freestream velocity, (m/s) = Cross-stream velocity (in y-irection), (m/s) = Spatially an temporally average jet velocity, (m/s) V = Prescribe velocity amplitue, (m/s) W x c z = With of the cavity, (m) = Center of the prescribe velocity profile, (m) ξ = Spanwise vorticity, (s - ) I. Introuction ERO-net mass-flux (ZNMF) actuators or synthetic jets have potential applications in the area of mixing Z enhancement, heat transfer,, mass transfer, jet vectoring, an active flow control of separation -8 an turbulence. 9 The ynamics an performance of these evices epen on several geometrical, structural an flow parameters. - When compare to the global omain, such as an airfoil, in which the actuator is imbee, the scales of the actuator are typically - - times smaller in size. Due to the range of scales involve, inclusion of a high-fielity moel of a ZNMF actuator within a macro-scale computational flow moel turns out to be an expensive, if not prohibitive, proposition. The esire to compute the flow physics associate with the ZNMF jet control makes reuce-orer moeling of these evices a practical necessity in these simulations. In the past these actuators have been represente as a reuce-orer moel in simulations in one form or another. Kral et al. moele the actuator via simplifie surface bounary conitions for velocity an pressure without the slot an cavity. For variations of the spatial istribution of the velocity profile, it was foun that the top-hat function provie the closest match to the experiments. Rizzetta et al. have use the recore flow fiel from a simulation of the isolate jet at the exit of a D perioic jet as a bounary conition to an external flow fiel. On the other han, Lockerby et al. have use a theoretical approach base on classic thin plate theory to moel the iaphragm eflection, while the slot is moele base on unsteay pipe-flow theory. A reuce-orer moel approximating D or D synthetic jets via quasi-d Euler equations was presente by Yamaleev & Carpenter. Filz et al. have moele D irecte synthetic jets using lumpe eterministic source terms (LDST) traine by a neural network. Rathnasingham & Breuer 8 have presente a semi-empirical analytical moel of ZNMF actuators using a system of couple nonlinear state equations escribing the structural an flui characteristics of the evice. An analytical lumpe element moel (LEM) of a piezoelectric-riven synthetic jet has been presente by Gallas et al. The moel represente the iniviual components of the synthetic jet as elements of an equivalent electric circuit using conjugate power variables. Sharma 9 has also use an alternate moel for LEM. More recently Tang et al. have compare the performance of the Dynamic Incompressible Flow moel (DI moel), Static Compressible Flow Moel (SC moel) an LEM in preicting the instantaneous space-average velocity at the orifice exit. Loworer moeling of two-imensional synthetic jets via Proper Orthogonal Decomposition (POD) was carrie out by Reiniotis et al. The ynamical moel base on Galerkin projection was erive from the flow for specific Reynols an Stokes numbers. Kihwan et al. have evelope a ynamical moel base on system ientification to ientify the interaction of synthetic jets with a laminar bounary layer with potential application to feeback control. The top-hat profile use by Kral et al. is the simplest moel of the jet an is referre to as the zeroth-orer moel in the current stuy. Such a representation eliminates the nee of simulating the cavity an the slot. A similar approach has also been aopte in some previous stuies., However, it is not clear to what extent such a bounary conition accurately moels the key features of the jet. Since the control mechanism of the actuator epens on transfer of the momentum or vorticity flux to the external bounary layer, it is essential that any such moel shoul be able to represent these quantities accurately. The purpose of the current stuy is to explore the issues of actuator moeling via numerical simulations. Base on these simulations we propose a new moel for ZNMF jets in grazing flows. The moel is compare with the ful cavity simulations an the zeroth-orer moel in attache flows. In aition, the effectiveness of the propose moel with the moifie bounary conition representation versus full cavity simulations has also been stuie for a separate bounary layer in a canonical flow configuration. American Institute of Aeronautics an Astronautics

3 II. Slot-Only ZNMF Actuator Moel The major reason for the shortcoming of the moifie bounary conition moels is the neglect of the flow physics associate with the slot. Our previous stuies show that flow in the slot tens to separate at the top an bottom lips, an the formation of seconary vortices near the exit tens to significantly alter the flow fiel. Thus the working hypothesis for the current reuce-orer moel is that the inclusion of just the slot with appropriate bounary conitions shoul significantly improve the fielity of the actuator moel. This hypothesis is supporte by past stuies, that have shown that the shape of the cavity has very little effect on the flow emanating from the jet provie that the incompressible flow assumption insie the actuator cavity is vali. The approach use for the current stuy explores the significance of the slot via simplifie bounary conitions an hence effectively reuces the computational complexity. An attempt is thus mae to replicate the flow physics insie the slot to further improve the fielity of the moel. Consier the flow insie the D slot of a typical ZNMF actuator as seen in Figure (a). The area change from the cavity to the slot causes the flow to turn near the bottom lip of the slot, an the flow pattern uring expulsion is similar to the one cause by a sink present somewhere along the slot center uring the expulsion phase. Hence the flow enters the slot raially at any given time uring expulsion an, in terms of its Cartesian components, will have two components of velocity, u an v. As seen in Figure (a), at the entrance to the slot the u -velocity will be significant, while near the slot center the v -velocity will ominate. Thus an approximation of this flow pattern prescribes a corresponing bounary conition profile for the u - velocity which varies linearly in the x-irection from a maximum positive value near the left wall to a minimum negative value near the right wall (see Table ). (a) (b) (c) Figure : Schematic showing (a) the typical flow pattern seen in a ZNMF jet slot for a full cavity () configuration an (b) moifie bounary conition (MBC) an (c) Propose slot-only (SO) moel. A. Moel Comparison Base on the above arguments, three separate flow configurations have been chosen for comparison. These inclue the full-cavity (), the slot-only (SO) an moifie bounary conition (MBC) configurations. As the name suggests, the moel consists of full-cavity simulations an is a complete representation of a D ZNMF jet, as seen in Figure (a). The river isplacement is represente as an oscillatory bounary conition at the bottom of the cavity an has been foun to be a fairly accurate representation when compare with experimental results. On American Institute of Aeronautics an Astronautics

4 the other han Figure (b) shows the MBC configuration, which consists of prescribing only the vertical component of the velocity over the exit plane of the jet, as mentione earlier. The propose SO moel is base on the configuration shown in Figure (c). For this configuration the slot is moele by consiering the full length of the slot, i.e., h=%, which is necessary for the replicating the sink type behavior of a configuration. Two variants of the slot-only moels are consiere. For the first, calle the SO- moel, only the sinusoial v-velocity is prescribe at the bottom while for the secon, the SO- moel, bounary conitions for both the u- & v-velocities uring expulsion have been use. Here the u-velocity is prescribe as a linearly varying sinusoial profile that satisfies no-slip conition at the bounary. Table lists the bounary conitions use for the ifferent configurations uring both expulsion an ingestion phases with the corresponing location. Table : Bounary conitions use for the ifferent configurations consiere for current work. y Expulsion Ingestion Moel location u vxt (, ) uxt (, ) vxt (, ) ( xt, ) w.r.t to y V sin( ωt) V sin( ω t) - MBC h+ H V sin( ωt) V sin( ω t) SO- H V sin( ωt) V sin( ω t) SO- H V sin( ωt) x U sin( ) ω t, for < x < V sin(, at x = & Care has been taken to ensure that the mass flow rate at the exit of the slot for the configuration an reuceorer moels are the same. The imensionless parameters of significance for a ZNMF jet uner quiescent conition are the Strouhal number, St, jet Reynols number, Re, velocity ratio, U an bounary-layer thickness to jet with ratio, δ. Table lists the cases along with the corresponing parameters use for the current stuy. It is known that the jet characteristics are etermine by functional epenence on the parameters as U δ fn St,,, Re j Vj. Hence for the current stuy three ifferent simulation sets were consiere for each of these parameters while keeping the rest constant. The velocity amplitue ratio U V = was maintaine for all the cases. Table : Parameters use for reuce-orer moeling simulations Case Re S U V δ Re δ St American Institute of Aeronautics an Astronautics

5 All the computations use a single gri of imensions 8 8 for the simulations an a omain size of 9. For the remaining configurations the gri resolution in the slot an external flow fiel is hel constant at all times. The irect numerical simulations are use to moel synthetic jets issuing from a cavity. Incompressible Navier-Stokes equations are solve in time using a secon-orer accurate fractional step metho. In a Cartesian framework, a secon-orer Aams-Bashforth scheme is employe for the convective terms, while the iffusion terms are iscretize using an implicit Crank-Nicolson scheme that eliminates the viscous stability constraint. The solver has been rigorously valiate by comparisons of several test cases against establishe experimental an computational ata incluing synthetic jets.,8 B. Canonical Separate Flow Configuration Whereas the above stuy examines the performance of the ifferent moes for an attache grazing flow, a relevant issue to be examine is to what extent the etails of the actuator moel moify the effect of the jet on a separation bubble. In orer to aress this issue we have also examine a case where a separation bubble is create ownstream of the jet by prescribing a blowing-suction bounary conition on the top bounary of the omain. A zero-vorticity conition, along the lines of Na & Moin 9 is applie on this bounary as follows: u G vxl (, y ) = Gx ( ), = () y x where Gx ( ) ( xl, y ) is the prescribe blowing an suction velocity profile efine in the form of: π( x xc ) Gx ( ) = Vtop sin e L β ( x x c ) α L where L is the length an xc is the center of the velocity profile. The parameters Vtop, α an β are set to U, an, respectively. Figure shows the configuration use to generate the separate flow. The, SO- an MBC moel configurations were teste uner these conitions. The simulations have been carrie out on the omain with size using a gri of imensions for case an same resolution for two other configurations. () Figure : Flow configuration use for simulating a canonical separate flow. III. Results A. Moel Performance in Attache Grazing Flow Holman et al. have shown that the jet formation criterion is etermine by the relation, Re S = St>K an is proportional to the non-imensional vorticity flux. Here, K is a constant foun to be O() for D jets an ~. for axisymmetric jets. Thus for the given set of conitions of St =.8, a jet is forme an is able to impart a significant amount of vorticity flux to external flow. Figure shows the comparison of the spanwise vorticity at the peak expulsion ( =9 ) an peak ingestion ( = ) for all moels at Re =, S =, U =. an δ =.. The flow tens to separate near the bottom lip of the slot for the full cavity simulations. The separating shear layers ten to roll up on either wall uring expulsion at the exit plane of the slot an a strong clockwise vortex American Institute of Aeronautics an Astronautics

6 is forme that is swept away by the incoming bounary layer. The MBC moel, seen in Figure (b), forms a smaller clockwise vortex near the lip of the slot in comparison to the moel. The SO moels are able to reprouce the flow physics of the baseline flow in the proximity of the slot exit, as seen in Figure (c) & (). Unlike the moel, insie the slot the bounary layer remains attache to the walls for SO- moel. On the other han the SO- oes show flow separation on the slot walls but is unable to capture roll up of the bounary layer on the walls. This ifference can be attribute to larger incoming angle of the flow from the cavity for the moel. The ingestion phase for the SO moels shows a near perfect match with the moel within the slot, proving the valiity of using a simple Neumann bounary conition for the u-velocity. The SO moels are able to capture even the seconary vortex prouce uring expulsion, although for the SO- moel the vortex loses its strength ownstream. ~9 ~ (a) (b) (c) () Figure : Instantaneous spanwise vorticity uring peak expulsion an peak ingestion for (a), (b) MBC moel (c) SO- moel an () SO- moel for Re =, S =, U =. an δ =.. Dashe lines compare the locations of vortex structures for the three moels with respect to configuration. American Institute of Aeronautics an Astronautics

7 The characteristics of the jet can be therefore efine on the basis of the integral quantities such as vorticity flux, momentum flux an kinetic energy flux. The unsteay variations of these quantities are compare for the representative cases in Figure. Note the efinitions for momentum flux an kinetic energy flux are ρ C () t v( x,) t x = ρ [ ] an ρ () ρ [ (,)] C t = v x t x an have been normalize by a suitable quantity. On the other han the vorticity flux is efine as () t ξ ( x,) t v( x, t) Ω = x. Figure shows that the accumulation of z vorticity from shear layer separation causes peaks in the integral quantities for the moel uring expulsion. Similar features are seen for the SO- moel an, although they are lower in magnitue, they appear to give a closer match to the moel than the SO- moel uring expulsion. On the other han the MBC moel is unable to preict the flow behavior in the slot uring most of the cycle. Hence although the flow fiel is significantly altere by seconary vortices, the integral measures can be reasonably approximate by using the SO- moel; although for a better match the incoming flow angle, i.e. U V ratio, might nee to be moifie further. SO- Moel SO- Moel SO- Moel - ρc /ρv ρc /ρv Figure : Comparison of the vorticity flux (left), momentum flux (mile) an kinetic energy flux (right) for all moels as a function of phase for Re =, S =, U =., δ =., St =.8. Simulations at higher Strouhal number ( St =.) are presente in Figure where Re =, S =, U V =. an δ =.. This correspons to the no-jet formation case ue to the high St, an uring peak expulsion the bounary layers near the bottom lips thicken to form a vena-contracta. At the exit, the shear layer rolls up to form a relatively small clockwise vortex which oes not travel far from the slot uring the suction stroke. Note that the expulsion stroke is efine as o < 8 o o an suction stroke as 8 < o. Interestingly it is foun that all the moels give a reasonable approximation to the moel at the slot exit uring the expulsion phase of the cycle. However, the MBC moel oes not prouce coherent vortex structures as seen for the rest of the configurations near the right exit lip which is apparent uring the suction stroke, as seen in Figure (b). The SO- moel captures the vena-contracta of the flow near the slot inlet which, although not seen for the SO- moel, oes not seem to significantly affect the external flowfiel. Both moels replicate the flow behavior of moel uring the suction stroke. American Institute of Aeronautics an Astronautics

8 ~9 ~ (a) (b) (c) () Figure : Instantaneous spanwise vorticity uring peak expulsion an peak ingestion for (a), (b) MBC moel (c) SO- moel an () SO- moel at Re =, S =, U =., δ =.. Dashe lines compare the locations of vortex structures for the three moels with respect to configuration. At higher Strouhal numbers, as is the case for Re =, S =, U V =. an δ =., seen in Figure, the integral measures o not show peaks uring the expulsion cycle. An increase in St tens to ecrease the contribution of vorticity flux of the jet to the external bounary layer. This also remains true for the momentum an kinetic energy fluxes. Both the SO moels compare reasonably well with the moel uring the whole cycle. On the other han, the MBC moel shows a slight phase shift uring the expulsion phase for the vorticity flux, while it prouces slightly larger eviations for the momentum an kinetic energy fluxes. 8 American Institute of Aeronautics an Astronautics

9 SO- Moel SO- Moel SO- Moel - ρc /ρv ρc /ρv Figure : Comparison of the vorticity flux (left), momentum flux (mile) an kinetic energy flux (right) for all moels as a function of phase for Re =, S =, U =., δ =., St =.. Similar analyses have been conucte for the other cases. Figure (a) an (b) compare the integral quantities for Re =, S =, U V =., δ =. an Re =8., S =, U =., δ =., respectively, where St =.8 for both the cases. When compare with Re =, S =, U =., δ =. shown in Figure, it can be seen that, although the Strouhal number remains the same, as Re increases the vorticity flux uring the expulsion phase also tens to increase an istort the sinusoial profile. These peaks correspon to the expulsion of vortices from the exit plane of the jet. On the other han the suction phase behavior remains the same for these cases. In terms of moel performance for Re =, S =, U =. an δ =. the SO- moel tens to slightly overpreict the integral measures versus the moel uring the expulsion phase, while the SO- moel slightly unerpreicts these values. Surprisingly the MBC moel also shows a reasonable match for the vorticity flux uring the expulsion phase but not for the other quantities. During ingestion the SO moels match the quantities at all times. However, this is not the case for the MBC moel. For Re =8., S =, U =. an δ =., similar behavior is observe with the SO- moel proviing the best approximation. At higher Strouhal numbers, as is the case for Re =, S =, U =. an δ =., seen in Figure (c) where St =., the integral measures show slightly lower values. Note that for this case the vorticity flux again shows several peaks uring the expulsion phase inicating that the jet is forme uner these conitions. The integral quantities of both SO moels compare reasonably well with the moel uring the whole cycle an are able to capture the small peaks in the vorticity flux. Note that for this case, the behavior of both moel is similar, inicating that uner these conitions the shear layer separation at the bottom lip is not as important as the cases with St =.8 where the SO- moel oes not perform as well as the SO- moel. This is apparent from Figure where both these moels show similar representation of the flow physics. In aition to the significance of Re an St, the effects of freestream to jet velocity ratio, U V an the bounary layer thickness to jet iameter ratio, δ are presente in Figure 8. The variation in the velocity ratio is presente in the Figure 8(a) an (b) for Re =, S =, U =., δ =. an Re =, S =, U V =., δ =. respectively. Due to the lower freestream velocity, the magnitues of the peaks are reuce. It can be observe that increasing the velocity ratio leas to an increase in the momentum an kinetic energy fluxes. For these cases the SO- moel is able to yiel a better approximation to the configuration, while both the SO- an MBC moels significantly unerpreict the moments uring the expulsion phase. On the other han, ecreasing the bounary layer thickness significantly alters the integral quantities an leas to an increase in the integral measures over a cycle with sharper peaks as seen in Figure 8(c) an () for two ifferent ratios δ =. an δ =., respectively. The preictive capabilities of the moels eteriorate as the bounary layer thickness is reuce, although the SO- moel consistently gives better performance over other moels. Overall the ingestion phase of both SO moels compares well with the moel. This implies that the inclusion of the slot is a significant factor for improvement of the moel, while accounting for the separation at the bottom lip of the slot increases the fielity of the moel still further. 9 American Institute of Aeronautics an Astronautics

10 SO- Moel SO- Moel SO- Moel - ρc /ρv ρc /ρv (a) (b) SO- Moel SO- Moel ρc /ρv ρc /ρv SO- Moel SO- Moel ρc /ρv ρc /ρv SO- Moel SO- Moel (c) Figure : Comparison of the vorticity flux (left), momentum flux (mile) an kinetic energy flux (right) for all moels as function of phase for (a) Re =, S =, U =., δ =., St =.8, (b) Re =8., S =, =., =., St =.8, an (c) Re =, S =, =., δ =., St =.. U δ U American Institute of Aeronautics an Astronautics

11 SO- Moel SO- Moel SO- Moel - ρc /ρv ρc /ρv (a) SO- Moel SO- Moel SO- Moel - ρc /ρv ρc /ρv (b) SO- Moel ρc /ρv SO- Moel ρc /ρv SO- Moel (c) SO- Moel SO- Moel 8 SO- Moel ρc /ρv ρc /ρv () Figure 8: Comparison of the vorticity flux (left), momentum flux (mile) an kinetic energy flux (right) for all moels as function of phase for (a) Re =, S =, U =., δ =., St =.8, (b) Re = S =, U δ U =., =., St =.8, (c) Re =, S =, =., δ =., St =.8 an () Re =, S =, U V =., δ =., St =.8. American Institute of Aeronautics an Astronautics

12 Figure 9: Spanwise vorticity an mean streamlines for the unforce separate flow. B. Moel Performance for a Canonical Separate Flow As seen in Figure 9, a separation bubble was create by imposing the bounary conition escribe in Eq. ()- () with length, L sep =. The effect of forcing of the ZNMF jet, place approximately. upstream of the point of separation, on the separation bubble has been examine for ifferent moel configurations. As seen in the previous section, since the SO- moel performs superior to other moels, it was teste along with an MBC configurations for six ifferent forcing frequencies base on shear layer frequency, such that f = ε f SL. The six values of ε were chosen as.,.,.,,. an.. For all the cases examine for the separate flow, the rest of the flow characteristics were hel fixe at Re =, U V = an δ =. Figure shows contours of instantaneous spanwise vorticity at peak expulsion an mean streamlines for these three cases at one particular forcing frequency ( f =. fsl ). By comparing the mean streamlines it can be seen that an SO- configuration yiel a very similar effect on the separation bubble. For both these cases the size of the separation bubble is reuce by a factor of.. On the other han the MBC moel, as expecte, is unable to match the moel, where the separation bubble size is reuce by a factor of. It is clear that the SO- moel is able to able to provie a better representation of the ZNMF actuator on a global scale. (a) (b) (c) Figure : Instantaneous spanwise vorticity an mean streamlines for the a), b) SO-, an c) MBC moels for f =. f. SL American Institute of Aeronautics an Astronautics

13 Figure compares the performance of three moels for ifferent forcing frequencies. The separation bubble frequency in this figure is efine, base on the separation bubble size, as f = U / L. It is foun that the SO- moel gives a goo approximation to the moel for separation bubble control at all forcing frequencies, especially for ε.. The maximum percentage ifference of separation bubble size seen between SO- an moel is approximately 8%. On the other han the MBC moel unerestimates the performance of on separation control for low frequencies ( f < fsl ) where control is most effective. Overall the percentage ifference of separation bubble size for the MBC moel varies from %. It is also clear that the separation bubble size is reuce by ecreasing the jet frequency, an maximum reuction in this particular problem an this frequency range is approximately %. sep sep Figure : Effect of forcing frequency on the separation bubble size for three moels. IV. Conclusions A simple reuce-orer moel is presente for the representation of ZNMF jets in large-scale flow control simulations. The moel inclues the slot of the actuator ue to its importance in governing the ynamics of the interaction process with a grazing flow. Two variants of the slot-only moel were consiere, an a parametric numerical stuy was carrie out to etermine the performance of each moel compare to a conventional moifie bounary conition. It was foun that, unlike the conventional MBC moel, the slot-only moels were able to capture more of the flow physics associate with full cavity simulations. A comparison of the integral measures of vorticity, momentum, an kinetic energy fluxes for these moels showe that the slot-only moel that assumes a sink-like flow at the slot inlet (i.e., the SO- moel) provies the best moeling approximation. The performance of the moel was also compare with full cavity simulations an the simple MBC moel for a canonical separate flow at ifferent forcing frequencies. It was foun that at lower forcing frequencies, the flow separation extent was reuce significantly, an the SO- moel was able to preict the separation bubble size reuction with improve accuracy. Future work in this irection will focus on etermining the appropriate value of U V, use in SO-, as a function of the imensionless flow parameters to stuy the effect of incoming flow angle at the bottom lip of the slot. Acknowlegments This work is supporte by grants from AFOSR an a NASA Cooperative Agreement NNXAD9A monitore by Dr. Brian Allan. American Institute of Aeronautics an Astronautics

14 References Chen, Y., Liang, S., Aung, K., Glezer, A., an agoa,, "Enhance Mixing in a Simulate Combustor Using Synthetic et Actuators." AIAA Paper 99-9, 999. Campbell,.S., Black, W.Z., Glezer, A., an Hartley,.G., "Thermal Management of a Laptop Computer with Synthetic Air Microjets." Intersociety Conference on Therm. Phenomenon, IEEE, 998, pp. -. Mahalingam, R., an Glezer, A., "Design an Thermal Characteristics of a Synthetic et Ejector Heat Sink." ournal of Electronic Packaging, Vol., No.,, pp. -. Trávníček, Z., an Tesař, V., "Annular Synthetic et Use for Impinging Flow Mass Transfer." International ournal of Heat an Mass Transfer, Vol., No.,, pp Smith, B.L., an Glezer, A., "et Vectoring Using Synthetic ets." ournal of Flui Mechanics, Vol. 8,, pp. -. Seifert, A., Bachar, T., Koss, D., Shepshelovich, M., an Wygnanski, I., "Oscillatory Blowing: A Tool to Delay Bounary-Layer Separation." AIAA ournal, Vol., No., 99, pp. -. Seifert, A., Darabi, A., an Wygnanski, I., "Delay of Airfoil Stall by Perioic Excitation." ournal of Aircraft, Vol., No., 99, pp Honohan, A.M., Amitay, M., an Glezer, A., "Aeroynamic Control Using Synthetic ets." AIAA Paper -,. 9 Rathnasingham, R., an Breuer, K.S., "System Ientification an Control of a Turbulent Bounary Layer." Physics of Fluis, Vol. 9, No., 99, pp Gallas, Q., Holman, R., Nishia, T., Carroll, B., Sheplak, M., an Cattafesta, L., "Lumpe Element Moeling of Piezoelectric-Driven Synthetic et Actuators." AIAA ournal, Vol., No.,, pp. -. Glezer, A., an Amitay, M., "Synthetic ets." Annual Review of Flui Mechanics, Vol.,, pp. -9. Holman, R., Utturkar, Y., Mittal, R., Smith, B.L., an Cattafesta, L., "Formation Criterion for Synthetic ets." AIAA ournal, Vol., No.,, pp. -. Kral, L.D., Donovan,.F., Cain, A.B., an Cary, A.W., "Numerical Simulation of Synthetic et Actuators." AIAA Paper 9-8, 99. Rizzetta, D.P., Visbal, M.R., an Stanek, M.., "Numerical Investigation of Synthetic-et Flow Fiels." AIAA ournal, Vol., No. 8, 999, pp Lockerby, D. A., Carpenter, P.W., an Davies, C., "Numerical Simulation of the Interaction of Microactuators an Bounary Layers." AIAA ournal, Vol., No.,, pp. -. Yamaleev, N.K., an Carpenter, M.H., "A Reuce-Orer Moel for Efficient Simulation of Synthetic et Actuators." NASA Technical Reports, NASA/TM--,. Filz, C., Lee, D., Orkwis, P.D., an Turner, M.G., "Moeling of Two Dimensional Directe Synthetic ets Using Neural Network-Base Detereministic Source Terms." AIAA Paper -,. 8 Rathnasingham, R., an Breuer, K.S., "Couple Flui-Structural Characteristics of Actuators for Flow Control." AIAA ournal, Vol., No., 99, pp Sharma, R. N., "Flui-Dynamics-Base Analytical Moel for Synthetic et Actuation." AIAA ournal, Vol., No. 8,, pp Tang, H., Zhong, S., abbal, M., Garcillan, L., Guo, F., Woo, N.., an Warsop, C., "Towars the Design of Synthetic-et Actuators for Full-Scale Flight Conitions. Part : Low-Dimensional Atuator Preiction Moels an Actuator Design Methos." Flow, Turbulence an Combustion, Vol. 8, No.,, pp Reiniotis, O.K., Ko,., an Kurila, A.., "Reuce Orer Nonlinear Navier-Stokes Moels for Synthetic ets." ournal of Fluis Engineering, Vol., No.,, pp. -. Kihwan, K., Beskok, A., an ayasuriya, S., "Nonlinear System Ientification for the Interaction of Synthetic ets with a Bounary Layer." American Control Conference,. Proceeings of the,. Rampunggoon, P., "Interaction of a Synthetic et with a Flat-Plate Bounary Layer." Ph Thesis, Department of Mechanical Engineering, University of Floria,. Ravi, B.R., "Numerical Stuy of Three Dimensional Synthetic ets in Quiescent an External Grazing Flows." DSc Thesis, Mechanical an Aerospace Engineering, The George Washington University,. Raju, R., Mittal, R., Gallas, Q., an Cattafesta, L., "Scaling of Vorticity Flux an Entrance Length Effects in Zero-Net Mass-Flux Devices." AIAA Paper -,. Kotapati, R.B., Mittal, R., an Cattafesta, L., "Numerical Stuy of Transitional Synthetic et in Quiescent External Flow." ournal of Flui Mechanics, Vol. 8,, pp. 8-. Utturkar, Y., Mittal, R., Rampunggoon, P., an Cattafesta, L., "Sensitivity of Synthetic ets to the Design of the et Cavity." AIAA Paper -,. 8 Mittal, R., Dong, H., Bozkurttas, M., Najjar, F., Vargas, A., an Loebbecke, A., "A Versatile Immerse Bounary Metho for Incompressible Flows with Complex Bounaries." ournal of Computational Physics (uner review),. 9 Na, Y., an Moin, P., "Direct Numerical Simulation of a Separate Turbulent Bounary Layer." ournal of Flui Mechanics, Vol., 998, pp. -. American Institute of Aeronautics an Astronautics

Simple Representations of Zero-Net Mass-Flux Jets in Grazing Flow for Flow-Control Simulations

Simple Representations of Zero-Net Mass-Flux Jets in Grazing Flow for Flow-Control Simulations Simple Representations of Zero-Net Mass-Flux ets in Grazing Flow for Flow-Control Simulations by Ehsan Aram, Rajat Mittal an Louis Cattafesta Reprinte from International ournal of Flow Control olume 2

More information

Toward Simple Boundary Condition Representations of Zero-Net Mass-Flux Actuators in Grazing Flow

Toward Simple Boundary Condition Representations of Zero-Net Mass-Flux Actuators in Grazing Flow 39th AIAA Fluid Dynamics Conference - 5 une 9, San Antonio, Texas AIAA 9-48 Toward Simple Boundary Condition Representations of Zero-Net Mass-Flux Actuators in Grazing Flow Ehsan Aram, Rajat Mittal Department

More information

Simple Models of Zero-Net Mass-Flux Jets for Flow Control Simulations

Simple Models of Zero-Net Mass-Flux Jets for Flow Control Simulations Simple Models of Zero-Net Mass-Flux Jets for Flow Control Simulations by Reni Raju, Ehsan Aram, Rajat Mittal and Louis Cattafesta Reprinted from International Journal of Flow Control Volume Number September

More information

Low Dimensional Modeling of Zero-Net Mass-Flux Actuators

Low Dimensional Modeling of Zero-Net Mass-Flux Actuators Low Dimensional Moeling of Zero-Net Mass-Flux Actuators Quentin Gallas * an Ryan olman * University of Floria, Gainesville, Fl, 368 Reni Rau * an Raat Mittal The George Washington University, Washington

More information

STUDY OF THREE-DIMENSIONAL SYNTHETIC JET FLOWFIELDS USING DIRECT NUMERICAL SIMULATION.

STUDY OF THREE-DIMENSIONAL SYNTHETIC JET FLOWFIELDS USING DIRECT NUMERICAL SIMULATION. 42 nd AIAA Aerospace Sciences Meeting and Exhibit 5-8 January 2004/Reno, NV STUDY OF THREE-DIMENSIONAL SYNTHETIC JET FLOWFIELDS USING DIRECT NUMERICAL SIMULATION. B.R.Ravi * and R. Mittal, Department of

More information

Numerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number

Numerical 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 information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE 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 information

A Jet Formation Criterion for Synthetic Jet Actuators

A Jet Formation Criterion for Synthetic Jet Actuators AIAA 3-636 A Jet Formation Criterion for Synthetic Jet Actuators Yogen Utturkar, Ryan Holman, Rajat Mittal, Bruce Carroll, Mark Sheplak, and Louis Cattafesta University of Florida Gainesville, FL The George

More information

LES Study of Synthetic Jet Frequency and Amplitude Effects on a Separated Flow

LES Study of Synthetic Jet Frequency and Amplitude Effects on a Separated Flow 1 LES Study of Synthetic Jet Frequency and Amplitude Effects on a Separated Flow ONERA Applied Aerodynamic Department Pierre-Yves Pamart, Julien Dandois, Eric Garnier, Pierre Sagaut 2 Introduction Separated

More information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow 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 information

The effect of nonvertical shear on turbulence in a stably stratified medium

The 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 information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

A SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS

A 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 information

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić

CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS. Jože Kutin *, Ivan Bajsić Sensors an Actuators A: Physical 168 (211) 149-154 oi: 1.116/j.sna.211..7 211 Elsevier B.V. CHARACTERISTICS OF A DYNAMIC PRESSURE GENERATOR BASED ON LOUDSPEAKERS Jože Kutin *, Ivan Bajsić Laboratory of

More information

LQG FLUTTER CONTROL OF WIND TUNNEL MODEL USING PIEZO-CERAMIC ACTUATOR

LQG 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 information

Numerical study of a transitional synthetic jet in quiescent external flow

Numerical study of a transitional synthetic jet in quiescent external flow J. Fluid Mech. (7), vol. 581, pp. 87 31. c 7 Cambridge University Press doi:1.117/s117564 Printed in the United Kingdom 87 Numerical study of a transitional synthetic jet in quiescent external flow RUPESH

More information

To understand how scrubbers work, we must first define some terms.

To 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 information

Evaluation of Column Breakpoint and Trajectory for a Plain Liquid Jet Injected into a Crossflow

Evaluation of Column Breakpoint and Trajectory for a Plain Liquid Jet Injected into a Crossflow ILASS Americas, 1 st Annual Conference on Liqui Atomization an Spray Systems, Orlano, Floria, May 008 Evaluation of Column Breakpoint an Trajectory for a Plain Liqui Jet Injecte into a Crossflow S.M. Thawley,

More information

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains

Hyperbolic Systems of Equations Posed on Erroneous Curved Domains Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (

More information

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes

Strength Analysis of CFRP Composite Material Considering Multiple Fracture Modes 5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation

More information

Stable and compact finite difference schemes

Stable and compact finite difference schemes Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long

More information

6. Friction and viscosity in gasses

6. 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 information

Influence the Nozzle Shape on Local Heat Transfer in Impinging Jet

Influence the Nozzle Shape on Local Heat Transfer in Impinging Jet Issue 6, Volume 6, 12 Influence the Nozzle Shape on Local Heat Transfer in Impinging Jet M. Attalla an M. S. Ahme 1 Abstract The local Nusselt number istributions of circular nozzle on a heate flat plate

More information

Table of Common Derivatives By David Abraham

Table 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 information

Dusty Plasma Void Dynamics in Unmoving and Moving Flows

Dusty Plasma Void Dynamics in Unmoving and Moving Flows 7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Dusty Plasma Voi Dynamics in Unmoving an Moving Flows O.V. Kravchenko*, O.A. Azarova**, an T.A. Lapushkina*** *Scientific an Technological

More information

Two- and Three-Dimensional Validation of Icing Model

Two- and Three-Dimensional Validation of Icing Model APCOM & ISCM 11-14 th December, 13, Singapore To- an Three-Dimensional Valiation of Icing Moel *Ryosuke Hayashi¹ an Makoto Yamamoto 1 Grauate School of Mechanical Engineering, Tokyo University of Science

More information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal 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 information

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite

More information

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.

More information

A new identification method of the supply hole discharge coefficient of gas bearings

A new identification method of the supply hole discharge coefficient of gas bearings Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico

More information

Problems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract

Problems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we

More information

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract

More information

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges

A Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

1 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.

1 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 information

05 The Continuum Limit and the Wave Equation

05 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 information

Study on aero-acoustic structural interactions in fan-ducted system

Study on aero-acoustic structural interactions in fan-ducted system Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,

More information

SIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES

SIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES SIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES Jan Macek, Miloš Polášek Czech Technical University in Prague, Josef Božek Research Center Introuction Improvement of emissions from reciprocating internal

More information

Evaporating droplets tracking by holographic high speed video in turbulent flow

Evaporating 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 information

Situation awareness of power system based on static voltage security region

Situation 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 information

Impurities in inelastic Maxwell models

Impurities in inelastic Maxwell models Impurities in inelastic Maxwell moels Vicente Garzó Departamento e Física, Universia e Extremaura, E-671-Baajoz, Spain Abstract. Transport properties of impurities immerse in a granular gas unergoing homogenous

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation

Electromagnet Gripping in Iron Foundry Automation Part II: Simulation www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper

More information

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section

The influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: 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 information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 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 information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers

Optimal Variable-Structure Control Tracking of Spacecraft Maneuvers Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The 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 information

V = Flow velocity, ft/sec

V = Flow velocity, ft/sec 1 Drag Coefficient Preiction Chapter 1 The ieal force acting on a surface positione perpenicular to the airflow is equal to a ynamic pressure, enote by q, times the area of that surface. Dynamic pressure

More information

Efficient Macro-Micro Scale Coupled Modeling of Batteries

Efficient 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 information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Hilton C. de M. Mello University of São Paulo School of Engineering, São Carlos São Carlos, SP.

Hilton C. de M. Mello University of São Paulo School of Engineering, São Carlos São Carlos, SP. Hilton C. de M. Mello et al Hilton C. de M. Mello hcmmello@sc.usp.br University of São Paulo School of Engineering, São Carlos 13566-590 São Carlos, SP. Brazil Fernando M. Catalano catalano@sc.usp.br Leandro

More information

Homework 7 Due 18 November at 6:00 pm

Homework 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 information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE 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 information

On Using Unstable Electrohydraulic Valves for Control

On Using Unstable Electrohydraulic Valves for Control Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves

More information

IPMSM Inductances Calculation Using FEA

IPMSM Inductances Calculation Using FEA X International Symposium on Inustrial Electronics INDEL 24, Banja Luka, November 68, 24 IPMSM Inuctances Calculation Using FEA Dejan G. Jerkan, Marko A. Gecić an Darko P. Marčetić Department for Power,

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining 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 information

EXPERIMENTAL INVESTIGATION ON PNEUMATIC COMPONENTS

EXPERIMENTAL INVESTIGATION ON PNEUMATIC COMPONENTS Conference on Moelling Flui Flow (CMFF 03) The 12 th International Conference on Flui Flow Technologies Buapest, Hungary, September 3-6, 2003 EXPERIMENTAL INVESTIGATION ON PNEUMATIC COMPONENTS Zoltán MÓZER,

More information

Damage detection of shear building structure based on FRF response variation

Damage detection of shear building structure based on FRF response variation , pp.18-5 http://x.oi.org/10.1457/astl.013.3.05 Damage etection of shear builing structure base on FRF response variation Hee-Chang Eun 1,*, Su-Yong Park 1, Rae-Jung im 1 1 Dept. of Architectural Engineering,

More information

Simulation of Angle Beam Ultrasonic Testing with a Personal Computer

Simulation 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 information

Quantum Mechanics in Three Dimensions

Quantum 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 information

Dynamics of the synchronous machine

Dynamics of the synchronous machine ELEC0047 - Power system ynamics, control an stability Dynamics of the synchronous machine Thierry Van Cutsem t.vancutsem@ulg.ac.be www.montefiore.ulg.ac.be/~vct These slies follow those presente in course

More information

Interaction force in a vertical dust chain inside a glass box

Interaction 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 information

MAE 210A FINAL EXAM SOLUTIONS

MAE 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 information

MULTI-SCALE METHODS FOR THE SOLUTION OF THE RADIATIVE TRANSFER EQUATION

MULTI-SCALE METHODS FOR THE SOLUTION OF THE RADIATIVE TRANSFER EQUATION MULTI-SCALE METHODS FOR THE SOLUTION OF THE RADIATIVE TRANSFER EQUATION Pero J. Coelho *,, Nicolas Crouseilles **, Pero Pereira * an Maxime Roger *** * LAETA, IDMEC, Dept. of Mechanical Engineering, Instituto

More information

SPE Copyright 1999, Society of Petroleum Engineers Inc.

SPE Copyright 1999, Society of Petroleum Engineers Inc. SPE 664 Effect of Flow Through a Choke Valve on Emulsion Stability M.J. van er Zane, SPE, K.R. van Heuven, J.H. Muntinga, SPE, an W.M.G.T. van en Broek, SPE, Delft University of Technology Copyright 1999,

More information

Domain Decomposition and Model Reduction of Systems with Local Nonlinearities

Domain Decomposition and Model Reduction of Systems with Local Nonlinearities Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Abstract The goal of this paper is to combine balance truncation

More information

Sparse Reconstruction of Systems of Ordinary Differential Equations

Sparse 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 information

Approaches for Predicting Collection Efficiency of Fibrous Filters

Approaches for Predicting Collection Efficiency of Fibrous Filters Volume 5, Issue, Summer006 Approaches for Preicting Collection Efficiency of Fibrous Filters Q. Wang, B. Maze, H. Vahei Tafreshi, an B. Poureyhimi Nonwovens Cooperative esearch Center, North Carolina State

More information

ECE 422 Power System Operations & Planning 7 Transient Stability

ECE 422 Power System Operations & Planning 7 Transient Stability ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of

More information

ON THE MODELING AND DESIGN OF ZERO-NET MASS FLUX ACTUATORS

ON THE MODELING AND DESIGN OF ZERO-NET MASS FLUX ACTUATORS ON THE MODELING AND DESIGN OF ZERO-NET MASS FLUX ACTUATORS By QUENTIN GALLAS A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR

More information

Transactions on Engineering Sciences vol 5, 1994 WIT Press, ISSN

Transactions on Engineering Sciences vol 5, 1994 WIT Press,   ISSN Generalize heat conuction transfer functions in two imensional linear systems of assigne geometry F. Marcotullio & A. Ponticiello Dipartimento i Energetica- Universit i L 'Aquila Monteluco i Roio - 67100

More information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic 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 information

This is a repository copy of PIV studies of coherent structures generated at the end of a stack of parallel plates in a standing wave acoustic field.

This is a repository copy of PIV studies of coherent structures generated at the end of a stack of parallel plates in a standing wave acoustic field. This is a repository copy of PIV stuies of coherent structures generate at the en of a stack of parallel plates in a staning wave acoustic fiel. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/79741/

More information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical 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 information

Vortex Shedding on Combined Bodies at Incidence to a Uniform Air Stream T. Yavuz x, Y. E. Akansu xx, M. Sarıo lu xxx, and M.

Vortex Shedding on Combined Bodies at Incidence to a Uniform Air Stream T. Yavuz x, Y. E. Akansu xx, M. Sarıo lu xxx, and M. Vorte Sheing on Combine Boies at Incience to a Uniform Air Stream T. Yavuz, Y. E. Akansu, M. Sarıo lu, an M. Özmert. Ba kent Universit, : Nige Universit,, : Karaeniz Technical Universit,Turke Abstract

More information

Sources and Sinks of Available Potential Energy in a Moist Atmosphere. Olivier Pauluis 1. Courant Institute of Mathematical Sciences

Sources 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 information

Chapter 14: Vibration Generations Mechanisms: Self Excited Vibration

Chapter 14: Vibration Generations Mechanisms: Self Excited Vibration hapter 14: ibration Generations echanisms: Self Ecite ibration Introuction: Self-ecite systems begin to vibrate of their own accor spontaneously, the amplitue increasing until some nonlinear effect limits

More information

ADIT 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 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 information

EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD

EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD 137 THERMAL SCIENCE, Year 1, Vol. 16, No. 5, pp. 137-1376 EVOLUTION OF PARTICLE SIZE DISTRIBUTION IN AIR IN THE RAINFALL PROCESS VIA THE MOMENT METHOD by Fu-Jun GAN a an Jian-Zhong LIN a,b * a Department

More information

Sensors & Transducers 2015 by IFSA Publishing, S. L.

Sensors & 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 information

Formation Criterion for Synthetic Jets

Formation Criterion for Synthetic Jets AIAA JOURNAL Vol. 43, No. 1, October 25 Formation Criterion for Synthetic Jets Ryan Holman and Yogen Utturkar University of Florida, Gainesville, Florida 32611-625 Rajat Mittal George Washington University,

More information

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Dissipative numerical methods for the Hunter-Saxton equation

Dissipative numerical methods for the Hunter-Saxton equation Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a

More information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear Dielectric Response of Periodic Composite Materials onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics

More information

Convective heat transfer

Convective heat transfer CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,

More information

Design and Aerodynamic Characterization of a Synthetic Jet for Boundary Layer Control

Design and Aerodynamic Characterization of a Synthetic Jet for Boundary Layer Control Design and Aerodynamic Characterization of a Synthetic Jet for Boundary Layer Control FRANCESCA SATTA, DANIELE SIMONI, MARINA UBALDI, PIETRO ZUNINO Department of Fluid Machines, Energy Systems, and Transportation

More information

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall.

Fluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall. Flui Mechanics EBS 189a. Winter quarter, 4 units, CRN 52984. Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 eihmeyer Hall. Course Description: xioms of flui mechanics, flui statics,

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - 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 information

Optimum design of tuned mass damper systems for seismic structures

Optimum design of tuned mass damper systems for seismic structures Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty

More information

An inductance lookup table application for analysis of reluctance stepper motor model

An inductance lookup table application for analysis of reluctance stepper motor model ARCHIVES OF ELECTRICAL ENGINEERING VOL. 60(), pp. 5- (0) DOI 0.478/ v07-0-000-y An inuctance lookup table application for analysis of reluctance stepper motor moel JAKUB BERNAT, JAKUB KOŁOTA, SŁAWOMIR

More information

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration Prouction of pions, kaons an protons in Xe Xe collisions at s =. ev arxiv:09.0v [hep-ex] Sep 0, for the ALICE Collaboration Università i Bologna an INFN (Bologna) E-mail: simone.ragoni@cern.ch In late

More information

Thermal Modulation of Rayleigh-Benard Convection

Thermal Modulation of Rayleigh-Benard Convection Thermal Moulation of Rayleigh-Benar Convection B. S. Bhaauria Department of Mathematics an Statistics, Jai Narain Vyas University, Johpur, Inia-3400 Reprint requests to Dr. B. S.; E-mail: bsbhaauria@reiffmail.com

More information

Fluid Flow Around and Heat Transfer from Elliptical Cylinders: Analytical Approach

Fluid Flow Around and Heat Transfer from Elliptical Cylinders: Analytical Approach JOURNAL OF THERMOPHYSICS AND HEAT TRANSFER Vol. 19, No., April June 5 Flui Flow Aroun an Heat Transfer from Elliptical Cyliners: Analytical Approach W. A. Khan, J. R. Culham, an M. M. Yovanovich University

More information

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method 1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse

More information

Heat and mass fluxes across a density interfaces submitted to a gridgenerated

Heat and mass fluxes across a density interfaces submitted to a gridgenerated Heat an mass fluxes across a ensity interfaces submitte to a grigenerate turbulence Pascal DUPONT 1&2 et Hassan PEERHOSSAINI 2 1 LGCGM, EA3913, INSA e Rennes, Campus Beaulieu, 35043 Rennes ceex 2 Thermocinétique,

More information