This study developed a numerical model to efficiently treat solid waste magnesium nitrate hydrate through multi-step chemical reactions.The model simulates two-phase flow,heat,and mass transfer processes in a pyrolysi...This study developed a numerical model to efficiently treat solid waste magnesium nitrate hydrate through multi-step chemical reactions.The model simulates two-phase flow,heat,and mass transfer processes in a pyrolysis furnace to improve the decomposition rate of magnesium nitrate.The performance of multi-nozzle and single-nozzle injection methods was evaluated,and the effects of primary and secondary nozzle flow ratios,velocity ratios,and secondary nozzle inclination angles on the decomposition rate were investigated.Results indicate that multi-nozzle injection has a higher conversion efficiency and decomposition rate than single-nozzle injection,with a 10.3%higher conversion rate under the design parameters.The decomposition rate is primarily dependent on the average residence time of particles,which can be increased by decreasing flow rate and velocity ratios and increasing the inclination angle of secondary nozzles.The optimal parameters are injection flow ratio of 40%,injection velocity ratio of 0.6,and secondary nozzle inclination of 30°,corresponding to a maximum decomposition rate of 99.33%.展开更多
Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile fac...Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile facilities such as embassies and petrochemical plants,this concern now extends to a wider array of infrastructures and facilities.Engineers and scholars increasingly prioritize structural safety against explosions,particularly to prevent disproportionate collapse and damage to nearby structures.Urbanization has further amplified the reliance on oil and gas pipelines,making them vital for urban life and prime targets for terrorist activities.Consequently,there is a growing imperative for computational engineering solutions to tackle blast loading on pipelines and mitigate associated risks to avert disasters.In this study,an empty pipe model was successfully validated under contact blast conditions using Abaqus software,a powerful tool in mechanical engineering for simulating blast effects on buried pipelines.Employing a Eulerian-Lagrangian computational fluid dynamics approach,the investigation extended to above-surface and below-surface blasts at standoff distances of 25 and 50 mm.Material descriptions in the numerical model relied on Abaqus’default mechanical models.Comparative analysis revealed varying pipe performance,with deformation decreasing as explosion-to-pipe distance increased.The explosion’s location relative to the pipe surface notably influenced deformation levels,a key finding highlighted in the study.Moreover,quantitative findings indicated varying ratios of plastic dissipation energy(PDE)for different blast scenarios compared to the contact blast(P0).Specifically,P1(25 mm subsurface blast)and P2(50 mm subsurface blast)showed approximately 24.07%and 14.77%of P0’s PDE,respectively,while P3(25 mm above-surface blast)and P4(50 mm above-surface blast)exhibited lower PDE values,accounting for about 18.08%and 9.67%of P0’s PDE,respectively.Utilising energy-absorbing materials such as thin coatings of ultra-high-strength concrete,metallic foams,carbon fiber-reinforced polymer wraps,and others on the pipeline to effectively mitigate blast damage is recommended.This research contributes to the advancement of mechanical engineering by providing insights and solutions crucial for enhancing the resilience and safety of underground pipelines in the face of blast events.展开更多
On the basis of computational fluid dynamics,the flow field characteristics of multi-trophic artificial reefs,including the flow field distribution features of a single reef under three different velocities and the ef...On the basis of computational fluid dynamics,the flow field characteristics of multi-trophic artificial reefs,including the flow field distribution features of a single reef under three different velocities and the effect of spacing between reefs on flow scale and the flow state,were analyzed.Results indicate upwelling,slow flow,and eddy around a single reef.Maximum velocity,height,and volume of upwelling in front of a single reef were positively correlated with inflow velocity.The length and volume of slow flow increased with the increase in inflow velocity.Eddies were present both inside and backward,and vorticity was positively correlated with inflow velocity.Space between reefs had a minor influence on the maximum velocity and height of upwelling.With the increase in space from 0.5 L to 1.5 L(L is the reef lehgth),the length of slow flow in the front and back of the combined reefs increased slightly.When the space was 2.0 L,the length of the slow flow decreased.In four different spaces,eddies were present inside and at the back of each reef.The maximum vorticity was negatively correlated with space from 0.5 L to 1.5 L,but under 2.0 L space,the maximum vorticity was close to the vorticity of a single reef under the same inflow velocity.展开更多
Computational fluid dynamics(CFD)provides a powerful tool for investigating complicated fluid flows.This paper aims to study the applicability of CFD in the preliminary design of linear and nonlinear fluid viscous dam...Computational fluid dynamics(CFD)provides a powerful tool for investigating complicated fluid flows.This paper aims to study the applicability of CFD in the preliminary design of linear and nonlinear fluid viscous dampers.Two fluid viscous dampers were designed based on CFD models.The first device was a linear viscous damper with straight orifices.The second was a nonlinear viscous damper containing a one-way pressure-responsive valve inside its orifices.Both dampers were detailed based on CFD simulations,and their internal fluid flows were investigated.Full-scale specimens of both dampers were manufactured and tested under dynamic loads.According to the tests results,both dampers demonstrate stable cyclic behaviors,and as expected,the nonlinear damper generally tends to dissipate more energy compared to its linear counterpart.Good compatibility was achieved between the experimentally measured damper force-velocity curves and those estimated from CFD analyses.Using a thermography camera,a rise in temperature of the dampers was measured during the tests.It was found that output force of the manufactured devices was virtually independent of temperature even during long duration loadings.Accordingly,temperature dependence can be ignored in CFD models,because a reliable temperature compensator mechanism was used(or intended to be used)by the damper manufacturer.展开更多
This paper presents a time-efficient numerical approach to modelling high explosive(HE)blastwave propagation using Computational Fluid Dynamics(CFD).One of the main issues of using conventional CFD modelling in high e...This paper presents a time-efficient numerical approach to modelling high explosive(HE)blastwave propagation using Computational Fluid Dynamics(CFD).One of the main issues of using conventional CFD modelling in high explosive simulations is the ability to accurately define the initial blastwave properties that arise from the ignition and consequent explosion.Specialised codes often employ Jones-Wilkins-Lee(JWL)or similar equation of state(EOS)to simulate blasts.However,most available CFD codes are limited in terms of EOS modelling.They are restrictive to the Ideal Gas Law(IGL)for compressible flows,which is generally unsuitable for blast simulations.To this end,this paper presents a numerical approach to simulate blastwave propagation for any generic CFD code using the IGL EOS.A new method known as the Input Cavity Method(ICM)is defined where input conditions of the high explosives are given in the form of pressure,velocity and temperature time-history curves.These time history curves are input at a certain distance from the centre of the charge.It is shown that the ICM numerical method can accurately predict over-pressure and impulse time history at measured locations for the incident,reflective and complex multiple reflection scenarios with high numerical accuracy compared to experimental measurements.The ICM is compared to the Pressure Bubble Method(PBM),a common approach to replicating initial conditions for a high explosive in Finite Volume modelling.It is shown that the ICM outperforms the PBM on multiple fronts,such as peak values and overall overpressure curve shape.Finally,the paper also presents the importance of choosing an appropriate solver between the Pressure Based Solver(PBS)and Density-Based Solver(DBS)and provides the advantages and disadvantages of either choice.In general,it is shown that the PBS can resolve and capture the interactions of blastwaves to a higher degree of resolution than the DBS.This is achieved at a much higher computational cost,showing that the DBS is much preferred for quick turnarounds.展开更多
Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the productio...Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.展开更多
Background:The assessment of renal function is important to the prognosis of patients needing Fontan palliation due to the reconstructed compromised circulation.To know the relationship between the kidney perfusion an...Background:The assessment of renal function is important to the prognosis of patients needing Fontan palliation due to the reconstructed compromised circulation.To know the relationship between the kidney perfusion and hemodynamic characteristics during surgical design could reduce the risk of acute kidney injury(AKI)and the postoperative complications.However,the issue is still unsolved because the current clinical evaluation methods are unable to predict the hemodynamic changes in renal artery(RA).Methods:We reconstructed a three-dimensional(3D)vascular model of a patient requiring Fontan palliation.The technique of computational fluid dynamics(CFD)was utilized to explore the changes of RA hemodynamics under different possible blood flow rates.The relationship between the kidney perfusion and hemodynamic characteristics was investigated.Results:The calculated results indicated the declined tendency of the pressure and pressure drop as the flow rate decreased.When the flow rate decreased to two-thirds of its baseline,both the pressure of left renal artery(LRA)and the pressure of right renal artery(RRA)dipped below 50%,and the pressure of RRA fell more quickly than that of LRA.Uneven distribution of WSS was observed on the trunk of RA,and the lowest WSS was found at the distal of RA.The average WSS in RA dropped to around 50%as the flow rate reached one-third of its baseline.Conclusions:As a promising approach,CFD can be utilized to quantitatively evaluate the hemodynamic characteristics of RA and contribute to offsetting the drawbacks of clinical assessments of renal function,to help realize better prognosis for the patients with Fontan palliation.展开更多
Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simpli...Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simplified physical model and classic Reynolds equation are always applied. While the application of the general computational fluid dynamics (CFD)-fluid structure interaction (FSI) techniques is more beneficial for analysis of the fluid field in a journal bearing when more detailed solutions are needed. This paper deals with the quasi-coupling calculation of transient fluid dynamics of oil film in journal bearings and rotor dynamics with CFD-FSI techniques. The fluid dynamics of oil film is calculated by applying the so-called "dynamic mesh" technique. A new mesh movement approacb is presented while the dynamic mesh models provided by FLUENT are not suitable for the transient oil flow in journal bearings. The proposed mesh movement approach is based on the structured mesh. When the joumal moves, the movement distance of every grid in the flow field of bearing can be calculated, and then the update of the volume mesh can be handled automatically by user defined function (UDF). The journal displacement at each time step is obtained by solving the moving equations of the rotor-bearing system under the known oil film force condition. A case study is carried out to calculate the locus of the journal center and pressure distribution of the journal in order to prove the feasibility of this method. The calculating results indicate that the proposed method can predict the transient flow field of a journal bearing in a rotor-bearing system where more realistic models are involved. The presented calculation method provides a basis for studying the nonlinear dynamic behavior of a general rotor-bearing system.展开更多
A three-dimensional model for gas-solid flow in a circulating fluidized bed(CFB) riser was developed based on computational particle fluid dynamics(CPFD).The model was used to simulate the gas-solid flow behavior ...A three-dimensional model for gas-solid flow in a circulating fluidized bed(CFB) riser was developed based on computational particle fluid dynamics(CPFD).The model was used to simulate the gas-solid flow behavior inside a circulating fluidized bed riser operating at various superficial gas velocities and solids mass fluxes in two fluidization regimes,a dilute phase transport(DPT) regime and a fast fluidization(FF) regime.The simulation results were evaluated based on comparison with experimental data of solids velocity and holdup,obtained from non-invasive automated radioactive particle tracking and gamma-ray tomography techniques,respectively.The agreement of the predicted solids velocity and holdup with experimental data validated the CPFD model for the CFB riser.The model predicted the main features of the gas-solid flows in the two regimes;the uniform dilute phase in the DPT regime,and the coexistence of the dilute phase in the upper region and the dense phase in the lower region in the FF regime.The clustering and solids back mixing in the FF regime were stronger than those in the DPT regime.展开更多
It is generally recognized that internal-loop reactors are well-developed mass and heat-transfer multiphase flow reactors. However, the internal flow field in the internal-loop reactor is influenced by the structure p...It is generally recognized that internal-loop reactors are well-developed mass and heat-transfer multiphase flow reactors. However, the internal flow field in the internal-loop reactor is influenced by the structure parameter of the reactor, which has a great effect on the reaction efficiency. In this study, the computational fluid dynamics simulation method was used to determine the influence of reactor structure on flow field, and a volume-offluid model was employed to simulate the gas–liquid, two-phase flow of the internal-loop micro-electrolysis reactor. Hydrodynamic factors were optimized when the height-to-diameter ratio was 4:1, diameter ratio was9:1, draft-tube axial height was 90 mm. Three-dimensional simulations for the water distributor were carried out, and the results suggested that the optimal conditions are as follows: the number of water distribution pipes was four, and an inhomogeneous water distribution was used. According to the results of the simulation,the suitable structure can be used to achieve good fluid mechanical properties, such as the good liquid circulation velocity and gas holdup, which provides a good theoretical foundation for the application of the reactor.展开更多
The effect of mixing on the precipitation of barium sulfate in a continuous stirred tank is simulated numerically with different feeding location, feed concentration, impeller speed and residence time through solving ...The effect of mixing on the precipitation of barium sulfate in a continuous stirred tank is simulated numerically with different feeding location, feed concentration, impeller speed and residence time through solving the standard momentum and mass transport equations in combination with the moment equations for crystal population balance. The numerical method was validated with the literature data. The simulation results including the distribution of the local supersaturation ratio distribution in the precipitator, mean crystal size and coefficient of variation under different operating conditions compared well with experimental data in the literature. The effect of the presence of a draft tube on precipitation were also investigated, and it is suggested that the installation of a draft tube increased the mean crystal size, in general agreement with experimental work in the literature.展开更多
Many fishes use undulatory fin to propel themselves in the underwater environment. These locomotor mechanisms have a popular interest to many researchers. In the present study, we perform a three-dimensional unsteady ...Many fishes use undulatory fin to propel themselves in the underwater environment. These locomotor mechanisms have a popular interest to many researchers. In the present study, we perform a three-dimensional unsteady computation of an undulatory mechanical fin that is driven by Shape Memory Alloy (SMA). The objective of the computation is to investigate the fluid dynamics of force production associated with the undulatory mechanical fin. An unstructured, grid-based, unsteady Navier-Stokes solver with automatic adaptive remeshing is used to compute the unsteady flow around the fin through five complete cycles. The pressure distribution on fin surface is computed and integrated to provide fin forces which are decomposed into lift and thrust. The velocity field is also computed throughout the swimming cycle. Finally, a comparison is conducted to reveal the dynamics of force generation according to the kinematic parameters of the undulatory fin (amplitude, frequency and wavelength).展开更多
Disturbance in wind regime and sand erosion deposition balance may lead to burial and eventual vanishing of a site.This study conducted 3D computational fluid dynamics(CFD)simulations to evaluate the effect of a propo...Disturbance in wind regime and sand erosion deposition balance may lead to burial and eventual vanishing of a site.This study conducted 3D computational fluid dynamics(CFD)simulations to evaluate the effect of a proposed city design on the wind environment of the Crescent Spring,a downwind natural heritage site located in Dunhuang,Northwestern China.Satellite terrain data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer(ASTER)Digital Elevation Model(DEM)were used to construct the solid surface model.Steady-state Reynolds Averaged Navier-Stokes equations(RANS)with shear stress transport(SST)k-ωturbulence model were then applied to solve the flow field problems.Land-use changes were modeled implicitly by dividing the underlying surface into different areas and by applying corresponding aerodynamic roughness lengths.Simulations were performed by using cases with different city areas and building heights.Results show that the selected model could capture the surface roughness changes and could adjust wind profile over a large area.Wind profiles varied over the greenfield to the north and over the Gobi land to the east of the spring.Therefore,different wind speed reduction effects were observed from various city construction scenarios.The current city design would lead to about 2 m/s of wind speed reduction at the downwind city edge and about 1 m/s of wind speed reduction at the north of the spring at 35-m height.Reducing the city height in the north greenfield area could efficiently eliminate the negative effects of wind spee.By contrast,restricting the city area worked better in the eastern Gobi area compared with other parts of the study area.Wind speed reduction in areas near the spring could be limited to 0.1 m/s by combining these two abatement strategies.The CFD method could be applied to simulate the wind environment affected by other land-use changes over a large terrain.展开更多
The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model...The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, and a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall; the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.展开更多
Computational fluid dynamics(CFD) simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds num...Computational fluid dynamics(CFD) simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.展开更多
As a result of environmental degradation,urban green space has become a key issue for urban sustainable development.This paper takes Liaoyang City in Northeast China as an example to develop green space planning using...As a result of environmental degradation,urban green space has become a key issue for urban sustainable development.This paper takes Liaoyang City in Northeast China as an example to develop green space planning using the computational fluid dynamics (CFD) model,landscape ecological principles and Geographical Information System (GIS).Based on the influencing factors of topography,building density and orientation,Shou Mountain,Longding Mountain and the Taizi River were selected as the urban ventilation paths to promote wind and oxygen circulation.Oxygen concentration around the green spaces gradually decreased with wind speed increase and wind direction change.There were obvious negative correlation relationships between the oxygen dispersion concentration and urban layout factors such as the building plot ratio and building density.Comparison with the field measurements found that there was significant correlation relationship between simulated oxygen concentration and field measurements (R 2=0.6415,p<0.001),moreover,simulation precision was higher than 92%,which indicated CFD model was effective for urban oxygen concentration simulation.Only less than 10% areas in Liaoyang City proper needed more green space urgently to improve oxygen concentration,mainly concentrated in Baitai and west Wensheng districts.Based on land-scape ecology principle,green space planning at different spatial scales were proposed to create a green space network system for Liaoyang City,including features such as green wedges,green belts and parks.Totally,about 2012 ha of green space need to be constructed as oxygen sources and ventilation paths.Compared with the current green space pattern,proposed green space planning could improve oxygen concentration obviously.The CFD model and research results in this paper could provide an effective way and theory support for sustainable development of urban green space.展开更多
Three approaches based on computational fluid dynamics(CFD) techniques have been assessed for their ability to describe the engineering flow environment in a miniaturized mechanically agitated bioreactor. The three a...Three approaches based on computational fluid dynamics(CFD) techniques have been assessed for their ability to describe the engineering flow environment in a miniaturized mechanically agitated bioreactor. The three approaches tested were the source-sink(SS), the multiple reference frames(MRF) and the sliding grids(SG). In all the cases, the predictions of the velocity components agree with reported experimental data. However, the analysis of the results of the turbulent intensities predicted by the three approaches indicates the MRF and the SG techniques under predicted turbulent intensities are comparable to both experimental measurements and the SS method. The predicted power number and pumping number based on the SS approach are closer to typical reported experimental values compared to those obtained from the MRF and SG methods.展开更多
Computational fluid dynamics (CFD) codes are being increasingly used in the simulation of submarine oil spills. This study focuses on the process of oil spills, from damaged submarine pipes, to the sea surface, usin...Computational fluid dynamics (CFD) codes are being increasingly used in the simulation of submarine oil spills. This study focuses on the process of oil spills, from damaged submarine pipes, to the sea surface, using numerical models. The underwater oil spill model is developed, and a description of the governing equations is proposed, along with modifications required for the particalization of the control volume. Available experimental data were introduced to evaluate the validity of the CFD predictions, the results of which proved to be in good agreement with the experimental data. The effects of oil leak rate, leak diameter, current velocity, and oil density are investigated, by the validated CFD model, to estimate the undersea leakage time, the lateral migration distance, and surface diffusion range when the oil reaches the sea surface. Results indicate that the leakage time and lateral migration distance increase with decreasing leak rates and leak diameter, and increase with increasing current velocity and oil density. On the other hand, a large leak diameter, high density, high leak rate, or fast currents result in a greater surface diffusion range. The findings and analysis presented here will provide practical predictions of oil spills, and guidance for emergency rescues.展开更多
The liquid flow on a single-pass sieve distillation tray is simulated with a three-dimensional computational fluid dynamics (CFD) program with the K-ε turbulence model. In the model, a source term SMi is formulatedin...The liquid flow on a single-pass sieve distillation tray is simulated with a three-dimensional computational fluid dynamics (CFD) program with the K-ε turbulence model. In the model, a source term SMi is formulatedin the Navier-Stokes equations to represent the interfacial momentum transfer and another term Sc is added to themass transfer equation as the source of interfacial mass transfer. The simulation provides the detailed informationof the three-dimensional distribution of liquid velocity on the tray, the circulation area and the concentration profilealong the height of liquid layer.展开更多
Lateral deformation of liquefiable soil is a cause of much damage during earthquakes, reportedly more than other forms of liquefaction-induced ground failures. Researchers have presented studies in which the liquefied...Lateral deformation of liquefiable soil is a cause of much damage during earthquakes, reportedly more than other forms of liquefaction-induced ground failures. Researchers have presented studies in which the liquefied soil is considered as viscous fluid. In this manner, the liquefied soil behaves as non-Newtonian fluid, whose viscosity decreases as the shear strain rate increases. The current study incorporates computational fluid dynamics to propose a simplified dynamic analysis for the liquefaction-induced lateral deformation of earth slopes. The numerical procedure involves a quasi-linear elastic model for small to moderate strains and a Bingham fluid model for large strain states during liquefaction. An iterative procedure is considered to estimate the strain-compatible shear stiffness of soil. The post-liquefaction residual strength of soil is considered as the initial Bingham viscosity. Performance of the numerical procedure is examined by using the results of centrifuge model and shaking table tests together with some field observations of lateral ground deformation. The results demonstrate that the proposed procedure predicts the time history of lateral ground deformation with a reasonable degree of precision.展开更多
基金the financial support for this work provided by the National Key R&D Program of China‘Technologies and Integrated Application of Magnesite Waste Utilization for High-Valued Chemicals and Materials’(2020YFC1909303)。
文摘This study developed a numerical model to efficiently treat solid waste magnesium nitrate hydrate through multi-step chemical reactions.The model simulates two-phase flow,heat,and mass transfer processes in a pyrolysis furnace to improve the decomposition rate of magnesium nitrate.The performance of multi-nozzle and single-nozzle injection methods was evaluated,and the effects of primary and secondary nozzle flow ratios,velocity ratios,and secondary nozzle inclination angles on the decomposition rate were investigated.Results indicate that multi-nozzle injection has a higher conversion efficiency and decomposition rate than single-nozzle injection,with a 10.3%higher conversion rate under the design parameters.The decomposition rate is primarily dependent on the average residence time of particles,which can be increased by decreasing flow rate and velocity ratios and increasing the inclination angle of secondary nozzles.The optimal parameters are injection flow ratio of 40%,injection velocity ratio of 0.6,and secondary nozzle inclination of 30°,corresponding to a maximum decomposition rate of 99.33%.
文摘Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile facilities such as embassies and petrochemical plants,this concern now extends to a wider array of infrastructures and facilities.Engineers and scholars increasingly prioritize structural safety against explosions,particularly to prevent disproportionate collapse and damage to nearby structures.Urbanization has further amplified the reliance on oil and gas pipelines,making them vital for urban life and prime targets for terrorist activities.Consequently,there is a growing imperative for computational engineering solutions to tackle blast loading on pipelines and mitigate associated risks to avert disasters.In this study,an empty pipe model was successfully validated under contact blast conditions using Abaqus software,a powerful tool in mechanical engineering for simulating blast effects on buried pipelines.Employing a Eulerian-Lagrangian computational fluid dynamics approach,the investigation extended to above-surface and below-surface blasts at standoff distances of 25 and 50 mm.Material descriptions in the numerical model relied on Abaqus’default mechanical models.Comparative analysis revealed varying pipe performance,with deformation decreasing as explosion-to-pipe distance increased.The explosion’s location relative to the pipe surface notably influenced deformation levels,a key finding highlighted in the study.Moreover,quantitative findings indicated varying ratios of plastic dissipation energy(PDE)for different blast scenarios compared to the contact blast(P0).Specifically,P1(25 mm subsurface blast)and P2(50 mm subsurface blast)showed approximately 24.07%and 14.77%of P0’s PDE,respectively,while P3(25 mm above-surface blast)and P4(50 mm above-surface blast)exhibited lower PDE values,accounting for about 18.08%and 9.67%of P0’s PDE,respectively.Utilising energy-absorbing materials such as thin coatings of ultra-high-strength concrete,metallic foams,carbon fiber-reinforced polymer wraps,and others on the pipeline to effectively mitigate blast damage is recommended.This research contributes to the advancement of mechanical engineering by providing insights and solutions crucial for enhancing the resilience and safety of underground pipelines in the face of blast events.
基金supported by the National Natural Science Foundation of China(No.32002442)the National Key R&D Program(No.2019YFD0902101).
文摘On the basis of computational fluid dynamics,the flow field characteristics of multi-trophic artificial reefs,including the flow field distribution features of a single reef under three different velocities and the effect of spacing between reefs on flow scale and the flow state,were analyzed.Results indicate upwelling,slow flow,and eddy around a single reef.Maximum velocity,height,and volume of upwelling in front of a single reef were positively correlated with inflow velocity.The length and volume of slow flow increased with the increase in inflow velocity.Eddies were present both inside and backward,and vorticity was positively correlated with inflow velocity.Space between reefs had a minor influence on the maximum velocity and height of upwelling.With the increase in space from 0.5 L to 1.5 L(L is the reef lehgth),the length of slow flow in the front and back of the combined reefs increased slightly.When the space was 2.0 L,the length of the slow flow decreased.In four different spaces,eddies were present inside and at the back of each reef.The maximum vorticity was negatively correlated with space from 0.5 L to 1.5 L,but under 2.0 L space,the maximum vorticity was close to the vorticity of a single reef under the same inflow velocity.
文摘Computational fluid dynamics(CFD)provides a powerful tool for investigating complicated fluid flows.This paper aims to study the applicability of CFD in the preliminary design of linear and nonlinear fluid viscous dampers.Two fluid viscous dampers were designed based on CFD models.The first device was a linear viscous damper with straight orifices.The second was a nonlinear viscous damper containing a one-way pressure-responsive valve inside its orifices.Both dampers were detailed based on CFD simulations,and their internal fluid flows were investigated.Full-scale specimens of both dampers were manufactured and tested under dynamic loads.According to the tests results,both dampers demonstrate stable cyclic behaviors,and as expected,the nonlinear damper generally tends to dissipate more energy compared to its linear counterpart.Good compatibility was achieved between the experimentally measured damper force-velocity curves and those estimated from CFD analyses.Using a thermography camera,a rise in temperature of the dampers was measured during the tests.It was found that output force of the manufactured devices was virtually independent of temperature even during long duration loadings.Accordingly,temperature dependence can be ignored in CFD models,because a reliable temperature compensator mechanism was used(or intended to be used)by the damper manufacturer.
文摘This paper presents a time-efficient numerical approach to modelling high explosive(HE)blastwave propagation using Computational Fluid Dynamics(CFD).One of the main issues of using conventional CFD modelling in high explosive simulations is the ability to accurately define the initial blastwave properties that arise from the ignition and consequent explosion.Specialised codes often employ Jones-Wilkins-Lee(JWL)or similar equation of state(EOS)to simulate blasts.However,most available CFD codes are limited in terms of EOS modelling.They are restrictive to the Ideal Gas Law(IGL)for compressible flows,which is generally unsuitable for blast simulations.To this end,this paper presents a numerical approach to simulate blastwave propagation for any generic CFD code using the IGL EOS.A new method known as the Input Cavity Method(ICM)is defined where input conditions of the high explosives are given in the form of pressure,velocity and temperature time-history curves.These time history curves are input at a certain distance from the centre of the charge.It is shown that the ICM numerical method can accurately predict over-pressure and impulse time history at measured locations for the incident,reflective and complex multiple reflection scenarios with high numerical accuracy compared to experimental measurements.The ICM is compared to the Pressure Bubble Method(PBM),a common approach to replicating initial conditions for a high explosive in Finite Volume modelling.It is shown that the ICM outperforms the PBM on multiple fronts,such as peak values and overall overpressure curve shape.Finally,the paper also presents the importance of choosing an appropriate solver between the Pressure Based Solver(PBS)and Density-Based Solver(DBS)and provides the advantages and disadvantages of either choice.In general,it is shown that the PBS can resolve and capture the interactions of blastwaves to a higher degree of resolution than the DBS.This is achieved at a much higher computational cost,showing that the DBS is much preferred for quick turnarounds.
基金National Natural Science Foundation of China(U2004176,22008055)Technology Research Project of Henan Province(232102240034)are gratefully acknowledged.
文摘Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.
基金Funding Statement:This study was supported by National Natural Science Foundation of China(No.81970439)Natural Science Foundation of Shanghai(No.19ZR1432700)+1 种基金Fund of the Shanghai Committee of Science and Technology(Nos.19411965400,17DZ2253100)the Development Fund of Shanghai Talents(No.2020114).
文摘Background:The assessment of renal function is important to the prognosis of patients needing Fontan palliation due to the reconstructed compromised circulation.To know the relationship between the kidney perfusion and hemodynamic characteristics during surgical design could reduce the risk of acute kidney injury(AKI)and the postoperative complications.However,the issue is still unsolved because the current clinical evaluation methods are unable to predict the hemodynamic changes in renal artery(RA).Methods:We reconstructed a three-dimensional(3D)vascular model of a patient requiring Fontan palliation.The technique of computational fluid dynamics(CFD)was utilized to explore the changes of RA hemodynamics under different possible blood flow rates.The relationship between the kidney perfusion and hemodynamic characteristics was investigated.Results:The calculated results indicated the declined tendency of the pressure and pressure drop as the flow rate decreased.When the flow rate decreased to two-thirds of its baseline,both the pressure of left renal artery(LRA)and the pressure of right renal artery(RRA)dipped below 50%,and the pressure of RRA fell more quickly than that of LRA.Uneven distribution of WSS was observed on the trunk of RA,and the lowest WSS was found at the distal of RA.The average WSS in RA dropped to around 50%as the flow rate reached one-third of its baseline.Conclusions:As a promising approach,CFD can be utilized to quantitatively evaluate the hemodynamic characteristics of RA and contribute to offsetting the drawbacks of clinical assessments of renal function,to help realize better prognosis for the patients with Fontan palliation.
基金supported by National Hi-tech Research and Development Program of China (863 Program, Grant No. 2009AA04Z413)Zhejiang Provincial Natural Science Foundation of China (Grant No. Y1110109)
文摘Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simplified physical model and classic Reynolds equation are always applied. While the application of the general computational fluid dynamics (CFD)-fluid structure interaction (FSI) techniques is more beneficial for analysis of the fluid field in a journal bearing when more detailed solutions are needed. This paper deals with the quasi-coupling calculation of transient fluid dynamics of oil film in journal bearings and rotor dynamics with CFD-FSI techniques. The fluid dynamics of oil film is calculated by applying the so-called "dynamic mesh" technique. A new mesh movement approacb is presented while the dynamic mesh models provided by FLUENT are not suitable for the transient oil flow in journal bearings. The proposed mesh movement approach is based on the structured mesh. When the joumal moves, the movement distance of every grid in the flow field of bearing can be calculated, and then the update of the volume mesh can be handled automatically by user defined function (UDF). The journal displacement at each time step is obtained by solving the moving equations of the rotor-bearing system under the known oil film force condition. A case study is carried out to calculate the locus of the journal center and pressure distribution of the journal in order to prove the feasibility of this method. The calculating results indicate that the proposed method can predict the transient flow field of a journal bearing in a rotor-bearing system where more realistic models are involved. The presented calculation method provides a basis for studying the nonlinear dynamic behavior of a general rotor-bearing system.
基金support by the National Basic Research Program (Grant No. 2010CB226906,and 2012CB215000)
文摘A three-dimensional model for gas-solid flow in a circulating fluidized bed(CFB) riser was developed based on computational particle fluid dynamics(CPFD).The model was used to simulate the gas-solid flow behavior inside a circulating fluidized bed riser operating at various superficial gas velocities and solids mass fluxes in two fluidization regimes,a dilute phase transport(DPT) regime and a fast fluidization(FF) regime.The simulation results were evaluated based on comparison with experimental data of solids velocity and holdup,obtained from non-invasive automated radioactive particle tracking and gamma-ray tomography techniques,respectively.The agreement of the predicted solids velocity and holdup with experimental data validated the CPFD model for the CFB riser.The model predicted the main features of the gas-solid flows in the two regimes;the uniform dilute phase in the DPT regime,and the coexistence of the dilute phase in the upper region and the dense phase in the lower region in the FF regime.The clustering and solids back mixing in the FF regime were stronger than those in the DPT regime.
基金Supported by the National Natural Science Foundation of China(21677018)Jointly Projects of Beijing Natural Science Foundation and Beijing Municipal Education Commission(KZ201810017024)
文摘It is generally recognized that internal-loop reactors are well-developed mass and heat-transfer multiphase flow reactors. However, the internal flow field in the internal-loop reactor is influenced by the structure parameter of the reactor, which has a great effect on the reaction efficiency. In this study, the computational fluid dynamics simulation method was used to determine the influence of reactor structure on flow field, and a volume-offluid model was employed to simulate the gas–liquid, two-phase flow of the internal-loop micro-electrolysis reactor. Hydrodynamic factors were optimized when the height-to-diameter ratio was 4:1, diameter ratio was9:1, draft-tube axial height was 90 mm. Three-dimensional simulations for the water distributor were carried out, and the results suggested that the optimal conditions are as follows: the number of water distribution pipes was four, and an inhomogeneous water distribution was used. According to the results of the simulation,the suitable structure can be used to achieve good fluid mechanical properties, such as the good liquid circulation velocity and gas holdup, which provides a good theoretical foundation for the application of the reactor.
基金Supported by the National Natural Science Foundation of China (Nos.20236050, 50134020) and the Special Funds for Major State Basic Research Program of China (973 Program, 2004CB217604).
文摘The effect of mixing on the precipitation of barium sulfate in a continuous stirred tank is simulated numerically with different feeding location, feed concentration, impeller speed and residence time through solving the standard momentum and mass transport equations in combination with the moment equations for crystal population balance. The numerical method was validated with the literature data. The simulation results including the distribution of the local supersaturation ratio distribution in the precipitator, mean crystal size and coefficient of variation under different operating conditions compared well with experimental data in the literature. The effect of the presence of a draft tube on precipitation were also investigated, and it is suggested that the installation of a draft tube increased the mean crystal size, in general agreement with experimental work in the literature.
文摘Many fishes use undulatory fin to propel themselves in the underwater environment. These locomotor mechanisms have a popular interest to many researchers. In the present study, we perform a three-dimensional unsteady computation of an undulatory mechanical fin that is driven by Shape Memory Alloy (SMA). The objective of the computation is to investigate the fluid dynamics of force production associated with the undulatory mechanical fin. An unstructured, grid-based, unsteady Navier-Stokes solver with automatic adaptive remeshing is used to compute the unsteady flow around the fin through five complete cycles. The pressure distribution on fin surface is computed and integrated to provide fin forces which are decomposed into lift and thrust. The velocity field is also computed throughout the swimming cycle. Finally, a comparison is conducted to reveal the dynamics of force generation according to the kinematic parameters of the undulatory fin (amplitude, frequency and wavelength).
基金supported by the National Basic Research Program of China(2012CB026105)the National Natural Science Foundation of China(41201003,41071009)the China Postdoctoral Science Foundation(2012M52819)
文摘Disturbance in wind regime and sand erosion deposition balance may lead to burial and eventual vanishing of a site.This study conducted 3D computational fluid dynamics(CFD)simulations to evaluate the effect of a proposed city design on the wind environment of the Crescent Spring,a downwind natural heritage site located in Dunhuang,Northwestern China.Satellite terrain data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer(ASTER)Digital Elevation Model(DEM)were used to construct the solid surface model.Steady-state Reynolds Averaged Navier-Stokes equations(RANS)with shear stress transport(SST)k-ωturbulence model were then applied to solve the flow field problems.Land-use changes were modeled implicitly by dividing the underlying surface into different areas and by applying corresponding aerodynamic roughness lengths.Simulations were performed by using cases with different city areas and building heights.Results show that the selected model could capture the surface roughness changes and could adjust wind profile over a large area.Wind profiles varied over the greenfield to the north and over the Gobi land to the east of the spring.Therefore,different wind speed reduction effects were observed from various city construction scenarios.The current city design would lead to about 2 m/s of wind speed reduction at the downwind city edge and about 1 m/s of wind speed reduction at the north of the spring at 35-m height.Reducing the city height in the north greenfield area could efficiently eliminate the negative effects of wind spee.By contrast,restricting the city area worked better in the eastern Gobi area compared with other parts of the study area.Wind speed reduction in areas near the spring could be limited to 0.1 m/s by combining these two abatement strategies.The CFD method could be applied to simulate the wind environment affected by other land-use changes over a large terrain.
基金Supported by the National 863 Project (2001AA642030-1) and Zhejiang Provincial Key Research Project (010007037).
文摘The flow field of gas and liquid in a φ150mm rotating-stream-tray (RST) scrubber is simulated by using computational fluid dynamic (CFD) method. The sismulation is based on the two-equation RNG κ-ε turbulence model, Eulerian multiphase model, and a real-shape 3D model with a huge number of meshes. The simulation results include detailed information about velocity, pressure, volume fraction and so on. Some features of the flow field are obtained: liquid is atomized in a thin annular zone; a high velocity air zone prevents water drops at the bottom from flying towards the wall; the pressure varies sharply at the end of blades and so on. The results will be helpful for structure optimization and engineering design.
基金National Natural Science Foundation of China (No.50435030)
文摘Computational fluid dynamics(CFD) simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.
基金Under the auspices of National Key Technology Research and Development Program of China (No.2008BAJ10B01-01)National Natural Science Foundation of China (No.40801069)
文摘As a result of environmental degradation,urban green space has become a key issue for urban sustainable development.This paper takes Liaoyang City in Northeast China as an example to develop green space planning using the computational fluid dynamics (CFD) model,landscape ecological principles and Geographical Information System (GIS).Based on the influencing factors of topography,building density and orientation,Shou Mountain,Longding Mountain and the Taizi River were selected as the urban ventilation paths to promote wind and oxygen circulation.Oxygen concentration around the green spaces gradually decreased with wind speed increase and wind direction change.There were obvious negative correlation relationships between the oxygen dispersion concentration and urban layout factors such as the building plot ratio and building density.Comparison with the field measurements found that there was significant correlation relationship between simulated oxygen concentration and field measurements (R 2=0.6415,p<0.001),moreover,simulation precision was higher than 92%,which indicated CFD model was effective for urban oxygen concentration simulation.Only less than 10% areas in Liaoyang City proper needed more green space urgently to improve oxygen concentration,mainly concentrated in Baitai and west Wensheng districts.Based on land-scape ecology principle,green space planning at different spatial scales were proposed to create a green space network system for Liaoyang City,including features such as green wedges,green belts and parks.Totally,about 2012 ha of green space need to be constructed as oxygen sources and ventilation paths.Compared with the current green space pattern,proposed green space planning could improve oxygen concentration obviously.The CFD model and research results in this paper could provide an effective way and theory support for sustainable development of urban green space.
基金Supported by the U CL ORS Award and KC Wong Scholarshi
文摘Three approaches based on computational fluid dynamics(CFD) techniques have been assessed for their ability to describe the engineering flow environment in a miniaturized mechanically agitated bioreactor. The three approaches tested were the source-sink(SS), the multiple reference frames(MRF) and the sliding grids(SG). In all the cases, the predictions of the velocity components agree with reported experimental data. However, the analysis of the results of the turbulent intensities predicted by the three approaches indicates the MRF and the SG techniques under predicted turbulent intensities are comparable to both experimental measurements and the SS method. The predicted power number and pumping number based on the SS approach are closer to typical reported experimental values compared to those obtained from the MRF and SG methods.
基金The National Basic Research Program(973 Program)under contract No.2014CB046803the National Natural Science Foundation of China under contract No.51239008the National Science and Technology Major Project under contract No.2016ZX05028005-004
文摘Computational fluid dynamics (CFD) codes are being increasingly used in the simulation of submarine oil spills. This study focuses on the process of oil spills, from damaged submarine pipes, to the sea surface, using numerical models. The underwater oil spill model is developed, and a description of the governing equations is proposed, along with modifications required for the particalization of the control volume. Available experimental data were introduced to evaluate the validity of the CFD predictions, the results of which proved to be in good agreement with the experimental data. The effects of oil leak rate, leak diameter, current velocity, and oil density are investigated, by the validated CFD model, to estimate the undersea leakage time, the lateral migration distance, and surface diffusion range when the oil reaches the sea surface. Results indicate that the leakage time and lateral migration distance increase with decreasing leak rates and leak diameter, and increase with increasing current velocity and oil density. On the other hand, a large leak diameter, high density, high leak rate, or fast currents result in a greater surface diffusion range. The findings and analysis presented here will provide practical predictions of oil spills, and guidance for emergency rescues.
基金Supported by Visiting Scholar Foundation of Key Lab in Tianjin University.
文摘The liquid flow on a single-pass sieve distillation tray is simulated with a three-dimensional computational fluid dynamics (CFD) program with the K-ε turbulence model. In the model, a source term SMi is formulatedin the Navier-Stokes equations to represent the interfacial momentum transfer and another term Sc is added to themass transfer equation as the source of interfacial mass transfer. The simulation provides the detailed informationof the three-dimensional distribution of liquid velocity on the tray, the circulation area and the concentration profilealong the height of liquid layer.
基金International Institute of Earthquake Engineering and Seismology(IIEES)under Grant No.AM-7/507-6 723-545
文摘Lateral deformation of liquefiable soil is a cause of much damage during earthquakes, reportedly more than other forms of liquefaction-induced ground failures. Researchers have presented studies in which the liquefied soil is considered as viscous fluid. In this manner, the liquefied soil behaves as non-Newtonian fluid, whose viscosity decreases as the shear strain rate increases. The current study incorporates computational fluid dynamics to propose a simplified dynamic analysis for the liquefaction-induced lateral deformation of earth slopes. The numerical procedure involves a quasi-linear elastic model for small to moderate strains and a Bingham fluid model for large strain states during liquefaction. An iterative procedure is considered to estimate the strain-compatible shear stiffness of soil. The post-liquefaction residual strength of soil is considered as the initial Bingham viscosity. Performance of the numerical procedure is examined by using the results of centrifuge model and shaking table tests together with some field observations of lateral ground deformation. The results demonstrate that the proposed procedure predicts the time history of lateral ground deformation with a reasonable degree of precision.