“流场法”在水库渗漏管涌通道探测中的应用

阐述了"流场法"的基本原理及探测数据的分析方法,通过对朝阳县北地水库大坝渗漏入水口的探测实例分析,认为"流场法"针对水库坝后渗漏出水点,可以查明坝体的渗漏进水区域,可为水库除险加固治理提供可靠依据。

Phase-field simulations of forced flow effect on dendritic growth perpendicular to flow

The effect of supercooled melt forced laminar flow at low Reynolds Number on dendritic growth perpendicular to melt flow direction was investigated with the phase-field method by incorporating melt convection and thermal noise under non-isothermal condition. By taking the dendritic growth of high pure succinonitrile （SCN） supercooled melt as an example, side-branching shape difference of melts with flow and without flow was analyzed. Relationships among supercooled melt inflow velocity, deflexion angle of dendritic arm and dendritic tip growth velocity were studied. Results show that the melt inflow velocity has few effects on the dendritic tip growth velocity. A formula of relationship between the velocity of the melt in front of primary dendritic tip and the dendritic growth time was deduced, and the calculated result was in quantitative agreement with the simulation result.

Phase-field simulation of forced flow effect on random preferred growth direction of multiple grains

The random distribution problem of dendrite preferred growth direction was settled by random grid method.This method was used to study the influence of forced laminar flow effect on multiple grains during solidification.Taking high pure succinonitrile （SCN） undercooled melt as an example,the forced laminar flow effect on multiple grains was studied by phase-field model of single grain which coupled with flow equations at non-isothermal condition.The simulation results show that the random grid method can reasonably settle the problem of random distribution and is more effective.When the solid fraction is relatively low,melt particles flow around the downstream side of dendrite,and the flow velocity between two dendrite arms becomes high.At the stage of solidification time less than 1800Δt,every dendrite grows freely;the upstream dendrites are stronger than the downstream ones.The higher the melt flow rate,the higher the solid fraction.However,when the solid fraction is relatively high,the dendrite arm intertwins and only a little residual melt which is not encapsulated can flow;the solid fraction will gradually tend to equal to solid fraction of melt without flow.

Effect of forced lamina flow on microsegregation simulated by phase field method quantitatively

The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.

TVD SCHEME WITH CHEMICAL REACTION FLOW AND ITS APPLICATION IN COMBUSTION GAS JETTING FLOW

By conjugating features of combustion gas jetting flows of the solid-rocket and using mathematical methods, a numerical scheme is systematically derived based on Harten′s standard TVD scheme, which fits for the flow with high temperature, pressure and velocity. The rational calculation formula of pressure partial derivation is also given out. By using the chemical kinetics knowledge, problems of multi-component and finite rate chemical reaction contained in combustion gas of the rocket flow field are discussed. The method for solving the mass source term of chemical reaction is clarified. Taking 9 reaction equations with 12 components as an example and utilizing the established calculation program, the free jetting flow field of the rocket is simulated. Numerical results show the correctness of the numerical scheme.

Analytical and Experimental Studies of Linear MHD Propulsor

A two dimensional mathematical model was developed to predict the performance characteristics for direct current, linear channel MHD propulsion system in a closed loop environment. The results of analytical and experimental studies of the linear channel MHD propulsor are described. Compared with the data of experiment, the correctness of the computation program is validated.

Flow and natural convection heat transfer characteristics of non-Newtonian nanofluid flow bounded by two infinite vertical flat plates in presence of magnetic field and thermal radiation using Galerkin method

The main goal of this paper is to investigate natural convective heat transfer and flow characteristics of non-Newtonian nanofluid streaming between two infinite vertical flat plates in the presence of magnetic field and thermal radiation.Initially,a similarity transformation is used to convert momentum and energy conservation equations in partial differential forms into non-linear ordinary differential equations (ODE) applying meaningful boundary conditions.In order to obtain the non-linear ODEs analytically,Galerkin method (GM) is employed.Subsequently,the ODEs are also solved by a reliable numerical solution.In order to test the accuracy,precision and reliability of the analytical method,results of the analytical analysis are compared with the numerical results.With respect to the comparisons,fairly good compatibilities with insignificant errors are observed.Eventually,the impacts of effective parameters including magnetic and radiation parameters and nanofluid volume fraction on the velocity,skin friction coefficient and Nusselt number distributions are comprehensively described.Based on the results,it is revealed that with increasing the role of magnetic force,velocity profile,skin friction coefficient and thermal performance descend.Radiation parameter has insignificant influence on velocity profile while it obviously has augmentative and decreasing effects on skin friction and Nusselt number,respectively.

Simulation of faceted dendrite growth of non-isothermal alloy in forced flow by phase field method

Numerical simulation based on a new regularized phase field model was presented to simulate the dendritic shape of a non-isothermal alloy with strong anisotropy in a forced flow. The simulation results show that a crystal nucleus grows into a symmetric dendrite in a free flow and into an asymmetry dendrite in a forced flow. As the forced flow velocity is increased, both of the promoting effect on the upstream arm and the inhibiting effects on the downstream and perpendicular arms are intensified, and the perpendicular arm tilts to the upstream direction. With increasing the anisotropy value to 0.14, all of the dendrite arms tip velocities are gradually stabilized and finally reach their relative saturation values. In addition, the effects of an undercooling parameter and a forced compound flow on the faceted dendrite growth were also investigated.

Influence of fins on tractor-type podded propulsor performance

A mathematical model of podded propulsors was established in order to investigate the influence of fins. The hydrodynamic performance of podded propulsors with and without fins was calculated, with interactions between propellers and pods and fins derived by iterative calculation. The differential equation based on velocity potential was adopted and hyperboloidal panels were used to avoid gaps between surface panels. The Newton-Raphson iterative procedure was used on the trailing edge to meet the pressure Kutta condition. The velocity distribution was calculated with the Yanagizawa method to eliminate the singularity caused by use of the numerical differential. Comparisons of the performance of podded propulsors with different fins showed that the thrust of propeller in a podded propulsor with fins is greater. The resistance of the pod is also reduced because of the thrust of the fin. The hydrodynamic performance of a podded propulsor with two fins is found to be best, the performance of a podded propulsor with one fin is not as good as two fins, and the performance of the common type is the worst.

Numerical simulation of transport phenomena during strip casting with EMBr in a single belt caster

A theoretical investigation of fluid flow,heat transfer and solidification(solidification transfer phenomena,STP)was presented which coupled with direct-current(DC)magnetic fields in a high-speed strip-casting metal delivery system.The bidirectional interaction between the STP and DC magnetic fields was simplified as a unilateral one,and the fully coupled solidification transport equations were numerically solved by the finite volume method(FVM).While the magnetic field contours for a localized DC magnetic field were calculated by software ANSYS and then incorporated into a three-dimensional(3-D)steady model of the liquid cavity in the mold by means of indirect coupling.A new FVM-based direct-SIMPLE algorithm was adopted to solve the iterations of pressure-velocity(P-V).The braking effects of DC magnetic fields with various configurations were evaluated and compared with those without static magnetic field(SMF).The results show that 0.6 T magnetic field with combination configuration contributes to forming an isokinetic feeding of melt,the re-circulation zone is shifted towards the back wall of reservoir,and the velocity difference on the direction of height decreases from 0.1 m/s to 0.Furthermore,the thickness of solidified skull increases uniformly from 0.45 mm to 1.36 mm on the chilled substrate(belt)near the exit.

After-market liquidity and IPOs underpricing： Evidence from Shari＇ah and non Shari＇ah-based firms

The underpricing of initial public offerings （IPOs） is generally explained with asymmetric information and risk. We complement these traditional explanations with a new theory proposed by Ellul and Pagano （2006） where investors worry also about the after-market illiquidity that may result from asymmetric information after the IPO. The less liquid the after-market is expected to be, the larger will be the IPO underpricing. The samples are the 41 IPOs carried out between 2001-2005. The samples are 7 Shari＇ah-based firms and 34 non Shariah-based firms. Shariah-based firms are those included in Jakarta Islamic Index （JII）, at least one period （one semester）. Regression results show that the relationship between after-market liquidity and underpricing is insignificant unless we use trading frequency as proxy for liquidity for non Shariah-based firms.

Optimal Design for Open MRI Superconducting Magnet with Active Shielding

The optimal design method for an open Magnetic Resonance Imaging （MRI） superconducting magnet with an active shielding configuration is proposed. Firstly, three pairs of current rings are employed as seed coils. By optimizing the homogeneity of Diameter Sphere Voltnne （DSV）, the positions and currents of the seed coils will be obtained. Secondly, according to the positions and currents of the seed coils, the current density of superconducting wires is determined, and then the original sections for the coils can be achieved. An optimization for the homogeneity based on the constrained nonlincar optimization method is employed to determine the coils with good homogeneity. Thirdly, the magnetic field generated by previous coils is set as the background field, then add two coils with reverse current, and optimize the stray field line of 5 Gauss in a certain scope. Finally, a further optimization for the homogeneity is used to get Final coils. This method can also be used in the design of other axisynmaetfic superconducting MRI magnets.

Numerical Prediction of Flow Patterns after Various Pipe Fittings

An accurate prediction of flows using CFD depends on a large number of factors. In addition to discretizing the flow region, the correct definition of boundary or initial conditions and the choice of suitable numerical methods, the applied turbulence model influences the results of the flow simulation to a great extent. Therefore, a validation of the results with the experimental data is of great importance for a correct selection of a turbulence model. It is the scope of this paper to assess different turbulence models for the simulation of pipe flows. The calculation results of pipe flows through a combination of 90~ elbows and a 1/3 segmental orifice are compared with experimental measurement results. This has the advantage that the suitability of the turbulence models for simulating both shear and swirl flows can be investigated. Thus, the k-ω, k-ε model and the Launder Reece Rodi Reynolds stress model are compared with each other and experimental results. Furthermore, this investigation is extended through including a much more c detached-eddy simulation. This model provides better prediction of the flow by resolving the large eddies and modeling the small ones. The experimental results originate from LDV measurements over the entire pipe cross-section. This measuring method provides velocity vectors over the measured surface.

The study of slip line field and upper bound method based on associated flow and non-associated flow rules

At present, associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line. It is proved that geotechnical materials do not abide by the associated flow rule. It is impossible for the stress characteristic line to conform to the velocity line. Generalized plastic mechanics theoretically proved that plastic potential surface intersects the Mohr-Coulomb yield surface with an angle, so that the velocity line must be studied by non-associated flow rule. According to limit analysis theory, the theory of slip line field is put forward in this paper, and then the ultimate beating capacity of strip footing is obtained based on the associated flow rule and the non-associated flow nile individually. These two results are identical since the ultimate bearing capacity is independent of flow role. On the contrary, the velocity fields of associated and non-associated flow rules are different which shows the velocity field based on the associat- ed flow rule is incorrect.

Characteristics of gas explosion flow fields in complex pipelines

The explosion flow field in five straight pipes with different diameters and one bending pipe selected from a domestic coal mine are studied by the method of numerical simulation. And the results show that,both in the straight and bending pipes, the pressure wave and velocity wave are accelerated by the rising of reaction rate. As the explosion progressed, with the temperature reaching approximately 3000 K, only one pressure wave and one reaction rate wave were observed, while several velocity waves were found.The larger diameter presented the highest relative pressure as well as the largest velocity increase and subsequent decrease inside the tube. The bent pipes caused both turbulence and kinetic energy to increase, resulting in the acceleration of the reaction rate. The burning time was 7.4% shorter than the burning time observed for the straight pipe. Based on these results, designing one explosion resistance device, and in the practical engineering applications, it was to be proved to meet the security requirements fully.

Simulation and experiment of starting transient flow field of the hydrostatic bearing based on dynamic mesh method

A new method is developed to assess and analyze the dynamic performance of hydrostatic bearing oil film by using an amulets-layer dynamic mesh technique. It is implemented using C Language to compile the UDF program of a single oil film of the hydrostatic bearing. The effects of key lubrication parameters of the hydrostatic bearing are evaluated and analyzed under various working conditions,i.e. under no-load,a load of 40 t,a full load of 160 t,and the rotation speed of 1r/min,2r/min,4r/min,8r/min,16r/min,32r/min. The transient data of oil film bearing capacity under different load and rotation speed are acquired for a total of 18 working conditions during the oil film thickness changing. It allows the effective prediction of dynamic performance of large size hydrostatic bearing. Experiments on hydrostatic bearing oil film have been performed and the results were used to define the boundary conditions for the numerical simulations and validate the developed numerical model. The results showed that the oil film thickness became thinner with the increase of the operating time of the hydrostatic bearing,both the oil film rigidity and the oil cavity pressure increased significantly,and the increase of the bearing capacity was inversely proportional to the cube of the change of the film thickness. Meanwhile,the effect of the load condition on carrying capacity of large size static bearing was more important than the speed condition. The error between the simulation value and the experimental value was 4.25%.

Installation position determination of wind speed sensors on steel pole along a high-speed railway

In order to consider the influence of steel pole on the measurement of wind speed sensors and determinate the installation position of wind speed sensors, the flow field around wind speed sensors was investigated. Based on the three-dimensional steady Reynolds-averaged Navier-Stokes equations and k-ε double equations turbulent model, the field flow around the wind speed sensor and the steel pole along a high-speed railway was simulated on an unstructured grid. The grid-independent validation was conducted and the accuracy of the present numerical simulation method was validated by experiments and simulations carried out by previous researchers. Results show that the steel pole has a significant influence on the measurement results of wind speed sensors. As the distance between two wind speed sensors is varied from 0.3 to 1.0 m, the impact angles are less than ±20°, it is proposed that the distance between two wind speed sensors is 0.8 m at least, and the interval between wind speed sensors and the steel pole is more than 1.0 m with the sensors located on the upstream side.

Phase-field simulation of dendritic growth in a forced liquid metal flow coupling with boundary heat flux

A phase-field model with forced liquid metal flow was employed to study the effect of boundary heat flux on the dendritic structure forming of a Ni-40.8%Cu alloy with liquid flow during solidification.The effect of the flow field coupling with boundary heat extractions on the morphology change and distributions of concentration and temperature fields was analyzed and discussed.The forced liquid flow could significantly affect the dendrite morphology,concentration and temperature distributions in the solidifying microstructure.And coupling with boundary heat extraction,the solute segregation and concentration diffusion were changed with different heat flux.The morphology,concentration and temperature distributions were significantly influenced by increasing the heat extraction,which could relatively make the effect of liquid flow constrained.With increasing the initial velocity of liquid flow,the lopsided rate of the primary dendrite arm was enlarged and the transition of developing manner of the secondary arms moved to the large heat extraction direction.It was the competition between heat flux and forced liquid flow that finally determined microstructure forming during solidification.

Numerical research on mixing characteristics of different injection schemes for supersonic transverse jet

A numerical method using AUSMDV scheme and k-ω SST turbulence model with an explicit compressibility correction was developed,and a 3-D numerical simulation of a supersonic flow field with a vertical sonic jet of hydrogen was performed.Good agreement between numerical results and experimental data validated the reliability of the numerical method.Whereafter,two parameters,mass-weighted average total pressure and mixing efficiency,were defined to evaluate the mixing performance of different injection schemes.Based on the numerical method and evaluation criterion,the mixing characteristics of different injection schemes were studied in detail.It was found that for the mixing field of supersonic transverse jet,the near-field mixing is controlled by convection transport while the far-field mixing is controlled by mass diffusion;the circular-hole injection causes a loss of total pressure comparable to the slot injection,but can induce a much higher mixing efficiency because of its 3-D flow characteristic;the variation of injection angle under circular-hole injection mainly affects the near-field mixing degree,and among the five injection angles studied in the present paper,angle 120° is the optimal one;with the increase of the ratio between injector space and diameter,the induced mixing efficiency increases while the caused loss of total pressure can grow greatly;the two-stage injection method designed through reducing the injector area to keep the same hydrogen mass flowrate can induce a much higher mixing efficiency while only a bit larger loss of total pressure when compared to the single-stage injection,and hence the two-stage injection is superior to the single-stage injection.The research results can direct the design of the fuel injection method in the combustor of scramjet engine.

Comparison of Reynolds Averaged Navier-Stokes Based Simulation and Large-eddy Simulation for One Isothermal Swirling Flow

The flow structure of one isothermal swirling case in the Sydney swirl flame database was studied using two numerical methods. Results from the Reynolds-averaged Navier-Stokes (RANS) approach and large eddy simulation (LES) were compared with experimental measurements. The simulations were applied in two different Cartesian grids which were investigated by a grid independence study for RANS and a post-estimator for LES. The RNG k-ε turbulence model was used in RANS and dynamic Smagorinsky-Lilly model was used as the sub-grid scale model in LES. A validation study and cross comparison of ensemble average and root mean square (RMS) results showed LES outperforms RANS statistic results. Flow field results indicated that both approaches could capture dominant flow structures, like vortex breakdown (VB), and precessing vortex core (PVC). Streamlines indicate that the formation mechanisms of VB deducted from the two methods were different. The vorticity field was also studied using a velocity gradient based method. This research gained in-depth understanding of isothermal swirling flow.