Search results

The purpose of this study is to examine theoretically the axisymmetric flow of a steady free-surface jet emerging from a tube for high inertia flow and moderate surface tension effect.

The method of matched asymptotic expansion is used to explore the rich dynamics near the exit where a stress singularity occurs. A boundary layer approach is also proposed to capture the flow further downstream where the free surface layer has grown significantly.

The jet is found to always contract near the tube exit. In contrast to existing numerical studies, the author explores the strength of upstream influence and the flow in the wall layer, resulting from jet contraction. This influence becomes particularly evident from the nonlinear pressure dependence on the upstream distance, as well as the pressure undershoot and overshoot at the exit for weak and strong gravity levels, respectively. The approach is validated against existing experimental and numerical data for the jet profile and centerline velocity where good agreement is obtained. Far from the exit, the author shows how the solution in the diffusive region can be matched to the inviscid far solution, providing the desired appropriate initial condition for the inviscid far flow solution. The location, at which the velocity becomes uniform across the jet, depends strongly on the gravity level and exhibits a non-monotonic behavior with respect to gravity and applied pressure gradient. The author finds that under weak gravity, surface tension has little influence on the final jet radius. The work is a crucial supplement to the existing numerical literature.

Given the presence of the stress singularity at the exit, the work constitutes a superior alternative to a computational approach where the singularity is typically and inaccurately smoothed over. In contrast, in the present study, the singularity is entirely circumvented. Moreover, the flow details are better elucidated, and the various scales involved in different regions are better identified.

A hybrid spectral/boundary element approach is proposed to examine the influence of Couette channel flow on transient coating of highly elastic fluids. The viscoelastic instability of one‐dimensional plane Couette flow is first determined for a large class of Oldroyd fluids with added viscosity, which typically represent polymer solutions composed of a Newtonian solvent and a polymeric solute. The Johnson‐Segalman equation is used as the constitutive model. The velocity profile inside the channel is taken as the exit profile for the emerging free‐surface flow. The flow is assumed to be Newtonian as it emerges from the channel. An estimate of the magnitude of the rate‐of‐strain tensor components in the free‐surface region reveals that they are generally smaller than the shear rate inside the channel. The evolution of the flow front is simulated using the boundary element method. For the channel flow, the problem is reduced to a non‐linear dynamical system using the Galerkin projection method. Stability analysis indicates that the channel velocity may be linear or non‐linear depending on the range of the Weissenberg number. The evolution of the coating flow at the exit is examined for steady as well as transient (monotonic and oscillatory) channel flow. It is found that adverse flow can exist as a result of fluid elasticity, which can hinder the process of blade coating.

The purpose of this paper is to investigate the thermal convection inside a spatially modulated domain.

The governing equations are mapped onto an infinite strip, allowing Fourier expansion of the flow and temperature in the streamwise direction.

Similar to Rayleigh‐Benard convection, conduction is lost to convection at a critical Rayleigh number, which depends strongly on both the modulation amplitude and the wavenumber. The effect of modulation is found to be destabilizing (stabilizing) for conduction for relatively large (small) modulation wavelength. Oscillatory convection sets in as the Rayleigh number is increased.

This paper presents novel results.

This paper gives a review of the finite element techniques (FE) applied in the area of material processing. The latest trends in metal forming, non‐metal forming, powder metallurgy and composite material processing are briefly discussed. The range of applications of finite elements on these subjects is extremely wide and cannot be presented in a single paper; therefore the aim of the paper is to give FE researchers/users only an encyclopaedic view of the different possibilities that exist today in the various fields mentioned above. An appendix included at the end of the paper presents a bibliography on finite element applications in material processing for 1994‐1996, where 1,370 references are listed. This bibliography is an updating of the paper written by Brannberg and Mackerle which has been published in Engineering Computations, Vol. 11 No. 5, 1994, pp. 413‐55.

This numerical study explores the influence of the time dependence of material fluid parameters on the transient temperature evolution during the growth of fluid shells. The shell is spherical, the fluid is Newtonian, and the flow is induced by a constant driving pressure. The coupled heat and flow equations are solved numerically using the cobody (Lagrangian) transformation and a central difference discretization in space. The range of material values is adjusted from existing experiments. It is generally found that the variation in viscosity, surface tension and specific heat can have a significant influence on both the growth rate and temperature evolution. Thermal conductivity is found to be of little influence.

The purpose of this paper is to determine the thermo-gravitational convective state of a non-Fourier fluid layer of the single-phase-lagging type, heated from below. Unlike existing methodologies, the spectral modes are not imposed arbitrarily. They are systematically identified by expanding the spectral coefficients in terms of the relative departure in the post-critical Rayleigh number (perturbation parameter). The number and type of modes is determined to each order in the expansion. Non-Fourier effects become important whenever the relaxation time (delay in the response of the heat flux with respect to the temperature gradient) is of the same order of magnitude as process time.

In the spectral method the flow and temperature fields are expanded periodically along the layer and orthonormal shape functions are used in the transverse direction. A perturbation approach is developed to solve the nonlinear spectral system in the post-critical range.

The Nusselt number increases with non-Fourier effect as suggested in experiments in microscale and nanofluid convection.

Unlike existing nonlinear formulations for RB thermal convection, the present combined spectral-perturbation approach provides a systematic method for mode selection.

The purpose of this study is to perform a detailed review on the numerical modeling of multiphase and multicomponent flows in microfluidic system using phase-field method. The phase-field method is of emerging importance in numerical computation of transport phenomena involving multiple phases and/or components. This method is not only used to model interfacial phenomena typical to multiphase flows encountered in engineering and nature but also turns out to be a promising tool in modeling the dynamics of complex fluid-fluid interfaces encountered in physiological systems such as dynamics of vesicles and red blood cells). Intrinsically, a priori unknown topological evolution of interfaces offers to be the most concerning challenge toward accurate modeling of moving boundary problems. However, the numerical difficulties can be tackled simultaneously with numerical convenience and thermodynamic rigor in the paradigm of the phase field method.

The phase-field method replaces the macroscopically sharp interfaces separating the fluids by a diffuse transition layer where the interfacial forces are smoothly distributed. As against the moving mesh methods (Lagrangian) for the explicit tracking of interfaces, the phase-field method implicitly captures the same through the evolution of a phase-field function (Eulerian). In contrast to the deployment of an artificially smoothing function for the interface as used in the volume of a fluid or level set method, however, the phase-field method uses mixing free energy for describing the interface. This needs the consideration of an additional equation for an order parameter. The dynamic evolution of the system (equation for order parameter) can be described by Allen–Cahn or Cahn–Hilliard formulation, which couples with the Navier–Stokes equation with the aid of a forcing function that depends on the chemical potential and the gradient of the order parameter.

In this review, first, the authors discuss the broad motivation and the fundamental theoretical foundation associated with phase-field modeling from the perspective of computational microfluidics. They subsequently pinpoint the outstanding numerical challenges, including estimations of the model-free parameters. They outline some numerical examples, including electrohydrodynamic flows, to demonstrate the efficacy of the method. Finally, they pinpoint various emerging issues and futuristic perspectives connecting the phase-field method and computational microfluidics.

This paper gives unique perspectives to future directions of research on this topic.

Lattice Boltzmann method (LBM) has made great success in computational fluid dynamics, and this paper aims to establish an efficient simulation model for the polymer injection molding process using the LBM. The study aims to validate the capacity of the model for accurately predicting the injection molding process, to demonstrate the superior numerical efficiency in comparison with the current model based on the finite volume method (FVM).

The study adopts the stable multi-relaxation-time scheme of LBM to model the non-Newtonian polymer flow during the filling process. The volume of fluid method is naturally integrated to track the movement of the melt front. Additionally, a novel fractional-step thermal LBM is used to solve the convection-diffusion equation of the temperature field evolution, which is of high Peclet number. Through various simulation cases, the accuracy and stability of the present model are validated, and the higher numerical efficiency verified in comparison with the current FVM-based model.

The paper provides an efficient alternative to the current models in the simulation of polymer injection molding. Through the test cases, the model presented in this paper accurately predicts the filling process and successfully reproduces several characteristic phenomena of injection molding. Moreover, compared with the popular FVM-based models, the present model shows superior numerical efficiency, more fit for the future trend of parallel computing.

Limited by the authors’ hardware resources, the programs of the present model and the FVM-based model are run on parallel up to 12 threads, which is adequate for most simulations of polymer injection molding. Through the tests, the present model has demonstrated the better numerical efficiency, and it is recommended for the researcher to investigate the parallel performance on even larger-scale parallel computing, with more threads.

To the authors’ knowledge, it is for the first time that the lattice Boltzmann method is applied in the simulation of injection molding, and the proposed model does obviously better in numerical efficiency than the current popular FVM-based models.

The post-critical convective state for Rayleigh-Benard (RB) convection is studied using a nonlinear spectral-amplitude-perturbation approach in a fluid layer heated from below. The paper aims to discuss these issues.

In the spectral method the flow and temperature fields are expanded periodically along the layer and orthonormal shape functions are used in the transverse direction. A combined amplitude-perturbation approach is developed to solve the nonlinear spectral system in the post-critical range, even far from the linear stability threshold. Also, to leading order, the Lorenz model is recovered.

It is found that very small Prandtl numbers (Pr < 0.1) can change the Nusselt number, when terms to O(ε5/2) and higher are considered. However, to lower orders the Prandtl number does not affect the results. Variation of the Nusselt number to different orders is found to be highly consistent. Comparison with experimental results is made and a very good qualitative agreement is observed, even far from the linear threshold.

Unlike existing nonlinear formulations for RB thermal convection, the present combined spectral-perturbation approach provides a systematic method for mode selection. The number and type of modes to be included are directly related to the post-critical Rayleigh number. The method is not limited to the weakly nonlinear range.

Considers the filling stage of an injection moulding process for polymer melts through a transient finite element method. Presents results for a thin‐walled cup‐shaped mould problem. Adopts a generalized Newtonian formulation for inelastic non‐isothermal incompressible flow under cooling conditions. Gives some attention to the transient aspects of the problem, including the motion of the melt front and the numerical deformation patterns that develop. Advocates a localized remeshing technique that tracks the moving front effectively.