The resistivity difference between oil and gas layers and the water layers in low contrast tight sandstone reservoirs is subtle. Fluid identification and saturation calculation based on conventional logging methods ar...The resistivity difference between oil and gas layers and the water layers in low contrast tight sandstone reservoirs is subtle. Fluid identification and saturation calculation based on conventional logging methods are facing challenges in such reservoirs. In this paper, a new method is proposed for fluid identification and saturation calculation in low contrast tight sandstone reservoirs. First, a model for calculating apparent formation water resistivity is constructed, which takes into account the influence of shale on the resistivity calculation and avoids apparent formation water resistivity abnormal values.Based on the distribution of the apparent formation water resistivity obtained by the new model, the water spectrum is determined for fluid identification in low contrast tight sandstone reservoirs.Following this, according to the average, standard deviation, and endpoints of the water spectrum, a new four-parameter model for calculating reservoir oil and gas saturation is built. The methods proposed in this paper are applied to the low contrast tight sandstone reservoirs in the Q4 formation of the X53 block and X70 block in the south of Songliao Basin, China. The results show that the water spectrum method can effectively distinguish oil-water layers and water layers in the study area. The standard deviation of the water spectrum in the oil-water layer is generally greater than that in the water layer. The new four-parameter model yields more accurate oil and gas saturation. These findings verify the effectiveness of the proposed methods.展开更多
Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have on...Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.展开更多
The Euler-Euler model is less effective in capturing the free surface of flow film in the spiral separator,and thus a Eulerian multi-fluid volume of fluid(VOF)model was first proposed to describe the particulate flow ...The Euler-Euler model is less effective in capturing the free surface of flow film in the spiral separator,and thus a Eulerian multi-fluid volume of fluid(VOF)model was first proposed to describe the particulate flow in spiral separators.In order to improve the applicability of the model in the high solid concentration system,the Bagnold effect was incorporated into the modelling framework.The capability of the proposed model in terms of predicting the flow film shape in a LD9 spiral separator was evaluated via comparison with measured flow film thicknesses reported in literature.Results showed that sharp air–water and air-pulp interfaces can be obtained using the proposed model,and the shapes of the predicted flow films before and after particle addition were reasonably consistent with the observations reported in literature.Furthermore,the experimental and numerical simulation of the separation of quartz and hematite were performed in a laboratory-scale spiral separator.When the Bagnold lift force model was considered,predictions of the grade of iron and solid concentration by mass for different trough lengths were more consistent with experimental data.In the initial development stage,the quartz particles at the bottom of the flow layer were more possible to be lifted due to the Bagnold force.Thus,a better predicted vertical stratification between quartz and hematite particles was obtained,which provided favorable conditions for subsequent radial segregation.展开更多
This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance ...This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.展开更多
Based on analysis of NMR T2 spectral characteristics,a new method for identifying fluid properties by decomposing T2 spectrum through signal analysis has been proposed.Because T2 spectrum satisfies lognormal distribut...Based on analysis of NMR T2 spectral characteristics,a new method for identifying fluid properties by decomposing T2 spectrum through signal analysis has been proposed.Because T2 spectrum satisfies lognormal distribution on transverse relaxation time axis,the T2 spectrum can be decomposed into 2 to 5 independent component spectra by fitting the T2 spectrum with Gauss functions.By analyzing the free relaxation response characteristics of crude oil and formation water,the dynamic response characteristics of the core mutual drive between oil and water,the petrophysical significance of each component spectrum is clarified.T2 spectrum can be decomposed into clay bound water component spectrum,capillary bound fluid component spectrum,micropores fluid component spectrum and macropores fluid component spectrum.According to the nature of crude oil in the target area,the distribution range of T2 component spectral peaks of oil-bearing reservoir is 165-500 ms on T2 time axis.This range can be used to accurately identify fluid properties.This method has high adaptability in identifying complex oil and water layers in low porosity and permeability reservoirs.展开更多
Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations wr...Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations written in a quasi conservative form, which is solved by a standard high resolution piecewise parabolic method (PPM) with multi fluid Riemann solver. The level set interface capturing method uses a narrow band ghost fluid method (GFM) with no numerical smearing. Several examples are presented and compared for one and two dimensions, which show the feasibility of the two methods applied to various multi fluid problems.展开更多
The ghost fluid method and the level set function is used to simulate the Richtmyer Meshkov instability of multi material flow interface in two dimensions. To keep the density profile from smearing out, the isobaric f...The ghost fluid method and the level set function is used to simulate the Richtmyer Meshkov instability of multi material flow interface in two dimensions. To keep the density profile from smearing out, the isobaric fix technique is used in numerical algorithms.展开更多
VOF(Volume of Fluid)方法由于其良好的守恒性和网格适应性,且具有计算资源需求相对较小等优点,成为船舶水动力学领域自由面流动CFD(Computational Fluid Dynamics)模拟的主流方法。但原始的VOF方法存在较为严重的界面扩散问题,导致模...VOF(Volume of Fluid)方法由于其良好的守恒性和网格适应性,且具有计算资源需求相对较小等优点,成为船舶水动力学领域自由面流动CFD(Computational Fluid Dynamics)模拟的主流方法。但原始的VOF方法存在较为严重的界面扩散问题,导致模拟的界面厚度过大、空间分辨率不够,进而影响流场相关变量的计算精度,这一问题在非定常自由面流动模拟中尤为明显。本文针对上述问题,通过在VOF控制方程中引入人工对流项以达到抑制界面扩散、压缩界面厚度的目的,并采用隐式离散人工对流项的方式提高计算稳定性,形成了反扩散VOF算法。经Zalesak和剪切场等经典算例在不同数量网格下的测试验证,表明反扩散VOF算法能够大幅压缩界面厚度,同时明显减小质量误差。随后的三维无障碍溃坝算例和破舱进水算例,进一步证明了反扩散VOF算法在实际非定常流动模拟中能够更好地捕捉自由面并提高计算精度。展开更多
A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-unifor...A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.展开更多
The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the bea...The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the beam channel for absorbing the divergence beam during the beam transmission process in the EAST-NBI system.The gas baffle entrance collimator(GBEC) is a typical high-heat-flux component located at the entrance of gas baffle. An efficient and accurate analysis of its thermodynamic performance is of great significance to explore the working limit and to ensure safe operation of the system under a high-parameter steady-state condition. Based on the thermo-fluid coupled method, thermodynamic analysis and simulation of GBEC is performed to get the working states and corresponding operating limits at different beam extraction conditions. This study provides a theoretical guidance for the next step to achieve long pulse with highpower experimental operation and has an important reference to ensure the safe operation of the system.展开更多
In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation...In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation mathematical model of fluid-structure interaction(FSI)was developed.Then,the transfer matrix method(TMM)was used to calculate the modal frequency,modal shape and frequency response.The results were compared with that in experiment to verify the correctness of the TMM and the results show that the fluid-structure coupling has a greater impact on the modal frequencies than the modal shape.Finally,the influence on the response spectrum of different damping ratios was studied and the results show that the natural frequency under different damping ratios has changed little but there is a big difference for the pressure spectrum.With the decreasing of damping ratio,the damping of the system on frequency spectrum is more and more significant and the dispersion and dissipation is more and more apparent.Therefore the appropriate damping ratio should be selected to minimize the effects of the vibration of the FSI.The results provide references for the theory research of FSI in the transient process.展开更多
A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as th...A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.展开更多
Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain M...Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.展开更多
Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The b...Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.展开更多
A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underex...A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.展开更多
基金funded by the National Natural Science Foundation of China (42174131)。
文摘The resistivity difference between oil and gas layers and the water layers in low contrast tight sandstone reservoirs is subtle. Fluid identification and saturation calculation based on conventional logging methods are facing challenges in such reservoirs. In this paper, a new method is proposed for fluid identification and saturation calculation in low contrast tight sandstone reservoirs. First, a model for calculating apparent formation water resistivity is constructed, which takes into account the influence of shale on the resistivity calculation and avoids apparent formation water resistivity abnormal values.Based on the distribution of the apparent formation water resistivity obtained by the new model, the water spectrum is determined for fluid identification in low contrast tight sandstone reservoirs.Following this, according to the average, standard deviation, and endpoints of the water spectrum, a new four-parameter model for calculating reservoir oil and gas saturation is built. The methods proposed in this paper are applied to the low contrast tight sandstone reservoirs in the Q4 formation of the X53 block and X70 block in the south of Songliao Basin, China. The results show that the water spectrum method can effectively distinguish oil-water layers and water layers in the study area. The standard deviation of the water spectrum in the oil-water layer is generally greater than that in the water layer. The new four-parameter model yields more accurate oil and gas saturation. These findings verify the effectiveness of the proposed methods.
基金supported by National Science Foundation of China (10576015)
文摘Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.
基金the National Natural Science Foundation of China(Nos.51974065 and 52274257)the Open Foundation of State Key Laboratory of Mineral Processing(No.BGRIMMKJSKL-2020-13)the Fundamental Research Funds for the Central Universities(Nos.N2201008 and N2201004).
文摘The Euler-Euler model is less effective in capturing the free surface of flow film in the spiral separator,and thus a Eulerian multi-fluid volume of fluid(VOF)model was first proposed to describe the particulate flow in spiral separators.In order to improve the applicability of the model in the high solid concentration system,the Bagnold effect was incorporated into the modelling framework.The capability of the proposed model in terms of predicting the flow film shape in a LD9 spiral separator was evaluated via comparison with measured flow film thicknesses reported in literature.Results showed that sharp air–water and air-pulp interfaces can be obtained using the proposed model,and the shapes of the predicted flow films before and after particle addition were reasonably consistent with the observations reported in literature.Furthermore,the experimental and numerical simulation of the separation of quartz and hematite were performed in a laboratory-scale spiral separator.When the Bagnold lift force model was considered,predictions of the grade of iron and solid concentration by mass for different trough lengths were more consistent with experimental data.In the initial development stage,the quartz particles at the bottom of the flow layer were more possible to be lifted due to the Bagnold force.Thus,a better predicted vertical stratification between quartz and hematite particles was obtained,which provided favorable conditions for subsequent radial segregation.
文摘This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.
基金Supported by the China National Science and Technology Major Project(2016ZX05050)
文摘Based on analysis of NMR T2 spectral characteristics,a new method for identifying fluid properties by decomposing T2 spectrum through signal analysis has been proposed.Because T2 spectrum satisfies lognormal distribution on transverse relaxation time axis,the T2 spectrum can be decomposed into 2 to 5 independent component spectra by fitting the T2 spectrum with Gauss functions.By analyzing the free relaxation response characteristics of crude oil and formation water,the dynamic response characteristics of the core mutual drive between oil and water,the petrophysical significance of each component spectrum is clarified.T2 spectrum can be decomposed into clay bound water component spectrum,capillary bound fluid component spectrum,micropores fluid component spectrum and macropores fluid component spectrum.According to the nature of crude oil in the target area,the distribution range of T2 component spectral peaks of oil-bearing reservoir is 165-500 ms on T2 time axis.This range can be used to accurately identify fluid properties.This method has high adaptability in identifying complex oil and water layers in low porosity and permeability reservoirs.
文摘Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations written in a quasi conservative form, which is solved by a standard high resolution piecewise parabolic method (PPM) with multi fluid Riemann solver. The level set interface capturing method uses a narrow band ghost fluid method (GFM) with no numerical smearing. Several examples are presented and compared for one and two dimensions, which show the feasibility of the two methods applied to various multi fluid problems.
文摘The ghost fluid method and the level set function is used to simulate the Richtmyer Meshkov instability of multi material flow interface in two dimensions. To keep the density profile from smearing out, the isobaric fix technique is used in numerical algorithms.
文摘VOF(Volume of Fluid)方法由于其良好的守恒性和网格适应性,且具有计算资源需求相对较小等优点,成为船舶水动力学领域自由面流动CFD(Computational Fluid Dynamics)模拟的主流方法。但原始的VOF方法存在较为严重的界面扩散问题,导致模拟的界面厚度过大、空间分辨率不够,进而影响流场相关变量的计算精度,这一问题在非定常自由面流动模拟中尤为明显。本文针对上述问题,通过在VOF控制方程中引入人工对流项以达到抑制界面扩散、压缩界面厚度的目的,并采用隐式离散人工对流项的方式提高计算稳定性,形成了反扩散VOF算法。经Zalesak和剪切场等经典算例在不同数量网格下的测试验证,表明反扩散VOF算法能够大幅压缩界面厚度,同时明显减小质量误差。随后的三维无障碍溃坝算例和破舱进水算例,进一步证明了反扩散VOF算法在实际非定常流动模拟中能够更好地捕捉自由面并提高计算精度。
基金Project supported by the National Natural Science Foundation of China (Grant No. 51075311)
文摘A two-phase wedge-sliding model is developed based on the micro-cellular structure and minimum entropy theory of a stable system, and it is used to describe the ingredient distribution of a mixed fluid in a non-uniform stress field and to analyse its phase drift phenomenon. In the model, the drift-inhibition angle and the expansion-inhibition angle are also deduced and used as evaluating indexes to describe the drifting trend of different ingredients among the mixed fluids. For solving above two indexes of the model, a new calculation method is developed and used to compute the phase distributions of multiphase fluid at peak stress and gradient area stress, respectively. As an example, the flow process of grease in a pipe is analysed by simulation method and used to verify the validity of the model.
基金supported by the National Natural Science Foundation of China(No.11605234)the Foundation of ASIPP(No.DSJJ-15-GC02)
文摘The world's first full Experimental Advanced Superconducting Tokamak(EAST) is designed with the auxiliary heating method of neutral beam injection(NBI)system. Beam collimators are arranged on both sides of the beam channel for absorbing the divergence beam during the beam transmission process in the EAST-NBI system.The gas baffle entrance collimator(GBEC) is a typical high-heat-flux component located at the entrance of gas baffle. An efficient and accurate analysis of its thermodynamic performance is of great significance to explore the working limit and to ensure safe operation of the system under a high-parameter steady-state condition. Based on the thermo-fluid coupled method, thermodynamic analysis and simulation of GBEC is performed to get the working states and corresponding operating limits at different beam extraction conditions. This study provides a theoretical guidance for the next step to achieve long pulse with highpower experimental operation and has an important reference to ensure the safe operation of the system.
文摘In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation mathematical model of fluid-structure interaction(FSI)was developed.Then,the transfer matrix method(TMM)was used to calculate the modal frequency,modal shape and frequency response.The results were compared with that in experiment to verify the correctness of the TMM and the results show that the fluid-structure coupling has a greater impact on the modal frequencies than the modal shape.Finally,the influence on the response spectrum of different damping ratios was studied and the results show that the natural frequency under different damping ratios has changed little but there is a big difference for the pressure spectrum.With the decreasing of damping ratio,the damping of the system on frequency spectrum is more and more significant and the dispersion and dissipation is more and more apparent.Therefore the appropriate damping ratio should be selected to minimize the effects of the vibration of the FSI.The results provide references for the theory research of FSI in the transient process.
文摘A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.
基金the sponsorship of the National Basic Research Program of China (973 Program,2013CB228604,2014CB239201)the National Oil and Gas Major Projects of China (2011ZX05014-001-010HZ,2011ZX05014-001-006-XY570) for their funding of this research
文摘Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.
文摘Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.
文摘A discontinuous Galerkin finite element method （DG-FEM） is developed for solving the axisymmetric Euler equations based on two-dimensional conservation laws. The method is used to simulate the unsteady-state underexpanded axisymmetric jet. Several flow property distributions along the jet axis, including density, pres- sure and Mach number are obtained and the qualitative flowfield structures of interest are well captured using the proposed method, including shock waves, slipstreams, traveling vortex ring and multiple Mach disks. Two Mach disk locations agree well with computational and experimental measurement results. It indicates that the method is robust and efficient for solving the unsteady-state underexpanded axisymmetric jet.
