Surrogate-Based Modeling of Cryogenic Turbulent Cavitating Flows
|
|
- August McKenzie
- 6 years ago
- Views:
Transcription
1 Proceedings of the 7 th Internationa Symposium on Caitation CAV009 Paper No. 77 August 17-, 009, Ann Arbor, Michigan, SA Surrogate-Based Modeing of Cryogenic Turbuent Caitating Fows Abstract The cryogenic caitation has critica impications on the performance and safety of iquid rocket engines. In this study, a systematic inestigation based on the surrogate modeing techniques is conducted to assess and improe the performance of a transport-based cryogenic caitation mode. Based on the surrogate mode, goba sensitiity anaysis is be conducted to assess the roe of mode parameters reguating the condensation and eaporation rates, and uncertainties in materia properties, specificay, the apor density and atent heat. The surrogate modes considered incude the response surface, radia basis neura network, Kriging, and a weighted aerage composite mode combining a surrogates. It is reeaed that the apor density and the mode parameter controing the eaporation rate are more critica than atent heat and the mode parameter controing the condensation rate. Based on the recommended mode parameter aues, better prediction of the cryogenic turbuent caitation can be attained. 1. Introduction For the fue deiery of iquid rockets, caitation occurs as the oca pressure is ower than apor pressure, which can seriousy compromise the engine performance and structura integrity [1,]. At cryogenic conditions, the iquid/apor density ratio is typicay ower than that under non-cryogenic conditions. Thus, to sustain a comparabe caity size, heat transfer caused by atent heat becomes more important due to higher mass transfer rates. Consequenty, the iquid surrounding the apor caity has more significant eaporatie cooing in the cryogenic iquids. In addition, other transport properties such as therma conductiity and thermodynamic properties such as apor pressure of cryogenic fuid are ery sensitie to temperature ariations. Substance Chien-Chou Tseng Department of Mechanica Engineering niersity of Michigan, Ann Arbor, MI, SA Heat capacity (J/Kg K) Density ratio Therma conductiity (W/mK) Latent heat (KJ/Kg) H O(98k) N (83k) H (0k) Tabe 1. Variation of physica properties for water (98k), iquid nitrogen (83k), and iquid hydrogen (k) on saturation cures. [5] Wei Shyy Department of Aerospace Engineering niersity of Michigan, Ann Arbor, MI, SA (a)vapor pressure s. temperature (b)liquid density s. temperature (c)density ratio (iq./apor) s. temperature (Soid ines represnt for water and use bottom and eft as x-axis and y-axis in each figure; dash ines represnt for iquid nitrogen and use top and right as x- axis and y-axis in each figure) Figure 1. Variation of physica properties for iquid nitrogen and water aong saturation cure. [5] Therefore, for cryogenic fuids, the energy transport equation needs to be incuded in the caitation mode [3,4]. Representatie aues of these quantities and the pressuretemperature saturation cures are summarized in Tabe 1 and Figure 1 [5]. The dynamic simiarity for isotherma case, such as water is goerned by the caitation number [1,,3]: σ = ( P P )0.5ρ (1) Where the caitation number σ is based on a constant apor pressure P at inet temperature T. P, ρ, and are reference pressure, iquid phase density and free stream eocity respectiey. For cryogenic caitation, the actua oca caitation number σ needs to be corrected according to the oca temperature T: σ = ( P P ( T )) / 0.5ρ () By the foowing first order approximation: 1 dp ρ ( σ σ ) = ( T T ) dt (3) dp T σ = σ ; T = T T < 0 dt 0.5ρ Equation (3) ceary shows the temperature dependency of caitation, and the oca temperature drop in cryogenic caitation wi produce a noticeabe rise for the oca caitation 1
2 number σ and hence suppress the caitation intensity [3,4]. The detai impact for the therma-sensibe materia properties to caitation mode wi be introduced ater. The numerica modeing of caitation argey foow two main categories: interface tracking methods with indiidua phases separatey treated [6,7], and homogeneous fow modes based on a singe-fuid framework with fuid properties estimated based on the iquid-apor mixture ratios [3,4,6,7,8,9,10,11,1,13,14,15]. Differences between the arious modes in the second category mosty come from the reation that defines the density fied. For oeriew of the arious modeing approaches, see, e.g., [3,4,8]. A homogeneous fow mode utiizing the framework of the transport-based equation (TEM) is adopted in the present study. In this method, the information of the apor oume/mass fraction distribution is obtained in a modeed transport equation based on the mass transfer between apor and iquid phases. This approach is we documented, see, e.g., [3,4,8]. To study the impact of the cryogenic mode parameters and fuid properties on the predictions, we conduct a mode improement exercise utiizing a goba sensitiity anaysis. We use the method deeoped by Sobo [16]. This method aows decomposition of a suitabe measure of prediction into the components of indiidua ariabes from which we can easiy cacuate the impact of each ariabe. To faciitate the framework which heps us probe the goba sensitiity of the caitation mode and fuid uncertainties, we wi first construct suitabe surrogate modes [17]. The practica utiity of surrogate modeing for design, optimization and sensitiity anaysis is we estabished [18,19]. There are many surrogate modes, but the mode that represents a particuar function the best is not known in adance. Thenthe predictions using different surrogate modes hae a certain amount of uncertainty. Goe et a [4,19] suggested that combinations of mutipe surrogate modes can be beneficia to quantify and to reduce uncertainties in predictions. They proposed a PRESS-based weighted aerage, namey PWS, (PRESS is the predicted residua sum of square) to reduce the mode uncertainties. Since the fideity of surrogate modes is critica in determining the success of the sensitiity anaysis, we wi use different surrogate modes to hep ascertain the performance measures. In this study, we use four surrogate modes, poynomia response surface approximation (PRS, [0]), Kriging (KRG, [1]), radia basis neura network (RBNN, []) and PRESS-based weighted aerage (PWS) surrogate mode constructed by using the preious three surrogates [4,17,19]. These surrogate modes are used to caibrate the mode parameters of the present transport based caitation mode for cryogenic caitation. The surrogate-based goba sensitiity anaysis can hep us inestigate the uncertainties from fuid properties and then identify the optima mode parameters to improe the prediction performance and robustness of cryogenic caitation modes. For turbuence cosure, the ensembe-aeraged modeing with a two-equation cosure [3] aong with a fiter-based mode (FBM) [4] is utiized. The approach reduces the infuence of the turbuent eddy iscosity based on the oca numerica resoution, essentiay bending direct numerica simuation (DNS) and conentiona turbuence mode in a singe framework. Specificay, the ee of the turbuent iscosity is corrected by comparing the turbuence ength scae computed from the turbuence cosure and the fiter size based on the oca mesh spacing. This approach can be aso categorized as Very Large Eddy Simuation (VLES) [5,6]. As documented in a preious study [7], the inet conditions of the turbuence mode can criticay affect the outcome of caitation structure. The fiter-based mode can hep significanty reduce the uncertainty in this regard. In this study, we re-examine the surrogate-based cryogenic modeing efforts preiousy reported by Goe et a. [4] to assess the predictie capabiity of turbuent caitating fows.. Goerning Equations and Numerica Approaches The set of goerning equations for cryogenic caitation under the homogeneous-fuid modeing consists of the conseratie form of the Fare-aeraged Naier-Stokes equations, the enthapy-based energy equation (for cryogenic caitation), the k-ε two-equation turbuence cosure, and a transport equation for the iquid oume fraction [3,4,8]. The continuity, momentum, enthapy, and caitation mode equations are gien beow. A computations presented beow are based on the steady-state equations. ( ρ u ) m = 0 x ( ρ uu ) p u u u m i i k = + [( µ + µ )( + δ )] L T i x x x x x 3 x i i k µ µ h m T [ ρ u ( h+ f L)] = [( + ) ] m x x Pr Pr x m T α + ( u ) = m& + m& (7) x where ρ m is the mixture density, u denotes the components of eocity, p is the pressure, µ and µ t are the mixture aminar and turbuent iscosities, respectiey, h is the sensibe enthapy, f is the apor mass fraction, L is the atent heat of aporization, Pr is the Prandt number, α is the iuid oume fraction, and m + and m are the source and sink terms for the caitation mode. The subscript t denotes turbuent properties, represents the iquid state, represents the apor state, and m denotes the mixture properties. The mixture property φ m and the apor mass fraction are, respectiey, expressed as φ = φα + φ ( 1 α ) (8) f m ρ ( 1 α ) = (9) ρ m The temperature can be interpoated based on enthapy in the data base [5]. We negect the effects of kinetic energy and iscous dissipation terms in Equation (4) (O(Re -0.5 ), Re is around 10 6 ) because the temperature fied is mainy contributed by the eaporatie cooing in cryogenic caitation. (4) (5) (6)
3 .1.Transport-Based Caitation Mode The source term m + and sink term m in Equation (7) represent for condensation and eaporation rates. They hae been deried from arious aspects, incuding dimensiona argument with empirica support [3,4,8,10,11,1], force baance based on the interfacia dynamics [3,8,9], and estimate of the bubbe growth rate through the Rayiegh-Pesset equation [13,14,15]. Numericay, [3,4,8,9,11,1] utiized pressure-based methods, and [10,1,16,13,14,15] empoyed the density-based methods. The iquid-apor eaporation and condensation rates for the present transport-based caitation mode [3,4,8,10,11,1] are respectiey shown as foowing: C α ρ min(0, p p ) dest m& = t ρ ρ (0.5 ) C (1 α ) max(0, p p ) + prod m& = (10) t (0.5 ρ ) where and are the empirica constants, is the reference eocity scae, and t is the reference time scae, which is the characteristic ength scae D diided by the reference eocity scae (t =D/ ). For non-cryogenic fuids ike water, the constants are specified =1 and =80 [3,8]. As for iquid nitrogen, the constants are chosen as =0.68 and =54.4, and for iquid hydrogen, =0.8 and =54.4 are suggested [3]. Further modifications are conducted to get a better prediction with =0.639 and =54.4 for iquid nitrogen, and =0.767 and =54.4 for iquid hydrogen [4]. For cryogenic caitation simuations, the temperature dependent properties are updated from a comprehensie data base [5] throughout the course of computations for eery iteration... Thermodynamics Effects The impact of therma effects in cryogenic caitation due to phase change on temperature prediction has been aready shown in Figure 1. These thermo-sensibe materia properties wi affect the energy equation in Equation (6) and caitation sink/source terms in Equation (10). First, the atent heat L in Equation (6) appears as a noninear source term and represents the atent heat transfer rate during the phase change. The spatia ariation of the thermodynamic properties together with the eaporatie cooing effect is embedded into this equation and causes a couping with the set of goerning equations [3,4]. As for the caitation sink/source terms in Equation (10), we can assess the impacts due to the thermo-sensibe materia properties by using Tayor s series and negect the higher order terms. We first consider the sink term m as the pressure is smaer than the apor pressure based on the oca temperature [4,7], or in other words, pressure coefficient C p is smaer than σ. Furthermore, the minus sign here in Equation (11) is for conenience to show a arger eaporation strength wi correspond to a arger positie aue of m : C α ρ min(0, p p ( T)) dest m& = ( ) = βr( T) min(0, C + σ) p t ρ 0.5ρ dr dp T = β( R( T ) + T+...)( C + σ +...) p dt dt 0.5ρ T dr dp T = β( R( T )( C + σ ) + ( C + σ ) T R( T ) +...) p p dt dt 0.5ρ T T T (11) where β is α /t and R is the temperature-dependent iquid/apor density ratio, Simiary for source term m + as the pressure is arger than the apor pressure (C p is arger than σ): C (1 α ) prod max(0, p p ( T)) + m& = = γ max(0, C + σ ) p t 0.5ρ dp T = γ (( C + σ ) +...) p dt 0.5ρ T (1) Where γ is (1-α )/t. Pease note that T=T-T <0 in both Equation (11) and (1), which is aso defined in Equation (3) as eaporatie cooing occurs. It can be concuded that the competing infuence of the therma effects in the caitation mode comes from two ways from Equation (11) as eaporatie cooing occurs: (1) therma rate of change of iquid/apor density ratio dr/dt which is negatie in Figure 1(c) together with C p +σ <0 and T<0 as eaporation occurs, wi tend to enhance the strength of m and () therma rate of change of apor pressure dp /dt which is positie in Figure 1 to suppress m. It is aso obious that the impacts of therma effects wi change significanty for different working temperature and pressure due to the non-inear ariation of materia properties from energy equation in Equation (6) and caitation sink/source terms in Equation (10)..3. Turbuence Mode The k-ε two-equation turbuence mode with a wa function treatment is presented as foows [3]: ( ρ u k) µ k m T = P ρ ε + [( µ + ) ] t m L x x σ x k ( ρ uε ) m ε ε µ ε T = C P C ρ + [( µ + ) ] ε t ε m L 1 x k k x σ x ε (13) (14) where the production term of turbuent kinetic energy (P t ) and the Reynods stress tensor are defined as: ui P = τ ; τ = ρ u u t i i m i x ρ kδ u u m i i ρ u u = µ ( + ) m i T 3 x x i (15) with C ε1 =1.44, C ε =1.9, σ ε =1.3,σ k =1.0. The turbuent eddy iscosity is defined as: ρ C k m µ µ =, C = 0.09 (16) T µ ε 3
4 As mentioned aboe, a fiter-based mode (FBM) [4] is aso adopted. This mode imits the infuence of the eddy iscosity based on the oca numerica resoution, essentiay forming a combined direct numerica simuation and RANS mode. Specificay, the ee of the turbuent iscosity is corrected by comparing the turbuence ength scae and the fiter size, which is based on the oca meshing spacing: k ε µ = 0.09ρ min(1, ) T m ε 3/ k (17) By imposing the fiter, the turbuence ength scae wi not be resoed if it is smaer than the fiter size. The fiter size is chosen to be comparabe to the maximum grid size: = max(, ) (18) present grid Thus if the grid resoution is significanty smaer than the turbuence ength scae in the entire fow fied, the soution wi approach that of a direct numerica simuation; for inadequatey resoed computations, the RANS mode is recoered. Simiar concepts can be found in studies of VLES [5,6]..4.Numerica method Detaied numerica procedures for the caitation mode and associated fuid dynamics equations adopted here utiize a modified pressure-based approach for arge density ump as we as therma effects, as reported in [8,9] The controed ariation scheme (CVS) [8] is appied to discretize the conection scheme, and centra difference is used for both pressure and diffusion terms. The CVS scheme can preent the osciations under sharp gradients caused by the phase change whie presering second-order accuracy esewhere. As for the boundary conditions, iquid oume fraction, eocity, temperature and turbuent quantities are specified at the inet. For the outet, pressure and other fow ariabes are extrapoated. On the was, pressure, iquid oume fraction, and turbuent quantities are extrapoated aong with no-sip and adiabatic conditions. Additionay, the pressure at the reference point (P ) in the upstream is aso fixed to define the caitation number σ [8,9]. Based on the eddy-to-aminar iscosity ratio at the inet, the inet turbuent quantities can be gien as foowing: 3 C ( ), µ k k = I ε = ν L ( µ T / µ L ) inet (19) where I is turbuence intensity (% for entire study), and eddyto-aminar iscosity ratio, is equa to1000 in a the cases during the current study. 3. Resuts and Discussion 3.1. Test Geometry We simuate fow oer a D quarter hydrofoi in Figure with the experimenta measurements by Hord [9]. Since we hae iustrated in Figure 1 and used Equation (11) and (1) to highight the impacts of therma effects wi change significanty for different temperatures, we seect both Case 90C and 96B to caibrate the roe of mode parameters and uncertainty of fuid properties. Figure and Tabe summarize the geometries, corresponding boundary conditions, and the fow conditions of the test cases seected to aid the mode aidation. Substance Case σ Re T Liquid nitrogen 90C K Liquid nitrogen 96B K Tabe. Summary of simuation setups and fuid properties Figure. Schematic of the geometries and the boundary conditions of the cases considered. 3.. Pressure and Temperature Predictions Figure 3 compares the predicted and experimentay measured pressure and temperature profies [9] on the hydrofoi surface with =0.68 and =54.4 suggested in [3]. Oera, the caitation and turbuence modes with fiter (FBM) can consistenty capture the main features of both pressure and therma profies. The temperature drop inside the caity in Figure 3(b) and (d) aso ceary demonstrates the eaporatie cooing resuting from cryogenic caitation. (a)pressure, 90C (c)pressure, 96B (b)temperature, 90C (d)temperature, 96B Figure 3. Pressure and temperature of cryogenic caitation aong surface. In Figure 4, we compare the present cryogenic mode soution with the isotherma soution for Case 90C, obtained by using the identica mode with =0.68 and =54.4 except that the energy equation is not inoked. Ceary, the therma fied does affect the caity structures. The caity size in Figure 4(b) is reduced due to the therma effect because the temperature drop inside the caity in Figure 4(a) decreases the oca apor pressure and hence increases the oca caitation number (pease refer to Equation (3)), resuting in a weaker 4
5 caitation intensity and higher oera iquid oume fraction in the caity (as shown in Figure 4(c)). Besides, the pressure inside the caity is steeper under the cryogenic condition than that under the isotherma condition in Figure 4(d) due to the ariation in therma rate of change of apor pressure dp /dt. The purpose of Figure 4 is to iustrate the therma effects coud significanty change the pressure fieds and caity structures under cryogenic conditions. (a) α with energy equation (b) α without energy equation (c) α aong surface (d)pressure aong surface Figure 4 Comparisons for Case 90C between resuts with/without energy equation 3.3 Surrogates-based Goba Sensitiity Assessment Since minor changes in fow enironment can ead to substantia changes in the predictions in cryogenic enironment, it is imperatie to appraise the roe of mode parameters and uncertainties in materia properties on the predictions. In this section, we characterize the parameters using surrogate-based goba sensitiity anaysis (GSA) and then caibrate the cryogenic caitation mode parameters. We use poynomia response surface approximation (PRS, [0]), Kriging (KRG, [1]), radia basis neura network (RBNN, []), and a weighted aerage surrogate mode (PWS) for approximation of response. [4,19]. We use ariance-based, non-parametric GSA method, proposed [16] to eauate the sensitiity of predictions with respect to mode parameters and materia properties. In this method, the obectie function is decomposed into additie functions of ariabes and interactions among ariabes. This aows the tota ariance (V) in the obectie function to be expressed as a combination of the main effect of each ariabe (V i ) and its interactions with other ariabes (V i ). The sensitiity of the obectie function with respect to any ariabe is measured by computing its sensitiity indices. The sensitiity indices of main effect (S i ) and tota effect ( S are gien as foows: tota i ) of a ariabe V i tota ( ), V i+ S V iz i = S V i = (0) V We choose,, ρ,, and L as design ariabes, whie hoding the Re and σ constant for the gien cases. The performance of predictions for the cryogenic caitation modes are represented by root mean square (RMS) aues of the differences between computed and experimenta aues aong hydrofoi surface for temperature ( ) and pressure ( ) as our obecties. The mode parameters, and, ary from to 0.68 and 46. to 54.4 respectiey. The materia properties ρ and L are perturbed within ±10% of the aue they assume from the NIST database [5]. These two empirica constants and directy contro the eaporation and condensation rate ia the caitation mode. Besides, ρ, as a fuid property which dominates the eaporatie cooing, aso appears directy in caitation sink term, and L wi determine the energy absorb or reease during the phase change. Therefore, these four mode parameters and fuid properties are seected as our design ariabes. To faciitate the deeopment of surrogate modes, 70 training points are seected using combined face-centered cubic composite design (FCCD, 5 points) and Latin hypercube samping (LHS, 45 points). Fie additiona test points which are not incuded in the 70 training points are used to aidate the surrogate modes for both case 90C and 96B. We eauate and for each data point using CFD simuations and construct PRS, Kriging, RBNN, and PWS modes of both obecties in normaized ariabe space. A ariabes and obecties are normaized such that 0 corresponds to the minimum aue and 1 corresponds to the maximum aue. Normaized ariabes and obecties are denoted by a superscript. We use second order poynomias for PRS and a spread coefficient=0.4~0.7 for RBNN. Reeant detais of the quaity of fit of surrogate modes are documented in [4]. The smaest root mean square of PRESS (PRESS is the predicted residua sum of square) in Appendix Tabe A1 indicates that the KRG mode has the best performance, whie RBNN has the worst oera performance. The contribution of different surrogate modes to the PWS mode is gien by the weights in Appendix Tabe A. Since the performance of KRG is the best, its weight is aso the greatest. There are additiona fie test points to aidate the surrogates. Appendix Tabe A3 shows the ocations of these points in the normaized design space for both Case 90C and 96B. The simuation resuts are presented to compare with the prediction of surrogates in Appendix Tabe A4 and A5 for Case 90C, and Appendix Tabe A6 and A7 for Case 96B. Due to the best performance in error estimate in Appendix Tabe A1 and tests for additiona fie sampes from Appendix Tabe A4 to A7, we wi use KRG to demonstrate the goba sensitiity anaysis. Appendix Tabe A8 and A9 show that the interactions among parameters are ery strong for pressure since the difference between main and tota sensitiity indices are obious for Case 90C and 96B. As for temperature, Case 96B has stronger interactions than those of Case 90C. From Figure 5 and 6 which eauate the weights of each ariabes ia goba sensitiity anaysis as pie-charts, and ρ, are ery important for, and whie and L don t hae noticeabe contribution in both cases. Besides, the weights of are ery simiar for Case 90C and 96B. As for, the 5
6 importance of L ceary increases, and een becomes the same important as and ρ in Case 96B, whie in Case 90C, the weight is ust around 6%. This indicates that the therma effects and importance of the thermodynamics property wi infuence the therma fied more significanty as temperature increases (Case 96B has higher inet temperature than Case 90C). ρ Figure 5. Pie-chart of goba sensitiity anaysis for Case 90C (KRG (Kriging), tota indices) ρ Figure 6. Pie-chart of goba sensitiity anaysis for Case 96B (KRG (Kriging), tota indices) is not important within this design space from the piechart, and this impies that the sensitiity of condensation term is not significant compared with eaporation term or. This shoud be because caitation initiates from eaporation, and this behaior aready decides how strong the caitation dynamics and how ow the iquid oume fraction wi be inside the caity. Besides, the condensation strength wi aso depend on the eaporation strength. Just imagine if we gie a ery ow strength on eaporation term, the condensation dynamics wi sti be weak een is arge. This is because the apor inside the caity wi not be sufficient for source term. Thus the weight of is much more important than in this design space. As we know, the caitation terms can be expressed as foowing in Equation (10): C α ρ min(0, p p ) dest m& = t ρ ρ L L (0.5 ) C (1 α ) max(0, p p ) + prod m& = t (0.5 ρ ) Therefore we can group /t /ρ and /t together to show the ariation in strength of caitation sink and source term under the combinations of design ariabes, and then ρ L L ρ normaize these aues in Figure 7 and 8. This is the direct impact which indicates the direct appearance of design ariabes in the caitation mode. A the normaized aues here are from the preious simuation resuts of the 70 training points. Figure 7 and 8 show the distributions of P diff s. ( /t /ρ ) are quite the same for both Case 90C and 96B. This is because both cases hae consistent pie-charts and weights in the pressure prediction in Figure 5 and 6. As, ( /t /ρ ) is sma in Figure 7 and 8, the sink term is not strong enough so that the caity size is too sma, and hence P diff wi be arge. When ( /t /ρ ) goes up to certain moderate aue, the corresponding eaporation term gies a more suitabe caity size with a smaer aue of P diff. For een arge aues of ( /t /ρ ), the caity sizes wi be too arge, and wi increase again. This ceary indicates that there exists a suitabe range for sink term or to obtain good pressure predictions. For the same aue of ( /t /ρ ), different P diff and T diff can sometimes be obtained. This is because different fuid properties, namey ρ or L, infuence the fow fieds. We haen t reay incuded this factor in Figure 7 and 8, and refer it as indirect impact. As for ( /t /ρ ) s. T diff, we sti can roughy see simiar trends as shown in ( /t /ρ ) s. P diff, but the distributions are not so consistent for both Case 90C in Figure 7 and 96B in Figure 8. Again, this is due to different pie-charts of for both cases as shown in Figure 5 and 6. This aso impies different impacts of therma effects at different temperatures which we hae aready stated before. For ρ, besides direct appearance in the sink term, it is aso part of fuid property. As a fuid property, it can change the aue of mixture density to exhibit different eaporatie cooing inside the caity, and hence affects the caity size and caitation dynamics. Therefore, it owns both direct and indirect impacts, and we can see its weight of from goba sensitiity anaysis is een arger than in Figure 5 and 6. As for ( /t ), there is reay no trend for P diff and. This is because from the pie-chart in Figure 5 and 6, we can see the importance of is insignificant, and a these distributions in Figure 7 and 8 woud be mainy due to the contributions of sink term. Therefore the random distributions for ( /t ) s. P diff and in Figure 7 and 8 hep us aidate the insignificant weight of in Figure 5 and 6 from another iewpoint. 3.4 Optimization of Mode Parameters In the preious section, we obsered that one of the mode parameters significanty infuences the performance of the present cryogenic caitation mode. Therefore we fix another mode parameter at 54.4 which is not reay infuentia on predictions, and assume the temperature-dependent materia properties ρ and L obtained from the NIST database [5] to be correct. Then we aow the mode parameter to ary between and 0.68, and try to find the optimization aues which gie and as ow as possibe. There are 11 training points within the design space with equa spacing. From the simuations shown in Figure 9, we can see for Case 90C, the trends of these two obecties are amost opposite, and for Case 96B, these two obecties are of the same trends. 6
7 . ( /t /ρ ) ( /t /ρ ) (a) 90C (b) 96B Figure 9. Location of points ( ) and corresponding responses ( is shown on eft y-axis, and is shown on right y-axis) used for caibration of the cryogenic caitation mode. Figure 7. Iustration of the direct impact of Case 90C. ( /t ) ( /t ) The function space and the ariabe space iustration of Pareto optima front (POF) obtained through different surrogate modes are shown in Figure 10 and 11. We obsere that different POFs obtained by using mutipe surrogate modes are cose to one another in both function and ariabe spaces (Besides PRS for in Case 90C). We note that the pressure fuctuations pay a more important roe in determining the caitation dynamics and the oadings on fuid machinery. Consequenty, more accurate pressure prediction is our primary obectie. We seect a soution on the POF for aidation, such that notabe reduction in can be reaized without incurring significant deterioration of.. Therefore the optima wi be 0.65 from Figure 10 and 11. It aso indicates that different therma effects wi not significanty affect the choice of optima. ( /t /ρ ) ( /t /ρ ) (a) Variabe space (b) Function space Figure 10. Pareto optima front and corresponding optima points for Case 90C. ( /t ) ( /t ) Figure 8. Iustration of the direct impact of Case 96B. To represent the responses and using surrogate modes, we sampe data using CFD simuations from Figure 9. and estimated by surrogate modes shown in Figure 10(a) and 11(a) ceary exhibit the same trend as simuations in Figure 9. As before, we construct PRS, KRG, RBNN, and PWS modes. The error estimates and the weights associated with different surrogates in PWS mode are summarized in Tabe A10 and A11. Besides PRS for in Case 90C, a the surrogates yied a good accuracy. (a) Variabe space (b) Function space Figure 11. Pareto optima front and corresponding optima points for Case 90C. (KRG,Kriging; PRS, Poynomia Response Surface; RBNN, Radio Basis Neura Network; PWS, PRESS-Based Weighted Aerage Surrogate.) 7
8 3.5 Optimization for Liquid hydrogen We hae repeated simiar mode caibration process for iquid nitrogen by considering Case 49D (σ =1.57, T =0.7K, Re = 10 7 ) and 55C (σ =1.49, T =.K, Re = ), the optima aue of wi be Improement by Fiter-based Turbuence Mode Compared with Goe. et a [4], our current resuts exhibit improement in both pressure and temperature prediction by utiizing fiter-based mode (FBM, [4]). This mode heps us reduce the uncertainties from inet turbuent quantities. The improement is summarized in Tabe 3. Pease note that our is 0.65 for iquid nitrogen and 0.78 for iquid hydrogen, whie Goe.et a use and respectiey. Case Fuid Current resuts Goe et a [4] Current resuts Goe et a [4] 90C Nitrogen B Nitrogen D Hydrogen C Hydrogen Tabe 3. Performance compasions between current study and Goe. et a [4]. We hae aso used =0.65 to other iquid nitrogen cases isted in the experiment [9], and the simuation resuts are aso ery consistent to experimenta data in Tabe 4. This exercise heps us to aidate our mode parameters are aso suitabe for other cases. Case σ T (K) 83B C A Tabe 4. =0.65 for other iquid nitrogen cases. 4. Concusions In this study, we choose,, ρ,, and L as design ariabes and use different surrogate modes, PRS, KRG, RBNN, and PWS, to construct the response of pressure prediction and temperature prediction for cryogenic caitation. As documented in our preious study [7], the inet conditions of the turbuence mode can criticay affect the outcome of caitation structure. The fiter-based mode can hep significanty reduce the uncertainty in this regard. We hae re-examined the surrogate-based cryogenic modeing efforts preiousy reported by Goe et a. [4] based on a more refined turbuence approach to better probe the combined turbuent caitating fow predictie capabiities. It is found that the performance of the current caitation mode is affected more by mode parameter associated with the eaporation term and the fuid property ρ which contros the eaporatie cooing. The condensation term ( ) is not important at a within this design space, and atent heat L becomes significant ony in temperature prediction. We hae aso obsered that there is a range of eaporation term or which can gie better accuracy of pressure prediction. We recommend =0.65 and 0.78 for iquid nitrogen and hydrogen respectiey, and is 54.4 for both fuids. 5. Acknowedgments The present efforts hae been supported by NASA Consteation niersity Institute Program (CIP). 6. Appendix Acronyms in this appendix: PRS Poynomia Response Surface [0] KRG Kriging [1] RBNN Radio Basis Neura Network [] PWS PRESS-Based Weighted Aerage Surrogate [4,17,19] PRESS The predicted residua sum of square R ad Coefficient of mutipe determination [0] (better fit as it is coser to 1) Surrogate 90C 96B 90C 96B PRESS rms of PRS 6.38% 11.90% 9.11% 10.0% R ad of PRS PRESS rms of KRG.97%.93%.48% 6.6% PRESS rms of RBNN 13.91% 11.67% 13.31% 19.03% PRESS rms of PWS 3.97% 5.44% 5.50% 9.0% Tabe A1 Error estimate for different cases and surrogates (70 training points) Surrogate 90C 96B 90C 96B PRS 44.% 7.5% 31.8% 36.4% KRG 47.4% 44.5% 45.0% 39.6% RBNN 8.4% 8.3% 3.% 4.0% Tabe A Weights associated with different surrogate modes (70 training points) No No No No No Tabe A3. Test points inside the normaized design space PRS KRG RBNN PWS CFD -.3% -.7% -8.0% -3.6% % -.7% -5.4% -5.0% % 3.0% 1.0% 14.8% % -5.3% -1.6% 4.3% % 0.% -3.7% -.7% Tabe A4. Predictions error of for case 90C ρ L 8
9 PRS KRG RBNN PWS CFD -0.% -0.5% 1.% 0.0% % -1.7% -9.1% -.4% % 1.8% -0.7%.7% % -1.1% 1.1% 0.8% % 0.1% -0.% -1.4% Tabe A5. Predictions error of for case 90C PRS KRG RBNN PWS CFD -15.9% -5.1% 1.%.8% % -3.6% 46.3% -0.9% % -6.6% -6.0% 0.4% % -6.8% 8.1% -4.6% % -4.6% 8.8% -.7%.6871 Tabe A6. Predictions error of for case 96B PRS KRG RBNN PWS CFD 0.6% 0.3% -.% 0.% % 6.4% 1.5% -4.9% % 6.7% 0.% -5.5% % 0.9% -1.6% -0.4% % 1.0% -1.3% 0.1% Tabe A7. Predictions error of for case 96B (The coumn for CFD in Appendix Tabe A4 to A7 denotes simuation resuts before normaized.) 0.176/ E-5/ / / /0.53 1E-4/ / /0.06 Tabe A8. Goba sensitiity for 90C with KRG (main/tota effect) 0.035/0.61 E-4/ / / / / / /0.404 Tabe A9. Goba sensitiity for 96B with KRG (main/tota effect) ρ ρ Surrogate 90C 96B 90C 96B PRESS rms of PRS 16.01% 4.46% 3.8% 6.00% R ad of PRS PRESS rms of KRG 7.01% 4.87% 1.89% 1.91% PRESS rms of RBNN 4.00% 6.74% 1.56% 0.73% PRESS rms of PWS 7.45% 5.10% 3.0% 1.83% Tabe A10 Error estimate for different cases and surrogates (11 training points) Surrogate 90C 96B 90C 96B PRS 0.8% 35.8% 3.6% 13.5% KRG 36.9% 34.5% 37.1% 39.5% RBNN 4.3% 9.6% 39.3% 47.0% Tabe A11 Weights associated with different surrogate modes (11 training points) L L 7. Nomencature σ Free stream caitation number σ Caitation number based on the oca temperature C ε1,c ε Coefficients of k-ε turbuence mode σ ε,σ k Coefficients of k-ε turbuence mode C p Pressure coefficient D Characteristic ength scae f Vapor mass fraction h Enthapy I Turbuence intensity K Turbuent kinetic energy L Latent heat m +, m - Source and sink terms in the caitation mode L norm between experiment and predicted pressure Pr Prandt number Pt Production term of turbuent kinetic energy P Saturation apor pressure Re Reynods number S Sensitiity indices T Temperature L norm between experiment and predicted temperature t Reference time scae, t =D/ Reference eocity u Veocity x Space ariabe V tota ariance α Liquid oume fraction ρ Density Dynamic iscosity Eddy-to-aminar iscosity ratio at the inet Mixture property ε Turbuent dissipation rate Fiter size in fiter-based mode Subscript i Interaction component Component Liquid L Laminar m Mixture property T Turbuent Vapor ω Free stream quantities Superscript Normaize aue tota Tota effect 9
10 8. References [1]Knapp RT Daiy JW, Hammitt FG. 1970, Caitation. McGraw-Hi, New York. []Brennen CE. 1995, Caitation and Bubbe Dynamics, Oxford niersity Press, New York. [3]tturkar Y, Wu J, Wang G, Shyy W. 005, Recent Progress in Modeing of Cryogenic Caitation for Liquid Rocket Propusion, Progress in Aerospace Sciences, 41(7), [4]Goe T, Zhao J, Thakur S, Haftka RT, Shyy W, Zhao J. 008, Surrogate Mode Based Strategy for Cryogenic Caitation Mode Vaidation and Sensitiity Eauation, Int. J. Numer. Meth. Fuids, 58, [5]Lemmon EW. McLinden M.O, Huber M.L. REFPROP: 00, Reference Fuid Thermodynamic and Transport Properties, NIST Standard Database 3, ersion 7.0. [6]Chen Y, Hesiter SD. 1994, A Numerica Treatment for Attached Caitaition, J Fuids Eng, 116, [7]Deshpande M, Feng J, Merke CL. 1997, Numerica Modeing of the Thermodynamic Effects of Caitation, J. Fuids Eng, 119, [8]Senocak I, Shyy W. 004, Interfacia Dynamics-Based Modeing of Turbuent Caitating Fows, Part-1: Mode Deeopment and Steady-State Computations, Int. J. Numer. Meth. Fuids, 44, [9]Senocak I, Shyy W. 004, Interfacia Dynamics-Based Modeing of Turbuent Caitating Fows, Part-: Time- Dependent Computations. Int. J. Numer. Meth. Fuids, 44, [10]Merke CL, Feng J, Bueow PEO. 1998, Computationa Modeing of Sheet Caitation, Proc. 3 rd Internationa Symposium on Caitation, Grenobe, France. [11]Singha Ak, Li H, Athaae MM, Jiang Y. 001, Mathematica Basis and Vaidation of the Fu Caitation Mode, ASME Paper FEDSM [1] Kunz RF, Boger DA, stinebring DR, Chyczewski TS, Lindau JW, Gibeing HJ, Venkateswaran S, Goindan TR. 000, A preconditioned Naier stokes method for two phase fows with appication to caitation prediction, Comput. Fuids, 9, [13]Hosangadi A, Ahua V. 006, A Numerica of Caitation in Cryogenic Fuids Part II; New nsteady Mode for Dense Coud, Caitation. 6 th Internationa Symposium on Caitation, Wageningen, Netherands. [14]Aiabadi S, Tu SZ, Watts MD. 005, Simuation of Hydrodynamic Caitating Fows sing Stabiized Finite Eement Method, 43 th AIAA Aerospace Science Meeting and Exhibit, Reno, Neada. [15]Giorgi MGD, Ficarea A, Chiara F, Laforgia D. 005, Experimenta and Numerica Inestigations of Caitating Fows, 35 th AIAA Fuids Dynamics Conference and Exhibit, Toronto, Ontario Canada. [16]Sobo, IM, 1993, Sensitiity Anaysis for Noninear Mathematica Modes, Mathematica Modeing & Computationa Experiment 1(4), 407. [17]Queipo, NV, Haftka, RT, Shyy, W, Goe, T, Vaidyanathan, R, and Tucker, PK, 005, Surrogate Based Anaysis and Optimization, Progress in Aerospace Sciences, 41, 1-8. [18]Li, W, Padua, S. 004, Approximation Methods for Conceptua Design of Compex Systems, Proceedings of 11 th Internationa Conference on Approximation Theory, (eds. Chui C, Neaumtu M, Schumaker L). [19]Goe, T, Haftka, RT, Shyy, W, Queipo, NV, 006, Ensembe of Mutipe Surrogates, accepted for pubication in Structura and Mutidiscipinary Optimization. [0]Myers, RH, Montgomery, DC, 1995, Response Surface Methodoogy, John Wiey, Sons, Inc. [1]Matheron, G, 1963, Principes of Geostatistics, Economic Geoogy, 58, []Orr, MJL, 1996, Introduction to Radia Basis Function Networks, Center for Cognitie Science, Edinburg niersity, EH 9LW, Scotand, K. [1]Matheron, G, 1963, Principes of Geostatistics, Economic Geoogy, 58, []Orr, MJL, 1996, Introduction to Radia Basis Function Networks, Center for Cognitie Science, Edinburg niersity, EH 9LW, Scotand, K. [3]Launder BE, Spading DB, 1974, The Numerica Computation of Turbuent Fow, Comp. Meth. App. Mech. Eng, 3, [4]Wu JY, Shyy W, 004, Fiter Based nsteady RANS Computationa, Internationa Journa of Heat and Fuid Fow, 5, [5]Ruprecht A, Hemrich T, Buntic I. 005, Very Large Eddy Simuation for the Prediction of nsteady Votex Motion, Conference on Modeing Fuid Fow (CMFF 05), The 1 th Internationa Conference on Fuid Fow Technoogies. Budapest, Hungary. [6]Liu NS, Shih TS, 006, Turbuence Modeing for Very Large-Eddy Simuation, AIAA J., 44(4), [7] Tseng C and Shyy W., 009, Turbuence Modeing for Isotherma and Cryogenic Caitation, AIAA Paper No In 47th AIAA Aerospace Science Meeting, Orando, FL [8]Shyy W, Thakur SS, Ouyang H, Liu J, Bosch E., 007, Computationa Techniques for Compex Transport Phenomenon. Cambridge niersity Press, New York. [9]Hord J, 1973, Caitation in Liquid Cryogens II-Hydrofoi. NASA CR-156.,. 10
Phase Change Equation of State for FSI Applications
15 th Internationa LS-DYNA Users Conference FSI / ALE Phase Change Equation of State for FSI Appications Mhamed Soui, Ramzi Messahe Lie Uniersity France Cyri Regan, Camie Ruiuc Ingeiance Technoogies, Agence
More informationNumerical Study on Subcooled Pool Boiling
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vo., pp.15-19 (011) ARTICLE Numerica Study on Subcooed Poo Boiing Yasuo OSE * and Tomoaki KUNUGI Kyoto Uniersity, Yoshida, Sakyo, Kyoto, 606-8501, Japan This
More informationCavitation Simulation on Conventional and Highly-Skewed Propellers in the Behind-Hull Condition
Second Internationa Symposium on Marine Propusors smp 11, Hamburg, Germany, June 2011 Caitation Simuation on Conentiona and Highy-Skewed Propeers in the Behind-Hu Condition Keun Woo Shin 1, Pou Andersen
More informationChapter 11: Two-Phase Flow and Heat Transfer Flow Models Flow Models
11.3 Fow Modes Chapter 11: Two-Phase Fow 11.3 Fow Modes The homogeneous fow mode proides an easier approach to determining fow properties and behaiors, but it underestimates the pressure drop, particuary
More information), enthalpy transport (i.e., the heat content that moves with the molecules)
Steady-state conseration statements for a composite of ces and airspace In steady state, conseration of moecues requires that the tota fux into a representatie oume of mesophy is equa to the fux out storage
More informationANALYSIS OF FLOW INSIDE THE FOCUSING TUBE OF THE ABRASIVE WATER JET CUTTING HEAD
7 American WJTA Conference and Expo August 9-, 7 Houston, Texas Paper ANALYSIS OF FLOW INSIDE THE FOCUSING TUBE OF THE ABRASIVE WATER JET CUTTING HEAD Viém Mádr, Jana Viiamsoá, Libor M. Haáč VŠB Technica
More information7.3 Filmwise Condensation Regimes of Filmwise Condensation
Adanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howe 7.3 Fimwise Condensation 7.3.1 Regimes of Fimwise Condensation Re 30 Figure 7.14 Fow regimes of fim condensate on a ertica wa.
More informationCOMPARISON OF HEAT TRANSFER CHARACTERISTICS IN SURFACE COOLING WITH BOILING MICROJETS OF WATER, ETHANOL AND HFE7100
4th Micro and Nano Fows Conference UCL, London, UK, 7-0 September 04 COMPARISON OF HEAT TRANSFER CHARACTERISTICS IN SURFACE COOLING WITH BOILING MICROJETS OF WATER, ETHANOL AND HFE700 Da0riusz Mikieewicz,
More informationCAV2009 Paper No. 9. Unsteady Dynamics of Cloud Cavitating Flows around a Hydrofoil
Proceedings of the 7 th Internationa Symposium on aitation AV009 August 17-, 009, Ann Arbor, Michigan, USA AV009 Paper No. 9 Unsteady Dynamics of oud aitating Fows around a Hydrofoi Guoyu Wang Emai: wangguoyu@bit.edu.cn
More informationEffect of Oxygen Injection into Argon Induction Plasmas on Chemically Non-Equilibrium Conditions
Proceedings of 17th Internationa Symposium on Pasma Chemistry, Toronto, Canada, August 7-12, 25 Effect of Oxygen Injection into Argon Induction Pasmas on Chemicay Non-Equiibrium Conditions Nobuhiko Atsuchi
More informationWilson Ogoh and Dominic Groulx Department of Mechanical Engineering Dalhousie University, Halifax, NS, Canada
Presented at the COMSOL Conference 2010 Boston Wison Ogoh and Dominic Groux Department of Mechanica Engineering Dahousie University, Haifax, NS, Canada Sensibe Heat Storage: A heat storage system that
More informationSimulations of Droplets falling on a solid surface Using Phase-Field Method
APCOM & ISCM 11-14 th December, 013, Singapore Simuations of Dropets faing on a soid surface Using Phase-Fied Method T. Sakakiabara¹, *T.Takaki 1, and M.Kurata 1 Graduate Schoo of Science and Technoogy,
More informationEffects of mass transfer time relaxation parameters on condensation in a thermosyphon
Journa of Mechanica Science and Technoogy 29 (12) (2015) 5497~5505 www.springerink.com/content/1738-494x(print)/1976-3824(onine) DOI 10.1007/s12206-015-1151-5 Effects of mass transfer time reaxation parameters
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationResearch on liquid sloshing performance in vane type tank under microgravity
IOP Conference Series: Materias Science and Engineering PAPER OPEN ACCESS Research on iquid soshing performance in vane type tan under microgravity Reated content - Numerica simuation of fuid fow in the
More informationThe Growth of Vapor Bubbles in the. Volume of Superheated Drops, Dispersed. in High-Boiling Liquid
Appied Mathematica Sciences, Vo. 8, 2014, no. 151, 7519-7528 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.49183 The Growth of Vapor Bubbes in the Voume of Superheated Drops, Dispersed
More informationDYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE
3 th Word Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 38 DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE Bo JIN SUMMARY The dynamic responses
More informationA SECOND ORDER TURBULENCE MODEL BASED ON A REYNOLDS STRESS APPROACH FOR TWO-PHASE BOILING FLOW AND APPLICATION TO FUEL ASSEMBLY ANALYSIS
A SEOND ORDER TURBULENE MODEL BASED ON A REYNOLDS STRESS APPROAH FOR TWO-PHASE BOILING FLOW AND APPLIATION TO FUEL ASSEMBLY ANALYSIS S. Mimouni 1, F. Archambeau 1, M. Boucer 1, J. Laieie 1,. More stephane.mimouni@edf.fr,
More information2.4 Volume-Averaged Models
Adanced Heat and Mass Transer by Amir Faghri, Yuwen Zhang, and John R. Howe 2.4 Voume-Aeraged Modes 2.4.1 Oeriew o Aeraging Approaches The objecties o the arious aeraging methods are twood: (1) to deine
More informationIn-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017
In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative
More informationMODELING OF VOID FORMATION DURING RESIN TRANSFER MOLDING
MODELIG OF VOID FORMAIO DURIG RESI RASFER MOLDIG Seong aek Lim Moon Koo Kang and Woo I Lee* Department of Mechanica Engineering Seou ationa Uniersity Seou 151-74 Korea SUMMARY: he oid content within RM
More informationVariation in the spatial distribution of fish length: a multi-annual geostatistics approach on anchovy in Biscay
Internationa Counci for the Exporation of the Sea ICES / CIEM ICES CM 2003/Q:15 Variation in the spatia distribution of fish ength: a muti-annua geostatistics approach on anchoy in Biscay 1983-2002. Pierre
More informationarxiv: v1 [physics.flu-dyn] 2 Nov 2007
A theoretica anaysis of the resoution due to diffusion and size-dispersion of partices in deterministic atera dispacement devices arxiv:7.347v [physics.fu-dyn] 2 Nov 27 Martin Heer and Henrik Bruus MIC
More informationBOILING FLOW SIMULATION IN NEPTUNE_CFD AND FLUENT CODES. L. Vyskocil, J. Macek
BOILING FLOW SIMULATION IN NEPTUNE_CFD AND FLUENT CODES L. Vyskoci, J. Macek Nucear Research Institute Rez (NRI), Dept. of Therma Hydrauic Anayses, 250 68 Rez, Czech Repubic Abstract This paper presents
More informationJJMIE Jordan Journal of Mechanical and Industrial Engineering
JJMIE Jordan Journa of Mechanica and Industria Engineering Voume 11 Number 4, December. 2017 ISSN 1995-6665 Pages 217-223 Investigate the Natura Convection Heat Transfer in A PCM Therma Storage System
More informationA Better Way to Pretrain Deep Boltzmann Machines
A Better Way to Pretrain Deep Botzmann Machines Rusan Saakhutdino Department of Statistics and Computer Science Uniersity of Toronto rsaakhu@cs.toronto.edu Geoffrey Hinton Department of Computer Science
More informationNumerical simulation of a high viscosity bubble column
20th Internationa Congress on Modeing and Simuation, Adeaide, Austraia, 1 6 December 2013 www.mssanz.org.au/modsim2013 Numerica simuation of a high viscosity bubbe coumn Danio Carvajaa, Victor Meendez-Vejara,
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationNumerical Simulation for Optimizing Temperature Gradients during Single Crystal Casting Process
ISIJ Internationa Vo 54 (2014) No 2 pp 254 258 Numerica Simuation for Optimizing Temperature Gradients during Singe Crysta Casting Process Aeksandr Aeksandrovich INOZEMTSEV 1) Aeksandra Sergeevna DUBROVSKAYA
More informationGenerality of Rotating Partial Cavitation in Two-Dimensional Cascades
Proceedings of the 7 th Internationa Symposium on Cavitation CAV2009 Paper No. 90 August 17-22, 2009, Ann Arbor, Michigan, USA Generaity of Partia Cavitation in Two-Dimensiona Cascades Byung-Jin AN Osaka
More informationAN IMPROVED CORRELATION FOR TWO-PHASE PRESSURE DROP OF R-22 AND R-410A IN 180 RETURN BENDS
Proceedings of the th Braziian Congress of Therma Sciences and Engineering -- ENCIT 006 Braz. Soc. of Mechanica Sciences and Engineering -- ABCM, Curitiba, Brazi,- Dec. 5-8, 006 Paper CIT06-5 AN IMPROVED
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationCrystallisation of a supercooled spherical nodule in a flow
EUROTHERM 69 Heat and Mass Transfer in Soid-Liquid Phase hange Processes June 25-27, 2003, Bistra caste, Ljubjana, Sovenia Eds.: B. Sarer, D. Gobin rystaisation of a supercooed spherica nodue in a fow
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationSimulation of Evaporator for Two-phase Flow in the New Plate-fin Desalination Unit
Research Journa of Appied Sciences, Engineering and echnoogy 5(13): 3554-3559, 2013 ISSN: 2040-7459; e-issn: 2040-7467 Maxwe Scientific Organization, 2013 Submitted: Juy 27, 2012 Accepted: September 17,
More informationChapter 11: Two-Phase Flow and Heat Transfer Forced Convective Boiling in Tubes
11.5 Forced Convective 11.5.1 Regimes in Horizonta and Vertica Tubes The typica sequence of fow regimes for upward fow forced convective boiing in a uniformy-heated vertica tube (q =const) is shown in
More informationThe cavitation behavior with short length blades in centrifugal pump
Journa of Mechanica Science and Technoogy 24 (10) (2010) 2007~2016 www.springerink.com/content/1738-494x DOI 10.1007/s12206-010-0705-9 The cavitation behavior with short ength bades in centrifuga pump
More informationFigure Transition boiling curves
Minimum Heat Fux 10.5 Transition Boiing and Minimum Heat Fux 10.5.1 Transition Boiing Figure 10.19 Transition boiing cures Transport Phenomena in Mutiphase Systems by A. Faghri & Y. Zhang 1 Minimum Heat
More informationCopyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU
Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water
More informationAdvanced Modeling of High Speed Turbulent Reacting Flows
50th AIAA Aerospace Sciences Meeting Incuding the New Horizons Forum and Aerospace Exposition Advanced Modeing of High Speed Turbuent Reacting Fows Journa: 50th AIAA Aerospace Sciences Meeting Incuding
More informationCFD MODELLING OF DIRECT CONTACT STEAM INJECTION
Fifth Internationa Conference on CFD in the Process Industries CSIRO, Meourne, Austraia 13-15 Decemer 006 CFD MODELLING OF DIRECT CONTACT STEAM INJECTION Curtis MARSH 1 and Denis WITHERS 1 CFD Design &
More information21th International Multidimensional Engine User s Meeting at the SAE Congress 2011, April, 11, 2011 DETOIT, MI
21th Internationa Mutidimensiona Engine User s Meeting at the SAE Congress 2011, Apri, 11, 2011 DETOIT, MI "A combined Euerian Lagrangian spray amortization (ELSA) in DI Diese combustion: Fuy couped Euerian
More informationModeling of Mass transfer and Chemical Reactions in a Bubble Column Reactor Using a Discrete Bubble Model
Modeing of Mass transfer and Chemica Reactions in a Bubbe Coumn Reactor Using a Discrete Bubbe Mode Dadan Darmana 1, Nies G. Deen 1, J.A.M. Kuipers 1 1: Facuty of Science and Technoogy, University of Twente,
More informationLarge Eddy Simulation of the Gas-Liquid Flow in a Square Cross-sectioned Bubble
Large Eddy Simuation of the Gas-Liquid Fow in a Square Cross-sectioned Bubbe Coumn 1 N.G. Deen 2, T. Soberg and B.H. Hjertager Chem. Eng. Lab., Aaborg Univ. Esbjerg; Nies Bohrs Vej 8; DK-6700 Esbjerg;
More informationASummaryofGaussianProcesses Coryn A.L. Bailer-Jones
ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe
More information17 Lecture 17: Recombination and Dark Matter Production
PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationQuantitative Phase Field Modeling of Boiling Phenomena
Quantitative Phase Fied Modeing of Boiing Phenomena Arnodo Badio Pau Scherrer Institute 5232 Viigen PSI, Schweiz arnodo.badio@psi.ch ABSTRACT A quantitative Phase Fied (PF) mode is deveoped for simuating
More informationStatistical Learning Theory: A Primer
Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationEffect of transport ratio on source term in determination of surface emission coefficient
Internationa Journa of heoretica & Appied Sciences, (): 74-78(9) ISSN : 975-78 Effect of transport ratio on source term in determination of surface emission coefficient Sanjeev Kumar and Apna Mishra epartment
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationHEAT EXCHANGER NETWORK SYNTHESIS CONSIDERING CHANGING PHASE STREAMS
HEAT EXCHANGER NETWORK SYNTHESIS CONSIDERING CHANGING PHASE STREAMS F. S. Liporace a, F. L. P. Pessoa b, and E. M. Queiroz b, a PETROBRAS/CENPES/EB/SAP Cidade Universitária - Iha do Fundão 21949-900, Rio
More informationResearch Article Solution of Point Reactor Neutron Kinetics Equations with Temperature Feedback by Singularly Perturbed Method
Science and Technoogy of Nucear Instaations Voume 213, Artice ID 261327, 6 pages http://dx.doi.org/1.1155/213/261327 Research Artice Soution of Point Reactor Neutron Kinetics Equations with Temperature
More informationThe influence of temperature of photovoltaic modules on performance of solar power plant
IOSR Journa of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 Vo. 05, Issue 04 (Apri. 2015), V1 PP 09-15 www.iosrjen.org The infuence of temperature of photovotaic modues on performance
More informationTIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART
VOL. 8, No. II, DECEMBER, 0, pp 66-76 TIME DEPENDENT TEMPERATURE DISTRIBUTION MODEL IN LAYERED HUMAN DERMAL PART Saraswati Acharya*, D. B. Gurung, V. P. Saxena Department of Natura Sciences (Mathematics),
More informationExperimental and Numerical Investigation of Thermosyphon Heat Pipe Performance at Various Inclination Angles
Journa of Adanced Research in Fuid Mechanics and Therma Sciences 44, Issue 1 (2018) 85-98 Journa of Adanced Research in Fuid Mechanics and Therma Sciences Journa homepage: www.akademiabaru.com/arfmts.htm
More informationTHE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP
More informationScienceDirect. Numerical modelling of debris bed water quenching
Avaiabe onine at www.sciencedirect.com ScienceDirect Procedia IUTAM 15 2015 64 71 IUTAM Symposium on Mutiphase fows with phase change: chaenges and opportunities, Hyderabad, India December 08 December
More informationHeat Transfer Analysis of Refrigerant Flow in an Evaporator Tube
Internationa OPEN ACCESS Journa Of Modern Engineering Research (IJMER) Heat Transfer Anaysis of Refrigerant Fow in an Evaporator Tube C. Rajasekhar 1, S. Suresh, R. T. Sarathbabu 3 1, Department of Mechanica
More informationMODELING AND DESIGN OF MICRO-GROOVED FLAT PLATE EVAPORATOR
Proceedings of ICMM005 3rd Internationa Conference on Microchannes and Minichannes June 13-15, 005, Toronto, Ontario, Canada ICMM005-75141 MODELING AND DESIGN OF MICRO-GROOVED FLAT PLATE EVAPORATOR Naoki
More informationThe persistence of wake vortices analyzed through the long-term evolution of the Crow instability
The persistence of wake vortices anayzed through the ong-term evoution of the Crow instabiity Hoy JOHNSON* Vincent BRION Laurent JACQUIN WakeNet-Europe Workshop 2015 1 Introduction Cam atmospheric conditions
More informationTWO- AND THREE-DIMENSIONAL SIMULATION OF A RISING BUBBLE AND FALLING DROPLET USING LEVEL SET METHOD
European Conference on Computationa Fuid Dynamics ECCOMAS CFD 2006 P. Wesseing, E. Oñate, J. Périaux (Eds) TU Deft, The Netherands, 2006 TWO- AND THREE-DIMENSIONAL SIMULATION OF A RISING BUBBLE AND FALLING
More informationA. Distribution of the test statistic
A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationD. Prémel, J.M. Decitre and G. Pichenot. CEA, LIST, F Gif-sur-Yvette, France
SIMULATION OF EDDY CURRENT INSPECTION INCLUDING MAGNETIC FIELD SENSOR SUCH AS A GIANT MAGNETO-RESISTANCE OVER PLANAR STRATIFIED MEDIA COMPONENTS WITH EMBEDDED FLAWS D. Préme, J.M. Decitre and G. Pichenot
More informationThe basic equation for the production of turbulent kinetic energy in clouds is. dz + g w
Turbuence in couds The basic equation for the production of turbuent kinetic energy in couds is de TKE dt = u 0 w 0 du v 0 w 0 dv + g w q 0 q 0 e The first two terms on the RHS are associated with shear
More informationLarge eddy simulation of the Gas Liquid ow in a square cross-sectioned bubble column
Chemica Engineering Science 56 (2001) 6341 6349 www.esevier.com/ocate/ces Large eddy simuation of the Gas Liquid ow in a square cross-sectioned bubbe coumn N. G. Deen, T. Soberg, B. H. Hjertager Chemica
More informationDISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE
DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,
More informationhttps://doi.org/ /epjconf/
HOW TO APPLY THE OPTIMAL ESTIMATION METHOD TO YOUR LIDAR MEASUREMENTS FOR IMPROVED RETRIEVALS OF TEMPERATURE AND COMPOSITION R. J. Sica 1,2,*, A. Haefee 2,1, A. Jaai 1, S. Gamage 1 and G. Farhani 1 1 Department
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationASSESSMENT OF THE HEAT TRANSFER MODEL AND TURBULENT WALL FUNCTIONS FOR TWO FLUID CFD SIMULATIONS OF SUBCOOLED AND SATURATED BOILING
ASSESSMENT OF THE HEAT TRANSFER MODEL AND TURBULENT WALL FUNCTIONS FOR TWO FLUID CFD SIMULATIONS OF SUBCOOLED AND SATURATED BOILING D. Prabhudharwadar 1,, M. Lopez de Bertodano 1,, J. Buchanan Jr. 1 Purdue
More informationDo Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix
VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationEPOXY MEDIATED ALIGNMENT OF SINGLE WALL CARBON NANOTUBES UNDER HIGH MAGNETIC FIELDS
EPOXY MEDIATED ALIGNMENT OF SINGLE WALL CARBON NANOTUBES UNDER HIGH MAGNETIC FIELDS H. Garmestani, M.S. A-Haika, D. Lia, K. Dahmen, Georgia Institute of Technoogy, Atanta, Georgia a Department of Mechanica
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationMelodic contour estimation with B-spline models using a MDL criterion
Meodic contour estimation with B-spine modes using a MDL criterion Damien Loive, Ney Barbot, Oivier Boeffard IRISA / University of Rennes 1 - ENSSAT 6 rue de Kerampont, B.P. 80518, F-305 Lannion Cedex
More informationBoiling heat transfer of HFO-1234yf flowing in a smooth small-diameter horizontal tube
Internationa Journa of Refrigeration 3 (11) 1-153 Boiing heat transfer of HFO-13yf fowing in a smooth sma-diameter horizonta tube Shizuo Saitoha, Chaobin Dangb, Yoshitaka Nakamurab, Eiji Hiharab a Department
More informationUniversity of Huddersfield Repository
University of Huddersfied Repository Hardaupas, Y., Souopouos, N., Stetsyuk, V. and Tayor, A. M. K. P. Experimenta Assessment of subgrid scae Probabiity Density Function Modes for Large Eddy Simuation
More informationHigh Efficiency Development of a Reciprocating Compressor by Clarification of Loss Generation in Bearings
Purdue University Purdue e-pubs Internationa Compressor Engineering Conference Schoo of Mechanica Engineering 2010 High Efficiency Deveopment of a Reciprocating Compressor by Carification of Loss Generation
More informationThermophoretic interaction of heat releasing particles
JOURNAL OF APPLIED PHYSICS VOLUME 9, NUMBER 7 1 APRIL 200 Thermophoretic interaction of heat reeasing partices Yu Doinsky a) and T Eperin b) Department of Mechanica Engineering, The Pearstone Center for
More informationRecursive Constructions of Parallel FIFO and LIFO Queues with Switched Delay Lines
Recursive Constructions of Parae FIFO and LIFO Queues with Switched Deay Lines Po-Kai Huang, Cheng-Shang Chang, Feow, IEEE, Jay Cheng, Member, IEEE, and Duan-Shin Lee, Senior Member, IEEE Abstract One
More informationSimplified analysis of EXAFS data and determination of bond lengths
Indian Journa of Pure & Appied Physics Vo. 49, January 0, pp. 5-9 Simpified anaysis of EXAFS data and determination of bond engths A Mishra, N Parsai & B D Shrivastava * Schoo of Physics, Devi Ahiya University,
More informationEXPERIMENT 5 MOLAR CONDUCTIVITIES OF AQUEOUS ELECTROLYTES
EXPERIMENT 5 MOLR CONDUCTIVITIES OF QUEOUS ELECTROLYTES Objective: To determine the conductivity of various acid and the dissociation constant, K for acetic acid a Theory. Eectrica conductivity in soutions
More informationSimulation of single bubble rising in liquid using front tracking method
Advances in Fuid Mechanics VI 79 Simuation o singe bubbe rising in iquid using ront tracking method J. Hua & J. Lou Institute o High Perormance Computing, #01-01 The Capricorn, Singapore Abstract Front
More informationModelling single-bubble sonoluminescence with chemical reactions and Coulomb interactions
Modeing singe-bubbe sonouminescence with chemica reactions and Couomb interactions Ping He, Li Yuan Institute of Computationa Mathematics and LSEC Academy of Mathematics and Systems Science Chinese Academy
More informationPerformance of heat pipes as capillary pumps: modelling and comparison with experimental results
Performance of heat pipes as capiary pumps: modeing and comparison with experimenta resuts Prof. Dionissios P. Margaris 1, Zisis G. Diamantis, Dionysios I. Photeinos, Prof. Demos T. Tsahais (corresponding
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationAdjustment of automatic control systems of production facilities at coal processing plants using multivariant physico- mathematical models
IO Conference Series: Earth and Environmenta Science AER OEN ACCESS Adjustment of automatic contro systems of production faciities at coa processing pants using mutivariant physico- mathematica modes To
More informationIncorporation of surface tension to interface energy balance in crystal growth
Cryst. Res. Techno. 42, No. 9, 914 919 (2007) / OI 10.1002/crat.200710927 Incorporation of surface tension to interface energy baance in crysta growth M. Yidiz and. ost* Crysta Growth aboratory, epartment
More informationNonlinear Gaussian Filtering via Radial Basis Function Approximation
51st IEEE Conference on Decision and Contro December 10-13 01 Maui Hawaii USA Noninear Gaussian Fitering via Radia Basis Function Approximation Huazhen Fang Jia Wang and Raymond A de Caafon Abstract This
More informationMVS. Multichannel Verification System
MVS Mutichanne Verification System Increase iquid handing quaity with easy, reiabe performance verification. Are your iquid handers transferring critica voumes accuratey and precisey? Woud you know if
More informationAerodynamic Design of Turbine Blades Using an Adjoint Equation Method
43rd AIAA Aerospace Sciences Meeting and Exhibit, January 10-13, 2005, Reno, Nevada Aerodynamic Design of Turbine ades Using an Adjoint Equation Method Hsiao-Yuan Wu, and Feng Liu Department of Mechanica
More informationConservation of Circulations in Turbulent Flow
(D) Conservation of Circuations in Turbuent Fow We have emphasized the importance of deveoping a better understanding of the dynamica & statistica origin of the positivity of vortex-stretching rate ω S
More informationExperimental Study and CFD Simulation of Mass Transfer Characteristics of a Gas-Induced Pulsating Flow Bubble Column
Z. ZHOU et a., Experimenta Study and CFD Simuation of Mass Transfer, Chem. Biochem. Eng. Q. 27 (2) 167 175 (2013) 167 Experimenta Study and CFD Simuation of Mass Transfer Characteristics of a Gas-Induced
More informationTHIN FILM EVAPORATION IN MICROCHANNEL MEMBRANE FOR SOLAR VAPOR GENERATION
Proceedings of the th Internationa Heat Transfer Conference, IHTC- August -,, Kyoto, Japan IHTC-9 THIN FILM EVAPORATION IN MICROCHANNEL MEMBRANE FOR SOLAR VAPOR GENERATION Ammar A. Asheghri, TieJun Zhang
More informationUniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete
Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity
More informationCalculating Alkalinity
Cacuating Akainity Overview How to cacuate acid:anion ratios at varying ionic strengths Use the Mixing Water cacuation as a ph titration apparatus Understand how weak acid chemistry affects akainity Compare
More informationThe mechanical energy balance equation used for the mh-b correlation 1 (2-6) sg u
Modified Haedron and Brown Method (mh-b) This is an empirica two-phase fow correation, the core of which is correation for iquid hod-up. Griffith correation is used for fow in the bubbe fow reion. The
More informationThe structure of turbulent flames in fractal- and regulargrid-generated
The structure of turbuent fames in fracta- and reguargrid-generated turbuence T Sponfedner 1, N Souopouos 1, F Beyrau 1, Y Hardaupas 1, AMKP Tayor 1,* and JC Vassiicos 2 1 Department of Mechanica Engineering,
More information