Parallel Adaptive Finite Element Methods for Problems in Natural Convection
|
|
- Megan Mitchell
- 5 years ago
- Views:
Transcription
1 Parallel Adaptive Finite Element Methods for Problems in Natural Convection John W. Peterson CFDLab Univ. of Texas at Austin April 18, 2008
2 Outline 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
3 Physical Problem 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
4 Physical Problem Buoyancy basic concepts Buoyancy: Warmer fluids are slightly less dense than cold fluids. g Cold Plate Slightly Cooler Fluid Slightly Warmer Fluid Hot Plate
5 Physical Problem Buoyancy basic concepts Rayleigh s Number: a non-dimensional parameter which tells the relative importance of buoyant forces to dissipative effects Ra = Buoyant Force Dissipative Effects Rayleigh s famous calculation predicted Ra crit = 1708 for onset of fluid motion in an infinite horizontal layer.
6 Physical Problem Surface tension basic concepts Fluid surface tension decreases with temperature g Slightly Cooler Fluid Slightly Warmer Fluid
7 Physical Problem Surface tension basic concepts The same dissipative effects (viscosity, diffusion) retard fluid motion on the surface Marangoni s Number: a non-dimensional ratio of surface tension strength to dissipative effects Ma = Surface Tension Dissipative Effects
8 Physical Problem Aspect ratio Free T. Surface n=bit d L T. n=0 Side Walls T=Hot Bottom Wall Aspect Ratio Γ := L d
9 Parallel FEM for Natural Convection Physical Problem Hexagonal Cells In large (Γ 1) aspect ratio fluid layers, hexagonal convection cells form Be nard, 1900 Velarde et. al (Album of Fluid Motion)
10 Physical Problem Triangular Container Medale & Cerisier, 2002 Peterson, 2004 At moderate Ma, presence of hexagonal cells is independent of container shape Γ := A/d = 22
11 Physical Problem Convection far from onset Square Cells Γ At (c), Ma 781! Exhibits hysteresis Schatz et. al, 1999
12 Physical Problem Small aspect ratio domains In smaller aspect ratio domains (Γ 10), hexagons do not form The number of convection cells depends on the aspect ratio and other parameters E. L. Koschmieder and S. A. Prahl, Surface-tension-driven Bénard convection in small containers, J. Fluid Mech
13 Physical Problem Small aspect ratio domains Γ = 6.18, Ma = 96 Γ = 6.36, Ma = 93 Koschmieder, 1990
14 Goals and Contributions 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
15 Goals and Contributions Goals Main Goal: Improve understanding of RBM convection in small aspect ratio containers. Hypothesis: In a square cross-section container of small aspect ratio Γ, there exist steady state cellular convection patterns which depend only on Γ.
16 Goals and Contributions Contributions Developed scaling (non-dimensionalization) for the RBM convection equations in which the aspect ratio, Γ, appears as an independent, continuous parameter. Developed a pseudo-arclength continuation algorithm for the RBM convection equations. Implemented generic continuation solver algorithm in LibMesh. Expanded catalog of spontaneous spatio-temporal pattern formation in RBM convection with high-fidelity numerical simulations.
17 Goals and Contributions Contributions Explicit second-order AB predictor implemented with continuation scheme. Performed linear stability analysis of solution branches; linked application code to state-of-the-art eigenanalaysis software package SLEPc. Performed exploratory AMR studies of combined Navier-slip thermocapillary convection. Several contributions to the LibMesh open-source library, and finite element analysis scientific community.
18 Governing Equations 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
19 Governing Equations Standard Non-Dimensional Equations ( ) 1 u + (u )u σ Pr t = Ra Tĝ u = 0 T t + u T + u T 0 = 2 T
20 Governing Equations Boundary Conditions Bottom: u = 0 T = 0
21 Governing Equations Boundary Conditions Bottom: u = 0 T = 0 Sides: u = 0 T n = 0
22 Governing Equations Boundary Conditions Bottom: u = 0 T = 0 Sides: u = 0 T n = 0 Free surface: σ n = Ma T T n = Bi T
23 Governing Equations Variational Statement: Conservation of Momentum & Mass Ω Find u, p, T such that: [ 1 ( u Pr t ) + (u )u Ra Ω Ω ] v + ε(u): v p ( v) Tĝ v dx + Ma ( u)q dx = 0 dx + Ω s ( T v ) ds = 0
24 Governing Equations Variational Statement: Conservation of Energy Ω ( ) T t + u T + u T 0 w + ( T w) dx + Bi Tw ds = 0 Ω s hold v, q, w V Q W.
25 Governing Equations Numerical Method From the variational statement, we construct an LBB-stable Galerkin finite element discretization in space.
26 Governing Equations Numerical Method From the variational statement, we construct an LBB-stable Galerkin finite element discretization in space. We then perform a standard θ-method discretization in time (plus adaptive step-doubling timestep selection).
27 Governing Equations Numerical Method From the variational statement, we construct an LBB-stable Galerkin finite element discretization in space. We then perform a standard θ-method discretization in time (plus adaptive step-doubling timestep selection). This leads to a discrete system of nonlinear equations which is linearized and solved via Newton s method.
28 Governing Equations Library Description Tumor Angiogenesis Compressible RBM NS Petsc Laspack DD LibMesh STL Metis Library Structure Basic libraries are LibMesh s roots Application branches built off the library trunk MPI
29 Unsteady Simulations 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
30 Unsteady Simulations Simulation Parameters Ma = 92 Ra = 30 (Ma/Ra 3.07, surf. ten. effects dominate) Pr = 880 (corresponds to a viscous silicon oil) Bi = 0.2 Values representative of those used by Koschmieder (1990) and Medale and Cerisier (2002) Direct quantitative comparisons with experiment not generally possible due to experimental uncertainty in material properties
31 Unsteady Simulations Adaptive Timestepping Adaptive step-doubling method with one-step θ scheme as base time discretization Steady state detection when the quantity u t un+1 u n t falls to approximately machine precision (10 13 to ) at t max. t max / t min = 10 2
32 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3
33 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3 1
34 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3 2a
35 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3 2b
36 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3 2c
37 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3 2d
38 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3 2e
39 Unsteady Simulations Representative Results, Γ = log t / t 0 log( U / t) Timestep 3 3
40 Unsteady Simulations Representative Results, Γ = 6 Movie: Time Evolution for Γ = 6
41 Unsteady Simulations Summary of Unsteady Results After collecting steady-state converged data for a range of aspect ratios, we plot it versus Nusselt number: Nu := ds Ω s T
42 Unsteady Simulations Summary of Unsteady Results After collecting steady-state converged data for a range of aspect ratios, we plot it versus Nusselt number: Nu := ds Ω s T The collected unsteady solution data does not provide the full story of the solution dependence on Γ.
43 Unsteady Simulations Summary of Unsteady Results After collecting steady-state converged data for a range of aspect ratios, we plot it versus Nusselt number: Nu := ds Ω s T The collected unsteady solution data does not provide the full story of the solution dependence on Γ. The branches of steady solutions appear to have discontinuous jumps near certain values of Γ.
44 Unsteady Simulations Summary of Unsteady Results Nu Cell 2 Cell 3 Cell 4 Cell Cell 6 Cell 7 Cell 8 Cell Aspect Ratio, Γ
45 Continuation Results 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
46 Continuation Results Variation of Aspect Ratio d Γ=1 g L Γ=2 Γ=3 Aspect Ratio Γ := L/d We consider containers of constant d and increasing L
47 Continuation Results Scaled Equations Cartesian (x 1, x 2, x 3 ) coordinate system, gravity vector in the x 3 direction. Non-dimensionalize with respect to L for horizontal scales, d for vertical scales x 1 := x 1 L x 2 := x 2 L The gradient operator is transformed as 1 (,, Γ ) L x1 x2 x3 x 3 := x 3 d := 1 L Γ
48 Continuation Results Scaled Equations ( ) 1 u Pr t + (u Γ)u Γ σ Γ = Γ 3 Ra Tĝ Γ u = 0 T t + u Γ T + u Γ T 0 = Γ T 2
49 Continuation Results Variational Statement: Conservation of Momentum & Mass Ω Find u, p, T such that: [ ( ) 1 u Pr t + (u Γ)u Γ 3 Ra Ω Ω ] v + ε Γ (u): Γ v p ( Γ v) Tĝ v dx + Γ 2 Ma ( Γ u)q dx = 0 dx + Ω s ( T v ) ds = 0
50 Continuation Results Variational Statement: Conservation of Energy Ω ( ) T t + u Γ T + u Γ T 0 w + ( Γ T Γ w) dx + Γ 2 Bi Tw ds = 0 Ω s hold v, q, w V Q W.
51 Continuation Results Continuation Basics We write the (discretized) steady form of the variational statement as G(u, λ) = 0 The Jacobian G u may be singular depending on u, λ A standard continuation algorithm will fail near such points (Euler-Newton) Generic continuation in arbitrary parameters now available in LibMesh
52 Continuation Results Continuation Basics The pseudo-arclength continuation technique augments G with a constraint N(u, λ, s) = 0 s is the arclength parameter, s > 0 Must solve augmented Jacobian systems [ Gu G λ N t u N λ ] k [ ] k+1 δu = δλ [ ] k G N The augmented matrix is of full rank even when G u is singular
53 Continuation Results Continuation Algorithm Obtain u 0, λ 0 from unsteady solve Choose next parameter value: λ 1 = λ 0 + λ ( λ small) Solve steady equations for u 1, use u 0 as initial guess Compute tangents u/ s 1, λ/ s 1 via finite difference Compute initial guess ũ 2, λ 2 using the tangent for n = 2 to n steps do Solve augmented Jacobian system for new values u n, λ n Update tangents u/ s n, λ/ s n with add l linear solve Compute initial guess ũ n+1, λ n+1 end for
54 Continuation Results Continuation Algorithm Prediction & Correction Steps u λ
55 Continuation Results Continuation Algorithm Prediction & Correction Steps u λ
56 Continuation Results Continuation Algorithm Prediction & Correction Steps u λ
57 Continuation Results Continuation Algorithm Prediction & Correction Steps u λ
58 Continuation Results Continuation Algorithm Prediction & Correction Steps u λ
59 Continuation Results Arclength Stepsize Selection Nu Proper adaptive stepsize selection is required to prevent stepping past important branch features. Ma
60 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 1
61 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 2
62 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 3
63 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 4
64 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 5
65 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 6
66 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 7
67 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 8
68 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 9
69 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 10
70 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ 11
71 Continuation Results Representative Solution Branch Nu Aspect Ratio, Γ Koschmieder, 1990
72 Continuation Results Summary of Continuation Results Nu Aspect Ratio, Γ
73 Parallel Adaptive Simulations 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
74 Parallel Adaptive Simulations Motivation Interested in controlling observed natural convection flow patterns and structure with non-invasive techniques such as surface treatments Shear-dependent slip, aka Navier slip, is a model for investigating these flow control ideas qualitatively and quantitatively
75 Parallel Adaptive Simulations Navier Slip BC Proposed by C. L. M. H. Navier, 1823 and independently by Maxwell, 1879 Amount of slip Rate of shear
76 Parallel Adaptive Simulations Navier Slip BC Proposed by C. L. M. H. Navier, 1823 and independently by Maxwell, 1879 Amount of slip Rate of shear Constant of proportionality is slip length, L
77 Parallel Adaptive Simulations Navier Slip BC Proposed by C. L. M. H. Navier, 1823 and independently by Maxwell, 1879 Amount of slip Rate of shear Constant of proportionality is slip length, L In the variational statement, we replace surface traction term (σ n) v = L 1 (u u wall ) v
78 Parallel Adaptive Simulations Navier Slip BC Proposed by C. L. M. H. Navier, 1823 and independently by Maxwell, 1879 Amount of slip Rate of shear Constant of proportionality is slip length, L In the variational statement, we replace surface traction term (σ n) v = L 1 (u u wall ) v Wall-normal and tangential slip lengths may be chosen independently
79 Parallel Adaptive Simulations Navier Slip Applications Moving contact line problem Flow near sharp corners/in small capillaries Electro-osmotically driven flow Hydrophobic channel flow Micro/nanoscale flows Surface polishing in semiconductor manufacturing
80 Parallel Adaptive Simulations Navier-slip in RBM Convection Performed a series of AMR simulations with unit Navier slip allowed on different quadrants of the container bottom.
81 Parallel Adaptive Simulations Navier-slip Quadrants and Reference Solution II III I IV Slip Quadrants Reference No-slip Solution
82 Parallel Adaptive Simulations Navier-slip Results Quadrant IV Quadrants I and III
83 Parallel Adaptive Simulations Navier-slip Results Quadrant III and IV Quadrants I, III, and IV
84 Parallel Adaptive Simulations Navier-slip Adapted Grids Temperature Contours Vel. Magnitude Contours
85 Parallel Adaptive Simulations Line Plots Through Slip Discontinuity x20x2 40x40x4 80x80x8 AMR(3) -4 u x 2
86 Parallel Adaptive Simulations Navier-slip Time Evolution Movies Movies: Top Surface Bottom Surface
87 Conclusion and Future Work 1 Physical Problem 2 Goals and Contributions 3 Governing Equations 4 Unsteady Simulations 5 Continuation Results 6 Parallel Adaptive Simulations 7 Conclusion and Future Work
88 Conclusion and Future Work Contributions Developed scaling (non-dimensionalization) for the RBM convection equations in which the aspect ratio, Γ, appears as an independent, continuous parameter. Developed a pseudo-arclength continuation algorithm for the RBM convection equations. Implemented generic continuation solver algorithm in LibMesh. Expanded catalog of spontaneous spatio-temporal pattern formation in RBM convection with high-fidelity numerical simulations.
89 Conclusion and Future Work Contributions Explicit second-order AB predictor implemented with continuation scheme. Performed linear stability analysis of solution branches; linked application code to state-of-the-art eigenanalaysis software package SLEPc. Performed exploratory AMR studies of combined Navier-slip thermocapillary convection. Several contributions to the LibMesh open-source library, and finite element analysis scientific community.
90 Conclusion and Future Work Future Work Additional parametric studies (Ma, Bi) including adaptive mesh refinement along solution branches Non-square (rectangular, circular, triangular... ) containers of small Γ Bifurcation tracking (via extended systems) Deforming free surface simulations Goal-oriented error indicators for driving AMR simulations to compute quantities of interest (Nu) Optimal shape design of containers to maximize heat transfer efficiency
91 Reference Slides Reference Material (Extra Slides)
92 Reference Slides Temperature Decomposition T T 0 T ~
93 Reference Slides Additional Remarks on the Scaling In the scaled variational statement, we obtained terms proportional to: Γ 3 Ra Γ 2 Ma Where Ra, Ma are defined in terms of the layer depth in the usual way Ra := ρ 0 β g T sd 3 αµ Ma := γ T T sd αµ
94 Reference Slides Additional Remarks on the Scaling In the scaled equations, the Γ 3 Ra term is constant with respect to d while the Γ 2 Ma term scales as d 1. Thus, the Ma term dominates for small d while the Ra term dominates for large d in both cases Ra term Ma term Classical d 3 d Scaled const d 1
95 Reference Slides Additional Remarks on the Scaling Alternatively, we can define analogous parameters Ra L := ρ 0 β g T sl 3 αµ Ma L := γ T T sl αµ in terms of the horizontal length scale L
96 Reference Slides Additional Remarks on the Scaling We then have Γ 3 Ra = Ra L Γ 2 Ma = = ( ) 2 L γ Ts d T d αµ ( ) L γ T T s L d αµ Let L be fixed = ΓMa L For small aspect ratio (Γ 0) layers, buoyant effects dominate Vice-versa for large aspect ratio layers
97 Reference Slides Surface Deflection General normal traction condition for curved surface σ n n = Kγ Kγ 0 Mean Curvature ( 1 K := + 1 ) R 1 R 2 Surface orientation n, t 1, t 2 is unknown; depends on solution {u, p, T}
98 Reference Slides Surface Deflection Surface deflections may be considered small if Kd Ca 1 Capillary number In the present work Ca := µα γ 0 d Ca
99 Reference Slides Surface Deflection Rahal et al. have recently (2007) employed interferometry to measure surface deflections The examined a related but slightly thicker (5 mm) fluid layer in a circular container of Γ := R/d = 6. Buoyant effects dominant (Ra/Ma 9) but thermocapillary flow still present. Rahal et al. observed convex cells (raised upwelling regions) on the order of a few ( 3) microns in height over distances of O(10) mm
100 Reference Slides Pseudo -Arclength Constraint We employ a linearized or pseudo arclength constraint as in Keller 1977, the LOCA library, and others N α (u, λ, s) := αi 2 (u u(s i )) t u s +(λ λ(s i )) λ si s (s s i ) si
101 Reference Slides Pseudo -Arclength Constraint We employ a linearized or pseudo arclength constraint as in Keller 1977, the LOCA library, and others N α (u, λ, s) := αi 2 (u u(s i )) t u s +(λ λ(s i )) λ si s (s s i ) si Solution and parameter tangents u si, λ s s si
102 Reference Slides Pseudo -Arclength Constraint We employ a linearized or pseudo arclength constraint as in Keller 1977, the LOCA library, and others N α (u, λ, s) := αi 2 (u u(s i )) t u s +(λ λ(s i )) λ si s (s s i ) si Solution and parameter tangents u si, λ s s si Scaling parameter α i ensures similar order-of-magnitude contributions from each constraint term
103 Reference Slides Schemes for Adaptively Selecting Arclength Stepsizes Ratio of successive tangents u λ si 1 s i+1 = s i u λ si
104 Reference Slides Schemes for Adaptively Selecting Arclength Stepsizes Most recent number of Newton iterations required [ ( ) ] 2 Nmax N i s i+1 = 1 + a s i N max 1
105 Reference Slides Generalized Eigenproblem The spatially-discretized equations can be written B u t = G(u, λ)
106 Reference Slides Generalized Eigenproblem The spatially-discretized equations can be written B u = G(u, λ) t B is a generalized mass matrix (linear operator)
107 Reference Slides Generalized Eigenproblem The spatially-discretized equations can be written B u = G(u, λ) t B is a generalized mass matrix (linear operator) u is a vector of FE coefficients
108 Reference Slides Generalized Eigenproblem The spatially-discretized equations can be written B u = G(u, λ) t B is a generalized mass matrix (linear operator) u is a vector of FE coefficients λ is a generalized parameter (Γ, Ma,... )
109 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ)
110 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ) Let η := u u 0
111 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ) Let η := u u 0 Applying B to η we obtain B η = B u Bu 0
112 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ) Let η := u u 0 Applying B to η we obtain B η = B u Bu 0 = G(u, λ) G(u 0, λ)
113 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ) Let η := u u 0 Applying B to η we obtain B η = B u Bu 0 = G(u, λ) G(u 0, λ) = G(u, λ)
114 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ) Let η := u u 0 Applying B to η we obtain B η = B u Bu 0 = G(u, λ) G(u 0, λ) = G(u, λ) = G(u 0 + η, λ)
115 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ) Let η := u u 0 Applying B to η we obtain B η = B u Bu 0 = G(u, λ) G(u 0, λ) = G(u, λ) = G(u 0 + η, λ) = G(u 0, λ) + G u (u 0, λ)η + O( η 2 )
116 Reference Slides Generalized Eigenproblem Let u 0 be given such that 0 = G(u 0, λ) Let η := u u 0 Applying B to η we obtain B η = B u Bu 0 = G(u, λ) G(u 0, λ) = G(u, λ) = G(u 0 + η, λ) = G(u 0, λ) + G u (u 0, λ)η + O( η 2 ) G u (u 0, λ)η
117 Reference Slides Generalized Eigenproblem Letting A := G u (u 0, λ)
118 Reference Slides Generalized Eigenproblem Letting A := G u (u 0, λ) And inserting the ansatz η = exp(σt)x yields B η = Aη
119 Reference Slides Generalized Eigenproblem Letting A := G u (u 0, λ) And inserting the ansatz η = exp(σt)x yields B η = Aη σb exp(σt)x = A exp(σt)x
120 Reference Slides Generalized Eigenproblem Letting A := G u (u 0, λ) And inserting the ansatz η = exp(σt)x yields B η = Aη σb exp(σt)x = A exp(σt)x σbx = Ax
121 Reference Slides Generalized Eigenproblem Letting A := G u (u 0, λ) And inserting the ansatz η = exp(σt)x yields B η = Aη σb exp(σt)x = A exp(σt)x σbx = Ax A generalized eigenproblem for σ, x
122 Reference Slides Eigensolver Details Direct calculation of eigenvalues is O(n 3 ) To determine stability, we only need the sign of the rightmost eigenvalue The Arnoldi iteration is a relatively efficient technique for computing eigenvalues largest in magnitude (not necessarily rightmost)
123 Reference Slides Eigensolver Details Arnoldi is related to the power iteration Computes eigenvalues of the orthogonal projection of A onto the Krylov subspace Already implemented in the SLEPc library, a collection of eigensolvers built on top of PETSc
124 Reference Slides Spectral Transformations Arnoldi converges to the eigenvalues largest in magnitude Our generalized mass matrix B has time-independent constraint rows which effectively lead to eigenvalues We can use a spectral transformation such as shift-and-invert to avoid converging to large negative eigenvalues.
125 Reference Slides Spectral Transformations σbx = Ax Shift-and-invert transform: start with generalized eigenproblem
126 Reference Slides Spectral Transformations σbx = Ax (σ s) Bx = (A sb) x Choose real-valued shift s σ, subtract sbx from both sides
127 Reference Slides Spectral Transformations σbx = Ax (σ s) Bx = (A sb) x (σ s) (A sb) 1 Bx = (A sb) 1 (A sb) x Left multiply both sides by (A sb) 1
128 Reference Slides Spectral Transformations σbx = Ax (σ s) Bx = (A sb) x (σ s) (A sb) 1 Bx = (A sb) 1 (A sb) x (σ s) (A sb) 1 Bx = x Left multiply both sides by (A sb) 1
129 Reference Slides Spectral Transformations σbx = Ax (σ s) Bx = (A sb) x (σ s) (A sb) 1 Bx = (A sb) 1 (A sb) x (σ s) (A sb) 1 Bx = x (A sb) 1 Bx = (σ s) 1 x Multiply both sides by (σ s) 1
130 Reference Slides Spectral Transformations σbx = Ax (σ s) Bx = (A sb) x (σ s) (A sb) 1 Bx = (A sb) 1 (A sb) x (σ s) (A sb) 1 Bx = x (A sb) 1 Bx = (σ s) 1 x Tx = γx We are left with a standard eigenproblem, where T := (A sb) 1 B, γ := (σ s) 1
131 Reference Slides Spectral Transformations Shift-and-invert maps eigenvalues far from σ to zero (small magnitude) Arnoldi does not converge to these small-magnitude eigenvalues Can also be used to converge quickly to a suspected eigenvalue σ s
132 Reference Slides Solution of Bordered Linear Systems The continuation algorithm requires solution of bordered linear systems [ ] [ ] [ ] Gu G λ δu r = δλ ρ N t u N λ
133 Reference Slides Solution of Bordered Linear Systems The continuation algorithm requires solution of bordered linear systems [ ] [ ] [ ] Gu G λ δu r = δλ ρ N t u N λ δu, δλ Newton updates or new tangent directions
134 Reference Slides Solution of Bordered Linear Systems The continuation algorithm requires solution of bordered linear systems [ ] [ ] [ ] Gu G λ δu r = δλ ρ N t u N λ δu, δλ Newton updates or new tangent directions Right-hand sides: [ ] r = ρ [ ] G N or [ 0 N s ]
135 Reference Slides Solution of Bordered Linear Systems Can solve the augmented system efficiently via iterative methods e.g. Barragy and Carey, 1988 Here we employ the decomposition approach described by Keller, 1977
136 Reference Slides Solution of Bordered Linear Systems Can solve the augmented system efficiently via iterative methods e.g. Barragy and Carey, 1988 (+) Single solve step Here we employ the decomposition approach described by Keller, 1977 ( ) Multi-solve step required
137 Reference Slides Solution of Bordered Linear Systems Can solve the augmented system efficiently via iterative methods e.g. Barragy and Carey, 1988 (+) Single solve step (+) Augmented matrix has full rank Here we employ the decomposition approach described by Keller, 1977 ( ) Multi-solve step required ( ) Solves may involve nearly singular G u
138 Reference Slides Solution of Bordered Linear Systems Can solve the augmented system efficiently via iterative methods e.g. Barragy and Carey, 1988 (+) Single solve step (+) Augmented matrix has full rank ( ) Changes the system matrix sparsity pattern Here we employ the decomposition approach described by Keller, 1977 ( ) Multi-solve step required ( ) Solves may involve nearly singular G u (+) System matrix sparsity unchanged
139 Reference Slides Solution of Bordered Linear Systems Can solve the augmented system efficiently via iterative methods e.g. Barragy and Carey, 1988 (+) Single solve step (+) Augmented matrix has full rank ( ) Changes the system matrix sparsity pattern ( ) Requires modification of existing code Here we employ the decomposition approach described by Keller, 1977 ( ) Multi-solve step required ( ) Solves may involve nearly singular G u (+) System matrix sparsity unchanged (+) Requires minimal modifications to existing code
140 Reference Slides Solution of Bordered Linear Systems Solution Algorithm [ Gu G λ N t u N λ ] [ ] δu δλ 1 Declare auxiliary vectors y, z = [ ] r ρ
141 Reference Slides Solution of Bordered Linear Systems Solution Algorithm [ Gu G λ N t u N λ ] [ ] δu δλ 1 Declare auxiliary vectors y, z 2 Solve #1: G u y = G λ for y = [ ] r ρ
142 Reference Slides Solution of Bordered Linear Systems Solution Algorithm [ Gu G λ N t u N λ ] [ ] δu δλ 1 Declare auxiliary vectors y, z 2 Solve #1: G u y = G λ for y 3 Solve #2: G u z = r for z = [ ] r ρ
143 Reference Slides Solution of Bordered Linear Systems Solution Algorithm [ Gu G λ N t u N λ ] [ ] δu δλ 1 Declare auxiliary vectors y, z 2 Solve #1: G u y = G λ for y 3 Solve #2: G u z = r for z 4 Set: = [ ] r ρ δλ = ρ Nt uz N λ N t uy
144 Reference Slides Solution of Bordered Linear Systems Solution Algorithm [ Gu G λ N t u N λ ] [ ] δu δλ 1 Declare auxiliary vectors y, z 2 Solve #1: G u y = G λ for y 3 Solve #2: G u z = r for z 4 Set: 5 Set: = [ ] r ρ δλ = ρ Nt uz N λ N t uy δu = z (δλ)y
145 Reference Slides The Stick-Navier Slip Singularity log 10 u Slope=1/ x20x2 40x40x x80x8 AMR(3) AMR(4) log 10 r
LibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu and Roy H. Stogner roystgnr@cfdlab.ae.utexas.edu Univ. of Texas at Austin September 9, 2008 1 Introduction 2 Weighted Residuals
More informationLibMesh Experience and Usage
LibMesh Experience and Usage John W. Peterson peterson@cfdlab.ae.utexas.edu Univ. of Texas at Austin January 12, 2007 1 Introduction 2 Weighted Residuals 3 Poisson Equation 4 Other Examples 5 Essential
More informationFinding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems
Finding Rightmost Eigenvalues of Large, Sparse, Nonsymmetric Parameterized Eigenvalue Problems AMSC 663-664 Final Report Minghao Wu AMSC Program mwu@math.umd.edu Dr. Howard Elman Department of Computer
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationMYcsvtu Notes HEAT TRANSFER BY CONVECTION
www.mycsvtunotes.in HEAT TRANSFER BY CONVECTION CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationNumerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 541 552 Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus M. Sankar
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationTHERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW
THERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW Luigi De Giorgi, Carlo Cima, Emilio Cafaro Dipartimento di Energetica, Politecnico di Torino, Torino, Italy Volfango Bertola School
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationAnalysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle
Nonlinear Analysis: Modelling and Control, 2011, Vol. 16, No. 1, 89 99 89 Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle S. Parvin,
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationNATURAL CONVECTION OF AIR IN TILTED SQUARE CAVITIES WITH DIFFERENTIALLY HEATED OPPOSITE WALLS
Proceedings of the International onference on Mechanical Engineering 0 (IME0 8-0 December 0, Dhaka, Bangladesh IME- NATURAL ONVETION OF AIR IN TILTED SQUARE AVITIES WIT DIFFERENTIALLY EATED OPPOSITE WALLS
More informationHEAT TRANSFER BY CONVECTION. Dr. Şaziye Balku 1
HEAT TRANSFER BY CONVECTION Dr. Şaziye Balku 1 CONDUCTION Mechanism of heat transfer through a solid or fluid in the absence any fluid motion. CONVECTION Mechanism of heat transfer through a fluid in the
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 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 informationTowards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below
Copyright 2008 Tech Science Press FDMP, vol.4, no.4, pp.263-269, 2008 Towards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below G. Accary 1, S. Meradji
More informationEffects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationNATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A FLUSH MOUNTED HEATER ON A SIDE WALL
Journal of Naval Architecture and Marine Engineering December, 2010 DOI: 10.3329/jname.v7i2.3292 http://www.banglajol.info NATURAL CONVECTION FLOW IN A SQUARE CAVITY WITH INTERNAL HEAT GENERATION AND A
More informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
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 informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationChapter 9 NATURAL CONVECTION
Heat and Mass Transfer: Fundamentals & Applications Fourth Edition in SI Units Yunus A. Cengel, Afshin J. Ghajar McGraw-Hill, 2011 Chapter 9 NATURAL CONVECTION PM Dr Mazlan Abdul Wahid Universiti Teknologi
More informationConvection Workshop. Academic Resource Center
Convection Workshop Academic Resource Center Presentation Outline Understanding the concepts Correlations External Convection (Chapter 7) Internal Convection (Chapter 8) Free Convection (Chapter 9) Solving
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationFINDING RIGHTMOST EIGENVALUES OF LARGE SPARSE NONSYMMETRIC PARAMETERIZED EIGENVALUE PROBLEMS
FINDING RIGHTMOST EIGENVALUES OF LARGE SPARSE NONSYMMETRIC PARAMETERIZED EIGENVALUE PROBLEMS Department of Mathematics University of Maryland, College Park Advisor: Dr. Howard Elman Department of Computer
More informationAnalysis of Natural Convection Flow in a Trapezoidal Cavity Containing a Rectangular Heated Body in Presence of External Oriented Magnetic Field
Publications Available Online J. Sci. Res. 10 (1), 11-23 (2018) JOURNAL OF SCIENTIFIC RESEARCH www.banglajol.info/index.php/jsr Analysis of Natural Convection Flow in a Trapezoidal Cavity Containing a
More informationAdvanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell
Laminar external natural convection on vertical and horizontal flat plates, over horizontal and vertical cylinders and sphere, as well as plumes, wakes and other types of free flow will be discussed in
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationRAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8
HEFAT01 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 01 Malta RAYLEIGH-BÉNARD CONVECTION IN A CYLINDER WITH AN ASPECT RATIO OF 8 Leong S.S. School of Mechanical
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More informationEffect of an adiabatic fin on natural convection heat transfer in a triangular enclosure
American Journal of Applied Mathematics 2013; 1(4): 78-83 Published online November 10, 2013 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.20130104.16 Effect of an adiabatic fin on
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationAnalysis of Turbulent Free Convection in a Rectangular Rayleigh-Bénard Cell
Proceedings of the 8 th International Symposium on Experimental and Computational Aerothermodynamics of Internal Flows Lyon, July 2007 Paper reference : ISAIF8-00130 Analysis of Turbulent Free Convection
More informationCode Verification of Multiphase Flow with MFIX
Code Verification of Multiphase Flow with MFIX A N I R U D D H A C H O U D H A R Y ( A N I R U D D @ V T. E D U ), C H R I S T O P H E R J. ROY ( C J R O Y @ V T. E D U ), ( V I R G I N I A T E C H ) 23
More informationAbstract. Introduction
Combined forced and natural convection in a square cavity - numerical solution and scale analysis A.T. Franco/ M.M. Ganzarolli'' "DAMEC, CEFET, PR 80230-901, Curitiba, PR Brasil >>DE, FEM, UNICAMP 13081-970,
More informationFINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
More informationLINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION
LINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION É. FONTANA 1, E. MANCUSI 1,2, A. A. ULSON DE SOUZA 1, S. M. A. GUELLI U. SOUZA 1 1 Universidade Federal
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationNatural Convection in Parabolic Enclosure Heated from Below
www.ccsenet.org/mas Modern Applied Science Vol. 5, No. 3; June 011 Natural Convection in Parabolic Enclosure Heated from Below Dr. Ahmed W. Mustafa (Corresponding auther) University of Tikrit, College
More informationHessian-based model reduction for large-scale systems with initial condition inputs
Hessian-based model reduction for large-scale systems with initial condition inputs O. Bashir 1, K. Willcox 1,, O. Ghattas 2, B. van Bloemen Waanders 3, J. Hill 3 1 Massachusetts Institute of Technology,
More informationOn two fractional step finite volume and finite element schemes for reactive low Mach number flows
Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes
More information(Refer Slide Time: 0:35)
Fluid Dynamics And Turbo Machines. Professor Dr Shamit Bakshi. Department Of Mechanical Engineering. Indian Institute Of Technology Madras. Part A. Module-1. Lecture-4. Tutorial. (Refer Slide Time: 0:35)
More informationAC : A STUDENT PROJECT ON RAYLEIGH-BENARD CONVECTION
AC 2008-945: A STUDENT PROJECT ON RAYLEIGH-BENARD CONVECTION John Matsson, Oral Roberts University O. JOHN E. MATSSON is an Associate Professor of Mechanical Engineering and Chair of the Engineering, Physics
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationINSTRUCTOR: PM DR MAZLAN ABDUL WAHID
SMJ 4463: HEAT TRANSFER INSTRUCTOR: PM ABDUL WAHID http://www.fkm.utm.my/~mazlan TEXT: Introduction to Heat Transfer by Incropera, DeWitt, Bergman, Lavine 5 th Edition, John Wiley and Sons Chapter 9 Natural
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationSELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE
Computational Thermal Sciences, 3 (1): 63 72 (2011) SELF-SUSTAINED OSCILLATIONS AND BIFURCATIONS OF MIXED CONVECTION IN A MULTIPLE VENTILATED ENCLOSURE M. Zhao, 1, M. Yang, 1 M. Lu, 1 & Y. W. Zhang 2 1
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationCahn-Hilliard Phase Decomposition and
Decomposition and Roy H. Stogner 1 Graham F. Carey 1 1 Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin 9th US National Congress on Computational
More informationPHYSICAL MECHANISM OF CONVECTION
Tue 8:54:24 AM Slide Nr. 0 of 33 Slides PHYSICAL MECHANISM OF CONVECTION Heat transfer through a fluid is by convection in the presence of bulk fluid motion and by conduction in the absence of it. Chapter
More informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationFREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER
Free Convective Heat Transfer From an Object at Low Rayleigh Number FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER Md. Golam Kader and Khandkar Aftab Hossain * Department of Mechanical
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationLAMINAR NATURAL CONVECTION IN VERTICAL 2D GLAZING CAVITIES
Mechanical and Industrial Engineering University of Massachusetts, Amherst AMINAR NATURA CONVECTION IN VERTICA 2D GAZING CAVITIES Bhaskar Adusumalli ABSTRACT Finite element predictions of natural convection
More informationTHE INFLUENCE OF HEIGHT RATIO ON RAYLEIGH-NUMBER SCALING AND STABILITY OF HORIZONTAL CONVECTION
Seventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 9-11 December 2009 THE INFLUENCE OF HEIGHT RATIO ON RAYLEIGH-NUMBER SCALING AND STABILITY OF HORIZONTAL
More informationMarangoni-Natural Convection in Liquid Metals in the Presence of a Tilted Magnetic Field
Copyright 2010 Tech Science Press FDMP, vol.6, no.4, pp.369-384, 2010 Marangoni-Natural Convection in Liquid Metals in the Presence of a Tilted Magnetic Field S. Hamimid 1 and A. Amroune 1 Abstract: The
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationBounds on Surface Stress Driven Flows
Bounds on Surface Stress Driven Flows George Hagstrom Charlie Doering January 1, 9 1 Introduction In the past fifteen years the background method of Constantin, Doering, and Hopf has been used to derive
More informationModeling of Thermoelastic Damping in MEMS Resonators
Modeling of Thermoelastic Damping in MEMS Resonators T. Koyama a, D. Bindel b, S. Govindjee a a Dept. of Civil Engineering b Computer Science Division 1 University of California, Berkeley MEMS Resonators
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationCombined Effect of Magnetic field and Internal Heat Generation on the Onset of Marangoni Convection
International Journal of Fluid Mechanics & Thermal Sciences 17; 3(4): 41-45 http://www.sciencepublishinggroup.com/j/ijfmts doi: 1.11648/j.ijfmts.1734.1 ISSN: 469-815 (Print); ISSN: 469-8113 (Online) ombined
More informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
More informationPHYSICAL MECHANISM OF NATURAL CONVECTION
1 NATURAL CONVECTION In this chapter, we consider natural convection, where any fluid motion occurs by natural means such as buoyancy. The fluid motion in forced convection is quite noticeable, since a
More informationChapter 6 Laminar External Flow
Chapter 6 aminar Eternal Flow Contents 1 Thermal Boundary ayer 1 2 Scale analysis 2 2.1 Case 1: δ t > δ (Thermal B.. is larger than the velocity B..) 3 2.2 Case 2: δ t < δ (Thermal B.. is smaller than
More informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
More informationFLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.
FLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.Bhosale 1 Research scholar, 2 Head of department & Asst professor Department
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 informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 5 NATURAL CONVECTION HEAT TRANSFER BASIC CONCEPTS MECHANISM OF NATURAL
More informationNumerical Study of the Moving Boundary Problem During Melting Process in a Rectangular Cavity Heated from Below
American Journal of Applied Sciences 4 (4): 25-256, 2007 ISSN 546-9239 2007 Science Publications Corresponding Author: Numerical Study of the Moving Boundary Problem During Melting Process in a Rectangular
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationNumerical Continuation of Bifurcations - An Introduction, Part I
Numerical Continuation of Bifurcations - An Introduction, Part I given at the London Dynamical Systems Group Graduate School 2005 Thomas Wagenknecht, Jan Sieber Bristol Centre for Applied Nonlinear Mathematics
More informationMaximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
More informationTurbulent Natural Convection in an Enclosure with Colliding Boundary Layers
Turbulent Natural Convection in an Enclosure with Colliding Boundary Layers Abstract Mutuguta John Wanau 1* 1. School of Pure and Applied Sciences, Murang a University of Technology, P.O box 75-10200,
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationComputation of Unsteady Flows With Moving Grids
Computation of Unsteady Flows With Moving Grids Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de Unsteady Flows With Moving Boundaries, I Unsteady flows
More informationMIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER
Suranaree J. Sci. Technol. Vol. 17 No. 2; April - June 2010 139 MIXED CONVECTION IN A SQUARE CAVITY WITH A HEAT-CONDUCTING HORIZONTAL SQUARE CYLINDER Md. Mustafizur Rahman 1 *, M. A. Alim 1 and Sumon Saha
More informationOn a Data Assimilation Method coupling Kalman Filtering, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model
On a Data Assimilation Method coupling, MCRE Concept and PGD Model Reduction for Real-Time Updating of Structural Mechanics Model 2016 SIAM Conference on Uncertainty Quantification Basile Marchand 1, Ludovic
More informationSIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
More informationINFLUENCE OF GEOMETRICAL ASPECT RATIO ON THE OSCILLATORY MARANGONI CONVECTION IN LIQUID BRIDGES R. MONTI, R. SAVINO, M. LAPPA *
AUTHOR S POST PRINT (Romeo Colour: Green) Acta Astronautica (ISSN: 94-5765) Vol. 47, No., pp. 753-76, 2. DOI.6/S94-5765()26- Publisher version available at http://www.sciencedirect.com/science/article/pii/s94576526
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
More informationMIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM
THERMAL SCIENCE, Year 015, Vol. 19, No. 1, pp. 119-18 119 MIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM by Gurminder SINGH *a and Oluwole Daniel MAKINDE
More informationFree and Forced Convection Heat Transfer Characteristics in an Opened Box with Parallel Heated Plates
American Journal of Energy and Power Engineering 2015; 2(1): 1-11 Published online February 20, 2015 (http://www.aascit.org/journal/ajepe) ISSN: 2375-3897 Free and Forced Convection Heat Transfer Characteristics
More informationTwo-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5
Two-Dimensional Unsteady Flow in a Lid Driven Cavity with Constant Density and Viscosity ME 412 Project 5 Jingwei Zhu May 14, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of
More informationThe dependence of the cross-sectional shape on the hydraulic resistance of microchannels
3-weeks course report, s0973 The dependence of the cross-sectional shape on the hydraulic resistance of microchannels Hatim Azzouz a Supervisor: Niels Asger Mortensen and Henrik Bruus MIC Department of
More informationApplications in Fluid Mechanics
CHAPTER 8 Applications in Fluid 8.1 INTRODUCTION The general topic of fluid mechanics encompasses a wide range of problems of interest in engineering applications. The most basic definition of a fluid
More informationLarge scale continuation using a block eigensolver
Universidad Central de Venezuela Facultad de Ciencias Escuela de Computación Lecturas en Ciencias de la Computación ISSN 1316-6239 Large scale continuation using a block eigensolver Z. Castillo RT 2006-03
More information