To investigate the influence of structural parameters on the performances and internal flow characteristics of partial flow pumps at a low specific speed of 10000 rpm,special attention was paid to the first and second...To investigate the influence of structural parameters on the performances and internal flow characteristics of partial flow pumps at a low specific speed of 10000 rpm,special attention was paid to the first and second stage impeller guide vanes.Moreover,the impeller blade outlet width,impeller inlet diameter,blade inclination angle,and number of blades were considered for orthogonal tests.Accordingly,nine groups of design solutions were formed,and then used as a basis for the execution of numerical simulations(CFD)aimed at obtaining the efficiency values and heads for each design solution group.The influence of impeller geometric parameters on the efficiency and head was explored,and the“weight”of each factor was obtained via a range analysis.Optimal structural parameters were finally chosen on the basis of the numerical simulation results,and the performances of the optimized model were verified accordingly(yet by means of CFD).Evidence is provided that the increase in the efficiency and head of the optimized model was 12.11%and 23.5 m,respectively,compared with those of the original model.展开更多
A self-priming pump is a centrifugal pump that has the ability to prime itself. Typically, its performance dependson the configuration of its reflux hole. In this study, the ANSYS FLUENT software is used to investigat...A self-priming pump is a centrifugal pump that has the ability to prime itself. Typically, its performance dependson the configuration of its reflux hole. In this study, the ANSYS FLUENT software is used to investigate the effectsof three different radial positions of the reflux hole on gas-liquid two-phase distribution, pressure pulsation, andimp during self-priming. The research results indicate that: (1) The effective channel size for the reflux liquid toenter the volute varies depending on the location of the reflux hole. The effect of the impeller rotation on thereflux liquid becomes more obvious as the setting distance of the reflux aperture decreases. (2) The position ofthe reflux hole significantly affects the gas phase mass fraction inside the impeller, resulting in a significant reductionin the time it takes for the mass fraction to exceed 80%. (3) The position of the reflux hole significantly affectsthe average pressure on each monitoring surface. (4) Placing the reflux hole at a excessively distant radial distancecan result in an excessive vertical component. (5) The self-priming performance of the pump can be improved tosome extent by placing the return hole at a small radial distance.展开更多
Plasma jet triggered gas gap switch has obvious advantages in fast control switch.The development of the plasma in the ambient medium is the key factor affecting the triggering conduction of the gas switch.However,the...Plasma jet triggered gas gap switch has obvious advantages in fast control switch.The development of the plasma in the ambient medium is the key factor affecting the triggering conduction of the gas switch.However,the plasma jet process and its characteristic parameters are complicated and the existing test methods cannot fully characterize its development laws.In this work,a two-dimensional transient fluid calculation model of the plasma jet process of the gas gap switch is established based on the renormalization-group k-εturbulence equation.The results show that the characteristic parameters and morphological evolution of the plasma jet are basically consistent with the experimental results,which verifies the accuracy of the simulation model calculation.The plasma jet is a long strip with an initial velocity of 1.0 km·s-1and develops in both axial and radial directions.The jet velocity fluctuates significantly with axial height.As the plasma jet enters the main gap,the pressure inside the trigger cavity drops by80%,resulting in a rapid drop in the jet velocity.When the plasma jet head interacts with the atmosphere,the two-phase fluid compresses each other,generating a forward-propelled pressure wave.The plasma jet heads flow at high velocity,a negative pressure zone is formed in the middle part of the jet,and the pressure peak decreases gradually with height.As the value of the inlet pressure increases,the characteristic parameters of the plasma jet increase.The entrainment phenomenon is evident,which leads to an increase in the pressure imbalance of the atmospheric gas medium,leading to a significant Coanda effect.Compared with air,the characteristic parameters of a plasma jet in SF6are lower,and the morphological evolution is significantly suppressed.The results of this study can provide some insight into the mechanism of action of the switch jet plasma development process.展开更多
Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, p...Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.展开更多
Based on the diurnal consecutively observed data in the offshore area of Jiaonan in 2005, the paper tries to make a preliminary analysis of the specificity of ocean currents, tidal current property and residual curren...Based on the diurnal consecutively observed data in the offshore area of Jiaonan in 2005, the paper tries to make a preliminary analysis of the specificity of ocean currents, tidal current property and residual current property in the area in observing dates. Then on the basis of observed data analysis and by employing the split-step method, the paper conducts a numerical simulation of the tidal current field, which can show the M2 tidal constituent tidal wave system, current ellipse distribution, maximum current velocity distribution and time-dependent current field. The calculated results agree well with the observed data, which can on the one hand reflect the basic specificities of temporal and spatial distribution of the M2 tidal constituent current field to some extent, and, on the other hand, offer more information about the hydrodynamic condition. So the paper would provide a scientific basis for the making of sea environment protection plans in the offshore area of Jiaonan under certain conditions.展开更多
The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their appl...The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their applications,focused on elasticity,heat conduction,electromagnetic field analysis,and fluid dynamics.The merits of the collocation method can be attributed to the need for element mesh,simple implementation,high computational efficiency,and ease in handling irregular domain problems since the collocation method is a type of node-based numerical method.Beginning with the fundamental principles of the collocation method,the discretization process in the continuous domain is elucidated,and how the collocation method approximation solutions for solving differential equations are explained.Delving into the historical development of the collocation methods,their earliest applications and key milestones are traced,thereby demonstrating their evolution within the realm of numerical computation.The mathematical foundations of collocation methods,encompassing the selection of interpolation functions,definition of weighting functions,and derivation of integration rules,are examined in detail,emphasizing their significance in comprehending the method’s effectiveness and stability.At last,the practical application of the collocation methods in engineering contexts is emphasized,including heat conduction simulations,electromagnetic coupled field analysis,and fluid dynamics simulations.These specific case studies can underscore collocation method’s broad applicability and effectiveness in addressing complex engineering challenges.In conclusion,this paper puts forward the future development trend of the collocation method through rigorous analysis and discussion,thereby facilitating further advancements in research and practical applications within these fields.展开更多
A combined experimental and numerical research study is conducted to investigate the complex relationship between the structure and the aerodynamic performances of an Archimedes spiral wind turbine(ASWT).Two ASWTs are...A combined experimental and numerical research study is conducted to investigate the complex relationship between the structure and the aerodynamic performances of an Archimedes spiral wind turbine(ASWT).Two ASWTs are considered,a prototypical version and an improved version.It is shown that the latter achieves the best aerodynamic performance when the spread angles at the three sets of blades areα_(1)=30°,α_(2)=55°,α3=60°,respectively and the blade thickness is 4 mm.For a velocity V=10 m/s,a tip speed ratio(TSR)=1.58 and 2,the maximum CP values are 0.223 and 0.263 for the prototypical ASWT and improved ASWT,respectively,and the maximum C_(P) enhancement is 17.93%.For V=10 m/s and TSR=2,the CP values of the prototypical ASWT and improved ASWT are 0.225 and 0.263,respectively,with an aerodynamic performance enhancement of 16.88%.Through mutual verification of the test outcomes and numerical results,it is concluded that the proposed approach can effectively lead to aerodynamic performance improvement.展开更多
In this paper, cold simulation experiments and numerical calculations are conducted to predict 3 D flow field aerodynamics for an oil furnace after being retrofitted due to its fuel variation. K ε model and SIMPLE ...In this paper, cold simulation experiments and numerical calculations are conducted to predict 3 D flow field aerodynamics for an oil furnace after being retrofitted due to its fuel variation. K ε model and SIMPLE program under body fit coordination (BFC) system, in which TTM non orthogonal method is used to control the irregular geometric boundary, are adopted to solve the control equations. Model tests are conducted to check the calculation results, showing that they are in agreement with each other. Three different alternatives with different side window locations are also calculated to optimize the designs. The field retrofitting results show that the combination of cold tests with numerical calculations has prosperous application in retrofitting or renewing medium and small boilers.展开更多
Natural convection heat transfer inside horizontal rectangular enclosure filled with the anisotropic porous media, with isothermally heated bottom and cooled top while the vertical walls are adiabatic, is analyzed num...Natural convection heat transfer inside horizontal rectangular enclosure filled with the anisotropic porous media, with isothermally heated bottom and cooled top while the vertical walls are adiabatic, is analyzed numerically by applying the Brinkman model-a modified form of Darcy model giving consideration to the viscous effect. The results show that: (1)the permeability ratio (K*=Ky/Kx) is an important factor affecting natural convection heat transfer in the porous media. As K' decreases, the circulation intensity of the natural convectioncells increase significantly, resulting in an enhancement of heat transfer coefficient; (2)the increase of Darcy number (aa=Ky/H2) implies that the viscous effect is more significant. As Da≥10-, there exists a certain difference between the Darcy model and the Brinkman model. It is more significant at a lower permeability ratio. In particalar, with K*≤0. 25, the Nusselt number for Da=10-3 would differ form that of Darcy model up to an amount of 30K. The Darcy flow as depicted by Darcy model is no longer existing and an analysis neglecting the viscous effect will inevitably be of considerable error.展开更多
Foundation settlement is of great significance for high-fill engineering in collapsible loess areas.To predict the construction settlement of Lüliang Airport located in Shanxi Province,China,a plane strain finite...Foundation settlement is of great significance for high-fill engineering in collapsible loess areas.To predict the construction settlement of Lüliang Airport located in Shanxi Province,China,a plane strain finite element method considering the linear variation in the modulus,was carried out in this paper based on the results of geotechnical tests.The stress and deformation of four typical sections caused by layered fill are simulated,and then the settlement of the high-fill airport is calculated and analyzed by inputting three sets of parameters.The relative soft parameters of loess geomaterials produce more settlement than the relatively hard parameters.The thicker the filling body is,the greater the settlement is.The filling body constrained by mountains on both sides produces less settlement than the filling body constrained by a mountain on only one side even the filling thickness is almost the same.The settlement caused by the original subbase accounts for 56%−77%of the total settlement,while the fill soils themselves accounts for 23%−44%of the total settlement,which is approximately consistent with the field monitoring results.It provides a good reference for predicting the settlement of similar high-fill engineering.展开更多
Numerical investigations on the flow field in Ti-Al melt during rectangular cold crucible directional solidification were carried out. Combined with the experimental results, 3-D finite element models for calculating ...Numerical investigations on the flow field in Ti-Al melt during rectangular cold crucible directional solidification were carried out. Combined with the experimental results, 3-D finite element models for calculating flow field inside melting pool were established, the characteristics of the flow under different power parameters were further studied. Numerical calculation results show that there is a complex circular flow in the melt, a rapid horizontal flow exists on the solid/liquid interface and those flows confluence in the center of the melting pool. The flow velocity v increases with the increase of current intensity, but the flow patterns remain unchanged. When the current is 1000 A, the vmax reaches 4 mm/s and the flow on the interface achieves 3 mm/s. Flow patterns are quite different when the frequency changes from 10 kHz to 100 kHz, the mechanism of the frequency influence on the flow pattern is analyzed, and there is an optimum frequency for cold crucible directional solidification.展开更多
Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties o...Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties of ship seismic wave, the numerical calculation of synthetic seismograms on seafloor aroused by a low frequency point sound source is carried out using a wave number integration technique combined with inverse Fourier transform. According to the numerical example of hard seafloor, the time series of seismic wave on seafloor are mostly composed of interface waves and normal mode waves. Each normal mode wave has a well defined low cut-off frequency, while the interface wave doesn't have. The frequency dispersion of normal mode wave is obvious when frequency is lower than 100Hz, while the interface wave is dispersive only in the infra-sound frequency range. The time series of seismic wave is dominated by the interface wave when the source frequency is less than the minimal cut-off frequency of normal mode wave.展开更多
The response characteristics of the warhead under thermal stimuli conditions are important to the safety improvement.The goal of this study is to obtain data on the warhead in the fast cook-off process.In this paper,a...The response characteristics of the warhead under thermal stimuli conditions are important to the safety improvement.The goal of this study is to obtain data on the warhead in the fast cook-off process.In this paper,a numerical calculation method is proposed,whose reliability is supported by comparison with experimental results.Through the numerical calculation,the temperature distribution,temperature change,and ignition time are acquired.The numerical results show that the ignition time is 76 s after the warhead started to burn and that the maximum temperature of the explosive’s outer surface is 238.3℃ at the ignition time.The fast cook-off experiment of the warhead is implemented so as to get the flame temperature and reaction grades that are not available through numerical calculation.The experimental results show that the overpressure fails to reach the preset minimumvalue which is equivalent to 6 kg of TNT and that the reaction grade is deflagration.The research results have reference value for the design of the warhead and the reduction of detonation risks.展开更多
Four ships,a twin-propeller naval ship,two single-propeller container ships,and a single-propeller very large crude carrier(VLCC),were studied to investigate the scale effect of the form factor.The viscous flow fields...Four ships,a twin-propeller naval ship,two single-propeller container ships,and a single-propeller very large crude carrier(VLCC),were studied to investigate the scale effect of the form factor.The viscous flow fields of the ships at different scales were solved numerically via the Reynolds-averaged Navier–Stokes method combined with the shear stress transport k–ωturbulence model.The numerical method was validated through comparisons with experimental data,and numerical uncertainty analysis was carried out based on the ITTC recommended procedure.On this basis,scale effects of the form factor were analyzed using different friction lines,and scale effects of flow fields and the mean axial wake fractions were further analyzed in details.The results showed that the form factor exhibited scale effects when adopting the ITTC-1957 line,and it increased with the increase in the Reynolds number.The scale effect of the form factor reduces the prediction precision of the full-scale ship resistance.The friction line has a significant effect on the form factor.The form factor exhibits little dependence on the Reynolds number when using the numerical friction line or the Katsui line,which is useful for full-scale ship resistance predictions.With the increasing Reynolds number,the boundary layer thickness becomes thinner and the axial velocity contour contracts toward the center plane,and there is nearly a linear relationship between the reciprocal of mean axial wake fraction on propeller disc and Reynolds number in logarithmic scale for the three types of ship forms.展开更多
The 3-D numerical computation of the flow and temperature fields for jet array impingement with initial crossflow investigates the effects of the jet-to-surface spacing, the impinging hole arrangement and the jet-to-c...The 3-D numerical computation of the flow and temperature fields for jet array impingement with initial crossflow investigates the effects of the jet-to-surface spacing, the impinging hole arrangement and the jet-to-crossflow mass flux ratio on heat transfer characteristics. The study shows that: (1) under the different jet-to-surface spacing, the impingement cooling with inline arrangement is better than that with staggered arrangement for a given jet-to-crossflow mass flux ratio;( 2 ) the value of jet-to-surface spacing impacts a complicated effect on the flow and heat transfer for jet array impingement; (3) as the ratio of crossflow-to-jet mass flux ratio increases, the cooling effectiveness on monotonous decrease for both inline and staggered arrangements at the same jet-to-surface spacing.展开更多
On the basis of a three-dimensional weakly nonliear theory of Lagrangian residual current in the Baroclinic shallow seas, a diagnostic numerical calculation of wind-driven, thermohaline and tide-induced Lagrangian res...On the basis of a three-dimensional weakly nonliear theory of Lagrangian residual current in the Baroclinic shallow seas, a diagnostic numerical calculation of wind-driven, thermohaline and tide-induced Lagrangian residual current in the Bohai Sea is made. The model involves the Richardson number in the eddy viscosity coefficient, wind, thcrmolialine and tidal effects in the focing terms. The runoff of the Huanghe River and a part of the Huanghai Warm Water coming from the Huanghai Sea through the Bohai Sea Strait is also considered. The velocity-splitting method is adopted. The wind-driven circu lation, thermohaline circulation and the tide-induced Lagrangian residual circulation are also obtained individually and analysed. The dynamics of the three main eddies in the Lagrangian mean circulation is discussed. Finally, the numerical result is partly verified with the observed data.展开更多
Based on the high order nonlinear and dispersive wave equation with a dissipative term, a numerical model for nonlinear waves is developed, It is suitable to calculate wave propagation in water areas with an arbitrari...Based on the high order nonlinear and dispersive wave equation with a dissipative term, a numerical model for nonlinear waves is developed, It is suitable to calculate wave propagation in water areas with an arbitrarily varying bottom slope and a relative depth h/L(0)less than or equal to1. By the application of the completely implicit stagger grid and central difference algorithm, discrete governing equations are obtained. Although the central difference algorithm of second-order accuracy both in time and space domains is used to yield the difference equations, the order of truncation error in the difference equation is the same as that of the third-order derivatives of the Boussinesq equation. In this paper, the correction to the first-order derivative is made, and the accuracy of the difference equation is improved. The verifications of accuracy show that the results of the numerical model are in good agreement with those of analytical Solutions and physical models.展开更多
Atomizing rainfall caused by flood discharge of high dams poses a great threat to the safety of powerhouse and ecological environment.As an indispensable means,numerical calculation is widely used in the safety design...Atomizing rainfall caused by flood discharge of high dams poses a great threat to the safety of powerhouse and ecological environment.As an indispensable means,numerical calculation is widely used in the safety design of discharge structures.The distribution of rainfall intensity is closely related to the trajectory nappe shape,jet trajectory distances,the splashed water droplet diameter and its velocity,and the spatial distribution of downstream nappe wind.In this paper,an experimental result is used to verify the improved stochastic splash mathematical model under different bucket types and discharge conditions,and the sensitivity of downstream rainfall intensity distribution to the shape of trajectory nappe,discharge flow,spatial distribution of downstream nappe wind,and the corresponding relationship between the droplet diameter and its splashing velocity is analyzed.The results show that the calculation accuracy of downstream rainfall intensity distribution is significantly improved when the above factors are taken into consideration.It is found that the bucket type and flood discharge rate play the greatest role in the rainfall intensity distribution,followed by the downstream nappe wind distribution,and finally the corresponding relationship between the diameter and velocity of splash droplets.Therefore,these factors should be considered comprehensively when the rainfall intensity distribution of flood discharge atomization is calculated.This study can help us to understand the influence factors of flood discharge atomization more deeply and predict the distribution of flood discharge atomization rainfall intensity more accurately.展开更多
Prediction of channel dredging volume is critical for project cost estimation. However, many proposed approximate methods are not accurate. This paper presents a novel numerical method to accurately calculate the dred...Prediction of channel dredging volume is critical for project cost estimation. However, many proposed approximate methods are not accurate. This paper presents a novel numerical method to accurately calculate the dredg- ing volume using a 3D stratum model (DSM) and a channel surface model. First, the 3D DSM is constructed rapidly yet accurately from non-uniform rational B-splines (NURBS) surfaces through Boolean operation between a physical terrain model and a stratum surfaces model. Then, a parametric channel surface model is built from cross-section data and a channel center line using code implemented in the VC++ programming language. Finally, the volumes of different types of physical stratums can be calculated automatically and hierarchically to determine the dredging volume. Practical application shows that the DSM method is more precise and faster compared to the section method, and that the implementation of the developed software provides an interactive graphical user interface and visual presentation.展开更多
Impulse thruster is a kind of actuator used for trajectory correction or attitude control of some in-flight munitions and vehicles.A simple mathematical model was set up to model interior characteristics of an impulse...Impulse thruster is a kind of actuator used for trajectory correction or attitude control of some in-flight munitions and vehicles.A simple mathematical model was set up to model interior characteristics of an impulse thruster.With this model,effects of some key parameters on interior performances of the impulse thruster were studied.Results show that action time is affected significantly by nozzle throat diameter and volume of combustion chamber,while output impulse is sensitive to charge particle diameter and nozzle throat diameter.Through numerical calculation,the ranges of the optimized values for some key parameters were obtained.展开更多
基金National Key R&D Program of China(Grant No.2020YFC1512404).
文摘To investigate the influence of structural parameters on the performances and internal flow characteristics of partial flow pumps at a low specific speed of 10000 rpm,special attention was paid to the first and second stage impeller guide vanes.Moreover,the impeller blade outlet width,impeller inlet diameter,blade inclination angle,and number of blades were considered for orthogonal tests.Accordingly,nine groups of design solutions were formed,and then used as a basis for the execution of numerical simulations(CFD)aimed at obtaining the efficiency values and heads for each design solution group.The influence of impeller geometric parameters on the efficiency and head was explored,and the“weight”of each factor was obtained via a range analysis.Optimal structural parameters were finally chosen on the basis of the numerical simulation results,and the performances of the optimized model were verified accordingly(yet by means of CFD).Evidence is provided that the increase in the efficiency and head of the optimized model was 12.11%and 23.5 m,respectively,compared with those of the original model.
基金the National Natural Science Foundation of China(Research Project No.52169018).
文摘A self-priming pump is a centrifugal pump that has the ability to prime itself. Typically, its performance dependson the configuration of its reflux hole. In this study, the ANSYS FLUENT software is used to investigate the effectsof three different radial positions of the reflux hole on gas-liquid two-phase distribution, pressure pulsation, andimp during self-priming. The research results indicate that: (1) The effective channel size for the reflux liquid toenter the volute varies depending on the location of the reflux hole. The effect of the impeller rotation on thereflux liquid becomes more obvious as the setting distance of the reflux aperture decreases. (2) The position ofthe reflux hole significantly affects the gas phase mass fraction inside the impeller, resulting in a significant reductionin the time it takes for the mass fraction to exceed 80%. (3) The position of the reflux hole significantly affectsthe average pressure on each monitoring surface. (4) Placing the reflux hole at a excessively distant radial distancecan result in an excessive vertical component. (5) The self-priming performance of the pump can be improved tosome extent by placing the return hole at a small radial distance.
基金supported by National Natural Science Foundation of China(No.52107142)。
文摘Plasma jet triggered gas gap switch has obvious advantages in fast control switch.The development of the plasma in the ambient medium is the key factor affecting the triggering conduction of the gas switch.However,the plasma jet process and its characteristic parameters are complicated and the existing test methods cannot fully characterize its development laws.In this work,a two-dimensional transient fluid calculation model of the plasma jet process of the gas gap switch is established based on the renormalization-group k-εturbulence equation.The results show that the characteristic parameters and morphological evolution of the plasma jet are basically consistent with the experimental results,which verifies the accuracy of the simulation model calculation.The plasma jet is a long strip with an initial velocity of 1.0 km·s-1and develops in both axial and radial directions.The jet velocity fluctuates significantly with axial height.As the plasma jet enters the main gap,the pressure inside the trigger cavity drops by80%,resulting in a rapid drop in the jet velocity.When the plasma jet head interacts with the atmosphere,the two-phase fluid compresses each other,generating a forward-propelled pressure wave.The plasma jet heads flow at high velocity,a negative pressure zone is formed in the middle part of the jet,and the pressure peak decreases gradually with height.As the value of the inlet pressure increases,the characteristic parameters of the plasma jet increase.The entrainment phenomenon is evident,which leads to an increase in the pressure imbalance of the atmospheric gas medium,leading to a significant Coanda effect.Compared with air,the characteristic parameters of a plasma jet in SF6are lower,and the morphological evolution is significantly suppressed.The results of this study can provide some insight into the mechanism of action of the switch jet plasma development process.
基金supported by the National Natural Science Foundation of China (Grant No.40874052)the Key Laboratory of Geo-detection (China University of Geosciences,Beijing),Ministry of Education
文摘Up to now, the primary method for studying critical porosity and porous media are experimental measurements and data analysis. There are few references on how to numerically calculate porosity at the critical point, pore fluid-related parameters, or framework-related parameters. So in this article, we provide a method for calculating these elastic parameters and use this method to analyze gas-bearing samples. We first derive three linear equations for numerical calculations. They are the equation of density p versus porosity Ф, density times the square of compressional wave velocity p Vp^2 versus porosity, and density times the square of shear wave velocity pVs^2 versus porosity. Here porosity is viewed as an independent variable and the other parameters are dependent variables. We elaborate on the calculation steps and provide some notes. Then we use our method to analyze gas-bearing sandstone samples. In the calculations, density and P- and S-velocities are input data and we calculate eleven relative parameters for porous fluid, framework, and critical point. In the end, by comparing our results with experiment measurements, we prove the viability of the method.
文摘Based on the diurnal consecutively observed data in the offshore area of Jiaonan in 2005, the paper tries to make a preliminary analysis of the specificity of ocean currents, tidal current property and residual current property in the area in observing dates. Then on the basis of observed data analysis and by employing the split-step method, the paper conducts a numerical simulation of the tidal current field, which can show the M2 tidal constituent tidal wave system, current ellipse distribution, maximum current velocity distribution and time-dependent current field. The calculated results agree well with the observed data, which can on the one hand reflect the basic specificities of temporal and spatial distribution of the M2 tidal constituent current field to some extent, and, on the other hand, offer more information about the hydrodynamic condition. So the paper would provide a scientific basis for the making of sea environment protection plans in the offshore area of Jiaonan under certain conditions.
基金the National Natural Science Foundation of China for financial support to this work under Grant NSFC No.12072064.
文摘The collocation method is a widely used numerical method for science and engineering problems governed by partial differential equations.This paper provides a comprehensive review of collocation methods and their applications,focused on elasticity,heat conduction,electromagnetic field analysis,and fluid dynamics.The merits of the collocation method can be attributed to the need for element mesh,simple implementation,high computational efficiency,and ease in handling irregular domain problems since the collocation method is a type of node-based numerical method.Beginning with the fundamental principles of the collocation method,the discretization process in the continuous domain is elucidated,and how the collocation method approximation solutions for solving differential equations are explained.Delving into the historical development of the collocation methods,their earliest applications and key milestones are traced,thereby demonstrating their evolution within the realm of numerical computation.The mathematical foundations of collocation methods,encompassing the selection of interpolation functions,definition of weighting functions,and derivation of integration rules,are examined in detail,emphasizing their significance in comprehending the method’s effectiveness and stability.At last,the practical application of the collocation methods in engineering contexts is emphasized,including heat conduction simulations,electromagnetic coupled field analysis,and fluid dynamics simulations.These specific case studies can underscore collocation method’s broad applicability and effectiveness in addressing complex engineering challenges.In conclusion,this paper puts forward the future development trend of the collocation method through rigorous analysis and discussion,thereby facilitating further advancements in research and practical applications within these fields.
基金supported by the National Natural Science Foundation of China.Project under Grant(Nos.51966018 and 51466015).
文摘A combined experimental and numerical research study is conducted to investigate the complex relationship between the structure and the aerodynamic performances of an Archimedes spiral wind turbine(ASWT).Two ASWTs are considered,a prototypical version and an improved version.It is shown that the latter achieves the best aerodynamic performance when the spread angles at the three sets of blades areα_(1)=30°,α_(2)=55°,α3=60°,respectively and the blade thickness is 4 mm.For a velocity V=10 m/s,a tip speed ratio(TSR)=1.58 and 2,the maximum CP values are 0.223 and 0.263 for the prototypical ASWT and improved ASWT,respectively,and the maximum C_(P) enhancement is 17.93%.For V=10 m/s and TSR=2,the CP values of the prototypical ASWT and improved ASWT are 0.225 and 0.263,respectively,with an aerodynamic performance enhancement of 16.88%.Through mutual verification of the test outcomes and numerical results,it is concluded that the proposed approach can effectively lead to aerodynamic performance improvement.
文摘In this paper, cold simulation experiments and numerical calculations are conducted to predict 3 D flow field aerodynamics for an oil furnace after being retrofitted due to its fuel variation. K ε model and SIMPLE program under body fit coordination (BFC) system, in which TTM non orthogonal method is used to control the irregular geometric boundary, are adopted to solve the control equations. Model tests are conducted to check the calculation results, showing that they are in agreement with each other. Three different alternatives with different side window locations are also calculated to optimize the designs. The field retrofitting results show that the combination of cold tests with numerical calculations has prosperous application in retrofitting or renewing medium and small boilers.
文摘Natural convection heat transfer inside horizontal rectangular enclosure filled with the anisotropic porous media, with isothermally heated bottom and cooled top while the vertical walls are adiabatic, is analyzed numerically by applying the Brinkman model-a modified form of Darcy model giving consideration to the viscous effect. The results show that: (1)the permeability ratio (K*=Ky/Kx) is an important factor affecting natural convection heat transfer in the porous media. As K' decreases, the circulation intensity of the natural convectioncells increase significantly, resulting in an enhancement of heat transfer coefficient; (2)the increase of Darcy number (aa=Ky/H2) implies that the viscous effect is more significant. As Da≥10-, there exists a certain difference between the Darcy model and the Brinkman model. It is more significant at a lower permeability ratio. In particalar, with K*≤0. 25, the Nusselt number for Da=10-3 would differ form that of Darcy model up to an amount of 30K. The Darcy flow as depicted by Darcy model is no longer existing and an analysis neglecting the viscous effect will inevitably be of considerable error.
基金Project(2020M670604)supported by the China Postdoctoral Science FoundationProject(41790434)supported by the National Natural Science Foundation of China。
文摘Foundation settlement is of great significance for high-fill engineering in collapsible loess areas.To predict the construction settlement of Lüliang Airport located in Shanxi Province,China,a plane strain finite element method considering the linear variation in the modulus,was carried out in this paper based on the results of geotechnical tests.The stress and deformation of four typical sections caused by layered fill are simulated,and then the settlement of the high-fill airport is calculated and analyzed by inputting three sets of parameters.The relative soft parameters of loess geomaterials produce more settlement than the relatively hard parameters.The thicker the filling body is,the greater the settlement is.The filling body constrained by mountains on both sides produces less settlement than the filling body constrained by a mountain on only one side even the filling thickness is almost the same.The settlement caused by the original subbase accounts for 56%−77%of the total settlement,while the fill soils themselves accounts for 23%−44%of the total settlement,which is approximately consistent with the field monitoring results.It provides a good reference for predicting the settlement of similar high-fill engineering.
基金Project (2011CB605504) supported by the National Basic Research Program of China
文摘Numerical investigations on the flow field in Ti-Al melt during rectangular cold crucible directional solidification were carried out. Combined with the experimental results, 3-D finite element models for calculating flow field inside melting pool were established, the characteristics of the flow under different power parameters were further studied. Numerical calculation results show that there is a complex circular flow in the melt, a rapid horizontal flow exists on the solid/liquid interface and those flows confluence in the center of the melting pool. The flow velocity v increases with the increase of current intensity, but the flow patterns remain unchanged. When the current is 1000 A, the vmax reaches 4 mm/s and the flow on the interface achieves 3 mm/s. Flow patterns are quite different when the frequency changes from 10 kHz to 100 kHz, the mechanism of the frequency influence on the flow pattern is analyzed, and there is an optimum frequency for cold crucible directional solidification.
基金Sponsored by National Nature Science Foundation of China ( 51179195)National Defense Foundation of China ( 513030203-02)
文摘Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties of ship seismic wave, the numerical calculation of synthetic seismograms on seafloor aroused by a low frequency point sound source is carried out using a wave number integration technique combined with inverse Fourier transform. According to the numerical example of hard seafloor, the time series of seismic wave on seafloor are mostly composed of interface waves and normal mode waves. Each normal mode wave has a well defined low cut-off frequency, while the interface wave doesn't have. The frequency dispersion of normal mode wave is obvious when frequency is lower than 100Hz, while the interface wave is dispersive only in the infra-sound frequency range. The time series of seismic wave is dominated by the interface wave when the source frequency is less than the minimal cut-off frequency of normal mode wave.
文摘The response characteristics of the warhead under thermal stimuli conditions are important to the safety improvement.The goal of this study is to obtain data on the warhead in the fast cook-off process.In this paper,a numerical calculation method is proposed,whose reliability is supported by comparison with experimental results.Through the numerical calculation,the temperature distribution,temperature change,and ignition time are acquired.The numerical results show that the ignition time is 76 s after the warhead started to burn and that the maximum temperature of the explosive’s outer surface is 238.3℃ at the ignition time.The fast cook-off experiment of the warhead is implemented so as to get the flame temperature and reaction grades that are not available through numerical calculation.The experimental results show that the overpressure fails to reach the preset minimumvalue which is equivalent to 6 kg of TNT and that the reaction grade is deflagration.The research results have reference value for the design of the warhead and the reduction of detonation risks.
基金financially supported by the National Natural Science Foundation of China(Grant No.51809273)。
文摘Four ships,a twin-propeller naval ship,two single-propeller container ships,and a single-propeller very large crude carrier(VLCC),were studied to investigate the scale effect of the form factor.The viscous flow fields of the ships at different scales were solved numerically via the Reynolds-averaged Navier–Stokes method combined with the shear stress transport k–ωturbulence model.The numerical method was validated through comparisons with experimental data,and numerical uncertainty analysis was carried out based on the ITTC recommended procedure.On this basis,scale effects of the form factor were analyzed using different friction lines,and scale effects of flow fields and the mean axial wake fractions were further analyzed in details.The results showed that the form factor exhibited scale effects when adopting the ITTC-1957 line,and it increased with the increase in the Reynolds number.The scale effect of the form factor reduces the prediction precision of the full-scale ship resistance.The friction line has a significant effect on the form factor.The form factor exhibits little dependence on the Reynolds number when using the numerical friction line or the Katsui line,which is useful for full-scale ship resistance predictions.With the increasing Reynolds number,the boundary layer thickness becomes thinner and the axial velocity contour contracts toward the center plane,and there is nearly a linear relationship between the reciprocal of mean axial wake fraction on propeller disc and Reynolds number in logarithmic scale for the three types of ship forms.
文摘The 3-D numerical computation of the flow and temperature fields for jet array impingement with initial crossflow investigates the effects of the jet-to-surface spacing, the impinging hole arrangement and the jet-to-crossflow mass flux ratio on heat transfer characteristics. The study shows that: (1) under the different jet-to-surface spacing, the impingement cooling with inline arrangement is better than that with staggered arrangement for a given jet-to-crossflow mass flux ratio;( 2 ) the value of jet-to-surface spacing impacts a complicated effect on the flow and heat transfer for jet array impingement; (3) as the ratio of crossflow-to-jet mass flux ratio increases, the cooling effectiveness on monotonous decrease for both inline and staggered arrangements at the same jet-to-surface spacing.
基金Project supported by the National Natural Science Foundation of China
文摘On the basis of a three-dimensional weakly nonliear theory of Lagrangian residual current in the Baroclinic shallow seas, a diagnostic numerical calculation of wind-driven, thermohaline and tide-induced Lagrangian residual current in the Bohai Sea is made. The model involves the Richardson number in the eddy viscosity coefficient, wind, thcrmolialine and tidal effects in the focing terms. The runoff of the Huanghe River and a part of the Huanghai Warm Water coming from the Huanghai Sea through the Bohai Sea Strait is also considered. The velocity-splitting method is adopted. The wind-driven circu lation, thermohaline circulation and the tide-induced Lagrangian residual circulation are also obtained individually and analysed. The dynamics of the three main eddies in the Lagrangian mean circulation is discussed. Finally, the numerical result is partly verified with the observed data.
基金This research was financially supported by China National Key Basic Research Project "Circulation Principal and Mathematic Model" (Grant No. 1999043810) Guangdong Science and Technology Innovation Project: "Disaster Diagnoses of Sea Walls" (99B07102G)
文摘Based on the high order nonlinear and dispersive wave equation with a dissipative term, a numerical model for nonlinear waves is developed, It is suitable to calculate wave propagation in water areas with an arbitrarily varying bottom slope and a relative depth h/L(0)less than or equal to1. By the application of the completely implicit stagger grid and central difference algorithm, discrete governing equations are obtained. Although the central difference algorithm of second-order accuracy both in time and space domains is used to yield the difference equations, the order of truncation error in the difference equation is the same as that of the third-order derivatives of the Boussinesq equation. In this paper, the correction to the first-order derivative is made, and the accuracy of the difference equation is improved. The verifications of accuracy show that the results of the numerical model are in good agreement with those of analytical Solutions and physical models.
基金supported by the National Natural Science Foundation of China(Grant Nos.U1765202,U20A20316,515779167)the Foundation for Innovative Research Groups of the Natural Science Foundation of Hebei Province of China(Grant No.E2020402074)。
文摘Atomizing rainfall caused by flood discharge of high dams poses a great threat to the safety of powerhouse and ecological environment.As an indispensable means,numerical calculation is widely used in the safety design of discharge structures.The distribution of rainfall intensity is closely related to the trajectory nappe shape,jet trajectory distances,the splashed water droplet diameter and its velocity,and the spatial distribution of downstream nappe wind.In this paper,an experimental result is used to verify the improved stochastic splash mathematical model under different bucket types and discharge conditions,and the sensitivity of downstream rainfall intensity distribution to the shape of trajectory nappe,discharge flow,spatial distribution of downstream nappe wind,and the corresponding relationship between the droplet diameter and its splashing velocity is analyzed.The results show that the calculation accuracy of downstream rainfall intensity distribution is significantly improved when the above factors are taken into consideration.It is found that the bucket type and flood discharge rate play the greatest role in the rainfall intensity distribution,followed by the downstream nappe wind distribution,and finally the corresponding relationship between the diameter and velocity of splash droplets.Therefore,these factors should be considered comprehensively when the rainfall intensity distribution of flood discharge atomization is calculated.This study can help us to understand the influence factors of flood discharge atomization more deeply and predict the distribution of flood discharge atomization rainfall intensity more accurately.
基金Supported by the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (No. 51021004)National Natural Science Foundation of China(No. 50879056)National Key Technologies R&D Program in the 12th Five-Year Plan of China(No. 2011BAB10B06)
文摘Prediction of channel dredging volume is critical for project cost estimation. However, many proposed approximate methods are not accurate. This paper presents a novel numerical method to accurately calculate the dredg- ing volume using a 3D stratum model (DSM) and a channel surface model. First, the 3D DSM is constructed rapidly yet accurately from non-uniform rational B-splines (NURBS) surfaces through Boolean operation between a physical terrain model and a stratum surfaces model. Then, a parametric channel surface model is built from cross-section data and a channel center line using code implemented in the VC++ programming language. Finally, the volumes of different types of physical stratums can be calculated automatically and hierarchically to determine the dredging volume. Practical application shows that the DSM method is more precise and faster compared to the section method, and that the implementation of the developed software provides an interactive graphical user interface and visual presentation.
基金Sponsored by the Ministerial Level Advanced Research Foundation(51305080302)
文摘Impulse thruster is a kind of actuator used for trajectory correction or attitude control of some in-flight munitions and vehicles.A simple mathematical model was set up to model interior characteristics of an impulse thruster.With this model,effects of some key parameters on interior performances of the impulse thruster were studied.Results show that action time is affected significantly by nozzle throat diameter and volume of combustion chamber,while output impulse is sensitive to charge particle diameter and nozzle throat diameter.Through numerical calculation,the ranges of the optimized values for some key parameters were obtained.