Search results

The purpose of this paper is to develop an effective numerical algorithm for a gas-melt two-phase flow and use it to simulate a polymer melt filling process. Moreover, the suggested algorithm can deal with the moving interface and discontinuities of unknowns across the interface.

The algebraic sub-grid scales-variational multi-scale (ASGS-VMS) finite element method is used to solve the polymer melt filling process. Meanwhile, the time is discretized using the Crank–Nicolson-based split fractional step algorithm to reduce the computational time. The improved level set method is used to capture the melt front interface, and the related equations are discretized by the second-order Taylor–Galerkin scheme in space and the third-order total variation diminishing Runge–Kutta scheme in time.

The numerical method is validated by the benchmark problem. Moreover, the viscoelastic polymer melt filling process is investigated in a rectangular cavity. The front interface, pressure field and flow-induced stresses of polymer melt during the filling process are predicted. Overall, this paper presents a VMS method for polymer injection molding. The present numerical method is extremely suitable for two free surface problems.

For the first time ever, the ASGS-VMS finite element method is performed for the two-phase flow of polymer melt filling process, and an effective numerical method is designed to catch the moving surface.

The purpose of this paper is to apply the variational multi-scale framework to the finite element approximation of the compressible Navier-Stokes equations written in conservation form. Even though this formulation is relatively well known, some particular features that have been applied with great success in other flow problems are incorporated.

The orthogonal subgrid scales, the non-linear tracking of these subscales, and their time evolution are applied. Moreover, a systematic way to design the matrix of algorithmic parameters from the perspective of a Fourier analysis is given, and the adjoint of the non-linear operator including the volumetric part of the convective term is defined. Because the subgrid stabilization method works in the streamline direction, an anisotropic shock capturing method that keeps the diffusion unaltered in the direction of the streamlines, but modifies the crosswind diffusion is implemented. The artificial shock capturing diffusivity is calculated by using the orthogonal projection onto the finite element space of the gradient of the solution, instead of the common residual definition. Temporal derivatives are integrated in an explicit fashion.

Subsonic and supersonic numerical experiments show that including the orthogonal, dynamic, and the non-linear subscales improve the accuracy of the compressible formulation. The non-linearity introduced by the anisotropic shock capturing method has less effect in the convergence behavior to the steady state.

A complete investigation of the stabilized formulation of the compressible problem is addressed.

This paper aims to present a finite element formulation to approximate systems of reaction–diffusion–advection equations, focusing on cases with nonlinear reaction. The formulation is based on the orthogonal sub-grid scale approach, with some simplifications that allow one to stabilize only the convective term, which is the source of potential instabilities. The space approximation is combined with finite difference time integration and a Newton–Raphson linearization of the reactive term. Some numerical examples show the accuracy of the resulting formulation. Applications using classical nonlinear reaction models in population dynamics are also provided, showing the robustness of the approach proposed.

A stabilized finite element method for advection–diffusion–reaction equations to the problem on nonlinear reaction is adapted. The formulation designed has been implemented in a computer code. Numerical examples are run to show the accuracy and robustness of the formulation.

The stabilized finite element method from which the authors depart can be adapted to problems with nonlinear reaction. The resulting method is very robust and accurate. The framework developed is applicable to several problems of interest by themselves, such as the predator–prey model.

A stabilized finite element method to problems with nonlinear reaction has been extended. Original contributions are the design of the stabilization parameters and the linearization of the problem. The application examples, apart from demonstrating the validity of the numerical model, help to get insight in the system of nonlinear equations being solved.

This paper aims to develop a robust set of advanced numerical tools to simulate multiphase flows under the superimposition of external uniform magnetic fields.

The flow has been simulated in a fully Eulerian framework by a {\it variational multi-scale} method, which allows to take into account the small-scale turbulence without explicitly model it. The multi-fluid problem has been solved through the convectively re-initialized level-set method to robustly deal with high density and viscosity ratio between the phases and the surface tension has been modelled implicitly in the level-set framework. The interaction with the magnetic field has been modelled through the classic induction equation for 2D problems and the time step computation is based on the electromagnetic interaction to guarantee convergence of the method. Anisotropic mesh adaptation is then used to adapt the mesh to the main problem’s variables and to reach good accuracy with a small number of degrees of freedom. Finally, the variational multiscale method leads to a natural stabilization of the finite elements algorithm, preventing numerical spurious oscillations in the solution of Navier–Stokes equations (fluid mechanics) and the transport equation (level-set convection).

The methodology has been validated, and it is shown to produce accurate results also with a low number of degrees of freedom. The physical effect of the external magnetic field on the multiphase flow has been analysed.

The dam-break benchmark case has been extended to include magnetically constrained flows.