Search results

This paper aims to construct positivity-preserving finite volume schemes for the three-dimensional convection–diffusion equation that are applicable to arbitrary polyhedral grids.

The cell vertices are used to define the auxiliary unknowns, and the primary unknowns are defined at cell centers. The diffusion flux is discretized by the classical nonlinear two-point flux approximation. To ensure the fully discrete scheme has positivity-preserving property, an improved discretization method for the convection flux was presented. Besides, a new positivity-preserving vertex interpolation method is derived from the linear reconstruction in the discretization of convection flux. Moreover, the Picard iteration method may have slow convergence in solving the nonlinear system. Thus, the Anderson acceleration of Picard iteration method is used to solve the nonlinear system. A condition number monitor of matrix is employed in the Anderson acceleration method to achieve better robustness.

The new scheme is applicable to arbitrary polyhedral grids and has a second-order accuracy. The results of numerical experiments also confirm the positivity-preserving of the discretization scheme.

1. This article presents a new positivity-preserving finite volume scheme for the 3D convection–diffusion equation. 2. The new discretization scheme of convection flux is constructed. 3. A new second-order interpolation algorithm is given to eliminate the auxiliary unknowns in flux expressions. 4. An improved Anderson acceleration method is applied to accelerate the convergence of Picard iterations. 5. This scheme can solve the convection–diffusion equation on the distorted meshes with second-order accuracy.

The purpose of this study is analysis of deformation and vibrations of turbine blades produced by high electrolyte pressure during electrochemical machining.

An experimental setup was designed, experiments were conducted and the obtained results were compared with the finite element results. The deformations were measured according to various flow rates of electrolyte. In finite element calculations, the pressure distribution created by the electrolyte on the blade surface was obtained in the ANSYS® (A finite element analysis software) Fluent software and transferred to the static structural where the deformation analysis was carried out. Three different parameters were examined, namely blade thickness, blade material and electrolyte pressure on blade disk caused by mass flow rate. The deformation results were compared with the gap distances between cathode and anode.

Large deformations were obtained at the free end of the blade and the most curved part of it. The appropriate pressure values for the electrolyte to be used in the production of blisk blades were proposed numerically. It has been determined that high pressure applications are not suitable for gap distance lower than 0.5 mm.

When the literature is examined, it is required that the high speed flow of the electrolyte is desired in order to remove the parts that are separated from the anode from the machining area during electrochemical machining. However, the electrolyte flowing at high speeds causes high pressure in the blisk blades, excessive deformation and vibration of the machined part, and as a result, contact of the anode with the cathode. This study provides important findings for smooth electro chemical machining at high electrolyte flows.

This work measures the performance characteristics of a hemispherical resonator gyroscope (HRG) and compares it with a numerical model.

This work we explore the optical and piezoelectric measurement methods to determine the resonant frequency of HRG. These experimental results are compared with their numerically obtained values. To explore the performance characteristics, the effect of varying actuation voltages on the sense mode displacement and the piezoelectric sensor output was studied in the absence of input angular rate. The structure was then subjected to range of angular rate signals, at a constant actuation voltage and the corresponding sensor response was analysed.

Experimental values of the resonant frequencies in drive and sense modes are found to be within 8% of the numerical results. The sensor output depicts a quadratic dependency on the applied angular rate, which is synchronous with the governing equations of the HRG. The experimental output is within 12% of that obtained numerically. The sensor is found to resolve upto 0.24 rad/s.

This work presents an in-house developed inexpensive measurement setup for static and dynamic characterization of mesoscale MEMS gyroscopes. The measurement setup can also be modified accordingly for measurement of other MEMS-based devices.

The purpose of this paper is to present a new approach to optimizing the design of 3D magnetic circuits. This approach is based on topology optimization, where derivative calculations are performed using the continuous adjoint method. Thus, the continuous adjoint method for magnetostatics has to be developed in 3D and has to be combined with penalization, filtering and homotopy approaches to provide an efficient optimization code.

To provide this new topology optimization code, this study starts from 2D magnetostatic results to perform the sensitivity analysis, and this approach is extended to 3D. From this sensitivity analysis, the continuous adjoint method is derived to compute the gradient of an objective function of a 3D topological optimization design problem. From this result, this design problem is discretized and can then be solved by finite element software. Thus, by adding the solid isotropic material with penalization (SIMP) penalization approach and developing a homotopy-based optimization algorithm, an interesting means for designing 3D magnetic circuits is provided.

In this paper, the 3D continuous adjoint method for magnetostatic problems involving an objective least-squares function is presented. Based on 2D results, new theoretical results for developing sensitivity analysis in 3D taking into account different parameters including the ferromagnetic material, the current density and the magnetization are provided. Then, by discretizing, filtering and penalizing using SIMP approaches, a topology optimization code has been derived to address only the ferromagnetic material parameters. Based on this efficient gradient computation method, a homotopy-based optimization algorithm for solving large-scale 3D design problems is developed.

In this paper, an approach based on topology optimization to solve 3D magnetostatic design problems when an objective least-squares function is involved is proposed. This approach is based on the continuous adjoint method derived for 3D magnetostatic design problems. The effectiveness of this topology optimization code is demonstrated by solving the design of a 3D magnetic circuit with up to 100,000 design variables.

As a new type of spinning machine, the jet spinning machine absorbs the carding system of the rotating cup spinning series and the nozzle part of the jet spinning. This paper aims to intends to introduce the double carding structure currently studied by the rotating cup spinning into the jet spinning machine, and analyze the influence of the nozzle characteristic number on the flow field in the double carding structure to verify the advantages of the double carding structure.

The simulation is used to evaluate the performance of single/double split jet spinning and nozzle feature number, verify the technical advantages of double split jet spinning and evaluate the influence of nozzle feature number on flow field. The influence of the nozzle characteristic number on the flow pattern in the four models is compared. The advantages and disadvantages of a conventional single comb and a double comb with a bypass channel on the longer side of the transport channel as an additional air supply channel are also evaluated.

At present, the double comb technology of rotary cup spinning is being studied at home and abroad to improve the spinning quality and improve the difficult problem of mixed yarn with large difference in processing fiber properties. At present, the jet spinning machine combines the advantages of rotary cup spinning and jet spinning, absorbing the comb system of rotary cup spinning series and the nozzle part of jet spinning. Therefore, it is found that the introduction of the double-split structure into the wool jet spinning has research value to improve the spinning quality.

The purpose of this paper is to refer to the previous research on the double comb structure in rotary spinning, and to apply the double comb structure in the new jet spinning machine to improve the spinning quality. The simulation is used to evaluate the performance of single/double split jet spinning and nozzle feature number, verify the technical advantages of double split jet spinning and evaluate the influence of nozzle feature number on flow field.

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.

The purpose of this study is to compare the thermal performance of two flow configurations in corrugated plate heat exchanger (CPHE): vertical flow configuration (CPHE_vert.) and diagonal flow configuration (CPHE_diag.). The study aims to determine the differences between these configurations and evaluate their respective thermal performance based on metrics such as heat transfer rates, pressure drop values and flow distribution.

The study compares the thermal performance of two flow arrangements of CPHE using identical geometrical dimensions and test conditions. Computational fluid dynamics (CFD) is employed, and a validated numerical model is used for the investigation. The comparison is based on analyzing the rate of heat transfer and pressure drop data between the two flow arrangements.

The findings indicate that the diagonal flow configuration in CPHEs offers improved flow distribution, enhanced heat transfer performance and lower pressure drop compared to the vertical flow configuration. However, the differences in general in the thermal performance of CPHE_vert. and CPHE_diag. are found to be minimal.

To the best of the author’s knowledge, this study represents the first attempt to investigate the impact of vertical and diagonal flow configurations on the thermal performance of the CPHE.

This study aims to investigate the erosion wear rate of a stainless steel automobile exhaust manifold, both computationally and physically.

The experiment was performed on a motorcycle exhaust manifold as well as on a 3D model, created using SolidWorks 2022 CAD software. The analysis was later achieved using ANSYS 19.2 simulation software using Fluent – code.

The analysis of solid particle erosion in the exhaust manifold revealed that erosion wear is concentrated predominantly at the extrados of the manifold, with the most significant wear occurring at the lowermost bend. The erosion wear rate increases with larger particulate sizes and varies among bends, with negligible wear observed in straight pipes. The SEM analysis further confirmed surface degradation, with rugged textures, pits and grooves indicating abrasive wear. Spine-like structures and fractured soot particles suggest erosive and abrasive forces caused by high-speed contact of exhaust gas compounds. Energy dispersive X-ray spectroscopy revealed significant carbon abundance, indicating carbonaceous compounds from fuel combustion, along with notable amounts of oxygen and iron, typical of oxidized metallic constituents. The discrete phase modeling (DPM) analysis highlighted peak particulate matter deposition at the first bend exit, with maximum concentrations observed at specific angles. This deposition is influenced by centrifugal force, leading to increased PM concentration at outer bend walls. Velocity magnitude contours showed asymmetrical flow profiles, with high turbulence levels and secondary flow induced by centrifugal effects in bend areas. Dynamic pressure contours revealed varying pressures at intrados and extrados, with maximum pressure observed at the intrados of the manifold’s bends. These findings provide valuable insights into erosion wear, particulate dispersion and flow dynamics within the exhaust manifold.

The study investigated an automobile exhaust manifold model using ANSYS Fluent code and DPM to analyze erosion wear rate phenomena and its various constituents. This analysis was conducted in comparison with a physically eroded sample. The study offers insights into the mechanism underlying the exhaust manifold of an automobile.

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 propose and numerically assess different ways of discretising a very weak formulation of the Poisson problem.

We use integration by parts twice to shift smoothness requirements to the test functions, thereby allowing low-regularity data and solutions.

Various conforming discretisations are presented and tested, with numerical results indicating good accuracy and stability in different types of problems.

This is one of the first articles to propose and test concrete discretisations for very weak variational formulations in primal form. The numerical results, which include a problem based on real MRI data, indicate the potential of very weak finite element methods for tackling problems with low regularity.