Cross-sectional Shape Research Articles

Numerical analysis on heat transfer enhancement of Al2O3 and CuO-water nanofluids in annular curved tubes

The present investigation presents numerically the thermo-hydraulic performance of annular curved tubes in the turbulent flow region. The three-dimensional (3D) computational fluid dynamic (CFD) model was developed by using a commercial ANSYS 14.5 package to get additional insights on the thermo-hydraulic performance on a level of detail that is not always available in experiment results. The turbulent k-ε realize model was employed to examine the effect of Al2O3-water and CuO-water nanofluids on the heat transfer and fluid flow characteristics at different key design parameters. Three coil geometry shapes (including helical, spiral, and conical), five annular cross-section shapes (including circular, elliptical, square, rectangular, and triangular), and three cross-section areas were also investigated in this study. The numerical simulations were carried out at Reynolds numbers of 4700 to 26,700 and nanofluid volume concentrations, φ of 1%, 3%, and 5%, with constant wall temperature (CWT) heat transfer boundary conditions. The result showed that the addition of Al2O3 (45 nm) and CuO (29 nm) nanoparticles to water improves the heat transfer coefficient by 48.8%, 25.7%, and 8.4% and 37.1%, 20%, and 8.7%, respectively, at φ of 5%, 3%, and 1%, while the penalty of pressure drop is approximately negligible. The heat transfer rate per unit pumping power for helical design achieved 21.6% and 5.34% higher than conical and spiral designs, respectively, for the same curvature ratio, cross-section area, and coil length. Also, the circular shape achieved a higher heat transfer enhancement than elliptical, square, rectangular, and triangular shapes by 4%, 14.6%, 17.8%, and 30.1%, respectively. The maximum thermo-hydraulic performance index reached 1.68 and 1.58 for both Al2O3-water and CuO-water nanofluids, respectively, at φ of 5%.

Dumbbell-shaped thrombectomy device for cerebral venous sinus thrombus removal with controllable axial and longitudinal maneuverability

Abstract Cerebral venous sinus thrombosis (CVST) is frequently observed in younger adults and features in large thrombus volume. Due to the triangular-like cross-sectional shape and large diameter of superior sagittal sinus, all the commercially available artery stent retrievers are not suitable for venous vessels. In this study, a dumbbell-like stent was designed and fabricated by 3D braided technology using NiTi wires, which was manually rotatable and stretchable with controlled length/diameter ratios (2.6 to 14.0) and reciprocating maneuverability. Computational modeling and an in vitro study were conducted to evaluate the mechanical properties of this device and its ability to trap and remove thrombi from occluded venous vessels was verified using a swine model. And a single-center retrospective clinical study of 10 patients using the Venus-TD to treat patients with CVST was conducted. Pre/postoperative thrombus volume in ten patients was quantitively analyzed (12855.3 ± 6417.1 vs. 2373.1 ± 2759.0 mm³, P < 0.001) with a high recanalization rate, yielding favorable clinical outcomes. This study offers a novel treatment option for patients with extensive CVST.

Influence of Processing Parameters on Flow Behaviour of Ultra-Large-Section Beam Blank Continuous Casting Mould

The complex cross-sectional shape of oversized beam blanks and the size effect of ultra-large-section beam blanks create severe issues related to the surface and internal quality of the castings. To ensure quality and control in the production of ultra-large-section beam blanks, a numerical and physical model of molten steel flow in the three-port submerged entrance nozzle (SEN) mould, with section dimensions of 1300 × 510 × 140 mm, was established. This model was created using numerical simulations and NSGA-II genetic algorithm optimisation, and the impact of the casting speed and SEN immersion depth on the mould’s flow behaviour was investigated. The results showed that a deeper SEN immersion depth resulted in, a greater impact depth of the molten steel, and the surface flow velocity decreased. Both the impact depth and the surface flow velocity of the molten steel increased with increasing casting speed. The physical simulation and numerical simulation of the molten steel flow form and flow velocity distribution in the mould were in good agreement with each other, thus verifying the accuracy of the numerical simulation. The process parameters derived from this study were all within an appropriate range, which can help to improve the quality of continuously cast beam blanks. This also provides guidance for selecting optimal parameters for actual continuous casting production processes.

Molecular dynamics simulation of bending behavior of B2-FeAl alloy nanowires with different crystallographic orientations

In nanosystems, the metallic nanowires are subjected to significant and cyclic bending deformation upon integration into stretchable and flexible nanoelectronic devices. The reliability and service life of these nanodevices depend fundamentally on the bending mechanical properties of the metallic nanowires that serve as the critical components. A deep understanding of the deformation behavior of the metallic nanowires under bending is not only essential but also imperative for design and manufacture of high-performance nanodevices. To explore the mechanism underlying the bending plasticity of the metallic nanowire, we have conducted a study on the bending deformation of B2-FeAl alloy nanowires with various crystallographic orientations, sizes and cross-sectional shapes by using molecular dynamics simulation. Our results show that the bending behavior of the B2-FeAl alloy nanowires is independent of the size and cross-sectional shape of the nanowire, but it is highly sensitive to its axial orientation. Specifically, both <111>- and <110>-oriented nanowires yield by dislocation nucleation upon bending, in which the <111>-oriented nanowire fails by brittle fracture soon after yielding, while the <110>-oriented nanowire exhibits good ductility due to homogeneous plastic flow raised by continuous nucleation and steady motion of dislocations. In contrast to the aforementioned two nanowires, the bending plasticity of the <001>-oriented nanowire is mediated by stress-induced transformation from B2 to L1<sub>0</sub> phases, which leads to excellent ductility and higher fracture strain. The orientation dependence of bending deformation can be understood by considering the Schmid factor. Moreover, the plastically bent nanowires with <110> and <001> orientations are able to recover to their original shape upon unloading, particularly, the plastic deformation in the <001>-oriented nanowire is recoverable completely via reverse transformation from L1<sub>0</sub> to B2 structures, exhibiting superelasticity. This work elucidates the deformation mechanism of the B2-FeAl alloy nanowire subjected to bending load, which provides a crucial insight for the design and optimization of flexible and stretchable nanodevices based on metallic nanowires.

Study on the influence of section shape of underground roadway on ventilation resistance

ABSTRACT Among the parameters that affect the ventilation capacity of an underground mine, the cross-sectional shape of the roadway plays a significant role. When the required airflow is constant, the roadways with different shapes but the same area evidently overcome the different wind resistance. In order to investigate the influence of roadway cross-sectional shape on airflow resistance, this study carried out experimental and Simulation Research on the influence of different cross-section shape and size changes on ventilation resistance. And the Relative Shape Factor (RSF) is established to quantitatively evaluate the influence of cross-section shape on ventilation resistance. Firstly, focusing on the four common types of cross-sectional in the roadway, including circular, rectangular, trapezoidal, and arched, the experimental platform of airflow resistance of the pipe is built, and the effect of the changes of cross-sectional shape on the ventilation resistance is studied. Furthermore, the law of airflow resistance with different scale changes under the same section shape is studied by numerical simulation. Subsequently, the scale-independent coefficient RSF is established through theoretical derivation to characterize the degree to which the section shape deviates from the circle. The further it deviates from the circular shape, the larger the RSF, indicating a greater impact of the cross-section shape on airflow resistance. Comparing the theoretical with the simulation results, the RSF errors of the rectangular, trapezoidal, and arched shapes are 7.3%, 4.01%, and 2.05%, respectively, indicating that the RSF can accurately describe the influence of the cross-section on ventilation resistance. Finally, the range of RSF for common roadway dimensions is summarized, the trapezoidal is 1.131 to 1.187, the arch-shaped is 1.0 to 1.160, and the rectangular is 1.128 to 1.152. The proposed RSF in this paper is intuitive and effective, which provides a reference for the engineering application and a new insight for the selection of the minimum ventilation resistance cross-section of roadway.

Sensitivity analysis of geometrical parameters of supercritical water in twisted spiral tubes

ABSTRACT The high thermal efficiency of supercritical water makes it a promising alternative to water cooled reactors. This study employs numerical analysis, utilizing the SST k-ω turbulence model, to investigate the heat transfer performance of supercritical water (SCW) in various tube configurations and fluid flow conditions across Reynolds numbers ranging from 8,000 to 20,000. This research examines how different geometrical parameters, such as the helical direction, number of lobes, and cross-sectional shape, impact the flow physics and heat transfer performance of different spiral tubes. The outcomes specify that increasing the lobed number in the tube improves the heat flux by about 5%–16%. Furthermore, introducing two direction changes in the twisted tube will cause a slight increase (about 20%) on heat transfer enhancement. Finally, the sensitivity analysis of heat flux and Nusselt number to each of the effective parameters in the heat transfer of supercritical water in twisted tubes has been accomplished and the Response Surface Method (RSM) utilizes the central composite design (face-centered) approach. According to the findings of these two studies, it has been established that the Reynolds number of the fluid flow is the most influential parameter in determining the extent of heat transfer. Specifically, it exerts an effectiveness of 26% and 76% on the Nusselt number and heat flux, respectively. Furthermore, the inscribed circle diameter of tube by 8% effectivity and minor axis length by 5% effectivity on heat flux are more effective than others.

Numerical Study on the Shear Behavior of a Late-Model Cold-Formed Stainless Steel C-Shaped Beam.

The failure mode of thin-walled C-channel beams typically manifests as premature local buckling of the compression flange, leading to insufficient utilization of material strength in both the flange and the web. To address this issue, this study adopts the approach of increasing the number of bends to reinforce the flange and adding V-shaped stiffeners in the middle of the web to reduce the width-to-thickness ratio of the plate elements, thereby delaying local buckling and allowing for greater plastic deformation. However, the challenge lies in the irregular cross-sectional shape and complex buckling patterns. Therefore, this paper aims to explore a suitable cross-sectional form to expand the application of stainless steel members. Subsequently, three-point bending tests were conducted on the optimally designed stainless C-channel beam with folded flanges and mid-web stiffeners. The finite element simulation results were compared and analyzed with the experimental results to validate the model's effectiveness. After verifying the correctness of the finite element model, this study conducted numerical parameterization research to investigate the effects of the shear span ratio, complex edge stiffeners, web height-thickness ratio, and V-shaped stiffener size on the shear performance of stainless steel folded flange C-beams. The results show that changing the shear span ratio has a significant impact on the shear capacity and vertical deflection deformation of components; increasing the web height-thickness ratio can enhance the shear capacity of the component; elevating the V-shaped stiffener size can slightly improve the shear capacity of components; and for the stainless steel C-shaped beam with folded flanges and intermediate stiffening webs, adding edge stiffeners cannot remarkably promote the shear capacity of the component.

Dynamic behavior and stability analysis of fluid-conveying nano-spiral shells: Modeling, simulation, and parametric study

This paper presents a study on a fluid-conveying nano-spiral shell’s free vibration and stability analysis. The nano-spiral shell with Archimedean cross-section shape is modeled according to the first-order shear deformation (FSDT) and modified coupled stress (MCST) theories. Van der Waals force is modeled based on the Leonard Jones potential for atomic interactions. The fluid-structure interaction equation is derived based on the assumption of slip boundary conditions for the fluid. The Rayleigh-Ritz method is used to obtain the natural frequencies of a macro-size spiral shell to be validated with the finite element method (FEM). The parametric vibration analysis is performed for a nano-spiral shell conveying fluid. The results show that the natural frequencies decrease as the distance between the layers of the nano-spiral shell increases. Additionally, an increase in length results in a decrease in natural frequencies. Stability analysis is conducted, and flutter and divergence bifurcations are accurately studied. The results show that the shell with clamped boundary conditions is more stable than the supported boundary condition due to its higher bending stiffness. As the distance between the layers of the fluid-conveying nano-spiral shell increases, instability occurs at a lower fluid velocity. Conversely, as the thickness of the shell increases, the nanostructure becomes unstable at higher fluid velocities.

Research on Multimodal Control Method for Prosthetic Hands Based on Visuo-Tactile and Arm Motion Measurement.

The realization of hand function reengineering using a manipulator is a research hotspot in the field of robotics. In this paper, we propose a multimodal perception and control method for a robotic hand to assist the disabled. The movement of the human hand can be divided into two parts: the coordination of the posture of the fingers, and the coordination of the timing of grasping and releasing objects. Therefore, we first used a pinhole camera to construct a visual device suitable for finger mounting, and preclassified the shape of the object based on YOLOv8; then, a filtering process using multi-frame synthesized point cloud data from miniature 2D Lidar, and DBSCAN algorithm clustering objects and the DTW algorithm, was proposed to further identify the cross-sectional shape and size of the grasped part of the object and realize control of the robot's grasping gesture; finally, a multimodal perception and control method for prosthetic hands was proposed. To control the grasping attitude, a fusion algorithm based on information of upper limb motion state, hand position, and lesser toe haptics was proposed to realize control of the robotic grasping process with a human in the ring. The device designed in this paper does not contact the human skin, does not produce discomfort, and the completion rate of the grasping process experiment reached 91.63%, which indicates that the proposed control method has feasibility and applicability.

Investigation and optimization of parameters for bidirectional composite vibratory finishing of gears

ABSTRACT To meet the core requirements for uniform and efficient finishing of high-precision gears, a bidirectional composite vibratory finishing (BCVF) process was introduced. The influence of motion parameters (gear rotational speed n, vibration frequency f x and f z , vibration amplitude A x and A z , vibration phase difference ψ), process parameters (media loading height H, gear installation position h), gear parameters and container geometric parameters (cross-sectional shape, size L x × L y ) on the finishing effect was investigated by kinematic analysis and discrete element simulation (DEM). On this basis, the parameters with significant influence were selected and further optimized by gray relational analysis (GRA). Numerical simulations were carried out by orthogonal design, taking the cumulative contact energy and distribution uniformity on gear surface as evaluation indicators for multi-objective optimization. The results indicated that the affecting significance on comprehensive performance was ranked as f x > f z > A x > A z . Besides, the optimal parameter combination was obtained.

Shear Force in RC Internal Beam-Column Connections for a Beam Loaded With a Transverse Force Occupying Different Possible Positions

Frame structures are vulnerable to seismic impacts. The frame joint is the element that is often responsible for compromising frame structures. In recent years, shear force in the beam-column connection has been pointed out as the main culprit for damage in the nodes, and from there to the entire structure. The complex nature of the stressed and strained state of the joint is due to the poor knowledge of the forces passing in the beam-beam or column-column direction. For their more accurate determination, one departs from the standard acceptance of the static scheme of the axis line of the structural elements and works with their actual dimensions. With a skillful selection of the support devices, a mathematical model of the beam is created, allowing the determination of the forces that arise along its height and that enter the joint. The magnitudes of the support reactions of the beam are applicable both to operation in the elastic stage and to the plasticization of the beam and the occurrence of a crack between the beam and the column on the face of the column. In the present work, a cantilever beam loaded with a transverse force, occupying different possible positions on the beam, is considered. The expressions for the support reactions are derived. The unfavorable position of the loading force, which results in the greatest shear force, was investigated. A comparison is made of the corresponding shear force in the beam-to-column connection, with that recommended in the literature. The results demonstrate differences of up to 18.20%. The main parameters of the cantilever beam that were monitored are crosssectional shape, modulus of elasticity of concrete and position of the loading force.

MATHEMATICAL MODELING OF THE DYNAMICS OF THE AEROELASTIC "PIPELINE – PRESSURE SENSOR" SYSTEM

This paper proposes mathematical models of mechanical systems "pipeline - pressure sensor" designed to control the pressure of the working medium in the combustion chambers of engines. In such systems, to mitigate the impact of vibration accelerations and high temperatures, the sensor is connected to the engine using a pipeline and is located at some distance from it. The movement of the working medium is described by linear models of fluid and gas mechanics. To describe the dynamics of an elastic sensitive element, linear models of the mechanics of a solid deformable body are used. Based on linear differential equations with partial derivatives, mathematical formulations of problems are proposed that correspond to three-dimensional models of pressure measurement systems in gas-liquid media for some pipeline cross-sectional shapes, namely, for a pipeline with a rectangular cross-section, with a section in the form of a sector and in the form of a ring. By introducing integral characteristics, the solution of problems is reduced to studying one-dimensional models. Equations have been obtained that make it possible to determine the pressure of the working medium in the combustion chamber at each moment of time by the value of deformation of the sensitive element of the sensor. Analytical and numerical-analytical methods for solving the corresponding initial-boundary value problems for systems of differential equations are proposed. In the analytical approach, the solution of the problem is reduced to solving a differential equation with a deviating argument. The numerical-analytical study of the problem is based on the application of the Galerkin method. Also, a numerical experiment was carried out and examples of calculating the deformation of the sensitive element of the sensor in the case of rigid fastening when specifying specific values of the mechanical parameters of the system are presented, also during setting the law of change in the excess pressure of the working medium in the engine.

Investigation of heat transfer and fluid flow in annular curved tubes in laminar flow using Al2O3–water nanofluid

In the present study, the combination of both the secondary flow and the nanofluids as a passive approach as a method to maximize the effectiveness of heat exchangers, and reducing the size was examined. The thermo-fluid behavior of annular curved tubes was investigated numerically ANSYS-FLUENT 14.5 CFD package. The effects of key design parameters, such as different annular tube designs (helical, spiral, and conical), different annular cross-sectional shapes (circular, elliptical, square, rectangular, and triangular), and different annular cross-sectional areas, were studied to evaluate the best design for annular curved heat exchanger. Furthermore, the effects of Al2O3–water nanofluids were presented to evaluate their superiority over the water. To achieve thermal and fluid characteristics, the coils are the same length and cross-sectional area. The numerical simulations were performed at Reynolds numbers, Re, of 850–6400 for various nanofluid concentrations, φ, of 1%, 3%, and 5%. The numerical findings showed that the annular helical tube design exhibited better performance than both the spiral and conical tube designs by 5.1% and 10.4%, respectively, while the pressure drop, ΔP, growth was almost negligible. A circular cross section achieved better heat transfer enhancement among the elliptical, square, rectangular, and triangular shapes by 13%, 30%, 40.5%, and 52.9%, respectively, for the same cross section area. At an Al2O3/water nanofluid concentration of 5%, the heat transfer augmentation was 31.7% greater than that of water, while the pressure drop was 2.37 times greater than that of water. The maximum thermo-hydraulic performance index, η, reaches 1.1% for Al2O3–water at a nanofluid concentration of 5%.Graphical abstract

Shell Ontogeny of Otoceras concavum Tozer (Ceratitida) and the Problem of the Origin of Otoceratids

As a result of a detailed study of cross sections of the most ancient (Late Permian) representatives of Otoceras Griesbach (O. concavum Tozer) from the base of the Nekuchan Formation in the upper reaches of the Vostochnaya Khandyga River in Southern Verkhoyanie, the features of ontogenetic shape changes were clarified. Trends in changes in the main parameters of the shell during its growth have been determined, the cyclical nature of the narrowing and expansion of the shell has been identified, and the sequence of changes in the shape of the whorl cross-section has been established. The formation of carinae and the emergence of a tricarinate form on the ventral side were traced. Identification of shell shape transformations in the ontogenesis of O. concavum contributes to the diagnosis of small-sized Otoceras and can serve as the basis for subsequent reconstruction of the morphogenetic development of Otoceratidae. The pentacarinate form of the outer whorl, established at the middle stage of ontogenesis, may be a character inherited from the ancestor, the genus Avushoceras Ruzhencev.

Aerodynamic Performance Evaluation of Subsonic Compressor Cascade Blade With Leading-Edge Damage

Abstract Material loss at the leading edge of in-service gas-turbine engine blades is detrimental to the aerodynamic performances of compressors. In this study, computational simulations and experiments using pneumatic probes were conducted to measure the total pressure loss of compressor cascade blades with material loss at their leading edges. The correlations between blade leading-edge damage geometries and cascade blade total pressure losses were analyzed. The geometrical features of blade leading-edge damage that most affect the cascade blade total pressure loss were distinguished, and the influences of leading-edge damage region cross-sectional shapes on the cascade blade total pressure loss were examined. Two methods for evaluating the total pressure loss of damaged cascade blades were developed using experimental and numerical data. In the first method, an attempt is made to correlate the cascade blade total pressure loss with a characteristic parameter that is defined using scale and angle parameters of the leading-edge damage region. The second method is based on the probability distribution function as well as the critical loss boundary of the geometrical parameters of the leading-edge damage region. The results of this study provide a better understanding of the influence of damaged leading edges on compressor cascade blade aerodynamic performance and can help establish low-cost methods for the estimation of compressor performance degradation.

Numerical Analysis of the Trapping Effect of Grooves with Various Cross-Sectional Shapes and Reynolds Numbers

Numerical Analysis of the Trapping Effect of Grooves with Various Cross-Sectional Shapes and Reynolds Numbers

A novel two-layer composite geomembrane lining structure to mitigate frost damage in cold-region canals: Model test and numerical simulation

The canal is crucial for water diversion projects, but it is susceptible to frost damage. To address this, a two-layer composite geomembrane lining structure (TLCGLS) was proposed that regulates the interaction between canal lining and frozen soil. Model tests were conducted to investigate its anti-frost heave effectiveness. Considering the interaction among the lining, two-layer composite geomembranes (TLCGs), and frozen soil, a canal frost heave model with heat-water-mechanical coupling was developed. The influence of canal cross-section shapes and TLCGs arrangements on anti-frost heave performance and mechanism of TLCGLS were discussed. Results show that TLCGLS reduces uneven frost heave degree and compressive/tensile strains of the lining by 35%, 29%, and 28% respectively. During melting, it rapidly reduces frost heave, tangential deformation, and strain with minimal residual effects. TLCGLS demonstrates strong resetting ability and excellent anti-frost heave performance. It is particular suitable for arc-bottomed trapezoidal canals. However, excessive reduction in friction between TLCGs weakens arching effect of the bottom lining, increasing tensile stress and safety risks. TLCGLS with geomembrane-geotextile contact exhibits superior anti-frost heave performance, mitigating compressive stress by over 50% while meeting design requirements for tensile stress. These findings provide a theoretical basis and technical solution for mitigating frost damage in canals.

An experimental study on flow induced motion and energy harvesting of cylinders with different cross sections

It has been known that the cross-sectional shape of a column oscillator significantly influences its vibrational characteristics and energy conversion capacity, and can alter the nature of flow-excited vibration (FIV). Whether the addition of appendages to oscillators with different cross-sectional shapes enhances energy conversion capacity remains uncertain. In this study, the vibration characteristics and energy capture capabilities of an elastically supported oscillator with a semicircular appendage, suitable for low-speed seafloor current environments, are investigated. Experiments were conducted at zero degrees incidence for Reynolds numbers ranging from 5.041 × 10³ to 7.562 × 104, resulting in turbulent wake conditions. The hydrodynamic properties of the oscillators were evaluated through statistical analysis, Proper Orthogonal Decomposition (POD), and vortex core identification of Particle Image Velocimetry (PIV) fields. The energy capture capability of the oscillator was assessed through statistical analysis of its vibration displacement, frequency, and amplitude. The study's results indicate that an oscillator with symmetric sharp attachments and without vortex reattachment is favorable for galloping with self-excitation. Under equivalent conditions, the Circular-T-shaped oscillator demonstrates superior energy conversion capacity compared to existing models, with the galloping branch being the most efficient for energy conversion; the peak efficiency is 24.5% (Ur = 14.5). This study provides some baseline data and optimization solutions for flow-induced motion power generation.

Effects of undulating printing paths on axial compressive behaviors of 3D-printed continuous fiber-reinforced multi-cell thin-walled structure

Thin-walled structures are widely used in passive safety systems of the vehicles because of their excellent lightweight and stable progressive deformation behavior. Recently, continuous fiber 3D printing technique has been successfully applied in the manufacturing of thin-walled structures, providing a promising path to meet increasing requirements for energy absorbers as vehicle speed increases. This paper proposed a range of novel undulating printing paths to fabricate two kinds of continuous Kevlar fiber reinforced multi-cell thin-walled structures (CFMSs) with different cross-section shapes, named MT1 and MT2, respectively. For each type of CFMS, three different printing paths (moving back and forth between two layers) were designed to manufacture CFMS with the wall thickness comprised of only a single deposited line. The axial compressive tests were conducted to investigate the compressive response of CFMSs. It was found that undulating printing paths affect both the mechanical behavior and collapse mode of CFMSs. During the compression process, MT1 and MT2 printed with path 3 showed progressively folding deformation modes (differed from the progressive crushing mode of composite laminated thin-walled structures). Thus, compared to MT2 printed with path 1 and path 2, which exhibit global buckling, MT1 and MT2 printed with path 3 exhibited more stable energy absorption performance. Besides, although the local failure mode in different CFMSs caused by the same path was similar, the same printing path had different influences on the axial compressive properties of different CFMSs.

Simulation of 3D concrete printing using an implicit formulation of the moving particle semi-implicit method

Additive construction techniques, such as 3D Concrete Printing (3DCP), offer significant advantages, as reduced material waste, optimized construction times, and the possibility to create complex shapes that would be difficult to achieve through conventional methods. Despite the advances in 3DCP, some challenges persist due to the still incipient development to predict the printing results. The fresh properties of the material play a crucial role. Therefore, in addition to experimental studies, numerical modeling emerges as a valuable tool to investigate the dynamic aspects of the extrusion process in 3DCP. However, the numerical modeling of 3DCP is challenging due to the complex rheological behavior of the self-supporting concrete, which typically exhibits high values of apparent viscosity at low deformation rates. To address these challenges, we adopted the moving particle semi-implicit (MPS) method, which is a particle-based simulation technique suitable for modeling free surface flows and materials undergoing large deformations. As a Lagrangian approach, MPS discretizes the computational domain into particles that move according to the governing equations of fluid motion. The conventional formulation of the MPS presents numerical restrictions when dealing with highly viscous fluids. To assure numerical stability, the semi-implicit approach requires very low time steps, resulting in high processing costs. In this way, this study adopts an implicit algorithm that allows for significantly larger time steps, resulting in reduced processing time compared to the conventional semi-implicit formulation. To validate the implicit algorithm, simulations of non-Newtonian flow between parallel plates were validated using analytical solutions, demonstrating better accuracy at lower Bingham numbers. The 3DCP simulations were conducted to describe the flow of fresh mortar using the Bingham-Papanastasiou constitutive model. The geometry of the 3D model includes the extrusion nozzle and a planar surface where concrete is printed. The numerical investigation involved simulations with different nozzle heights (Hn), printing speeds (V), and extrusion volumetric fluxes (U). The cross-sectional shapes of the extruded layers were compared with experimental data, showing qualitative agreement. The results demonstrate the potential of the implicit MPS method for modeling 3DCP processes.