The space time CE/SE method for solving one-dimensional batch crystallization model with fines dissolution

This article is concerned with the numerical investigation of one-dimensional population balance models for batch crystallization process with fines dissolution.In batch crystallization,dissolution of smaller unwanted nuclei below some critical size is of vital importance as it improves the quality of product.The crystal growth rates for both size-independent and size-dependent cases are considered.A delay in recycle pipe is also included in the model.The space–time conservation element and solution element method,originally derived for non-reacting flows,is used to solve the model.This scheme has already been applied to a range of PDEs,mainly in the area of fluid mechanics.The numerical results are compared with those obtained from the Koren scheme,showing that the proposed scheme is more efficient.

Finite Element Harmonic Solution of the Coupled Rotor-bearing System

Fluid-solid interaction problems have been studied q uite extensively in the past years. Rotor-bearing system is a typical example. Fluid field is changed under the exciting of rotor vibration. On the same ti me, a net force caused by fluid pressure exerts on rotor, which will change roto r vibration. So, the fluid-solid coupled analysis method must be used. Traditionally, numerical difference method was used to solve fluid problems. The coupled fluid-solid equation could not be set up based on the method. It is no t until finite element method was used in fluid dynamics area then can the coupl ed dynamics be researched. Recently many experimental, analytical and numerical studies have been used in the area . But in these investigations, it is a ssumed that the solid vibration could not be influenced by fluid. In the other w ords, the force exerted on solid from fluid was neglected in the papers. So, the models built were some kinds of semi-coupled model only. In this paper, the Galerkin finite-element method, two-dimension vibration equ ation of rigid body and Navier-Stokes equations are used to build a full-coupl ed fluid-solid model in rotor-bearing system. Some assumptions are taken: 1) In fluid equation, the nonlinear terms are relatively small and neglected. 2) The gravity takes no effect on this system. 3) The bearing and the rotor are long. Flow and leakage along the axis is neglec ted. 4) The fluid is a kind of Newtonian incondensable viscous fluid. 5) The rotor is considered to be a rigid body. Using the model established, we calculated all the examples given by paper , results show the error are less than 7%. So the full-coupled model is built c orrectly. Examples are given in the end of the paper. After analyzing the examples, we get some conclusions: 1) In rotor-bearing system, while being taken under two conditions that whether coupled method is taken or not, difference of pressure and vibration amplitude could reach 76% and 120%. Therefore coupled method must be taken to investigate fluid-solid system. 1) Amplitude of fluid pressure can be more or less influenced by rotor unbalance , gap, eccentricity and other factors. 2) By using coupling method, results show that the amplitudes of vibration and p ressure are greater than ignoring the method. It should be paid more attention t o.

First-principles Investigation on Solution Hardening of Interstitial Impurity Elements (O,N) inγ-TiAl

We have calculated the electronic Structures of O-doped and N-doped r-TiAl using the firstprinciples discrete variational method (DVM) with the aim to understand the solution hardening effects of oxygen and nitrogen in r-TiAl. Our combination analysis on the electronic density, density of states (DOS) and the local environment total bond orders (LTBO) will show that, X atom (X is O or N) can strongly bind with its six surrounding atoms via electronic hybridizations of Ti-3d/X-2p and Al-3p/X-2p. As a sequence, there forms a 'hard' cohesive region around the impurity atom. A pinning model based on the calculations is proposed to explain the hardening effects. The consistent results are obtained between the present calculation and formal test experiments.

THE ANALYTICAL SOLUTIONS BASED ON THE CONCEPT OF FINITE ELEMENT METHODS

On the basis of the concept of finite element methods, the rigorous analytical solutions of structural response in terms of the design variables are researched in this paper. The spatial trusses are taken as an example for the solution of the analytical expressions of the explicit displacements which are proved mathematically; then some conclusions are reached that are useful to structural sensitivity analysis and optimization. In the third part of the paper, a generalized geometric programming method is sugguested for the optimal model with the explicit displacement. Finally, the analytical solutions of the displacements of three trusses are given as examples.

The role of solutes in grain refinement of hypoeutectic magnesium and aluminum alloys

The effect of solutes on grain formation has been studied over a century but is still under debating,simply because it is a very complex topic.This article focuses on the effect of dissolved solute on the growth,fragmentation of a dendrite and the resultant grain size.Experimental data on grain size in magnesium and aluminum alloys with various solute concentrations are collected and analyzed using phase diagram variables including Q,P,andΔT.The physical meaning of each phase diagram variable is discussed.Curve fitting of the recently proposed two-parameter models with experimental data suggests that there is a clear correlation between the measured grain size and the solidification range of the alloy over the entire range of the hypoeutectic composition.Such a trend of grain size vs.ΔT is closely related to the grain refining mechanisms operating under relevant experimental conditions.A critical review of the grain refining mechanisms indicates that the mechanisms associated with the dendrite fragmentation are operating under the conditions where convection is substantial in the melt during its solidification.

A multiscale 3D finite element analysis of fluid/solute transport in mechanically loaded bone

The transport of fluid, nutrients, and signaling molecules in the bone lacunar-canalicular system （LCS） is critical for osteocyte survival and function. We have applied the fluorescence recovery after photobleaching （FRAP） approach to quantify load-induced fluid and solute transport in the LCS in situ, but the measurements were limited to cortical regions 30-50 μm underneath the periosteum due to the constrains of laser penetration. With this work, we aimed to expand our understanding of load-induced fluid and solute transport in both trabecular and cortical bone using a multiscaled image-based finite element analysis （FEA） approach. An intact murine tibia was first re-constructed from microCT images into a three-dimensional （3D） linear elastic FEA model, and the matrix deformations at various locations were calculated under axial loading. A segment of the above 3D model was then imported to the biphasic poroelasticity analysis platform （FEBio） to predict load-induced fluid pressure fields, and interstitial solute/fluid flows through LCS in both cortical and trabecular regions. Further, secondary flow effects such as the shear stress and/or drag force acting on osteocytes, the presumed mechano-sensors in bone, were derived using the previously developed ultrastructural model of Brinkman flow in the canaliculi. The material properties assumed in the FEA models were validated against previously obtained strain and FRAP transport data measured on the cortical cortex. Our results demonstrated the feasibility of this computational approach in estimating the fluid flux in the LCS and the cellular stimulation forces （shear and drag forces） for osteocytes in any cortical and trabecular bone locations, allowing further studies of how the activation of osteocytes correlates with in vivo functional bone formation. The study provides a promising platform to reveal potential cellular mechanisms underlying the anabolic power of exercises and physical activities in treating patients with skeletal deficiencies.

A three-dimensional Eulerian method for the numerical simulation of high-velocity impact problems

In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.

A PARAMETER FOR EVALUATING ABILITY OF GRAIN REFINEMENT OF SOLUTE ELEMENTS

The gradient of the equilibrium melting temperature distribution at the advancing liquid-solid interface of liquid side was derived.From which,a parameter,H,used for evaluating the abili- ty of grain refinement of solute elements was proposed.Experiments on Al-based alloys show that the grain size decreases as the value of parameter,H,increases.Such a dependance indi- cates that the parameter H can be used for choosing grain refiner.

Boundary element analysis for elastic and elastoplastic problems of 2D orthotropic media with stress concentration

Both the orthotropy and the stress concentration are common issues in modem structural engineering. This paper introduces the boundary element method （BEM） into the elastic and elastoplastic analyses for 2D orthotropic media with stress concentration. The discretized boundary element formulations are established, and the stress formulae as well as the fundamental solutions are derived in matrix notations. The numerical procedures are proposed to analyze both elastic and elastoplastic problems of 2D orthotropic me- dia with stress concentration. To obtain more precise stress values with fewer elements, the quadratic isoparametric element formulation is adopted in the boundary discretization and numerical procedures. Numerical examples show that there are significant stress concentrations and different elastoplastic behaviors in some orthotropic media, and some of the computational results are compared with other solutions. Good agreements are also observed, which demonstrates the efficiency and reliability of the present BEM in the stress concentration analysis for orthotropic media.

Enforcing the Discrete Maximum Principle for Linear Finite Element Solutions of Second-Order Elliptic Problems

The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems.The preservation of the qualitative characteristics,such as the maximum principle,in discrete model is one of the key requirements.It is well known that standard linear finite element solution does not satisfy maximum principle on general triangular meshes in 2D.In this paper we consider how to enforce discrete maximum principle for linear finite element solutions for the linear second-order self-adjoint elliptic equation.First approach is based on repair technique,which is a posteriori correction of the discrete solution.Second method is based on constrained optimization.Numerical tests that include anisotropic cases demonstrate how our method works for problems for which the standard finite element methods produce numerical solutions that violate the discrete maximum principle.

TEMPORAL AND SPATIAL DISCRETIZATION ON QUASI-3-D GROUNDWATER FINITE ELEMENT MODELLING TO AVOID SPURIOUS OSCILLATION

In this article, the fmite element solution of quasi-three-dimensional （quasi-3-D） groundwater flow was mathematically analyzed. The research shows that the spurious oscillation solution to the Finite Element Model （FEM） is the results choosing the small time step △t or the large element size L and using the non-diagonal storage matrix. The mechanism for this phenomenon is explained by the negative weighting factor of implicit part in the discretized equations. To avoid spurious oscillation solution, the criteria on the selection of △t and L for quasi-3-D groundwater flow simulations were identified. An application example of quasi-3-D groundwater flow simulation was presented to verify the criteria. The results indicate that temporal discretization scale has significant impact on the spurious oscillations in the finite-element solutions, and the spurious oscillations can be avoided in solving practical quasi-3-D groundwater flow problems if the criteria are satisfied.

Modeling the Earth's magnetosphere under the influence of solar wind with due northward IMF by the AMR-CESE-MHD model

We present a newly developed global magnetohydrodynamic(MHD) model to study the responses of the Earth's magnetosphere to the solar wind. The model is established by using the space-time conservation element and solution element(CESE) method in general curvilinear coordinates on a six-component grid system. As a preliminary study, this paper is to present the model's numerical results of the quasi-steady state and the dynamics of the Earth's magnetosphere under steady solar wind flow with due northward interplanetary magnetic field(IMF). The model results are found to be in good agreement with those published by other numerical magnetospheric models.

Effect of Solute Elements and Cooling Rate on Strain in Brittle Temperature Range of Continuously Cast Strand

A micro-segregation model of solute elements in mushy zone with δ/γ transformation during solidification was established based on the regular hexagon transverse cross section of dendrite shape proposed by finite difference method under the non-equilibrium solidification condition. The model was used to calculate the non-equilibrium pseudo binary Fe-C phase diagram and the strain of steels induced by variation of temperature in brittle temperature range. On the basis of the phase diagram and the strain, the strain curve in brittle temperature range as a function of carbon content for continuously fast strand was introduced and obtained, Solute elements change the position of the strain curve. And cooling rate changes the position arid the shape of the strain curve. The comprehensive formula of the strain as functions of solute elements and cooling rate in brittle temperature range has been obtained by nonlinear fitting program.