Search results

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.

Understanding and tailoring the solidification characteristics and microstructure evolution in as-built parts fabricated by laser powder bed fusion (LPBF) is crucial as they influence the final properties. Experimental approaches to address this issue are time and capital-intensive. This study aims to develop an efficient numerical modeling approach to develop the process–structure (P-S) linkage for LPBF-processed Inconel 718.

In this study, a numerical approach based on the finite element method and cellular automata was used to model the multilayer, multitrack LPBF build for predicting the solidification characteristics (thermal gradient G and solidification rate R) and the average grain size. Validations from published experimental studies were also carried out to ensure the reliability of the proposed numerical approach. Furthermore, microstructure simulations were used to develop P-S linkage by evaluating the effects of key LPBF process parameters on G × R, G/R and average grain size. A solidification or G-R map was also developed to comprehend the P-S linkage.

It was concluded from the developed G-R map that low laser power and high scan speed will result in a finer microstructure due to an increase in G × R, but due to a decrease in G/R, columnar characteristics are also reduced. Moreover, increasing the layer thickness and decreasing the hatch spacing lowers the G × R, raises the G/R and generates a coarse columnar microstructure.

The proposed numerical modeling approach was used to parametrically investigate the effect of LPBF parameters on the resulting microstructure. A G-R map was also developed that enables the tailoring of the as-built LPBF microstructure through solidification characteristics by tuning the process parameters.

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.

This study aims to investigate the possible defects and their root causes in a soft-termination multilayered ceramic capacitor (MLCC) when subjected to a thermal reflow process.

Specimens of the capacitor assembly were subjected to JEDEC level 1 preconditioning (85 °C/85%RH/168 h) with 5× reflow at 270°C peak temperature. Then, they were inspected using a 2 µm scanning electron microscope to investigate the evidence of defects. The reliability test was also numerically simulated and analyzed using the extended finite element method implemented in ABAQUS.

Excellent agreements were observed between the SEM inspections and the simulation results. The findings showed evidence of discontinuities along the Cu and the Cu-epoxy layers and interfacial delamination crack at the Cu/Cu-epoxy interface. The possible root causes are thermal mismatch between the Cu and Cu-epoxy layers, moisture contamination and weak Cu/Cu-epoxy interface. The maximum crack length observed in the experimentally reflowed capacitor was measured as 75 µm, a 2.59% difference compared to the numerical prediction of 77.2 µm.

This work's contribution is expected to reduce the additional manufacturing cost and lead time in investigating reliability issues in MLCCs.

Despite the significant number of works on the reliability assessment of surface mount capacitors, work on crack growth in soft-termination MLCC is limited. Also, the combined experimental and numerical investigation of reflow-induced reliability issues in soft-termination MLCC is limited. These cited gaps are the novelties of this study.

Present study aims to extend the lattice Boltzmann method (LBM) to simulate radiation in geometries with curved boundaries, as the first step to simulate radiation in complex porous media. In recent years, researchers have increasingly explored the use of porous media to improve the heat transfer processes. The lattice Boltzmann method (LBM) is one of the most effective techniques for simulating heat transfer in such media. However, the application of the LBM to study radiation in complex geometries that contain curved boundaries, as found in many porous media, has been limited.

The numerical evaluation of the effect of the radiation-conduction parameter and extinction coefficient on temperature and incident radiation distributions demonstrates that the proposed LBM algorithm provides highly accurate results across all cases, compared to those found in the literature or those obtained using the finite volume method (FVM) with the discrete ordinates method (DOM) for radiative information.

For the case with a conduction-radiation parameter equal to 0.01, the maximum relative error is 1.9% in predicting temperature along vertical central line. The accuracy improves with an increase in the conduction-radiation parameter. Furthermore, the comparison between computational performances of two approaches reveals that the LBM-LBM approach performs significantly faster than the FVM-DOM solver.

The difficulty of radiative modeling in combined problems involving irregular boundaries has led to alternative approaches that generally increase the computational expense to obtain necessary radiative details. To address the limitations of existing methods, this study presents a new approach involving a coupled lattice Boltzmann and first-order blocked-off technique to efficiently model conductive-radiative heat transfer in complex geometries with participating media. This algorithm has been developed using the parallel lattice Boltzmann solver.

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.

This study aims to investigate the morphing concepts able to manipulate the dynamics of the downstream unsteadiness in the separated shear layers and, in the wake, be able to modify the upstream shock–boundary layer interaction (SBLI) around an A320 morphing prototype to control these instabilities, with emphasis to the attenuation or even suppression of the transonic buffet. The modification of the aerodynamic performances according to a large parametric study carried out at Reynolds number of 4.5 × 10⁶, Mach number of 0.78 and various angles of attack in the range of (0, 2.4)° according to two morphing concepts (travelling waves and trailing edge vibration) are discussed, and the final benefits in aerodynamic performance increase are evaluated.

This article examines through high fidelity (Hi-Fi) numerical simulation the effects of the trailing edge (TE) actuation and of travelling waves along a specific area of the suction side starting from practically the most downstream position of the shock wave motion according to the buffet and extending up to nearly the TE. The present paper studies through spectral analysis the coherent structures development in the near wake and the comparison of the aerodynamic forces to the non-actuated case. Thus, the physical mechanisms of the morphing leading to the increase of the lift-to-drag ratio and the drag and noise sources reduction are identified.

This study investigates the influence of shear-layer and near-wake vortices on the SBLI around an A320 aerofoil and attenuation of the related instabilities thanks to novel morphing: travelling waves generated along the suction side and trailing-edge vibration. A drag reduction of 14% and a lift-to-drag increase in the order of 8% are obtained. The morphing has shown a lift increase in the range of (1.8, 2.5)% for angle of attack of 1.8° and 2.4°, where a significant lift increase of 7.7% is obtained for the angle of incidence of 0° with a drag reduction of 3.66% yielding an aerodynamic efficiency of 11.8%.

This paper presents results of morphing A320 aerofoil, with a chord of 70cm and subjected to two actuation kinds, original in the state of the art at M = 0.78 and Re = 4.5 million. These Hi-Fi simulations are rather rare; a majority of existing ones concern smaller dimensions. This study showed for the first time a modified buffet mode, displaying periodic high-lift “plateaus” interspersed by shorter lift-decrease intervals. Through trailing-edge vibration, this pattern is modified towards a sinusoidal-like buffet, with a considerable amplitude decrease. Lock-in of buffet frequency to the actuation is obtained, leading to this amplitude reduction and a drastic aerodynamic performance increase.

Realistic composite vehicles with 2, 3, 5 and 9 axles, consisting of a truck with one or two trailers, are addressed in this paper by computational models for vehicle–bridge interaction analysis.

The vehicle–bridge interaction (VBI) models are formed by sets of 2-D rigid blocks interconnected by mass, damping and stiffness elements to simulate their suspension system. The passage of the vehicles is performed at different speeds. Several rolling surface profiles are admitted, considering the maintenance grade of the pavement. The spectral density functions are generated from an experimental database to form the longitudinal surface irregularity profiles. A computational code written in Phyton based on the finite element method was developed considering the Euler–Bernoulli beam model.

Several models of composite heavy vehicles are presented as manufactured and currently travel on major roads. Dynamic amplification factors are presented for each type of composite vehicle.

The VBI models for compound heavy vehicles are 2-D.

This work contributes to improving the safety and lifetime of the bridges, as well as the stability and comfort of the vehicles when passing over a bridge.

The structural response of the bridge is affected by the type and size of the compound vehicles, their speed and the conservative grade of the pavement. Moreover, one axle produces vibrations that can be superposed by the vibrations of the other axles. This effect can generate not usual dynamic responses.

Abnormal vibrations often occur in the liquid oxygen kerosene transmission pipelines of rocket engines, which seriously threaten their safety. Improper handling can result in failed rocket launches and significant economic losses. Therefore, this paper aims to examine vibrations in transmission pipelines.

In this study, a three-dimensional high-pressure pipeline model composed of corrugated pipes, multi-section bent pipes, and other auxiliary structures was established. The fluid–solid coupling method was used to analyse vibration characteristics of the pipeline under various external excitations. The simulation results were visualised using MATLAB, and their validity was verified via a thermal test.

In this study, the vibration mechanism of a complex high-pressure pipeline was examined via a visualisation method. The results showed that the low-frequency vibration of the pipe was caused by fluid self-excited pressure pulsation, whereas the vibration of the engine system caused a high-frequency vibration of the pipeline. The excitation of external pressure pulses did not significantly affect the vibrations of the pipelines. The visualisation results indicated that the severe vibration position of the pipeline thermal test is mainly concentrated between the inlet and outlet and between the two bellows.

The results of this study aid in understanding the causes of abnormal vibrations in rocket engine pipelines.

The causes of different vibration frequencies in the complex pipelines of rocket engines and the propagation characteristics of external vibration excitation were obtained.

This paper aims to investigate the effect of frictional heat on the wear of high-speed rotary lip seals in engines.

In this research paper, the authors focus on the high-speed rotating lip seal of aircraft engines. Using the hybrid lubrication theory, a thermal-fluid-solid coupled numerical simulation model is established to investigate the influence of parameters such as contact pressure distribution, temperature rise and leakage rate on the sealing performance under different operating conditions. By incorporating the Rhee wear theory and combining simulation results with experimental data, a method for predicting the wear of the rotating seal lip profile is proposed. Experimental validation is conducted using a high-speed rotating test rig.

The results indicate that as the speed increases, the rise in frictional heat leads to a decrease in the sealing performance of the lip seal contact region. The experimental results show a similar trend to the numerical simulation results, and considering the effect of frictional heat, the predicted wear of the lip seal profile aligns more closely with the actual wear curve. This highlights the importance of considering the influence of frictional heat in the analysis of rotating seal mechanisms.

This study provides a reference for the prediction of wear profiles of engine high-speed rotary lip seals.