Numerical Simulation of Heat Transfer in Materials with Anisotropic Thermal Conductivity: A Finite Volume Scheme to Handle Complex Geometries
|
|
- Kory Trevor Shaw
- 5 years ago
- Views:
Transcription
1 Numerical Simulation of Heat Transfer in Materials with Anisotropic Thermal Conductivity: A Finite Volume Scheme to Handle Complex Geometries Olaf Klein, Jürgen Geiser, Peter Philip 2 Weierstrass Institute for Applied Analysis and Stochastics (WIAS) Berlin, Germany 2 University of Minnesota Institute for Mathematics and its Applications (IMA) Minneapolis, USA IMA Workshop: New Paradigms in Computation Minneapolis, March 28 30, 2005
2 Complex Sample Domain from Crystal Growth Blind hole (for cooling and measurements). Ω 3 Ω Ω Ω 2 Ω 2 Ω 4 Ω 5 Ω : Insulation (often anisotropic). Ω 2 : Graphite crucible. Ω 3 : SiC crystal seed. Ω 4 : Gas enclosure. Ω 5 : SiC powder source. Ω 6 : Quartz. Ω = 6 S m= Ω m. Ω 6 Figure : Axisymmetric domain representing a growth apparatus used in silicon carbide single crystal growth by physical vapor transport (PVT). The geometry is a modified version of K. Semmelroth et al., J. Phys.-Condes. Matter 6 (2004).
3 Model for Stationary Anisotropic Heat Conduction div(k m (θ) θ) = f m in Ω m (m M), θ: absolute temperature, K m : symmetric and positive definite tensor of thermal conductivity, f m : heat sources, Ω m : domain of material m. Assumed form of K m : K m (θ) = ą κ m i,j (θ)ć, where κ m i,j (θ) = 8 < α m i κ m iso (θ) for i = j, : 0 for i j. Interface Conditions on Ω m Ω m2 : ą Km (θ) θ ć Ωm n m = ą K m2 (θ) θ ć Ωm2 n m. : restriction, n m : outer unit normal vector to material m. Boundary Conditions: Dirichlet, Robin θ = θ Dir on Γ Dir, ą K m (θ) θ ć n m = ξ m (θ θ ext ) on Γ Rob Ω m, m M,
4 Finite Volume Discretization Σ m = (σ m,i ) i Im conforming triangulation of Ω m satisfying the constrained Delaunay property: If γ is an interior edge of Σ m, α and β the angles opposite to γ, then α + β π. If γ Ω m is a boundary edge of Σ m, α the angle opposite γ, then α π/2. V (σ m,i ) = ľ v m i,j : j {, 2, 3}ł : Set of vertices. V := S m M, i I m V (σ m,i ). ω v := ľ x Ω : x v 2 < x z 2 for each z V \ {v} ł, ω m,v := ω v Ω m, V m := {z V : ω m,z }. A m = (a m i,j ), am i,j := 8 < α m i for i = j, : 0 for i j. w ω,w σ = conv{v, w, u }, σ 2 = conv{v, w, u 2 } Ω = σ σ 2 u ω,u σ ω,u σ 2 γ,v,u σ γ,v,u σ 2 ω,v γ,v,w γ,w,u2 γ,v,u2 ω,u2 u 2 v Figure 2: Illustration of the space discretization.
5 Approximation of Anisotropic Terms φ σ,v : σ [0, ]: Affine coordinates on triangle σ w.r.t. v V (σ). For each edge [v, w] of some σ Σ m : Letting Σ m,v,w := ľ σ Σ m : {v, w} V (σ) ł. Σ γm,v,w := ľ σ Σ m,v,w : λ (H v,w,σ γ m,v,w ) 0 ł, decompose γ m,v,w : γ m,v,w = [ σ Σ γ m,v,w σ γ m,v,w. Approximation: (A m θ) σ n ωv γm,v,w X θ(ṽ) (A m φ σ,ṽ ) ṽ V (σ) w v w v 2. Finite Volume Scheme Find (θ v ) v V satisfying: θ v = θ Dir (v) for each v V Dir, 0 = X m M X m M X m M ξ m ą θv θ ext (v) ć λ ( ω m,v Γ Rob ) X σ Σ γ m,v,w X ṽ V (σ) 2 ą κ m iso (θ v ) + κ m iso (θ w) ć θṽ (A m φ σ,ṽ ) w v w v 2 λ (H v,w,σ γ m,v,w ) f m,v λ 2 (ω m,v ) for each v V Dir = V \ V Dir.
6 Comparison with Closed-Form Solution Axisymmetric Single-Material Domain Ω = {(r, z) : 0 < r < 0.2, 0.2 < z < 0.2}: ţ ű θ r α r ţ ű θ α z = 0 r r r z z in Ω, (a) θ Dir (r, z) := 2 Solution: θ(r, z) = 2 Numerical solution θ 0 num r 2 z 2 on Ω. (b) α r α z r 2 z 2 on Ω. α r α z Exact solution θ Exact Stationary Solution z = 20 T_min=-0.04 T_max=0.002 Levels: : 0.0 2: : : : z = -20 r = 0 r = 20 Figure 3: Solution of (): Numerical θnum, 0 37 triangles (left); exact θ (right). Isolevel difference: Discrete L -error: ɛ l L := X vol(ω v ) θnum l (v) θ(v), v V l v V l : vertices, vol(ω v ): r-weighted area of Voronoï cell. Numerical convergence rate: ρ l L := (ln(ɛ l L ) ln(ɛ l L ))/(ln(h l ) ln(h l )), h l : upper bound for triangle area of level l.
7 Axisymmetric Multi-Material Domain r r 00 Ω = {(r, z) : 0 < r < r 0, 0 < z < z 0 }, Ω 2 = {(r, z) : r 0 < r < r max, 0 < z < z 0 }, Ω 3 = {(r, z) : 0 < r < r 0, z 0 < z < z max }, Ω 4 = {(r, z) : r 0 < r < r max, z 0 < z < z max }, ţ θ r α m,r α m,r 0 0 α m,z 00 α m,r 0 0 α m,z Solution: ű z A θ Ωm ţ ű θ α m,z z A n m A θ Ω m = f m in Ω m, (2a) A n m on Ω m Ω m, (2b) θ Dir,m (r, z) := a m r 2 + b m z 2 + c m on Ω Ω m. (2c) θ(r, z) := a m r 2 + b m z 2 + c m on Ω m, θ Dir,m (r, z) := a m r 2 + b m z 2 + c m on Ω Ω m, where r 0 = z 0 = 0., r max = z max = 0.2, α,r = 2, α 2,r =, α 3,r = 4, α 4,r = 2, α,z =, α 2,z = 2, α 3,z = 3, α 4,z = 6, a =, a 2 = 2, a 3 =, a 4 = 2, b =, b 2 =, b 3 = /3, b 4 = /3, c = 0, c 2 = /00, c 3 = 2/300, c 4 = /300, f = 0, f 2 = 2.0, f 3 = 8.0, f 4 = 20.0
8 0 Numerical solution θnum Exact solution θ Levels: : Ω3 Ω Ω 2: : : : Ω2 0 Figure 4: Solution of (2): Numerical θnum, 37 triangles (left); exact θ (right). Isolevel difference: level number l of triangles max area hl L -error ²lL numerical convergence rate ρll Table : L -error and numerical convergence rate for the numerical solution of () with anisotropy (αr, αz ) = (0, ). level number l of triangles max area hl L -error ²lL numerical convergence rate ρll Table 2: L -error and numerical convergence rate for the numerical solution of (2).
9 Stationary heat field Stationary temperature field Heat source 580 K αr αz Maximal Temperature [K] 220 K 820 K Figure 5: Left: Domain of heat sources highlighted. Right: T -field for isotropic insulation, i.e. αr = αz =. Stationary temperature field Stationary temperature field 660 K Stationary temperature field 580 K 580 K 740 K 740 K 700 K 900 K 900 K 580 K Figure 6: T -field for anisotropic insulation with αr = 000 (left), with αr = 000 for sides, αz = 000 for top and bottom (middle), αz = 000 (right).
10 Publications P. PHILIP: Transient Numerical Simulation of Sublimation Growth of SiC Bulk Single Crystals. Modeling, Finite Volume Method, Results, Thesis, Department of Mathematics, Humboldt University of Berlin, Germany, 2003 Report No. 22, Weierstrass Institute for Applied Analysis and Stochastics, Berlin. J. GEISER, O. KLEIN, P. PHILIP: Numerical simulation of heat transfer in materials with anisotropic thermal conductivity: A finite volume scheme to handle complex geometries. In preparation. J. GEISER, O. KLEIN, P. PHILIP: Influence of anisotropic thermal conductivity in the apparatus insulation for sublimation growth of SiC: Numerical investigation of heat transfer. In preparation. Funding: Supported by the DFG Research Center Matheon: Mathematics for Key Technologies in Berlin, by the IMA in Minneapolis, and by the WIAS in Berlin.
Institute for Mathematics and its Applications
Institute for Mathematics and its Applications University of Minnesota 400 Lind Hall, 207 Church St. SE, Minneapolis, MN 55455 Numerical Simulation and Control of Sublimation Growth of SiC Bulk Single
More informationWIAS-HITNIHS: A SOFTWARE TOOL FOR SIMULATION IN SUBLIMATION GROWTH OF SIC SINGLE CRYSTALS: APPLICATIONS AND METHODS.
WIAS-HITNIHS: A SOFTWARE TOOL FOR SIMULATION IN SUBLIMATION GROWTH OF SIC SINGLE CRYSTALS: APPLICATIONS AND METHODS. JÜRGEN GEISER, OLAF KLEIN, AND PETER PHILIP Abstract. The numerous technical applications
More informationInstitute for Mathematics and its Applications
Institute for Mathematics and its Applications University of Minnesota 400 Lind Hall, 207 Church St. SE, Minneapolis, MN 55455 Numerical Simulation and Control of Sublimation Growth of SiC Bulk Single
More informationNumerical simulations of the influence of a traveling magnetic field, generated by an internal heater-magnet module, on Czochralski crystal growth
International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, October 7-9, 8 Numerical simulations of the influence of a traveling magnetic field, generated by an internal heater-magnet
More informationIMA Preprint Series # 2002
OPTIMAL CONTROL OF A SEMILINEAR PDE WITH NONLOCAL RADIATION INTERFACE CONDITIONS By C. Meyer P. Philip and F. Tröltzsch IMA Preprint Series # ( October 4 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Report ISSN
Weierstraß-Institut für Angewandte Analysis und Stochasti im Forschungsverbund Berlin e.v. Report ISSN 0946 8838 Transient Numerical Simulation of Sublimation Growth of SiC Bul Single Crystals Modeling,
More informationUniform exponential decay of the free energy for Voronoi finite volume discretized reaction-diffusion systems
Weierstrass Institute for Applied Analysis and Stochastics Special Session 36 Reaction Diffusion Systems 8th AIMS Conference on Dynamical Systems, Differential Equations & Applications Dresden University
More informationFinite element analysis of the temperature field in a vertical MOCVD reactor by induction heating
Vol. 30, No. 11 Journal of Semiconductors November 2009 Finite element analysis of the temperature field in a vertical MOCVD reactor by induction heating Li Zhiming( ), Xu Shengrui( ), Zhang Jincheng(
More informationETIKA V PROFESII PSYCHOLÓGA
P r a ž s k á v y s o k á š k o l a p s y c h o s o c i á l n í c h s t u d i í ETIKA V PROFESII PSYCHOLÓGA N a t á l i a S l o b o d n í k o v á v e d ú c i p r á c e : P h D r. M a r t i n S t r o u
More informationHp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme
Hp-Adaptivity on Anisotropic Meshes for Hybridized Discontinuous Galerkin Scheme Aravind Balan, Michael Woopen and Georg May AICES Graduate School, RWTH Aachen University, Germany 22nd AIAA Computational
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationMeridional Flow, Differential Rotation, and the Solar Dynamo
Meridional Flow, Differential Rotation, and the Solar Dynamo Manfred Küker 1 1 Leibniz Institut für Astrophysik Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany Abstract. Mean field models of rotating
More informationDielectric Boundary in Biomolecular Solvation Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego
Dielectric Boundary in Biomolecular Solvation Bo Li Department of Mathematics and Center for Theoretical Biological Physics UC San Diego International Conference on Free Boundary Problems Newton Institute,
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationThree remarks on anisotropic finite elements
Three remarks on anisotropic finite elements Thomas Apel Universität der Bundeswehr München Workshop Numerical Analysis for Singularly Perturbed Problems dedicated to the 60th birthday of Martin Stynes
More informationApproximating Riemannian Voronoi Diagrams and Delaunay triangulations.
Approximating Riemannian Voronoi Diagrams and Delaunay triangulations. Jean-Daniel Boissonnat, Mael Rouxel-Labbé and Mathijs Wintraecken Boissonnat, Rouxel-Labbé, Wintraecken Approximating Voronoi Diagrams
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 informationarxiv: v1 [cond-mat.other] 11 Sep 2008
arxiv:0809.1990v1 [cond-mat.other] 11 Sep 2008 Momentum deficit in quantum glasses A.F. Andreev Kapitza Institute for Physical Problems, Russian Academy of Sciences, Kosygin Str. 2, Moscow, 119334 Russia
More informationA Semi-Analytical Thermal Elastic Model for Directional Crystal Growth with Weakly Anisotropy
A Semi-Analytical Thermal Elastic Model for Directional Crystal Growth with Weakly Anisotropy School of Mathematical Sciences Peking University with Sean Bohun, Penn State University Huaxiong Huang, York
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationQueen s University at Kingston. Faculty of Arts and Science. Department of Physics PHYSICS 106. Final Examination.
Page 1 of 5 Queen s University at Kingston Faculty of Arts and Science Department of Physics PHYSICS 106 Final Examination April 16th, 2009 Professor: A. B. McLean Time allowed: 3 HOURS Instructions This
More informationString-representation of hysteresis operators acting on vector-valued, left-continuous and piecewise monotaffine and continuous functions
Weierstrass Institute for Applied Analysis and Stochastics String-representation of hysteresis operators acting on vector-valued, left-continuous and piecewise monotaffine and continuous functions Olaf
More informationarxiv: v2 [cond-mat.mes-hall] 9 Nov 2017
arxiv:1601.01962v2 [cond-mat.mes-hall] 9 Nov 2017 Zeszyty Naukowe WWSI, No 13, Vol. 9, 2015, pp. 103-126 Behavioral modeling of stressed MOSFET Zenon Gniazdowski * Warsaw School of Computer Science Abstract
More informationSupporting Online Material for
www.sciencemag.org/cgi/content/full/314/586/1757/dc1 Supporting Online Material for Quantum Spin Hall Effect and Topological Phase Transition in HgTe Quantum Wells B. Andrei Bernevig, Taylor L. Hughes,
More informationRemote Asymmetric Induction in an Intramolecular Ionic Diels-Alder Reaction: Application to the Total Synthesis of (+)-Dihydrocompactin
Page S16 Remote Asymmetric Induction in an Intramolecular Ionic Diels-Alder Reaction: Application to the Total Synthesis of (+)-Dihydrocompactin Tarek Sammakia,* Deidre M. Johns, Ganghyeok Kim, and Martin
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationCompare and Contrast electrical resistivity measurements in linear 4 probe and Montgomery
Compare and Contrast electrical resistivity measurements in linear 4 probe and Montgomery Rabindra Bhattarai Department of Physics Iowa State University What is Electrical Resistivity Electrical resistivity
More informationA note on scenario reduction for two-stage stochastic programs
A note on scenario reduction for two-stage stochastic programs Holger Heitsch a and Werner Römisch a a Humboldt-University Berlin, Institute of Mathematics, 199 Berlin, Germany We extend earlier work on
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 Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
More informationFrame Dragging Anomalies for Rotating Bodies
General Relativity and Gravitation, Vol. 36, No. 5, May 2004 ( C 2004) LETTER Frame Dragging Anomalies for Rotating Bodies Peter Collas 1 and David Klein 2 Received October 7, 2003 Examples of axially
More informationMathematical Modeling for a PEM Fuel Cell.
Mathematical Modeling for a PEM Fuel Cell. Mentor: Dr. Christopher Raymond, NJIT. Group Jutta Bikowski, Colorado State University Aranya Chakrabortty, RPI Kamyar Hazaveh, Georgia Institute of Technology
More informationAn Estimator for statistical anisotropy from the CMB. CMB bispectrum
An Estimator for statistical anisotropy from the CMB bispectrum 09/29/2012 1 2 3 4 5 6 ...based on: N. Bartolo, E. D., M. Liguori, S. Matarrese, A. Riotto JCAP 1201:029 N. Bartolo, E. D., S. Matarrese,
More informationA Space-Time Boundary Element Method for the Wave Equation
W I S S E N T E C H N I K L E I D E N S C H A F T A Space-Time Boundary Element Method for the Wave Equation Marco Zank and Olaf Steinbach Institut für Numerische Mathematik Space-Time Methods for PDEs,
More informationREGULARITY AND MULTIGRID ANALYSIS FOR LAPLACE-TYPE AXISYMMETRIC EQUATIONS. Hengguang Li. IMA Preprint Series #2417.
REGULARITY AND MULTIGRID ANALYSIS FOR LAPLACE-TYPE AXISYMMETRIC EQUATIONS By Hengguang Li IMA Preprint Series #417 February 013) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF MINNESOTA 400
More informationGetting The Spin Right In Black-Hole Binaries
Getting The Spin Right In Black-Hole Binaries Gregory B. Cook Wake Forest University July 25, 2005 Abstract We will take a detailed look at the issues involved in setting the spin of a black hole during
More informationCIMPA Summer School on Current Research in Finite Element Methods
CIMPA Summer School on Current Research in Finite Element Methods Local zooming techniques for the finite element method Alexei Lozinski Laboratoire de mathématiques de Besançon Université de Franche-Comté,
More informationBlack Hole-Neutron Star Binaries in General Relativity. Thomas Baumgarte Bowdoin College
Black Hole-Neutron Star Binaries in General Relativity Thomas Baumgarte Bowdoin College 1 Why do we care? Compact binaries (containing neutron stars and/or black holes) are promising sources of gravitational
More informationStructure Report for J. Reibenspies
X-ray Diffraction Laboratory Center for Chemical Characterization and Analysis Department of Chemistry Texas A & M University Structure Report for J. Reibenspies Project Name: Sucrose Date: January 29,
More informationChemistry 532 Practice Final Exam Fall 2012 Solutions
Chemistry 53 Practice Final Exam Fall Solutions x e ax dx π a 3/ ; π sin 3 xdx 4 3 π cos nx dx π; sin θ cos θ + K x n e ax dx n! a n+ ; r r r r ˆL h r ˆL z h i φ ˆL x i hsin φ + cot θ cos φ θ φ ) ˆLy i
More informationAnalysis of a high order trace finite element method for PDEs on level set surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 7 Analysis of a high order trace finite element method for PDEs on level set surfaces Jörg Grande *, Christoph Lehrenfeld and Arnold Reusken * Institut für Geometrie
More informationA Unified Theorem on SDP Rank Reduction. yyye
SDP Rank Reduction Yinyu Ye, EURO XXII 1 A Unified Theorem on SDP Rank Reduction Yinyu Ye Department of Management Science and Engineering and Institute of Computational and Mathematical Engineering Stanford
More informationMathematical Modeling of Displacements and Thermal Stresses in Anisotropic Materials (Sapphire) in Cooling
Mathematical Modeling of Displacements and Thermal Stresses in Anisotropic Materials (Sapphire) in Cooling Timo Tiihonen & Tero Tuovinen September 11, 2015 European Study Group with Industry, ESGI 112,
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationAnisotropic turbulence in rotating magnetoconvection
Anisotropic turbulence in rotating magnetoconvection André Giesecke Astrophysikalisches Institut Potsdam An der Sternwarte 16 14482 Potsdam MHD-Group seminar, 2006 André Giesecke (AIP) Anisotropic turbulence
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
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 informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationBachelor Thesis. General Relativity: An alternative derivation of the Kruskal-Schwarzschild solution
Bachelor Thesis General Relativity: An alternative derivation of the Kruskal-Schwarzschild solution Author: Samo Jordan Supervisor: Prof. Markus Heusler Institute of Theoretical Physics, University of
More informationSmall Molecule Crystallography Lab Department of Chemistry and Biochemistry University of Oklahoma 101 Stephenson Parkway Norman, OK
Small Molecule Crystallography Lab Department of Chemistry and Biochemistry University of Oklahoma 101 Stephenson Parkway Norman, OK 73019-5251 Sample: KP-XI-cinnamyl-chiral alcohol Lab ID: 12040 User:
More informationA posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model
A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP
More informationThe QHE as laboratory system to study quantum phase transitions
The QHE as laboratory system to study quantum phase transitions Anne de Visser, WZI-UvA Quantum Hall effect Critical behavior and scaling Magnetotransport at the plateau-insulator transition Scaling functions
More informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
More informationA variational approach to the macroscopic. Electrodynamics of hard superconductors
A variational approach to the macroscopic electrodynamics of hard superconductors Dept. of Mathematics, Univ. of Rome La Sapienza Torino, July 4, 2006, joint meeting U.M.I. S.M.F. Joint work with Annalisa
More informationWeierstraß Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 0946 8633 Discrete Sobolev-Poincaré inequalities for Voronoi finite volume approximations Annegret
More informationComparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems
Comparative Analysis of Mesh Generators and MIC(0) Preconditioning of FEM Elasticity Systems Nikola Kosturski and Svetozar Margenov Institute for Parallel Processing, Bulgarian Academy of Sciences Abstract.
More informationE&M. 1 Capacitors. January 2009
E&M January 2009 1 Capacitors Consider a spherical capacitor which has the space between its plates filled with a dielectric of permittivity ɛ. The inner sphere has radius r 1 and the outer sphere has
More informationNumerical construction of Wannier functions
July 12, 2017 Internship Tutor: A. Levitt School Tutor: É. Cancès 1/27 Introduction Context: Describe electrical properties of crystals (insulator, conductor, semi-conductor). Applications in electronics
More informationSpin-orbit coupling fields in Fe/GaAs heterostructures
Spin-orbit coupling fields in Fe/GaAs heterostructures Outline motivation a simplified model of the Fe/GaAs heterostructure extracting spin-orbit coupling parameters spin-orbit coupling field conclusions
More informationMotivation. Confined acoustics phonons. Modification of phonon lifetimes Antisymmetric Bulk. Symmetric. 10 nm
Motivation Confined acoustics phonons Modification of phonon lifetimes 0 0 Symmetric Antisymmetric Bulk 0 nm A. Balandin et al, PRB 58(998) 544 Effect of native oxide on dispersion relation Heat transport
More informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More informationOn Surface Meshes Induced by Level Set Functions
On Surface Meshes Induced by Level Set Functions Maxim A. Olshanskii, Arnold Reusken, and Xianmin Xu Bericht Nr. 347 Oktober 01 Key words: surface finite elements, level set function, surface triangulation,
More informationA Variational Analysis of a Gauged Nonlinear Schrödinger Equation
A Variational Analysis of a Gauged Nonlinear Schrödinger Equation Alessio Pomponio, joint work with David Ruiz Dipartimento di Meccanica, Matematica e Management, Politecnico di Bari Variational and Topological
More informationKerr black hole and rotating wormhole
Kerr Fest (Christchurch, August 26-28, 2004) Kerr black hole and rotating wormhole Sung-Won Kim(Ewha Womans Univ.) August 27, 2004 INTRODUCTION STATIC WORMHOLE ROTATING WORMHOLE KERR METRIC SUMMARY AND
More informationLectures on magnetism at the Fudan University, Shanghai October 2005
Lectures on magnetism at the Fudan University, Shanghai 10. 26. October 2005 Klaus Baberschke Institut für Experimentalphysik Freie Universität Berlin Arnimallee 14 D-14195 D Berlin-Dahlem Germany 1 Introduction
More informationCharge cloud simulations, overview and very first results
W eierstraß -Institut für Angew andte Analysis und Stochastik K. Gärtner Charge cloud simulations, overview and very first results SRH generation 6M pairs, 7.5ns Berlin, July 28, 2008 charge cloud sim.
More informationThe anisotropic Strang-Fix conditions
The anisotropic Strang-Fix conditions Ronny Bergmann Mathematical Image Processing and Data Analysis Department of Mathematics University of Kaiserslautern March 19, 2015 Mecklenburg Workshop Approximation
More informationJacco Snoeijer PHYSICS OF FLUIDS
Jacco Snoeijer PHYSICS OF FLUIDS dynamics dynamics freezing dynamics freezing microscopics of capillarity Menu 1. surface tension: thermodynamics & microscopics 2. wetting (statics): thermodynamics & microscopics
More informationCharge cloud simulations, status report
W eierstraß-institut für Angew andte Analysis und Stochastik K. Gärtner Charge cloud simulations, status report SRH generation 6M pairs, 7.5ns Hamburg, Oct. 14, 2008 XDAC meeting Outline - Short recapitulation
More informationShear Flow of a Nematic Liquid Crystal near a Charged Surface
Physics of the Solid State, Vol. 45, No. 6, 00, pp. 9 96. Translated from Fizika Tverdogo Tela, Vol. 45, No. 6, 00, pp. 5 40. Original Russian Text Copyright 00 by Zakharov, Vakulenko. POLYMERS AND LIQUID
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationOn the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods
On the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods Marlis Hochbruck, Jonas Köhler CRC Preprint 2017/34 (revised), March 2018 KARLSRUHE INSTITUTE OF TECHNOLOGY KIT The Research
More informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
More informationModeling of electronic properties of interfaces in solar cells
Modeling of electronic properties of interfaces in solar cells Annegret Glitzky and Alexander Mielke Funding period 06/2010-05/2014 Requested funding: 1 position DFG Research Center MATHEON Mathematics
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationLarge, non-saturating magnetoresistance in WTe 2
Large, non-saturating magnetoresistance in WTe 2 M. N. Ali et al., Nature 514, 205-208 (2014) Team 5 Alina Kononov, Sangjun Lee, Mao Lin November 21, 2014 Outline Introduction Summary of the paper Comparison
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationarxiv:gr-qc/ v1 19 Feb 2004
On the construction of global models describing isolated rotating charged bodies; uniqueness of the exterior gravitational field Raül Vera Dublin City University, Ireland. arxiv:gr-qc/0402086v1 19 Feb
More informationMorphological evolution of single-crystal ultrathin solid films
Western Kentucky University From the SelectedWorks of Mikhail Khenner March 29, 2010 Morphological evolution of single-crystal ultrathin solid films Mikhail Khenner, Western Kentucky University Available
More informationLIQUID CRYSTAL ELECTRO-OSMOSIS
LIQUID CRYSTAL ELECTRO-OSMOSIS Chris Conklin and Jorge Viñals School of Physics and Astronomy University of Minnesota Acknowledgments: Chenhui Peng and Oleg Lavrentovich (LCI, Kent State). Carme Calderer
More informationTHE EDDY CURRENT PROBLEM WITH TEMPERATURE DEPENDENT PERMEABILITY
Electronic Journal of Differential Equations, Vol. 003003, No. 91, pp. 1 5. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu login: ftp THE EDDY CURRENT PROBLEM
More informationBlack-hole binary inspiral and merger in scalar-tensor theory of gravity
Black-hole binary inspiral and merger in scalar-tensor theory of gravity U. Sperhake DAMTP, University of Cambridge General Relativity Seminar, DAMTP, University of Cambridge 24 th January 2014 U. Sperhake
More informationNumerical computation of the genus of an irreducible curve within an algebraic set
Numerical computation of the genus of an irreducible curve within an algebraic set Dan Bates Chris Peterson Andrew J. Sommese Charles W. Wampler March 21, 2007 Abstract The common zero locus of a set of
More informationf(r) theory, Cosmology and Beyond
f(r) theory, Cosmology and Beyond Li-Fang Li The institute of theoretical physics, Chinese Academy of Sciences 2011-7-30 1 / 17 Outline 1 Brief review of f(r) theory 2 Higher dimensional f(r) theory: case
More informationSupporting information for: Adjusting the thermoelectric properties in. copper(i) oxide - graphite - polymer pastes
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Supporting information for: Adjusting the thermoelectric properties in copper(i)
More informationSupplementary Information: Dispersions of ellipsoidal particles in a nematic liquid crystal
Supplementary Information: Dispersions of ellipsoidal particles in a nematic liquid crystal Mykola Tasinkevych, 1, 2, Frédéric Mondiot, 3 Olivier Mondain-Monval, 3 and Jean-Christophe Loudet 3, 1 Max-Planck-Institut
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationLinear stability of MHD configurations
Linear stability of MHD configurations Rony Keppens Centre for mathematical Plasma Astrophysics KU Leuven Rony Keppens (KU Leuven) Linear MHD stability CHARM@ROB 2017 1 / 18 Ideal MHD configurations Interested
More informationIntrinsic finite element modeling of a linear membrane shell problem
arxiv:3.39v [math.na] 5 Mar Intrinsic finite element modeling of a linear membrane shell problem Peter Hansbo Mats G. Larson Abstract A Galerkin finite element method for the membrane elasticity problem
More informationSmall Molecule Crystallography Lab Department of Chemistry and Biochemistry University of Oklahoma 101 Stephenson Parkway Norman, OK
Small Molecule Crystallography Lab Department of Chemistry and Biochemistry University of Oklahoma 101 Stephenson Parkway Norman, OK 73019-5251 Sample: KP-XI-furan-enzymatic alcohol Lab ID: 12042 User:
More informationWeierstraß-Institut. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 0946 8633 Transient conductive-radiative heat transfer: Discrete existence and uniqueness for a
More informationAutomatic optimization of the cooling of injection mold based on the boundary element method
Automatic optimization of the cooling of injection mold based on the boundary element method E. Mathey, L. Penazzi, F.M. Schmidt and F. Rondé-Oustau Research Center on Tools, Materials and Forming Processes,
More informationStorm Open Library 3.0
S 50% off! 3 O L Storm Open Library 3.0 Amor Sans, Amor Serif, Andulka, Baskerville, John Sans, Metron, Ozdoby,, Regent, Sebastian, Serapion, Splendid Quartett, Vida & Walbaum. d 50% f summer j sale n
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationFinite Element Simulation of the Liquid Silicon Oriented Crystallization in a Graphite Mold
Finite Element Simulation of the Liquid Silicon Oriented Crystallization in a Graphite Mold Alexey I. Borovkov Vadim I. Suprun Computational Mechanics Laboratory, St. Petersburg State Technical University,
More informationRBF-FD Approximation to Solve Poisson Equation in 3D
RBF-FD Approximation to Solve Poisson Equation in 3D Jagadeeswaran.R March 14, 2014 1 / 28 Overview Problem Setup Generalized finite difference method. Uses numerical differentiations generated by Gaussian
More informationOn the parameter choice in grad-div stabilization for the Stokes equations
Noname manuscript No. (will be inserted by the editor) On the parameter choice in grad-div stabilization for the Stokes equations Eleanor W. Jenkins Volker John Alexander Linke Leo G. Rebholz Abstract
More informationA Continuum Theory for the Natural Vibrations of Spherical Viral Capsids
Viral Capsids A Continuum Theory for the Natural Vibrations of Spherical Viral Capsids Francesco Bonaldi Master s Thesis in Mathematical Engineering Advisors Prof. Paolo Podio-Guidugli Università di Roma
More information