𝔖 Bobbio Scriptorium
✦   LIBER   ✦

A non-physical enthalpy method for the numerical solution of isothermal solidification

✍ Scribed by K. Davey; R. Mondragon


Publisher
John Wiley and Sons
Year
2010
Tongue
English
Weight
473 KB
Volume
84
Category
Article
ISSN
0029-5981

No coin nor oath required. For personal study only.

✦ Synopsis


Abstract

In this paper, a new formulation for the solution of the discontinuous isothermal solidification problem is presented. The formulation has similarities with the now classical capacitance and source methods traditionally used in commercial software. However, the new approach focuses on the solution of the governing enthalpy‐transport equation rather than the governing parabolic partial differential heat equation. The advantage is that discontinuous physics can be accounted for without approximation and the arbitrariness common to classic approaches is avoided. Also introduced is the concept of non‐physical enthalpy, which unlike physical enthalpy has numerical values that are not moving‐frame invariant. Understanding the behaviour of the non‐physical enthalpy is central to the successful treatment of discontinuities. A particular drawback is that non‐physical enthalpy is non‐intuitive and new mathematical constructs are required to describe its behaviour. This involves the introduction of transport equations, which provide the new concept of relative moving‐frame invariance for the non‐physical enthalpy. The principal advantage, however, is that a unified methodology is established for the treatment of discontinuities. This is shown to establish real rigour and in many respects the formulation highlights the erroneous choices made with established classical approaches and casts in a totally new light a somewhat traditional problem. The new methodology is applied to a range of simple problems not only to provide an in‐depth treatment and for ease of understanding but also to best describe the behaviour of the non‐intuitive non‐physical enthalpy. Demonstrated in the paper is the methods' remarkable accuracy and stability. Copyright © 2010 John Wiley & Sons, Ltd.


📜 SIMILAR VOLUMES


Existence of solutions to a phase-field
✍ J. Rappaz; J. F. Scheid 📂 Article 📅 2000 🏛 John Wiley and Sons 🌐 English ⚖ 180 KB 👁 2 views

We investigate the well-posedness of a phase-"eld model for the isothermal solidi"cation of a binary alloy due to Warren}Boettinger [12]. Existence of weak solution as well as regularity and uniqueness results are established under Lipschitz and boundedness assumptions for the non-linearities. A max

A Numerical Method for Solution of the G
✍ J. Candy 📂 Article 📅 1996 🏛 Elsevier Science 🌐 English ⚖ 354 KB

dq N dp N ϭ const, (3a) A numerical method for the time evolution of systems described by Liouville-type equations is derived. The algorithm uses a lattice of numerical markers, which follow exactly Hamiltonian trajectories, to represent the operator d/dt in moving (i.e., Lagrangian) coordinates. H

A Wavelet Collocation Method for the Num
✍ S. Bertoluzza; G. Naldi 📂 Article 📅 1996 🏛 Elsevier Science 🌐 English ⚖ 355 KB

We describe a wavelet collocation method for the numerical solution of partial differential equations which is based on the use of the autocorrelation functions of Daubechie's compactly supported wavelets. For such a method we discuss the application of wavelet based preconditioning techniques along

A numerical method for the solution of t
✍ Mailand R. Strunk; Frank F. Tao 📂 Article 📅 1964 🏛 American Institute of Chemical Engineers 🌐 English ⚖ 509 KB

T = pressure, atm. u = collision diameter, A (a = azimuthal angle between the axes of the two di-(p(r) = Stockmayer potential, Equation (6) f i ( l J ) \* [ T ~] = reduced collision integral for the Lennardfi(2,2)\* [ T N ] = reduced collision integral for the Lennard-O ( l J ) \* [ TN, S o ] = redu

A Conservative Finite Difference Method
✍ Phillip Colella; Milo R Dorr; Daniel D Wake 📂 Article 📅 1999 🏛 Elsevier Science 🌐 English ⚖ 328 KB

This paper describes a numerical method for the solution of a system of plasma fluid equations. The fluid model is similar to those employed in the simulation of high-density, low-pressure plasmas used in semiconductor processing. The governing equations consist of a drift-diffusion model of the ele