3D Numerical Analysis Research Articles

The Combination of X-ray Ct and Multi-particle Finite Element Method for Powder Compaction

Abstract The macroscopic mechanical properties of granular materials are closely related to their microscopic properties, and studying the compression problem of powder is helpful for establishing a constitutive law of powder compaction process. This study utilized X-ray CT technology to investigate the 3D numerical analysis model of loose particle stacking and compaction based on MPFEM, and investigated the mechanical behavior of powder under compression. The results of MPFEM were experimentally verified through powder uniaxial compression experiments, and the results showed that the calculated results of MPFEM were in good agreement with the experimental results. The method of X-ray CT can effectively capture the geometric characteristics of powder particles, and the MPFEM can obtain the true characteristics of plastic deformation of powder bodies, which is an important means of developing powder compaction constitutive models.

Thermal Hydraulic Performance and Characteristics of a Microchannel Heat Exchanger: Experimental and Numerical Investigations

Abstract This paper presents extensive fluid flow and Heat Transfer studies conducted through an experimental setup followed by a detailed three-dimensional (3D) numerical analysis of the same setup using a commercial package for computational fluid dynamics (CFD), known as cfd-ace® for additive-manufactured counterflow AlSi10 Mg microchannel heat exchangers (MCHEs). A detailed 3D computational model of the experimentally tested MCHEs was built and analyzed using the commercial software cfd-ace® for the same experimentally tested operating conditions. The computational model results are in good agreement with experimental data of tested MCHE within +2% to +7% and ∼0% to −13.5% variation for cold and hot fluids for the entire set of design of experiments (DoEs). This percentage disagreement may be due to various factors, such as manufacturing deviation within tolerance, longitudinal conduction, variation in the thermal conductivity of the material after heat treatment, variation in environmental temperature, sensor deviation, and surface roughness of internal channels. Instead of Stainless steel (SST), AlSi10 Mg was used because of its lower manufacturing cost because AlSi10 Mg was lighter than SST, though its thermal conductivity is almost ∼8–10 times more than that of SST. A higher thermal conductivity is not good for MCHEs because it leads to higher longitudinal conduction, which eventually degrades the performance of MCHEs in terms of effectiveness. MCHE effectiveness is also reduced by ∼12% to 18% owing to longitudinal conduction from ideal effectiveness.

Three-Dimensional CFD Analysis of a Hot Water Storage Tank with Various Inlet/Outlet Configurations

This study presents a comprehensive 3D numerical analysis of thermal stratification, fluid dynamics, and heat transfer efficiency across six hot water storage tank configurations, identified as Tank-1 through Tank-6. The objective is to determine the most effective design for achieving uniform temperature distribution, stable stratification, and efficient heat retention in sensible heat storage systems, with potential for integration with phase change materials (PCMs). Using COMSOL Multiphysics 5.6, simulations were conducted to evaluate key performance indicators, including the Richardson number, capacity ratio, and exergy efficiency. Among the tanks, Tank-1 demonstrated the highest efficiency, with a capacity ratio of 84.6% and an exergy efficiency of 72.5%, while Tank-3, which achieved a capacity ratio of 70.2% and exergy efficiency of 50.5%, was identified as the most practical for real-world applications due to its balanced heat distribution and feasibility for PCM integration. Calculated dimensionless numbers (Reynolds number: 635, Prandtl number: 4.5, and Peclet number: 2858) indicated laminar flow and dominant convective heat transfer across all the configurations. These findings provide valuable insights into the design of efficient thermal storage systems, with Tank-3’s configuration offering a practical balance of thermal performance and operational feasibility. Future work will explore the inclusion of PCM containers within Tank-3, as well as applications for heat pump and solar water heaters, and high-temperature heat storage with various working fluids.

Three-Dimensional Numerical Analysis of Seismic Response of Steel Frame–Core Wall Structure with Basement Considering Soil–Structure Interaction Effects

In recent years, numerous studies highlighted the crucial role of the soil–structure interaction (SSI) in the seismic performance of basement structures. However, there remains a limited understanding of how this interaction affects buildings with basement structures under varying site conditions. Based on the three-dimensional (3D) numerical analysis method, the influence of the SSI on the seismic response of high-rise steel frame–core wall (SFCW) structures situated on shallow-box foundations were investigated in this study. To further investigate the effects of the SSI and site conditions, three types of soil profiles—soft, medium, and hard—were considered, along with a fixed-foundation model. The results were compared in terms of the maximum lateral displacement, inter-story drift ratio (IDR), acceleration amplification coefficient, and tensile damage for the SFCW structure under different site conditions, with both fixed-base and shallow-box foundation configurations. The findings highlight that the site conditions significantly affected the seismic performance of the SFCW structure, particularly in the soft soil, which increased the lateral deflection and inter-story drift. Moreover, compared with non-pulse-like ground motion, pulse-like ground motion resulted in a higher acceleration amplification coefficient and greater structural response in the SFCW structure. The RC core wall–basement slab junction was a critical region of stress concentration that exhibited a high sensitivity to the site conditions. Additionally, the maximum IDRs showed a more significant variation at incidence angles between 20 and 30 degrees, with a more pronounced effect at a seismic input intensity of 0.3 g than at 0.2 g.

3D Numerical Analysis of a Geogrid-Reinforced Piled Embankment: High-Speed Railway

3D Numerical Analysis of a Geogrid-Reinforced Piled Embankment: High-Speed Railway

Do gravity data justify a rifted “Liguro-Provençal Basin”?

The geodynamic evolution of the Liguro-Provençal Basin and its crust and upper mantle structure remain debated, especially regarding the role of rifting in continental break-up and seafloor spreading. Our study incorporates updated datasets, including new gravity maps from the AlpArray Gravity Working Group (complete Bouguer, free air, and isostatic anomalies) for 3D modeling and gravity field analysis, seismic data from Lobster offshore campaigns for direct comparison, and geodynamic models, supplemented by seismic profiles from previous French and Italian campaigns to constrain the interpretation. We used GFZ’s IGMAS + software for interactive 3D modeling, creating a density model extending to 300 km depth that includes crustal and upper mantle inhomogeneities based on prior geodynamic models. This hybrid approach, with polygonal structures for the crust and voxels for the upper mantle, clarifies individual contributions to the gravity field. Extending initial gravity modeling from the SPP MB4D project INTEGRATE, our work provides a consistent 3D density model for the Alps and Ligurian Basin. The constrained 3D modeling and numerical analyses (terracing, clustering, filtering, curvature), along with vertical stress and gravitational potential energy calculations, suggest that rifting has significantly influenced the basin’s geological evolution.

Stability Assessment of the Maltravieso Cave (Caceres, Spain) Through Engineering Rock Mass Classification, Empirical, Numerical and Remote Techniques

Caves have long fascinated humanity, serving as shelters, canvases for artistic expression and now significant attractions in the realm of tourism. Among these remarkable geological formations, the Maltravieso cave in Extremadura, Spain, stands out for its rich archaeological and paleontological heritage, particularly its collection of Paleolithic rock art. Despite its cultural significance, there is a notable dearth of studies addressing the stability of the cave from an engineering perspective. This article presents a pioneering study aimed at assessing the stability of the Maltravieso cave through a multidisciplinary approach: using empirical geomechanical classifications such as the Q Index, Rock Mass Rating (RMR) and the recently formulated Cave Geomechanical Index (CGI), alongside other techniques like Structure from Motion (SfM), 2D numerical modeling and 3D wedge analysis. This research aims to fill the gap in our opinion of cave stability assessment. By combining field data collection with sophisticated analysis methods, this study seeks to provide valuable insights into the geomechanical properties of the Maltravieso cave and validate a simple yet effective methodology for evaluating the stability of natural caves. This work not only contributes to the body of knowledge regarding cave geomechanics but also underscores the importance of preserving these invaluable cultural and geological treasures for future generations.

Experimental and 3D numerical analysis of embankment on soft soil improved with cement bottom ash columns

Experimental and 3D numerical analysis of embankment on soft soil improved with cement bottom ash columns

Mechanical Characteristics of Deep Excavation Support Structure with Asymmetric Load on Ground Surface

The purpose of this paper is to capture the mechanical response of the support structure of deep excavation subject asymmetric load. A two-dimensional (2D) numerical analysis model was established by taking a pipe gallery deep excavation subject to asymmetric load as an example. The numerical analysis results were in good agreement with the measured data, thus verified the validity of the numerical model. On this basis, the stress and displacement of support structure caused by the change in foundation asymmetric load were studied. According to the numerical results, horizontal displacement of the diaphragm wall (DW) was dominant, and the maximum horizontal displacement of the DW was 7.54 mm when the deep excavation was completed. With the increase in asymmetric load, the left wall displacement continued to increase, while the displacement of the right DW continued to decrease, and the maximum horizontal wall displacement occurred near the excavation face. The DW was the main bending component, and the maximum wall bending moment when the deep excavation was completed was 173.5 kN·m. The maximum wall bending moment increased with the increase in asymmetric load, and the maximum wall bending moment on the left of the deep excavation was greater than that on the right. The inner support sustained the main component of axial force, with the axial force peaking at 1051.8 kN when the deep excavation was completed. The axial force of the inner support increased with increasing the asymmetric load, and the axial force of the second inner support was obviously greater than that of the first inner support. This research has a positive effect on the design and optimization of deep excavation support structure subject to asymmetric load on ground surface.

Comprehensive study of vortices interaction and blades height effect in a Darrieus vertical axis wind turbine with J‐type blades

AbstractThere is a growing demand to improve the performance of vertical axis wind turbines to facilitate their commercialization for application in urban areas. This study utilizes a 3D numerical analysis to examine the influence of different vortices generated on turbine efficiency with straight and J‐type blades. The numerical simulation of this study employs the Reynolds‐Averaged Navier–Stokes equations and ‎sliding ‎mesh techniques ‎to more accurately model the rotational motion of blades about the turbine axis in relation to the ‎wind. Comparing the output ‎torque and the flow field at different span‐wise sections, the J‐type blades achieve better ‎performance at mid‐spans where the effect of stall vortices is dominant. Conversely, the lower ‎performance of J‐type blades is seen at tip spans due to stronger tip vortices. Investigations also ‎reveal the criticality of the downwind region on the overall performance at high tip speed ratios. It is observed that by ‎increasing the height, the tip vortices are limited to the tip sections, and stall vortices expand further ‎along the blade. At TSR = 1, the improvement by J‐type blades rises from 10% at a height of 0.8 m to 44% ‎at 3 m. The growth in height at lower wind speeds becomes more beneficial. Compared to the straight blades, the self‐starting ‎generated torque by J‐type blades for heights of 0.8, 1.2, and 1.6 m, are improved by 15.6%, ‎‎26.9%, and 34.7%, respectively. Overall, it is concluded that by increasing the blade height, the superiority ‎of the J‐type blade becomes more noticeable as the blade mainly contributes to suppressing the stall ‎vortices effect where the tip vortices effect is not presented. ‎

Impact of trapezoidal ribs on the performance of solar air collector: A numerical solution with optimized rib dimensions for better performance

In the present study, a 2D numerical analysis of the solar air collector (SAC) of an indirect solar dryer having trapezoidal corrugations on the absorber plate was performed. Corrugation pitch, p (twelve values ranged from 20 to 160 mm) and height, e (six values ranged from 1 to 10 mm) were varied and analyzed for six values of Reynolds numbers (Re). The output characteristics such as Nusselt number (Nu), friction factor (f) and thermo-hydraulic performance index (Thp) were calculated for different p, e and Re. The total work was categorized into two parts (part-I for optimizing p and part-II for optimizing e). 18 domains (twelve for part-I and six for part-II simulations) were generated and 108 simulations were executed to find the optimum dimensions (p, e and corrugation angle, α) of the corrugation. ANSYS Fluent-v15 was used to solve the problem. The maximum Nu for the corrugated sheet was 2.663 times greater than the flat absorber plate. The maximum Thp range was from 1.435 to 1.699 and obtained at the optimal values of p = 140 mm, e = 4 mm and α = 38.66° The numerical results were compared with the existing literature.

Parametric analysis of innovative corrugated profile of soil‐steel composite bridge

AbstractComposite soil‐steel corrugated structures are well recognized and widely used for construction of culverts and pedestrian or animal crossings. Moreover, for the past decade corrugated soil‐steel structures because of their large span are also recognized and chosen for bridges and tunnels engineering. Record holder composite corrugated structure of 32.40 m span encourages new research in the development of innovations. This paper investigates innovative the deepest corrugation cross‐section strengthened with circular hollow section steel pipes influence on plate utilization of large span soil‐steel composite bridge. The numerical 2D model of 17.5 m span structure was developed for investigation. Parametric analysis indicated that introduction of circular pipes will reduce corrugated profile steel thickness because of reduction of buckling length of straight region of corrugation. Moreover, 3D numerical analysis of buckling shapes allowed to conclude that local buckling of the deepest corrugation of 500 × 237 mm could be controlled strengthening the plates with circular steel pipes.

3D modeling of deep borehole electromagnetic measurements with energized casing source for fracture mapping at the Utah Frontier Observatory for Research in Geothermal Energy

AbstractWe present a 3D numerical modelling analysis evaluating the deployment of a borehole electromagnetic measurement tool to detect and image a stimulated zone at the Utah Frontier Observatory for Research in Geothermal Energy geothermal site. As the depth to the geothermal reservoir is several kilometres and the size of the stimulated zone is limited to several 100 m, surface‐based controlled‐source electromagnetic measurements lack the sensitivity for detecting changes in electrical resistivity caused by the stimulation. To overcome the limitation, the study evaluates the feasibility of using a three‐component borehole magnetic receiver system at the Frontier Observatory for Research in Geothermal Energy site. To provide sufficient currents inside and around the enhanced geothermal reservoir, we use an injection well as an energized casing source. To efficiently simulate energizing the injection well in a realistic 3D resistivity model, we introduce a novel modelling workflow that leverages the strengths of both 3D cylindrical‐mesh‐based electromagnetic modelling code and 3D tetrahedral‐mesh‐based electromagnetic modelling code. The former is particularly well‐suited for modelling hollow cylindrical objects like casings, whereas the latter excels at representing more complex 3D geological structures. In this workflow, our initial step involves computing current densities along a vertical steel‐cased well using a 3D cylindrical electromagnetic modelling code. Subsequently, we distribute a series of equivalent current sources along the well's trajectory within a complex 3D resistivity model. We then discretize this model using a tetrahedral mesh and simulate the borehole electromagnetic responses excited by the casing source using a 3D finite‐element electromagnetic code. This multi‐step approach enables us to simulate 3D casing source electromagnetic responses within a complex 3D resistivity model, without the need for explicit discretization of the well using an excessive number of fine cells. We discuss the applicability and limitations of this proposed workflow within an electromagnetic modelling scenario where an energized well is deviated, such as at the Frontier Observatory for Research in Geothermal Energy site. Using the workflow, we demonstrate that the combined use of the energized casing source and the borehole electromagnetic receiver system offer measurable magnetic field amplitudes and sensitivity to the deep localized stimulated zone. The measurements can also distinguish between parallel‐fracture anisotropic reservoirs and isotropic cases, providing valuable insights into the fracture system of the stimulated zone. Besides the magnetic field measurements, vertical electric field measurements in the open well sections are also highly sensitive to the stimulated zone and can be used as additional data for detecting and imaging the target. We can also acquire additional multiple‐source data by grounding the surface electrode at various locations and repeating borehole electromagnetic measurements. This approach can increase the number of monitoring data by several factors, providing a more comprehensive dataset for analysing the deep‐localized stimulated zone. The numerical analysis indicates that it is feasible to use the combination of the energized casing and downhole electromagnetic measurements in monitoring localized stimulated zone at large depths.

Definition and validation of a valley amplification factor for seismic linear response of 2D homogeneous alluvial basins

The paper presents findings from a parametric study analysing geometric (e.g., shape ratio, edge inclination) and stratigraphic factors (e.g. impedance ratio) influencing ground motion in trapezoidal valleys. The study involved 2160 visco-elastic analyses, considering 180 2D models with diverse shapes and soil properties, undergoing 12 synthetic input motions. Analyses results showed that the motion at the valley centre increases with both shape and impedance ratios, while it is independent of the edge slope; on the other hand, the maximum amplification at the edges depends on their inclination and on the impedance ratio, while it is independent of the valley shape. The position and size of the zone of maximum amplification at the edges depend on all the previous parameters. A valley amplification factor (VAF) is introduced to quantify spectral acceleration increase due to 2D effects. Closed-form equations are proposed to evaluate VAF based on valley properties. The proposed VAF is then applied to predict seismic amplification in two central Italian valleys, providing results well-comparable to those obtained from 2D numerical analyses. The described approach can be easily implemented into codes of practice as a conservative design tool to estimate 2D amplification along the surface of ‘shallow valleys’ subjected to moderate seismic actions.

3D transient numerical analysis of PCM dendritic cylindrical heat sinks designed by constructal theory

The phase change material (PCM) heat sink represents an important cooling tool for many applications. The presence of a heat sink in devices requires additional space to occupy it. This led to modeling this space to obtain the highest thermal performance with the smallest size. Among the modern modeling methods is constructal theory. This work presents a numerical modeling of the performance of cylindrical phase change material heat sinks with branched heat transfer fluid tubes in the manner of arterial and venous according to constructal theory. The aim of numerical simulations is to demonstrate the effect of growth in design by increasing the number of branches on thermal performance and melting time and the extent to which this is reflected in pressure drops. Several branching configurations are studied for comparison purpose. The governing continuity, momentum and energy equations were solved using the computational fluid dynamics method. Three-dimensional transient simulations of the melting process of phase change materials were performed using an enthalpy and porosity model. Four combinations of constructal dendritic tubes were considered. The results show that the proposed dendritic PCM heat sink significantly reduces the melting time of PCM compared to the traditional heat sink consisting of a shell and single tube up to 69.1 %. The drop in heat transfer fluid temperature at the outlet decreased by up to 51.7 %. The results proved also the heat transfer enhancement is governed with pressure drop increase because of changing branches configuration. Finally, this paper demonstrated the ability of constructal theory to provide promising solutions for cylindrical heat sinks of phase change materials in various applications.

Multidisciplinary design and manufacturing of a Tesla pump prototype

To widen the range of hydraulic efficiencies of boundary layer pumps, a full design methodology has been proposed in order to identify critical issues for their performance and manufacturing. The methodology integrated a 2D numerical code, CFD and FEM analyses, coupled with manufacturing assessments as feedback mechanism. Considering budget constraints and in-house machining capabilities, a quick first prototype was produced. Analyses of the design are pointing out that the volute design initially chosen will not help to achieve an increase in the overall efficiency. The curves of head achieved with 2D and CFD are in agreement, but the latter determines the losses with larger accuracy, thus achieving lower values of head. The 2D model shows limits in the determination of the efficiency, effectively corrected by the CFD analysis. Critical parameters as disc thickness and gap between discs will require a more sophisticated assembly process and materials outsource. The proposed methodology could be used as a reference for the design and performance evaluation of this kind of turbomachinery in the future. The procedure lead to a prototype design, whose optimal efficiency slightly lower than 30 % was achieved at 5000 rpm with 0.3 mm disks gap.

Proposed formulas for pillar stress estimation in a regular room-and-pillar pattern

Room-and-pillar mining is one of the oldest and most widely employed underground mining methods worldwide, while in recent years it has also been used for civil engineering projects and applications. The pillar design process is quite straightforward, requiring the assessment of two key elements, namely the pillar's strength and the anticipated loading imposed on it. Lots of research work has been elaborated for the estimation of pillars' strength and the main parameters associated with it, while, on the other hand, in terms of the loading the Tributary Area Theory (TAT) has been traditionally and successfully used as the dominant methodology to assess the post-mining pillar induced stresses. This paper focuses primarily on decoding the pillar loading regime by analyzing the identifying the stresses imposed to them through the use of modern computational tools offered by 2D and 3D numerical analyses, over a series of pillars' type (square, ribs), configurations and initial virgin stress fields. This allows for the development of a benchmarking framework that showcase the differences in the stress conditions and reveal the over-conservative behavior of TAT. More importantly though, it allowed for the expression of two (2) general analytical formulae for the direct estimation of the average vertical elastic stresses on pillars, both for the case of rib & square pillars layout; the most typical patterns applied in this method. They aim at providing an easy to use and more accurate estimation of the pillar stress - as compared to TAT - in alignment with the results obtained from numerical analyses. Thus, with the proposed formulas, the time-consuming numerical process needed for preliminary pillars' dimensioning is tackled in a swift and effective manner.

2D Numerical Analysis of Prefabricated Vertical Drains Using Different Matching Methods

2D Numerical Analysis of Prefabricated Vertical Drains Using Different Matching Methods

Performance improvement of novel latent thermal energy storage system with two-layer metal foam porosities

A 2D numerical analysis was conducted to examine the melting performance of a shell-and-tube-type phase change material with a two-layer copper foam having higher and lower porosity levels. The design parameters for these two layers of metal foam consisted of three different volume fill ratios (VFR) (30 %/70 %, 50 %/50 %, and 70 %/30 %) and three different shapes (circular, semi-circular, and elliptical) for the lower porosity layer. The results showed that a semicircular-shaped low porosity layer with a 30 % VFR had the least melting time of 475 s. When the VFR increased to 50 %, the elliptical and circular shapes had the least melting time of 394 s. However, at a VFR of 70 %, only the elliptical shape had the least melting time of 323 s. Furthermore, the results were compared with the uniform porosity metal foam system, and deduced that the two-layer system did not always have the lower melting time compared with that of the uniform porosity system because it depended on the VFR and shape of the different porosities filled in the different layers. Additionally, for an average porosity of the metal foam, an optimal combination of two different porosities was observed for the melting performance.

Study on significance of pillar dimensions on stability of caverns in rocks with high overburden

The present study aims to examine the influence of pillar widths on the stability of caverns. Case study considering two caverns viz. Powerhouse Cavern (PHC) and Transformerhall Cavern (TC), in a major hydro-electric project in the eastern Himalayas is considered. 2D and 3D numerical analysis was carried out for w/B ratios 1.5 and 3.0 respectively. Primary aspects like major principal stress and development of plastic zone were investigated for the two pillar widths. An optimum pillar width was observed that resulted in reduced stress acting along the cavern periphery, a better stress distribution, and no overlap of plastic zones between the caverns. Further, the optimum pillar width resulted in a better stress-redistribution with progress of excavation and the in-situ stress became constant at an earlier stage of excavation. Observations from comparative analysis revealed that a pillar width nearly equal to the largest dimension or twice the width of the larger of the caverns in the group resulted in a better stability and hence can be considered as the optimum width. Furthermore, the analysis suggests that along the pillar width, maximum stress was observed at mid-height, and it is more in the vicinity of the face of the caverns.