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.

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.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

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.

Hybrid Methods for Computing the Streamfunction and Velocity Potential for Complex Flow Fields over Mesoscale Domains

Three types of previously used numerical methods are revisited for computing the streamfunctionψand velocity potentialχfrom the horizontal velocity v in limited domains.The first type,called the SOR-based method,uses a classical successive over-relaxation(SOR)scheme to computeψ(orχ)first with an arbitrary boundary condition(BC)and thenχ(orψ)with the BC derived from v.The second type,called the spectral method,uses spectral formulations to construct the inner part of(ψ,χ)-the inversion of(vorticity,divergence)with a homogeneous BC,and then the remaining harmonic part of(ψ,χ)with BCs from v.The third type,called the integral method,uses integral formulas to compute the internally induced(ψ,χ)-the inversion of domain-internal(vorticity,divergence)using the free-space Greenꞌs function without BCs and then the remaining harmonicψ(orχ)with BCs from v minus the internally-induced part.Although these methods have previously been successfully applied to flows in large-scale and synoptic-scale domains,their accuracy is compromised when applied to complex flows over mesoscale domains,as shown in this paper.To resolve this problem,two hybrid approaches,the integral-SOR method and the integral-spectral method,are developed by combining the first step of the integral method with the second step adopted from the SOR-based and spectral methods,respectively.Upon testing these methods on real-case complex flows,the integral-SOR method is significantly more accurate than the integral-spectral method,noting that the latter is still generally more accurate than the three previously-used methods.The integral-SOR method is recommended for future applications and diagnostic studies of complex flows.

Combined Application of Wide-Field Electromagnetic Method and Flow Field Fitting Method for High-Resolution Exploration： A Case Study of the Anjialing No. 1 Coal Mine

The Anjialing No. 1 Coal Mine in Shanxi Province, China, contains a complicated old goaf and an unknown water distribution that hold high potential for serious water hazards. Due to poor detection resolution, previous attempts have failed to determine the scope of the old goal and the water distribution in the mine by separate use of various exploration methods such as seismic method, direct current resistivity, audio magnetotellurics, controlled-source audio-frequency magnetotellurics, and transient electromag-netics. To solve this difficult problem, a combination of the wide-field electromagnetic method and the flow field fitting method with three-dimensional resistivity data inversion was applied to determine the precise scope of the goal and the locations where water is present, and to identify the hydraulic con- nection between the water layers so as to provide reliable technical support for safe coal production. Reasonable results were achieved, with all these goals being met. As a result, a mining area of nearly 4 km^2 has been released for operation.

Effect of fluid flow on solidified equiaxed dendrite morphology evolution based on phase field-lattice Boltzmann method

Fluid flow can significantly change the evolution of microstructural morphology. However, relatively little is known how the fluid flow, concentration and microstructure affect each other quantitatively, which is essential to optimize processing parameters. A quantitative simulation study of Al-Cu solidified equiaxed dendrite evolution under forced flow based on phase field-lattice Boltzmann method(PF-LBM) is carried out. Results obtained are validated by Gibbs-Thomson relation at the dendrite tip. Compared with the equiaxed dendrite evolution without flow, the upstream dendrite arm is enhanced while the downstream arm is inhibited. Besides, as the inlet flow rate increases, the secondary arms attached onto the upstream primary arm and the upstream side of the primary arm normal to the inflow velocity has been well developed. Results show that sidewise instabilities of the primary dendrite arm and onset of secondary arm is caused by the local concentration perturbation and will be enhanced or inhabited by the flow. It is believed that the coupled PF-LBM method is able to handle dendrite evolution under forced flow quantitatively, which helps in investigating the solidified dendrite morphology evolution.

Simulation of Temperature and Flow Field by Three Dimension Finite ElementMethod for Castex Process of AS Wire

Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with the theory of heat-transfer and hydrodynamics by means of 3-dimensional finite element method. From simulation results, it is found that the variation of temperature and velocityis mainly influenced by the casting temperature of aluminum, rotating speed of Castex wheel and flow of cooling water. Among theseinfluencing factors, the casting temperature distributes most to the length of liquid phase metal. Moreover, the faster the metal solidifies,the higher the metal there moves with the overall trend of descending from the bottom of the wheel to the shoe wall as well as from sidewalls to the center of wheel groove. In comparison with the practical value, the simulation is reliable.

Simulation of Flow Field of Molten Salt in Neodymium Metal Electrolytic Cell Using Vortex-Flow Function Method

With the applications of Nd-Fe-B material extending in recent years, the materials of neodymium metal and other rare earth metal alloy confront the increased demand and the high quality request at the same time.These factors stimulated greatly to perfect the producing craft of RE metals and improve the equipments.The rare earth electrolysis cell is developing towards large-scale way.Notwithstanding the present electrolysis cell of Nd metal, include 6 kA and 10 kA cell, exists some insurmountable problems during operation and these problems lead to lower electric efficiency and higher operating costs.So it is significant to study the physical fields of rare earth electrolysis cell.In this paper,a numerical flow mode is established using vortex- flowing function method and the fluid flow field of 3000A Nd electrolysis cell is computed using MATLAB.The results of the study will be important reference in theory for improving and enlarging rare earth fluoride system cell.

Flow field fusion simulation method based on model features and its application in CRDM

The control rod drive mechanism(CRDM)is an essential part of the control and safety protection system of pressurized water reactors.Current CRDM simulations are mostly performed collectively using a single method,ignoring the influence of multiple motion units and the differences in various features among them,which strongly affect the efficiency and accuracy of the simulations.In this study,we constructed a flow field fusion simulation method based on model features by combining key motion unit analysis and various simulation methods and then applied the method to the CRDM simulation process.CRDM performs motion unit decomposition through the structural hierarchy of function-movement-action method,and the key meta-actions are identified as the nodes in the flow field simulation.We established a fused feature-based multimethod simulation process and processed the simulation methods and data according to the features of the fluid domain space and the structural complexity to obtain the fusion simulation results.Compared to traditional simulation methods and real measurements,the simulation method provides advantages in terms of simulation efficiency and accuracy.

RECONSTRUCTION OF PART OF AN ACTUAL BLAST-WAVE FLOW FIELD TO AGREE WITH EXPERIMENTAL DATA BY USING NUMERICAL METHOD WITH HIGH IDENTIFICATION

In this paper, on the basis of experimental data of two kinds of chemical explosions, the piston-pushing model of spherical blast-waves and the second-order Godunov-type scheme of finite difference methods with high identification to discontinuity are used to the numerical reconstruction of part of an actual hemispherical blast-wave flow field by properly adjusting the moving bounary conditions of a piston. This method is simple and reliable. It is suitable to the evaluation of effects of the blast-wave flow field away from the explosion center.

Internal Flow Field Measurement of Gas Turbine Based on Optical Flow Method

The in-cylinder flow field of the internal combustion engine is an important factor affecting the quality and combustion quality of the fuel mixture in the cylinder. In order to calculate the high-precision flow field, the paper presents a flow field calculation method based on the optical flow algorithm. The motion of the point was calculated using the change in pixel intensity within two temporally adjacent frame images. The results show the high accuracy and resolution of the flow field at small displacement conditions.

Finite Element Numerical Simulation and PIV Measurement of Flow Field inside Metering-in Spool Valve

The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.

Numerical study on the characteristics of flow field and wave propagation near submerged breakwater on slope

In this study, characteristics of flow field and wave propagation near submerged breakwater on a sloping bed are investigated with numerical model. The governing equations of the vertical two-dimensional model are Reynolds Averaged Navier Stokes equations. The Reynolds stress terms are closed by a nonlinear k - ε turbulence transportation model. The free surface is traced through the PILC-VOF method. The proposed numerical model is verified with experimental results. The numerical result shows that the wave profile may become more asymmetrical when wave propagates over breakwater. When wave crest propagates over breakwater, the anticlockwise vortex may generate. On the contrary, when wave hollow propagates over breakwater, the clockwise vortex may generate. Meanwhile, the influenced zone of vortex created by wave crest is larger than that created by wave hollow. All the maximum values of the turbulent kinetic energy, turbulent dissipation and eddy viscosity occur on the top of breakwater. Both the turbulent dissipation and eddy viscosity increase as the turbulent kinetic energy increases. Wave energy may rapidly decrease near the breakwater because turbulent dissipation increases and energy in lower harmonics is transferred into higher harmonics.

Effects of operating parameters on the flow field in slab continuous casting molds with narrow widths

Computational simulations and high-temperature measurements of velocities near the surface of a mold were carried out by using the rod deflection method to study the effects of various operating parameters on the flow field in slab continuous casting(CC)molds with narrow widths for the production of automobile exposed panels.Reasonable agreement between the calculated results and measured subsurface velocities of liquid steel was obtained under different operating parameters of the CC process.The simulation results reveal that the flow field in the horizontal plane located 50 mm from the meniscus can be used as the characteristic flow field to optimize the flow field of molten steel in the mold.Increases in casting speed can increase the subsurface velocity of molten steel and shift the position of the vortex core downward in the downward circulation zone.The flow field of liquid steel in a 1040 mm-wide slab CC mold can be improved by an Ar gas flow rate of 7 L·min^−1 and casting speed of 1.7 m·min^−1.Under the present experimental conditions,the double-roll flow pattern is generally stable at a submerged entry nozzle immersion depth of 170 mm.

Density and temperature reconstruction of a flame-induced distorted flow field based on background-oriented schlieren(BOS) technique

An experimental system based on the background-oriented schlieren（BOS） technique is built to reconstruct the density and temperature distribution of a flame-induced distorted flow field which has a density gradient. The cross-correlation algorithm with sub-pixel accuracy is introduced and used to calculate the background-element displacement of a disturbed image and a fourth-order difference scheme is also developed to solve the Poisson equation. An experiment for a disturbed flow field caused by a burning candle is performed to validate the built BOS system and the results indicate that density and temperature distribution of the disturbed flow field can be reconstructed accurately. A notable conclusion is that in order to make the reconstructed results have a satisfactory accuracy, the inquiry step length should be less than the size of the interrogation window.

Flow Field Analysis of Submerged Horizontal Plate Type Breakwater

Submerged horizontal plate can be considered as a new concept breakwater. In order to reveal the wave elimination mechanism of this type breakwater, boundary element method is utilized to investigate the velocity field around plate carefully. The flow field analysis shows that the interaction between incident wave and reverse flow caused by submerged plate will lead to the formation of wave elimination area around both sides of the plate. The velocity magnitude of flow field has been reduced and this is the main reason of wave elimination.

Effects of harmonic electric field frequency on electromagnetic and flow fields of Ti-Al melt in a rectangular mold based on numerical simulation

A model was established based on Maxwell's equations and Navier-Stokes' equations to numerically simulate the electromagnetic field and flow field in a rectangular mold with sectional aspect ratio of 5:1. The FEM (Finite Element Method) and APDL (ANSYS Parametric Design Language) were employed for the model to execute the modeling, meshing, load applying and solving. The Ti-Al alloy melt was selected to illustrate and validate the effects of the harmonic field frequency on the distribution of the physical fields in the mold. The simulated results demonstrate that with an increasing frequency the electric current forms an ellipsoid cavity where it becomes much weaker, and that the melt flows more intensely with low frequency (less than 5 kHz) than with high frequency (more than 5 kHz). The melt is pinched from the central part in the mold to bipolar parts in which it forms two vortexes in each side. The maximum value of fluid velocity exists near the bipolar zone.

Numerical simulation of a falling droplet surrounding by air under electric field using VOF method:A CFD study

This is a numerical study of a falling droplet surrounding by air under the electric field modeled with finite volume method by means of CFD.The VOF method has been employed to model the two-phase flow of the present study.Various capillary numbers are investigated to analyze the effects of electric field intensity on the falling droplet deformation.Also,the effects of electric potential on the heat transfer coefficient have been examined.The obtained results showed that by applying the electric field at a capillary number of 0.2 the droplet tends to retain its primitive shape as time goes by,with a subtle deformation to an oblate form.Intensifying the electric field to a capillary number of 0.8 droplet deformation is almost insignificantwith time progressing;however,further enhancement in capillary number to 2 causes the droplet to deform as a prolate shape and higher values of this number intensify the prolate form deformation of the droplet and result in pinch-off phenomenon.Ultimately,it is showed that as the electric potential augments the heat transfer coefficient increases in which for electric potential values higher than 2400 V the heat transfer coefficient enhances significantly.