Search results

This paper aims to establish a lattice Boltzmann (LB) method for solid-liquid phase transition (SLPT) from the pore scale to the representative elementary volume (REV) scale. By applying this method, detailed information about heat transfer and phase change processes within the pores can be obtained, while also enabling the calculation of larger-scale SLPT problems, such as shell-and-tube phase change heat storage systems.

Three-dimensional (3D) pore-scale enthalpy-based LB model is developed. The computational input parameters at the REV scale are derived from calculations at the pore scale, ensuring consistency between the two scales. The approaches to reconstruct the 3D porous structure and determine the REV of metal foam were discussed. The implementation of conjugate heat transfer between the solid matrix and the solid−liquid phase change material (SLPCM) for the proposed model is developed. A simple REV-scale LB model under the local thermal nonequilibrium condition is presented. The method of bridging the gap between the pore-scale and REV-scale enthalpy-based LB models by the REV is given.

This coupled method facilitates detailed simulations of flow, heat transfer and phase change within pores. The approach holds promise for multiscale calculations in latent heat storage devices with porous structures. The SLPT of the heat sinks for electronic device thermal control was simulated as a case, demonstrating the efficiency of the present models in designing and optimizing SLPT devices.

A coupled pore-scale and REV-scale LB method as a numerical tool for investigating phase change in porous materials was developed. This innovative approach allows for the capture of details within pores while addressing computations over a large domain. The LB method for simulating SLPT from the pore scale to the REV scale was given. The proposed method addresses the conjugate heat transfer between the SLPCM and the solid matrix in the enthalpy-based LB model.

The purpose of this study is to explore numerical scrutinization of micropolar and Walters-B non-Newtonian fluids motion under the influence of thermal radiation and chemical reaction.

The two fluids micropolar and Walters-B liquid are considered to start flowing from the slot to the stretching sheet. A magnetic field of constant strength is imposed on their flow transversely. The problems on heat and mass transport are set up with thermal, chemical reaction, heat generation, etc. to form partial differential equations. These equations were simplified into a dimensionless form and solved using spectral homotopy analysis method (SHAM). SHAM uses the basic concept of both Chebyshev pseudospectral method and homotopy analysis method to obtain numerical computations of the problem.

The outcomes for encountered flow parameters for temperature, velocity and concentration are presented with the aid of figures. It is observed that both the velocity and angular velocity of micropolar and Walters-B and thermal boundary layers increase with increase in the thermal radiation parameter. The decrease in velocity and decrease in angular velocity occurred are a result of increase in chemical reaction. It is hoped that the present study will enhance the understanding of boundary layer flow of micropolar and Walters-B non-Newtonian fluid under the influences of thermal radiation, thermal conductivity and chemical reaction as applied in various engineering processes.

All results are presented graphically and all physical quantities are computed and tabulated.

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.

This review aims to give a critical view of the wheel/rail high frequency vibration-induced vibration fatigue in railway bogie.

Vibration fatigue of railway bogie arising from the wheel/rail high frequency vibration has become the main concern of railway operators. Previous reviews usually focused on the formation mechanism of wheel/rail high frequency vibration. This paper thus gives a critical review of the vibration fatigue of railway bogie owing to the short-pitch irregularities-induced high frequency vibration, including a brief introduction of short-pitch irregularities, associated high frequency vibration in railway bogie, typical vibration fatigue failure cases of railway bogie and methodologies used for the assessment of vibration fatigue and research gaps.

The results showed that the resulting excitation frequencies of short-pitch irregularity vary substantially due to different track types and formation mechanisms. The axle box-mounted components are much more vulnerable to vibration fatigue compared with other components. The wheel polygonal wear and rail corrugation-induced high frequency vibration is the main driving force of fatigue failure, and the fatigue crack usually initiates from the defect of the weld seam. Vibration spectrum for attachments of railway bogie defined in the standard underestimates the vibration level arising from the short-pitch irregularities. The current investigations on vibration fatigue mainly focus on the methods to improve the accuracy of fatigue damage assessment, and a systematical design method for vibration fatigue remains a huge gap to improve the survival probability when the rail vehicle is subjected to vibration fatigue.

The research can facilitate the development of a new methodology to improve the fatigue life of railway vehicles when subjected to wheel/rail high frequency vibration.

This study establishes a method for predicting ground vibrations caused by railway tunnels in unsaturated soils with spatial variability.

First, an improved 2.5D finite-element-method-perfect-matching-layer (FEM-PML) model is proposed. The Galerkin method is used to derive the finite element expression in the ub-pl-pg format for unsaturated soil. Unlike the ub-v-w format, which has nine degrees of freedom per node, the ub-pl-pg format has only five degrees of freedom per node; this significantly enhances the calculation efficiency. The stretching function of the PML is adopted to handle the unlimited boundary domain. Additionally, the 2.5D FEM-PML model couples the tunnel, vehicle and track structures. Next, the spatial variability of the soil parameters is simulated by random fields using the Monte Carlo method. By incorporating random fields of soil parameters into the 2.5D FEM-PML model, the effect of soil spatial variability on ground vibrations is demonstrated using a case study.

The spatial variability of the soil parameters primarily affected the vibration acceleration amplitude but had a minor effect on its spatial distribution and attenuation over time. In addition, ground vibration acceleration was more affected by the spatial variability of the soil bulk modulus of compressibility than by that of saturation.

Using the 2.5D FEM-PML model in the ub-pl-pg format of unsaturated soil enhances the computational efficiency. On this basis, with the random fields established by Monte Carlo simulation, the model can calculate the reliability of soil dynamics, which was rarely considered by previous models.

The susceptible-infectious-susceptible (SIS) infectious disease models without spatial heterogeneity have limited applications, and the numerical simulation without considering models’ global existence and uniqueness of classical solutions might converge to an impractical solution. This paper aims to develop a robust and reliable numerical approach to the SIS epidemic model with spatial heterogeneity, which characterizes the horizontal and vertical transmission of the disease.

This study used stability analysis methods from nonlinear dynamics to evaluate the stability of SIS epidemic models. Additionally, the authors applied numerical solution methods from diffusion equations and heat conduction equations in fluid mechanics to infectious disease transmission models with spatial heterogeneity, which can guarantee a robustly stable and highly reliable numerical process. The findings revealed that this interdisciplinary approach not only provides a more comprehensive understanding of the propagation patterns of infectious diseases across various spatial environments but also offers new application directions in the fields of fluid mechanics and heat flow. The results of this study are highly significant for developing effective control strategies against infectious diseases while offering new ideas and methods for related fields of research.

Through theoretical analysis and numerical simulation, the distribution of infected persons in heterogeneous environments is closely related to the location parameters. The finding is suitable for clinical use.

The theoretical analysis of the stability theorem and the threshold dynamics guarantee robust stability and fast convergence of the numerical solution. It opens up a new window for a robust and reliable numerical study.

The paper aims to present an experimental and numerical investigation of the load–settlement behavior of soil reinforced by stone column, as well as to evaluate the plane strain unit cell model for the analysis of stone columns.

The numerical analysis was done using both axisymmetric and plane strain models. The elastic perfectly plastic behavior of Mohr–Coulomb was adopted for both soil and column material. The numerical results of this study were validated by the comparison with the in-situ measurements of a full-scale loading test on a stone column. This study also evaluated the effect of different parameters involved in the design of a stone column, including Young’s modulus of the column material, column diameter, spacing between the stone columns and Poisson’s ratio of the column material.

After the numerical simulation, the results from both axisymmetric and plane strain models are quite comparable. In addition, the numerical results revealed that the stone column with low spacing, a large diameter and a high Young’s modulus indicated better behavior against the settlement.

The axisymmetric unit cell model was used in many numerical studies on the behavior of stone columns. In the present work, a field load test on stone column was simulated using a plane strain unit cell model. This research adds that the plane strain unit cell model can be used to predict the settlement of reinforced soil with stone columns.

The purpose of this paper is to investigate the impact of near-wall treatment approaches, which are crucial parameters in predicting the flow characteristics of open channels, and the influence of different vegetation covers in different layers.

Ansys Fluent, a computational fluid dynamics software, was used to calculate the flow and turbulence characteristics using a three-dimensional, turbulent (k-e realizable), incompressible and steady-flow assumption, along with various near-wall treatment approaches (standard, scalable, non-equilibrium and enhanced) in the vegetated channel. The numerical study was validated concerning an experimental study conducted in the existing literature.

The numerical model successfully predicted experimental results with relative error rates below 10%. It was determined that nonequilibrium wall functions exhibited the highest predictive success in experiment Run 1, standard wall functions in experiment Run 2 and enhanced wall treatments in experiment Run 3. This study has found that plant growth significantly alters open channel flow. In the contact zones, the velocities and the eddy viscosity are low, while in the free zones they are high. On the other hand, the turbulence kinetic energy and turbulence eddy dissipation are maximum at the solid–liquid interface, while they are minimum at free zones.

This is the first study, to the best of the author’s knowledge, concerning the performance of different near-wall treatment approaches on the prediction of vegetation-covered open channel flow characteristics. And this study provides valuable insights to improve the hydraulic performance of open-channel systems.

This study aims to investigate the behavior of different types of stone columns, including the short and floating columns, as well as the ordinary and the geosynthetic encased stone columns (OSC and GESC). The effectiveness of the geosynthetic encasement and the impact of the installation using the lateral expansion method on the column performance is evaluated through a three-dimensional (3D) unit cell numerical analysis.

A full 3D numerical analysis is carried out using the explicit finite element code PLAXIS 3D to examine the installation influence on settlement reduction (ß), lateral displacement (U_x) and vertical displacement (U_z) relative to different values of lateral expansion of the column (0% to 15%).

The findings demonstrate the superior performance of GESC, particularly short columns outperforming floating counterparts. This enhanced performance is attributed to the combined effects of geosynthetic encasement and increased lateral expansion. Notably, these strategies contribute significantly to decreasing lateral displacement (U_x) at the column’s edge and reducing vertical displacement (U_z) under the rigid footing.

In contrast to previous studies that examined the installation effect of OSC contexts, this paper presents a comprehensive investigation into the effect of geosynthetic encasement and the installation effects using the lateral expansion method in very soft soil, using 3D numerical simulation. The study emphasizes the significance of the consideration of geosynthetic encasement and lateral expansion of the column during the design process to enhance column performance.

In recent times, there has been a growing interest in buoyancy-induced heat transfer within confined enclosures due to its frequent occurrence in heat transfer processes across diverse engineering disciplines, including electronic cooling, solar technologies, nuclear reactor systems, heat exchangers and energy storage systems. Moreover, the reduction of entropy generation holds significant importance in engineering applications, as it contributes to enhancing thermal system performance. This study, a numerical investigation, aims to analyze entropy generation and natural convection flow in an inclined square enclosure filled with Ag–MgO/water and Ag–TiO2/water hybrid nanofluids under the influence of a magnetic field. The enclosure features heated slits along its bottom and left walls. Following the Boussinesq approximation, the convective flow arises from a horizontal temperature difference between the partially heated walls and the cold right wall.

The governing equations for laminar unsteady natural convection flow in a Newtonian, incompressible mixture is solved using a Marker-and-Cell-based finite difference method within a customized MATLAB code. The hybrid nanofluid’s effective thermal conductivity and viscosity are determined using spherical nanoparticle correlations.

The numerical investigations cover various parameters, including nanoparticle volume concentration, Hartmann number, Rayleigh number, heat source/sink effects and inclination angle. As the Hartmann and Rayleigh numbers increase, there is a significant enhancement in entropy generation. The average Nusselt number experiences a substantial increase at extremely high values of the Rayleigh number and inclination.

This numerical investigation explores advanced applications involving various combinations of influential parameters, different nanoparticles, enclosure inclinations and improved designs. The goal is to control fluid flow and enhance heat transfer rates to meet the demands of the Fourth Industrial Revolution.

In a 90° tilted enclosure, the addition of 5% hybrid nanoparticles to the base fluid resulted in a 17.139% increase in the heat transfer rate for Ag–MgO nanoparticles and a 16.4185% increase for Ag–TiO2 nanoparticles compared to the base fluid. It is observed that a 5% nanoparticle volume fraction results in an increased heat transfer rate, influenced by variations in both the Darcy and Rayleigh numbers. The study demonstrates that the Ag–MgO hybrid nanofluid exhibits superior heat transfer and fluid transport performance compared to the Ag–TiO2 hybrid nanofluid. The simulations pertain to the use of hybrid magnetic nanofluids in fuel cells, solar cavity receivers and the processing of electromagnetic nanomaterials in enclosed environments.