Sharp-Interface Simulations of Heat-Driven Flows with Marangoni Effect
|
|
- Avice Greene
- 5 years ago
- Views:
Transcription
1 th International Smposium on Turbulence and Shear Flow Phenomena (TSFP), Chicago, USA, Jul, 27 Sharp-Interface Simulations of Heat-Driven Flows with Marangoni Effect Johannes L. Pieper Felix S. Schranner Gurvan Hermange Nikolaus A. Adams Abstract A method to account for the effect of interfacial surfactant concentration on multi-phase flows has been proposed b Schranner & Adams (26). The therein validated method builds up on the the sharp-interface method of Hu et al. (26). This work advances the method to simulate and analse configurations of industrial applicabilit. In the majorit of these configurations temperature gradients, commonl due to local heat sources or sinks, are the dominant cause of the Marangoni flow. In order to account for the Marangoni effect properl, the temperature distribution within the two-phase flows under consideration has to be known. In this work, two numerical approaches for the calculation of the temperature are presented and evaluated. In the first approach the temperature transport equation, realized as a scalar transport equation, is coupled to the finite volume algorithm of Schranner & Adams (26). The underling method for the second approach is the weakl compressible high-resolution approach for incompressible flows (Schranner et al. (23)), which is generalized to compressible and incompressible flows. The extension to multi-phase flows with discontinuous phases is formulated in this work. The approaches are evaluated and compared on basis of the test flow proposed b Pendse & Esmaeeli (2). It is sufficientl complex in its evolution, et simple in its set-up and permits comparison to an analtical reference solution. The application of the methods to further thermocapillar flows is pursued. Both methods are robust and correctl predict the evolution of the Marangoni flows and temperature fields. Numerical experiments indicate that the passive scalar method is faster than the energ equation method. Yet, for two-phase flows with fluids of ver diverse properties the energ equation method is more accurate. Introduction In man classical two-phase flows, such as viscous fingering (Matar & Troian (999)), drop break-up and coalescence (Leal (24)), tip-streaming (Bruijn (993)), and buoanc-driven bubble-motion (Fdhila & Duineveld (996)), the effects of locall variable surface tension are significant. The causes are inhomogeneousl distributed interfacial surfactants and temperature gradients. The accurate and concurrentl robust two-phase flow simulation of immiscible fluids with large densit and viscosit ratios in combination with resolved interface dnamics can be numericall challenging. Hu et al. (26) have developed a full conservative, Eulerian, level-set based sharp-interface method (SIM) to simulate immiscible, compressible flows. The sharp-interface framework has demonstrated to accuratel predict two-phase flows, with each phase obeing different equations of state, and large densit and viscosit ratios (Lauer et al. (22); Schranner et al. (26)). A viscous extension of the SIM for incompressible two-fluid flows is due to Luo et al. (25). A robust, more accurate, simpler and eas to implement further-development of the SIM that incorporates inviscid, viscous, capillar and Marangoni interfacial stresses is due to Schranner & Adams (26). Generall, the SIM has demonstrated to be a reliable and simple approach to capture multi-fluid interfaces robustl even for complex interface topolog changes (Luo et al. (25); Hu et al. (26); Lauer et al. (22); Schranner et al. (26)). To account for interfacial surfactant transport a method has been developed and applied for the diffuse interface approach b Teigen et al. (2), Adalsteinsson & Sethian (23), and Xu & Zhao (23). It has been improved in terms of accurac, surfactant mass conservation and computational efficienc b Schranner & Adams (26), emploing the sharp-interface methodolog. To allow for the simulation of temperature-gradient driven Marangoni flows, the transport of mass, momentum and energ within each phase and across the interface needs to be considered appropriatel. Based on the accurate interfacial temperature distribution, the local capillar and Marangoni stresses as well as their contribution to the momentum and energ balance of the involved fluids can be computed straightforwardl. This generalization is pursued within the scope of this work. Therefore, two numerical approaches for advancing the temperature within the phases and exchanging energ across the interface are presented. The are evaluated on basis of a selection of test flows. Furthermore, the industrial application are sought. Basic Concepts The finite-volume algorithm of Schranner et al. (23) for weakl compressible turbulent and non-turbulent single-phase flows, further-developed to multi-phase configurations with capillar and Marangoni effects Schranner & Adams (26) is the underling numerical method. For incompressible flows the energ equation is decoupled from the transport of mass and momentum Panton (984) and becomes obsolete. Yet, it is not obsolete if internal energ or temperature distributions are of interest.
2 Passive scalar temperature transport equation (PS) For incompressible multi-phase flows with temperature gradients and negligible temporal pressure fluctuations the temperature equation ( ) T ρc p t + u T = λ T + λ T () 2 3 µ ( u)2 + 2µ (( u) : ( u)) + βtu p + ω holds within the two-fluid domain Ω = Ω ξ Ω ξ 2. The heat capacit c p, coefficient of thermal expansion β, thermal conductivit λ, dnamic viscosit µ = ρν are material properties of the fluids ξ and ξ 2. Numericall, within the multi-fluid cells, which are those cells containing both fluids, the material properties are determined from the volume fractions ζ ξ i(t) = V dv of the fluid phases Ω Ω ξ i ξ i and the individual fluid properties: e.g. c p = (ζ c p ) ξ + (ζ c p ) ξ 2. The derivatives of Eq. () are approximated with 4 th -order central (LHS: st term, RHS: st, 2 nd, 3 rd, 4 th term) and 5 th -order upwindbiased (LHS: 2 nd term, RHS: 5 th term) finite difference schemes, respectivel. Note, however, from numerical experiments we find that the order of convergence is independent from the choice of the order for the central differences, thus due to the order of accurac of the scheme for fluid transport. The 6 th term of Eq. () defines volume and area sources. To consider an area-distributed heat source (AHS) the power of a radiative heat source P is considered to act uniforml on a certain area AIF tot. The product with the ratio of interface-area AIF [i, j,k] to volume V [i, j,k] of a cell [i, j,k] is the heat added to each cell at the interface: ω = P AIF tot AIF [i, j,k] V [i, j,k]. Energ equation method (EEM) The method of Schranner & Adams (26) is extended to incorporate the energ equation for compressible as well as weakl compressible flows, see also Schranner (27). Hence, within each fluid phase the transport equations for the conservative states U = (ρ,ρu,ρe t ) are solved seperatel. e t = e kin + e i. B considering the interface-interaction flux in each multi-fluid cell X ξ i = {,X M,ξ i,x E,ξ i } (2) the two subdomains Ω ξ i containing one of the fluids ξ i are coupled and the sstem is globall conservative. The momentum flux is: X M,ξ i = {X,ξ i + X,ξ i} (Schranner & Adams (26)). In contrast to Schranner (27) conductive heat transfer is also considered in this work. Thus, the interfacial energ flux is X E,ξ i = X M,ξi u ξ i I + q I n ξ i I. (3) Thereb, the interface heat flux is q I = q I,α = λ I,α D(T ) I,α (4) with α =,2,3 denoting the Cartesian directions. Numericall, the temperature gradient across the interface in direction α, denoted as D(T ) I,α, is approximated with a 4 th -order central difference. The stencil is constructed of T = {(ζ T ) ξ +(ζ T ) ξ 2}. To account for the discontinuous λ within interface cells, let λ I,α = ). (5) D( xα λ(u) du The order of the approximation of D of Eq. (5) is consistent with the order of D(T ) I,α. The interface temperature T I is obtained from e i = c p T I. Temperature dependent material properties Temperature dependencies of the material properties are considered. For incompressible or weakl compressible fluids, the Boussinesq approximation is used to account for a temperature dependent densit in the buoant force, modeled as a source term: ρ = ρ re f ( β T ( T Tre f )), (6) where β T is the coefficient of thermal expansion. The temperature dependent surface-tension σ is modeled with a linear Langmuir constitutive equation: σ(t ) = σ re f + σ ( ) T Tre f. (7) T σ re f is the surface tension coefficient at a reference temperature T re f. These as well as the temperature gradient σ are constant T and material-specific. To consider solidification and melting phase changes the dnamic viscosit is considered constant and relativel high for temperatures lower than a specific the melting temperature T S. For higher temperatures, the dnamic viscosit obes an Arrhenius function. Continuous transition is achieved b means of a Heaviside function. µ max, T < T S 2 δt (µ max µ(t S )) µ(t ) = ( H (T T S + δt,δt)), T S 2 δt T S (8) µ exp( E T T S ), T T S. The temperature dependenc of the specific heat capacit c p includes a normall distributed latent heat h S around the melting point, the width is T : ( ) ) 2 exp ( T TS T c p (T ) = c p + h S. (9) π T Thermocapillar Flow To validate the two models and verif their implementation a thermocapillar driven creeping flow with two superimposed planar fluids is considered, see Pendse & Esmaeeli (2). The analtical solutions for the temperature T field, and the velocities u and v are available. The geometrical set-up is shown in figure. Figure. Geometrical Set-up for thermocapillar flow. δ = δ 2 =. and l = 4δ are used. The boundar condition are periodic in x -direction. In x 2 -directions a no-slip condition is assumed. The uniform temperature T (x,δ 2 ) = T c is imposed at the top boundar. Along the bottom boundar a sinusoidal temperature distribution with T (x, δ ) = T h + T cos( 2π l x ) is assumed. T h = 2, T = 4, T c =. The stead-state solution ma be found in Pendse & Esmaeeli (2) or Schranner & Adams (26). The kinematic viscosities are ν ξ = ν ξ 2 =.2. Both fluids obe the Tait equation of state, with γ =.. The speed of sound depends on the densit. For fluids with ρ =. the speed of sound is a =, for fluids with ρ = the speed of sound is a =.. Moreover, σ T = 5 4 and σ =.3 at T re f =. Four different cases are considered. For the upper fluid ρ ξ 2 =., c ξ 2 = 2/3 and λ ξ 2 =.2. The densit, heat capacit and the thermal conductivit of the lower fluid are case-dependent:
3 Case : ρ ξ = ρ ξ 2, c ξ = c ξ 2 and λ ξ = λ ξ 2, Case 2: ρ ξ = ρ ξ 2, c ξ = 5c ξ 2 and λ ξ = 5λ ξ 2, Case 3: ρ ξ = ρ ξ 2, c ξ = c ξ 2/ and λ ξ = λ ξ 2, Case 4: ρ ξ = ρ ξ 2, c ξ = c ξ 25/ and λ ξ = 5λ ξ 2. The reference length and velocit are L re f = δ =. and U re f = σ T T δ µ ξ l = 4, respectivel. Thus, one finds that the Renolds, Marangoni, Peclet, capillar, and Weber numbers are Re= U re f L re f =.25, Ma= U re f ρ ξ 2 c ξ 2 l = ν ξ 2 λ ξ 2 3, Pe= U re f L re f ρ ξ 2 c ξ 2 = λ ξ 2 2, Ca= ρξ 2U re f ν ξ 2 σ = 6 We= ρξ 2U re 2 f L re f σ = 48. These numbers compl Re<<, Ca<< and Ma<<, hence the analtical solution of Pendse & Esmaeeli (2) hold. To evaluate the simulation, the L - L - and L 2 - absolut errors are evaluated. The absolute error norms are calculated for u, v and T L = max M ( ε ), L = V ε dv { ε V } [i, j], M ( ( ) /2 L 2 = ε 2 dv V ) /2 { ε 2 V } [i, j]. M The is reduced in five steps from.2 to.625. (a) (b) (c) Case The first case serves as a basic verification and validation configuration. The temperature and flow field is depicted in Fig. 2. Figure 3 shows the L errors for the EEM and PS approaches. For both methods the results converge to the analtical solution independent from the chosen approach. Regarding T, L is below.49 for PS and.266 for EEM. The maximum error is lower than.2 for both methods at the largest resolution. The largest deviation of u and v, found on the interface, is (PS) and (EEM),.8 5 (PS) and.74 5 (EEM), respectivel. With respect to the maximum velocities of u max = 3 4 and v max = 6 5 these values indicate good results even for the lowest resolution. Spatial refinement leads to a reduction in the errors. The convergence order is between. and.3 for the velocities, and. < L (T ) <.6 for T x Figure 2. Results (PS) for Case with x 2 at stead state, t = 2. Case 2 The second case simulates two fluids with the same densit and different thermal conductivities and capacities. This corresponds to a highl simplified model of a laser beam welding process, as air and metal have a large difference in their thermal conductivities and capacities. Figure 4 depicts the L errors for the EEM and PS. L 2 and L are not shown for clarit. The simulation results for both methods are in a good agreement with the analtical solution. EEM leads to more accurate results than the results of PS. This is found from evaluating all three error norms. A possible explanation ma be that while the EEM calculates the T L x. x.6 x.3 x. x x.8.5 L,u,PS L,u,EEM L,v,PS L,v,EEM L,T,PS L,T,EEM Figure 3. L error norms for u, v and T for Case with PS and EEM at stead state, t = 2. energ for both fluids and models the heat flux across the sharp interface, the PS calculates the temperature for the entire domain without explicitl considering the interface. The deviations from the analtical solution are larger with the PS. Consequentl, larger deviations in the temperature gradient are found which lead to different results of the flow velocities. The maximum velocities of the analtical solution are u max,anl and v max,anl.2 4. The simulated maximum velocities are u max,ps , v max,ps.24 4, u max,eem , v max,eem.6 4 found at a mesh width of x.5. The accurac of the first case is not achieved for both methods. The temperature order of convergence is about.2 for PS and. for EEM. L x. x.2 x.3 x.3 x x..5 L,u,PS L,u,EEM L,v,PS L,v,EEM L,T,PS L,T,EEM Figure 4. L error norms for velocit u, v and temperature T for Case 2 with PS and EEM at stead state, t = 2. Case 3 This case extends Case b considering a densit of. Figure 5 depicts the L errors for the temperature at different resolutions. The velocit contours are depicted in figure 6 The temperature distribution approaches the analtical solution well. The maximum errors, e.g. L, are identical to Case for both methods and same resolutions. In the upper fluid, Ω ξ 2, which has a lower densit, significantl higher velocities are observed in the stationar state than in Ω ξ
4 L 2 3 x. x.6 L,T,PS L,T,EEM Case 4 Case 4 is a combination of Cases 2 and 3. It is a configuration that resembles laser beam welding better than the previous cases. Figure 7 depicts the L temperature error norms for when emploing EEM and PS. The temperature fields as well as the resolution dependent error norms are found to be identical to those of Case 2. The velocities are found to differ significantl from the analtical solution. A possible explanation is given for Case 3, which also holds for Case 4. x Figure 5. L error norms for T for Case 3 with PS and EEM at stead state, t = 2. L x. L,T,PS L,T,EEM independent from when emploing EEM or PS According to Pendse & Esmaeeli (2) the densit has no influence on the solution as long as Re<<, Ca<< and Ma<<. The ratio of the two dnamic viscosities, which is reduced from (Case ) to, has onl an influence on the maximum velocit in the analtical solution. Pendse & Esmaeeli (2) do not impose an restriction to the densit ratios. The question is raised, at this venue, whether the analtical solution does or does not appl to combinations of fluids with ver large densit ratios, although the creeping flow conditions are fulfilled for one of the two fluid. Within ξ (the lower fluid) Ca = 6 for Case 3, for Cases and 2 it is Ca = 6. Hence, ξ cannot be considered a creeping flow, even though ξ 2 can. Figure 7. L error norms for T for Case 4 with PS and EEM at stead state, t = 2. Thermocapillar motion of a Drop 5E 6 5E 6.5 5E 6 5E 6.5E 5.5E 5 5E 7 5E 7.5 E 7 E 7 E 8 E x 5E 7.5 5E 6 E 8 (a) u-velocit field 5E 7 3E 63E 6 3E 6 5E 6 E 7 E 7 E 8 E 8 E 8.5 E 8 5E 7 3E 6 E 7 E x (b) v-velocit field Figure 6. Velocities u and v with a densit-ratio of at stead state, t = 2. Figure 8. Principal sketch of the Set-up for Rising Bubble. Another test case for the validation of the implemented equations is provided b Young et al. (959), who have investigated the velocit of a two-dimensional rising drop. For this purpose, a drop with the radius R is set in a clinder with the height L and the radius 5R (see figure 8). Furthermore, a linear temperature field is specified over the height of the clinder. The temperature at the bottom of the clinder is T ( = ) = T H and at the top T ( = L) = T C. This results in a linear temperature profile. The surface tension changes linearl with the height according to the temperature distribution. For validation, the average ascent rate is calculated and compared to the ascent rate u Y GB according to Young et al. (959). The empiricall evaluated ascent rate is ) 2( σ T TH T C L R ρgr 2 µ +µ 2 µ u Y GB =. () 9µ + 6µ 2
5 The average ascent rate of the simulated drop is determined b ū 2 = u 2 dv = V 2 Ω 2 V 2 u 2 ζ ξ 2 V, (2) M where M is the number of cells, V 2 is the volume of the drop, and ζ ξ 2 is the volume fraction of the disperse phase. Furthermore, R =.5 and L = 7.5R. For the fluids, ρ ξ = ρ ξ 2 =.2, a ξ = 5., a ξ 2 = 4. The Tait equation of state holds for both fluids. The initial pressure for both fluids is equal to p =. With the dnamic viscosit µ ξ = µ ξ 2 =., g =, σ T = and the temperatures T C = 3 and T H = 3.5, a theoretical ascent rate of u Y GB = 45 4 is obtained. Thus, the Renolds number is Re= This configuration is simulated with EEM and PS. The results are nearl identical. For clarit, we onl provide results for PS in Figs. 9,. A convergence analsis is performed with a spatial resolution ranging from to 256, four different refinement levels are simulated. For a resolution of x = the deviation from the theoretical ascent rate is alread less than %. The results can be considered to be converged for a resolution of x = as these do not differ for a grid of x = The predicted ascent rates of the drop and the pressures agree well with the empirical solutions. ū2 uy GB pressure p Figure time t Normalised ascent velocit for different resolutions theoretical solution simulation 256 Figure. Simulated pressure at x =.25 and t = compared to theoreticall expected pressure Application Example - Laserdeep-Welding Based on Case 4, a simple set-up to simulate heat conduction welding is developed. Copper is chosen to be the welding material and air as the ambient gas. The geometrical set-up is l =. and δ = δ 2 =.25, comparable to figure. Moreover, ρ ξ 2 =.9, c ξ 2 = 4, µ ξ 2 = and λ ξ 2 =.262, ρ ξ = 896, c ξ = 385 and λ ξ = 394. The latent heat of copper is h ξ S = 27, it is distributed normall with T = 25 around the melting point T ξ S = 357, see Eq. (9). The dnamic viscosit is modeled according to Eq. (8), µ ξ max =, µ ξ (T S ) =.3 and E ξ = The heat source is a laser beam, modeled as an AHS. It acts within a small region with an focus diameter of d f =.3 and a powerdensit of 3 9. To account for thermal convection, ρ(t ) is modeled according to Eq. (6). To model the temperature dependent surface tension Eq. (7) holds. The set-up is shown in Figure. The feed Figure. Geometrical setup of the laser welding configuration. rate is zero in order to examine the convection in the molten pool. We simulate the initial heat-up process. The flow within the molten metal pool and air are visualized in Figs. 2 and 3, respectivel. Note that the molten metal pool is the area enclosed b the T = T S - line, the temperature within the enclosed area is higher than T S. The two-dimensional model with the selected parameters shows that the melting temperature is reached after s. A Marangoni flow develops, driving the liquid phase. The maximum velocities are observed at the interface where the melting temperature T S = 357 is reached. The comparabl higher velocities in the air are due to the comparabl low dnamic viscosit of the air. Figure 2. Fluid flows in molten metal pool at t =.5 3. Conclusion Two methods to simulate two-phase immiscible, incompressible flows with temperature dependent material properties have been presented. Both methods are applicable to simulate Marangoni flows with a sharp-interface. In the first method temperature is propagated as a passive scalar, multifluid cells are considered b mixing fluid properties, the transport equations of mass and momentum appl for the fluids in motion, the exchange of momentum across the interface couples the two phases. The second method emplos the transport equations for mass, momentum and energ and couples the immiscible phases b means of interface-fluxes, modeling the
6 Figure 3. Fluid flow in ambient gas at t =.5 3. exchange of momentum and energ across the interface. B comparing the results to analtical solutions of a selection of test flows the methods are verified and their implementations validated. Both methods are robust and correctl predict the evolution of the Marangoni flows and temperature fields. The simulation results of a thermocapillar flow are in excellent agreement with analtical data for both methods. Numerical experiments indicate that the passive scalar method is faster than the energ equation method. Yet, for two-phase flows with fluids of ver diverse properties the energ equation method is more accurate. As an industrial application example a two-dimensional base configuration of a laser beam welding process has been simulated. Acknowledgment We acknowledge the Deutsche Forschungsgemeinschaft (DFG) for funding. Felix S. Schranner is member of the TUM Graduate School. The Munich Centre of Advanced Computing (MAC) has provided the computational resources. REFERENCES Adalsteinsson, David & Sethian, J.A. 23 Transport and diffusion of material quantities on propagating interfaces via level set methods. Journal of Computational Phsics 85 (), Bruijn, R.A. De 993 Tipstreaming of drops in simple shear flows. Chemical Engineering Science 48 (2), Fdhila, R. Bel & Duineveld, P.C. 996 The effect of surfactant on the rise of a spherical bubble at high renolds and peclet numbers. Phsics of Fluids 8 (2), Hu, XY, Khoo, BC, Adams, Nikolaus A & Huang, FL 26 A conservative interface method for compressible flows. Journal of Computational Phsics 29 (2), Lauer, E., Hu, X. Y., Hickel, S. & Adams, N. A. 22 Numerical investigation of collapsing cavit arras. Phsics of Fluids 24 (5). Leal, L. G. 24 Flow induced coalescence of drops in a viscous fluid. Phsics of Fluids 6 (6), Luo, J., Hu, X.Y. & Adams, N.A. 25 A conservative sharp interface method for incompressible multiphase flows. Journal of Computational Phsics 284 (), Matar, Omar K. & Troian, Sandra M. 999 The development of transient fingering patterns during the spreading of surfactant coated films. Phsics of Fluids (), Panton, R.L. 984 Incompressible Flow. John Wile & Sons Australia, Limited. Pendse, Bhushan & Esmaeeli, Asghar 2 An analtical solution for thermocapillar-driven convection of superimposed fluids at zero renolds and marangoni numbers. International Journal of Thermal Sciences 49 (7), Schranner, Felix Sebastian 27 Weakl compressible models for complex flows. PhD thesis,. Schranner, Felix S. & Adams, Nikolaus A. 26 A conservative interface-interaction model with insoluble surfactant. Journal of Computational Phsics 327, Schranner, Felix. S., Hu, Xiangu & Adams, Nikolaus A. 23 A phsicall consistent weakl compressible high-resolution approach to underresolved simulations of incompressible flows. Computers & Fluids 86, Schranner, Felix S., Hu, Xiangu & Adams, Nikolaus A. 26 On the convergence of the weakl compressible sharp-interface method for two-phase flows. Journal of Computational Phsics 324, Teigen, Knut Erik, Song, Peng, Lowengrub, John & Voigt, Axel 2 A diffuse-interface method for two-phase flows with soluble surfactants. Journal of Computational Phsics 23 (2), Xu, Jian-Jun & Zhao, Hong-Kai 23 An eulerian formulation for solving partial differential equations along a moving interface. Journal of Scientific Computing 9 (-3), Young, N. O., Goldstein, J. S. & Block, M. J. 959 The motion of bubbles in a vertical temperature gradient. Journal of Fluid Mechanics 6 (3),
Chapter 2 Basic Conservation Equations for Laminar Convection
Chapter Basic Conservation Equations for Laminar Convection Abstract In this chapter, the basic conservation equations related to laminar fluid flow conservation equations are introduced. On this basis,
More information7. TURBULENCE SPRING 2019
7. TRBLENCE SPRING 2019 7.1 What is turbulence? 7.2 Momentum transfer in laminar and turbulent flow 7.3 Turbulence notation 7.4 Effect of turbulence on the mean flow 7.5 Turbulence generation and transport
More informationModelling Axial Flow Generation and Profile Evolution in the CSDX Linear Device R. J. Hajjar, P. H. Diamond, G. R. Tynan, R. Hong, S.
Modelling Axial Flow Generation and Profile Evolution in the CSDX Linear Device R. J. Hajjar, P. H. Diamond, G. R. Tnan, R. Hong, S. Thakur This material is based upon work supported b the U.S. Department
More informationInternational Journal of Applied Mathematics and Physics, 3(2), July-December 2011, pp Global Research Publications, India
International Journal of Applied Mathematics and Phsics, 3(), Jul-December 0, pp. 55-67 Global Research Publications, India Effects of Chemical Reaction with Heat and Mass Transfer on Peristaltic Flow
More informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationVelocity Limit in DPD Simulations of Wall-Bounded Flows
Velocit Limit in DPD Simulations of Wall-Bounded Flows Dmitr A. Fedosov, Igor V. Pivkin and George Em Karniadakis Division of Applied Mathematics, Brown Universit, Providence, RI 2912 USA Abstract Dissipative
More informationComputational Fluid Dynamics (CFD, CHD)*
1 / 1 Computational Fluid Dnamics (CFD, CHD)* PDE (Shocks 1st); Part I: Basics, Part II: Vorticit Fields Rubin H Landau Sall Haerer, Producer-Director Based on A Surve of Computational Phsics b Landau,
More informationChapter 4 Transport of Pollutants
4- Introduction Phs. 645: Environmental Phsics Phsics Department Yarmouk Universit hapter 4 Transport of Pollutants - e cannot avoid the production of pollutants. hat can we do? - Transform pollutants
More informationCONSERVATION OF ENERGY FOR ACONTINUUM
Chapter 6 CONSERVATION OF ENERGY FOR ACONTINUUM Figure 6.1: 6.1 Conservation of Energ In order to define conservation of energ, we will follow a derivation similar to those in previous chapters, using
More informationA FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS NUMBER FLUID FLOWS
Mathematical Models and Methods in Applied Sciences Vol. 12, No. 4 (2002) 525 539 c World Scientific Publishing Compan A FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS
More informationScale-Adaptive Simulation (SAS) Turbulence Modeling. F.R. Menter, ANSYS Germany GmbH
Scale-Adaptive Simulation (SAS) Turbulence Modeling F.. Menter, ANSYS German GmbH 1 Unstead ANS Based Models UANS (Unstead enolds averaged Navier Stoes) Methods UANS gives unphsical single mode unstead
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at.sciencedirect.com ScienceDirect Procedia Engineering 9 (14 383 388 1th International Conference on Mechanical Engineering, ICME 13 Effects of volumetric heat source and temperature
More informationUnsteady free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion *
CHAPTER 4 Unstead free MHD convection flow past a vertical porous plate in slip-flow regime under fluctuating thermal and mass diffusion * 4. Nomenclature A : Suction parameter C : Dimensionless species
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationJournal of Colloid and Interface Science
Accepted Manuscript On the phase-field modelling of a miscible liquid/liquid boundar Ruilin Xie, Anatoli Vorobev PII: S21-9797(15)3336-2 DOI: http://dx.doi.org/1.116/j.jcis.215.11.26 Reference: YJCIS 2877
More informationInternational Journal of Scientific Research and Reviews
Research article Available online www.ijsrr.org ISSN: 2279 543 International Journal of Scientific Research and Reviews Soret effect on Magneto Hdro Dnamic convective immiscible Fluid flow in a Horizontal
More informationEntropy ISSN
Entrop 5, 7[4], 3 37 3 Entrop ISSN 99-43 www.mdpi.org/entrop/ Full Paper Second law analsis of laminar flow in a channel filled with saturated porous media: a numerical solution Kamel Hooman, Arash Ejlali
More informationJoule Heating Effects on MHD Natural Convection Flows in Presence of Pressure Stress Work and Viscous Dissipation from a Horizontal Circular Cylinder
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 7-3, 04. Available online at www.jafmonline.net, ISSN 735-357, EISSN 735-3645. Joule Heating Effects on MHD Natural Convection Flows in Presence of
More informationLarge Eddy Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes
Int. J Heat Mass Transfer, Vol. 76 pp. 432-446, 214 Large Edd Simulation of the flow and heat transfer in a half-corrugated channel with various wave amplitudes Maram Mirzaei, Ahmad Sohankar, Lars Davidson,
More informationFilm thickness Hydrodynamic pressure Liquid saturation pressure or dissolved gases saturation pressure. dy. Mass flow rates per unit length
NOTES DERITION OF THE CLSSICL REYNOLDS EQTION FOR THIN FIL FLOWS Te lecture presents te derivation of te Renolds equation of classical lubrication teor. Consider a liquid flowing troug a tin film region
More informationIntroduction and Fundamental Concepts (Lectures 1-7)
Introduction and Fundamental Concepts (Lectures -7) Q. Choose the crect answer (i) A fluid is a substance that (a) has the same shear stress at a point regardless of its motion (b) is practicall incompressible
More informationA NEW APPROACH TO TREAT THE RANS-LES
ETMM: th International ERCOFTAC Smposium on Engineering, Marbella, 7-9 September 4 A NEW APPROACH TO TREAT THE RANS-LES INTERFACE IN PANS Lars Davidson Div. of Fluid Dnamics, Dept. of Applied Mechanics
More informationInvestigation of Cavitation Bubble Cloud with discrete Lagrangian Tracking
Investigation of Cavitation Bubble Cloud with discrete Lagrangian Tracking Xiuxiu Lyu, Xiangyu Hu, Nikolaus A. Adams Chair of Aerodynamics and Fluid Mechanics, Department of Mechanical Engineering, Technical
More informationIncompressible Viscous Flows
Incompressible Viscous Flows Q. Choose the crect answer (i) The maximum velocit of a one-dimensional incompressible full developed viscous flow between two fixed parallel plates is 6m/s. The mean velocit
More informationESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundary Layer
ESCI 485 Air/sea Interaction Lesson 5 Oceanic Boundar Laer References: Descriptive Phsical Oceanograph, Pickard and Emer Introductor Dnamical Oceanograph, Pond and Pickard Principles of Ocean Phsics, Apel
More informationPHYSICAL REVIEW E 78, Lattice Boltzmann method for simulation of compressible flows on standard lattices
PHYSICAL REVIEW E 7, 016704 00 Lattice Boltzmann method for simulation of compressible flows on standard lattices Nikolaos I. Prasianakis 1, * and Ilia V. Karlin 1,, 1 Aerothermochemistr and Combustion
More informationVisualizing Non-Equilibrium Flow Simulations using 3-D Velocity Distribution Functions
Visualizing Non-Equilibrium Flow Simulations using 3-D Velocit Distribution Functions A. Venkattraman and A.A. Alexeenko School of Aeronautics & Astronautics, Purdue Universit, West Lafaette IN 797 Abstract.
More informationA Direct Derivation of the Griffith-Irwin Relationship using a Crack tip Unloading Stress Wave Model.
A Direct Derivation of the Griffith-Irwin Relationship using a Crack tip Unloading Stress Wave Model. C.E. Neal-Sturgess. Emeritus Professor of Mechanical Engineering, The Universit of Birmingham, UK.
More informationA first contact with STAR-CCM+ Comparing analytical and finite volume solutions with STAR-CCM+ simulations
A first contact with STAR-CCM+ Comparing analtical and finite volume solutions with STAR-CCM+ simulations Michael Heer Analtical Finite volumes STAR-CCM+ What is ParisTech? ParisTech is a consortium of
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationA high order Padé ADI method for unsteady convection-diffusion equations
Center for Turbulence Research Annual Research Briefs 2005 85 A high order Padé ADI method for unstead convection-diffusion equations B D. You 1. Motivation and objectives The unstead convection-diffusion
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationVibration Analysis of Isotropic and Orthotropic Plates with Mixed Boundary Conditions
Tamkang Journal of Science and Engineering, Vol. 6, No. 4, pp. 7-6 (003) 7 Vibration Analsis of Isotropic and Orthotropic Plates with Mied Boundar Conditions Ming-Hung Hsu epartment of Electronic Engineering
More informationBubble curtain system for prevention of seawater intrusion in coastal aquifers
Available online at www.pelagiaresearchlibrar.com Advances in Applied Science Research, 013, 4(4):81-90 ISSN: 0976-8610 CODEN (USA): AASRFC Bubble curtain sstem for prevention of seawater intrusion in
More informationMathematical Modeling on MHD Free Convective Couette Flow between Two Vertical Porous Plates with Chemical Reaction and Soret Effect
International Journal of Advanced search in Phsical Science (IJARPS) Volume, Issue 5, September 4, PP 43-63 ISSN 349-7874 (Print) & ISSN 349-788 (Online) www.arcjournals.org Mathematical Modeling on MHD
More informationModule 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction
1 Module 1 : The equation of continuit Lecture 4: Fourier s Law of Heat Conduction NPTEL, IIT Kharagpur, Prof. Saikat Chakrabort, Department of Chemical Engineering Fourier s Law of Heat Conduction According
More informationChE 385M Surface Phenomena University of Texas at Austin. Marangoni-Driven Finger Formation at a Two Fluid Interface. James Stiehl
ChE 385M Surface Phenomena University of Texas at Austin Marangoni-Driven Finger Formation at a Two Fluid Interface James Stiehl Introduction Marangoni phenomena are driven by gradients in surface tension
More informationVISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL CHANNEL THROUGH A POROUS MEDIUM
Volume 2, No. 1, Januar 214 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info VISCO-ELASTIC FLUID FLOW WITH HEAT AND MASS TRASNFER IN A VERTICAL
More informationOscillatory heating in a microchannel at arbitrary oscillation frequency in the whole range of the Knudsen number
Journal of Phsics: Conference Series Oscillator heating in a microchannel at arbitrar oscillation freuenc in the whole range of the Knudsen number To cite this article: O Buchina and D Valougeorgis J.
More informationDirect Numerical Simulation of Two-Phase Flows with phase change
Direct Numerical Simulation of Two-Phase Flows with phase change Sébastien Tanguy, Lucia Rueda Villegas, Catherine Colin, Grégory Huber Institut de Mécanique des Fluides de Toulouse Maître de conférence
More informationHigh resolution of 2D natural convection in a cavity by the DQ method
Journal of Computational and Applied Mathematics 7) 9 6 www.elsevier.com/locate/cam High resolution of D natural convection in a cavit b the DQ method D.C. Lo a,, D.L. Young b, C.C. Tsai c a Department
More informationSimulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundary Element Method
10 th International LS-DYNA Users Conference Simulation Technolog (2) Simulation of Acoustic and Vibro-Acoustic Problems in LS-DYNA using Boundar Element Method Yun Huang Livermore Software Technolog Corporation
More informationExperimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper
Tamkang Journal of Science and Engineering, Vol. 8, No 4, pp. 319-35 (005) 319 Experimental Calibration and Head Loss Prediction of Tuned Liquid Column Damper Jong-Cheng Wu Department of Civil Engineering,
More informationDNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT
10 th International Symposium on Turbulence and Shear Flow Phenomena (TSFP10), Chicago, USA, July, 2017 DNS STUDY OF TURBULENT HEAT TRANSFER IN A SPANWISE ROTATING SQUARE DUCT Bing-Chen Wang Department
More informationUNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX METHOD
Proceedings of the 3rd ASME/JSME Joint Fluids Engineering Conference Jul 8-23, 999, San Francisco, California FEDSM99-8 UNSTEADY LOW REYNOLDS NUMBER FLOW PAST TWO ROTATING CIRCULAR CYLINDERS BY A VORTEX
More informationSTUDY OF HEAT TRANSFER MECHANISMS DURING THE LENS TM PROCESS
STUDY OF HEAT TRANSFER MECHANISMS DURING THE LENS TM PROCESS Liang Wang 1 and Sergio Felicelli 1. Center for Advanced Vehicular Systems, Mississippi State University, Mississippi State, MS 3976, USA; email:
More informationNumerical Analysis of Heat Transfer in the Unsteady Flow of a non- Newtonian Fluid over a Rotating cylinder
Bulletin of Environment, Pharmacolog and Life Sciences Bull. Env.Pharmacol. Life Sci., Vol 4 [Spl issue 1] 215: 318-323 214 Academ for Environment and Life Sciences, India Online ISSN 2277-188 Journal
More informationA surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant
A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant Ashley J. James Department of Aerospace Engineering and Mechanics, University of Minnesota John Lowengrub
More informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More informationME615 Project Presentation Aeroacoustic Simulations using Lattice Boltzmann Method
ME615 Project Presentation Aeroacoustic Simulations using Lattice Boltzmann Method Kameswararao Anupindi Graduate Research Assistant School of Mechanical Engineering Purdue Universit December 11, 11 Outline...
More informationMagnetic Field and Chemical Reaction Effects on Convective Flow of
Communications in Applied Sciences ISSN 221-7372 Volume 1, Number 1, 213, 161-187 Magnetic Field and Chemical Reaction Effects on Convective Flow of Dust Viscous Fluid P. Mohan Krishna 1, Dr.V.Sugunamma
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationDepartment of Mathematics, Usmanu Danfodiyo University, Sokoto, Nigeria. DOI: / Page
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. II (Jan - Feb. 5), PP 4-58 www.iosrjournals.org Finite difference solutions of magneto hdrodnamic free convective
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationChemical Reaction, Radiation and Dufour Effects on Casson Magneto Hydro Dynamics Fluid Flow over A Vertical Plate with Heat Source/Sink
Global Journal of Pure and Applied Mathematics. ISSN 973-768 Volume, Number (6), pp. 9- Research India Publications http://www.ripublication.com Chemical Reaction, Radiation and Dufour Effects on Casson
More informationFigure 11.1: A fluid jet extruded where we define the dimensionless groups
11. Fluid Jets 11.1 The shape of a falling fluid jet Consider a circular orifice of a radius a ejecting a flux Q of fluid density ρ and kinematic viscosity ν (see Fig. 11.1). The resulting jet accelerates
More informationNumerical Investigation of Cavitation Bubble Collapsing Behavior
First International Smposium on Marine Propulsors SMP 09, Trondheim, Norwa, June 009 Numerical Investigation of Cavitation Bubble Collapsing Behavior Beong Rog Shin 1 1 Deptartment of Mechanical Engineering,
More informationCHAPTER-III CONVECTION IN A POROUS MEDIUM WITH EFFECT OF MAGNETIC FIELD, VARIABLE FLUID PROPERTIES AND VARYING WALL TEMPERATURE
CHAPER-III CONVECION IN A POROUS MEDIUM WIH EFFEC OF MAGNEIC FIELD, VARIABLE FLUID PROPERIES AND VARYING WALL EMPERAURE 3.1. INRODUCION Heat transer studies in porous media ind applications in several
More informationVERTICAL BRIDGMAN CONFIGURATION HEATED FROM BELOW: 3D BIFURCATION AND STABILITY ANALYSIS
Third International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia -2 December 23 VERTICAL BRIDGMAN CONFIGURATION HEATED FROM BELOW: 3D BIFURCATION AND STABILITY ANALYSIS
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More informationLift Enhancement on Unconventional Airfoils
Lift Enhancement on nconventional Airfoils W.W.H. Yeung School of Mechanical and Aerospace Engineering Nanang Technological niversit, North Spine (N3), Level 2 50 Nanang Avenue, Singapore 639798 mwheung@ntu.edu.sg
More informationRANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM 1
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM Jirka Roubal Vladimír Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Karlovo náměstí
More informationRadiation and Soret Effects to MHD Flow in Vertical Surface with Chemical reaction and Heat generation through a Porous Medium
ISSN (e): 5 35 Vol, 4 Issue, 7 Jul 4 International Journal of Computational Engineering Research (IJCER) Radiation and Soret Effects to MHD Flow in Vertical Surface with Chemical reaction and Heat generation
More informationRESEARCH ARTICLE. Zonal PANS: evaluation of different treatments of the RANS-LES interface
March 3, 26 Journal of Turbulence paper Journal of Turbulence Vol., No., 2, 32 Volume 7, Issue 3, pp. 274-37, 26 RESEARCH ARTICLE 2 3 4 6 7 8 9 Zonal PANS: evaluation of different treatments of the RANS-LES
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationBuoyancy Effects on Unsteady MHD Flow of a Reactive Third-Grade Fluid with Asymmetric Convective Cooling
Journal of Applied Fluid Mechanics, Vol. 8, No. 4, pp. 931-941, 215. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. DOI: 1.18869/acadpub.jafm.73.238.22865 Buoanc Effects on Unstead
More informationNumerical investigation of cavitation-regimes in a converging-diverging nozzle
Numerical investigation of cavitation-regimes in a converging-diverging nozzle 1 Polina Gorkh, 1 Steffen J. Schmidt, and 1 Nikolaus A. Adams 1 Institute of Aerodynamics and Fluid Mechanics, Department
More informationEFFECTS OF THERMAL RADIATION AND HEAT TRANSFER OVER AN UNSTEADY STRETCHING SURFACE EMBEDDED IN A POROUS MEDIUM IN THE PRESENCE OF HEAT SOURCE OR SINK
THERMAL SCIENCE, Year 011, Vol. 15, No., pp. 477-485 477 EFFECTS OF THERMAL RADIATION AND HEAT TRANSFER OVER AN UNSTEADY STRETCHING SURFACE EMBEDDED IN A POROUS MEDIUM IN THE PRESENCE OF HEAT SOURCE OR
More informationNUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A
More informationPlanar Geometry Ferrofluid Flows in Spatially Uniform Sinusoidally Time-varying Magnetic Fields
Presented at the 11 COMSOL Conference in Boston Planar Geometry Ferrofluid Flows in Spatially Uniform Sinusoidally Time-varying Magnetic Fields Shahriar Khushrushahi, Alexander Weddemann, Young Sun Kim
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationNumerical Solution of Non-Darcian Effects on Natural Convection in a Rectangular Porous Enclosure with Heated Walls
International Journal of Advanced Mechanical Engineering. ISSN 2250-3234 Volume 8, Number 1 (18), pp. 71-86 Research India Publications http://www.ripublication.com Numerical Solution of Non-Darcian Effects
More informationNUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR
NUMERICAL COMPUTATION OF PARTICLE TRANSPORT IN TURBULENT FLOW FIELD IN ELECTROSTATIC PRECIPITATOR Jenő SUDA, István BERTA, Gergel KRISTÓF +, Tamás LAJOS # Ph.D. Student, + Assistant, # Professor, Department
More informationHeat and Mass Transfer Flow of MHD Viscous Dissipative Fluid in a Channel with a Stretching and Porous Plate
Appl. Math. Inf. Sci. Lett. 5, No. 3, 8-87 (27) 8 Applied Mathematics & Information Sciences Letters An International Journal http://dx.doi.org/.8576/amisl/53 Heat and Mass Transfer Flow of MHD Viscous
More informationVALIDATION OF CALCULATIONS OF THE VISCOUS FLOW AROUND A SHIP IN OBLIQUE MOTION
The First MARIN-NMRI Workshop in Toko (Oktober 25-26, 24) VALIDATION OF CALCULATIONS OF THE VISCOUS FLOW AROUND A SHIP IN OBLIQUE MOTION Serge TOXOPEUS Maritime Research Institute Netherlands (MARIN),
More informationFinite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport
Finite Element Techniques for the Numerical Simulation of Two-Phase Flows with Mass Transport Christoph Lehrenfeld and Arnold Reusken Preprint No. 413 December 2014 Key words: Two-phase flow, mass transport,
More informationarxiv: v1 [physics.flu-dyn] 30 Oct 2017
Passive scalar transport in rotating turbulent channel flow Geert Brethouwer 1 1 Linné FLOW Centre, KTH Mechanics, SE-144 Stockholm, Sweden (Dated: October 31, 217) arxiv:171.1931v1 [phsics.flu-dn] 3 Oct
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationL O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S
L O N G I T U D I N A L V O R T I C E S I N G R A N U L A R F L O W S Yoël FORTERRE, Olivier POULIQUEN Institut Universitaire des Sstèmes Thermiques et Industriels (UMR 6595 CNRS) 5, rue Enrico Fermi 13453
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationSupplementary Information: Manipulating Acoustic Wavefront by Inhomogeneous Impedance and Steerable Extraordinary Reflection
Supplementar Information: Manipulating Acoustic Wavefront b Inhomogeneous Impedance and Steerable Extraordinar Reflection Jiajun Zhao 1,, Baowen Li,3, Zhining Chen 1, and Cheng-Wei Qiu 1, 1 Department
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationSimulation of Electric Fields in Small Size Divertor Tokamak Plasma Edge
Energ and Power Engineering, 21, 39-45 doi:1.4236/epe.21.217 Published Online Februar 21 (http://www.scirp.org/journal/epe) Simulation of Electric Fields in Small Size Divertor Tokamak Plasma Edge Plasma
More informationARTICLE IN PRESS. Numerical investigation of the effect of insoluble surfactants on drop deformation and breakup in simple shear flow
S0021-9797(05)01265-8/FLA AID:11865 Vol. ( ) [+model] P.1 (1-26) YJCIS:m5+ v 1.50 Prn:10/01/2006; 14:37 yjcis11865 by:vik p. 1 Journal of Colloid and Interface Science ( ) www.elsevier.com/locate/jcis
More informationGrid-independent large-eddy simulation of compressible turbulent flows using explicit filtering
Center for Turbulence Researc Proceedings of te Summer Program 2010 203 Grid-independent large-edd simulation of compressible turbulent flows using explicit filtering B D. You, S. T. Bose AND P. Moin Te
More informationTransient Heat Transfer Experiment. ME 331 Introduction to Heat Transfer. June 1 st, 2017
Transient Heat Transfer Experiment ME 331 Introduction to Heat Transfer June 1 st, 2017 Abstract The lumped capacitance assumption for transient conduction was tested for three heated spheres; a gold plated
More informationThe Gas-assisted Expelled Fluid Flow in the Front of a Long Bubble in a Channel
C. H. Hsu, P. C. Chen,. Y. ung, G. C. uo The Gas-assisted Epelled Fluid Flow in the Front of a Long Bubble in a Channel C.H. Hsu 1, P.C. Chen 3,.Y. ung, and G.C. uo 1 1 Department of Mechanical Engineering
More informationFLOW ASSURANCE: DROP COALESCENCE IN THE PRESENCE OF SURFACTANTS
FLOW ASSURANCE: DROP COALESCENCE IN THE PRESENCE OF SURFACTANTS Vishrut Garg and Osman A. Basaran Davidson School of Chemical Engineering Purdue University With special thanks to: Krish Sambath (now at
More informationIntroduction to Heat and Mass Transfer. Week 12
Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationPhysics of Tsunamis. This is an article from my home page: Ole Witt-Hansen 2011 (2016)
Phsics of Tsunamis This is an article from m home page: www.olewitthansen.dk Ole Witt-Hansen 11 (16) Contents 1. Constructing a differential equation for the Tsunami wave...1. The velocit of propagation
More information