Search results

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.

The paper aims to expand the scope of application of the lattice Boltzmann method (LBM), especially in the field of aircraft engineering. The traditional LBM is usually applied to incompressible flows at a low Reynolds number, which is not sufficient to satisfy the needs of aircraft engineering. Devoted to tackling the defect, the paper proposes a developed LBM combining the subgrid model and the multiple relaxation time (MRT) approach. A multilayer adaptive Cartesian grid method to improve the computing efficiency of the traditional LBM is also employed.

The subgrid model and the multilayer adaptive Cartesian grid are introduced into MRT-LBM for simulations of incompressible flows at a high Reynolds number. Validated by several typical flow simulations, the numerical methods in this paper can efficiently study the flows under high Reynolds numbers.

Some numerical simulations for the lid-driven flow of cavity, flow around iced GLC305, LB606b and ONERA-M6 are completed. The paper presents the investigation results, indicating that the methods are accurate and effective for the separated flow after icing.

LBM is developed with the addition of the subgrid model and the MRT method. A numerical strategy is proposed using a multilayer adaptive Cartesian grid method and its treatment of boundary conditions. The paper refers to innovative algorithm developments and applications to the aircraft engineering, especially for iced wing simulations with flow separations.

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.

Identifying the flow regime is a prerequisite for accurately modeling two-phase flow. This paper aims to introduce a comprehensive data-driven workflow for flow regime identification.

A numerical two-phase flow model was validated against experimental data and was used to generate dynamic pressure signals for three different flow regimes. First, four distinct methods were used for feature extraction: discrete wavelet transform (DWT), empirical mode decomposition, power spectral density and the time series analysis method. Kernel Fisher discriminant analysis (KFDA) was used to simultaneously perform dimensionality reduction and machine learning (ML) classification for each set of features. Finally, the Shapley additive explanations (SHAP) method was applied to make the workflow explainable.

The results highlighted that the DWT + KFDA method exhibited the highest testing and training accuracy at 95.2% and 88.8%, respectively. Results also include a virtual flow regime map to facilitate the visualization of features in two dimension. Finally, SHAP analysis showed that minimum and maximum values extracted at the fourth and second signal decomposition levels of DWT are the best flow-distinguishing features.

This workflow can be applied to opaque pipes fitted with pressure sensors to achieve flow assurance and automatic monitoring of two-phase flow occurring in many process industries.

This paper presents a novel flow regime identification method by fusing dynamic pressure measurements with ML techniques. The authors’ novel DWT + KFDA method demonstrates superior performance for flow regime identification with explainability.

Gas metal arc-based directed energy deposition (GMA-DED) process experiences residual stress (RS) developed due to heat accumulation during successive layer deposition as a significant challenge. To address that, monitoring of transient temperature distribution concerning time is a critical input. Finite element analysis (FEA) is considered a decisive engineering tool in quantifying temperature and RS in all manufacturing processes. However, computational time and prediction accuracy has always been a matter of concern for FEA-based prediction of responses in the GMA-DED process. Therefore, this study aims to investigate the effect of finite element mesh variations on the developed RS in the GMA-DED process.

The variation in the element shape functions, i.e. linear- and quadratic-interpolation elements, has been used to model a single-track 10-layered thin-walled component in Ansys parametric design language. Two cases have been proposed in this study: Case 1 has been meshed with the linear-interpolation elements and Case 2 has been meshed with the combination of linear- and quadratic-interpolation elements. Furthermore, the modelled responses are authenticated with the experimental results measured through the data acquisition system for temperature and RS.

A good agreement of temperature and RS profile has been observed between predicted and experimental values. Considering similar parameters, Case 1 produced an average error of 4.13%, whereas Case 2 produced an average error of 23.45% in temperature prediction. Besides, comparing the longitudinal stress in the transverse direction for Cases 1 and 2 produced an error of 8.282% and 12.796%, respectively.

To avoid the costly and time-taking experimental approach, the experts have suggested the utilization of numerical methods in the design optimization of engineering problems. The FEA approach, however, is a subtle tool, still, it faces high computational cost and low accuracy based on the choice of selected element technology. This research can serve as a basis for the choice of element technology which can predict better responses in the thermo-mechanical modelling of the GMA-DED process.

The purpose of this research is to efficiently separate incident terahertz (THz) waves into distinct transmission and reflection channels by minimizing the absorption ratio. So, the optical systems operating within the THz frequency range can developed. To achieve a multi-band response, four different periodic arrays of graphene patterns are used. These arrays are strategically stacked on both sides of three SU-8 photoresists, serving as dielectric materials. Consequently, each layer exhibits a unique influence on the device's response, and by applying four external bias voltages, the behavior of the device can be precisely controlled and adjusted.

A novel optoelectronic device operating in the THz frequency range is introduced, using periodic arrays of graphene patterns and SU-8 photoresist dielectrics. The design of this device is based on meta-surface principles, using both the equivalent circuit model (ECM) and transmission line concept. The output of the device is a THz coupler implemented by analyzing the reflection and transmission channels. The structure is characterized using the ECM and validated through comprehensive full-wave simulations. By representing the electromagnetic phenomenon with passive circuit elements, enabling the calculation of absorption, reflection and transmission through the application of the theory of maximum power transfer.

Based on simulation results and theoretical analysis, the proposed device exhibits sensitivity to gate biasing, enabling efficient reflection and transmission of THz waves. The device achieves reflection and transmission peaks exceeding across the five distinct THz bands 90%, and its behavior can be tuned by external gate biasing. Moreover, the device's sensitivity to variations in geometrical parameters and chemical potentials demonstrates its reliable performance. With its outstanding performance, this high-performance meta-surface emerges as an ideal candidate for fundamental building blocks in larger optical systems, including sensors and detectors, operating within the THz frequency band.

The proposed device covers a significant portion of the THz gap through the provision of five adjustable peaks for reflection and transmission channels. Additionally, the ECM and impedance matching concept offers a simplified and time-efficient approach to designing the meta-surface. Leveraging this approach, the proposed device is effectively represented using passive circuit elements such as inductors, capacitors and resistors, while its performance is validated through the utilization of the finite element method (FEM) as a full-wave simulation tool. This combination of circuit modeling and FEM simulation contributes to the robustness and accuracy of the device's performance evaluation.

The aerodynamic differences between the head car (HC) and tail car (TC) of a high-speed maglev train are significant, resulting in control difficulties and safety challenges in operation. The arch structure has a significant effect on the improvement of the aerodynamic lift of the HC and TC of the maglev train. Therefore, this study aims to investigate the effect of a streamlined arch structure on the aerodynamic performance of a 600 km/h maglev train.

Three typical streamlined arch structures for maglev trains are selected, i.e. single-arch, double-arch and triple-arch maglev trains. The vortex structure, pressure of train surface, boundary layer, slipstream and aerodynamic forces of the maglev trains with different arch structures are compared by adopting improved delayed detached eddy simulation numerical calculation method. The effects of the arch structures on the aerodynamic performance of the maglev train are analyzed.

The dynamic topological structure of the wake flow shows that a change in arch structure can reduce the vortex size in the wake region; the vortex size with double-arch and triple-arch maglev trains is reduced by 15.9% and 23%, respectively, compared with a single-arch maglev train. The peak slipstream decreases with an increase in arch structures; double-arch and triple-arch maglev trains reduce it by 8.89% and 16.67%, respectively, compared with a single-arch maglev train. The aerodynamic force indicates that arch structures improve the lift imbalance between the HC and TC of a maglev train; double-arch and triple-arch maglev trains improve it by 22.4% and 36.8%, respectively, compared to a single-arch maglev train.

This study compares the effects of a streamlined arch structure on a maglev train and its surrounding flow field. The results of the study provide data support for the design and safe operation of high-speed maglev trains.

This article aims to present a direct solution to handle linear constraints in finite element (FE) analysis without penalties or the Lagrange multipliers introduced.

First, the system of linear equations corresponding to the linear constraints is solved for the leading variables in terms of the free variables and the constants. Then, the reduced system of equilibrium equations with respect to the free variables is derived from the finite-dimensional virtual work equation. Finally, the algorithm is designed.

The proposed procedure is promising in three typical cases: (1) to enforce displacement constraints in any direction; (2) to implement local refinements by allowing hanging nodes from element subdivision and (3) to treat non-matching grids of distinct parts of the problem domain. The procedure is general and suitable for 3D non-linear analyses.

The algorithm is fitted only to the Galerkin-based numerical methods.

The proposed procedure does not need Lagrange multipliers or penalties. The tangential stiffness matrix of the reduced system of equilibrium equations reserves positive definiteness and symmetry. Besides, many contemporary Galerkin-based numerical methods need to tackle the enforcement of the essential conditions, whose weak forms reduce to linear constraints. As a result, the proposed procedure is quite promising.

The purpose of this paper is to design a simple and accurate a-posteriori Lagrangian-based error estimator is developed for the class of backward differentiation formula (BDF) algorithms with variable time step size, and the adaptive time-stepping in BDF algorithms is demonstrated for efficient time-dependent simulations in fluid flow and heat transfer.

The Lagrange interpolation polynomial is used to predict the time derivative, and then the accurate primary result is obtained by the Gauss integral, which is applied to evaluate the local error. Not only the generalized formula of the proposed error estimator is presented but also the specific expression for the widely applied BDF1/2/3 is illustrated. Two essential executable MATLAB functions to implement the proposed error estimator are appended for practical applications. Then, the adaptive time-stepping is demonstrated based on the newly proposed error estimator for BDF algorithms.

The validation tests show that the newly proposed error estimator is accurate such that the effectivity index is always close to unity for both linear and nonlinear problems, and it avoids under/overestimation of the exact local error. The applications for fluid dynamics and coupled fluid flow and heat transfer problems depict the advantage of adaptive time-stepping based on the proposed error estimator for time-dependent simulations.

In contrast to existing error estimators for BDF algorithms, the present work is more accurate for the local error estimation, and it can be readily extended to practical applications in engineering with a few changes to existing codes, contributing to efficient time-dependent simulations in fluid flow and heat transfer.

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.