Magnetic convection heat transfer in an open ended enclosure filled with paramagnetic fluids

Numerical simulations of thermomagnetic convection of paramagnetic fluids placed in a micro-gravity condition (g nearly 0) and under a uniform vertical gradient magnetic field in an open ended square enclosure with ramp heating temperature condition applied on a vertical wall is investigated in this study. In presence of the strong magnetic gradient field thermal convection of the paramagnetic fluid might take place even in a zero-gravity environment as a direct consequence of temperature differences occurring within the fluid. The thermal boundary layer develops adjacent to the hot wall as soon as the ramp temperature condition is applied on it. There are two scenario that can be observed based on the ramp heating time. The steady state of the thermal boundary layer can be reached before the ramp time is finished or vice versa. If the ramp time is larger than the quasi-steady time then the thermal boundary layer is in a quasi-steady mode with convection balancing conduction after the quasi-steady time. Further increase of the heat input simply accelerates the flow to maintain the proper thermal balance. Finally, the boundary layer becomes completely steady state when the ramp time is finished. Effects of magnetic Rayleigh number, Prandtl number and paramagnetic fluid parameter on the flow pattern and heat transfer are presented.

Impact of the piston secondary motion on its slap force

A theoretical model is developed for the analysis of piston secondary motion. Based on this model, the slap force of a specific L6 diesel engine was compared when considering different boundary conditions, such as lubricating oil on cylinder liner, surface roughness, deformation of cylinder liner and piston skirt. It is concluded that it is necessary to consider the secondary motion of piston in the analysis of the inner excitation for an internal combustion engine. A more comprehensive consideration of the boundary condition (i.e., more close to the actual condition) will lead to a smaller maximum slap force, and among all boundary conditions considered in this paper, the structural deformation of the piston skirt and cylinder liner is the most influential factor. The theoretical model developed and findings obtained in this study will benefit the future analysis and design of advanced internal combustion engine structures.

Agent-based application on different boarding strategies

In this paper we implemented six different boarding strategies (Wilma, Steffen, Reverse Pyramid, Random, Blocks and By letter) in order to minimize boarding time and turnaround time for Boeing 777 and Airbus 380 aircrafts by using Agent-based modelling approach. In the simulation, we divided passengers into six different categories which are group size more than 5 people, passengers with child, gold members, first class passengers, business class passengers and economy class passengers. Results from the simulation demonstrates Reverse Pyramid method is the best boarding method for Boeing 777 and Steffen method is the best boarding method for Airbus 380.

Tuneable resonance properties of graphene by nitrogen-dopant

Doping as one of the popular methods to manipulate the properties of nanomaterials has received extensive application in deriving different types of graphene derivates, while the understanding of the resonance properties of dopant graphene is still lacking in literature. Based on the large-scale molecular dynamics simulation, reactive empirical bond order potential, as well as the tersoff potential, the resonance properties of N-doped graphene were studied. The studied samples were established according to previous experiments with the N atom’s percentage ranging from 0.43%-2.98%, including three types of N dopant locations, i.e., graphitic N, pyrrolic N and pyridinic N. It is found that different percentages of N-dopant exert different influence to the resonance properties of the graphene, while the amount of N-dopant is not the only factor that determines its impact. For all the considered cases, a relative large percentage of N-dopant (2.98% graphitic N-dopant) is observed to introduce significant influence to the profile of the external energy, and thus lead to an extremely low Q-factor comparing with that of the pristine graphene. The most striking finding is that, the natural frequency of the defective graphene with N-dopant appears uniformly larger than that of the pristine defective graphene. While for the perfect graphene, the N-dopant shows less influence to its natural frequency. This study will enrich the current understanding of the influence of dopants on graphene, which will eventually shed lights on the design of different molecules-doped graphene sheet.

Optimal Bayesian experimental design for models with intractable likelihoods using indirect inference applied to biological process models

This paper addresses the problem of determining optimal designs for biological process models with intractable likelihoods, with the goal of parameter inference. The Bayesian approach is to choose a design that maximises the mean of a utility, and the utility is a function of the posterior distribution. Therefore, its estimation requires likelihood evaluations. However, many problems in experimental design involve models with intractable likelihoods, that is, likelihoods that are neither analytic nor can be computed in a reasonable amount of time. We propose a novel solution using indirect inference (II), a well established method in the literature, and the Markov chain Monte Carlo (MCMC) algorithm of Müller et al. (2004). Indirect inference employs an auxiliary model with a tractable likelihood in conjunction with the generative model, the assumed true model of interest, which has an intractable likelihood. Our approach is to estimate a map between the parameters of the generative and auxiliary models, using simulations from the generative model. An II posterior distribution is formed to expedite utility estimation. We also present a modification to the utility that allows the Müller algorithm to sample from a substantially sharpened utility surface, with little computational effort. Unlike competing methods, the II approach can handle complex design problems for models with intractable likelihoods on a continuous design space, with possible extension to many observations. The methodology is demonstrated using two stochastic models; a simple tractable death process used to validate the approach, and a motivating stochastic model for the population evolution of macroparasites.

The renormalization group: A perturbation method for the graduate curriculum

In this paper the renormalization group (RG) method of Chen, Goldenfeld, and Oono [Phys. Rev. Lett., 73 (1994), pp.1311-1315; Phys. Rev. E, 54 (1996), pp.376-394] is presented in a pedagogical way to increase its visibility in applied mathematics and to argue favorably for its incorporation into the corresponding graduate curriculum.The method is illustrated by some linear and nonlinear singular perturbation problems. Key word. © 2012 Society for Industrial and Applied Mathematics.

A survey in mathematics for industry: Two-timing and matched asymptotic expansions for singular perturbation problems

Following the derivation of amplitude equations through a new two-time-scale method [O'Malley, R. E., Jr. & Kirkinis, E (2010) A combined renormalization group-multiple scale method for singularly perturbed problems. Stud. Appl. Math. 124, 383-410], we show that a multi-scale method may often be preferable for solving singularly perturbed problems than the method of matched asymptotic expansions. We illustrate this approach with 10 singularly perturbed ordinary and partial differential equations. © 2011 Cambridge University Press.

On helical shear of a compressible elastic circular cylindrical tube

In this paper we examine the combined azimuthal and axial shear of a compressible isotropic elastic circular cylindrical tube of finite extent, otherwise referred to as helical shear (which is an isochoric deformation). The equilibrium equations are formulated in terms of the principal stretches, and explicit necessary and sufficient conditions on the strain-energy function for the material to support this deformation are obtained and compared with those obtained previously for this problem. Several classes of strain-energy functions are derived and in some general cases complete solutions of the equilibrium equations are obtained. Existing results are recovered as special cases and some new results for the strain-energy functions derived are determined and discussed.

On extension and torsion of a compressible elastic circular cylinder

In this paper we examine the combined extension and torsion of a compressible isotropic elastic cylinder of finite extent. The equilibrium equations are formulated in terms of the principal stretches and then applied to the special case of pure torsion superimposed on a uniform extension (an isochoric deformation). Explicit necessary and sufficient conditions on the strain-energy function for the material to support this deformation with vanishing traction on the lateral surfaces of the cylinder are obtained. Some strain-energy functions satisfying these conditions are considered, existing results are recovered as special cases and new results are obtained. We also point out how the strain-energy functions generated from the considered isochoric deformation considered (of a compressible material) can be used to generate energy functions and corresponding solutions for the incompressible theory.

Surface science of plasma exposed surfaces : a challenge for applied plasma science

This article introduces a deterministic approach to using low-temperature, thermally non-equilibrium plasmas to synthesize delicate low-dimensional nanostructures of a small number of atoms on plasma exposed surfaces. This approach is based on a set of plasma-related strategies to control elementary surface processes, an area traditionally covered by surface science. Major issues related to balanced delivery and consumption of building units, appropriate choice of process conditions, and account of plasma-related electric fields, electric charges and polarization effects are identified and discussed in the quantum dot nanoarray context. Examples of a suitable plasma-aided nanofabrication facility and specific effects of a plasma-based environment on self-organized growth of size- and position-uniform nanodot arrays are shown. These results suggest a very positive outlook for using low-temperature plasma-based nanotools in high-precision nanofabrication of self-assembled nanostructures and elements of nanodevices, one of the areas of continuously rising demand from academia and industry.

Corner and finger formation in Hele-Shaw flow with kinetic undercooling regularisation

We examine the effect of a kinetic undercooling condition on the evolution of a free boundary in Hele--Shaw flow, in both bubble and channel geometries. We present analytical and numerical evidence that the bubble boundary is unstable and may develop one or more corners in finite time, for both expansion and contraction cases. This loss of regularity is interesting because it occurs regardless of whether the less viscous fluid is displacing the more viscous fluid, or vice versa. We show that small contracting bubbles are described to leading order by a well-studied geometric flow rule. Exact solutions to this asymptotic problem continue past the corner formation until the bubble contracts to a point as a slit in the limit. Lastly, we consider the evolving boundary with kinetic undercooling in a Saffman--Taylor channel geometry. The boundary may either form corners in finite time, or evolve to a single long finger travelling at constant speed, depending on the strength of kinetic undercooling. We demonstrate these two different behaviours numerically. For the travelling finger, we present results of a numerical solution method similar to that used to demonstrate the selection of discrete fingers by surface tension. With kinetic undercooling, a continuum of corner-free travelling fingers exists for any finger width above a critical value, which goes to zero as the kinetic undercooling vanishes. We have not been able to compute the discrete family of analytic solutions, predicted by previous asymptotic analysis, because the numerical scheme cannot distinguish between solutions characterised by analytic fingers and those which are corner-free but non-analytic.

Stress relaxation analysis of single chondrocytes using porohyperelastic model based on AFM experiments

The aim of this study is to investigate the stress relaxation behavior of single chondrocytes using the Porohyperelastic (PHE) model and inverse Finite Element Analysis (FEA). Firstly, based on Atomic Force Microscopy (AFM) technique, we have found that the chondrocytes exhibited stress relaxation behavior. We explored the mechanism of this stress relaxation behavior and concluded that the intracellular fluid exuding out from the cells during deformation plays the most important role in the stress relaxation. Next, we have applied the inverse FEA technique to determine necessary material parameters for PHE model to simulate this stress relaxation behavior as this model is proven capable of capturing the non-linear behavior and the fluid-solid interaction during the stress relaxation of the single chondrocytes. It is observed that this PHE model can precisely capture the stress relaxation behavior of single chondrocytes and would be a suitable model for cell biomechanics.

Determination of strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes using atomic force microscopy and inverse finite element analysis

The aim of this paper is to determine the strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes, in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of these single cells at four different strain-rates. These results were then employed in inverse finite element analysis (FEA) using Modified Standard neo-Hookean Solid (MSnHS) idealization of these cells to determine their mechanical properties. In addition, a FEA model with a newly developed spring element was employed to accurately simulate AFM evaluation in this study. We report that both cytoskeleton (CSK) and intracellular fluid govern the strain-rate-dependent mechanical property of living cells whereas intracellular fluid plays a predominant role on fixed cells’ behavior. In addition, through the comparisons, it can be concluded that osteocytes are stiffer than chondrocytes at all strain-rates tested indicating that the cells could be the biomarker of their tissue origin. Finally, we report that MSnHS is able to capture the strain-rate-dependent mechanical behavior of osteocyte and chondrocyte for both living and fixed cells. Therefore, we concluded that the MSnHS is a good model for exploration of mechanical deformation responses of single osteocytes and chondrocytes. This study could open a new avenue for analysis of mechanical behavior of osteocytes and chondrocytes as well as other similar types of cells.