Numerical study of two ill-posed one-phase Stefan problems

We treat two related moving boundary problems. The first is the ill-posed Stefan problem for melting a superheated solid in one Cartesian coordinate. Mathematically, this is the same problem as that for freezing a supercooled liquid, with applications to crystal growth. By applying a front-fixing technique with finite differences, we reproduce existing numerical results in the literature, concentrating on solutions that break down in finite time. This sort of finite-time blow-up is characterised by the speed of the moving boundary becoming unbounded in the blow-up limit. The second problem, which is an extension of the first, is proposed to simulate aspects of a particular two-phase Stefan problem with surface tension. We study this novel moving boundary problem numerically, and provide results that support the hypothesis that it exhibits a similar type of finite-time blow-up as the more complicated two-phase problem. The results are unusual in the sense that it appears the addition of surface tension transforms a well-posed problem into an ill-posed one.

Advanced numerical study of impact forces on railway sleepers under wheel flat

Wheel-rail interaction is one of the most important research topics in railway engineering. It includes track vibration, track impact response and safety of the track. Track structure failures caused by impact forces can lead to significant economic loss for track owners through damage to rails and to the sleepers beneath. The wheel-rail impact forces occur because of imperfections on the wheels or rails such as wheel flats, irregular wheel profile, rail corrugation and differences in the height of rails connected at a welded joint. In this paper, a finite element model for the wheel flat study is developed by use of the FEA software package ANSYS. The effect of the wheel flat to impact force on sleepers is investigated. It has found that the wheel flat significantly increases impact forces and maximum Von Mises stress, and also delays the peak position of dynamic variation for impact forces on both rail and sleeper.

Numerical study of cold-formed steel beam subject to lateral-torsional buckling

The use of cold-formed steel members as structural columns and beams in residential, industrial and commercial buildings has increased significantly in recent times. This study is focused on the use of cold-formed steel sections as flexural members subject to lateral-torsional buckling. For this purpose a finite element model of a simply supported lipped channel beam under uniform bending was developed, validated using available numerical and experimental results, and used in a detailed parametric study. The moment capacity results were then compared with the predictions from the current ambient temperature design rules in the cold-formed steel structures codes of Australia, New Zealand, North America and Europe. European design rules were found to be conservative while Australian and American design rules were unsafe. This paper presents the results of the numerical study, the comparison with the current design rules and the new proposed design rules.

Coupling of energy from quantum emitters to the plasmonic mode of V groove waveguides: A numerical study

This work is a theoretical investigation into the coupling of a single excited quantum emitter to the plasmon mode of a V groove waveguide. The V groove waveguide consists of a triangular channel milled in gold and the emitter is modeled as a dipole emitter, and could represent a quantum dot, nitrogen vacancy in diamond, or similar. In this work the dependence of coupling efficiency of emitter to plasmon mode is determined for various geometrical parameters of the emitter-waveguide system. Using the finite element method, the effect on coupling efficiency of the emitter position and orientation, groove angle, groove depth, and tip radius, is studied in detail. We demonstrate that all parameters, with the exception of groove depth, have a significant impact on the attainable coupling efficiency. Understanding the effect of various geometrical parameters on the coupling between emitters and the plasmonic mode of the waveguide is essential for the design and optimization of quantum dot–V groove devices.

Numerical study of impact forces on railway sleepers under wheel flat

Wheel-rail interaction is one of the most important research topics in railway engineering. It includes track vibration, track impact response and safety of the track. Track structure failures caused by impact forces can lead to significant economic loss for track owners through damage to rails and to the sleepers beneath. The wheel-rail impact forces occur because of imperfections on the wheels or rails such as wheel flats, irregular wheel profile, rail corrugation and differences in the height of rails connected at a welded joint. The vehicle speed and static wheel load are important factors of the track design, because they are related to the impact forces under wheel-rail defects. In this paper, a 3-Dimensional finite element model for the study of wheel flat impact is developed by use of the FEA software package ANSYS. The effects of the wheel flat to impact force on sleepers with various speeds and static wheel loads under a critical wheel flat size are investigated. It has found that both wheel-rail impact force and impact force on sleeper induced by wheel flat are varying nonlinearly by increasing the vehicle speed; both impact forces are nonlinearly and monotonically increasing by increasing the static wheel load. The relationships between both of impact forces induced by wheel flat and vehicles speed or static load are important to the track engineers to improve the design and maintenance methods in railway industry.

Simulation of phosphine flow in vertical grain storage : a preliminary numerical study

To fumigate grain stored in a silo, phosphine gas is distributed by a combination of diffusion and fan-forced advection. This initial study of the problem mainly focuses on the advection, numerically modelled as fluid flow in a porous medium. We find satisfactory agreement between the flow predictions of two Computational Fluid Dynamics packages, Comsol and Fluent. The flow predictions demonstrate that the highest velocity (>0.1 m/s) occurs less than 0.2m from the inlet and reduces drastically over one metre of silo height, with the flow elsewhere less than 0.002 m/s or 1% of the velocity injection. The flow predictions are examined to identify silo regions where phosphine dosage levels are likely to be too low for effective grain fumigation.

Experimental and numerical study of water-filled road safety barriers

Portable water-filled road barriers (PWFB) are roadside structures placed on temporary construction zones to separate work site from moving traffic. Recent changes in governing standards require PWFB to adhere to strict compliance in terms of lateral displacement of the road barriers and vehicle redirectionality. Actual road safety barrier test can be very costly, thus researchers resort to Finite Element Analysis (FEA) in the initial designs phase prior to real vehicle test. There has been many research conducted on concrete barriers and flexible steel barriers using FEA, however not many is done pertaining to PWFB. This research probes a new method to model joint mechanism in PWFB. Two methods to model the joining mechanism are presented and discussed in relation to its practicality and accuracy to real work applications. Moreover, the study of the physical gap and mass of the barrier was investigated. Outcome from this research will benefit PWFB research and allow road barrier designers better knowledge in developing the next generation of road safety structures.

Numerical study of lid-driven square cavity with heat generation using LBM

In this study, the mixed convection heat transfer and fluid flow behaviors in a lid-driven square cavity filled with high Prandtl number fluid (Pr = 5400, ν = 1.2×10-4 m2/s) at low Reynolds number is studied using thermal Lattice Boltzmann method (TLBM) where ν is the viscosity of the fluid. The LBM has built up on the D2Q9 model and the single relaxation time method called the Lattice-BGK (Bhatnagar-Gross-Krook) model. The effects of the variations of non dimensional mixed convection parameter called Richardson number(Ri) with and without heat generating source on the thermal and flow behavior of the fluid inside the cavity are investigated. The results are presented as velocity and temperature profiles as well as stream function and temperature contours for Ri ranging from 0.1 to 5.0 with other controlling parameters that present in this study. It is found that LBM has good potential to simulate mixed convection heat transfer and fluid flow problem. Finally the simulation results have been compared with the previous numerical and experimental results and it is found to be in good agreement.

Experimental and numerical study of polymeric foam efficacy in portable water filled barriers

Portable, water filled road safety barriers are used to provide protection and reduce the potential hazard due to errant vehicles in areas where the road conditions change frequently (e.g. near road work sites). As part of an effort to reduce excessive working widths typical of these systems, a study was conducted to assess the effectiveness of introducing polymeric foam filled panels into the design. Surrogate impact tests of a design typical of such as barrier system were conducted utilising a pneumatically powered horizontal impact testing machine up to impact energies of 7.40 kJ. Results of these tests are utilised to examine the barrier behaviour, in addition to being used to validate a couple FE/SPH model of the barrier system. Once validated, the FE/SPH model it utilised as the basis for a parametric study into the efficacy and effects of the inclusion of polymeric foam filled panels on the performance of portable water filled road safety barriers. It was found that extruded polystyrene foam functioned well, with a greater thickness of the foam panel significantly reducing the impacting body velocity as the barrier began to translate.

Numerical study of turbulent fluid flow and heat transfer in lateral perforated extended surfaces

Numerical study has been performed in this study to investigate the turbulent convection heat transfer on a rectangular plate mounted over a flat surface. Thermal and fluid dynamic performances of extended surfaces having various types of lateral perforations with square, circular, triangular and hexagonal cross sections are investigated. RANS (Reynolds averaged Navier–Stokes) based modified k–ω turbulence model is used to calculate the fluid flow and heat transfer parameters. Numerical results are compared with the results of previously published experimental data and obtained results are in reasonable agreement. Flow and heat transfer parameters are presented for Reynolds numbers from 2000 to 5000 based on the fin thickness.

Numerical study of LSF floors made of hollow flange channels in fire

Light gauge steel frame (LSF) floor systems are generally made of lipped channel section joists and lined with gypsum plasterboards to provide adequate fire resistance rating under fire conditions. Recently a new LSF floor system made of welded hollow flange channel (HFC) section was developed and its fire performance was investigated using full scale fire tests. The new floor systems gave higher fire resistance ratings in comparison to conventional LSF floor systems. To avoid expensive and time consuming full scale fire tests, finite element analyses were also performed to simulate the fire performance of LSF floors made of HFC joists using both steady and transient state methods. This paper presents the details of the developed finite element models of HFC joists to simulate the structural fire performance of the LSF floor systems under standard fire conditions. Finite element analyses were performed using the measured time–temperature profiles of the failed joists from the fire tests, and their failure times, temperatures and modes, and deflection versus time curves were obtained. The developed finite element models successfully predicted the structural performance of LSF floors made of HFC joists under fire conditions. They were able to simulate the complex behaviour of thin cold-formed steel joists subjected to non-uniform temperature distributions, and local buckling and yielding effects. This study also confirmed the superior fire performance of the newly developed LSF floors made of HFC joists.

A 3D numerical study of the collateral capacity of the circle of Willis with anatomical variation in the posterior circulation

Background The Circle of Willis (CoW) is the most important collateral pathway of the cerebral artery. The present study aims to investigate the collateral capacity of CoW with anatomical variation when unilateral internalcarotid artery (ICA) is occluded. Methods Basing on MRI data, we have reconstructed eight 3D models with variations in the posterior circulation of the CoW and set four different degrees of stenosis in the right ICA, namely 24%, 43%, 64% and 79%, respectively. Finally, a total of 40 models are performed with computational fluid dynamics simulations. All of the simulations share the same boundary condition with static pressure and the volume flow rate (VFR) are obtained to evaluate their collateral capacity. Results As for the middle cerebral artery (MCA) and the anterior cerebral artery (ACA), the transitional-type model possesses the best collateral capacity. But for the posterior cerebral artery (PCA), unilateral stenosis of ICA has the weakest influence on the unilateral posterior communicating artery (PCoA) absent model. We also find that the full fetal-type posterior circle of Willis is an utmost dangerous variation which must be paid more attention. Conclusion The results demonstrate that different models have different collateral capacities in coping stenosis of unilateral ICA and these differences can be reflected by different outlets. The study could be used as a reference for neurosurgeon in choosing the best treatment strategy.

A numerical study of thermocapillary flow in a rectangular cavity during laser melting

The surface tension gradient driven flow that occurs during laser melting has been studied. The vorticity-streamfunction form of the Navier-Stokes equations and the energy equation has been solved by the ‘Alternative Direction Implicit’ method. It has been shown that the inertia forces in the melt strongly influence the flow pattern in the melt. The convection in the melt modifies the isotherms in the melt at high surface tension Reynolds number and high Prandtl number. The buoyancy driven flow has been shown to be negligible compared to the surface tension gradient driven flow in laser melting.