Coupled Numerical Simulation of Electromagnetic and Flow Fields in a Magnetohydrodynamic Induction Pump

Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.

Flow field, sedimentation, and erosion characteristics around folded linear HDPE sheet sand fence: Numerical simulation study

Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.

Microwave-induced thermoacoustic elastic imaging:A simulation study

Microwave-induced thermoacoustic imaging(MTI)has the advantages of high resolution,high contrast,non-ionization,and non-invasive.Recently,MTI was used in the¯eld of breast cancer screening.In this paper,based on the¯nite element method(FEM)and COMSOL Multiphysics software,a three-dimensional breast cancer model suitable for exploring the MTI process is proposed to investigate the in°uence of Young's modulus(YM)of breast cancer tissue on MTI.It is found that the process of electromagnetic heating and initial pressure generation of the entire breast tissue is earlier in time than the thermal expansion process.Besides,compared with normal breast tissue,tumor tissue has a greater temperature rise,displacement,and pressure rise.In particular,YM of the tumor is related to the speed of thermal expansion.In particular,the larger the YM of the tumor is,the higher the heating and contraction frequency is,and the greater the maximum pressure is.Di®erent Young's moduli correspond to di®erent thermoacoustic signal spectra.In MTI,this study can be used to judge di®erent degrees of breast cancer based on elastic imaging.In addition,this study is helpful in exploring the possibility of microwave-induced thermoacoustic elastic imaging(MTAE).

An internal ballistic model of electromagnetic railgun based on PFN coupled with multi-physical field and experimental validation

To accelerate the practicality of electromagnetic railguns,it is necessary to use a combination of threedimensional numerical simulation and experiments to study the mechanism of bore damage.In this paper,a three-dimensional numerical model of the augmented railgun with four parallel unconventional rails is introduced to simulate the internal ballistic process and realize the multi-physics field coupling calculation of the rail gun,and a test experiment of a medium-caliber electromagnetic launcher powered by pulse formation network(PFN)is carried out.Various test methods such as spectrometer,fiber grating and high-speed camera are used to test several parameters such as muzzle initial velocity,transient magnetic field strength and stress-strain of rail.Combining the simulation results and experimental data,the damage condition of the contact surface is analyzed.

Effect of traveling-wave magnetic field on dendrite growth of high-strength steel slab: Industrial trials and numerical simulation

The dendrite growth behavior of high-strength steel during slab continuous casting with a traveling-wave magnetic field was studied in this paper. The morphology of the solidification structure and composition distribution were analyzed. Results showed that the columnar crystals could deflect and break when the traveling-wave magnetic field had low current intensity. With the increase in current intensity, the secondary dendrite arm spacing and solute permeability decreased, and the columnar crystal transformed into an equiaxed crystal. The electromagnetic force caused by the traveling-wave magnetic field changed the temperature gradient and velocity magnitude and promoted the breaking and fusing of dendrites. Dendrite compactness and composition uniformity were arranged in descending order as follows:columnar-toequiaxed transition (high current intensity), columnar crystal zone (low current intensity), columnar-to-equiaxed transition (low current intensity), and equiaxed crystal zone (high current intensity). Verified numerical simulation results combined with the boundary layer theory of solidification front and dendrite breaking–fusing model revealed the dendrite deflection mechanism and growth process. When thermal stress is not considered, and no narrow segment can be found in the dendrite, the velocity magnitude on the solidification front of liquid steel can reach up to 0.041 m/s before the dendrites break.

Full-Scale Numerical Simulation of the Local Scour Under Combined Current and Wave Conditions Based on Field Data

The monopile is the most common foundation to support offshore wind turbines.In the marine environment,local scour due to combined currents and waves is a significant issue that must be considered in the design of wind turbine foundations.In this paper,a full-scale numerical model was developed and validated based on field data from Rudong,China.The scour development around monopiles was investigated,and the effects of waves and the Reynolds number Re were analyzed.Several formulas for predicting the scour depth in the literature have been evaluated.It is found that waves can accelerate scour development even if the KC number is small(0.78<KC<1.57).The formula obtained from small-scale model tests may be unsafe or wasteful when it is applied in practical design due to the scale effect.A new equation for predicting the scour depth based on the average pile Reynolds number(Rea)is proposed and validated with field data.The equilibrium scour depth predicted using the proposed equation is evaluated and compared with those from nine equations in the literature.It is demonstrated that the values predicted from the proposed equation and from the S/M(Sheppard/Melville)equation are closer to the field data.

Numerical Simulation of 3D Flow Field and Flow-Induced Noise Characteristics in a T-Shaped Reducing Tee Junction

The so-called T-shaped reducing tees are typically used to divide,change and control(to a certain extent)the flow direction in pipe networks.In this study,the Ffowcs Williams–Hawkings(FW-H)equation and the Large Eddy Simulation(LES)methods are used to simulate the flow-induced noise related to T-shaped reducing tees under different inlet flow velocities and for different pipe diameter ratios.The results show that the maximum flow velocity,average flow velocity,and vorticity in the branch pipe increase gradually as the related diameter decreases.Strong vorticity and secondary flows are also observed in the branch pipe,and the associated violent pressure fluctuations are found to be the main sources of flow-induced noise.In particular,as the pipe diameter ratio decreases from 1 to 0.45,the Total Sound Pressure Level(TSPL)increases by 6.8,6.26,and 7.43 dB for values of the inlet flow velocity of 1,2,and 3 m/s,respectively.The distribution characteristics of the flow-induced noise in the frequency domain follow similar trends for different pipe diameter ratios.

Simulation for MSS-2 low-perigee elliptical orbit satellites:an example of lithospheric magnetic field modelling

A future constellation of at least four geomagnetic satellites(designated Macao Scientific Satellite-1(MSS-1)and Macao Scientific Satellite-2(MSS-2))was recently proposed,to continue high-quality geomagnetic observations in the post-Swarm period,focusing especially on collecting data that will provide a global,three-dimensional survey of the geomagnetic field.In this paper,we present a simulation of two years of orbits(2020.01.01-2022.01.01)of two satellites(tentatively denoted as MSS-2)that are constellated in elliptical(200×5,300 km)low-perigee orbits.By comparing error variances of Gauss coefficients,we investigate the sensitivity of lithospheric magnetic field modelling to data collected from various satellite orbits,including a near circular reference orbit of 300×350km,and elliptical orbit of 180×5,300 km,220×5,300 km,200×3,000 km and 200×1,500 km.We find that in two years the two MSS-2 satellites can collect 35,000 observations at altitude below 250 km,data that will be useful in advancing the quality of lithospheric magnetic field modelling;this number of observations reflects the fact that only 4.5%of the flight time of these satellites will be below250 km(just 6.4%of their flight time below 300 km).By combining observations from the MSS-2 satellites’elliptical orbits of 200×5,300km with observations from a circular reference orbit,the variance of the geomagnetic model can be reduced by a factor of 285 at spherical harmonic degree n=200 and by a factor of 1,300 at n=250.The planned lower perigee of their orbits allows the new satellites to collect data at unprecedentedly lower altitudes,thus dramatically improving the spatial resolution of satellite-derived lithospheric field models,(up to 80%at n=150).In addition,lowering the apogee increases the time interval during which the satellites fly at near-Earth altitudes,thus improving the model predictions at all spherical harmonic degrees(around 52%-62%at n=150).The upper limit of the expected improvement to the field model at the orbital apogee is not as good as at the perigee.However,data from the MSS-1 orbit can help fill the gap between data from the MSS-2 orbits and from the circular reference orbit for the low-degree part of the model.The feasibility of even lower-altitude flight requires further discussion with satellite engineers.

Phase field simulation of grain refinement in silver-based filler metal

Numerical simulation is one of the important auxiliary methods for studying materials-related problems. In this study, phase field simulation was employed to investigate the refinement behavior of BAg55CuZn-x B brazing alloys. Simulation and experimental studies were conducted for B contents ranging from 0 wt.% to 0.2 wt.%. The results demonstrated that the addition of 0.05 wt.% B in the brazing alloy leads to a significant refinement effect. As the B content increases, the grain size further reduces, and a refinement stagnation phenomenon occurs after exceeding 0.15 wt.%. The solidification process of brazing alloys with different B content was predicted by simulation, and the simulation results showed that with the increase of B content, the initial number of nucleation increased, and the radius of the dendrite tip decreased. The simulation results are in good agreement with the experimental findings, providing further evidence of the refining effect of the B element and the reliable predictive capability of the phase field model.

COUPLED SIMULATION OF 3D ELECTRO-MAGNETO-FLOW FIELD IN HALL-HEROULT CELLS USING FINITE ELEMENT METHOD

Two full 3D steady mathematical models are developed by finite element method （FEM） to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.

Numerical Simulation of Multi-track and Multi-layer Temperature Field on Laser Direct Metal Shaping

To improve the mechanical properties of the parts fabricated by Laser Direct Metal Shaping(LDMS),it is of great significance to understand the distribution regularities of transient temperature field during LDMS process.Based on the“el- ement birth and death”technique of finite element method,a three-dimensional multi-track and multi-layer model for the transient temperature field analysis of LDMS is developed by ANSYS Parametric Design Language(APDL)for the first time.In the fab- ricated modal,X-direction parallel reciprocating scanning paths is introduced.Using the same process parameters,the simulation results show good agreement with the microstructure features of samples which fabricated by LDMS.

2-D Current Field Numerical Simulation Integrating Yangtze Estuary with Hangzhou Bay

In this paper, integrating the Yangtze Estuary with the Hangzhou Bay, a 2-D velocity field model is established. In the model, fine self-adaptive grids are employed to adapt to the complicated coastal shape. The hydrodynamic equations satisfied by two contravariant components of velocity vector and surface elevation in non-orthogonal curvilinear coordinates are used. In each momentum equation the coefficients before the two partial derivatives of surface elevation with respect to variables of alternative direction coordinates have different orders of magnitude, i. e., the derivative with the larger coefficient may play a more important role than that with the smaller one. With this advantage, the ADI scheme can be easily employed. The hydrodynamic factors include tidal current, river runoff and wind-induced current. In terms of tidal current, seven main constituents in the area are considered in the open boundaries. The verifications of surface elevation process and current velocity process in the spring tide and in the neap tide show that the model can preferably reflect current fields in the area. Through the simulation of Lagrangian residual current fields in summer and in winter, the paths of the exchange of water and sediment between the Yangtze Estuary and the Hangzhou Bay are elementarily discussed.

COUPLED NUMERICAL SIMULATION ON COLD ROLLER'S TEMPERATURE FIFLD-PHASE TRANSFORMATION - STRESS FIELD DURING ITS QUENCHING PSOCESS

Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short time, and these variables take a contact interactions in the whole process. In this paper, a three dimensional non - linear mathematical model for queeching process has been founded and the numerical simulation on temperature field,microstructre and stress field has been realized.In the FEM analysis, the incremental iteration method is used to deal with such complicated nonlinear as boundary nonlinear, physical property nonlinear,transformation nonlinear etc.The effect of stress on transformation kinetics has been considered in the calculation of microstructure. In the stress field anal- ysis,a thermo- elasto - plastic model has been founded, which considers such factors as transforma- tion strain,transformation plastic strain, themal strain and the effect of temperature and transforma- tion on mechanical propertier etc. The transient temperature field, microstructure distribution and stress field of the roller on any time can be displayed vividly,and the cooling curve and the changes of stress on any position can also be given.

Numerical simulation of the temperature fields of stainless steel with different roller parameters during twin-roll strip casting

The temperature field of stainless steel during twin-roll strip casting was simulated by experiment and a finite element （FE） model. By comparing the measured result with the simulated values, it is found that they fit close to each other, which indicates this FE model is effective. Based on this model, the effects of roll gap （t） and roll radius （R） on solidification were simulated. The simulated results give the relationship between t or R and the position of the freezing point. The larger the t is and the smaller the R is, the closer the position of the freezing point is to the exit.

Simulation of Paleotectonic Stress Fields and Distribution Prediction of Tectonic Fractures at the Hudi Coal Mine, Qinshui Basin

Study on tectonic fractures based on the inversion of tectonic stress fields is an effective method. In this study, a geological model was set up based on geological data from the Hudi Coal Mine, Qinshui Basin, a mechanical model was established under the condition of rock mechanics and geostress, and the finite element method was used to simulate the paleotectonic stress field. Based on the Griffith and Mohr-Coulomb criterion, the distribution of tectonic fractures in the Shanxi Formation during the Indosinian, Yanshanian, and Himalayan period can be predicted with the index of comprehensive rupture rate. The results show that the acting force of the Pacific Plate and the India Plate to the North China Plate formed the direction of principal stress is N-S, NW - SE, and NE - SW, respectively, in different periods in the study area. Changes in the direction and strength of the acting force led to the regional gradients of tectonic stress magnitude, which resulted in an asymmetrical distribution state of the stress conditions in different periods. It is suggested that the low-stress areas are mainly located in the fault zones and extend along the direction of the fault zones. Furthermore, the high-stress areas are located in the junction of fold belts and the binding site of multiple folds. The development of tectonic fractures was affected by the distribution of stress intensity and the tectonic position of folds and faults, which resulted in some developed areas with level I and II. There are obvious differences in the development of tectonic fractures in the fold and fault zones and the anticline and syncline structure at the same fold zones. The tectonic fractures of the Shanxi Formation during the Himalayan period are more developed than those during the Indosinian and Yanshanian period due to the superposition of the late tectonic movement to the early tectonic movement and the differences in the magnitude and direction of stress intensity.

SIMULATION ON TEMPERATURE FIELD OF FRICTION STIR WELDED JOINTS OF 2024-T4 Al

The thermal model of FSW based on the thermal elastic-plastic finite element method, and the transient temperature distribution of FS welded joints of 2024-T4 Al was simulated by using this model, which provides useful information for the investigation of FSW process. Simulation resuits show that the temperature distribution of the weld gradually decreases toward periphery in a radiate format, whose center is the probe, and the highest temperature in the weld can reach about 400℃. The initial terminal of the weld is a zone, in which the temperature gradient is great, and defects of the welding are easily produced in this zone. Temperature change at the end of the welded joint is as layer variation, the local serious defects are not easy to produce in this zone.

Physical Simulation of Mold-Filling Processing of Thin-Walled Castings under Traveling Magnetic Field

Mold-filling process of thin-walled castings under the condition of traveling magnetic field has been studied by physical simulation method using gallium melt and fast speed photography. Flow morphology and its formation mechanism were obtained and discussed for thin-walled casting. The influences of magnetic field density on the filling ability, filling velocity and mold filling time have been studied. The differences in filling capability between gravity casting and casting under the traveling magnetic field have been compared. The results indicate that the mold filling ability of the gallium melt increases greatly under the condition of traveling magnetic field; the filling time is shortened from 18 s under gravity field to 3 s under the traveling magnetic field and average flow rate of the melt increases from 1.6 to 8.68 cm3/s; the change law of the cross-section morphology of the gallium melt during the mold filling is that at first, the cross-section area does not change, then it decreases gradually. When the front of the melt reaches the end of the mold cavity, the front melt will backfill the mold; the wider the width of mold cavity, the better the mold filling ability. The mold filling ability of gallium melt in mold with upper magnetic conductor is better than that without upper magnetic conductor.

Coherency matrix-based proper orthogonal decomposition with application to wind field simulation

Proper Orthogonal Decomposition （POD） provides a powerful modal transformation tool for stochastic dynamics. In this paper, coherency matrix-based proper orthogonal decomposition （CPOD） is presented as an innovative form of the POD based on cross power spectral density matrices. By introducing a discretizing scheme, the CPOD-based spectral representation method is obtained for use in stochastic simulation. Moreover, some criteria are proposed that allow the truncation order of CPOD to be conveniently determined. A numerical example to illustrate the application of the proposed method for the simulation of a wind velocity field is provided.

Phase field simulation of single bubble behavior under an electric field

Based on the Cahn-Hilliard phase field model, a three-dimensional multiple-field coupling model for simulating the motion characteristics of a rising bubble in a liquid is established in a gas-liquid two-phase flow. The gas-liquid interface motion is simulated by using a phase-field method, and the effect of the electric field intensity on bubble dynamics is studied without electric field, or with vertical electric field or horizontal electric field. Through the coupling effect of electric field and flow field, the deformation of a single rising bubble and the formation of wake vortices under the action of gravity and electric field force are studied in detail. The correctness of the results is verified by mass conservation, and the influences of different electric field directions and different voltages on the movement of bubbles in liquid are considered. The results show that the ratio of the length to axis is proportional to the strength of the electric field when the air bubble is stretched into an ellipsoid along the electric field line under the action of electrostatic gravity and surface tension. In addition, the bubble rising speed is affected by the electric field, the vertical electric field accelerates the bubble rise, and the horizontal direction slows it down.

Numerical simulation on temperature field for resistance spot welding of non-equal thickness stainless steel

An axisymmetric finite element model is developed to simulate the temperature field of resistant spot welding according to the process characters of nugget formation of non equal stainless steel sheets. A simulation method of the interaction of electrical and thermal factors is presented. The spot welding process of nugget formation is simulated using hard and soft welding technique norms. The heating characters of soft and hard norms determine the differences in the process of nugget formation and determine the finally shape and offset of nugget. Experimental verification shows that the model prediction agrees well with the practical.