Search results

The purpose of this paper is to simulate two macrosegregation benchmarks with a newly developed stabilized finite element algorithm based on a semi-implicit pressure correction scheme.

A streamline-upwind/Petrov–Galerkin (SUPG) stabilized finite element algorithm is developed for the coupled conservation equations of mass, momentum, energy and species. A semi-implicit pressure correction method combined with SUPG stabilization technique is proposed to solve the convection flow during solidification. An analytically derived enthalpy method is adopted to solve the energy conservation equation. The nonlinearities of the energy and species equations are tackled by Newton–Raphson method. Two macrosegregation benchmarks considering the solidification of an Al-4.5 per cent Cu alloy and a Sn-10 per cent Pb alloy are simulated.

A very good agreement is achieved by comparison with the classical finite volume method and a novel meshless method for the Al-4.5 per cent Cu alloy solidification benchmark. Moreover, a unique reference numerical solution has been obtained. Besides, it is demonstrated that the stabilized finite element algorithm can capture the flow instability and channel segregation during solidification of the Sn-10 per cent Pb alloy.

A semi-implicit pressure correction method combined with SUPG stabilization technique is adopted to develop robust stabilized finite element algorithm for the macrosegregation model. A new enthalpy formulation for heat transfer problems with phase change is derived analytically.

The purpose of this paper is to investigate the feasibility of solving turbulent flows based on smoothed finite element method (S-FEM). Then, the differences between S-FEM and finite element method (FEM) in dealing with turbulent flows are compared.

The stabilization scheme, the streamline-upwind/Petrov-Galerkin stabilization is coupled with stabilized pressure gradient projection in the fractional step framework. The Reynolds-averaged Navier-Stokes equations with standard k-epsilon model are selected to solve turbulent flows based on S-FEM and FEM. Standard wall functions are applied to predict boundary layer profiles.

This paper explores a completely new application of S-FEM on turbulent flows. The adopted stabilization scheme presents a good performance on stabilizing the flows, especially for very high Reynolds numbers flows. An advantage of S-FEM is found in applying wall functions comparing with FEM. The differences between S-FEM and FEM have been investigated.

The research in this work is limited to the two-dimensional incompressible turbulent flow.

The verification and validation of a new combination are conducted by several numerical examples. The new combination could be used to deal with more complicated turbulent flows.

The applications of the new combination to study basic and complex turbulent flow are also presented, which demonstrates its potential to solve more turbulent flows in nature and engineering.

This work carries out a great extension of S-FEM in simulations of fluid dynamics. The new combination is verified to be very effective in handling turbulent flows. The performances of S-FEM and FEM on turbulent flows were analyzed by several numerical examples. Superior results were found compared with existing results and experiments. Meanwhile, S-FEM has an advantage of accuracy in predicting boundary layer profile.

The purpose of this paper is to reveal the solute segregation behavior in the molten and solidified regions during direct chill (DC) casting of a large-size magnesium alloy slab under no magnetic field (NMF), harmonic magnetic field (HMF), pulsed magnetic field (PMF) and two types of out-of-phase pulsed magnetic field (OPMF).

A 3-D multiphysical coupling mathematical model is used to evaluate the corresponding physical fields. The coupling issue is addressed using the method of separating step and result inheritance.

The results suggest that the solute deficiency tends to occur in the central part, while the solute-enriched area appears near the fillet in the molten and solidified regions. Applying magnetic field could greatly homogenize the solute field in the two-phase region. The variance of relative segregation level in the solidified cross-section under NMF is 38.1%, while it is 21.9%, 18.6%, 16.4% and 12.4% under OPMF2 (the current phase in the upper coil is ahead of the lower coil), HMF, PMF and OPMF1 (the current phase in the upper coil lags behind the lower coil), respectively, indicating that OPMF1 is more effective to reduce the macrosegregation level.

There are few reports on the solute segregation degree in rectangle slab under magnetic field, especially for magnesium alloy slab. This paper can act a reference to make clear the solute transport behavior and help reduce the macrosegregation level during DC casting.

We introduce a second‐order accurate time‐marching pressure‐correction algorithm to accommodate weakly‐compressible highly‐viscous liquid flows at low Mach number. As the incompressible limit is approached (Ma ≈ 0), the consistency of the compressible scheme is highlighted in recovering equivalent incompressible solutions. In the viscous‐dominated regime of low Reynolds number (zone of interest), the algorithm treats the viscous part of the equations in a semi‐implicit form. Two discrete representations are proposed to interpolate density: a piecewise‐constant form with gradient recovery and a linear interpolation form, akin to that on pressure. Numerical performance is considered on a number of classical benchmark problems for highly viscous liquid flows to highlight consistency, accuracy and stability properties. Validation bears out the high quality of performance of both compressible flow implementations, at low to vanishing Mach number. Neither linear nor constant density interpolations schemes degrade the second‐order accuracy of the original incompressible fractional‐staged pressure‐correction scheme. The piecewise‐constant interpolation scheme is advocated as a viable method of choice, with its advantages of order retention, yet efficiency in implementation.

Several constitutive models are available in the literature to describe the mechanical behaviour of cement stabilized soils. However, difficulties in implementing such models within commercial finite element programs have hindered their application to solve related boundary value problems. Therefore, the aim of this study is to implement a constitutive model, which has the capability to simulate cement stabilized soil behaviour, into the finite element program ABAQUS through the user material subroutine UMAT.

After a detailed review of existing constitutive models for cement stabilized soils, a model based on the elasto‐plastic theory and the extended critical state concept with an associated flow rule is selected for the finite element implementation. A semi‐implicit integration method (cutting plane algorithm) with a continuum elasto‐plastic modulus and path dependent stress prediction strategy has been used in the implementation. The performance of the new finite element formulation of the constitutive model is verified by simulating triaxial test data using the finite element program with the new implementation and predictions from constitutive equations as well as experimental data.

The paper provides the implementation procedure of the constitutive model into ABAQUS but this method is useful for the implementation of any other constitutive model into ABAQUS or any other finite element program. Simulated results for the volumetric deformation of cement stabilized soils show that the cement stabilized soils do not obey the associated flow rule at high confining pressures. The parametric study shows that the influence of cementation increases the brittle nature and the bearing capacity of treated clay. In addition the results show that proposed finite element implementation has the ability to illustrate key features of the cement stabilized clay.

This paper presents an implementation of an elasto‐plastic constitutive model, based on the extended critical state concept, for cement stabilized soils into a finite element programme, which has been identified as an important and challenging topic in computational geomechanics. This implementation is useful in solving boundary value problems in geomechanics involving cement stabilized soils, incorporating key characteristics of these soils.

This paper aims to propose a characteristic stabilized finite element method for non-stationary conduction-convection problems.

To avoid difficulty caused by the trilinear term, the authors use the characteristic method to deal with the time derivative term and the advection term. The space discretization adopts the low-order triples (i.e. P₁-P₁-P₁ and P₁-P₀-P₁ triples). As low-order triples do not satisfy inf-sup condition, the authors use the stability technique to overcome this flaw.

The stability and the convergence analysis shows that the method is stable and has optimal-order error estimates.

Numerical experiments confirm the theoretical analysis and illustrate that the authors’ method is highly effective and reliable, and consumes less CPU time.

This paper aims to propose two parallel two-step finite element algorithms based on fully overlapping domain decomposition for solving the 2D/3D time-dependent natural convection problem.

The first-order implicit Euler formula and second-order Crank–Nicolson formula are used to time discretization respectively. Each processor of the algorithms computes a stabilized solution in its own global composite mesh in parallel. These algorithms compute a nonlinear system for the velocity, pressure and temperature based on a lower-order element pair (P₁b-P₁-P₁) and solve a linear approximation based on a higher-order element pair (P₂-P₁-P₂) on the same mesh, which shows that the new algorithms have the same convergence rate as the two-step finite element methods. What is more, the stability analysis of the proposed algorithms is derived. Finally, numerical experiments are presented to demonstrate the efficacy and accuracy of the proposed algorithms.

Finally, numerical experiments are presented to demonstrate the efficacy and accuracy of the proposed algorithms.

The novel parallel two-step algorithms for incompressible natural convection problem are proposed. The rigorous analysis of the stability is given for the proposed parallel two-step algorithms. Extensive 2D/3D numerical tests demonstrate that the parallel two-step algorithms can deal with the incompressible natural convection problem for high Rayleigh number well.

The purpose of this paper is to consider double‐diffusive convection in a heated porous medium saturated with a fluid. Of particular interest is the case where the fluid has a stabilizing concentration gradient and small diffusivity.

A fully‐coupled stabilized finite element scheme and adaptive mesh refinement (AMR) methodology are introduced to solve the resulting coupled multiphysics application and resolve fine scale solution features. The code is written on top of the open source finite element library LibMesh, and is suitable for parallel, high‐performance simulations of large‐scale problems.

The stabilized adaptive finite element scheme is used to compute steady and unsteady onset of convection in a generalized Horton‐Rogers‐Lapwood problem in both two and three‐dimensional domains. A detailed study confirming the applicability of AMR in obtaining the predicted dependence of solutal Nusselt number on Lewis number is given. A semi‐permeable barrier version of the generalized HRL problem is also studied and is believed to present an interesting benchmark for AMR codes owing to the different boundary and internal layers present in the problem. Finally, some representative adaptive results in a complex 3D heated‐pipe geometry are presented.

This work demonstrates the feasibility of stabilized, adaptive finite element schemes for computing simple double‐diffusive flow models, and it represents an easily‐generalizable starting point for more complex calculations since it is based on a highly‐general finite element library. The complementary nature of h‐adaptivity and stabilized finite element techniques for this class of problem is demonstrated using particularly simple error indicators and stabilization parameters. Finally, an interesting double‐diffusive convection benchmark problem having a semi‐permeable barrier is suggested.

The purpose of this paper is to describe a finite element formulation to approximate thermally coupled flows using both the Boussinesq and the low Mach number models with particular emphasis on the numerical implementation of the algorithm developed.

The formulation, that allows us to consider convection dominated problems using equal order interpolation for all the valuables of the problem, is based on the subgrid scale concept. The full Newton linearization strategy gives rise to monolithic treatment of the coupling of variables whereas some fixed point schemes permit the segregated treatment of velocity‐pressure and temperature. A relaxation scheme based on the Armijo rule has also been developed.

A full Newtown linearization turns out to be very efficient for steady‐state problems and very robust when it is combined with a line search strategy. A segregated treatment of velocity‐pressure and temperature happens to be more appropriate for transient problems.

A fractional step scheme, splitting also momentum and continuity equations, could be further analysed.

The results presented in the paper are useful to decide the solution strategy for a given problem.

The numerical implementation of a stabilized finite element approximation of thermally coupled flows is described. The implementation algorithm is developed considering several possibilities for the solution of the discrete nonlinear problem.

Introduces papers from this area of expertise from the ISEF 1999 Proceedings. States the goal herein is one of identifying devices or systems able to provide prescribed performance. Notes that 18 papers from the Symposium are grouped in the area of automated optimal design. Describes the main challenges that condition computational electromagnetism’s future development. Concludes by itemizing the range of applications from small activators to optimization of induction heating systems in this third chapter.