Numerical Study of Fluid Dynamics and Heat Transfer Induced by Plasma Discharges

A numerical investigation is conducted to explore the evolution of a plasma discharge and its interaction with the fluid flow based on a self-consistent fluid model which couples the discharge dynamics with the fluid dynamics.The effects of the applied voltage on the distribution of velocity and temperature in initially static air are parainetrically studied.Furthermore,the spatial structure of plasma discharge and the resulting force contours in streamwise and normal directions are discussed in detail.The result shows that the plasma actuator produces a net force that should always be directed away from the exposed electrode,which results in an ionic wind pushing particles into a jet downstream of the actuator.When the energy added by the plasma is taken into account,the ambient air temperature is increased slightly around the electrode,but the velocity is almost not affected.Therefore it is unlikely that the induced flow is buoyancy driven.For the operating voltages considered in this paper,the maximum induced velocity is found to follow a power law,i.e.,it is proportional to the applied voltage to the 3.5 power.This promises an efficient application in the flow control with plasma actuators.

Study on Numerical Simulation of Mold Filling and Heat Transfer in Die Casting Process

A 3-D mathematical model considering turbulence phenomena has been established based on a computational fluid dynamics technique, so called 3-D SOLA-VOF (Solution Algorithm-Volume of Fluid), to simulate the fluid flow of mold filling process of die casting. In addition, the mathematical model for simulating the heat transfer in die casting process has also been established. The computation program has been developed by the authors with the finite difference method (FDM) recently. As verification, the mold filling process of a S-shaped die casting has been simulated and the simulation results coincide with that of the benchmark test. Finally, as a practical application, the gating design of a motorcycle component was modified by the mold filling simulation and the dies design of another motorcycle component was optimized by the heat transfer simulation. All the optimized designs were verified by the production practice.

Numerical Simulation on Subcooled Boiling Heat Transfer Characteristics of Water-Cooled W/Cu Divertors

In order to realize safe and stable operation of a water-cooled W/Cu divertor under high heating condition,the exact knowledge of its subcooled boiling heat transfer characteristics under different design parameters is crucial.In this paper,subcooled boiling heat transfer in a water-cooled W/Cu divertor was numerically investigated based on computational fluid dynamic（CFD）.The boiling heat transfer was simulated based on the Euler homogeneous phase model,and local differences of liquid physical properties were considered under one-sided high heating conditions.The calculated wall temperature was in good agreement with experimental results,with the maximum error of 5%only.On this basis,the void fraction distribution,flow field and heat transfer coefficient（HTC）distribution were obtained.The effects of heat flux,inlet velocity and inlet temperature on temperature distribution and pressure drop of a water-cooled W/Cu divertor were also investigated.These results provide a valuable reference for the thermal-hydraulic design of a water-cooled W/Cu divertor.

NUMERICAL MODELING OF RADIATIVE TRANSFER FOR MICROWAVE REMOTE SENSING

An overall vector radiative transfer theory was developed for numerical modeling, in both active and passive microwave remote sensing. The Theory and approaches are briefly summerized.To quantitatively understand scattering and thermal emission from targets in active and passive remote sensing, we have developed an overall vector radiative transfer theory for a set of theoretical models of discrete scatterer and continuous random media for the earth terrain (wet soil, vegetation, snow, sea-ice, etc.) and atmosphere, and numerical approaches for simulation, data analysis, and parameter sensitivity test. Our numerical results favorably agreed with experimental data in microwave re mote sensing of various earth surfaces. Main approaches are briefly summerized herewith.

Numerical and Experimental Investigation on the In-Flight Melting Behaviour of Granulated Powders in Induction Thermal Plasmas

An innovative in-flight glass melting technology with thermal plasmas was developed for the purpose of energy conservation and environment protection. In this study, modelling and experiments of argon-oxygen induction thermal plasmas were conducted to investigate the melting behaviour of granulated soda-lime glass powders injected into the plasma. A two-dimensional local thermodynamic equilibrium （LTE） model was performed to simulate the heat and momentum transfer between plasma and particle. Results showed that the particle temperature was strongly affected by the flow rate of carrier gas and the particle size of raw material. A higher flow rate of carrier gas led to lower particle temperature and less energy transferred to particles which resulted in lower vitrification. The incomplete melting of large particles was attributed to the lower central temperature of the particle caused by a larger heat capacity. The numerical analysis explained well the experimental results, which can provide valuable practical guidelines for the process control in the melting process for the glass industry.

An Improved Coupling of Numerical and Physical Models for Simulating Wave Propagation

An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model （FEM） based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level.

Real-time 3-D space numerical shake prediction for earthquake early warning

In earthquake early warning systems, real-time shake prediction through wave propagation simulation is a promising approach. Compared with traditional methods, it does not suffer from the inaccurate estimation of source parameters. For computation efficiency, wave direction is assumed to propagate on the 2-D surface of the earth in these methods. In fact, since the seismic wave propagates in the 3-D sphere of the earth, the 2-D space modeling of wave direction results in inaccurate wave estimation. In this paper, we propose a 3-D space numerical shake pre- diction method, which simulates the wave propagation in 3-D space using radiative transfer theory, and incorporate data assimilation technique to estimate the distribution of wave energy. 2011 Tohoku earthquake is studied as an example to show the validity of the proposed model. 2-D space model and 3-D space model are compared in this article, and the prediction results show that numerical shake prediction based on 3-D space model can estimate the real-time ground motion precisely, and overprediction is alleviated when using 3-D space model.

Real-time numerical shake prediction and updating for earthquake early warning

Ground motion prediction is important for earthquake early warning systems, because the region＇s peak ground motion indicates the potential disaster. In order to predict the peak ground motion quickly and pre- cisely with limited station wave records, we propose a real- time numerical shake prediction and updating method. Our method first predicts the ground motion based on the ground motion prediction equation after P waves detection of several stations, denoted as the initial prediction. In order to correct the prediction error of the initial prediction, an updating scheme based on real-time simulation of wave propagation is designed. Data assimilation technique is incorporated to predict the distribution of seismic wave energy precisely. Radiative transfer theory and Monte Carlo simulation are used for modeling wave propagation in 2-D space, and the peak ground motion is calculated as quickly as possible. Our method has potential to predict shakemap, making the potential disaster be predicted before the real disaster happens. 2008 Ms8.0 Wenchuan earthquake is studied as an example to show the validity of the proposed method.

Numerical simulation of flow and heat transfer in a random packed bed

Random packed beds have more complex interior structure than structured beds and are widely used in industry and engineering. CFD simulation was carried out to investigate and analyze the local flow and heat transfer in a 120-sphere random packed bed. 3D Navier-Stokes equation was solved with a finite volume formulation based on the Chimera meshing technique. Investigation was focused on low Reynolds number flow （Re=4.6-56.2）, which typically occurs in packed bed reactors in bio-chemical fields. Detailed temperature field information was obtained. Inhomogeneity of flow and heat transfer due to the non-uniform distribution of void fraction was discussed and analyzed.

A comparison of approaches in fitting continuum SEDs

We present a detailed comparison of two approaches, the use of a precalculated database and simulated annealing （SA）, for fitting the continuum spectral energy distribution （SED） of astrophysical objects whose appearance is dominated by surrounding dust. While pre-calculated databases are commonly used to model SED data, only a few studies to date employed SA due to its unclear accuracy and convergence time for this specific problem. From a methodological point of view, different approaches lead to different fitting quality, demand on computational resources and calculation time. We compare the fitting quality and computational costs of these two approaches for the task of SED fitting to provide a guide to the practitioner to find a compromise between desired accuracy and available resources. To reduce uncertainties inherent to real datasets, we introduce a reference model resembling a typical circumstellar system with 10 free parameters. We derive the SED of the reference model with our code MC3 D at 78 logarithmically distributed wavelengths in the range [0.3 μm, 1.3 mini and use this setup to simulate SEDs for the database and SA. Our result directly demonstrates the applicability of SA in the field of SED modeling, since the algorithm regularly finds better solutions to the optimization problem than a precalculated database. As both methods have advantages and shortcomings, a hybrid approach is preferable. While the database provides an approximate fit and overall probability distributions for all parameters deduced using Bayesian analysis, SA can be used to improve upon the results returned by the model grid.

A Numerical Methodology for Enforcing Maximum Principles and the Non-Negative Constraint for Transient Diffusion Equations

Transient diffusion equations arise in many branches of engineering and applied sciences(e.g.,heat transfer and mass transfer),and are parabolic partial differential equations.It is well-known that these equations satisfy important mathematical properties like maximum principles and the non-negative constraint,which have implications in mathematical modeling.However,existing numerical formulations for these types of equations do not,in general,satisfy maximum principles and the nonnegative constraint.In this paper,we present a methodology for enforcing maximum principles and the non-negative constraint for transient anisotropic diffusion equation.The proposed methodology is based on the method of horizontal lines in which the time is discretized first.This results in solving steady anisotropic diffusion equation with decay equation at every discrete time-level.We also present other plausible temporal discretizations,and illustrate their shortcomings in meeting maximum principles and the non-negative constraint.The proposedmethodology can handle general computational grids with no additional restrictions on the time-step.We illustrate the performance and accuracy of the proposed methodology using representative numerical examples.We also perform a numerical convergence analysis of the proposed methodology.For comparison,we also present the results from the standard singlefield semi-discrete formulation and the results froma popular software package,which all will violate maximum principles and the non-negative constraint.

Numerical Study and Performance Analyses of Mems-Based Particle Velocity Sensor with Combined Inclined Rib Pair

The combined inclined rib pair(CIRP)is the first time proposed to improve the sensor performance of particle velocity sensor(PVS)by using a three-dimensional numerical method.The method is verified by the experiment results in the literature.The optimal plain channel parameters are determined as the basic sensor structure.In comparison of plain channel,both heat transfer characteristics and sensor performance are enhanced effectively by arranging the CIRP.The reason is that the high flow rate region caused by the CIRP can maintain strongly in the whole fluid field if there are enough rib pairs.Furthermore,the produced longitudinal vortex pair can get a better fluid mix,which is more conductive to heat transfer.The increasing height and number of the CIRP can improve the heat transfer characteristics,but the flow resistance will increase as well.For the purpose of finding the best overall performance,the effects of the parameters including the geometric sizes and the position of the CIRP have been investigated.The results show that PVS will get the best sensitivity when the rib length and width are 0.2 mm and 0.03 mm respectively,and the distance between rib pair and between ribs in the same pair are 0.15 mm and 0.3 mm respectively.Besides,the most suitable crossing angle is 45°.Thus,the performance of PVS can be significantly improved by this novel structure.

Simulation of steel beam under ceiling jet based on a wind–fire–structure coupling model

For localized fires, it is necessary to consider the thermal and mechanical responses of building elements subject to uneven heating under the influence of wind. In this paper, the thermomechanical phenomena experienced by a ceiling jet and I-beam in a structural fire were simulated. Instead of applying the concept of adiabatic surface temperature (AST) to achieve fluid–structure coupling, this paper proposes a new computational fluid dynamics–finite element method numerical simulation that combines wind, fire, thermal, and structural analyses. First, to analyze the velocity and temperature distributions, the results of the numerical model and experiment were compared in windless conditions, showing good agreement. Vortices were found in the local area formed by the upper and lower flanges of the I-beam and the web, generating a local high-temperature zone and enhancing the heat transfer of convection. In an incoming-flow scenario, the flame was blown askew significantly;the wall temperature was bimodally distributed in the axial direction. The first temperature peak was mainly caused by radiative heat transfer, while the second resulted from convective heat transfer. In terms of mechanical response, the yield strength degradation in the highest-temperature region in windless conditions was found to be significant, thus explaining the stress distribution of steel beams in the fire field. The mechanical response of the overall elements considering the incoming flows was essentially elastic.

Impacts of Previous Winter Kuroshio SSTA on Summer Rainfall in China

The relationship of Kuroshio sea surface temperature anomaly (SSTA) in theprevious winter and summer rainfall in China was analyzed based on observational studies andnumerical simulations. Observational results indicate that there is a close relation betweenKuroshio SSTA and precipitation in China. When Kuroshio SSTA is positive, the western Pacificsubtropical high will be stronger and extend farther westward in, the following summer, with Asiansummer monsoon weaker and the frontal precipitation further southward. As a result, summerprecipitation increases (decreases) in the mid-lower reaches of the Yangtze River (in North andNortheast China), and vice versa. Conclusions drawn by NCAR-CCM3 testify the observational results.

Multiphase Flow and Thermo-Mechanical Behaviors of Solidifying Shell in Continuous Casting Mold

The metallurgical phenomena occurring in the continuous casting mold have a significant influence on the performance and the quality of steel product.The multiphase flow phenomena of molten steel,steel/slag interface and gas bubbles in the slab continuous casting mold were described by numerical simulation,and the effect of electromagnetic brake（EMBR） and argon gas blowing on the process were investigated.The relationship between wavy fluctuation height near meniscus and the level fluctuation index F,which reflects the situation of mold flux entrapment,was clarified.Moreover,based on a microsegregation model of solute elements in mushy zone with δ/γ transformation and a thermo-mechanical coupling finite element model of shell solidification,the thermal and mechanical behaviors of solidifying shell including the dynamic distribution laws of air gap and mold flux,temperature and stress of shell in slab continuous casting mold were described.

Computational fluid dynamics simulation of friction stir welding:A comparative study on different frictional boundary conditions

Numerical simulation based on computational fluid dynamics （CFD） is a useful approach for quantitatively investigating the underlying thermal-mechanical conditions during FSW, such as temperature field and material deformation field. One of the critical issues in CFD simulation of FSW is the use of the frictional boundary condition, which represents the friction between the welding tool and the workpiece in the numerical models. In this study, three-dimensional numerical simulation is conducted to analyze the heat transfer and plastic deformation behaviors during the FSW of AA2024. For comparison purposes, both the boundary velocity （BV） models and the boundary shear stress （BSS） models are employed in order to assess their performances in predicting the temperature and material deformation in FSW. It is interesting to note that different boundary conditions yield similar predictions on temperature, but quite different predictions on material deformation. The numerical predictions are compared with the experimental results. The predicted deformation zone geometry by the BSS model is consistent with the experimental results while there is large difference between the predictions by the BV models and the experimental measurements. The fact that the BSS model yields more reasonable predictions on the deformation zone geometry is attributed to its capacity to automatically adjust the contact state at the tool/workpiece interface. Based on the favorable predictions on both the temperature field and the material deformation field, the BSS model is suggested to have a better performance in numerical simulation of FSW than the BV model.