Alloy composition design usually contributes to eliminating cracking in Ni-based superalloys during addi-tive manufacturing(AM).However,a detailed understanding of each solid solution element in the crack-ing suscepti...Alloy composition design usually contributes to eliminating cracking in Ni-based superalloys during addi-tive manufacturing(AM).However,a detailed understanding of each solid solution element in the crack-ing susceptibility of Ni-based superalloys during AM is still missing.Thirteen newly designed alloys are considered to investigate the combined effect of solid solution elements on cracking susceptibility.The behaviors of solidification cracking,liquation cracking,and solid-state cracking were analyzed by the microstructural characterization and thermodynamic calculations.The results showed that W and Mo cause the formation of high melting-point carbides at grain boundaries(GBs),which increase solidifica-tion cracking susceptibility.Moreover,W and Mo lead to a slightly higher solidification cracking index(SCI)compared to Co,Cr,and Re.In the successive solidification and remelting process,the borides en-riched in W,Mo,and B around GBs will cause grain boundary segregation and liquation cracking.W and Re extend the freezing range(FR)and exacerbate the segregation of Al and Ti in the inter-dendritic re-gions,contributing to the formation of eutectics.Similarly,complete or partial melting of the eutectic can induce liquation cracking during the thermal cycling in AM.The solid-state cracking susceptibility can be reduced by solid solution elements,especially Re and Co.In summary,compared to Co,Cr,and Re,W and Mo exacerbate the cracking susceptibility.展开更多
In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. ...In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.展开更多
Truck frames should be designed and fabricated with enough rigidity to avoid excessive deflections. Finite element analysis (FEA) plays an important role in all stages of frame designs. While being accurate, 3D solid ...Truck frames should be designed and fabricated with enough rigidity to avoid excessive deflections. Finite element analysis (FEA) plays an important role in all stages of frame designs. While being accurate, 3D solid element FEA models are built upon frame configuration details which are not feasible in the preliminary design stage, partially because of limited available design data of frames and heavy computation costs. This research develops 1D beam element FEA models for simulating frame structures. In this paper, the CAD model of a truck frame is first created. The solid element FEA analysis, which is adopted as the baseline in this study, is subsequently conducted for the stiffness of the frame, Next, beam element FEA analysis is performed for validating the feasibility of the beam element FEA model by comparing the results from the solid and beam element FEA models. It is found that the beam element FEA model can predict the frame stiffness with acceptable accuracy and reduce the computation cost significantly.展开更多
Significant springback occurs after tube rotary-draw-bending （RDB）, especially for a high-strength Ti-3A1-2.5V tube （HSTT） due to its high ratio of yield strength to Young＇s modulus. The combination scheme of exp...Significant springback occurs after tube rotary-draw-bending （RDB）, especially for a high-strength Ti-3A1-2.5V tube （HSTT） due to its high ratio of yield strength to Young＇s modulus. The combination scheme of explicit and implicit is preferred to predict the springback. This simulation strategy includes several numerical parameters, such as element type, number of elements through thickness （NEL）, element size, etc. However, the influences of these parameters on spring- back prediction accuracy are not fully understood. Thus, taking the geometrical specification 9.525 mm × 0.508 mm ofa HSTT as the objective, the effects of numerical parameters on prediction accuracy and computation efficiency of springback simulation of HSTT RDB are investigated. The simulated springback results are compared with experimental ones. The main results are： （1） solid and continuum-shell elements predict the experimental results well; （2） for C3DSR elements, NEL of at least 3 is required to obtain reliable results and a relative error of 29% can occur as NEL is varied in the range of 1-3; （3） specifying damping factor typically works well in Abaqus/Emplicit simulation of springback and the springback results are sensitive to the magnitude of damping factor. In addition, the explanations of the effect rules are given and a guideline is added.展开更多
In vibration active control of composite structures, piezoelectricsensors/actuators are usually bonded to the surface of a host structure. Debonding of piezoelectricsensors/actuators can result in significant changes ...In vibration active control of composite structures, piezoelectricsensors/actuators are usually bonded to the surface of a host structure. Debonding of piezoelectricsensors/actuators can result in significant changes to the static and dynamic response. In thepresent paper, an novel Enhanced Assumed Strain(EAS) piezoelectric solid element formulation isdeveloped for vibration active control of laminated structures bonded with piezoelectric sensors andactuators. Unlike the conventional brick elements, the present formulation is very reliable, moreaccurate, and computationally efficient and can be used to model the response of shell structuresbesides thin plates. Delaminations are modeled by pairs of nodes with the same coordinates butdifferent node numbers, and numerical results demonstrate the performance of the element and theglobal and local effects of debonding sensors/actuators on the dynamics of the adaptive laminates.展开更多
The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix ar...The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demon- strated in our calculation results.展开更多
Conductor and suction anchor are the key equipment providing bearing capacity in the field of deep-water drilling or offshore engineering,which have the advantages of high operation efficiency and short construction p...Conductor and suction anchor are the key equipment providing bearing capacity in the field of deep-water drilling or offshore engineering,which have the advantages of high operation efficiency and short construction period.In order to drill a horizontal well in the shallow hydrate reservoir in the deep water,the suction anchor wellhead assembly is employed to undertake the main vertical bearing capacity in the second round of hydrate trial production project,so as to reduce the conductor running depth and heighten the kick-off point position.However,the deformation law of the deep-water suction anchor wellhead assembly under the moving load of the riser is not clear,and it is necessary to understand the lateral bearing characteristics to guide the design of its structural scheme.Based on 3D solid finite element method,the solid finite element model of the suction anchor wellhead assembly is established.In the model,the seabed soil is divided into seven layers,the contact between the wellhead assembly and the soil is simulated,and the vertical load and bending moment are applied to the wellhead node to simulate the riser movement when working in the deep water.The lateral bearing stability of conventional wellhead assembly and suction anchor wellhead assembly under the influence of wellhead load is discussed.The analysis results show that the bending moment is the main factor affecting the lateral deformation of the wellhead string;the anti-bending performance from increasing the outer conductor diameter is better than that from increasing the conductor wall thickness;for the subsea wellhead,the suction anchor obviously improves the lateral bearing capacity and reduces the lateral deformation.The conduct of the suction anchor wellhead assembly still needs to be lowered to a certain depth that below the maximum disturbed depth to ensure the lateral bearing stability,Thus,a method for the minimum conductor running depth for the suction anchor wellhead assembly is developed.The field implementations show that compared with the first round of hydrate trial production project,the conductor running depth is increased by 9.42 m,and there is no risk of wellhead overturning during the trial production.The method for determining the minimum conductor running depth in this paper is feasible and will still play an important role in the subsequent hydrate exploration and development.展开更多
Brick masonry constructions are very common in many areas in the world and their failure in earthquakes has been the cause of many deaths. Since human safety is main issue of disaster management, people are more conce...Brick masonry constructions are very common in many areas in the world and their failure in earthquakes has been the cause of many deaths. Since human safety is main issue of disaster management, people are more concerned about the structural assessment and strengthening of those constructions. One historical brick masonry house located in Kathmandu world heritage site is modeled by FEM （finite element method） and analyzed in E1 Centro earthquake ground motions. Bricks are modeled as solid elements and the interfaces between the brick units are modeled as zero thickness joint elements. Then, non-linear analyses of the house are applied satisfying the famous Mohr-Coulomb failure criterion. The result shows that the house is very weak and sustains large deformation in El Centro 1940 Earthquake. A strengthening solution modifying the connections of existing elements and adding wooden frame inside the house can reduce the deformations significantly.展开更多
Morphology evolution of prior β grains of laser solid forming （LSF） Ti-xAl-yV （x 11,y 20） alloys from blended elemental powders is investigated. The formation mechanism of grain morphology is revealed by incorpor...Morphology evolution of prior β grains of laser solid forming （LSF） Ti-xAl-yV （x 11,y 20） alloys from blended elemental powders is investigated. The formation mechanism of grain morphology is revealed by incorporating columnar to equiaxed transition （CET） mechanism during solidification. The morphology of prior β grains of LSF Ti-6Al-yV changes from columnar to equiaxed grains with increasing element V content from 4 to 20 wt.-%. This agrees well with CET theoretical prediction. Likewise, the grain morphology of LSF Ti-xAl-2V from blended elemental powders changes from large columnar to small equiaxed with increasing Al content from 2 to 11 wt.-%. The macro-morphologies of LSF Ti-8Al-2V and Ti-11Al-2V from blended elemental powders do not agree with CET predictions. This is caused by the increased disturbance effects of mixing enthalpy with increasing Al content, generated in the alloying process of Ti, Al, and V in the molten pool.展开更多
The CFD-DEM model was developed to simulate solid exchange behavior between two half beds in a bench-scale two-dimensional dual-leg fluidized bed （DL-FB）. Power spectrum density （PSD） analysis was applied to obtai...The CFD-DEM model was developed to simulate solid exchange behavior between two half beds in a bench-scale two-dimensional dual-leg fluidized bed （DL-FB）. Power spectrum density （PSD） analysis was applied to obtain the dominant frequency （F） of the simulated differential particle number （APLR） between the two half beds. Effects of fluidization velocity （u） and bed material inventory （H） on the solid exchange behavior were studied using the CFD-DEM model. Not only snapshots of the simulated particle flow patterns using the OpenGL code but also the dominant frequency of APLR was similar to the experimental results. The simulation results show that higher fluidization velocity assists the exchange of more particles between the two half beds, but the dispersion of clusters on the bed surface into single particles decreases the cluster exchange frequency. A greater bed material inventory results in more intense cluster exchange. The cluster exchange frequency decreases with an increase of the bed material inventory.展开更多
Fluidized beds with multiple jets have widespread industrial applications. The objective of this paper is to investigate the jet interactions and hydrodynamics of a fluidized bed with multiple jets. Discrete element m...Fluidized beds with multiple jets have widespread industrial applications. The objective of this paper is to investigate the jet interactions and hydrodynamics of a fluidized bed with multiple jets. Discrete element modeling coupled with in-house CFD code GenlDLEST has been used to simulate a bed with nine jets. The results are compared with published experiments. Mono dispersed particles of size 550 ~m are used with 1.4 times the minimum fluidization velocity of the particles. Both two and three dimensional computations have been performed. To the best of our knowledge, the results presented in this paper are the first full 3D simulations of a fluidized bed performed with multiple jets. Discrepancies between the experiment and simulations are discussed in the context of the dimensionality of the simulations. The 2D solid fraction profile compares well with the experiment close to the distributor plate. At higher heights, the 2D simulation over-predicts the solid fraction profiles near the walls. The 3D simulation on the other hand is better able to capture the solid fraction profile higher up in the bed compared to that near the distributor plate. Similarly, the normalized particle velocities and the particle fluxes compare well with the experiment closer to the distributor plate for the 2D simulation and the freeboard for the 3D simulation, respectively. A lower expanded bed height is predicted in the 2D simulation compared to the 3D simulation and the experiment. The results obtained from DEM computations show that a 2D simulation can be used to capture essential jetting trends near the distributor plate regions, whereas a full scale 3D simulation is needed to capture the bubbles near the freeboard regions. These serve as validations for the experiment and help us understand the complex jet interaction and solid circulation patterns in a multiple jet fluidized bed system.展开更多
基金supported by the National Science and Technology Major Project,China(Project No.Y2019-Ⅶ-0011-0151)and the Science Center for Gas Turbine Project(Project No.P2022-C-Ⅳ-002-001).
文摘Alloy composition design usually contributes to eliminating cracking in Ni-based superalloys during addi-tive manufacturing(AM).However,a detailed understanding of each solid solution element in the crack-ing susceptibility of Ni-based superalloys during AM is still missing.Thirteen newly designed alloys are considered to investigate the combined effect of solid solution elements on cracking susceptibility.The behaviors of solidification cracking,liquation cracking,and solid-state cracking were analyzed by the microstructural characterization and thermodynamic calculations.The results showed that W and Mo cause the formation of high melting-point carbides at grain boundaries(GBs),which increase solidifica-tion cracking susceptibility.Moreover,W and Mo lead to a slightly higher solidification cracking index(SCI)compared to Co,Cr,and Re.In the successive solidification and remelting process,the borides en-riched in W,Mo,and B around GBs will cause grain boundary segregation and liquation cracking.W and Re extend the freezing range(FR)and exacerbate the segregation of Al and Ti in the inter-dendritic re-gions,contributing to the formation of eutectics.Similarly,complete or partial melting of the eutectic can induce liquation cracking during the thermal cycling in AM.The solid-state cracking susceptibility can be reduced by solid solution elements,especially Re and Co.In summary,compared to Co,Cr,and Re,W and Mo exacerbate the cracking susceptibility.
文摘In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.
文摘Truck frames should be designed and fabricated with enough rigidity to avoid excessive deflections. Finite element analysis (FEA) plays an important role in all stages of frame designs. While being accurate, 3D solid element FEA models are built upon frame configuration details which are not feasible in the preliminary design stage, partially because of limited available design data of frames and heavy computation costs. This research develops 1D beam element FEA models for simulating frame structures. In this paper, the CAD model of a truck frame is first created. The solid element FEA analysis, which is adopted as the baseline in this study, is subsequently conducted for the stiffness of the frame, Next, beam element FEA analysis is performed for validating the feasibility of the beam element FEA model by comparing the results from the solid and beam element FEA models. It is found that the beam element FEA model can predict the frame stiffness with acceptable accuracy and reduce the computation cost significantly.
基金the National Natural Science Foundation of China (No.51275415)Program for New Century Excellent Talents in University+1 种基金the fund of the State Key Laboratory of Solidifcation Processing in NWPUNatural Science Basic Research Plan in Shaanxi Province (No.2011JQ6004),and the 111 Project (No.B08040) for the support
文摘Significant springback occurs after tube rotary-draw-bending （RDB）, especially for a high-strength Ti-3A1-2.5V tube （HSTT） due to its high ratio of yield strength to Young＇s modulus. The combination scheme of explicit and implicit is preferred to predict the springback. This simulation strategy includes several numerical parameters, such as element type, number of elements through thickness （NEL）, element size, etc. However, the influences of these parameters on spring- back prediction accuracy are not fully understood. Thus, taking the geometrical specification 9.525 mm × 0.508 mm ofa HSTT as the objective, the effects of numerical parameters on prediction accuracy and computation efficiency of springback simulation of HSTT RDB are investigated. The simulated springback results are compared with experimental ones. The main results are： （1） solid and continuum-shell elements predict the experimental results well; （2） for C3DSR elements, NEL of at least 3 is required to obtain reliable results and a relative error of 29% can occur as NEL is varied in the range of 1-3; （3） specifying damping factor typically works well in Abaqus/Emplicit simulation of springback and the springback results are sensitive to the magnitude of damping factor. In addition, the explanations of the effect rules are given and a guideline is added.
文摘In vibration active control of composite structures, piezoelectricsensors/actuators are usually bonded to the surface of a host structure. Debonding of piezoelectricsensors/actuators can result in significant changes to the static and dynamic response. In thepresent paper, an novel Enhanced Assumed Strain(EAS) piezoelectric solid element formulation isdeveloped for vibration active control of laminated structures bonded with piezoelectric sensors andactuators. Unlike the conventional brick elements, the present formulation is very reliable, moreaccurate, and computationally efficient and can be used to model the response of shell structuresbesides thin plates. Delaminations are modeled by pairs of nodes with the same coordinates butdifferent node numbers, and numerical results demonstrate the performance of the element and theglobal and local effects of debonding sensors/actuators on the dynamics of the adaptive laminates.
基金supported by the Funds for Creative Research Groups of China (50921001)the State Key Development Program for Basic Research of China (2010CB832704)
文摘The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demon- strated in our calculation results.
基金This research was jointly supported by the National Key R&D Program of China(2021YFC2800801)the Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory(Guangzhou)(GML2019ZD0501)+3 种基金the Guangzhou Basic and Applied Basic Foundation(202102020611)the project of Guangzhou Marine Geological Survey of the China Geology Survey(DD20221700)the Key-Area Research and Development Program of Guangdong Province(2020B1111030003)the High-tech Ship Research Project of the Ministry of Industry and Information Technology(CJ05N20).
文摘Conductor and suction anchor are the key equipment providing bearing capacity in the field of deep-water drilling or offshore engineering,which have the advantages of high operation efficiency and short construction period.In order to drill a horizontal well in the shallow hydrate reservoir in the deep water,the suction anchor wellhead assembly is employed to undertake the main vertical bearing capacity in the second round of hydrate trial production project,so as to reduce the conductor running depth and heighten the kick-off point position.However,the deformation law of the deep-water suction anchor wellhead assembly under the moving load of the riser is not clear,and it is necessary to understand the lateral bearing characteristics to guide the design of its structural scheme.Based on 3D solid finite element method,the solid finite element model of the suction anchor wellhead assembly is established.In the model,the seabed soil is divided into seven layers,the contact between the wellhead assembly and the soil is simulated,and the vertical load and bending moment are applied to the wellhead node to simulate the riser movement when working in the deep water.The lateral bearing stability of conventional wellhead assembly and suction anchor wellhead assembly under the influence of wellhead load is discussed.The analysis results show that the bending moment is the main factor affecting the lateral deformation of the wellhead string;the anti-bending performance from increasing the outer conductor diameter is better than that from increasing the conductor wall thickness;for the subsea wellhead,the suction anchor obviously improves the lateral bearing capacity and reduces the lateral deformation.The conduct of the suction anchor wellhead assembly still needs to be lowered to a certain depth that below the maximum disturbed depth to ensure the lateral bearing stability,Thus,a method for the minimum conductor running depth for the suction anchor wellhead assembly is developed.The field implementations show that compared with the first round of hydrate trial production project,the conductor running depth is increased by 9.42 m,and there is no risk of wellhead overturning during the trial production.The method for determining the minimum conductor running depth in this paper is feasible and will still play an important role in the subsequent hydrate exploration and development.
文摘Brick masonry constructions are very common in many areas in the world and their failure in earthquakes has been the cause of many deaths. Since human safety is main issue of disaster management, people are more concerned about the structural assessment and strengthening of those constructions. One historical brick masonry house located in Kathmandu world heritage site is modeled by FEM （finite element method） and analyzed in E1 Centro earthquake ground motions. Bricks are modeled as solid elements and the interfaces between the brick units are modeled as zero thickness joint elements. Then, non-linear analyses of the house are applied satisfying the famous Mohr-Coulomb failure criterion. The result shows that the house is very weak and sustains large deformation in El Centro 1940 Earthquake. A strengthening solution modifying the connections of existing elements and adding wooden frame inside the house can reduce the deformations significantly.
基金supported by the State Key Laboratory of Solidification Processing in NWPU (Nos.SKLSP201102 and 06-BZ-2010)Lthe China Postdoc-toral Science Foundation (No.20100470040)the National Natural Science Foundation of China (No.50871089)
文摘Morphology evolution of prior β grains of laser solid forming （LSF） Ti-xAl-yV （x 11,y 20） alloys from blended elemental powders is investigated. The formation mechanism of grain morphology is revealed by incorporating columnar to equiaxed transition （CET） mechanism during solidification. The morphology of prior β grains of LSF Ti-6Al-yV changes from columnar to equiaxed grains with increasing element V content from 4 to 20 wt.-%. This agrees well with CET theoretical prediction. Likewise, the grain morphology of LSF Ti-xAl-2V from blended elemental powders changes from large columnar to small equiaxed with increasing Al content from 2 to 11 wt.-%. The macro-morphologies of LSF Ti-8Al-2V and Ti-11Al-2V from blended elemental powders do not agree with CET predictions. This is caused by the increased disturbance effects of mixing enthalpy with increasing Al content, generated in the alloying process of Ti, Al, and V in the molten pool.
基金the support provided by the National Science and Technology Support Program of China(No.2012BAA02B00)
文摘The CFD-DEM model was developed to simulate solid exchange behavior between two half beds in a bench-scale two-dimensional dual-leg fluidized bed （DL-FB）. Power spectrum density （PSD） analysis was applied to obtain the dominant frequency （F） of the simulated differential particle number （APLR） between the two half beds. Effects of fluidization velocity （u） and bed material inventory （H） on the solid exchange behavior were studied using the CFD-DEM model. Not only snapshots of the simulated particle flow patterns using the OpenGL code but also the dominant frequency of APLR was similar to the experimental results. The simulation results show that higher fluidization velocity assists the exchange of more particles between the two half beds, but the dispersion of clusters on the bed surface into single particles decreases the cluster exchange frequency. A greater bed material inventory results in more intense cluster exchange. The cluster exchange frequency decreases with an increase of the bed material inventory.
文摘Fluidized beds with multiple jets have widespread industrial applications. The objective of this paper is to investigate the jet interactions and hydrodynamics of a fluidized bed with multiple jets. Discrete element modeling coupled with in-house CFD code GenlDLEST has been used to simulate a bed with nine jets. The results are compared with published experiments. Mono dispersed particles of size 550 ~m are used with 1.4 times the minimum fluidization velocity of the particles. Both two and three dimensional computations have been performed. To the best of our knowledge, the results presented in this paper are the first full 3D simulations of a fluidized bed performed with multiple jets. Discrepancies between the experiment and simulations are discussed in the context of the dimensionality of the simulations. The 2D solid fraction profile compares well with the experiment close to the distributor plate. At higher heights, the 2D simulation over-predicts the solid fraction profiles near the walls. The 3D simulation on the other hand is better able to capture the solid fraction profile higher up in the bed compared to that near the distributor plate. Similarly, the normalized particle velocities and the particle fluxes compare well with the experiment closer to the distributor plate for the 2D simulation and the freeboard for the 3D simulation, respectively. A lower expanded bed height is predicted in the 2D simulation compared to the 3D simulation and the experiment. The results obtained from DEM computations show that a 2D simulation can be used to capture essential jetting trends near the distributor plate regions, whereas a full scale 3D simulation is needed to capture the bubbles near the freeboard regions. These serve as validations for the experiment and help us understand the complex jet interaction and solid circulation patterns in a multiple jet fluidized bed system.
