scholarly journals Precise Numerical Differentiation of Thermodynamic Functions with Multicomplex Variables

Author(s):  
Ulrich K. Deiters ◽  
Ian H. Bell

The multicomplex finite-step method for numerical differentiation is an extension of the popular Squire–Trapp method, which uses complex arithmetics to compute first-order derivatives with almost machine precision. In contrast to this, the multicomplex method can be applied to higher-order derivatives. Furthermore, it can be applied to functions of more than one variable and obtain mixed derivatives. It is possible to compute various derivatives at the same time. This work demonstrates the numerical differentiation with multicomplex variables for some thermodynamic problems. The method can be easily implemented into existing computer programs, applied to equations of state of arbitrary complexity, and achieves almost machine precision for the derivatives. Alternative methods based on complex integration are discussed, too.

2021 ◽  
Vol 5 (4) ◽  
pp. 203
Author(s):  
Suzan Cival Buranay ◽  
Nouman Arshad ◽  
Ahmed Hersi Matan

We give fourth-order accurate implicit methods for the computation of the first-order spatial derivatives and second-order mixed derivatives involving the time derivative of the solution of first type boundary value problem of two dimensional heat equation. The methods are constructed based on two stages: At the first stage of the methods, the solution and its derivative with respect to time variable are approximated by using the implicit scheme in Buranay and Arshad in 2020. Therefore, Oh4+τ of convergence on constructed hexagonal grids is obtained that the step sizes in the space variables x1, x2 and in time variable are indicated by h, 32h and τ, respectively. Special difference boundary value problems on hexagonal grids are constructed at the second stages to approximate the first order spatial derivatives and the second order mixed derivatives of the solution. Further, Oh4+τ order of uniform convergence of these schemes are shown for r=ωτh2≥116,ω>0. Additionally, the methods are applied on two sample problems.


Author(s):  
Vasily I. Repnikov ◽  
Boris V. Faleichik ◽  
Andrew V. Moisa

In this work we present explicit Adams-type multi-step methods with extended stability intervals, which are analogous to the stabilised Chebyshev Runge – Kutta methods. It is proved that for any k ≥ 1 there exists an explicit k-step Adams-type method of order one with stability interval of length 2k. The first order methods have remarkably simple expressions for their coefficients and error constant. A damped modification of these methods is derived. In the general case, to construct a k-step method of order p it is necessary to solve a constrained optimisation problem in which the objective function and p constraints are second degree polynomials in k variables. We calculate higher-order methods up to order six numerically and perform some numerical experiments to confirm the accuracy and stability of the methods.


Author(s):  
Greg M. Anderson ◽  
David A. Crerar

The procedures described in Chapter 15 are well suited to solid and liquid solutions and could also be applied to gases, but in fact, other approaches are generally used. The main reason for this is partly historical; much work was done early in the history of physical chemistry on the behavior of gases, and these methods have continued to evolve to the present day. We have also just seen that the Margules equations become very unwieldy with multi-component systems. Because true gases are completely miscible, natural gases often contain many different components, so the Margules approach is not very suitable. Unfortunately, the most successful alternative methods described in this section are also quite unwieldy; however, they do not become much more complicated for multi-component gases than they are for the pure gases themselves, and this is a definite advantage. We have seen that with real, non-ideal gases, all the thermodynamic properties are described if we know the T, P, and the fugacity coefficient. For gaseous solutions, the fugacity coefficient for each component generally depends on the concentrations and types of other gaseous species in the same mixture. All gases, whether pure or multi-component, should approach ideality at higher T and lower P; conversely, non-ideality is most pronounced in dense, low-temperature gases where intermolecular forces are strongest. The challenge here is to find an equation of state that can adequately cover this range of conditions for gases of many different constituents. In the following discussion we first briefly outline some of the equations of state used to describe pure gases. We will introduce these from the molecular point of view since this helps understand the physical basis (and limitations) of each model. Each of these equations of state can then be applied to mixtures of gases using a set of rules which we describe at the end of this section.


1992 ◽  
Vol 03 (05) ◽  
pp. 1071-1082 ◽  
Author(s):  
D.B. ABRAHAM ◽  
P.J. UPTON

Problems associated with analyticity of thermodynamic functions close to first-order phase transitions are briefly reviewed. The bubble model for correlation functions is then applied to planar Ising-like models at subcritical temperatures (T<Tc) with a bulk magnetic field h. The fluctuation sum is used to calculate the susceptibility χ(h) from the bubble correlation function. We show that χ(h), calculated this way, must contain an essential singularity at h=0 i.e. at the first-order phase boundary. This has important implications to metastability, where we demonstrate that if the ensemble is restricted such that the magnetization stays positive when h goes negative, χ(h) has an infinite number of poles along the negative real axis with a limit-point at h=0. For an unrestricted ensemble, a Yang-Lee circle theorem is derived.


2006 ◽  
Vol 03 (01) ◽  
pp. 27-52 ◽  
Author(s):  
FRANÇOIS BOUCHUT ◽  
HERMANO FRID

We propose finite difference schemes for multidimensional quasilinear parabolic systems whose main feature is the introduction of correctors which control the second-order terms with mixed derivatives. We show that with these correctors the schemes inherit physically relevant properties present at the continuous level, such as the existence of invariant domains and/or the nonincrease of the total amount of entropy. The analysis is performed with some general tools that could be used also in the analysis of finite volume methods based on flux vector splitting for first-order hyperbolic problems on unstructured meshes. Applications to the compressible Navier–Stokes system are given.


Sign in / Sign up

Export Citation Format

Share Document