An enthalpy-based model of dendritic growth in a convecting binary alloy melt

Purpose-In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues. Design/methodology/approach-The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed. Findings-Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated. Originality/value-As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.

Effect of shrinkage induced flow on binary alloy dendrite growth: An equivalent undercooling model

This paper presents a theoretical model for studying the effects of shrinkage induced flow on the growth rate of binary alloy dendrites. An equivalent undercooling of the melt is defined in terms of ratio of the phase densities to represent the change in dendrite growth rate due to variation in solutal and thermal transport resulting from shrinkage induced flow. Subsequently, results for dendrite growth rate predicted by the equivalent undercooling model is compared with the corresponding predictions obtained using an enthalpy based numerical method for dendrite growth with shrinkage. The agreement is found to be good. Published by Elsevier Ltd.

Comparison of a priori and a posteriori meshes for singularly perturbed nonlinear parameterized problems

This paper deals with the adaptive mesh generation for singularly perturbed nonlinear parameterized problems with a comparative research study on them. We propose an a posteriori error estimate for singularly perturbed parameterized problems by moving mesh methods with fixed number of mesh points. The well known a priori meshes are compared with the proposed one. The comparison results show that the proposed numerical method is highly effective for the generation of layer adapted a posteriori meshes. A numerical experiment of the error behavior on different meshes is carried out to highlight the comparison of the approximated solutions. (C) 2015 Elsevier B.V. All rights reserved.

Evaluation of the Grounding Resistance of Conductors Buried in Steplike Terrain

Availability of land for conventional air-insulated substations is becoming increasingly difficult not only in urban but also in semiurban areas. When the land made available is highly uneven, the associated technoeconomic factors favors the erection of substations on a steplike-formed ground surface and such constructions are in service for more than ten years in some parts of southern India. Noting that the literature on the performance of ground grids in such a construction is rather scarce, the present work was taken up. Evaluation of the performance of earthing elements in steplike ground forms the main goal of the present work. For the numerical evaluation, a suitable boundary-based methodology is employed. This method retains the classical Galerkin approach for the conductors, while the interfaces are replaced by equivalent fictitious surface sources defined over unstructured mesh. Details of the implementation of this numerical method, along with special measures to minimize the computation, are presented. The performance of basic earthing elements, such as the driven rod, counterpoise, and simple grids buried in steplike ground, are analyzed and compared with that for the case with uniform soil surface. It is shown that more than the earthing resistances, the step potentials can get significantly affected.

Vaporization of heated materials into discharge plasmas

The vaporization of condensed materials in contact with high-current discharge plasmas is considered. A kinetic numerical method named direct simulation Monte Carlo (DSMC) and analytical kinetic approaches based on the bimodal distribution function approximation are employed. The solution of the kinetic layer problem depends upon the velocity at the outer boundary of the kinetic layer which varies from very small, corresponding to the high-density plasma near the evaporated surface, up to the sound speed, corresponding to evaporation into vacuum. The heavy particles density and temperature at the kinetic and hydrodynamic layer interface were obtained by the analytical method while DSMC calculation makes it possible to obtain the evolution of the particle distribution function within the kinetic layer and the layer thickness.

Elastic Waves in a Hybrid Multilayered Piezoelectric Plate

An analytical-numerical method is presented for analyzing dispersion and characteristic surface of waves in a hybrid multilayered piezoelectric plate. In this method, the multilayered piezoelectric plate is divided into a number of layered elements with three-nodal-lines in the wall thickness, the coupling between the elastic field and the electric field is considered in each element. The associated frequency dispersion equation is developed and the phase velocity and slowness, as well as the group velocity and slowness are established in terms of the Rayleigh quotient. Six characteristic wave surfaces are introduced to visualize the effects of anisotropy and piezoelectricity on wave propagation. Examples provide a full understanding for the complex phenomena of elastic waves in hybrid multilayered piezoelectric media.

Assessment of the contribution of surface roughness to airframe noise

The generation of sound by turbulent boundary-layer flow at low Mach number over a rough wall is investigated by applying a theoretical model that describes the scattering of the turbulence near field into sound by roughness elements. Attention is focused on the numerical method to approximately quantify the absolute level of far-field radiated roughness noise. Models for the source statistics are obtained by scaling smooth-wall data by the increased skin friction velocity and boundary-layer thickness for a rough surface. Numerical integration is performed to determine the roughness noise, and it reproduces the spectral characteristics of the available empirical formula and experimental data. Experiments are conducted to measure the radiated sound from two rough plates in an open jet The measured noise spectra of the rough plates are above that of a smooth plate in 1-2.5 kHz frequency and exhibit reasonable agreement with the predicted level. Estimates of the roughness noise for a Boeing 757 sized aircraft wing with idealized levels of surface roughness show that hi the high-frequency region the sound radiated from surface roughness may exceed that from the trailing edge, and higher overall sound pressure levels are observed for the roughness noise. The trailing edge noise is also enhanced by surface roughness somewhat A parametric study indicates that roughness height and roughness density significantly affect the roughness noise with roughness height having the dominant effect The roughness noise directivity varies with different levels of surface roughness. Copyright © 2007 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Effects of surface roughness on airframe noise

The generation of sound by turbulent boundary layer flow at low Mach number over a rough wall is investigated by applying the theoretical model which describes the scattering of the turbulence near field into sound by roughness elements. Attention is focused on the numerical method to approximately quantify the absolute level of the roughness noise radiated to far field. Empirical models for the source statistics are obtained by scaling smooth-wall data through increased skin friction velocity and boundary layer thickness for the rough surface. Numerical integration is performed to determine the roughness noise, and it reproduces the spectral characteristics of the available empirical formula and experimental data. Experiments are conducted to measure the radiated sound from two rough plates in an open jet by four 1/2'' free-field condenser microphones. The measured noise spectra of the rough plates are above that of a smooth plate in 1-2.5 kHz frequency and exhibits encouraging agreement with the predicted spectra. Also, a phased microphone array is utilized to localize the sound source, and it confirms that the rough plates generate higher source strengthes in this frequency range. A parametric study illustrates that the roughness height and roughness density significantly affect the far-field radiated roughness noise with the roughness height having the dominant effect. The estimates of the roughness noise for a Boeing 757 sized aircraft wing show that in high frequency region the sound radiated from surface roughness may exceed that from the trailing edge, and higher overall sound pressure levels for the roughness noise are also observed.

Relation between Stabilization Energy, Crystal Field Coefficient and the Magnetic Exchange Interaction for Tb3+ Ion

Based on a single ion model, Hamiltonian of the simplest form about magnetocrystalline anisotropy for Tb3+ ion was solved by using the numerical method. The relation between the stabilization energy, crystal field coefficient B-2(0) and the magnetic exchange interaction was studied as temperature approaches to 0 K. The results show that the stabilization energy contributed by Tb3+ is linear with crystal field coefficient B-2(0) approximately, but it is insensitive to the change of magnetic exchange interaction for the strong magnetic substances such as TbCo5, Tb2Co17 and Tb2Fe14B compounds.

Transition waves and nonlinear interactions in the near wake of a circular cylinder

Transition waves and interactions between two kinds of instability-vortex shedding and transition wave in the near wake of a circular cylinder in the Reynolds number range 3 000-10 000 are studied by a domain decomposition hybrid numerical method. Based on high resolution power spectral analyses for velocity new results on the Reynolds-number dependence of the transition wave frequency, i.e. f(t)/f(s) similar to Re-0.87 are obtained. The new predictions are in good agreement with the experimental results of Wei and Smith but different from Braza's prediction and some early experimental results f(t)/f(s) similar to Re-0.5 given by Bloor et nl. The multi-interactions between two kinds of vortex are clearly visualized numerically. The strong nonlinear interactions between the two independent frequencies (f(t), f(s)) leading to spectra broadening to form the coupling mf(s) +/- nf(t) are predicted and analyzed numerically, and the characteristics of the transition are described. Longitudinal variations of the transition wave and its coupling are reported. Detailed mechanism of the flow transition in the near wake before occurrence of the three-dimensional evolution is provided.

Wave pattern characteristics of a two-phase nozzle flow by shock propagation

In this paper, the wave pattern characteristics of shock-induced two-phase nozzle Hows with the quiescent or moving dusty gas ahead of the incident-shock front has been investigated by using high-resolution numerical method. As compared with the corresponding results in single-phase nozzle flows of the pure gas, obvious differences between these two kinds of flows can be obtained.

The mixed-type integral equations for transient elastodynamic plane crack problems

In the present paper, by use of the boundary integral equation method and the techniques of Green fundamental solution and singularity analysis, the dynamic infinite plane crack problem is investigated. For the first time, the problem is reduced to solving a system of mixed-typed integral equations in Laplace transform domain. The equations consist of ordinary boundary integral equations along the outer boundary and Cauchy singular integral equations along the crack line. The equations obtained are strictly proved to be equivalent with the dual integral equations obtained by Sih in the special case of dynamic Griffith crack problem. The mixed-type integral equations can be solved by combining the numerical method of singular integral equation with the ordinary boundary element method. Further use the numerical method for Laplace transform, several typical examples are calculated and their dynamic stress intensity factors are obtained. The results show that the method proposed is successful and can be used to solve more complicated problems.

A solution to the upward migration of an intrusion

The evolution of the upward migration of the magma is a nonlinear and unstable problem in mathematics. It is difficult to solve it. And using the numerical method, the solution is relatively tedious and time-consuming. This paper introduces a method of the instantaneous point source to solve the linear and unstable heat conduction equation during the infinite period of time instead of the solution of the nonlinear and unstable heat conduction equation. The results obtained by this method coincide with those by the numerical method, meaning that this method offers a simple way to solve the nonlinear and unstable heat conduction equation.

The wave model of mesothermal plasma near wakes and Korteweg-de Vries equation

The stationary two-dimensional (x, z) near wakes behind a flat-based projectile which moves at a constant mesothermal speed (V∞) along a z-axis in a rarefied, fully ionized, plasma is studied using the wave model previously proposed by one of the authors (VCL). One-fluid theory is used to depict the free expansion of ambient plasma into the vacuum produced behind a fast-moving projectile. This nonstationary, one-dimensional (x, t) flow which is approximated by the K-dV equation can be transformed, through substitution, t=z/V∞, into a stationary two-dimensional (x, z) near wake flow seen by an observer moving with the body velocity (V∞). The initial value problem of the K-dV equation in (x, t) variables is solved by a specially devised numerical method. Comparisons of the present numerical solution for the asymptotically small and large times with available analytical solutions are made and found in satisfactory agreements.

地下结构震动响应数值方法研究

随着计算机技术的飞速发展，数值计算已经成为模拟在各种介质中的爆炸过程及其近场波动效应的重要研究手段。针对目前爆炸技术在工程的重要应用课题之一，即地下结构在爆炸载荷作用下的震动响应数值方法进行了较深入的研究，研究了适用于松动爆破的爆破震动的等效模型，该模型简单可行，基本反映了爆破震动过程的物理特征。利用该模型对爆炸载荷作用下地下结构的震动过程进行了模拟，分析爆破在介质中激发的瞬态扰动及其传播的特征，如频率特征、振幅变化规律等，通过与实验数据进行对比，解释了近场爆炸现象并为实验提供理论依据。通过场地模拟实验的数值模拟验证了方法的有效性，为后续工作打下了良好的基础。