Multiscale Computation for Incompressible Flow and Transport Problems
|
|
- Arabella Allison
- 6 years ago
- Views:
Transcription
1 Multiscale Computation for Incompressible Flow and Transport Problems Thomas Y. Hou Applied Mathematics, Caltech Collaborators: Y. Efenidev (TAMU), V. Ginting (Colorado), T. Strinopolous (Caltech), Danping Yang (Shandong Univ.), H. Ran (Caltech) Multiscale Computation for Incompressible Flow and Transport Problems p.1/66
2 Introduction Subsurface flows and transport are affected by heterogeneities at multiple scales (pore scale, core scale, field scale). Because of wide range of scales direct numerical simulations are not affordable. Upscaling of flow and transport parameters is commonly used in practice. Multiscale approaches are developed as an alternative to perform upscaling in the solution space. Multiscale Computation for Incompressible Flow and Transport Problems p.2/66
3 Outline Multiscale finite element methods (MsFEM) as a subgrid modeling technique; Applications of multiscale finite element methods to porous media flows; Upscaling of two-phase flow in flow-based coordinate system; Multiscale computation for the Navier-Stokes equations; Multiscale Computation for Incompressible Flow and Transport Problems p.3/66
4 Darcy s law and permeability Darcy s empirical law, 1856: The volumetric flux u(x, t) (Darcy velocity) is proportional to the pressure gradient u = k p = K p, µ where k(x) is the measured permeability of the rock, µ is the fluid viscosity, p(x) is the fluid pressure, u(x) is the Darcy velocity. We obtain the second order elliptic system u = K p in Q Darcy s Law div(u) = f in Q conservation Sand Water Multiscale Computation for Incompressible Flow and Transport Problems p.4/66
5 A natural porous media Rock Facies (1 1000m) Fault (1e 3m) Vugs (0.01m) Fracture (1e 4 m) Rock Facies Shale (0.01m) (1 1000m) Rock Grains and Voids (1e 5 m) Multiscale Computation for Incompressible Flow and Transport Problems p.5/66
6 Requirements/Challenges Accuracy and Robustness Valid for different types of subsurface heterogeneity Left: two-point variogram based. Middle: synthetic channelized. Right: channelized (North Sea). Applicable for varying flow scenarios Multiscale Computation for Incompressible Flow and Transport Problems p.6/66
7 Requirements/Challenges Accuracy and Robustness Valid for different types of subsurface heterogeneity Left: two-point variogram based. Middle: synthetic channelized. Right: channelized (North Sea). Applicable for varying flow scenarios Multiscale Computation for Incompressible Flow and Transport Problems p.6/66
8 Requirements/Challenges Accuracy and Robustness Valid for different types of subsurface heterogeneity Left: two-point variogram based. Middle: synthetic channelized. Right: channelized (North Sea). Applicable for varying flow scenarios Multiscale Computation for Incompressible Flow and Transport Problems p.6/66
9 Multiscale Finite Element Methods Babuska-Osborn (83,94), Hou and Wu (1997). Consider div(k ǫ (x) p ǫ ) = f, where ǫ is a small parameter. The central idea is to incorporate the small scale information into the finite element bases Basis functions are constructed by solving the leading order homogeneous equation in an element K (coarse grid or Representative Elementary Volume (RVE)) div(k ǫ (x) φ i ) = 0 in K It is through the basis functions that we capture the local small scale information of the differential operator. Multiscale Computation for Incompressible Flow and Transport Problems p.7/66
10 Multiscale Finite Element Methods Boundary conditions? φ i = linear function on K, φ i (x j ) = δ ij Coarse grid Fine grid Multiscale Computation for Incompressible Flow and Transport Problems p.8/66
11 Basis functions 40*(2.1+sin(60*(x y))+sin(60*y)) K bilinear 0 φ i = φ i 0 on K, where φi 0 are standard bilinear basis functions. Multiscale Computation for Incompressible Flow and Transport Problems p.9/66
12 Multiscale Finite Element Methods Except for the multiscale basis functions, MsFEM is the same as the traditional FEM (finite element method). Find p h ǫ V h = {φ i } such that k(p h ǫ, vh ) = f(v h ) v h V h, where k(u, v) = Z Q kij ǫ u v (x) dx, f(v) = x i x j Z Q fvdx The coupling of the small scales is through the variational formulation Solution has the form pǫ p 0 (x) + ǫn k (x/ǫ) x k p 0 (x), if k(x/ǫ) is periodic. If standard finite element method is used (linear basis functions): p ǫ p h ǫ H 1 (Q) Chα p ǫ H 1+α (Q) = O(hα ǫ α ). Thus, h ǫ, which is not affordable in practice. Multiscale Computation for Incompressible Flow and Transport Problems p.10/66
13 Relation to other approaches Subgrid modeling (by I. Babuska, T. Hughes, G. Allaire, T. Arbogast, and others) Subgrid stabilization (by F. Brezzi, L Franco, J.L. Guermond, T. Hughes, A. Russo, and others). Multiscale finite element methods (Hou and Wu, Allaire and Brizzi,...) Multiscale finite volume methods (Jenny, Lee and Tchelepi) Multiscale finite element methods based on two-scale convergence (C. Schwab,...) Variational multiscale method and residual free bubble methods (T. Hughes, F. Brezzi, T. Arbogast, G. Sangalli,...) Partition of Unity Methods (Babuska, Osborn, Fish,...) Wavelet based homogenization (Dorobantu, Engquist,...) Heterogeneous multiscale methods (E and Engquist...) Multiscale Computation for Incompressible Flow and Transport Problems p.11/66
14 Convergence property of MsFEM Consider k ǫ (x) = k(x/ǫ), where k(y) is periodic in y. h - computational mesh size. Theorem Denote p h ǫ the numerical solution obtained by MsFEM, and p ǫ the solution of the original problem. Then, If h >> ǫ, r ǫ p ǫ p h ǫ 1,Q C(h + h ) This theorem shows that MsFEM converges to the correct solution as ǫ 0 The ratio ǫ/h reflects two intrinsic scales. We call ǫ/h the resonance error The theorem shows that there is a scale resonance when h ǫ. Numerical experiments confirm the scale resonance. Multiscale Computation for Incompressible Flow and Transport Problems p.12/66
15 Resonance errors For problems with scale separation, we can choose h ǫ in order to avoid the resonance, but for problems with continuous spectrum of scales, we cannot avoid this resonance. To demonstrate the influence of the boundary condition of the basis function on the overall accuracy of the method we perform multiscale expansion of φ i Multiscale expansion of φ i φ i = φ 0 (x) + ǫφ 1 (x, x/ǫ) + ǫθ ǫ +..., φ1 (x, x/ǫ) = N k (x/ǫ) x k φ 0, where N k (x/ǫ) is a periodic function which depends on k(x/ǫ). θǫ satisfies div(k ǫ θ ǫ ) = 0 in K, θ i ǫ = φ 1(x, x/ǫ) + (φ i φ 0 )/ǫ on K Oscillations near the boundaries (in ǫ vicinity) of θ i ǫ lead to the resonance error Multiscale Computation for Incompressible Flow and Transport Problems p.13/66
16 Illustration of θ ǫ θ (without oversampling) 4 θ 0 y x Multiscale Computation for Incompressible Flow and Transport Problems p.14/66
17 Oversampling technique To capture more accurately the small scale information of the problem, the effect of θ ǫ needs to be moderated Since the boundary layer of θǫ is thin (O(ǫ)) we can sample in a domain with size larger than h + ǫ and use only interior sampled information to construct the basis functions. Let ψ k be the functions in the domain S, div(k ǫ (x) ψ k ) = 0 in S, ψ k = linear function on S, ψ k (s i ) = δ ik. S K Oversampled domain Coarse grid Fine grid Multiscale Computation for Incompressible Flow and Transport Problems p.15/66
18 Oversampling technique The base functions in a domain K S constructed as φ i K = X c ij ψ j K, φ i (x k ) = δ ik The method is non-conforming. The derivation of the convergence rate uses the homogenization method combined with the techniques of non-conforming finite element method (Efendiev et al., SIAM Num. Anal. 1999) By a correct choice of the boundary condition of the basis functions we can reduce the effects of the boundary layer in θ. Multiscale Computation for Incompressible Flow and Transport Problems p.16/66
19 Illustration of θ with oversampling θ (with oversampling) 2.8 θ 1.8 y x Multiscale Computation for Incompressible Flow and Transport Problems p.17/66
20 Numerical Results U h ǫ U h 0 l 2, ǫ/h = 0.64 h MsFEM MsFEM-os Resolved FEM l 2 rate l 2 rate h fine l 2 1/ e e-5 1/ e-4 1/ e e / e-4 1/ e e / e-4 1/ e e / e-4 Multiscale Computation for Incompressible Flow and Transport Problems p.18/66
21 MsFEM for problems with scale separation For periodic problems or problems with scale separation, multiscale finite element methods can take an advantage of scale separation. Local problems can be solved in RVE div(k φ i ) = 0 φ i = φ i 0 on RV E. Basis functions can be also approximated φ i = φ i 0 + N ǫ φ i 0, where φ 0 i is linear basis functions and N is the periodic solution of auxiliary problem in ǫ-size period div(k ǫ (x)( N + I)) = 0. (cf. Durlofsky 1981, etc.). Note, the above procedure works when homogenization by periodization is applicable (e.g., random homogeneous case). Multiscale Computation for Incompressible Flow and Transport Problems p.19/66
22 Various generalizations Once basis functions are constructed, various global formulation (mixed, control volume finite element, DG and etc) can be used to couple the subgrid effects. Control volume finite element: Find ph V h such that Z V z k(x) p h n dl = Z V z q dx V z Q, where V z is control volume. Mixed finite element: In each coarse block K, we construct basis functions for the velocity field div(k(x) wi K ) = 1 in K K ( 1 k(x) wi K on nk = e K e K i i 0 else. For the pressure, the basis functions are taken to be constants. Multiscale Computation for Incompressible Flow and Transport Problems p.20/66
23 Applications of MsFEM to subsurface flow simulations Two-phase flow model. Darcy s law for each phase v i = kk i(s i ) µ i p i, i=1,2. Here k - permeability field representing the heterogeneities (micro-level information), p i - the pressure, v i - velocity, k i - relative permeability, S i -saturation, µ i - viscosity phase 1 phase 2 Multiscale Computation for Incompressible Flow and Transport Problems p.21/66
24 Applications of MsFEM At least two way one can apply MsFEM 1) Solve the pressure equation on the coarse-grid and solve the saturation equation on the fine-grid div(λ(s)k p) = 0 S + v f(s) = 0, t where v = λ(s)k p. Basis functions are updated only near sharp fronts Coarse grid Fine grid 2) The coarse-scale equation for the transport is obtained and coupled with MsFEM for pressure equation. Multiscale Computation for Incompressible Flow and Transport Problems p.22/66
25 MsFVEM applied to two-phase flow problem Given S 0, for n = 1,2,3,, do the following: find p n 1 h V h such that Z (IM)plicit (P)ressure (E)xplicit (S)aturation: V z λ(s n 1 )k(x) p n 1 h n dl = compute v n 1 = λ(s n 1 )k(x) p n 1 h time march on the saturation equation: Z V z q dx V z Q Z `Sn S n 1 Z Z dx + t n 1 f(s n 1 )v n 1 n dl = t n 1 c z c z c z q S dx Multiscale Computation for Incompressible Flow and Transport Problems p.23/66
26 Numerical Setting Rectangular domain is considered. The permeability field is generated using geostatistical libraries. The boundary conditions: no flow on top and bottom boundaries, a fixed pressure and saturation (S = 1) at the inlet (left edge), fixed pressure at the outlet (right edge). The production rate F = q0 /q, where q 0 the volumetric flow rate of oil produced at the outlet edge and q the volumetric flow rate of the total fluid produced at the outlet edge. The dimensionless time is defined as PV I = qt/v p, where t is time, V p is the total pore volume of the system. F(t) no flow S=1 p=1 p= no flow Multiscale Computation for Incompressible Flow and Transport Problems p.24/66
27 A random porosity field with layered structure Multiscale Computation for Incompressible Flow and Transport Problems p.25/66
28 Recovery of fine scale velocity Left: Fine grid horizontal velocity field, N = Right: Recovered horizontal velocity field from the coarse grid N = using multiscale bases. Multiscale Computation for Incompressible Flow and Transport Problems p.26/66
29 Recovery of fine scale saturation Left: Fine grid saturation, N = Right: Saturation computed from the recovered fine grid velocity field. Multiscale Computation for Incompressible Flow and Transport Problems p.27/66
30 Fractional flow and total flow curves fine modified MsFVM standard MsFVM fine modified MsFVM standard MsFVM F PVI Q PVI Fractional flow and total flow for a realization of permeability field with exponential variogram and l x = 0.4, l z = 0.02, σ = 1.5. Multiscale Computation for Incompressible Flow and Transport Problems p.28/66
31 Channelized permeability fields Benchmark tests: SPE 10 Comparative Project Multiscale Computation for Incompressible Flow and Transport Problems p.29/66
32 Channelized reservoir 1 fine scale model standard MsFVEM 350 fine scale model standard MsFVEM F Q PVI PVI Comparison of upscaled quantities (Layer 43) Multiscale Computation for Incompressible Flow and Transport Problems p.30/66
33 Channelized reservoir z 0.5 z x 0 x 0 Comparison of saturation profile at PVI=0.5: (left) fine-scale model, (right) standard MsFVEM Multiscale Computation for Incompressible Flow and Transport Problems p.31/66
34 MsFVEM utilizing limited global information The numerical tests using strongly channelized permeability fields (such as SPE 10 Comparative) show that local basis functions can not accurately capture the long-range information. There is a need to incorporate a global information. The main idea is to use the solution of the fine-scale problem at time zero, p 0, to determine the boundary conditions for the multiscale basis formulation. φ ( x )=1 i i x i φ ( x i+1 )=0 i x i+1 φ ( x)=0 i x i 1 φ ( x i 1 )=0 i φ i ( x)=0 These approach is different from oversampling technique. Previous related work: J. Aarnes; L. Durlofsky et al. Multiscale Computation for Incompressible Flow and Transport Problems p.32/66
35 Channelized reservoir fine scale model modified MsFVEM standard MsFVEM fine scale model modified MsFVEM standard MsFVEM F Q PVI PVI Comparison of upscaled quantities Multiscale Computation for Incompressible Flow and Transport Problems p.33/66
36 Channelized reservoir z 0.5 z x 0 x 0 Comparison of saturation profile at PVI=0.5: (left) fine-scale model, (right) modified MsFVEM Multiscale Computation for Incompressible Flow and Transport Problems p.34/66
37 Channelized reservoir fine modified MsFVM fine modified MsFVM F 0.5 Q PVI PVI Comparison of upscaled quantities (Layer 43, changing boundary conditions) Multiscale Computation for Incompressible Flow and Transport Problems p.35/66
38 Channelized reservoir fine scale saturation plot at PVI=0.7 1 saturation plot at PVI=0.5 using modified MsFVEM Comparison of saturation profile at PVI=0.5: (left) fine-scale model (right) modified MsFVEM (changing boundary condition) Multiscale Computation for Incompressible Flow and Transport Problems p.36/66
39 Brief Analysis Main goal is to show that time-varying pressure is strongly influenced by the initial pressure field. Use the streamline-pressure coordinates: ψ/ x 1 = v 2, ψ/ x 2 = v 1 Set η = ψ(x, t = 0) and ζ = p(x, t = 0) and transform as follows: y x η ζ Ο(δ) Ο(δ) high flow channel Multiscale Computation for Incompressible Flow and Transport Problems p.37/66
40 Brief Analysis The transformed pressure equation: η k 2 λ(s) p «+ η ζ λ(s) p «ζ = 0 The transformed saturation equation: S t + (v η) f(s) η + (v ζ) f(s) ζ = 0 k 2 λ(s) = k 0 2 λ 0 (ζ, t)1 Q1 δ + k 1 2 λ 1 (η, ζ, t)1 Qδ, λ(s) = λ 0 (ζ, t)1 Q1 δ + λ 1 (η, ζ, t)1 Qδ. The pressure has the following expansion: p(η, ζ, t) = p 0 (ζ, t) + δp 1 (η, ζ, t) +..., λ ζ 0 (ζ, t) p 0 = 0. ζ Modified basis functions can exactly recover the initial pressure. Multiscale Computation for Incompressible Flow and Transport Problems p.38/66
41 Analysis Assumption G. There exists a sufficiently smooth scalar valued function G(η) (G C 3 ), such that p G(p sp ) 1,Q Cδ, where δ is sufficiently small. Under Assumption G and p sp W 1,s (Q) (s > 2), multiscale finite element method converges with the rate given by p p h 1,Q Cδ + Ch 1 2/s p sp W 1,s (Q) Cδ + Ch1 2/s. Multiscale Computation for Incompressible Flow and Transport Problems p.39/66
42 Extensions Assume there exists a sufficiently smooth scalar valued function G(η), η R N (G C 3 ), such that p G(u 1,..., u N ) 1,Q Cδ, where δ is sufficiently small. Let ω i be a patch, and define φ 0 i to be piecewise linear basis function in patch ω i, such that φ 0 i (x j) = δ ij. For simplicity of notation, denote u 1 = 1. Then, the multiscale finite element method for each patch ω i is constructed by ψ ij = φ 0 i u j where j = 1,.., N and i is the index of nodes. First, we note that in each K, P n i=1 ψ ij = u j is the desired single-phase flow solution. Theorem. Assume u i W 1,s (Q), s > 2, i = 1,..., N. Then p p h 1,Q Cδ + Ch 1 2/s. Multiscale Computation for Incompressible Flow and Transport Problems p.40/66
43 An approach to general 2nd order pdes (H. Owhadi and L. Zhang, CPAM, 2006) div(λ(x)k(x) p) = 0, where λ(x) is smooth function, while k(x) is rough (e.g., k(x) = k(x/ǫ)). Take u 1 and u 2 that satisfy div(k(x) u i ) = 0 in Q, u i = x i on Q. Then, p(u 1, u 2 ) W 2,p because it satisfies a ij 2 p u i u j 0. (e.g., p = p 0 + ǫ 2 N(x/ǫ)...). Owhadi and Zhang showed that the non-conforming method with basis functions that span u 1 and u 2 converge. Multiscale Computation for Incompressible Flow and Transport Problems p.41/66
44 Two-phase flow equations in flow-based coor. ψ k 2 λ(s) P «+ ψ p λ(s) P p «= 0. S t + (v ψ) f(s) ψ + (v p) f(s) p = 0. Consider λ(s) = 1. Homogenization of hyperbolic equations. S ǫ t + vǫ 0 f(sǫ ) p = 0 S(p, ψ, t = 0) = S 0, v ǫ 0 (p) = v 0(p, p ǫ ). Multiscale Computation for Incompressible Flow and Transport Problems p.42/66
45 Homogenization of transport Then, for each ψ, it can be shown that S ǫ (p, ψ, t) S(p, ψ, t) in L 1 ((0,1) (0, T)), where S satisfies S t + ṽ 0 f( S) p = 0, where ṽ 0 is harmonic average of v ǫ 0, i.e., 1 v ǫ 0 1 ṽ 0 weak in L (0,1). Theorem. S ǫ S n Gǫ 1/n. Note. S can be considered as an upscaled S ǫ along streamlines. Can we average across streamlines? Multiscale Computation for Incompressible Flow and Transport Problems p.43/66
46 Numerical results Permeability fields used in the simulations. Left - permeability field with exponential variogram, middle - synthetic channelized permeability field, right - layer 36 of SPE comparative project Multiscale Computation for Incompressible Flow and Transport Problems p.44/66
47 Numerical results Saturation snapshots for variogram based permeability field (top) and synthetic channelized permeability field (bottom). Linear flux is used. Left figures represent the upscaled saturation plots and the right figures represent the fine-scale saturation plots. Multiscale Computation for Incompressible Flow and Transport Problems p.45/66
48 Numerical results Saturation snapshots for variogram based permeability field (top) and synthetic channelized permeability field (bottom). Nonlinear flux is used. Left figures represent the upscaled saturation plots and the right figures represent the fine-scale saturation plots. Multiscale Computation for Incompressible Flow and Transport Problems p.46/66
49 Numerical results Upscaling error for permeability generated using two-point geostatistics LINEAR FLUX 25x25 50x50 100x x200 L 1 error of S e e e 5 L 1 error of S with macrodispersion L 1 error of S fine without macrodispersion NONLINEAR FLUX 25x25 50x50 100x x200 L 1 error of S e e e 4 L 1 error of S with macrodispersion L 1 error of S fine without macrodispersion Multiscale Computation for Incompressible Flow and Transport Problems p.47/66
50 Numerical results Upscaling error for for synthetic channelized permeability field LINEAR FLUX 25x25 50x50 100x x200 L 1 error of S L 1 error of S with macrodispersion L 1 error of S fine without macrodispersion NONLINEAR FLUX 25x25 50x50 100x x200 L 1 error of S L 1 error of S with macrodispersion L 1 error of S fine without macrodispersion Multiscale Computation for Incompressible Flow and Transport Problems p.48/66
51 Numerical results Upscaling error for SPE 10, layer 36 LINEAR FLUX 25x25 50x50 100x x200 L 1 error of S L 1 error of S with macrodispersion L 1 error of S fine without macrodispersion NONLINEAR FLUX 25x25 50x50 100x x200 L 1 error of S L 1 error of S with macrodispersion L 1 error of S fine without macrodispersion Multiscale Computation for Incompressible Flow and Transport Problems p.49/66
52 Numerical results Computational cost fine x.y fine p, ψ S S layered, linear flux layered, nonlinear flux percolation, linear flux percolation, nonlinear flux SPE10 36, linear flux SPE10 36, nonlinear flux Multiscale Computation for Incompressible Flow and Transport Problems p.50/66
53 Numerical results Left: Pressure and streamline function at time t = 0.4 in Cartesian frame. Right: pressure and streamline function at time t = 0.4 in initial pressure-streamline frame. For two-phase flow, equations are upscaled on flow-based coordinate system. MsFEM using limited global information is equivalent to standard MsFEM. Multiscale Computation for Incompressible Flow and Transport Problems p.51/66
54 Multiscale analysis for Navier-Stokes Equations t u ǫ + (u ǫ )u ǫ = p ǫ + ν u ǫ + f, u ǫ = 0, u ǫ t=0 = U(x) + W(x, x ǫ ). Motivated by a multiscale analysis for the Euler equations, we look for (u ǫ, p ǫ ) of form as u ǫ = ū(t,x, τ) + eu(t, θ, τ,z), p ǫ = p(t,x, τ) + ep(t, θ, τ,z), where eu and ep are periodic in z with zero mean, z = θ/ǫ. Multiscale Computation for Incompressible Flow and Transport Problems p.52/66
55 Main Results. The main results can be summarized as follows: Effective equation for ū: t ū + (ū x )ū + x w w = x p + ν ū + f, x ū = 0, ū t=0 = U(x). Leading order small scale velocity field w: τ w + z w x θw + x θ z q = ν ǫ z ( x θ x θ z w), with ` x θ z w = 0 and, with w τ=t=0 = W. Effective equation for flow map θ: t θ + (ū x ) θ + ǫ x Θ w = 0, θ t=0 = x. t Θ + (I + z Θ) x θw = 0, with Θ t=0 = 0. Multiscale Computation for Incompressible Flow and Transport Problems p.53/66
56 Numerical Experiments for 2D NSE We use direct numerical simulation (DNS) to check the accuracy of our multiscale analysis. In computation, we choose ν = The DNS is performed using 512 by 512 mesh. The multiscale computations use 64 by 64 on the coarse grid, and 32 by 32 for the subgrid. The initial vorticity is generated in Fourier space, satisfying bω k O( k 1 ) with random phases. Our main objective is to test how accuate we can capture the inverse cascade of 2D turbulence. Multiscale Computation for Incompressible Flow and Transport Problems p.54/66
57 Numerical results in 2D Vorticity contour at t = 0. Multiscale Computation for Incompressible Flow and Transport Problems p.55/66
58 Numerical results in 2D Vorticity contour, DNS (left) and Homogenization (right) at t=5.0. Multiscale Computation for Incompressible Flow and Transport Problems p.56/66
59 Numerical results in 2D Vorticity contour, DNS (left) and Homogenization (right) at t=10.0. Multiscale Computation for Incompressible Flow and Transport Problems p.57/66
60 Numerical results in 2D Vorticity contour, DNS (left) and Homogenization (right) at t=20.0. Multiscale Computation for Incompressible Flow and Transport Problems p.58/66
61 Numerical results in 2D Spectrum of velocity component u 1.(dashed), t = 0.0; (solid), DNS, t = 20.0; (chndot), homogenization, t = Multiscale Computation for Incompressible Flow and Transport Problems p.59/66
62 Numerical Experiments for 3D forced NSE DNS is performed using mesh in a periodic box to check accuracy of the multiscale method. The Taylor Reynolds number is 223 in our computation (equvilent viscosity is ). The specific form of the forcing is given by ˆf k = δ N û k pûk û k (-12) N is the number of wave modes that are forced. The above form of forcing ensures that the energy injection rate, P ˆf û, is a constant which is equal to δ. We chose δ = 0.1 for all of our runs. Multiscale Computation for Incompressible Flow and Transport Problems p.60/66
63 3D Results e k t Temporal evolution of total kinetic energy. Multiscale Computation for Incompressible Flow and Transport Problems p.61/66
64 3D Results Iso-surface plot of vorticity. Multiscale Computation for Incompressible Flow and Transport Problems p.62/66
65 3D Results Comparison of energy spectrum at t = (blue), DNS. (red), N=128; (green), N=64; (yellow), N=32. The computations are performed for 30 eddy turn-over times. Multiscale Computation for Incompressible Flow and Transport Problems p.63/66
66 3D Results PDF of ω e 1 (solid), ω e 2 (dashed line), ω e 3 (dotted line). DNS (left) vs LES (right), at t=20. Multiscale Computation for Incompressible Flow and Transport Problems p.64/66
67 3D Results PDF of e 3 t 1 (solid), e 2 t 2 (dashed line), e 3 t 1 (dotted line). DNS (left) vs LES (right), at t=20. Here t j are eigenvectors of the SGS tensor. Multiscale Computation for Incompressible Flow and Transport Problems p.65/66
68 Conclusions Multiscale finite element methods provide an effective multiscale computational method for flow and transport in heterogeneous porous media. Resonance error can be effectively removed by the over-sampling technique. Long range scale interaction can be captured by incorporating limited global information into the multiscale bases. Flow-based adaptive coordinates offer a natural physical coordinate to capture the long range interaction in a channelized media. The hybrid Eulerian-Lagrangian coordinates offer a natural multiscale frame to capture the nonlinear nonlocal scale interactions for the Navier-Stokes equations. Multiscale Computation for Incompressible Flow and Transport Problems p.66/66
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationA MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS
A MULTISCALE METHOD FOR MODELING TRANSPORT IN POROUS MEDIA ON UNSTRUCTURED CORNER-POINT GRIDS JØRG E. AARNES AND YALCHIN EFENDIEV Abstract. methods are currently under active investigation for the simulation
More informationMultiscale Finite Element Methods. Theory and
Yalchin Efendiev and Thomas Y. Hou Multiscale Finite Element Methods. Theory and Applications. Multiscale Finite Element Methods. Theory and Applications. October 14, 2008 Springer Dedicated to my parents,
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationNumerical Simulation of Flows in Highly Heterogeneous Porous Media
Numerical Simulation of Flows in Highly Heterogeneous Porous Media R. Lazarov, Y. Efendiev, J. Galvis, K. Shi, J. Willems The Second International Conference on Engineering and Computational Mathematics
More informationA GENERALIZED CONVECTION-DIFFUSION MODEL FOR SUBGRID TRANSPORT IN POROUS MEDIA
MULTISCALE MODEL. SIMUL. Vol. 1, No. 3, pp. 504 526 c 2003 Society for Industrial and Applied Mathematics A GENERALIZED CONVECTION-DIFFUSION MODEL FOR SUBGRID TRANSPORT IN POROUS MEDIA Y. EFENDIEV AND
More informationXingye Yue. Soochow University, Suzhou, China.
Relations between the multiscale methods p. 1/30 Relationships beteween multiscale methods for elliptic homogenization problems Xingye Yue Soochow University, Suzhou, China xyyue@suda.edu.cn Relations
More informationMULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationMMsFEM and Streamlines
MMsFEM and Streamlines Combining a mixed multiscale FEM with streamline simulation for enhanced reservoir performance prediction on large grid models. Jørg E. Aarnes, Vegard Kippe and Knut-Andreas Lie
More informationMULTISCALE ANALYSIS AND COMPUTATION FOR THE THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER STOKES EQUATIONS
MULTISCALE MODEL. SIMUL. Vol. 6, No. 4, pp. 1317 1346 c 2008 Society for Industrial and Applied Mathematics MULTISCALE ANALYSIS AND COMPUTATION FOR THE THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER STOKES EQUATIONS
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationParallel Simulation of Subsurface Fluid Flow
Parallel Simulation of Subsurface Fluid Flow Scientific Achievement A new mortar domain decomposition method was devised to compute accurate velocities of underground fluids efficiently using massively
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationA Framework for Modeling Subgrid Effects for Two-Phase Flows in Porous Media
A Framework for Modeling Subgrid Effects for Two-Phase Flows in Porous Media Thomas Y. Hou Andrew Westhead Danping Yang Abstract In this paper, we study upscaling for two-phase flows in strongly heterogeneous
More informationPRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS
PRECONDITIONING MARKOV CHAIN MONTE CARLO SIMULATIONS USING COARSE-SCALE MODELS Y. EFENDIEV, T. HOU, AND W. LUO Abstract. We study the preconditioning of Markov Chain Monte Carlo (MCMC) methods using coarse-scale
More informationMULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION
MULTISCALE FINITE ELEMENT METHODS FOR STOCHASTIC POROUS MEDIA FLOW EQUATIONS AND APPLICATION TO UNCERTAINTY QUANTIFICATION P. DOSTERT, Y. EFENDIEV, AND T.Y. HOU Abstract. In this paper, we study multiscale
More informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationMixed Multiscale Methods for Heterogeneous Elliptic Problems
Mixed Multiscale Methods for Heterogeneous Elliptic Problems Todd Arbogast Abstract We consider a second order elliptic problem written in mixed form, i.e., as a system of two first order equations. Such
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationModeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids
Modeling of two-phase flow in fractured porous media on unstructured non-uniform coarse grids Jørg Espen Aarnes and Vera Louise Hauge SINTEF ICT, Deptartment of Applied Mathematics Applied Mathematics
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationarxiv: v1 [math.na] 4 Jun 2014
MIXE GENERALIZE MULTISCALE FINITE ELEMENT METHOS AN APPLICATIONS ERIC T. CHUNG, YALCHIN EFENIEV, AN CHA SHING LEE arxiv:1406.0950v1 [math.na] 4 Jun 2014 Abstract. In this paper, we present a mixed Generalized
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationSeminar zu aktuellen Themen der Numerik im Wintersemester 2010/2011
Seminar zu aktuellen Themen der Numerik im Wintersemester 2010/2011 Modeling Two-Phase-Two-Component Processes in Porous Media Using an Upscaling-Multi-Scale Method Elin Solberg 27 January 2011 Contents
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationComparison of the Effects of k-ϵ, k-ω, and Zero Equation Models on Characterization of Turbulent Permeability of Porous Media
Comparison of the Effects of k-ϵ, k-ω, and Zero Equation Models on Characterization of Turbulent Permeability of Porous Media Charlie Matsubara *, Tim Kuo, Helen Wu PIM Consulting, No. 16, Gonyequ 1 st
More informationModeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures
Modeling seismic wave propagation during fluid injection in a fractured network: Effects of pore fluid pressure on time-lapse seismic signatures ENRU LIU, SERAFEIM VLASTOS, and XIANG-YANG LI, Edinburgh
More informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More informationUniversal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis.
Universal Scalable Robust Solvers from Computational Information Games and fast eigenspace adapted Multiresolution Analysis Houman Owhadi Joint work with Clint Scovel IPAM Apr 3, 2017 DARPA EQUiPS / AFOSR
More informationarxiv: v1 [math.na] 25 Aug 2016
A conservative local multiscale model reduction technique for Stoes flows in heterogeneous perforated domains Eric T. Chung Maria Vasilyeva Yating Wang arxiv:1608.07268v1 [math.na] 25 Aug 2016 October
More informationWeak Galerkin Finite Element Methods and Applications
Weak Galerkin Finite Element Methods and Applications Lin Mu mul1@ornl.gov Computational and Applied Mathematics Computationa Science and Mathematics Division Oak Ridge National Laboratory Georgia Institute
More informationMultiscale Computations for Flow and Transport in Porous Media
Multiscale Computations for Flow and Transport in Porous Media Thomas Y. Hou Applied Mathematics, 217-5, Caltech, Pasadena, CA 91125, USA. Email: hou@ama.caltech.edu. Abstract. Many problems of fundamental
More informationICES REPORT Adaptive Numerical Homogenization for Non-Linear Multiphase Flow and Transport
ICES REPORT 17-13 June 217 Adaptive Numerical Homogenization for Non-Linear Multiphase Flow and Transport by Gurpreet Singh, Yerlan Amanbek, and Mary F. Wheeler The Institute for Computational Engineering
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
More informationEddy viscosity of cellular flows by upscaling
Eddy viscosity of cellular flows by upscaling Alexei Novikov a a California Institute of Technology, Applied & Computational Mathematics 1200 E. California Boulevard, MC 217-50, Pasadena, CA 91125, USA
More informationMULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS
COMM. MATH. SCI. Vol. 2, No. 4, pp. 553 589 c 2004 International Press MULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS Y. EFENDIEV, T. HOU, AND V. GINTING Abstract. In this
More informationTwo-Scale Wave Equation Modeling for Seismic Inversion
Two-Scale Wave Equation Modeling for Seismic Inversion Susan E. Minkoff Department of Mathematics and Statistics University of Maryland Baltimore County Baltimore, MD 21250, USA RICAM Workshop 3: Wave
More informationModeling two-phase flow in strongly heterogeneous porous media
Presented at the COMSOL Conference 2010 China COMSOL 2010 I«^rc Modeling two-phase flow in strongly heterogeneous porous media Zhaoqin Huang Research Center for Oil & Gas Flow in Reservoir, School of Petroleum
More informationTENSOR RELATIVE PERMEABILITIES: ORIGINS, MODELING AND NUMERICAL DISCRETIZATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 9, Number 3, Pages 701 724 c 2012 Institute for Scientific Computing and Information TENSOR RELATIVE PERMEABILITIES: ORIGINS, MODELING AND
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationOnline basis construction for goal-oriented adaptivity in the Generalized Multiscale Finite Element Method
Online basis construction for goal-oriented adaptivity in the Generalized Multiscale Finite Element Method Eric Chung, Sara Pollock, Sai-Mang Pun arxiv:1812.02290v1 [math.na] 6 Dec 2018 December 7, 2018
More informationUpscaling for Two-phase Flows in Porous Media
Upscaling for Two-phase Flows in Porous Media Thesis by Andrew Westhead In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California
More informationModeling of 1D Anomalous Diffusion In Fractured Nanoporous Media
LowPerm2015 Colorado School of Mines Low Permeability Media and Nanoporous Materials from Characterisation to Modelling: Can We Do It Better? IFPEN / Rueil-Malmaison - 9-11 June 2015 CSM Modeling of 1D
More informationEddy viscosity of cellular flows by upscaling
Eddy viscosity of cellular flows by upscaling Alexei Novikov a a California Institute of Technology, Applied & Computational Mathematics 200 E. California Boulevard, MC 27-50, Pasadena, CA 925, USA E-mail:
More informationREDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION
REDUCING ORDER METHODS APPLIED TO RESERVOIR SIMULATION Lindaura Maria Steffens Dara Liandra Lanznaster lindaura.steffens@udesc.br daraliandra@gmail.com Santa Catarina State University CESFI, Av. Central,
More informationOperator Upscaling for the Wave Equation
Operator Upscaling for the Wave Equation Tetyana Vdovina Susan E. Minkoff UMBC), Oksana Korostyshevskaya Department of Computational and Applied Mathematics Rice University, Houston TX vdovina@caam.rice.edu
More informationMultiscale Methods for Subsurface Flow. SINTEF ICT, Dept. of Applied Mathematics
Multiscale Methods for Subsurface Flow Jørg Aarnes, Knut Andreas Lie, Stein Krogstad, and Vegard Kippe SINTEF ICT, Dept. of Applied Mathematics... and for saving our planet Applied Mathematics 1/89 Subsurface
More informationReactive Transport in Porous Media
Reactive Transport in Porous Media Daniel M. Tartakovsky Department of Mechanical & Aerospace Engineering University of California, San Diego Alberto Guadagnini, Politecnico di Milano, Italy Peter C. Lichtner,
More informationMULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS
MULTISCALE FINITE ELEMENT METHODS FOR NONLINEAR PROBLEMS AND THEIR APPLICATIONS Y. EFENDIEV, T. HOU, AND V. GINTING Abstract. In this paper we propose a generalization of multiscale finite element methods
More informationOptimizing calculation costs of tubulent flows with RANS/LES methods
Optimizing calculation costs of tubulent flows with RANS/LES methods Investigation in separated flows C. Friess, R. Manceau Dpt. Fluid Flow, Heat Transfer, Combustion Institute PPrime, CNRS University
More informationA NUMERICAL APPROXIMATION OF NONFICKIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXX, 1(21, pp. 75 84 Proceedings of Algoritmy 2 75 A NUMERICAL APPROXIMATION OF NONFICIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA R. E. EWING, Y. LIN and J. WANG
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
More informationFracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood
Fracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood R. Sanaee *, G. F. Oluyemi, M. Hossain, and M. B. Oyeneyin Robert Gordon University *Corresponding
More informationNUMERICAL METHODS FOR MULTISCALE ELLIPTIC PROBLEMS
NUMERICAL METODS FOR MULTISCALE ELLIPTIC PROBLEMS PINGBING MING AND XINGYE YUE Abstract. We present an overview of the recent development on numerical methods for elliptic problems with multiscale coefficients.
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 informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationA FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS
A FRONT-TRACKING METHOD FOR HYPERBOLIC THREE-PHASE MODELS Ruben Juanes 1 and Knut-Andreas Lie 2 1 Stanford University, Dept. Petroleum Engineering, USA 2 SINTEF IKT, Dept., Norway ECMOR IX, August 30 September
More informationA NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS
June 30 - July 3, 2015 Melbourne, Australia 9 7B-4 A NOVEL VLES MODEL FOR TURBULENT FLOW SIMULATIONS C.-Y. Chang, S. Jakirlić, B. Krumbein and C. Tropea Institute of Fluid Mechanics and Aerodynamics /
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationOversampling for the partition of unity parallel finite element algorithm
International conference of Computational Mathematics and its application Urumqi, China, July 25-27, 203 Oversampling for the partition of unity parallel finite element algorithm HAIBIAO ZHENG School of
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationOn the (multi)scale nature of fluid turbulence
On the (multi)scale nature of fluid turbulence Kolmogorov axiomatics Laurent Chevillard Laboratoire de Physique, ENS Lyon, CNRS, France Laurent Chevillard, Laboratoire de Physique, ENS Lyon, CNRS, France
More informationFlow-based Upscaling of Reservoir Models
Flow-based Upscaling of Reservoir Models Knut Andreas Lie Department of Mathematics and Cybernetics, SINTEF Digital/ Department of Mathematical Sciences, NTNU, Norway Multiscale Methods Summer School June
More informationWeierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, Berlin, Germany,
Volker John On the numerical simulation of population balance systems Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstr. 39, 10117 Berlin, Germany, Free University of Berlin, Department
More informationSplitting methods in the design of coupled flow and mechanics simulators
Splitting methods in the design of coupled flow and mechanics simulators Sílvia Barbeiro CMUC, Department of Mathematics, University of Coimbra PhD Program in Mathematics Coimbra, November 17, 2010 Sílvia
More informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationMultiscale Methods for Multiphase Flow in Porous Media
Multiscale Methods for Multiphase Flow in Porous Media Jan M. Nordbotten Department of Mathematics, University of Bergenjan.nordbotten@math.uib.no Department of Civil and Environmental Engineering, Princeton
More informationHETEROGENOUS CARBONATES INTEGRATING PLUG AND WHOLE CORE DATA USING ROCK TYPES
SCA2012-12 1/12 HETEROGENOUS CARBONATES INTEGRATING PLUG AND WHOLE CORE DATA USING ROCK TYPES Mark Skalinski, Rafael Salazar, Gerry LaTorraca, Zheng Yang, and John Urbach Chevron ETC This paper was prepared
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationNUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 NUMERICAL MODELING OF FLOW THROUGH DOMAINS WITH SIMPLE VEGETATION-LIKE OBSTACLES Steven
More informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationTodd Arbogast. Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and Sciences (ICES)
CONSERVATIVE CHARACTERISTIC METHODS FOR LINEAR TRANSPORT PROBLEMS Todd Arbogast Department of Mathematics and, (ICES) The University of Texas at Austin Chieh-Sen (Jason) Huang Department of Applied Mathematics
More informationUncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media
Uncertainty Quantification of Two-Phase Flow in Heterogeneous Porous Media M.Köppel, C.Rohde Institute for Applied Analysis and Numerical Simulation Inria, Nov 15th, 2016 Porous Media Examples: sponge,
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
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 informationChapter Seven. For ideal gases, the ideal gas law provides a precise relationship between density and pressure:
Chapter Seven Horizontal, steady-state flow of an ideal gas This case is presented for compressible gases, and their properties, especially density, vary appreciably with pressure. The conditions of the
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationHOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM
Homogenization of a heat transfer problem 1 G. Allaire HOMOGENIZATION OF A CONVECTIVE, CONDUCTIVE AND RADIATIVE HEAT TRANSFER PROBLEM Work partially supported by CEA Grégoire ALLAIRE, Ecole Polytechnique
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationImproved discontinuity-capturing finite element techniques for reaction effects in turbulence computation
Comput. Mech. (26) 38: 356 364 DOI.7/s466-6-45- ORIGINAL PAPER A. Corsini F. Rispoli A. Santoriello T.E. Tezduyar Improved discontinuity-capturing finite element techniques for reaction effects in turbulence
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationA posteriori error estimation in the FEM
A posteriori error estimation in the FEM Plan 1. Introduction 2. Goal-oriented error estimates 3. Residual error estimates 3.1 Explicit 3.2 Subdomain error estimate 3.3 Self-equilibrated residuals 3.4
More informationUpscaling non-darcy flow
Transport in Porous Media manuscript No. (will be inserted by the editor) Upscaling non-darcy flow C.R. Garibotti M.Peszyńska Received: date / Accepted: date Abstract We consider upscaling of non-darcy
More information