Search results

The purpose of this paper is to overview successful approaches to the computational simulation of real fluid (R744 – carbon dioxide (CO₂)) flow within an ejector is presented. Important issues such as the ejector geometry and its optimisation, the adapted equations of state and the proposed models of the process, fluid parameters, etc., are examined and critically discussed. Whenever possible, the discussed models are experimentally validated. In the conclusion, some trends in future research are pointed out.

Flow within CO₂ ejector is generally transcritical and compressible. Models existing in the literature are shortly described and critically compared. Whenever possible, those models were validated against the experimental data. In a model validation process, the primary and secondary mass flow rates as well as the pressures at the selected points in the mixing section and diffuser were compared, showing a satisfactory agreement between experimental and computational results.

Developed CO₂ ejector flow models are tested in few industrial applications. All these initiatives bring solutions which are interesting and very promising from technological point of view.

This is an extensive overview of successful approaches to computational simulation of the real fluid (R744 – CO₂) flow within ejector. It brings many useful information.

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.

A phase‐change temperature‐based formulation including general latent heat effects is presented. These effects are taken into account by means of an explicit “phase‐change function” (or liquid fraction‐temperature relationship in a more specific context) defined analytically or based on experimental measurements. The behaviour of different functions is studied and compared. The finite element equations of this formulation are also described. Finally, a numerical example is analysed to illustrate the performance of the proposed methodology.

The homographic approximation, in which the Heaviside step function is replaced by a continuous smooth curve, is applied to the enthalpy method for heat transfer problems with isothermal phase change. Both the finite difference and finite element implementations, based on the basic enthalpy, the apparent heat capacity and the source term formulations, are considered. A 1‐D Stefan problem of melting a solid is used as a test problem. The accuracy of the numerical solutions is measured globally using L2 error norms and comparison is made between the solutions using homographic approximation and those using linear approximation. The advantages of using homographic approximation are examined.

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 present original methods related to the modeling of material deposit and associated heat sources for finite element simulation of gas metal arc welding (GMAW).

The filler deposition results from high-frequency impingements of melted droplets. The present modeling approach consists of a time-averaged source term in the mass equation for selected finite elements in the fusion zone. The associated expansion of the mesh is controlled by means of adaptive remeshing. The heat input includes a volume source corresponding to the droplets energy, for which a model from the literature is expressed in coherency with mass supply. Finally, an inverse technique has been developed to identify different model parameters. The objective function includes the differences between calculations and experiments in terms of temperature, but also shape of the fusion zone.

The proposed approach for the modeling of metal deposition results in a direct calculation of the formation of the weld bead, without any a priori definition of its shape. Application is shown on GMAW of steel 316LN, for which parameters of the model have been identified by the inverse method. They are in agreement with literature and simulation results are found quite close to experimental measurements.

The proposed algorithm for material deposit offers an alternative to the element activation techniques that are commonly used to simulate the deposition of filler metal. The proposed inverse method for parameter identification is original in that it encompasses an efficient and convenient technique to take into account the shape of the fusion zone.

The freezing phenomenon of a falling droplet is a frequently encountered phenomenon in various applications, such as spray crystallization, hail formation and artificial snowmaking. Therefore, this paper aims to understand the freezing processes of a falling droplet without and with initial horizontal velocity in a cold space.

The freezing processes of a falling droplet were characterized using a modified enthalpy-based lattice Boltzmann method.

The temperature field, streamlines and freezing process of the falling droplet were investigated and analyzed. The lower part of the droplet was frozen earlier than the upper part. The freezing trend slowed down in the later stage of the freezing process. The droplet shape was related to the initial vertical velocity, nucleation temperature and initial horizontal velocity.

A modified enthalpy-based lattice Boltzmann method is proposed. In the model, the improved pseudo-potential model is used and the radiation is considered. This method was firstly used to simulate the freezing process of a falling droplet. By examining these freezing processes in detail, the freezing trend and the effect factors of droplet deformation and freezing time were obtained, respectively.

This paper aims to propose a semi-analytical benchmarking framework for enthalpy-based methods used in problems involving phase change with latent heat. The benchmark is based on a class of semi-analytical solutions of spatially symmetric Stefan problems in an arbitrary spatial dimension. Via a public repository this study provides a finite element numerical code based on the FEniCS computational platform, which can be used to test and compare any method of choice with the (semi-)analytical solutions. As a particular demonstration, this paper uses the benchmark to test several standard temperature-based implementations of the enthalpy method and assesses their accuracy and stability with respect to the discretization parameters.

The class of spatially symmetric semi-analytical self-similar solutions to the Stefan problem is found for an arbitrary spatial dimension, connecting some of the known results in a unified manner, while providing the solutions’ existence and uniqueness. For two chosen standard semi-implicit temperature-based enthalpy methods, the numerical error assessment of the implementations is carried out in the finite element formulation of the problem. This paper compares the numerical approximations to the semi-analytical solutions and analyzes the influence of discretization parameters, as well as their interdependence. This study also compares accuracy of these methods with other traditional approach based on time-explicit treatment of the effective heat capacity with and without iterative correction.

This study shows that the quantitative comparison between the semi-analytical and numerical solutions of the symmetric Stefan problems can serve as a robust tool for identifying the optimal values of discretization parameters, both in terms of accuracy and stability. Moreover, this study concludes that, from the performance point of view, both of the semi-implicit implementations studied are equivalent, for optimal choice of discretization parameters, they outperform the effective heat capacity method with iterative correction in terms of accuracy, but, by contrast, they lose stability for subcritical thickness of the mushy region.

The proposed benchmark provides a versatile, accessible test bed for computational methods approximating multidimensional phase change problems. The supplemented numerical code can be directly used to test any method of choice against the semi-analytical solutions.

While the solutions of the symmetric Stefan problems for individual spatial dimensions can be found scattered across the literature, the unifying perspective on their derivation presented here has, to the best of the authors’ knowledge, been missing. The unified formulation in a general dimension can be used for the systematic construction of well-posed, reliable and genuinely multidimensional benchmark experiments.

The purpose of this paper is to develop a coupled lattice Boltzmann model for the simulation of the freezing process in unsaturated porous media.

In the developed model, the porous structure with complexity and disorder was generated by using a stochastic growth method, and then the Shan-Chen multiphase model and enthalpy-based phase change model were coupled by introducing a freezing interface force to describe the variation of phase interface. The pore size of porous media in freezing process was considered as an influential factor to phase transition temperature, and the variation of the interfacial force formed with phase change on the interface was described.

The larger porosity (0.2 and 0.8) will enlarge the unfrozen area from 42 mm to 70 mm, and the rest space of porous medium was occupied by the solid particles. The larger specific surface area (0.168 and 0.315) has a more fluctuated volume fraction distribution.

The concept of interfacial force was first introduced in the solid–liquid phase transition to describe the freezing process of frozen soil, enabling the formulation of a distribution equation based on enthalpy to depict the changes in the water film. The increased interfacial force serves to diminish ice formation and effectively absorb air during the freezing process. A greater surface area enhances the ability to counteract liquid migration.

In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues.

The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed.

Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated.

As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.