C程序中的内存泄漏机制分析与检测方法设计

C语言作为安全关键软件的主要实现语言,其存在的内存泄漏缺陷具有很高的隐蔽性和危害性,如何保证内存泄漏检测的准确性和高效性是一大挑战。静态分析具有直接分析源码、能够较早发现软件错误,从而降低修复代价的优势。基于静态分析技术,提出了一种基于路径敏感的值流分析的内存泄漏检测方法,首先进行指针分析生成精确指向信息;然后基于指向信息构建值流约束,执行可达性分析以识别程序中的泄漏路径;最后借助指针与内存地址的有效生命周期进行验证。在典型基准C程序上的实验结果分析表明,本文方法与现有技术相比在效率和精度上都具有一定优势。

基于摘要的内存泄露分析方法改进

为了更加精确高效地分析源代码中的内存泄露,针对内存泄露属性分析建模,通过多种途径对现有的方法进行改进。使用函数摘要的方法提高分析效率,并使用对库函数建模、检测虚假路径的方法降低误报,同时尝试检测动态内存通过参数逃逸的情况。在crystal基础上实现了静态分析内存泄露的工具SMD,取得了良好的实验结果。

反例引导的C代码空间流模型检测方法

软件验证一直是确保软件正确性和安全性的热点研究问题.然而,由于程序语言复杂的语法语义特性,应用形式化方法验证程序的正确性存在准确度低和效率差的问题.其中,由指针操作带来的地址空间的状态变化使得现有模型检测方法的检测准确度难以得到保证.为此,通过结合模型检测与稀疏值流分析方法,设计了一种空间流模型,实现了对C程序在符号变量层面和地址空间层面的状态行为的有效描述,并提出了一种反例引导的抽象细化和稀疏值流强更新算法(CEGAS),实现了C程序指向信息敏感的形式化验证.建立了包含多种指针操作的C代码基准库,并基于该基准库进行了对比实验.实验结果表明:所提出的模型检测算法CEGAS在分析含有多种C代码特性的任务中,与现有模型检测工具相比均能取得突出的结果,其检测准确度为92.9%,每行代码的平均检测时间为2.58 ms,优于现有检测工具.

Numerical modeling of the 2D time-domain transient electromagnetic secondary field of the line source of the current excitation

To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.

VIBRATION NUMERICAL ANALYSIS OF COUNTER-ROTATING TURBINE WITH WAKE-FLOW USING FLUID-STRUCTURE INTERACTION METHOD

The design of counter-rotating turbine is one of new techniques to improve the thrust-weight ratio of jet propulsion engines.Numerical analysis of a low pressure（LP）counter-rotating turbine rotor blade is presented by using ANSYS/CFX software.Interaction of aerodynamics and solid mechanics coupling in the computation is applied.In some rating of turbine,stress distribution and vibration characteristics of low pressure turbine（LPT）blade are computed.The wake aerodynamic forces and LPT blade vibration are transformed in frequency domain using fast Fourier transform（FFT）method.The results show that under wake aerodynamic force excitation,the first order modal vibration is more easily aroused and the higher order response cannot be ignored.Moreover,with different temperature fields,the vibration responses of blade are also different.

Numerical Simulation and Analysis of Gas Flow Field in Serrated Valve Column

A three-dimensional computational fluid dynamics （CFD） model for gas flow through a serrated valve tray was presented. The flow field, as well as the dry-pressure drop of the valve under the full-opening condition was simulated based on the proposed model by using FLUENT 6.0 software. Compared with the values of dry-pressure dro.p in different turbulent models, the.simulated.results using RNG κ-ε model are in reasonable agreement with experimental data, indicating that RNG κ-ε model is suitable in simulating gas flow through the serrated valve tray. Then the CFD model combining RNG κ-ε model was used to study the three-dimensional gas flow through the considered serrated valve tray. The simulated results showed that various eddies existed on the serrated valve tray, and both the eddy and the non-eddy areas were nearly equal. The existence of addendum can decrease the eddy area caused by gas passing through the lateral outlet slots. The size of eddies can be reduced by optimizing the distance between valves.

Influence of water-rich tunnel by shield tunneling on existing bridge pile foundation in layered soils

At present,shield tunneling often needs to pass through a large number of bridge pile foundations.However,there are few studies on the influence of shield tunneling on adjacent pile foundations by combining with groundwater seepage.Based on Winkler model,the calculation equations of shield tunneling on vertical and horizontal displacement of adjacent bridge pile are derived.Meanwhile,full and part three-dimensional finite element models are established to analyze the trend of bridge pier settlement,ground surface settlement trough,vertical and horizontal displacement of the pile and pile stress under three calculation conditions,i.e.,not considering groundwater effect,considering stable groundwater effect and fluid-soil interaction.The results show that the calculated value is small when the effect of groundwater is not considered;the seepage velocity of the soil above the excavation face is faster than that of the surrounding soil under fluid-soil interaction,and after the shield passing,the groundwater on both sides shows a flow trend of“U”shape on the ground surface supplying to the upper part of the tunnel;the vertical displacement of the pile body is bounded by the horizontal position of the top of the tunnel,the upper pile body settles,and the lower pile body deforms upward.The horizontal displacement of pile body presents a continuous“S”shape distribution,causing stress concentration near the tunnel.The calculated results of fluid-soil interaction are in good agreement with the field measured data and accord with the actual situation.

Experimental measurement and numerical analysis of fused taper shape for optical fiber coupler

To find out the effect of the shape of fused taper region on the optical fiber coupler, the fiber couplers were fabricated at different drawing speeds with a six-axes fiber coupler machine. The results, which were obtained fi＇om the shape of fused taper region measured with microscope, show that there is a close correlation between the cone angle and optical performance of fiber coupler. High-performance fiber coupler cannot be obtained until rheological shape is controlled accurately. The numerical analysis model, which was built based on generalized Maxwell viscoelastic theory, is resolved with ANSYS software. The calculated results accord with the experimental data. It can apply a theoretic basis for forecasting the shape of fiber coupler fabricated under the conditions of different technological parameters.

Prediction of multi-borehole undermine coalbed gas drainage

By analyzing the flow character of a single drainage borehole in its effectingtime and the correlative theory introduced,the reason for 'inflexion' appearance in theflow character curve of the single draining borehole in a multi-borehole was studied.Takingthe theory of permeation fluid mechanics and so on as basis,the coalbed gas flowmodel was set up,and the numerical simulation analyzer was built for undermine gasproducts.With the results from the analyzer,the gas capacity could be calculated underdifferent conditions and comparisons made with the site measurement data.

Numerical analysis on coal-breaking process under high pressure water jet

Based on the theory of nonlinear dynamic finite element,the control equation ofcoal and water jet was acquired in the coal breaking process under a water jet.The calculationmodel of coal breaking under a water jet was established;the fluid-structure couplingof water jet and coal was implemented by penalty function and convection calculation.The dynamic process of coal breaking under a water jet was simulated and analyzed bycombining the united fracture criteria of the maximum tensile strain and the maximal shearstrain in the two cases of damage to coal and damage failure to coal.

Rheological numerical simulation for thermo-hydro-mechanical coupling analysis for rock mass

Under the environment of seepage field, stress field and temperature field interaction and influence, the three fields will not only produce coupling effect, but also have deformation with time due to the rheological behavior of rock mass. In the paper, based on the fundamental theories of rock mass coupling theory and rheological mechanics, the rheological model for fully coupled thermo-hydro-mechanical analysis for rock mass was set up, and the corresponding constitutive relationship, the conservation equation of mass and the conservation equation of energy were given, and the finite element formulas were derived for coupling analysis of rock mass. During establishing governing equations, rock mass was assumed approximately as macro-equivalent continuum medium. The obtained rheological numerical model for fully coupled thermo-hydro-mechanical analysis can be used for analyzing and predicting the long-term stability of underground caverns and slope engineering under the condition of thermo-hydro-mechanical coupling with rheological deformation.

Computational analysis of the flow of pseudoplastic power-law fluids in a microchannel plate

The flow of pseudoplastic power-law fluids with different flow indexes at a microchannel plate was studied using computational fluid dynamic simulation.The velocity distribution along the microchannel plate and especially in the microchannel slits,flow pattern along the outlet arc and the pressure drop through the whole of microchannel plate were investigated at different power-law flow indexes.The results showed that the velocity profile in the microchannel slits for low flow index fluids was similar to the plug flow and had uniform pattern.Also the power-law fluids with lower flow indexes had lower stagnation zones near the outlet of the microchannel plate.The pressure drop through the microchannel plate showed huge differences between the fluids.The most interesting result was that the pressure drops for power-law fluids were very smaller than that of Newtonian fluids.In addition,the heat transfer of the fluids through the microchannel with different channel numbers in a wide range of Reynolds number was investigated.For power-law fluid with flow index(n=0.4),the Nusselt number increases continuously as the number of channels increases.The results highlight the potential use of using pseudoplastic fluids in the microheat exchangers which can lower the pressure drop and increase the heat transfer efficiency.

Numerical analysis of behavior of active layer in rotary kilns by discrete element method

The behavior of the active layer of material bed within rotary kilns plays a key role in industrial applications.To obtain its influences on industrial process,different regimes of particle motion have been simulated by discrete element method(DEM) in three dimensions under variant rotation speeds,filling degree,based on the background of induration process of iron ore pellets.The influences of the mentioned factors on the maximum thickness of the active layer and the average velocity of particles have been investigated.The average velocity of particles increases with Froude number following the power function over a wide range,and the maximum thickness rises with increasing rotation speed in a way of logarithm.The influence of the filling degree f on the maximum thickness exhibits a good linearity under two classic regimes,but the increasing of the average velocity of the active layer is limited at f=0.4.This basic research highlights the impact of the active layer within rotary kilns,and lays a good foundation for the further investigation in mixing and heat transfer within the particle bed inside rotary kilns.

Effects of Shape and Quantity of Helical Baffle on the Shell-side Heat Transfer and Flow Performance of Heat Exchangers

Shape and quantity of helical baffles have great impact on the shell-side performance of helical baffled heat exchangers （HBHE）. In this work, three physical models of HBHE with baffles of different shape （trisection, quadrant and sextant sector） were investigated. Numerical simulations were performed on HBHE at three helix an- gles （10°, 25° and 40°） by the software ANSYS CFX. Analyses of numerical results indicate that the sextant HBHE shows relatively better fluid flow performance because the leakage flow in the triangle area is evidently reduced and the fluid streamline appears much closer to an ideal spiral flow, while the trisection and quadrant HBHE show more scattered and disordered streamline distributions. The convective heat transfer coefficient and pressure drop in three types of HBHE were presented. Further investigations on the shell side performance with different helical baf- fles were implemented by the field synergy theory. Both theoretical and numerical analyses gave support on the re- lations between helical baffle shape and shell-side performance. This paper may provide useful reference for the selection of baffle shade and auantitv in HBHE.

Numerical Analysis of Adaptation-to-capacity Length for Fluvial Sediment Transport

Over the last several decades,various sediment transport capacity formulations have been used by geomorphologists and engineers to calculate fluvial morphological changes.However,it remains poorly understood if the adaptation to capacity could be fulfilled instantly in response to differing inflow discharges and sediment supplies,and thus if the calculation of morphological changes in rivers based on the assumed capacity status is fully justified.Here we present a numerical investigation on this issue.The distance required for sediment transport to adapt to capacity(i.e.,adaptation-to-capacity length) of both bed load and suspended sediment transport is computationally studied using a coupled shallow water hydrodynamic model,in line with varied inlet sediment concentrations.It is found that the adaptation-to-capacity length generally decreases as the Rouse number increases,irrespective of whether the inlet sediment concentration increases or reduces.For cases with vanishing inlet sediment concentration a unified relationship is found between the adaptation-to-capacity length and the Rouse number.Quantitatively,the adaptation-to-capacity length of bed load sediment is limited to tens of times of the flow depth,whilst that of suspended sediment increases substantially with decreasing Rouse number and can be up to hundreds of times of the flow depth.The present finding concurs that bed load sediment transport can adapt to capacity much more rapidly than suspended sediment transport,and it facilitates a quantitative criterion on which the applicability of bed load or suspended sediment transport capacity for natural rivers can be readily assessed.

Numerical Analysis of Emergency River Restoration Scheme for Qingping Mega Debris Flow

The mega debris flow occurred on August 13 th 2010 in Qingping town,China(hereafter called '8.13' Debris Flow) have done great damage to the local habitants as well as to the re-construction projects in the quake-hit areas,and the channel-fill deposit problem caused by the debris flow was the most destructive.Moreover,it is of high possibility that an even severe deposit problem would reappear and result in worse consequences.In order to maximize risk reduction of this problem,relevant departments of the government established a series of emergency river restoration schemes,for which the numerical analysis is an important procedure to evaluate and determine the optimized one.This study presents a numerical analysis by applying a twodimensional debris flow model combined with a relevant water-sediment model to simulate the deposit during the progress of the debris flow,and to calculate and analyze the river flow field under both the present condition and different restoration conditions.The results show that the debris flow model,which takes the confluence of the Wenjia Gully to the main river into account,could simulate the deposit process quite well.In the reproduced debris flow from the simulation of the '8.13' Debris Flow,the original river flow path has switched to a relatively lower place just along the right bank with a high speed of near 7m.s-1 after being blocked by the deposit,which is highly hazardous.To prevent this hazard,a recommended scheme is derived through inter-comparison of different restoration conditions.It shows that the recommended scheme is able to reduce the water level and as well to regulate the flow path.Based on the given conditions of the mainstream and the tributary confluence for the simulated '8.13' Debris Flow,when encountering a debris flow with deposit volume less than 0.5 million m3,the river channel can endure a 20-year return flood;however,when the deposit volume increases to 2 million m3,the flood capacity of the river will be greatly impacted and the scheme becomes invalid.The recommended scheme supported by the present study has been applied to the emergency river restoration after this mega-debris flow.

Numerical analysis on car-following traffic flow models with delay time

Effects of the speed relaxation time on the optimal velocity car-following model (OVM) with delay time due to driver reaction time proposed by Bando et al.(1995) were studied by numerical methods. Results showed that the OVM including the delay is not physically sensitive to the speed relaxation times. A modified car-following model is proposed to overcome the deficiency. Analyses of the linear stability of the modified model were conducted. It is shown that coexisting flows appear if the initial homogeneous headway of the traffic flow is between critical values. In addition, phase transitions occur on varying the initially homogeneous headway.

Numerical analysis of heat transfer intensity from twin slot vertical jet impingement on strip surface after hot rolling

The flow field and heat transfer of the strip surface due to the twin slot vertical jet impingement were investigated using the ANSYS FLUENT.The RNG k-ε model was carried out in the turbulent calculation.Systematic parametric research was conducted by varying the jet velocity of nozzle exit(V=5 m/s,7.5 m/s,10 m/s),the temperature of cooling water(T_w=280 K,300 K),the normalized spacing from the nozzle to the strip surface(H=10,15,20,33),and the normalized spacing from the nozzle to nozzle centerline(W=0,15,30).The velocity streamline of the flow domain and the general trend of the distribution of the local Nusselt number on the impingement surface of strip were obtained.The result indicate that,the average Nusselt number increases by about70%(90%) as the jet velocity is increased from 5 m/s to 7.5 m/s(from 7.5 m/s to 10 m/s),and T_w,Hand//have minimal effect on it.While the valley Nusselt number decreases by about 10%-43%with the increase of H and W.The functional relationship between the average Nusselt number and the systematic parameters is derived by the least square regression method.

NUMERICAL MANIFOLD METHOD AND ITS APPLICATION IN UNDERGROUND OPENINGS

A brief introduction is made for the Numerical Manifold Method and its analysing process in rock mechanics. Some aspects of the manifold method are improved in implementing process according to the practice of excavating underground openings. Corresponding formulas are given and a computer program of the Numerical Manifold Method has been completed in this paper.

Numerical Analysis of the Nonlinear Parameterization of Waves in Currents over a Submerged Sill with a Non-Hydrostatic Model

An investigation of the effects of a uniform current strength direction(following or opposing wave propagation) on the nonlinear transformation of irregular waves over a submerged trapezoidal sill is carried out using SWASH,a non-hydrostatic numerical wave model.The nonlinear parameters(i.e.,asymmetry,skewness,and kurtosis) are calculated,and the empirical formulas for these parameters are presented as a function of the local Ursell number based on the present numerical data measured.In the shoaling area of the submerged sill,the nonlinear characteristics of waves are more obvious when waves propagate in the same direction as the currents than when waves propagate in the opposite direction.Whereas nonlinear parameters grow with the strengthening of the following currents over the crest,they tend to decrease as the adverse current velocity increases over the crest area of the submerged sill.