A MIXED FINITE ELEMENT AND CHARACTERISTIC MIXED FINITE ELEMENT FOR INCOMPRESSIBLE MISCIBLE DARCY-FORCHHEIMER DISPLACEMENT AND NUMERICAL ANALYSIS

In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order.The concentration equation is solved by the method of mixed finite element,where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method.The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation.In actual computations the characteristics adopts a large time step without any loss of accuracy.The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux.In order to derive the optimal 3/2-order error estimate in L^(2) norm,a post-processing technique is included in the approximation to the scalar unknowns.Numerical experiments are illustrated finally to validate theoretical analysis and efficiency.This method can be used to solve such an important problem.

Dynamic Modelling of a Hybrid Variable Reluctance Machine Using the 3D Finite Element Method

This paper presents the work carried out to evaluate the dynamic performance of the Hybrid Variable Reluctance Motor (HVRM). The fourth-order Runge-Kutta integration algorithm was employed to solve the equations of the dynamic model, in conjunction with the three-dimensional finite element method. The 3D numerical data was calculated using Comsol Multiphysics, which accounts for the nonlinearity of the ferromagnetic material and the 3D nature of the HVRM. The outcomes of this study are precise and accurately predict the dynamic behaviour of the HVRM in terms of rotor position response, rotational speed and torque. The distinctive contribution of this work lies in the 3D numerical modelling of the HVRM and the subsequent evaluation and analysis of its dynamic operation. Analytical and numerical 2D studies are less resource-intensive and time-consuming, and are more straightforward and rapid to analyse and interpret. However, they are constrained in their capacity to examine spatial, volumetric interactions and intricate dynamics such as flux studies where three 3D effects cannot be disregarded, winding end effects and the configuration and positioning of the interposed permanent magnet.

Numerical modeling dynamic process of multi-feed microwave heating of industrial solution media

The exothermic efficiency of microwave heating an electrolyte/water solution is remarkably high due to the dielectric heating by orientation polarization of water and resistance heating by the Joule process occurred simultaneously compared with pure water.A three-dimensional finite element numerical model of multi-feed microwave heating industrial liquids continuously flowing in a meter-scale circular tube is presented.The temperature field inside the applicator tube in the cavity is solved by COMSOL Multiphysics and professional programming to describe the momentum,energy and Maxwell's equations.The evaluations of the electromagnetic field,the temperature distribution and the velocity field are simulated for the fluids dynamically heated by singleand multi-feed microwave system,respectively.Both the pilot experimental investigations and numerical results of microwave with single-feed heating for fluids with different effective permittivity and flow rates show that the presented numerical modeling makes it possible to analyze dynamic process of multi-feed microwave heating the industrial liquid.The study aids in enhancing the understanding and optimizing of dynamic process in the use of multi-feed microwave heating industrial continuous flow for a variety of material properties and technical parameters.

Numerical Simulation of Flow Field in Tidal Estuary Using Body-Fitted Coordinate System

On the basis of the physical mechanism, a body-fitted coordinate system is developed. By using this system the boundaries in simulation and in real are fitted well, and simulation with great accuracy is achieved. A computation example indicates that compared to traditional two-dimensional computation methods, the body-fitted simulation has an advantange of better coincidence with the real and can be adopted in simulating flow fields in tidal estuaries.

Numerical Analysis of Magnetic Fluid Sealing Performance

The sealing performance of magnetic fluid is related to the magnetic fluid itself. Many factors can influence the magnetic field and the seal pressure differences of magnetic fluid seals, such as the sealing gap, the shaft eccentricity, the shaft diameter, the volume of the magnetic fluid and the centrifugal force. These factors are analyzed by numerical computation. When the seal material and structure are the same, the seal pressure difference is directly proportional to the magnetic field intensity and the saturation magnetization of the magnetic fluid. The sealing performance of the magnetic fluid will reduce with the increase of the sealing gap and shaft eccentricity. The sealing performance will increase with the volume of the magnetic fluid and decrease with the increase of the shaft diameter taking gravity into account. The increase of the shaft diameter is the same as the reduction of the volume of the magnetic fluid. The magnetic fluid cross-section can change because of the centrifugal force. Some improvements can reducc the influence of the centrifugal force. The centrifugal force can be utilized to improve the sealing performance.

Sheltering effect of punched steel plate sand fences for controlling blown sand hazards along the Golmud-Korla Railway:Field observation and numerical simulation studies

Sand fences made of punched steel plate(PSP)have recently been applied to control wind-blown sand in desertified and Gobi areas due to their strong wind resistance and convenient in situ construction.However,few studies have assessed the protective effect of PSP sand fences,especially through field observations.This study analyzes the effects of double-row PSP sand fences on wind and sand resistance using field observations and a computational fluid dynamics(CFD)numerical simulation.The results of field observations showed that the average windproof efficiencies of the first-row and second-row sand fences were 79.8%and 70.8%,respectively.Moreover,the average windproof efficiencies of the numerical simulation behind the first-row and second-row sand fences were 89.8%and 81.1%,respectively.The sand-resistance efficiency of the double-row PSP sand fences was 65.4%.Sand deposition occurred close to the first-row sand fence;however,there was relatively little sand on the leeward side of the second-row sand fence.The length of sand accumulation near PSP sand fences obtained by numerical simulation was basically consistent with that through field observations,indicating that field observations combined with numerical simulation can provide insight into the complex wind-blown sand field over PSP sand fences.This study indicates that the protection efficiency of the double-row PSP sand fences is sufficient for effective control of sand hazards associated with extremely strong wind in the Gobi areas.The output of this work is expected to improve the future application of PSP sand fences.

Numerical Computation in Airflow Field of Compact Spinning with Pneumatic Groove

In compact spinning with pneumatic groove,the computational fluid dynamic model,computed with parallel technologies & Fluent 6.3,is developed to simulate the flow field in condensing zone with 3D computational fluid dynamic (CFD) technology.Flowing state,distribution rules of static pressure,and velocity in condensing zone are characterized and analyzed.The results show that the fiber bundle in compact spinning with pneumatic groove is compacted by airflow and the shape of the pneumatic groove,and the static pressure in condensing zone is negative,as well as the velocity of airflow in condensing zone is not zero.The fluctuation frequencies of the static pressure and velocity near the bottom of the pneumatic groove are relatively higher,and the number of the fluctuation is equal to that of the round holes in condensing zone.

Efficiency analysis of numerical integrations for finite element substructure in real-time hybrid simulation

Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time（TET） decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method（CDM）, the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ（λ = 4） method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.

Pulse cleaning flow models and numerical computation of candle ceramic filters

Analytical and numerical computed models are developed for reverse pulse cleaning system of candle ceramic filters. A standard turbulent model is demonstrated suitably to the designing computation of reverse pulse cleaning system from the experimental and one dimensional computational result. The computed results can be used to guide the designing of reverse pulse cleaning system, which is optimum Venturi geometry. From the computed results, the general conclusions and the designing methods are obtained.

Ultrasonic prediction of crack density using machine learning:A numerical investigation

Cracks are accounted as the most destructive discontinuity in rock, soil, and concrete. Enhancing our knowledge from their properties such as crack distribution, density, and/or aspect ratio is crucial in geo-systems. The most well-known mechanical parameter for such an evaluation is wave velocity through which one can qualitatively or quantitatively characterize the porous media. In small scales, such information is obtained using the ultrasonic pulse velocity(UPV) technique as a non-destructive test. In large-scale geo-systems, however, it is inverted from seismic data. In this paper, we take advantage of the recent advancements in machine learning(ML) for analyzing wave signals and predict rock properties such as crack density(CD) – the number of cracks per unit volume. To this end, we designed numerical models with different CDs and, using the rotated staggered finite-difference grid(RSG) technique, simulated wave propagation. Two ML networks, namely Convolutional Neural Networks(CNN) and Long Short-Term Memory(LSTM), are then used to predict CD values. Results show that, by selecting an optimum value for wavelength to crack length ratio, the accuracy of predictions of test data can reach R2> 96% with mean square error(MSE) < 25e-4(normalized values). Overall, we found that:(i) performance of both CNN and LSTM is highly promising,(ii) accuracy of the transmitted signals is slightly higher than the reflected signals,(iii) accuracy of 2D signals is marginally higher than 1D signals,(iv)accuracy of horizontal and vertical component signals are comparable,(v) accuracy of coda signals is less when the whole signals are used. Our results, thus, reveal that the ML methods can provide rapid solutions and estimations for crack density, without the necessity of further modeling.

Wind Tunnel Test and Numerical Computation on Ice Accretion on Blade Airfoil for Straight-bladed VAWT

To invest the condition of ice accretion on the blade used for straight-bladed vertical axis wind turbine （SB-VAWT）, wind tunnel tests were carried out on a blade with NACA0015 airfoil by using a small simple icing wind tunnel. Tests were carried out at some typical attack angles under different wind speeds and flow discharges of a water spray with wind. The icing shape and area on blade surface were recorded and measured, Then the numerical computation was carded out to calculate the lift and drag coefficients of the blade before and after ice accretion according to the experiment result, the effect of icing on the aerodynamic characteristics of blade were discussed.

Numerical Computation of Figure-eight Solution for 3-body Problems

The main goal of this paper is to compute the Figure-eight solutions for the planar Newtonian 3-body problem with equal masses by finding the critical points of the functional associated with the motion equations of 3-body in plane R2. The algorithm adopted here is the steepest descent method, which is simple but very valid for our problem.

Design of a Computational Heuristic to Solve the Nonlinear Liénard Differential Model

In this study,the design of a computational heuristic based on the nonlinear Liénard model is presented using the efficiency of artificial neural networks(ANNs)along with the hybridization procedures of global and local search approaches.The global search genetic algorithm(GA)and local search sequential quadratic programming scheme(SQPS)are implemented to solve the nonlinear Liénard model.An objective function using the differential model and boundary conditions is designed and optimized by the hybrid computing strength of the GA-SQPS.The motivation of the ANN procedures along with GA-SQPS comes to present reliable,feasible and precise frameworks to tackle stiff and highly nonlinear differentialmodels.The designed procedures of ANNs along with GA-SQPS are applied for three highly nonlinear differential models.The achieved numerical outcomes on multiple trials using the designed procedures are compared to authenticate the correctness,viability and efficacy.Moreover,statistical performances based on different measures are also provided to check the reliability of the ANN along with GASQPS.

Computation and Validation of Parameter Effects on Lobed Mixer-Ejector Performances

Three-dimensional numerical computation of the flow fields and pumping performances for the lobed mixer-ejector are conducted using full Navier-Stokes equations. In the computation, the inlet of the primary flow uses the mass flowrate boundary condition. The inlet of the second flow and the outlet of the mixing flow use the pressure boundary condition. Compared with the relative experimental resuits, it is shown that the present calculation is reasonable. And a series of numerical studies is performed to obtain the effects of area ratio and length-to-diameter ratio of mixing duct on pumping coefficient and thermal mixing efficiency of a lobed mixer-ejector.

Simulation of moving boundary of the reaction shaft in a flash smelting furnace

As to solve the online monitoring of the inner temperatur e and freezing profile of the reaction shaft of flash smelting furnace, simulation of the wall in the reaction shaft in a flash smelting furnace was made through numerical computation. The computational results are very near the data got in s ite. The error of the moving boundary is approximately 3%, and that of the tempe rature is less than 5%. It is proved that the simulation software is applicable for practice to monitor the temperature and moving boundary inside the hearth on line. Based on a large number of the data computed, the relation between the cha nge of the moving boundary and inner temperature is summarized, and the great in fluence of the cooling system on the forming and stability of the moving boundar y inside the hearth is emphasized, which provide the theoretical bases for optim izing the flash smelting operation.

Multiply accumulate operations in memristor crossbar arrays for analog computing

Memristors are now becoming a prominent candidate to serve as the building blocks of non-von Neumann inmemory computing architectures.By mapping analog numerical matrices into memristor crossbar arrays,efficient multiply accumulate operations can be performed in a massively parallel fashion using the physics mechanisms of Ohm’s law and Kirchhoff’s law.In this brief review,we present the recent progress in two niche applications:neural network accelerators and numerical computing units,mainly focusing on the advances in hardware demonstrations.The former one is regarded as soft computing since it can tolerant some degree of the device and array imperfections.The acceleration of multiple layer perceptrons,convolutional neural networks,generative adversarial networks,and long short-term memory neural networks are described.The latter one is hard computing because the solving of numerical problems requires high-precision devices.Several breakthroughs in memristive equation solvers with improved computation accuracies are highlighted.Besides,other nonvolatile devices with the capability of analog computing are also briefly introduced.Finally,we conclude the review with discussions on the challenges and opportunities for future research toward realizing memristive analog computing machines.

Research and design of the main equipments and structure of Xiangjiaba shiplift

The type of pinion and rack vertical shiplifts has been developed in recent a couple of years in the construction of dams.But the design methods and methodologies have rarely been discussed in literature.The Xiangjiaba shiplift is the second shiplift of this type following the Three Gorges shiplift.Being aimed at the technological rationality of the design in synthetically considering security,economy and applicability,this paper presents the research results of some vital issues relating the design of the Xiangjiaba shiplift,including the determination of design water depth of ship chamber based on fluid numeral computation and physical model test,the optimum design of general layout of main equipments and the civil structure of the Xiangjiaba shiplift,the finite element method(FEM) analysis of stress,vibration modes and the buckling of ship chamber,antiseismic research and the design of structures and mechanisms of the shiplift and the optimum design of driving mechanisms.This research provides the theoretical basis for the design of the Xiangjiaba shiplift.The design principles and research methods are valuable for the design of the same type of shiplifts.

Critical condition for the transformation from Taylor cone to cone-jet

An energy method is proposed to investigate the critical transformation condition from a Taylor cone to a cone-jet. Based on the kinetic theorem, the system power allocation and the electrohydrodynamics stability are discussed. The numerical results indicate that the energy of the liquid cone tip experiences a maximum value during the transformation. With the proposed jetting energy, we give the critical transformation condition under which the derivative of jetting energy with respect to the surface area is greater than or equal to the energy required to form a unit of new liquid surface.

RESEARCH ON STABILITY OF MOVING JET CONTAINING DENSE SUSPENDED SOLID PARTICLES

The spatial stability equation of moving jet containing dense suspended solid particles was derived out by means of the continuum phase-coupled model. The stability curves of moving jet far different downstream distances, Reynolds number of flow-field, particle properties and velocities of jetting device are got by the finite difference method based on the asymptotic method and the Eulerian conservative difference scheme. Founded on the analysis of the obtained stability curves it is found that the positive velocity of jetting device widens the unstable frequency range of flow-field hut the effect of the negative one is contrary. In addition, particles existing in the flow-field curb the instability of flow-field and the effect enhances with the decrease of Reynolds number of flow-field. These conclusions benefit learning the development of moving two-phase jet.

TRANSIENT PRESSURE OF PERCOLATION THROUGH ONE DIMENSION POROUS MEDIA WITH THRESHOLD PRESSURE GRADIENT

This paper studies the transient pressure of percolation during one production and one shutting in one dimension porous media with threshold pressure gradient, the differential equations are derived and solved with numerical computation. Basing on numerical solution, it is analyzed that: 1. the relation between the steady pressure at well bore (or at endpoint) and threshold pressure gradient, shut-in time, and the corresponding formulae are derived; 2, the regulation of transient pressure peak. The result is very useful and will help experiments and applications in the development of low permeability reservoirs with threshold pressure gradient.