A three-dimensional numerical model is developed to study the behaviour of an argon-nitrogen plasma arc inside a non-transferred torch. In this model, both the entire cathode and anode nozzle are considered to simulat...A three-dimensional numerical model is developed to study the behaviour of an argon-nitrogen plasma arc inside a non-transferred torch. In this model, both the entire cathode and anode nozzle are considered to simulate the plasma arc. The argon-nitrogen plasma arc is simulated for different arc currents and gas flow rates of argon. Various combinations of arc core radius and arc length, which correspond to a given torch power, are predicted. A most feasible combination of the same, which corresponds to an actual physical situation of the arc inside the torch, is identified using the thermodynamic principle of minimum entropy production for a particular torch power. The effect of the arc current and gas flow rate on the plasma arc characteristics and torch efficiency is explained. The effect of the nitrogen content in the plasma gas on the torch power and efficiency is clearly detected. Predicted torch efficiencies are comparable to the measured ones and the effect of the arc current and gas flow rate on predicted and measured efficiencies is almost similar. The efficiency of the torch, cathode and anode losses and core temperature and velocity at the nozzle exit are reported for five different cases.展开更多
The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance...The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.展开更多
A three-dimensional fluid model is developed to investigate the radio-frequency inductively coupled H2 plasma in a reactor with a rectangular expansion chamber and a cylindrical driver chamber,for neutral beam injecti...A three-dimensional fluid model is developed to investigate the radio-frequency inductively coupled H2 plasma in a reactor with a rectangular expansion chamber and a cylindrical driver chamber,for neutral beam injection system in CFETR.In this model,the electron effective collision frequency and the ion mobility at high E-fields are employed,for accurate simulation of discharges at low pressures(0.3 Pa-2 Pa)and high powers(40 kW-100 kW).The results indicate that when the high E-field ion mobility is taken into account,the electron density is about four times higher than the value in the low E-field case.In addition,the influences of the magnetic field,pressure and power on the electron density and electron temperature are demonstrated.It is found that the electron density and electron temperature in the xz-plane along permanent magnet side become much more asymmetric when magnetic field enhances.However,the plasma parameters in the yz-plane without permanent magnet side are symmetric no matter the magnetic field is applied or not.Besides,the maximum of the electron density first increases and then decreases with magnetic field,while the electron temperature at the bottom of the expansion region first decreases and then almost keeps constant.As the pressure increases from 0.3 Pa to 2 Pa,the electron density becomes higher,with the maximum moving upwards to the driver region,and the symmetry of the electron temperature in the xz-plane becomes much better.As power increases,the electron density rises,whereas the spatial distribution is similar.It can be summarized that the magnetic field and gas pressure have great influence on the symmetry of the plasma parameters,while the power only has little effect.展开更多
This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare ...This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.展开更多
Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and effic...Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and efficiency of the conventional numerical modeling, I develop a new seismic modeling method by combining the FD scheme with the numerical dispersion suppression neural network(NDSNN). This method involves the following steps. First, a training data set composed of a small number of wavefield snapshots is generated. The wavefield snapshots with the low-accuracy wavefield data and the high-accuracy wavefield data are paired, and the low-accuracy wavefield snapshots involve the obvious numerical dispersion including both the temporal and spatial dispersion. Second, the NDSNN is trained until the network converges to simultaneously suppress the temporal and spatial dispersion.Third, the entire set of low-accuracy wavefield data is computed quickly using FD modeling with the large time step and the coarse grid. Fourth, the NDSNN is applied to the entire set of low-accuracy wavefield data to suppress the numerical dispersion including the temporal and spatial dispersion.Numerical modeling examples verify the effectiveness of my proposed method in improving the computational accuracy and efficiency.展开更多
Aerocapture is one of the key technologies for low-cost transportation,with high demands of autonomy,accuracy,and robustness of guidance and control,due to its high reliability requirements for only one chance of tryi...Aerocapture is one of the key technologies for low-cost transportation,with high demands of autonomy,accuracy,and robustness of guidance and control,due to its high reliability requirements for only one chance of trying.A unified numerical predictor-corrector guidance method based on characteristic models for aerocapture is proposed.The numerical predictor-corrector guidance method is used to achieve autonomy and high accuracy,and the characteristic model control method is introduced to achieve robustness.At the same time,by transforming path constraints,characteristic model equations including apogee deviation and altitude differentiation are established.Based on the characteristic model equations,a unified guidance law which can satisfy path constraints and guidance objectives simultaneously is designed.In guidance problems,guidance deviation is not directly obtained from the output of the dynamics at present,but is calculated through integral and algebraic equations.Therefore,the method of directly discretizing differential equations cannot be used to establish characteristic models,which brings great difficulty to characteristic modeling.A method for characteristic modeling of guidance problems is proposed,and convergence analysis of the proposed guidance law is also provided.Finally,a joint numerical simulation of guidance and control considering navigation deviation and various uncertainties is conducted to verify the effectiveness of the proposed method.The proposed unified method can be extended to general aerodynamic entry guidance designs,providing theoretical and methodological support for them.展开更多
Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a ...Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.展开更多
To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyur...To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.展开更多
With drilling and seismic data of Transtensional(strike-slip)Fault System in the Ziyang area of the central Sichuan Basin,SW China plane-section integrated structural interpretation,3-D fault framework model building,...With drilling and seismic data of Transtensional(strike-slip)Fault System in the Ziyang area of the central Sichuan Basin,SW China plane-section integrated structural interpretation,3-D fault framework model building,fault throw analyzing,and balanced profile restoration,it is pointed out that the transtensional fault system in the Ziyang 3-D seismic survey consists of the northeast-trending F_(I)19 and F_(I)20 fault zones dominated by extensional deformation,as well as 3 sets of northwest-trending en echelon normal faults experienced dextral shear deformation.Among them,the F_(I)19 and F_(I)20 fault zones cut through the Neoproterozoic to Lower Triassic Jialingjiang Formation,presenting a 3-D structure of an“S”-shaped ribbon.And before Permian and during the Early Triassic,the F_(I)19 and F_(I)20 fault zones underwent at least two periods of structural superimposition.Besides,the 3 sets of northwest-trending en echelon normal faults are composed of small normal faults arranged in pairs,with opposite dip directions and partially left-stepped arrangement.And before Permian,they had formed almost,restricting the eastward growth and propagation of the F_(I)19 fault zone.The F_(I)19 and F_(I)20 fault zones communicate multiple sets of source rocks and reservoirs from deep to shallow,and the timing of fault activity matches well with oil and gas generation peaks.If there were favorable Cambrian-Triassic sedimentary facies and reservoirs developing on the local anticlinal belts of both sides of the F_(I)19 and F_(I)20 fault zones,the major reservoirs in this area are expected to achieve breakthroughs in oil and gas exploration.展开更多
In an arcjet thruster,the cathode and constrictor degrade with time,and the electrical arc discharge may become unsymmetrical.In this work,a three-dimensional numerical model of a hydrogen plasma arcjet is developed a...In an arcjet thruster,the cathode and constrictor degrade with time,and the electrical arc discharge may become unsymmetrical.In this work,a three-dimensional numerical model of a hydrogen plasma arcjet is developed and validated to study the effect of unsymmetrical electric arc discharge on thruster performance.The unsymmetrical arc discharge is realized by introducing a radial shift of the cathode so that the cathode tip offset is 80μm(25%of the constrictor radius).Simulations are conducted for both axially centered cathode(coaxial)and off-centered cathode(non-coaxial)configurations with identical propellant flow rates and input current.Simulations show asymmetrical arc discharge in the non-coaxial cathode configuration,resulting in azimuthally asymmetric Joule heating,species concentrations,and velocity field.This asymmetry continues as the plasma expands in the divergent section of the nozzle.Temperature,species concentrations,and axial velocity exhibit asymmetric radial distribution at the nozzle exit.The computed Joule heating was found to reduce with cathode shift,and consequently,the thrust and specific impulse of the thruster was decreased by about 6.6%.In the case of the non-coaxial cathode,geometric asymmetry also induces a small side thrust.展开更多
Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were ex...Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.展开更多
A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-co...A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.展开更多
In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solut...In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.展开更多
A mathematical model of resin flow and temperature variation in the filling stage of the resin transfer molding (RTM) is developed based on the control volume/finite element method (CV/FEM). The effects of the heat tr...A mathematical model of resin flow and temperature variation in the filling stage of the resin transfer molding (RTM) is developed based on the control volume/finite element method (CV/FEM). The effects of the heat transfer and chemical reaction of the resin on the flow and temperature are considered. The numerical algorithm of the resin flow and temperature variation in the process of RTM are studied. Its accuracy and convergence are analyzed. The comparison of temperature variations between experimental results and model predictions is carried out for two RTM cases. Result shows that the model is efficient for evaluating the flow and temperature variation in the filling stage of RTM and there is a good coincidence between theory and experiment.展开更多
Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes ...Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.展开更多
In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, re...In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group （RNG） model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.展开更多
Coal mine longwall gateroads are subject to changing loading conditions induced by the advancing longwall face. The ground response and support requirements are closely related to the magnitude and orientation of the ...Coal mine longwall gateroads are subject to changing loading conditions induced by the advancing longwall face. The ground response and support requirements are closely related to the magnitude and orientation of the stress changes, as well as the local geology. This paper presents the monitoring results of gateroad response and support performance at two longwall mines at a 180-m and 600-m depth of cover.At the first mine, a three-entry gateroad layout was used. The second mine used a four-entry, yieldabutment-yield gateroad pillar system. Local ground deformation and support response were monitored at both sites. The monitoring period started during the development stage and continued during first panel retreat and up to second panel retreat. The two data sets were used to compare the response of the entries in two very different geotechnical settings and different gateroad layouts. The monitoring results were used to validate numerical models that simulate the loading conditions and entry response for these widely differing conditions. The validated models were used to compare the load path and ground response at the two mines. This paper demonstrates the potential for numerical models to assist mine engineers in optimizing longwall layouts and gateroad support systems.展开更多
Quantifying the rockburst consequence is of critical importance to reduce the hazards with preventative measures in underground mines and deep tunnels. Contours of energy components within a pillar model are plotted a...Quantifying the rockburst consequence is of critical importance to reduce the hazards with preventative measures in underground mines and deep tunnels. Contours of energy components within a pillar model are plotted at different rockmass damage stages, and plastic strain work and released energy are proposed as indicators of rockmass damage consequence. One pillar model under different loading stiffness is simulated to assess indicators of pillar burst and the resulting damages. The results show the rockmass damage under soft loading stiffness has larger magnitude of plastic strain work and released energy than that which is under stiff loading stiffness, indicating the rockburst consequence can be quantified with plastic strain work and released energy in numerical models. With the quantified rockburst consequence,preventative measures can be taken to avoid severe hazards to mine safety.展开更多
Mine or longwall panel layout is a 3D structure with highly non-uniform stress distribution. Recognition of such fact will facilitate underground problem identification/investigation and solving by numerical modeling ...Mine or longwall panel layout is a 3D structure with highly non-uniform stress distribution. Recognition of such fact will facilitate underground problem identification/investigation and solving by numerical modeling through proper model construction. Due to its versatility, numerical modeling is the most popular method for ground control design and problem solving. However numerical modeling results require highly experienced professionals to interpret its validity/applicability to actual mining operations due to complicated mining and geological conditions. Underground ground control monitoring is routinely performed to predict roof behavior such as weighting and weighting interval without matching observation of face mining condition while the mining pressures are being monitored, resulting in unrealistic interpretation of the obtained data on mining pressure. The importance of ground control pressure monitoring and simultaneous observation of mining and geological conditions is illustrated by an example of shield leg pressure monitoring and interpretation in an U.S. longwall coal mine: it was found that the roof strata act like a plate, not an individual block of the size of a shield dimension, as commonly assumed by all researchers and shield capacity is not a fixed property for a longwall panel or a mine or a coal seam. A new mechanism on the interaction between shield's hydraulic leg pressure and roof strata for shield loading is proposed.展开更多
A numerical model and transmission characteristic analysis of DPSK （differential phase shift keying） pressure signals in mud channels is introduced. With the control logic analysis of the rotary valve mud telemetry,...A numerical model and transmission characteristic analysis of DPSK （differential phase shift keying） pressure signals in mud channels is introduced. With the control logic analysis of the rotary valve mud telemetry, a logical control signal is built from a Gate function sequence according to the binary symbols of transmitted data and a phase-shift function is obtained by integrating the logical control signal. A mathematical model of the DPSK pressure signal is built based on principles of communications by modulating carrier phase with the phase-shift function and a numerical simulation of the pressure wave is implemented with the mathematical model by MATLAB programming. Considering drillpipe pressure and drilling fluid temperature profile along drillpipes, the drillpipe of a vertical well is divided into a number of sections. With water-based drilling fluids, the impacts of travel distance, carrier frequency, drillpipe size, and drilling fluids on the signal transmission were studied by signal transmission characteristic analysis for all the sections. Numerical calculation results indicate that the influences of the viscosity of drilling fluids and volume fraction of gas in drilling fluids on the DPSK signal transmission are more notable than the others and the signal will distort in waveform with differential attenuations of the signal frequent component.展开更多
文摘A three-dimensional numerical model is developed to study the behaviour of an argon-nitrogen plasma arc inside a non-transferred torch. In this model, both the entire cathode and anode nozzle are considered to simulate the plasma arc. The argon-nitrogen plasma arc is simulated for different arc currents and gas flow rates of argon. Various combinations of arc core radius and arc length, which correspond to a given torch power, are predicted. A most feasible combination of the same, which corresponds to an actual physical situation of the arc inside the torch, is identified using the thermodynamic principle of minimum entropy production for a particular torch power. The effect of the arc current and gas flow rate on the plasma arc characteristics and torch efficiency is explained. The effect of the nitrogen content in the plasma gas on the torch power and efficiency is clearly detected. Predicted torch efficiencies are comparable to the measured ones and the effect of the arc current and gas flow rate on predicted and measured efficiencies is almost similar. The efficiency of the torch, cathode and anode losses and core temperature and velocity at the nozzle exit are reported for five different cases.
基金Project supported by the Major Basic Research Project of National Security of China(Grant No.613157)the Excellent Young Scientists Fund of China(Grant No.51222904)
文摘The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.
基金Project supported by the National Key Research and Development Program of China(Grant No.2017YFE0300106)the National Natural Science Foundation of China(Grant No.12075049)the Fundamental Research Funds for the Central Universities,China(Grant Nos.DUT20LAB201 and DUT21LAB110).
文摘A three-dimensional fluid model is developed to investigate the radio-frequency inductively coupled H2 plasma in a reactor with a rectangular expansion chamber and a cylindrical driver chamber,for neutral beam injection system in CFETR.In this model,the electron effective collision frequency and the ion mobility at high E-fields are employed,for accurate simulation of discharges at low pressures(0.3 Pa-2 Pa)and high powers(40 kW-100 kW).The results indicate that when the high E-field ion mobility is taken into account,the electron density is about four times higher than the value in the low E-field case.In addition,the influences of the magnetic field,pressure and power on the electron density and electron temperature are demonstrated.It is found that the electron density and electron temperature in the xz-plane along permanent magnet side become much more asymmetric when magnetic field enhances.However,the plasma parameters in the yz-plane without permanent magnet side are symmetric no matter the magnetic field is applied or not.Besides,the maximum of the electron density first increases and then decreases with magnetic field,while the electron temperature at the bottom of the expansion region first decreases and then almost keeps constant.As the pressure increases from 0.3 Pa to 2 Pa,the electron density becomes higher,with the maximum moving upwards to the driver region,and the symmetry of the electron temperature in the xz-plane becomes much better.As power increases,the electron density rises,whereas the spatial distribution is similar.It can be summarized that the magnetic field and gas pressure have great influence on the symmetry of the plasma parameters,while the power only has little effect.
基金This research has been conducted under SEGTRANS project,funded by the Centre for Industrial Technological Development(CDTI,Government of Spain).
文摘This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.
基金supported by the National Natural Science Foundation of China (grant numbers: 41874160 and 92055213)。
文摘Seismic finite-difference(FD) modeling suffers from numerical dispersion including both the temporal and spatial dispersion, which can decrease the accuracy of the numerical modeling. To improve the accuracy and efficiency of the conventional numerical modeling, I develop a new seismic modeling method by combining the FD scheme with the numerical dispersion suppression neural network(NDSNN). This method involves the following steps. First, a training data set composed of a small number of wavefield snapshots is generated. The wavefield snapshots with the low-accuracy wavefield data and the high-accuracy wavefield data are paired, and the low-accuracy wavefield snapshots involve the obvious numerical dispersion including both the temporal and spatial dispersion. Second, the NDSNN is trained until the network converges to simultaneously suppress the temporal and spatial dispersion.Third, the entire set of low-accuracy wavefield data is computed quickly using FD modeling with the large time step and the coarse grid. Fourth, the NDSNN is applied to the entire set of low-accuracy wavefield data to suppress the numerical dispersion including the temporal and spatial dispersion.Numerical modeling examples verify the effectiveness of my proposed method in improving the computational accuracy and efficiency.
基金The National Key R&D Program of China(2018YFA0703800)。
文摘Aerocapture is one of the key technologies for low-cost transportation,with high demands of autonomy,accuracy,and robustness of guidance and control,due to its high reliability requirements for only one chance of trying.A unified numerical predictor-corrector guidance method based on characteristic models for aerocapture is proposed.The numerical predictor-corrector guidance method is used to achieve autonomy and high accuracy,and the characteristic model control method is introduced to achieve robustness.At the same time,by transforming path constraints,characteristic model equations including apogee deviation and altitude differentiation are established.Based on the characteristic model equations,a unified guidance law which can satisfy path constraints and guidance objectives simultaneously is designed.In guidance problems,guidance deviation is not directly obtained from the output of the dynamics at present,but is calculated through integral and algebraic equations.Therefore,the method of directly discretizing differential equations cannot be used to establish characteristic models,which brings great difficulty to characteristic modeling.A method for characteristic modeling of guidance problems is proposed,and convergence analysis of the proposed guidance law is also provided.Finally,a joint numerical simulation of guidance and control considering navigation deviation and various uncertainties is conducted to verify the effectiveness of the proposed method.The proposed unified method can be extended to general aerodynamic entry guidance designs,providing theoretical and methodological support for them.
基金Engineering and Physical Sciences Research Council (EPSRC) is also acknowledged for funding this work under Grant Number EP/N009207/1.
文摘Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.
基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200494)China Postdoctoral Science Foundation(Grant No.2021M701725)+3 种基金Jiangsu Postdoctoral Research Funding Program(Grant No.2021K522C)Fundamental Research Funds for the Central Universities(Grant No.30919011246)National Natural Science Foundation of China(Grant No.52278188)Natural Science Foundation of Jiangsu Province(Grant No.BK20211196)。
文摘To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.
基金Supported by the Key Project of National Natural Science Foundation of China(42330810).
文摘With drilling and seismic data of Transtensional(strike-slip)Fault System in the Ziyang area of the central Sichuan Basin,SW China plane-section integrated structural interpretation,3-D fault framework model building,fault throw analyzing,and balanced profile restoration,it is pointed out that the transtensional fault system in the Ziyang 3-D seismic survey consists of the northeast-trending F_(I)19 and F_(I)20 fault zones dominated by extensional deformation,as well as 3 sets of northwest-trending en echelon normal faults experienced dextral shear deformation.Among them,the F_(I)19 and F_(I)20 fault zones cut through the Neoproterozoic to Lower Triassic Jialingjiang Formation,presenting a 3-D structure of an“S”-shaped ribbon.And before Permian and during the Early Triassic,the F_(I)19 and F_(I)20 fault zones underwent at least two periods of structural superimposition.Besides,the 3 sets of northwest-trending en echelon normal faults are composed of small normal faults arranged in pairs,with opposite dip directions and partially left-stepped arrangement.And before Permian,they had formed almost,restricting the eastward growth and propagation of the F_(I)19 fault zone.The F_(I)19 and F_(I)20 fault zones communicate multiple sets of source rocks and reservoirs from deep to shallow,and the timing of fault activity matches well with oil and gas generation peaks.If there were favorable Cambrian-Triassic sedimentary facies and reservoirs developing on the local anticlinal belts of both sides of the F_(I)19 and F_(I)20 fault zones,the major reservoirs in this area are expected to achieve breakthroughs in oil and gas exploration.
基金the Indian Space Research Organization(VSSC-ISRO)for funding this research through ISRO-IITM Cell。
文摘In an arcjet thruster,the cathode and constrictor degrade with time,and the electrical arc discharge may become unsymmetrical.In this work,a three-dimensional numerical model of a hydrogen plasma arcjet is developed and validated to study the effect of unsymmetrical electric arc discharge on thruster performance.The unsymmetrical arc discharge is realized by introducing a radial shift of the cathode so that the cathode tip offset is 80μm(25%of the constrictor radius).Simulations are conducted for both axially centered cathode(coaxial)and off-centered cathode(non-coaxial)configurations with identical propellant flow rates and input current.Simulations show asymmetrical arc discharge in the non-coaxial cathode configuration,resulting in azimuthally asymmetric Joule heating,species concentrations,and velocity field.This asymmetry continues as the plasma expands in the divergent section of the nozzle.Temperature,species concentrations,and axial velocity exhibit asymmetric radial distribution at the nozzle exit.The computed Joule heating was found to reduce with cathode shift,and consequently,the thrust and specific impulse of the thruster was decreased by about 6.6%.In the case of the non-coaxial cathode,geometric asymmetry also induces a small side thrust.
基金Supported by the China National Science and Technology Major Project(2016ZX05033002,2016ZX05033001).
文摘Taking the Paleogene salt strata in the west of Kuqa foreland thrust belt as study object, the deformation features of salt structure in the compression direction and perpendicular to the compression direction were examined to find out the control factors and formation mechanisms of the salt structures. By using the three-dimensional discrete element numerical simulation method, the formation mechanisms of typical salt structures of western Kuqa foreland thrust belt in Keshen and Dabei work areas were comprehensively analyzed. The simulation results show that the salt deformation in Keshen and Dabei work areas is of forward spread type, with deformation concentrated in the piedmont zone;the salt deformation is affected by the early uplift near the compression end, pre-existing basement faults, synsedimentary process and the initial salt depocenter;in the direction perpendicular to the compression direction, salt rocks near the compression end have strong lateral mobility with the velocity component moving towards the middle part, and the closer to the middle, the larger the velocity will be, so that salt rocks will aggregate towards the middle and deform intensely, forming complex folds and separation of salt structures from salt source, and local outcrop with thrust faults. Compared with 2 D simulation, 3 D simulation can analyze salt structures in the principal stress direction and direction perpendicular to the principal stress, give us a full view of the formation mechanisms of salt structures, and guide the exploration of oil and gas reservoirs related to salt structures.
文摘A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept (EDC) model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10732010,10972010,and 11332002)
文摘In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.
文摘A mathematical model of resin flow and temperature variation in the filling stage of the resin transfer molding (RTM) is developed based on the control volume/finite element method (CV/FEM). The effects of the heat transfer and chemical reaction of the resin on the flow and temperature are considered. The numerical algorithm of the resin flow and temperature variation in the process of RTM are studied. Its accuracy and convergence are analyzed. The comparison of temperature variations between experimental results and model predictions is carried out for two RTM cases. Result shows that the model is efficient for evaluating the flow and temperature variation in the filling stage of RTM and there is a good coincidence between theory and experiment.
基金funded by a grant from Natural Sciences and Engineering Research Council of Canada (NSERC)the authors would like to acknowledge the Niobec mine (Saint-Honoré, QuébecQuébec)
文摘Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.
基金Project supported by the National Natural Science Foundation of China （Grant Nos 40233030, 40405004 and 40405014） and the Special Program of the Scientific and Social Practices for Graduate Students in Chinese Academy of Sciences, China.
文摘In this study, the Reynolds-averaged Navier-Stokes （RANS） method is employed to simulate the flow within and over an intersection model with three kinds of k-ε turbulence closure schemes, namely, standard model, renormalization group （RNG） model and realizable k-ε model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizble k-ε model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable k-ε model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.
文摘Coal mine longwall gateroads are subject to changing loading conditions induced by the advancing longwall face. The ground response and support requirements are closely related to the magnitude and orientation of the stress changes, as well as the local geology. This paper presents the monitoring results of gateroad response and support performance at two longwall mines at a 180-m and 600-m depth of cover.At the first mine, a three-entry gateroad layout was used. The second mine used a four-entry, yieldabutment-yield gateroad pillar system. Local ground deformation and support response were monitored at both sites. The monitoring period started during the development stage and continued during first panel retreat and up to second panel retreat. The two data sets were used to compare the response of the entries in two very different geotechnical settings and different gateroad layouts. The monitoring results were used to validate numerical models that simulate the loading conditions and entry response for these widely differing conditions. The validated models were used to compare the load path and ground response at the two mines. This paper demonstrates the potential for numerical models to assist mine engineers in optimizing longwall layouts and gateroad support systems.
基金funded by the National Institute of Occupational Health and Science (NIOSH) under Grant Number 200-2016-90154
文摘Quantifying the rockburst consequence is of critical importance to reduce the hazards with preventative measures in underground mines and deep tunnels. Contours of energy components within a pillar model are plotted at different rockmass damage stages, and plastic strain work and released energy are proposed as indicators of rockmass damage consequence. One pillar model under different loading stiffness is simulated to assess indicators of pillar burst and the resulting damages. The results show the rockmass damage under soft loading stiffness has larger magnitude of plastic strain work and released energy than that which is under stiff loading stiffness, indicating the rockburst consequence can be quantified with plastic strain work and released energy in numerical models. With the quantified rockburst consequence,preventative measures can be taken to avoid severe hazards to mine safety.
基金supported by the National Natural Science Foundation of China (Nos. 51604267 and 51704095)
文摘Mine or longwall panel layout is a 3D structure with highly non-uniform stress distribution. Recognition of such fact will facilitate underground problem identification/investigation and solving by numerical modeling through proper model construction. Due to its versatility, numerical modeling is the most popular method for ground control design and problem solving. However numerical modeling results require highly experienced professionals to interpret its validity/applicability to actual mining operations due to complicated mining and geological conditions. Underground ground control monitoring is routinely performed to predict roof behavior such as weighting and weighting interval without matching observation of face mining condition while the mining pressures are being monitored, resulting in unrealistic interpretation of the obtained data on mining pressure. The importance of ground control pressure monitoring and simultaneous observation of mining and geological conditions is illustrated by an example of shield leg pressure monitoring and interpretation in an U.S. longwall coal mine: it was found that the roof strata act like a plate, not an individual block of the size of a shield dimension, as commonly assumed by all researchers and shield capacity is not a fixed property for a longwall panel or a mine or a coal seam. A new mechanism on the interaction between shield's hydraulic leg pressure and roof strata for shield loading is proposed.
基金supported by High Technology Research and Development Program of China(No.2006AA06A101).
文摘A numerical model and transmission characteristic analysis of DPSK （differential phase shift keying） pressure signals in mud channels is introduced. With the control logic analysis of the rotary valve mud telemetry, a logical control signal is built from a Gate function sequence according to the binary symbols of transmitted data and a phase-shift function is obtained by integrating the logical control signal. A mathematical model of the DPSK pressure signal is built based on principles of communications by modulating carrier phase with the phase-shift function and a numerical simulation of the pressure wave is implemented with the mathematical model by MATLAB programming. Considering drillpipe pressure and drilling fluid temperature profile along drillpipes, the drillpipe of a vertical well is divided into a number of sections. With water-based drilling fluids, the impacts of travel distance, carrier frequency, drillpipe size, and drilling fluids on the signal transmission were studied by signal transmission characteristic analysis for all the sections. Numerical calculation results indicate that the influences of the viscosity of drilling fluids and volume fraction of gas in drilling fluids on the DPSK signal transmission are more notable than the others and the signal will distort in waveform with differential attenuations of the signal frequent component.
