Search results

This study aims to analyse the non-linear losses of a porous media (stack) composed by parallel plates and inserted in a resonator tube in oscillatory flows by proposing numerical correlations between pressure gradient and velocity.

The numerical correlations origin from computational fluid dynamics simulations, conducted at the microscopic scale, in which three fluid channels representing the porous media are taken into account. More specifically, for a specific frequency and stack porosity, the oscillating pressure input is varied, and the velocity and the pressure-drop are post-processed in the frequency domain (Fast Fourier Transform analysis).

It emerges that the viscous component of pressure drop follows a quadratic trend with respect to velocity inside the stack, while the inertial component is linear also at high-velocity regimes. Furthermore, the non-linear coefficient b of the correlation ax + bx² (related to the Forchheimer coefficient) is discovered to be dependent on frequency. The largest value of the b is found at low frequencies as the fluid particle displacement is comparable to the stack length. Furthermore, the lower the porosity the higher the Forchheimer term because the velocity gradients at the stack geometrical discontinuities are more pronounced.

The main novelty of this work is that, for the first time, non-linear losses of a parallel plate stack are investigated from a macroscopic point of view and summarised into a non-linear correlation, similar to the steady-state and well-known Darcy–Forchheimer law. The main difference is that it considers the frequency dependence of both Darcy and Forchheimer terms. The results can be used to enhance the analysis and design of thermoacoustic devices, which use the kind of stacks studied in the present work.

The purpose of this study is to propose a precise and standardized strategy for numerically simulating vehicle aerodynamics.

Error sources in computational fluid dynamics were analyzed. Additionally, controllable experiential and discretization errors, which significantly influence the calculated results, are expounded upon. Considering the airflow mechanism around a vehicle, the computational efficiency and accuracy of each solution strategy were compared and analyzed through numerous computational cases. Finally, the most suitable numerical strategy, including the turbulence model, simplified vehicle model, calculation domain, boundary conditions, grids and discretization scheme, was identified. Two simplified vehicle models were introduced, and relevant wind tunnel tests were performed to validate the selected strategy.

Errors in vehicle computational aerodynamics mainly stem from the unreasonable simplification of the vehicle model, calculation domain, definite solution conditions, grid strategy and discretization schemes. Using the proposed standardized numerical strategy, the simulated steady and transient aerodynamic characteristics agreed well with the experimental results.

Building upon the modified Low-Reynolds Number k-e model and Scale Adaptive Simulation model, to the best of the authors’ knowledge, a precise and standardized numerical simulation strategy for vehicle aerodynamics is proposed for the first time, which can be integrated into vehicle research and design.

This study aims to develop an experimentally validated three-dimensional numerical model for predicting different flow patterns produced with a gas dynamic virtual nozzle (GDVN).

The physical model is posed in the mixture formulation and copes with the unsteady, incompressible, isothermal, Newtonian, low turbulent two-phase flow. The computational fluid dynamics numerical solution is based on the half-space finite volume discretisation. The geo-reconstruct volume-of-fluid scheme tracks the interphase boundary between the gas and the liquid. To ensure numerical stability in the transition regime and adequately account for turbulent behaviour, the k-ω shear stress transport turbulence model is used. The model is validated by comparison with the experimental measurements on a vertical, downward-positioned GDVN configuration. Three different combinations of air and water volumetric flow rates have been solved numerically in the range of Reynolds numbers for airflow 1,009–2,596 and water 61–133, respectively, at Weber numbers 1.2–6.2.

The half-space symmetry allows the numerical reconstruction of the dripping, jetting and indication of the whipping mode. The kinetic energy transfer from the gas to the liquid is analysed, and locations with locally increased gas kinetic energy are observed. The calculated jet shapes reasonably well match the experimentally obtained high-speed camera videos.

The model is used for the virtual studies of new GDVN nozzle designs and optimisation of their operation.

To the best of the authors’ knowledge, the developed model numerically reconstructs all three GDVN flow regimes for the first time.

1) The OE-MLPG penalty meshfree method is developed to solve cracked structure.(2) Smartening the numerical meshfree method by combining the particle swarm optimization (PSO) optimization algorithms and Voronoi computational geometric algorithm. (3). Selection of base functions, finding optimal penalty factor and distribution of appropriate nodal points to the accuracy of calculation in the meshless local Petrov–Galekrin (MLPG) meshless method.

Using appropriate shape functions and distribution of nodal points in local domains and sub-domains and choosing an approximation or interpolation method has an effective role in the application of meshless methods for the analysis of computational fracture mechanics problems, especially problems with geometric discontinuity and cracks. In this research, computational geometry technique, based on the Voronoi diagram (VD) and Delaunay triangulation and PSO algorithm, are used to distribute nodal points in the sub-domain of analysis (crack line and around it on the crack plane).

By doing this process, the problems caused by too closeness of nodal points in computationally sensitive areas that exist in general methods of nodal point distribution are also solved. Comparing the effect of the number of sentences of basic functions and their order in the definition of shape functions, performing the mono-objective PSO algorithm to find the penalty factor, the coefficient, convergence, arrangement of nodal points during the three stages of VD implementation and the accuracy of the answers found indicates, the efficiency of V-E-MLPG method with Ns = 7 and ß = 0.0037–0.0075 to estimation of 3D-stress intensity factors (3D-SIFs) in computational fracture mechanics.

The present manuscript is a continuation of the studies (Ref. [33]) carried out by the authors, about; feasibility assessment, improvement and solution of challenges, introduction of more capacities and capabilities of the numerical MLPG method have been used. In order to validate the modeling and accuracy of calculations, the results have been compared with the findings of reference article [34] and [35].

Water jet propulsion is widely used in various military and civilian fields due to its advantages of simple structure and high propulsion efficiency. The process of mooring involves utilizing specially designed equipment to secure a ship at a designated berth. During the process of water jet propulsion, the single propeller operates within a complex and turbulent three-dimensional flow. Hence, studying the coupling between the water jet propeller and the hull is critical to comprehending the characteristics of the device and the distribution of the flow field in detail.

Firstly, we conducted computational fluid dynamics (CFD)-based self-propulsion calculations to evaluate the interaction between the hull and the propeller. We subsequently analyzed the propeller's performance and the forces acting on the hull to understand how the presence or absence of the hull influenced the water jet propeller. Finally, we performed calculations and analysis of the cavitation characteristics of the coupling between the hull and the water jet propeller, considering different rotational speeds and water depths at the bottom of the pool.

The study demonstrated that the presence of the hull boundary layer under the hull-propeller coupling condition led to reduced uniformity of propeller inlet flow and lower efficiency of the propulsion pump. However, it also increased the bias toward low-flow conditions. Additionally, increasing the impeller speed led to a gradual increase in the cavitation volume within the water jet propeller, resulting in a gradual decrease in the propeller's performance.

This research provides the technical support required for effective design and operation of water jet propulsion systems. This paper involves studying and analyzing the performance and flow field of the coupling between the hull and propeller under mooring conditions with a specified hull model.

To improve the speed and accuracy of turbine blade film cooling design process, the most advanced deep learning models were introduced into this study to investigate the most suitable define for prediction work. This paper aims to create a generative surrogate model that can be applied on multi-objective optimization problems.

The latest backbone in the field of computer vision (Swin-Transformer, 2021) was introduced and improved as the surrogate function for prediction of the multi-physics field distribution (film cooling effectiveness, pressure, density and velocity). The basic samples were generated by Latin hypercube sampling method and the numerical method adopt for the calculation was validated experimentally at first. The training and testing samples were calculated at experimental conditions. At last, the surrogate model predicted results were verified by experiment in a linear cascade.

The results indicated that comparing with the Multi-Scale Pix2Pix Model, the Swin-Transformer U-Net model presented higher accuracy and computing speed on the prediction of contour results. The computation time for each step of the Swin-Transformer U-Net model is one-third of the original model, especially in the case of multi-physics field prediction. The correlation index reached more than 99.2% and the first-order error was lower than 0.3% for multi-physics field. The predictions of the data-driven surrogate model are consistent with the predictions of the computational fluid dynamics results, and both are very close to the experimental results. The application of the Swin-Transformer model on enlarging the different structure samples will reduce the cost of numerical calculations as well as experiments.

The number of U-Net layers and sample scales has a proper relationship according to equation (8). Too many layers of U-Net will lead to unnecessary nonlinear variation, whereas too few layers will lead to insufficient feature extraction. In the case of Swin-Transformer U-Net model, incorrect number of U-Net layer will reduce the prediction accuracy. The multi-scale Pix2Pix model owns higher accuracy in predicting a single physical field, but the calculation speed is too slow. The Swin-Transformer model is fast in prediction and training (nearly three times faster than multi Pix2Pix model), but the predicted contours have more noise. The neural network predicted results and numerical calculations are consistent with the experimental distribution.

This paper creates a generative surrogate model that can be applied on multi-objective optimization problems. The generative adversarial networks using new backbone is chosen to adjust the output from single contour to multi-physics fields, which will generate more results simultaneously than traditional surrogate models and reduce the time-cost. And it is more applicable to multi-objective spatial optimization algorithms. The Swin-Transformer surrogate model is three times faster to computation speed than the Multi Pix2Pix model. In the prediction results of multi-physics fields, the prediction results of the Swin-Transformer model are more accurate.

This paper aims to investigate the responses of single piles and pile groups due to tunneling-induced ground movements in a two-layered soil system. The analyses mainly focus on the additional single pile responses in terms of bending moment, lateral deflection, axial force, shaft resistance and pile settlement. Subsequently, a series of parametric studies were carried out to better understand the responses of single piles induced by tunneling. To give further understanding regarding the pile groups, a 2 × 2 pile group with two different pile head conditions, namely, free and capped, was considered.

Using the PLAXIS three-dimensional (3D) software, a full 3D numerical modeling is performed to investigate the effects of ground movements caused by tunneling on adjacent pile foundations. The numerical model was validated using centrifuge test data found in the literature. The relevance of the 3D model is also judged by comparison with the 2D plane strain model using the PLAXIS 2D code.

The numerical test results reveal that tunneling induces significant displacements and internal forces in nearby piles. The magnitude and distribution of internal forces depend mainly on the position of the pile toe relative to the tunnel depth and the distance between the pile and the vertical axis of the tunnel. As the volume loss increases from 1% to 3%, the apparent loss of pile capacity increases from 11% to 20%. By increasing the pile length from 0.5 to 1.5 times, the tunnel depth, the maximum pile settlement and lateral deflection decrease by about 63% and 18%, respectively. On the other hand, the maximum bending moment and axial load increase by about 7 and 13 times, respectively. When the pile is located at a distance of 2.5 times the tunnel diameter (Dt), the additional pile responses become insignificant. It was found that an increase in tunnel depth from 1.5Dt to 2.5Dt (with a pile length of 3Dt) increases the maximum lateral deflection by about 420%. Regarding the interaction between tunneling and group of piles, a positive group effect was observed with a significant reduction of the internal forces in rear piles. The maximum bending moment of the front piles was found to be higher than that of the rear piles by about 47%.

Soil is a complex material that shows differently in primary loading, unloading and reloading with stress-dependent stiffness. This general behavior was not possibly being accounted for in simple elastic perfectly plastic Mohr–Coulomb model which is often used to predict the behavior of soils. Thus, in the present study, the more advanced hardening soil model with small-strain stiffness (HSsmall) is used to model the non-linear stress–strain soil behavior. Moreover, unlike previous studies THAT are usually based on the assumption that the soil is homogeneous and using numerical methods by decoupled loadings under plane strain conditions; in this study, the pile responses have been exhaustively investigated in a two-layered soil system using a fully coupled 3D numerical analysis that takes into account the real interactions between tunneling and pile foundations. The paper presents a distinctive set of findings and insights that provide valuable guidance for the design and construction of shield tunnels passing through pile foundations.

The restoration and strengthening of QT600 is an industry bottleneck challenge. The Co-12 cladding layer has great wear and corrosion resistance. The purpose of this paper is to quantitatively reveal the transient evolution law of the corrosion process of Co-12 cladding layer on QT600 surface.

In this paper, a three-dimensional numerical model of the corrosion process of Co-12 cladding layer by QT600 laser cladding is established. The interaction between pitting pits and corrosion medium is considered to reveal the transient evolution of ion concentration, electrode potential, pH and corrosion rate at different locations.

The calculation shows that the ion concentration in pitting pit changes Cl⁻>Co²⁺>Na⁺, pH value decreases from top to bottom and corrosion rate at bottom is greater than that at top. The electrochemical corrosion test of Co-12 cladding layer was carried out. It is shown that the current density of QT600 increases by an order of magnitude compared to the Co-12 cladding layer, and the corrosion rate is 4.862 times higher than that of the cladding layer.

The results show that Co-12 cladding layer has great corrosion resistance, which provides an effective way for QT600 protection.

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.

A model for calculating the global overpressure time history of a single cloud detonation from overpressure time history of discrete positions in the range of single cloud detonation is to be proposed and verified. The overpressure distribution produced by multiple cloud detonation and the influence of cloud spacing and fuel mass of every cloud on the overpressure distribution are to be studied.

A calculation method is used to obtain the global overpressure field distribution after single cloud detonation from the overpressure time history of discrete distance to detonation center after single cloud detonation. On this basis, the overpressure distribution produced by multi-cloud under different cloud spacing and different fuel mass conditions is obtained.

The results show that for 150 kg fuel, when the spacing of three clouds is 40 m, 50 m, respectively, the overpressure range of larger than 0.1 MPa is 5496.48 mˆ2 and 6235.2 mˆ2, which is 2.89 times and 3.28 times of that of single cloud detonation. The superposition effect can be ignored when the spacing between the three clouds is greater than 60 m. In the case of fixed cloud spacing, once the overpressure forms continuous effective superposition, the marginal utility of fuel decreases.

A model for calculating the global overpressure time history of a single cloud detonation from overpressure time history of discrete positions in the range of single cloud detonation is proposed and verified. Based on this method, the global overpressure field of single cloud detonation is reconstructed, and the superimposed overpressure distribution characteristics of three cloud detonation are calculated and analyzed.