Semi-Lagrangian kinetic and gyrokinetic simulations
|
|
- Harvey Dennis
- 5 years ago
- Views:
Transcription
1 Semi-Lagrangian kinetic and gyrokinetic simulations M. Mehrenberger IRMA, University of Strasbourg and TONUS project (INRIA) Marseille, Equations cinétiques, Collaborators : B. Afeyan, A. Back, F. Casas, E. Chacon Golcher, A. Dodhy, N. Crouseilles, E. Faou, P. Glanc, V. Grandgirard, A. Hamiaz, P. Helluy, S. Hirstoaga, G. Latu, E. Madaule, L. Mendoza, P. Navaro, E. Sonnendrücker, M. Ottaviani, C. Prouveur, F. Rozar, C. Steiner M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 1 / 40
2 Outline Overview of recent works around kinetic and gyrokinetic simulations 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) 2. 2D conservative Semi-Lagrangian scheme (11-13) 3. Curvilinear and hexagonal meshes (14-17) 4. Gyroaverage operator on polar grids (18-33) 5. Field aligned interpolation (34-37) 6. Conclusion/perspectives/references (38-40) M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 2 / 40
3 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) Strang splitting Since Cheng-Knorr (1976), most Vlasov codes use Strang splitting. Transport in x over t f (t, x, v)+v@ x f (t, x, v) =0, t! t + Update of E through the Poisson equation Transport in v over t t 2 Transport in x over t f (t, x, v)+e(t, x)@ v f (t, x, v) t f (t, x, v)+v@ x f (t, x, v) =0, t! t + Splitting is parametrized by s = 3, begin_with_x = true, a 1 = 1/2, a 2 = 1, a 3 = 1/2. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 3 / 40 t 2
4 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) Abstract form Hamiltonian is defined as H(f )=T(f)+U(f), with Z v 2 Z 1 H(f ) = f (x, v)dxdv E(f )(x) 2 dx. (0,2 ) R Vlasov-Poisson can be written as (0,2 t f { H f (f ), f } = 0, with Poisson bracket {f, g} x v v x g Equivalent to d 8 G, G(f )=[H, G](f ) dt with Poisson bracket for functionals Z H [H, G] = f (f ){ G (f ), f }dxdv = f (0,2 ) R [G, H]. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 4 / 40
5 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) Relation between T and U We have the relation (in 1D 1D) [[T, U], U] =2U This implies RKN condition [[[T, U], U], U] =0. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 5 / 40
6 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) Derivation of high order methods We write p = e b1 U e a1 T e b 2 U e b2 U e a1 T e b 1 U Using for example the Baker-Campbell-Haussdorff formula, we get where (T +U)+R( ) p = e R( ) = 2 p 21 [T, U]+ 3 (p 31 [[T, U], T ]+p 32 [[T, U], U]) + 4 (p 41 [[[T, U], T ], T ]+p 42 [[[T, U], U], T ] +p 43 [[[T, U], U], U]) + O( 5 ) p ij are polynomials in the parameters a i, b i. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 6 / 40
7 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) Use of modified potentials We can incorporate the flow associated to [[T, U], U]. p = e b 1 U+c 1 3 [[T,U],U] e a 1 T e a1 T e b 1 U+c 1 3 [[T,U],U] ) Achievement of given order with less stages further reduction of e Ci, when [[T, U], U] =2U, with C i = b i U + c i 2 [[T, U], U]+d i 4 W 5,1 + e i 6 W 7,1 = (b i + 2c i 2 + 4d i 4 8e i 6 )U with W 5,1 = [U, [U, [T, [T, U]]]] W 7,1 = [U, [T, [U, [U, [T, [T, U]]]]]]. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 7 / 40
8 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) 6-th order time splitting s = 11, a 1 = t 2 ( ), a 2 = , a 3 = t 2 ( ) +4 t 4 ( ), a 4 = , a 5 = t 2 ( ) +4 t 4 ( ) 8 t 6 ( ), a 6 = , together with a 6+i = a 6 i, i = 1,...,5 and begin_with_x = false. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 8 / 40
9 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) Application to a KEEN wave test problem FIGURE: Comparison of the time evolution of the first 5 harmonics t = (Strang) comparable to t = 0.5 (Sixth-order scheme). gain of factor 10 at least M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM 14 9 / 40
10 1. Sixth order splitting scheme for Vlasov-Poisson (3-10) Application to Landau damping test problem t = 2 M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
11 2. 2D conservative Semi-Lagrangian scheme (11-13) CSL2D : Conservative Semi-Lagrangian Vlasov solver Typical 2D advection block, here in @ t r g g = 0. r r Unknowns are mean values ḡ i,j = 1 A i,j RA i,j g h Degree 2 reconstruction of g h with finite difference of order p for edge and point values Update of mean values by intersection of backward volume f method for radial boundary conditions M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
12 2. 2D conservative Semi-Lagrangian scheme (11-13) CSL2D : diocotron instability test case FIGURE: Time evolution of relative total mass (left) and of energy (right) M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
13 2. 2D conservative Semi-Lagrangian scheme (11-13) CSL2D : drift kinetic simulation Better mass (left) and energy (right) conservation w.r.t classical BSL2D (Backward semi-lagrangian) Here on grid, with t = 8. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
14 3. Curvilinear and hexagonal mesh (14-17) Semi-Lagrangian scheme on curvilinear grids FIGURE: Representation of a mapped mesh in two dimensions M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
15 3. Curvilinear and hexagonal mesh (14-17) Guiding center simulation on curvilinear mesh Examples of meshes : Colella mesh, D-shaped mesh M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
16 3. Curvilinear and hexagonal mesh (14-17) Semi-Lagrangian schemes on hexagonal mesh r 3 r 2 r 1 FIGURE: Representation of an hexagonal mesh M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
17 3. Curvilinear and hexagonal mesh (14-17) Circular advection test problem M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
18 4. Gyroaverage operator on polar grids (18-33) Gyroaverage operator definition The gyroaverage J (f ) of a polar fonction f is defined as J (f )(r, )= 1 Z 2 f (~x 2 G + ~ )d. 0 where ~x G = r(cos( ), sin( )) and ~ = (cos( ), sin( )). This gyroaverage process consists in computing an average on the Larmor circle. It tends to damp any fluctuation which develops at sub-larmor scales. Let consider a uniform polar mesh (r, ) 2 [r min, r max ] [0, 2 [ with N r N cells, our goal is to approximate the operator f j,k 2 R (Nr +1) N 7! J (f ) j,k 2 R (Nr +1) N. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
19 4. Gyroaverage operator on polar grids (18-33) Method based on Padé approximation FIGURE: The zero-th order Bessel function J 0 (k ) compare to its Padé approximation 1/ 1 +(k ) 2 /4. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
20 4. Gyroaverage operator on polar grids (18-33) Method based on interpolation M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
21 4. Gyroaverage operator on polar grids (18-33) Benchmark with the classical 5D Cyclone DIII-D case FIGURE: Linear growth rate versus k i for the CYCLONE DIII-D base case : Padé vs interpolation M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
22 4. Gyroaverage operator on polar grids (18-33) Derivation of the quasi-neutrality equation Electrons follow a Boltzmann distribution : ( h i) n e = n 0 exp n 0 1 T e 1 ( h i). T e The ion density n i can be written at the first gyrokinetic order : Z n i (x) = J (f + g)(x, v z ) d dv z, = p 2µ with the correction : g(x, ):=@ µ F M (r, )[ (x) J ( )(x)]. Under the assumption of quasi-neutral limit n i = n e : Z n 0 ( J 2 n p 0 ( )) exp( µ/t T i R + 2µ i )dµ ( h i) = n i n 0. T e M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
23 4. Gyroaverage operator on polar grids (18-33) QN solver using Padé := 1 Z Jp T i ZR 2 ( )e µ/t i dµ = Jp 2 ( )e µ dµ. + 2µ R + 2µTi In Fourier space : \J ( )(k) =J 0 ( k ) b (k) ) b (k) = 0 ( k 2 T i ) b (k) where the function 0 is defined by : Z 0(k 2 ):= exp( R + A development of 0 leads to x 2 /2)J 2 0 (kx)xdx. r? (n 0 r? ) + n 0 T e ( h i) = n i n 0. This equation can be solved by FFT in and finite differences in r. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
24 4. Gyroaverage operator on polar grids (18-33) QN solver using Padé in a simplified model We also consider a simplified case of the quasi-neutrality equation where the integrals in µ are replaced by the evaluation at one value µ = µ 0 : n 0 T i ( J 2 p 2µ0 T i ( )) + n 0 T e ( h i) = In Fourier space, \ [ J 2 p 2µ0 T i ( )](k) = Z J p 2µ 0 (f )(x, v z )dv z n 0. 1 J 2 0 ( k p 2µ 0 T i ) b(k). The development : 1 J0 2 ( k p 2µ 0 T i )= k 2 µ 0 T i + O( k 4 ). leads to the equation : µ 0? + 1 ( h i) = 1 Z J T e n p (f )(x, v 2µ z )dv z 0 0 n 0. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
25 4. Gyroaverage operator on polar grids (18-33) QN solver using interpolation Matrix formulation of the linear operator : 7! n Z 0 ( Jp 2 ( )) exp( µ/t T i R + 2µ i )dµ. Use of the method with interpolation for the computation of the gyroaverage Quadrature in µ using Gauss-Legendre method Precomputation in Fourier space. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
26 4. Gyroaverage operator on polar grids (18-33) Numerical comparison with analytical solutions Homogeneous Dirichlet conditions on r min and r max. m,` (r, )= J m ( m,`)y m r Y m ( m,`)j m r r max where m,` is the `th zero of y 7! J m (y)y m y r min verifies Then : Z +1 0 r max p2µ J p 2µ ( )(r m,` 0, 0 )=J 0 r max ( J 2 p 2µ ( ))e µ dµ + = m,` r max e im Y m (y)j m y r min r max (r 0, 0 ). 2 m,` rmax 2!!. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
27 4. Gyroaverage operator on polar grids (18-33) Numerical comparison with analytical solutions FIGURE: Comparison between the analytical solution of the quasi-neutrality equation (in red) and the solution obtained with Padé (left - in green) and with the interpolation method (right - in green) for ` = 1, m = 1. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
28 4. Gyroaverage operator on polar grids (18-33) Numerical comparison with analytical solutions FIGURE: Comparison between the analytical solution of the quasi-neutrality equation (in red) and the solution obtained with Padé (left - in green) and with the interpolation method (right - in green) for ` = 1, m = 5. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
29 4. Gyroaverage operator on polar grids (18-33) Application to gyrokinetic simulations The distribution function f (t, r,,z, v) is solution of J 2µ r J p! 0 2µ t r f f + r r v@ z z J 2µ0 v f = 0. for (r,,z, v) 2 [r min, r max ] [0, 2 ] [0, L] [ v max, v max ]. Coupling with : the quasi-neutrality equation solved with the interpolation method : 1 T i ( J 2 p 2µ0 T i ( ))+ 1 T e ( h i) = 1 n 0 Z J p 2µ 0 (f )(x, v z )dv z or its Padé approximation solved by FFT in and finite diff. in r : µ 0? + 1 ( h i) = 1 Z J T e n p (f )(x, v 2µ z )dv z n n 0 M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
30 4. Gyroaverage operator on polar grids (18-33) Numerical results FIGURE: f (r,,z = 0, v k = 0), , t = 2, interpolation (left), Padé (right), top : µ = 0.1 and T = 3000, bottom : µ = 0.2 and T = M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
31 4. Gyroaverage operator on polar grids (18-33) Numerical results FIGURE: f (r,,z = 0, v k = 0), , t = 2, µ = 0.7, T = 5000 (top), T = 7000 (bottom), with interpolation method (left) and Padé (right). M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
32 4. Gyroaverage operator on polar grids (18-33) Numerical results F IGURE: f (r,, z = 0, vk = 0), t = 2, µ = 1, T = 7000, interpolation (top), Padé (bottom). Mesh (from left to right) : ; ; and M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
33 4. Gyroaverage operator on polar grids (18-33) Numerical results FIGURE: Time evolution of R r max R 2 r min (r,,0)rdrd, mesh : , 0 comparison Padé/interpolation for different µ (0, 0.1, 0.2, 0.5, 0.7, 1). M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
34 5. Field aligned interpolation (34-37) Field aligned drift kinetic model Vlasov equation reads using v = v k for shorter @ r f + + v f + vb z f rb 0 rb 0 r r + b v f = 0, together with quasi-neutral equation (without zonal flow) r 2 + r r n r + 1 n 0 r + 1Te = 1n0 Z f f eq dv. = b /(2 r) b z /L, L = z max z min = 2 R. M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
35 Numerical method 5. Field aligned interpolation (34-37) Oblic constant advection with oblic interpolation constant advection in on a stencil around feet of characteristic for a given z = z k Lagrange interpolation of degree 2d + 1 in field aligned direction M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
36 5. Field aligned interpolation (34-37) Problem of visualization vs for representing the mode (30, 23) M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
37 5. Field aligned interpolation (34-37) = 0.8, m = 30, n = 23, T = 4000, 6000, 8000 M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
38 Conclusion Conclusion/perspectives/references (38-40) High order splitting scheme for Vlasov Poisson can outperform Strang splitting scheme Better energy (and mass) conservation property for conservative semi-lagrangian scheme Curvilinear grid for dealing with different geometries is an option Hexagonal grid is promising Interpolation method for gyroaverage operator leads to better accuracy First comparisons for quasi neutrality equation Validity of field aligned interpolation approach M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
39 Perspectives Conclusion/perspectives/references (38-40) Robust and efficient numerical methods for the poloidal plane Gyroaverage New Padé approximants Comprehension of continuous model and properties field aligned interpolation Test with shear safety factor Extension to 2D Hermite Explore "spaghetti" interpolation M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
40 References Conclusion/perspectives/references (38-40) Selalib : semi-lagrangian library. Afeyan, Casas, Crouseilles, Dodhy, Faou, M., Sonnendrücker, EPJD Vlasovia Steiner, M., Crouseilles, Grandgirard, Latu, Rozar, EPJD Vlasovia Crouseilles, Glanc, Hirstoaga, Madaule, M., Pétri, EPJD Vlasovia Mendoza, Prouveur, M., Sonnendrücker, Selhex project, CEMRACS Hamiaz, M., Back, Inria report, Hariri-Ottaviani, CPC M. Mehrenberger (UDS) Semi-Lagrangian (gyro)kinetic simulations Equations cinétiques, CIRM / 40
Advective and conservative semi-lagrangian schemes on uniform and non-uniform grids
Advective and conservative semi-lagrangian schemes on uniform and non-uniform grids M. Mehrenberger Université de Strasbourg and Max-Planck Institut für Plasmaphysik 5 September 2013 M. Mehrenberger (UDS
More informationConservative Semi-Lagrangian solvers on mapped meshes
Conservative Semi-Lagrangian solvers on mapped meshes M. Mehrenberger, M. Bergot, V. Grandgirard, G. Latu, H. Sellama, E. Sonnendrücker Université de Strasbourg, INRIA Grand Est, Projet CALVI, ANR GYPSI
More informationNew formulations of the semi-lagrangian method for Vlasov-type equations
New formulations of the semi-lagrangian method for Vlasov-type equations Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 17 September 2008 In collaboration
More informationSimulation gyrocinétique des plasmas de fusion magnétique
Simulation gyrocinétique des plasmas de fusion magnétique IRMA, Université de Strasbourg & CNRS Calvi project-team INRIA Nancy Grand Est Collaborators: J.P. Braeunig, N. Crouseilles, M. Mehrenberger, T.
More informationCONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS
CONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS CHANG YANG AND FRANCIS FILBET Abstract. We develop weighted essentially
More informationSolving the guiding-center model on a regular hexagonal mesh
Solving the guiding-center model on a regular hexagonal mesh Michel Mehrenberger, Laura S. Mendoza, Charles Prouveur, Eric Sonnendrücker To cite this version: Michel Mehrenberger, Laura S. Mendoza, Charles
More informationCONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS
CONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS CHANG YANG AND FRANCIS FILBET Abstract. We introduce a WENO reconstruction
More informationA Block-Based Parallel Adaptive Scheme for Solving the 4D Vlasov Equation
A Block-Based Parallel Adaptive Scheme for Solving the 4D Vlasov Equation Olivier Hoenen and Eric Violard Laboratoire LSIIT - ICPS, UMR CNRS 7005 Université Louis Pasteur, Strasbourg {hoenen,violard}@icps.u-strasbg.fr
More informationAsynchronous OpenCL/MPI numerical simulations of conservation laws
Asynchronous OpenCL/MPI numerical simulations of conservation laws Philippe HELLUY 1,3, Thomas STRUB 2. 1 IRMA, Université de Strasbourg, 2 AxesSim, 3 Inria Tonus, France IWOCL 2015, Stanford Conservation
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationMichel Mehrenberger 1, Laura S. Mendoza 2, 3, Charles Prouveur 4 and Eric Sonnendrücker 2,3
ESAIM: PROCEEDINGS AND SURVEYS, March 2016, Vol. 53, p. 149-176 M. Campos Pinto and F. Charles, Editors SOLVING THE GUIDING-CENTER MODEL ON A REGULAR HEXAGONAL MESH Michel Mehrenberger 1, Laura S. Mendoza
More informationApplication of High Dimensional B-Spline Interpolation in Solving the Gyro-Kinetic Vlasov Equation Based on Semi-Lagrangian Method
Commun Comput Phys doi: 1428/cicpOA-216-92 Vol 22, No 3, pp 789-82 September 217 Application of High Dimensional B-Spline Interpolation in Solving the Gyro-Kinetic Vlasov Equation Based on Semi-Lagrangian
More informationVerification and Validation for Seismic Wave Propagation Problems
Chapter 26 Verification and Validation for Seismic Wave Propagation Problems (1989-2-24-25-28-29-21-211-217-) (In collaboration with Dr. Nima Tafazzoli, Dr. Federico Pisanò, Mr. Kohei Watanabe and Mr.
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationLong time integrations of a convective PDE on the sphere by RBF collocation
Long time integrations of a convective PDE on the sphere by RBF collocation Bengt Fornberg and Natasha Flyer University of Colorado NCAR Department of Applied Mathematics Institute for Mathematics Applied
More informationAlmost Curvature Continuous Fitting of B-Spline Surfaces
Journal for Geometry and Graphics Volume 2 (1998), No. 1, 33 43 Almost Curvature Continuous Fitting of B-Spline Surfaces Márta Szilvási-Nagy Department of Geometry, Mathematical Institute, Technical University
More informationA mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates
A mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates Peter Hjort Lauritzen Atmospheric Modeling & Predictability Section National Center for Atmospheric
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More information= f (a, b) + (hf x + kf y ) (a,b) +
Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals
More informationConservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow. Abstract
Conservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow Jing-Mei Qiu 1 and Chi-Wang Shu Abstract In this paper, we propose a semi-lagrangian finite difference
More informationThin Film Simulation on a Rotating Wafer. B. Gschaider, D. Prieling, H. Steiner, P. Vita
Thin Film Simulation on a Rotating Wafer B. Gschaider, D. Prieling, H. Steiner, P. Vita Topics Motivation Finite Area Method Thin Film Model Impinging Jet Polydual Mesh Comparison with 3D Solution Conclusion
More informationVOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL
VOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL The following is a description of the model setups and input/output parameters for benchmarking analytical volcanic deformation models
More informationATM 298, Spring 2013 Lecture 4 Numerical Methods: Horizontal DiscreDzaDons April 10, Paul A. Ullrich (HH 251)
ATM 298, Spring 2013 Lecture 4 Numerical Methods: Horizontal DiscreDzaDons April 10, 2013 Paul A. Ullrich (HH 251) paullrich@ucdavis.edu Outline 1. Introduction / Motivation 2. Finite Difference Methods
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More informationStress due to surface load
Stress due to surface load To analyze problems such as compressibility of soils, bearing capacity of foundations, stability of embankments, and lateral pressure on earth retaining structures, we need to
More informationSPH: Why and what for?
SPH: Why and what for? 4 th SPHERIC training day David Le Touzé, Fluid Mechanics Laboratory, Ecole Centrale de Nantes / CNRS SPH What for and why? How it works? Why not for everything? Duality of SPH SPH
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationPROPERTIES OF NATURAL ELEMENT COORDINATES ON ANY POLYHEDRON
PROPRTIS OF NATURAL LMNT COORDINATS ON ANY POLYHDRON P. Milbradt and T. Fröbel Institute of Computer Science in Civil ngineering, Univercity of Hanover, 3067, Hanover, Germany; PH (+49) 5-76-5757; FAX
More informationApplication of optimal sampling lattices on CT image reconstruction and segmentation or three dimensional printing
Application of optimal sampling lattices on CT image reconstruction and segmentation or three dimensional printing XIQIANG ZHENG Division of Health and Natural Sciences, Voorhees College, Denmark, SC 29042
More informationIntroduction to Computer Graphics. Modeling (3) April 27, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (3) April 27, 2017 Kenshi Takayama Solid modeling 2 Solid models Thin shapes represented by single polygons Unorientable Clear definition of inside & outside
More informationA parametric ROM for the linear frequency domain approach to time-accurate CFD
Platzhalter für Bild, Bild auf itelfolie hinter das Logo einsetzen A parametric ROM for the linear frequency domain approach to time-accurate CFD Ralf Zimmermann, AG Numerik, U-BS Model Reduction of Complex
More informationHigh-Order Finite Difference Schemes for computational MHD
High-Order Finite Difference Schemes for computational MHD A. Mignone 1, P. Tzeferacos 1 and G. Bodo 2 [1] Dipartimento di Fisica Generale, Turin University, ITALY [2] INAF Astronomic Observatory of Turin,,
More informationParallelization of an advection-diffusion problem arising in edge plasma physics using hybrid MPI/OpenMP programming
Parallelization of an advection-diffusion problem arising in edge plasma physics using hybrid MPI/OpenMP programming Matthieu Kuhn 1, Guillaume Latu 2, Nicolas Crouseilles 3, and Stéphane Genaud 1 1 ICube,
More informationOverview of Traditional Surface Tracking Methods
Liquid Simulation With Mesh-Based Surface Tracking Overview of Traditional Surface Tracking Methods Matthias Müller Introduction Research lead of NVIDIA PhysX team PhysX GPU acc. Game physics engine www.nvidia.com\physx
More informationMeshless Modeling, Animating, and Simulating Point-Based Geometry
Meshless Modeling, Animating, and Simulating Point-Based Geometry Xiaohu Guo SUNY @ Stony Brook Email: xguo@cs.sunysb.edu http://www.cs.sunysb.edu/~xguo Graphics Primitives - Points The emergence of points
More informationSAMCEF for ROTORS. Chapter 3.2: Rotor modeling. This document is the property of SAMTECH S.A. MEF A, Page 1
SAMCEF for ROTORS Chapter 3.2: Rotor modeling This document is the property of SAMTECH S.A. MEF 101-03-2-A, Page 1 Table of contents Introduction Introduction 1D Model 2D Model 3D Model 1D Models: Beam-Spring-
More informationModeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures
University of California, Davis Modeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures N. Sukumar UC Davis Rutgers University,
More informationRegional Cooperation for Limited Area Modeling in Central Europe. Dynamics in LACE. Petra Smolíková thanks to many colleagues
Dynamics in LACE Petra Smolíková thanks to many colleagues Outline 1 1. VFE NH Design of vertical finite elements scheme for NH version of the model 2. SL scheme Application of ENO technique to SL interpolations
More informationMathematical Tools in Computer Graphics with C# Implementations Table of Contents
Mathematical Tools in Computer Graphics with C# Implementations by Hardy Alexandre, Willi-Hans Steeb, World Scientific Publishing Company, Incorporated, 2008 Table of Contents List of Figures Notation
More informationNumerical Aspects of Special Functions
Numerical Aspects of Special Functions Nico M. Temme In collaboration with Amparo Gil and Javier Segura, Santander, Spain. Nico.Temme@cwi.nl Centrum voor Wiskunde en Informatica (CWI), Amsterdam Numerics
More informationCharacteristic Aspects of SPH Solutions
Characteristic Aspects of SPH Solutions for Free Surface Problems: Source and Possible Treatment of High Frequency Numerical Oscillations of Local Loads. A. Colagrossi*, D. Le Touzé & G.Colicchio* *INSEAN
More informationParticle Tracing Module
Particle Tracing Module Particle Tracing Module Released with version 4.2a in October 2011 Add-on to COMSOL Multiphysics Combines with any COMSOL Multiphysics Module Particle Tracing Particle tracing
More informationFlow Structures of Jupiter s Great Red Spot Extracted by Using Optical Flow Method
Flow Structures of Jupiter s Great Red Spot Extracted by Using Optical Flow Method Tianshu Liu & Bo Wang Department of Mechanical & Aeronautical Engineering Western Michigan University, Kalamazoo, MI 49008
More informationFlow Structures Extracted from Visualization Images: Vector Fields and Topology
Flow Structures Extracted from Visualization Images: Vector Fields and Topology Tianshu Liu Department of Mechanical & Aerospace Engineering Western Michigan University, Kalamazoo, MI 49008, USA We live
More informationChallenges in adapting Particle-In-Cell codes to GPUs and many-core platforms
Challenges in adapting Particle-In-Cell codes to GPUs and many-core platforms L. Villard, T.M. Tran, F. Hariri *, E. Lanti, N. Ohana, S. Brunner Swiss Plasma Center, EPFL, Lausanne A. Jocksch, C. Gheller
More informationPlane wave solutions for the radiative transport equation
Plane wave solutions for the radiative transport equation Arnold D. Kim adkim@ucmerced.edu School of Natural Sciences University of California, Merced Overview Introduction to plane wave solutions Green
More informationAnimation, Teeth and Skeletons...
Animation, Teeth and Skeletons... (using Discrete Differential Geometry!) June 7, 2010 Animation Animation 3-dimensional Animation relies on computer graphics Animation 3-dimensional Animation relies on
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationEfficiency of adaptive mesh algorithms
Efficiency of adaptive mesh algorithms 23.11.2012 Jörn Behrens KlimaCampus, Universität Hamburg http://www.katrina.noaa.gov/satellite/images/katrina-08-28-2005-1545z.jpg Model for adaptive efficiency 10
More informationDirections: 1) Delete this text box 2) Insert desired picture here
Directions: 1) Delete this text box 2) Insert desired picture here Multi-Disciplinary Applications using Overset Grid Technology in STAR-CCM+ CD-adapco Dmitry Pinaev, Frank Schäfer, Eberhard Schreck Outline
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 58 63 (2007) Published online 22 May 2007 in Wiley InterScience (www.interscience.wiley.com).151 An explicit and conservative remapping strategy for semi-lagrangian
More informationKnow it. Control points. B Spline surfaces. Implicit surfaces
Know it 15 B Spline Cur 14 13 12 11 Parametric curves Catmull clark subdivision Parametric surfaces Interpolating curves 10 9 8 7 6 5 4 3 2 Control points B Spline surfaces Implicit surfaces Bezier surfaces
More informationQUIZ 4 (CHAPTER 17) SOLUTIONS MATH 252 FALL 2008 KUNIYUKI SCORED OUT OF 125 POINTS MULTIPLIED BY % POSSIBLE
QUIZ 4 (CHAPTER 17) SOLUTIONS MATH 5 FALL 8 KUNIYUKI SCORED OUT OF 15 POINTS MULTIPLIED BY.84 15% POSSIBLE 1) Reverse the order of integration, and evaluate the resulting double integral: 16 y dx dy. Give
More informationSLICE-3D: A three-dimensional conservative semi-lagrangian scheme for transport problems. Mohamed Zerroukat. Met Office, Exeter, UK
SLICE-D: A three-dimensional conservative semi-lagrangian scheme for transport problems Mohamed Zerroukat (Nigel wood & Andrew Staniforth) Met Office, Exeter, UK SRNWP-NT, Zagreb, c Crown Copyright Content
More informationFFD, mesh morphing & reduced order models: Enablers for efficient aerodynamic shape optimization
FFD, mesh morphing & reduced order models: Enablers for efficient aerodynamic shape optimization F Salmoiraghi, G Rozza SISSA mathlab A Scardigli, H Telib OPTIMAD engineering srl Outline 1. Prac'cal)problems)in)shape)op'miza'on)
More informationMetafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction
ALE simulations ua sus using Metafor eao 1. Introduction 2. Operator split 3. Convection schemes 4. Rezoning methods 5. Contact with friction 1 Introduction EULERIAN FORMALISM Undistorted mesh Ideal for
More informationAn added mass partitioned algorithm for rigid bodies and incompressible flows
An added mass partitioned algorithm for rigid bodies and incompressible flows Jeff Banks Rensselaer Polytechnic Institute Overset Grid Symposium Mukilteo, WA October 19, 216 Collaborators Bill Henshaw,
More informationCHARMS: A Simple Framework for Adaptive Simulation SIGGRAPH Presented by Jose Guerra
CHARMS: A Simple Framework for Adaptive Simulation SIGGRAPH 2002 Eitan Grinspun Caltech Petr Krysl UCSD Peter Schröder Caltech Presented by Jose Guerra 1 Outline Background Motivation (Element vs. Basis
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationIntroduction to Multigrid and its Parallelization
Introduction to Multigrid and its Parallelization! Thomas D. Economon Lecture 14a May 28, 2014 Announcements 2 HW 1 & 2 have been returned. Any questions? Final projects are due June 11, 5 pm. If you are
More informationDrop Impact Simulation with a Velocity-Dependent Contact Angle
ILASS Americas 19th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Drop Impact Simulation with a Velocity-Dependent Contact Angle S. Afkhami and M. Bussmann Department
More informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
More informationDiscrete representations of geometric objects: Features, data structures and adequacy for dynamic simulation. Part I : Solid geometry
Discrete representations of geometric objects: Features, data structures and adequacy for dynamic simulation. Surfaces Part I : Solid geometry hachar Fleishman Tel Aviv University David Levin Claudio T.
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 28 / 46
Polar Coordinates Polar Coordinates: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ)
More informationAdarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs)
OBJECTIVE FLUID SIMULATIONS Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs) The basic objective of the project is the implementation of the paper Stable Fluids (Jos Stam, SIGGRAPH 99). The final
More informationName: Class: Date: 1. Use Lagrange multipliers to find the maximum and minimum values of the function subject to the given constraint.
. Use Lagrange multipliers to find the maximum and minimum values of the function subject to the given constraint. f (x, y) = x y, x + y = 8. Set up the triple integral of an arbitrary continuous function
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationLS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.
12 th International LS-DYNA Users Conference FSI/ALE(1) LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA Part 1 Facundo Del
More informationMath 225 Scientific Computing II Outline of Lectures
Math 225 Scientific Computing II Outline of Lectures Spring Semester 2003 I. Interpolating polynomials Lagrange formulation of interpolating polynomial Uniqueness of interpolating polynomial of degree
More informationIntroduction to Computer Graphics. Image Processing (1) June 8, 2017 Kenshi Takayama
Introduction to Computer Graphics Image Processing (1) June 8, 2017 Kenshi Takayama Today s topics Edge-aware image processing Gradient-domain image processing 2 Image smoothing using Gaussian Filter Smoothness
More informationMulti-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech
Multi-Mesh CFD Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech cjroy@vt.edu May 21, 2014 CCAS Program Review, Columbus, OH 1 Motivation Automated
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 27 / 45
: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ) Chapter 10: Parametric Equations
More informationGEOMETRIC TOOLS FOR COMPUTER GRAPHICS
GEOMETRIC TOOLS FOR COMPUTER GRAPHICS PHILIP J. SCHNEIDER DAVID H. EBERLY MORGAN KAUFMANN PUBLISHERS A N I M P R I N T O F E L S E V I E R S C I E N C E A M S T E R D A M B O S T O N L O N D O N N E W
More information3D Finite Element Software for Cracks. Version 3.2. Benchmarks and Validation
3D Finite Element Software for Cracks Version 3.2 Benchmarks and Validation October 217 1965 57 th Court North, Suite 1 Boulder, CO 831 Main: (33) 415-1475 www.questintegrity.com http://www.questintegrity.com/software-products/feacrack
More informationDesign Optimization of a Weather Radar Antenna using Finite Element Analysis (FEA) and Computational Fluid Dynamics (CFD)
Design Optimization of a Weather Radar Antenna using Finite Element Analysis (FEA) and Computational Fluid Dynamics (CFD) Fernando Prevedello Regis Ataídes Nícolas Spogis Wagner Ortega Guedes Fabiano Armellini
More informationExploring unstructured Poisson solvers for FDS
Exploring unstructured Poisson solvers for FDS Dr. Susanne Kilian hhpberlin - Ingenieure für Brandschutz 10245 Berlin - Germany Agenda 1 Discretization of Poisson- Löser 2 Solvers for 3 Numerical Tests
More informationAuto Injector Syringe. A Fluent Dynamic Mesh 1DOF Tutorial
Auto Injector Syringe A Fluent Dynamic Mesh 1DOF Tutorial 1 2015 ANSYS, Inc. June 26, 2015 Prerequisites This tutorial is written with the assumption that You have attended the Introduction to ANSYS Fluent
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationOverview and Recent Developments of Dynamic Mesh Capabilities
Overview and Recent Developments of Dynamic Mesh Capabilities Henrik Rusche and Hrvoje Jasak h.rusche@wikki-gmbh.de and h.jasak@wikki.co.uk Wikki Gmbh, Germany Wikki Ltd, United Kingdom 6th OpenFOAM Workshop,
More informationMath background. 2D Geometric Transformations. Implicit representations. Explicit representations. Read: CS 4620 Lecture 6
Math background 2D Geometric Transformations CS 4620 Lecture 6 Read: Chapter 2: Miscellaneous Math Chapter 5: Linear Algebra Notation for sets, functions, mappings Linear transformations Matrices Matrix-vector
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationAudio acoustic modeling using full-wave methods
Acoustics 8 Paris Audio acoustic modeling using full-wave methods Timo Lahivaara, Tomi Huttunen and Simo-Pekka Simonaho University of Kuopio, P.O.Box 1627, 7211 Kuopio, Finland simo-pekka.simonaho@uku.fi
More informationA Comparison of Some Numerical Methods for the Advection-Diffusion Equation
Res Lett Inf Math Sci, 26, Vol1, pp49-62 Available online at http://iimsmasseyacnz/research/letters/ 49 A Comparison of Some Numerical Methods for the Advection-Diffusion Equation M Thongmoon 1 & R McKibbin
More informationA brief description of the particle finite element method (PFEM2). Extensions to free surface
A brief description of the particle finite element method (PFEM2). Extensions to free surface flows. Juan M. Gimenez, L.M. González, CIMEC Universidad Nacional del Litoral (UNL) Santa Fe, Argentina Universidad
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationParallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes Stefan Vater 1 Kaveh Rahnema 2 Jörn Behrens 1 Michael Bader 2 1 Universität Hamburg 2014 PDES Workshop 2 TU München Partial
More informationAcknowledgment. I would like to thank several people who have helped me to complete this
Acknowledgment I would like to thank several people who have helped me to complete this thesis. First of all, I need to mention my advisor, Professor Chris Lacor, who has given me the opportunity to do
More informationA complete methodology for the implementation of XFEM models. Team: - Carlos Hernán Villanueva - Kai Yu
A complete methodology for the implementation of XFEM models Team: - Carlos Hernán Villanueva - Kai Yu Problem Statement QUAD4 DISCONTINUITY Problem Statement Traditional FEM Mesh generates discrete representation
More informationFederal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, Berlin, Germany
Jannis Bulling 1, Jens Prager 1, Fabian Krome 1 1 Federal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, 12205 Berlin, Germany Abstract: This paper addresses the computation of
More informationDiscrete Geometry Processing
Non Convex Boundary Convex boundary creates significant distortion Free boundary is better Some slides from the Mesh Parameterization Course (Siggraph Asia 008) 1 Fixed vs Free Boundary Fixed vs Free Boundary
More informationBubble Dynamics using Free Surfaces in a VOF framework
Bubble Dynamics using Free Surfaces in a VOF framework Sevilla Meeting: Numerical Challenges in two-phase flows 27 October 2014 Léon MALAN Supervisor: Prof. Stéphane ZALESKI 1 Contents Introduction: Simulation
More informationixcube 4-10 Brief introduction for membrane and cable systems.
ixcube 4-10 Brief introduction for membrane and cable systems. ixcube is the evolution of 20 years of R&D in the field of membrane structures so it takes a while to understand the basic features. You must
More informationA Brief Introduction to Ray Tracing and Ionospheric Models
ABriefIntroductiontoRayTracingand Ionospheric Models Stephen R. Kaeppler AD0AE skaeppl@clemson.edu Clemson University March 2, 2018 Motivation - why should you care? Prediction tool for understanding HF
More informationCS-184: Computer Graphics Lecture #21: Fluid Simulation II
CS-184: Computer Graphics Lecture #21: Fluid Simulation II Rahul Narain University of California, Berkeley Nov. 18 19, 2013 Grid-based fluid simulation Recap: Eulerian viewpoint Grid is fixed, fluid moves
More information