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 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
- DOI
- 10.1002/nme.2896
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
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
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
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
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