639 results about "Numerical computing" patented technology

The embodiment of the invention discloses a numerical calculation method for water outlet icebreaking of a navigation body and related equipment, and the method comprises the steps: generating a floating ice field geometric model according to a working condition, and constructing a navigation body model of the navigation body; grid division is carried out on the floating ice field geometric model and the navigation body model, and a water outlet icebreaking calculation model is obtained; fluid-solid coupling relation definition is conducted on the effluent icebreaking calculation model, and an optimized effluent icebreaking calculation model is obtained; the optimized water outlet icebreaking calculation model is used for simulating the water outlet icebreaking conditions of the navigation body at different moments, and water outlet icebreaking numerical values of the navigation body at different moments are obtained. The fluid-solid coupling relation definition can more accurately describe the interaction between fluid and solid, and in the water outlet icebreaking process, collision and extrusion between a navigation body and an ice layer and subsequent icebreaking behaviors all relate to a complex fluid-solid coupling phenomenon, so that the fluid-solid coupling relation definition is carried out on the water outlet icebreaking calculation model, and the calculation accuracy is improved. Therefore, an accurate value of water outlet icebreaking of the navigation body is obtained.

The invention provides a numerical calculation method for double-explosion-source underwater explosion impact cavitation loads, and belongs to the technical field of computer data processing-based explosion impact wave calculation, and the method comprises the following steps: solving a propagation process of far-field underwater double-explosion-source explosion impact waves in a two-dimensional acoustic flow field by using a local discontinuous Galerkin (LDG) method; and the cavitation characteristics of the flow field under the action of the double-explosion-source underwater explosion shock waves are researched. A flow field dynamic pressure result calculated by the method is compared with a result calculated by an acoustic finite element method. The result shows that the flow field dynamic pressure result calculated by the method is well matched with the calculation result data of the acoustic finite element method, and the effectiveness and the accuracy of the method are verified.

The invention discloses a coupling iterative algorithm library and a reactor multi-physical field coupling calculation program construction method. The coupling iterative algorithm library comprises an open source library layer used for integrating a plurality of numerical calculation libraries and providing an input / output interface for a reactor multi-physics field coupling framework and each numerical calculation library; the core algorithm layer is used for integrating a plurality of coupling iteration algorithms and a plurality of acceleration algorithms related to each coupling iteration algorithm; the interface layer is used for providing a data interaction interface between a nuclear reactor multi-physics field coupling calculation program and the core algorithm layer; and the iterative algorithm template layer is used for forming an iterative solution process control template function based on the open source library layer, the core algorithm layer and the interface layer, and the iterative solution process control template function is used for constructing a nuclear reactor multi-physics field coupling calculation program based on the process control parameters set in the input file. According to the coupling iterative algorithm library, the development efficiency and the development quality of a nuclear reactor multi-physics field coupling calculation program can be improved.

The invention belongs to the technical field of weak compressible large deformation of soft materials, provides a dynamic analysis explicit substance point method for large deformation of weak compressible materials, and provides a brand new numerical calculation method for large deformation research of the weak compressible materials. According to the method, a weak-form control equation is discretized in a total Lagrangian substance point method framework, and challenges caused by grid distortion when a traditional method is used for solving a large deformation problem are avoided. Besides, the invention provides a volume self-locking processing technology based on an F-bar method, the problem of volume locking of the weakly compressible substance in numerical simulation can be simply, conveniently and effectively solved, and the volume self-locking processing technology can be conveniently embedded into an existing substance point method program. Meanwhile, a construction method of a smooth kernel interpolation function is developed, and the problem that the solving precision of a traditional substance point method in a boundary region is reduced can be effectively solved. According to the method, the problem of large deformation of the weak compressible material can be solved in a stable and efficient manner, and a new technical path is provided for research and application in related fields.

The invention relates to the field of reactor data processing, and particularly discloses a numerical calculation method for solving a neutron transport equation. The method comprises the following steps: 1, describing the geometry of a reactor core of a nuclear reactor by utilizing a construction entity geometric modeling method, and subdividing a curved surface grid along a curved surface structure; 2, directly expanding a neutron flux distribution function in each curved surface grid by using a multivariate polynomial; 3, deriving a neutron transport equation and an integral residual equation of a boundary condition under each curved surface grid; 4, solving geometric integral parameters of each curved surface grid; step 5, establishing an adaptive linear system; and step 6, iteratively calculating and solving a maximum characteristic value and an unknown expansion coefficient.

The invention relates to the field of theoretical calculation and numerical simulation of ship structure protection and cavitation, in particular to a bubble pulsation initial condition inversion method and system based on a bubble energy loss target. Inputting characteristic parameters and energy parameters of the bubbles, and determining a target variable value of inversion calculation; taking a bubble maximum volume moment as an initial calculation moment, randomly giving a bubble internal pressure value at the bubble maximum volume moment, and determining a bubble initial condition in inversion calculation; reversely carrying out time integration on the basis of a compressible spherical bubble pulsation equation to obtain bubble characteristic parameters at zero moment; calculating a time calendar evolution process of a bubble radius by taking a bubble characteristic parameter obtained by reverse time integration as a bubble pulsation initial condition, and judging whether the energy loss of the bubble reaches a target value or not; and adjusting the bubble internal pressure value at the moment of the maximum volume of the bubbles for re-calculation, solving a bubble radius time calendar process, and recording bubble characteristic parameters after reaching an expected target as a new bubble pulsation initial condition. The method is used for solving the problem that the initial bubble pulsation condition is difficult to determine in the existing bubble dynamics theory and numerical calculation.

The invention discloses a plasma sheath electromagnetic characteristic simulation method and system for a hypersonic flight vehicle, and the method comprises the steps: carrying out the fluid simulation according to the geometric shape and flight parameters of the hypersonic flight vehicle, so as to obtain the information of a flow field containing a plasma sheath region; according to the flow field information, simulation parameters of all positions of the plasma sheath space are determined, and the simulation parameters at least comprise electron density distribution and electron collision frequency distribution; based on the simulation parameters, a three-dimensional model of the plasma sheath is established, the three-dimensional model is layered according to electron density distribution, a layered three-dimensional model composed of a plurality of plasma layers is formed, and each plasma layer has the electron density within a preset range; based on the layered three-dimensional model, establishing a non-uniformly distributed electromagnetic simulation model of the plasma sheath; and carrying out grid division on the non-uniformly distributed electromagnetic simulation model, and carrying out numerical calculation by adopting a numerical solver so as to obtain the electromagnetic characteristics of the plasma sheath. According to the method, the problem of deviation between a simulation result and an actual situation of a traditional simulation model can be solved.

The invention discloses a comprehensive intelligent monitoring decision-making system and method for an overlapping shield tunnel system. The system comprises an information acquisition system, an iterative optimization state real-time evaluation system and an intelligent maintenance decision-making system. Firstly, a preset model is trained through data such as numerical calculation; operating measured data to train a feedback model, and weighting output results of the two models to obtain a tunnel deformation predicted value; then establishing an experience overlapping tunnel safety evaluation system and an iterative optimization evaluation system, performing deep association on the tunnel state, damage and maintenance decision through association degree intelligent real-time analysis, realizing tunnel state real-time monitoring, tunnel damage real-time capturing and tunnel maintenance decision real-time updating, performing intelligent matching on maintenance measures, and performing real-time monitoring on the tunnel state, the tunnel damage and the maintenance decision. And the prediction model and the safety evaluation system are iteratively optimized according to operation effect feedback. According to the invention, real-time perception, intelligent diagnosis and maintenance of the whole life cycle overall service performance of the overlapping shield system are realized, and the long-term service quality of the overlapping shield tunnel structure is improved.

The invention provides a computing power metering method and device of an intelligent computing center cloud platform based on a multi-dimensional performance index, and relates to the technical field of intelligent computing centers, intelligent computing centers, computing power infrastructures and intelligent computing clouds, and the method comprises the steps: S1, obtaining a first task type of a first computing power operation task; s2, determining a plurality of resource indexes corresponding to the first task type and a weight coefficient corresponding to each resource index; s3, collecting a parameter value corresponding to each resource index in each of the N periods; s4, calculating a computing power value corresponding to the first computing power operation task based on the weight coefficient corresponding to each resource index and the parameter value corresponding to each resource index in N periods; and S5, generating a first verification data packet based on the first task identifier of the first computing power operation task, the weight coefficient corresponding to each resource index and the parameter value corresponding to each resource index in N periods. The computing power credibility of the intelligent computing center cloud platform can be greatly improved.

The invention relates to the technical field of underground engineering safety monitoring and disaster simulation, in particular to a digital twinborn deduction platform for underground engineering disaster chain evolution simulation, which comprises a multi-disaster coupling numerical simulation module used for acquiring multi-source data including design drawing data, geological data and engineering monitoring data, and based on the multi-source data, establishing an initial stress field model by using a numerical calculation method combining a finite element and a finite difference, and carrying out leakage-settlement-structure failure multi-disaster coupling numerical simulation. According to the method, the partial differential equation for describing the evolution law of the underground engineering is used as a constraint term to be embedded into the deduction model, and online identification and dynamic extrapolation of parameters are executed by fusing real-time monitoring data, so that the accuracy of disaster evolution simulation under complex working conditions is remarkably improved; and the problem that the traditional numerical calculation method is difficult to meet the real-time requirement of digital twinning is solved.

The invention discloses a data center temperature field and flow field rapid order reduction simulation method and system, and the method comprises the following steps: S1, carrying out the CFD numerical calculation of a data center, and obtaining a CFD data set; s2, performing dimensionality reduction on the CFD data set by adopting a weighted intrinsic orthogonal decomposition (POD) method to obtain a low-dimensional mode; s3, optimizing the low-dimensional mode by adopting a machine learning model to obtain a reduced-order model; s4, performing error evaluation and physical constraint processing on the reduced-order model, and performing verification; and S5, generating a temperature cloud chart and a streamline chart by adopting a reduced-order model, and performing hotspot optimization based on a visualization result. Temperature field and flow field simulation is carried out under the guidance of the thermodynamic theory, and the actual requirement for efficient heat management of the data center is met.

The invention discloses a flow field prediction-oriented full-process automatic joint reduced-order modeling method, which comprises the following steps of: specifying a target physical field parameter space, and randomly generating a sample space according to a Latin hypercube sampling method; constructing a full-process automatic simulation tool chain, driving target physical field numerical calculation and generating a training data set; carrying out singular value decomposition-based intrinsic orthogonal decomposition on the output physical field data, and only retaining first r main feature components to construct a reduced-order data set; constructing a multi-input multi-output full-connection feedforward neural network, and modeling and training a nonlinear mapping relation between input parameters and reduced-order features; new working condition parameters are input, reduced-order features are predicted through the trained neural network, distribution of a target physical field is reconstructed according to a singular value decomposition reduction matrix, and more flexible and reliable technical support is provided for reducing the training cost of a reduced-order model and improving simulation efficiency.

The invention discloses a rock fracture prediction method and system based on a resistivity data assimilation algorithm, and the method comprises the steps: building a numerical model for simulating the asymptotic fracture process of a rock mass under the load effect according to the boundary conditions and initial conditions of the rock mass; simulating the stress damage process of the rock mass by using a numerical calculation method to obtain prediction data of stress field and fracture damage field distribution of the rock mass; monitoring the damage process of the rock mass sample in real time by using a high-density electrical method to obtain observation data; rock mass fracture damage distribution is obtained based on observation data; an assimilation algorithm is adopted, information of observation data and information of prediction data are combined, and initial conditions or parameters of the numerical model are corrected; and then corrected rock fracture prediction data is obtained. Through the data assimilation algorithm, different data and model simulation results are fused in the power frame of rock fracture, the model track is automatically adjusted continuously depending on observation, the analysis result with higher precision and more consistency is obtained, and the prediction precision and predictability of the model are improved.

The invention provides a sparse approximate inverse Shenwei parallel preprocessing method and system for electromagnetic scattering, and relates to the technical field of processor parallel computing, and the method comprises the steps: carrying out the performance hotspot analysis of a sparse approximate inverse preprocessor in an electromagnetic scattering simulation process, and determining a hotspot function of the sparse approximate inverse preprocessor; extracting a sparse matrix vector multiplication operation of a hotspot function, performing master-slave core parallel calculation on the sparse matrix vector multiplication operation, migrating the sparse matrix vector multiplication operation to slave core calculation, and accelerating hotspot calculation by using a heterogeneous parallel sparse approximate inverse algorithm; parallel computing of master-slave core intensive numerical computing tasks is achieved; in the parallel computing process, the slave cores communicate with the master core in a direct memory access mode, and during the period, the master core is in a waiting state to ensure data consistency and communication synchronization until all the slave cores compute distributed tasks.

The invention relates to the technical field of tunnel deformation prediction, in particular to a BPNN parameter inversion and LSTM tunnel deformation prediction method. On the basis of the orthogonal test design principle, multiple parameter combinations are designed for the M-C constitutive model and the HSS constitutive model, numerical calculation is carried out, and vertical displacement data of the tunnel measuring points are obtained through numerical calculation; constructing a BPNN model, and training the BPNN model by taking the vertical displacement data as input and the corresponding constitutive model parameters as output; inputting field actual measurement tunnel settlement data into the trained BPNN model, and performing inversion to obtain optimal constitutive parameters of M-C and HSS constitutive models; performing forward numerical simulation on the whole process of pipe jacking construction by adopting the optimal constitutive parameters obtained by inversion to obtain a tunnel vertical displacement evolution rule; and inputting a vertical displacement monitoring value of the measuring point of the downlink tunnel based on the long-short term memory neural network to obtain a prediction result of future settlement of the tunnel. The accuracy of the predicted value is improved.

The invention belongs to the technical field of computational electromagnetics, and discloses a numerical calculation method for simulating electromagnetic wave propagation characteristics in time-varying plasma, which is a three-dimensional Newton-HIE-FDTD method for efficiently simulating propagation of electromagnetic waves in time-varying plasma based on a hybrid explicit-implicit finite difference time domain method and Newton's second law of motion. The method inherits the advantages of the Newton-ADE-FDTD method, the time-varying electron density is accurately represented through the time derivative, and the abundant dynamic characteristics of the time-varying plasma can be accurately captured. On the basis, the method also fuses the advantages of an HIE-FDTD method, so that the time step length of the algorithm is only restricted by the sizes of grids in two directions, and compared with other simulation methods based on time-varying plasma, the method provided by the invention is more excellent in calculation efficiency and is more suitable for large-scale popularization and application. The method is especially suitable for analysis and simulation of electromagnetic wave propagation characteristics in time-varying plasma.

The invention relates to the field of multiphase numerical simulation, and discloses a large-scale particle flow numerical calculation method based on multi-ball filling, and the method comprises the steps: obtaining a fluid domain grid size, a fluid density and a large particle radius; calculating the radius of the small filling ball according to the size of the fluid domain grid; based on golden angle spiral distribution and minimum distance coverage judgment, determining the number and the positions of small filling balls, and generating a large particle model of multi-ball filling; obtaining the background flow field speed of the large particles and the filling small balls in the current time step, and calculating the background flow field speed and drag force of the large particles; updating through iterative calculation to obtain the speed and the position of the large particles and the small filling balls in the next time step; and adding drag force as a source item into a momentum equation to solve fluid motion, and obtaining a fluid computational domain grid speed in the next time step. Accurate numerical simulation of large-scale particles is realized.

The invention belongs to the technical field of nuclear reactor physical numerical calculation, and particularly relates to a method for improving precision and convergence of solving a diffusion equation based on a PINN, and the method comprises the following steps: S1, building a neural network algorithm model combining a physical information neural network PINN and source iteration; s2, aiming at a multi-region smooth transition neutron diffusion problem with small cross section difference between materials, optimizing the neural network algorithm model by adopting an optimization strategy of a training level; and S3, aiming at the multi-region neutron diffusion problem of large neutron flux gradient at the boundary due to large cross section difference between materials, improving the neural network algorithm model by adopting a domain-divided neural network architecture. According to the method, the problem that the calculation cost is exponentially increased in a high-dimensional or complex geometric scene due to the fact that a traditional finite difference method and a finite element method depend on grid division is avoided.

The invention discloses a micro-channel cooling turbine blade inner cavity anti-deposition bionic structure design method and system, and relates to the field of turbine blade cooling of aero-engines. The problems of optimizing the wall cooling effect, flow resistance and calculation efficiency of the micro-channel structure are solved. The method comprises the steps that parameterization design is conducted on a bionic turbulent flow structure, and characteristic geometric parameters are obtained; establishing a micro-channel three-dimensional model of the bionic fish scale structure; screening the characteristic geometric parameters as design variables, and designing upper and lower intervals of the variables; taking heat exchange, pressure drop and particle deposition of the micro-channel as optimization targets, adjusting a design variable interval, and constructing a test sample space; performing numerical calculation by using the channel model to obtain a multi-objective optimization parameter value; establishing an agent model between the optimization parameter values and the geometric characteristic parameters; and optimizing by utilizing an agent model and a multi-objective optimization algorithm to obtain the optimal geometric parameters of the bionic fish scale turbulent flow structure, so that the efficient design of the micro-channel structure and the efficient utilization of the cold air consumption are realized.

The present invention relates to a method and system for predicting non - linear wave loads in ship structure design. When constructing the non - linear wave load prediction model STLN, a WLIP network model based on a neural network is constructed by using a numerical calculation data set of the ship to be designed; taking the model test regular wave data set as the target domain, knowledge transfer learning of non - linear wave loads is carried out to obtain a short - term prediction model FTLN for regular waves; for the short - term prediction model FTLN of regular waves, taking the model test irregular wave data set as the target domain, knowledge is transferred again for non - linear wave loads to obtain the non - linear wave load prediction model STLN; when predicting non - linear wave loads, for an irregular wave condition of the ship, the non - linear wave load prediction model STLN is used to obtain the wave load corresponding to the irregular wave condition. The present invention can quickly and accurately predict non - linear wave loads and meet the requirements in the ship structure design stage.

The invention discloses a comprehensive pipe gallery fire positioning method based on modal recognition and a Bayesian algorithm, and the method comprises the steps: firstly, determining a comprehensive pipe gallery disaster condition for a space scene, and building a comprehensive pipe gallery fire numerical calculation model based on a computational fluid mechanics method; secondly, a physical field modal reduced-order model is constructed based on time slice data, and nonlinear modal recognition of physical information such as temperature, carbon dioxide and carbon monoxide in the utility tunnel during a fire disaster is achieved; thirdly, performing nonlinear modal feature recognition on features of physical information such as temperature, carbon dioxide and carbon monoxide by using a neural network driven by physical information; and finally, based on a Bayesian reasoning algorithm method, carrying out fire position rapid positioning. According to the method, the machine learning theory is fused, and the accuracy, interpretability and timeliness of the model are improved. The method has an important practical application value for rapid prediction and disaster source positioning of the comprehensive pipe gallery fire, and especially can rapidly respond in an emergency, thereby providing support for disaster management.

The invention provides a method and a device for dynamically metering computing power based on computing power operation tasks of an intelligent computing center cloud platform, and relates to the technical field of intelligent computing centers, intelligent computing centers, computing power infrastructures and intelligent computing cloud. The method for dynamically metering the computing power of the intelligent computing center cloud platform based on the computing power operation task comprises the following steps: S1, acquiring a parameter value of first operation data of a first GPU in a first time period, wherein the first GPU is a GPU for operating a first computing power operation task; s2, determining a first task type of the first computing power operation task based on a parameter value of the first operation data; s3, determining a weighting coefficient corresponding to the first task type; and S4, based on the weighting coefficient and the parameter value of the first operation data, calculating to obtain the calculation power value of the first calculation power operation task in the first time period. According to the method, the computing power metering accuracy of the computing power operation task in the intelligent computing center cloud platform can be greatly improved.

The invention discloses a semiconductor device critical breakdown characteristic simulation method based on a hybrid finite volume method. The method comprises the following steps: firstly, carrying out geometric modeling on a device, and carrying out mesh generation; then, the quality of the grids is judged, an FBM method is adopted for the grids with the good quality, and a CVFEM method is adopted for the grids with the poor quality to discretize a drift-diffusion equation set; and then, solving the discrete drift-diffusion equation set system by adopting a full-coupling Newton iteration method, and realizing the optimization of the boost numerical simulation process by adopting an automatic boost strategy and an initial value extrapolation method. And finally, obtaining physical field information of a power semiconductor device (LDMOS) under different external bias voltages, calculating terminal current through a post-processing method, and obtaining an I-V characteristic curve of the device after the avalanche breakdown effect is considered. The method has good numerical calculation precision, numerical stability and grid adaptability, and the complete benchmarking with the solution precision of an international advanced semiconductor simulation tool Sentaur TCAD is realized.

The invention belongs to the technical field of heat exchangers, and discloses a step-by-step iteration multi-objective optimization design method and system for a heat exchanger, numerical calculation verification is carried out on performance indexes of decision points, the maximum error of a result obtained through first optimization design is 25.3%, and the maximum error of a result obtained through second optimization design is only 0.45%. Calculation shows that the step-by-step iteration optimization strategy effectively screens out key optimization variables and gradually reduces the optimization range, and multi-parameter, wide-range, low-cost, high-precision and multi-objective optimization design is achieved. From the three performance indexes of the unit volume heat exchange amount, the average Nusselt number and the average friction coefficient, the performance of the large flow and the small fins is better; the heat exchange amount of unit volume can be improved by reducing the transverse spacing of the fins, and the friction coefficient can be reduced by increasing the transverse spacing of the fins. The method and the conclusion provide a new design calculation method and theoretical support for design and optimization of the heat exchanger.

The invention discloses a fracturing fracture extension simulation method based on a self-adaptive multi-scale phase field model, and relates to the field of oil and gas reservoir renovation, and the method is characterized in that reservoir physical property and rock mechanics parameters are obtained, and a reservoir geologic model is established according to the parameters; dividing the reservoir geometric model into a plurality of finite element grids, and assigning grid attributes of the finite element grids by using the parameters to obtain preset hydraulic fractures in the finite element grids and initial variable values on nodes of the finite element grids; according to a preset position relation between the hydraulic fracture and the finite element mesh, constructing a multi-scale adaptive mesh for dynamic expansion of the hydraulic fracture; establishing a coupling mathematical model of a stress field, a seepage field and a phase field of hydraulic fracturing, and constructing a numerical calculation expression of the coupling mathematical model according to the multi-scale adaptive grid; and based on the numerical calculation expression, solving the coupling mathematical model by taking parameters of reservoir physical properties and rock mechanics and boundary conditions of the multi-scale adaptive grid as constraints to obtain hydraulic fracturing fracture extension simulation results under different reservoir and construction conditions.

The invention discloses a seismic wave propagation numerical simulation method and a seismic wave propagation numerical simulation system suitable for a tunnel peripheral frozen soil layer, and relates to the field of seismic wave simulation propagation. According to the method, a layered model is constructed to describe physical property distribution of the frozen soil layer at different depths, and a viscoelastic model is combined to integrate the layered model into a finite element model of the tunnel peripheral frozen soil layer; dividing the finite element model by adopting an initialized uniform grid, then preliminarily simulating a seismic wave propagation process, and by comparing the spatial change rate of the wave velocity, identifying a region with a remarkably changed wave velocity gradient, and extracting the wave velocity change amplitude and directivity data to determine the boundary and shape of the region; and planning the density distribution of the updated grid in combination with the geological parameters of the gradient change region. According to the method, the boundary and internal characteristics of the wave velocity gradient region are identified, the grid shape and density are planned according to the spatial change rate and distribution characteristics of the gradient, it is ensured that the calculation region can reflect the multi-scale characteristics of the frozen soil layer, and the seismic wave propagation path and the numerical calculation of energy distribution can be optimized.

The invention provides an iterative solver optimization method and system based on a Shenwei architecture, relates to the technical field of electronic information, and aims to solve the problems that an AztecOO suitable for an SW26010Pro architecture in a new generation of Shenwei super computer at present still adopts a double-precision algorithm, the storage and calculation cost is high, and an AztecOO solver is low in calculation efficiency and low in calculation precision. The method comprises the following steps: acquiring a solving task of a sparse linear equation set, and issuing the solving task to a main core; the main core divides a solving task and allocates a memory for the solving task; the master core calls a slave core starting function and executes calculation; and the slave core feeds back a calculation result to the master core to obtain a solution of the sparse linear equation set. According to the method, the iterative solver is optimized through a designed precision selection strategy, a matrix segmentation algorithm and a mixed precision algorithm, so that the calculation speed of the algorithm is improved and memory resources are saved under the condition that sufficient calculation precision is ensured, the whole solving process of the solver is accelerated, and more efficient development efficiency is provided for researchers in the aspect of numerical calculation.

The present invention relates to the technical field of fluid mechanics. Disclosed is a variable-density low-speed flow mesoscopic numerical calculation method, comprising: determining a computational domain, and dividing or reading grid information, fluid parameters and control parameters in the computational domain; on the basis of the read information, allocating arrays, initializing corresponding parameters of the computational domain, and entering a time step loop upon completion; setting boundary conditions, and assigning all boundary macroscopic quantities in a flow field; under a finite volume framework, deriving a numerical flux of the flow field by means of a BGK equation to solve a Navier-Stokes equation of the flow field; calculating pressure and velocity fluxes at cell interfaces of the flow field, and updating boundary macroscopic quantity information on the basis of flux information; determining whether there is an external force, and if so, applying the external force by means of splitting, and updating the boundary macroscopic quantities; determining whether a time step meets a loop condition; if so, re-entering setting boundary conditions, and assigning all boundary macroscopic quantities in the flow field; and if not, outputting flow field data, and ending simulation. The present invention has the advantages of high algorithm accuracy and simple programming.

The invention relates to a numerical calculation method and device for simulating slippage of a duct piece and a built-in steel bar. The numerical calculation method for simulating slippage of the duct piece and the built-in steel bar comprises the following steps: performing grid division on a three-dimensional shield duct piece according to reinforcement information, endowing the grid three-dimensional shield duct piece with basic physical parameters of concrete, and constructing cable units arranged along actual reinforcement positions, coupling with the concrete model in combination with the calibrated bonding parameters to obtain a concrete model with slippage characteristics; a load path conforming to an engineering scene is applied to the concrete model with the slippage characteristic, nonlinear increment iteration solving is adopted, and simulation time history data is obtained; and finally, performing index calculation and evaluation on the simulation time history data. According to the numerical calculation method for simulating slippage of the duct piece and the built-in steel bar, the problem of overestimation of rigidity and bearing capacity caused by assumption of complete bonding is effectively avoided.

The invention provides a method for calculating a non-center-guide multi-arch tunnel structure under blasting vibration velocity control, and belongs to the field of tunnel engineering design and construction, and the method comprises the following steps: 1, obtaining stratum physical mechanics and blasting excavation design parameters of a designed tunnel; 2, establishing a stratum-structure model of the non-middle-guide multi-arch tunnel, performing numerical simulation calculation of tunnel blasting load action, and obtaining a blasting excavation surrounding rock plastic zone and blasting vibration velocity distribution; 3, the relation between the surrounding rock peak blasting vibration velocity and surrounding rock rigidity reduction is established, and a surrounding rock rigidity reduction coefficient is obtained; 4, based on surrounding rock rigidity reduction, a load-structure model of the non-middle-guide multi-arch tunnel is established, and structural stress safety is analyzed; and step 5, based on a numerical calculation result and a safety discriminant of the non-center-guiding multi-arch tunnel lining, evaluating the mechanical safety reliability of the tunnel structure, and feeding back and optimizing the tunnel structure design according to the mechanical safety reliability of the tunnel structure, including structure enhancement or blasting vibration velocity control target adjustment.