Search results

Thermal solidification processes are an important concern in today’s manufacturing technology. Because of the complex geometric nature of real‐world problems, analytical techniques with closed‐form solutions are scarce and/or not feasible. As a consequence, various numerical techniques have been employed for the numerical simulations. Of interest in the present paper are thermal solidification problems involving single or multiple arbitary phases. In order to effectively handle such problems, the finite element method is employed in conjunction with adaptive time stepping approaches to accurately and effectively track the various phase fronts and describe the physics of phase front interactions and thermal behaviour. In conjunction with the enthalpy method which is employed to handle the latent heat release, a fixed‐grid finite element technique and an automatic time stepping approach which uses the norm of the temperature distribution differences between adjacent time step levels to control the error are employed with the scale of the norm being automatically selected. Several numerical examples, including single and multiple phase change problems, are described.

The purpose of this paper is to describe a novel implementation of a multispatial method, multitime-scheme subdomain differential algebraic equation (DAE) framework allowing a mix of different space discretization methods and different time schemes by a robust generalized single step single solve (GS4) family of linear multistep (LMS) algorithms on a single body analysis for the first-order nonlinear transient systems.

This proposed method allows the coupling of different numerical methods, such as the finite element method and particle methods, and different implicit and/or explicit algorithms in each subdomain into a single analysis with the GS4 framework. The DAE, which constrains both space and time in multi-subdomain analysis, combined with the GS4 framework ensures the second-order time accuracy in all primary variables and Lagrange multiplier. With the appropriate GS4 parameters, the algorithmic temperature rate variable shift can be matched for all time steps using the DAE. The proposed method is used to solve various combinations of spatial methods and time schemes between subdomains in a single analysis of nonlinear first-order system problems.

The proposed method is capable of coupling different spatial methods for multiple subdomains and different implicit/explicit time integration schemes in the GS4 framework while sustaining second-order time accuracy.

Traditional approaches do not permit such robust and flexible coupling features. The proposed framework encompasses most of the LMS methods that are second-order time accurate and unconditionally stable.

The time‐discretization process of transient equation systems is an important concern in computational heat transfer applications. As such, the present paper describes a formal basis towards providing the theoretical concepts, evolution and development, and characterization of a wide class of time discretized operators for transient heat transfer computations. Therein, emanating from a common family tree and explained via a generalized time weighted philosophy, the paper addresses the development and evolution of time integral operators [IO], and leading to integration operators [InO] in time encompassing single‐step integration operators [SSInO], multi‐step integration operators [MSInO], and a class of finite element in time integration operators [FETInO] including the relationships and the resulting consequences. Also depicted are those termed as discrete numerically assigned [DNA] algorithmic markers essentially comprising of both: the weighted time fields, and the corresponding conditions imposed upon the dependent variable approximation, to uniquely characterize a wide class of transient algorithms. Thereby, providing a plausible standardized formal ideology when referring to and/or relating time discretized operators applicable to transient heat transfer computations.

The present paper describes the applicability of hybrid transfinite element modelling/analysis formulations for non‐linear heat conduction problems involving phase change. The methodology is based on application of transform approaches and classical Galerkin schemes with finite element formulations to maintain the modelling versatility and numerical features for computational analysis. In addition, in conjunction with the above, the effects due to latent heat are modelled using enthalpy formulations to enable a physically realistic approximation to be effectively dealt computationally for materials exhibiting phase change within a narrow band of temperatures. Pertinent details of the approach and computational scheme adapted are described in technical detail. Numerical test cases of comparative nature are presented to demonstrate the applicability of the proposed formulations for numerical modelling/analysis of non‐linear heat conduction problems involving phase change.

The purpose of this study is to further advance the multiple space/time subdomain framework with model reduction. Existing linear multistep (LMS) methods that are second-order time accurate, and useful for practical applications, have a significant limitation. They do not account for separable controllable numerical dissipation of the primary variables. Furthermore, they have little or no significant choices of altogether different algorithms that can be integrated in a single analysis to mitigate numerical oscillations that may occur. In lieu of such limitations, under the generalized single-step single-solve (GS4) umbrella, several of the deficiencies are circumvented.

The GS4 framework encompasses a wide variety of LMS schemes that are all second-order time accurate and offers controllable numerical dissipation. Unlike existing state-of-art, the present framework permits implicit–implicit and implicit–explicit coupling of algorithms via differential algebraic equations (DAE). As further advancement, this study embeds proper orthogonal decomposition (POD) to further reduce model sizes. This study also uses an iterative convergence check in acquiring sufficient snapshot data to adequately capture the physics to prescribed accuracy requirements. Simple linear/nonlinear transient numerical examples are presented to provide proof of concept.

The present DAE-GS4-POD framework has the flexibility of using different spatial methods and different time integration algorithms in altogether different subdomains in conjunction with the POD to advance and improve the computational efficiency.

The novelty of this paper is the addition of reduced order modeling features, how it applies to the previous DAE-GS4 framework and the improvement of the computational efficiency. The proposed framework/tool kit provides all the needed flexibility, robustness and adaptability for engineering computations.

The purpose of this paper is to describe a novel universal error estimator and the adaptive time-stepping process in the generalized single-step single-solve (GS4-1) computational framework, applied for the fluid dynamics with illustrations to incompressible Navier–Stokes equations.

The proposed error estimator is universal and versatile that it works for the entire subsets of the GS4-1 framework, encompassing the nondissipative Crank–Nicolson method, the most dissipative backward differential formula and anything in between. It is new and novel that the cumbersome design work of error estimation for specific time integration algorithms can be avoided. Regarding the numerical implementation, the local error estimation has a compact representation that it is determined by the time derivative variables at four successive time levels and only involves vector operations, which is simple for numerical implementation. Additionally, the adaptive time-stepping is further illustrated by the proposed error estimator and is used to solve the benchmark problems of lid-driven cavity and flow past a cylinder.

The proposed computational procedure is capable of eliminating the nonphysical oscillations in GS4-1(1,1)/Crank–Nicolson method; being CPU-efficient in both dissipative and nondissipative schemes with better solution accuracy; and detecting the complex physics and hence selecting a suitable time step according to the user-defined error threshold.

To the best of the authors’ knowledge, for the first time, this study applies the general purpose GS4-1 family of time integration algorithms for transient simulations of incompressible Navier–Stokes equations in fluid dynamics with constant and adaptive time steps via a novel and universal error estimator. The proposed computational framework is simple for numerical implementation and the time step selection based on the proposed error estimation is efficient, benefiting to the computational expense for transient simulations.

The purpose of this paper is to provide an overview and some recent advances in the models, analysis and simulation of thermal transport of phonons as related to the field of microscale/macroscale heat conduction in solids. The efforts focus upon a fairly comprehensive overview of the subject matter from a unified standpoint highlighting the various approximations inherent in the thermal models. Subsequently, the numerical formulations and illustrations using the current state-of-the-art are provided.

This paper is dedicated to the approximate solution to the relaxation time phonon Boltzmann equation (BE). While original contributions are pointed out and addressed appropriately, the efforts and contributions will be focussed on a relatively complete overview highlighting the field from one unified standpoint and clearly stating all assumptions that go into the approximations inherent to existing models. The contents will be divided as follows: In the first section the authors will give an overview of semi-classical phonon transport physics. Then the authors will discuss the equation of phonon radiative transport (EPRT) and its approximations—the ballistic-diffusive approximation (BDA) and the new heat equation (NHE). Next the authors derive and discuss the C-F model. A numerical discretization method valid for all models is then presented followed by results to numerical simulations and discussion.

From a unified treatment based on the introduction of an energy distribution function, the authors have derived the EPRT and its two well-known approximations: BDA and NHE. For completeness and to provide a vehicle for a general numerical discretization approach, the authors have also included analysis of the C-F model and the parabolic and hyperbolic descriptions of heat transfer along with it. The approximation of angular dependence of phonons in radiation-like descriptions of transport has been given special attention. The assumption of isotropy was found to be of paramount importance in the formulation of position space models for phononic thermal transport. For the thin film problem considered here, the NHE along with the proper boundary condition appears to be the best choice to approximate the phonon BE. Not only does it provide predictions that are in excellent agreement with EPRT, it does not require the discretization of phase space making it far more computationally efficient.

The authors hope this work will help dispel the idea that since Fourier’s law describes diffusion (under limiting assumptions) and it has shown to be ineffective in describing heat transfer for very thin films, that diffusion cannot describe heat transfer in thin films and one should look to a radiative description instead. If one considers diffusion in the sense of random motion, as invisaged by the original builders of the subject (Smoluchowski, Einstein, Ornstein et al.), instead of a temperature gradient, the idea that diffusion can govern thermal transport at this scale is not surprising. Indeed, the NHE is essentially a diffusion equation that describes the motion of particles up to the point of true randomness (isotropy) as well as thereafter.

We first provide an overview of some predominant theoretical methods currently used for predicting thermal conductivity of thin dielectric films: the equation of radiative transfer, the temperature‐dependent thermal conductivity theories based on the Callaway model, and the molecular dynamics simulation. This overview also highlights temporal and spatial scale issues by looking at a unified theory that bridges physical issues presented in the Fourier and Cattaneo models. This newly developed unified theory is the so‐called C‐ and F‐processes constitutive model. This model introduces the notion of a new dimensionless heat conduction model number, which is the ratio of the thermal conductivity of the fast heat carrier F‐processes to the total thermal conductivity comprised of both the fast heat carriers F‐processes and the slow heat carriers C‐processes. Illustrative numerical examples for prediction of thermal conductivity in thin films are primarily presented.

This paper aims to propose a gradient-based shape optimization framework in which traditional time-consuming conversions between computer-aided design and computer-aided engineering and the mesh update procedure are avoided/eliminated. The scheme is general so that it can be used in all cases as a black box, no matter what the objective and/or design variables are, whilst the efficiency and accuracy are guaranteed.

The authors integrated CAD and CAE by using isogeometric analysis (IGA), enabling the present methodology to be robust and accurate. To overcome the difficulty in evaluating the sensitivities of objective and/or constraint functions by analytic method in some cases, the authors adopt the finite difference method to calculate these sensitivities, thereby providing a universal approach. Moreover, to further eliminate the inefficiency caused by the finite difference method, the authors advance the exact reanalysis method, the indirect factorization updating (IFU), to exactly and efficiently calculate functions and their sensitivities, which guarantees its generality and efficiency at the same time.

The proposed isogeometric gradient-based shape optimization using our IFU approach is reliable and accurate, as well as general and efficient.

The authors proposed a gradient-based shape optimization framework in which they first integrate IGA and the proposed exact reanalysis method for applicability to structural response and sensitivity analysis.

The purpose of this paper is to design a simple and accurate a-posteriori Lagrangian-based error estimator is developed for the class of backward differentiation formula (BDF) algorithms with variable time step size, and the adaptive time-stepping in BDF algorithms is demonstrated for efficient time-dependent simulations in fluid flow and heat transfer.

The Lagrange interpolation polynomial is used to predict the time derivative, and then the accurate primary result is obtained by the Gauss integral, which is applied to evaluate the local error. Not only the generalized formula of the proposed error estimator is presented but also the specific expression for the widely applied BDF1/2/3 is illustrated. Two essential executable MATLAB functions to implement the proposed error estimator are appended for practical applications. Then, the adaptive time-stepping is demonstrated based on the newly proposed error estimator for BDF algorithms.

The validation tests show that the newly proposed error estimator is accurate such that the effectivity index is always close to unity for both linear and nonlinear problems, and it avoids under/overestimation of the exact local error. The applications for fluid dynamics and coupled fluid flow and heat transfer problems depict the advantage of adaptive time-stepping based on the proposed error estimator for time-dependent simulations.

In contrast to existing error estimators for BDF algorithms, the present work is more accurate for the local error estimation, and it can be readily extended to practical applications in engineering with a few changes to existing codes, contributing to efficient time-dependent simulations in fluid flow and heat transfer.