Search results

The purpose of this paper is to present an automatic approach for mesh sizing field generation of complicated  computer-aided design (CAD) models.

In this paper, the authors present an automatic approach for mesh sizing field generation. First, a source point extraction algorithm is applied to capture curvature and proximity features of CAD models. Second, according to the distribution of feature source points, an octree background mesh is constructed for storing element size value. Third, mesh size value on each node of background mesh is calculated by interpolating the local feature size of the nearby source points, and then, an initial mesh sizing field is obtained. Finally, a theoretically guaranteed smoothing algorithm is developed to restrict the gradient of the mesh sizing field.

To achieve high performance, the proposed approach has been implemented in multithreaded parallel using OpenMP. Numerical results demonstrate that the proposed approach is remarkably efficient to construct reasonable mesh sizing field for complicated CAD models and applicable for generating geometrically adaptive triangle/tetrahedral meshes. Moreover, since the mesh sizing field is defined on an octree background mesh, high-efficiency query of local size value could be achieved in the following mesh generation procedure.

How to determine a reasonable mesh size for complicated CAD models is often a bottleneck of mesh generation. For the complicated models with thousands or even ten thousands of geometric entities, it is time-consuming to construct an appropriate mesh sizing field for generating high-quality mesh. A parallel algorithm of mesh sizing field generation with low computational complexity is presented in this paper, and its usability and efficiency have been verified.

Hexahedral meshing is one of the most important steps in performing an accurate simulation using the finite element analysis (FEA). However, the current hexahedral meshing method in the industry is a nonautomatic and inefficient method, i.e. manually decomposing the model into suitable blocks and obtaining the hexahedral mesh from these blocks by mapping or sweeping algorithms. The purpose of this paper is to propose an almost automatic decomposition algorithm based on the 3D frame field and model features to replace the traditional time-consuming and laborious manual decomposition method.

The proposed algorithm is based on the 3D frame field and features, where features are used to construct feature-cutting surfaces and the 3D frame field is used to construct singular-cutting surfaces. The feature-cutting surfaces constructed from concave features first reduce the complexity of the model and decompose it into some coarse blocks. Then, an improved 3D frame field algorithm is performed on these coarse blocks to extract the singular structure and construct singular-cutting surfaces to further decompose the coarse blocks. In most modeling examples, the proposed algorithm uses both types of cutting surfaces to decompose models fully automatically. In a few examples with special requirements for hexahedral meshes, the algorithm requires manual input of some user-defined cutting surfaces and constructs different singular-cutting surfaces to ensure the effectiveness of the decomposition.

Benefiting from the feature decomposition and the 3D frame field algorithm, the output blocks of the proposed algorithm have no inner singular structure and are suitable for the mapping or sweeping algorithm. The introduction of internal constraints makes 3D frame field generation more robust in this paper, and it can automatically correct some invalid 3–5 singular structures. In a few examples with special requirements, the proposed algorithm successfully generates valid blocks even though the singular structure of the model is modified by user-defined cutting surfaces.

The proposed algorithm takes the advantage of feature decomposition and the 3D frame field to generate suitable blocks for a mapping or sweeping algorithm, which saves a lot of simulation time and requires less experience. The user-defined cutting surfaces enable the creation of special hexahedral meshes, which was difficult with previous algorithms. An improved 3D frame field generation method is proposed to correct some invalid singular structures and improve the robustness of the previous methods.

To present the detailed implementation processes of the IDEAL algorithm for two-dimensional compressible flows based on Delaunay triangular mesh, and compare the performance of the SIMPLE and IDEAL algorithms for solving compressible problems. What’s more, the implementation processes of Delaunay mesh generation and derivation of the pressure correction equation are also introduced.

Programming completely in C++.

Five compressible examples are used to test the SIMPLE and IDEAL algorithms, and the comparison with measurement data shows good agreement. The IDEAL algorithm has much better performance in both convergence rate and stability over the SIMPLE algorithm.

The detail solution procedure of implementing the IDEAL algorithm for compressible flows based on Delaunay triangular mesh is presented in this work, seemingly first in the literature.

This study aims to investigate the noise-inducing characteristics during the start-up process of a mixed-flow pump and the impact of different start-up schemes on pump noise.

This study conducted numerical simulations on the mixed-flow pump under different start-up schemes and investigated the flow characteristics and noise distribution under these schemes.

The results reveal that the dipole noise is mainly caused by pressure fluctuations, while the quadrupole noise is mainly generated by the generation, development and breakdown of vortices. Additionally, the noise evolution characteristics during the start-up process of the mixed-flow pump can be divided into the initial stage, stable growth stage, impulse stage and stable operation stage.

The findings of this study can provide a theoretical basis for the selection of start-up schemes for mixed-flow pumps, reducing flow noise and improving the operational stability of mixed-flow pumps.

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.

The purpose of this study is to investigate the influence of enclosure shape on magnetohydrodynamic (MHD) nanofluidic flow, heat transfer and irreversibility in square, trapezoidal and triangular thermal systems under fluid volume constraints, with the aim of optimizing thermal behavior in diverse applications.

The study uses numerical simulations based on a finite element-based technique to analyze the effects of the Rayleigh number (Ra), Hartmann number (Ha), magnetic field orientation (γ) and nanoparticle concentration (ζ) on heat transfer characteristics and thermodynamic entropy production.

The key findings reveal that the geometrical design significantly influences fluid velocity, heat transfer and irreversibility. Trapezoidal thermal systems outperform square systems, while triangular systems achieve optimal enhancement. Nanoparticle concentration enhances heat transfer and flow strength at higher Rayleigh numbers. The magnetic field intensity has a significant impact on fluid flow and heat transport in natural convection, with higher Hartmann numbers resulting in reduced flow strength and heat transfer. The study also highlights the influence of various parameters on thermodynamic entropy production.

Further research can explore additional geometries, parameters and boundary conditions to expand the understanding of enclosure shape effects on MHD nanofluidic flow and heat transfer. Experimental validation can complement the numerical simulations presented in this study.

This study provides valuable insights into the impact of enclosure shape on heat transfer performance in MHD nanofluid flow systems. The findings contribute to the optimization of thermal behavior in applications such as electronics cooling and energy systems. The comparison of different enclosure shapes and the analysis of thermodynamic entropy production add novelty to the study.

Topology optimization of structures under self-weight loading is a challenging problem which has received increasing attention in the past years. The use of standard formulations based on compliance minimization under volume constraint suffers from numerous difficulties for self-weight dominant scenarios, such as non-monotonic behaviour of the compliance, possible unconstrained character of the optimum and parasitic effects for low densities in density-based approaches. This paper aims to propose an alternative approach for dealing with topology design optimization of structures into three spatial dimensions subject to self-weight loading.

In order to overcome the above first two issues, a regularized formulation of the classical compliance minimization problem under volume constraint is adopted, which enjoys two important features: (a) it allows for imposing any feasible volume constraint and (b) the standard (original) formulation is recovered once the regularizing parameter vanishes. The resulting topology optimization problem is solved with the help of the topological derivative method, which naturally overcomes the above last issue since no intermediate densities (grey-scale) approach is necessary.

A novel and simple approach for dealing with topology design optimization of structures into three spatial dimensions subject to self-weight loading is proposed. A set of benchmark examples is presented, showing not only the effectiveness of the proposed approach but also highlighting the role of the self-weight loading in the final design, which are: (1) a bridge structure is subject to pure self-weight loading; (2) a truss-like structure is submitted to an external horizontal force (free of self-weight loading) and also to the combination of self-weight and the external horizontal loading; and (3) a tower structure is under dominant self-weight loading.

An alternative regularized formulation of the compliance minimization problem that naturally overcomes the difficulties of dealing with self-weight dominant scenarios; a rigorous derivation of the associated topological derivative; computational aspects of a simple FreeFEM implementation; and three-dimensional numerical benchmarks of bridge, truss-like and tower structures.

This study aims to investigate the fabrication of solid ankle foot orthoses (SAFOs) using fused deposition modeling (FDM) printing technology. It emphasizes cost-effective 3D scanning with the Kinect sensor and conducts a comparative analysis of SAFO durability with varying thicknesses and materials, including polylactic acid (PLA) and carbon fiber-reinforced (PLA-C), to address research gaps from prior studies.

In this study, the methodology comprises key components: data capture using a cost-effective Microsoft Kinect® Xbox 360 scanner to obtain precise leg dimensions for SAFOs. SAFOs are designed using CAD tools with varying thicknesses (3, 4, and 5 mm) while maintaining consistent geometry, allowing controlled thickness impact investigation. Fabrication uses PLA and PLA-C materials via FDM 3D printing, providing insights into material suitability. Mechanical analysis uses dual finite element analysis to assess force–displacement curves and fracture behavior, which were validated through experimental testing.

The results indicate that the precision of the scanned leg dimensions, compared to actual anthropometric data, exhibits a deviation of less than 5%, confirming the accuracy of the cost-effective scanning approach. Additionally, the research identifies optimal thicknesses for SAFOs, recommending a 4 and 5 mm thickness for PLA-C-based SAFOs and an only 5 mm thickness for PLA-based SAFOs. This optimization enhances the overall performance and effectiveness of these orthotic solutions.

This study’s innovation lies in its holistic approach, combining low-cost 3D scanning, 3D printing and computational simulations to optimize SAFO materials and thickness. These findings advance the creation of cost-effective and efficient orthotic solutions.

Constructing porous wind barriers is one of the most effective approaches to increase the running safety of trains on viaducts in crosswinds. This paper aims to further improve the wind-sheltering performance of the porous wind barriers.

Improved delayed detached eddy simulations based on the k-ω turbulence model were carried out, and the results were validated with wind tunnel tests. The effects of the hole diameter on the flow characteristics and wind-sheltering performance were studied by comparing the wind barriers with the porosity of 21.6% and the hole diameters of 60 mm–360 mm. The flow characteristics above the windward and leeward tracks were analyzed, and the wind-sheltering performance of the wind barriers was assessed using the wind speed reduction coefficients.

The hole diameters affected the jet behind the wind barriers and the recirculation region above the tracks. Below the top of the wind barriers, the time-averaged velocity first decreased and then increased with the increase in the hole diameter. The wind barrier with the hole diameter of 120 mm had the best wind-sheltering performance for the windward track, but such barrier might lead to overprotection on the leeward track. The wind-sheltering performance of the wind barriers with the hole diameters of 240 mm and 360 mm was significantly degraded, especially above the windward track.

The effects of the hole diameters on the wake and wind-sheltering performance of the wind barriers were studied, by which the theoretical basis is provided for a better design of the porous wind barrier.

The purpose of this paper is to calculate the meshing stiffness of nutation face gear considering the roughness, establish the calculation method of time-varying meshing stiffness of rough tooth surface and analyze the influence of roughness, load and other factors on the meshing stiffness of tooth surface.

The Weierstrass–Mandelbrot (W-M) function in the Majumdar–Bhushan model is used to characterize the rough contact line of the tooth surface, the normal height and radius of the micro convex body are calculated and the contact flexibility of the contact point of the tooth surface is obtained. The contact flexibility and the bending shear deformation flexibility obtained previously are substituted into the improved deformation compatibility equation for iterative calculation, and the time-varying meshing stiffness of the nutation face gear considering the roughness is obtained.

Compared with ABAQUS finite element simulation results, it is found that the meshing stiffness curve of rough tooth surface is more gentle than that of smooth tooth surface, the meshing stiffness value is smaller and the meshing stiffness change is smaller at the position where the number of gear teeth coincide changes.

In the process of calculating contact deformation, the fractal theory W-M function is used to characterize the contact line of the rough nutation face gear, and the deformation coordination condition considering roughness is improved. Therefore, the method of time-varying meshing stiffness considering roughness can obtain more accurate results, which provides theory and data for the subsequent dynamics analysis of the nutation face gear transmission.