Integral approach and numerical improvement to calculate carbon concentration profiles in carburising

The carbon diffusion in steel, where the carbon diffusivity varies with the carbon content, was solved with the integral methods under the third boundary condition. The variation of carbon diffusivity in steel with the carbon content was described with two different functions, linear dependence and exponential dependence. The integral approximation for both cases was improved with the numerical computation to more accurately predict the carbon profiles. The integral solution is more accurate than the formulation based on the assumption of a constant diffusivity or those based on the assumption of a constant diffusivity and/or constant carbon content at part surface. It is also more easily used in practice than the numerical method to describe the carburising process and predict the carbon content at steel surface and carbon profiles in treated layer.

Optimization of bypass graft using FSI numerical method

It is well documented in literature that the coronary artery bypass graft is normally fail after a short period of time, due to the development of plaque known as intimal hyperplasia within the graft. Various in vivo and in vitro studies have linked the development of intimal hyperplasia to the abnormal hemodynamics and compliance mismatch. Therefore, it is essential to fully understand the relationship between the hemodynamics inside the coronary artery bypass and its mechanical and geometrical characteristics under the correct physiological conditions. In this work, hemodynamic of the bypass graft is studied numerically. The effect of the host and graft diameters ratio, the angle of anastomosis and the graft configuration on the local flow patterns and the distribution of wall shear stress are examined. The pulsatile waveforms boundary conditions are adopted from in vivo measurement data to study the hemodynamics of composite grafts namely Consequence and Y grafting in terms temporal and spatial distributions of the blood flows. Moreover, various non-Newtonian and Newtonian models of blood have been carried out to examine the numerical simulation of blood flow in stenosis artery. The results are presented and discussed for various operating conditions.

Simulation and analysis of a sub-wavelength grating based multilayer surface plasmon resonance biosensor

This paper presents a subwavelength grating based multilayer surface plasmon resonance biosensor (SPRB) which includes a periodic array of subwavelength grating on top of a layer of graphene sheet in the biosensor. The proposed biosensor is named grating-graphene SPRB (GG-SPRB). The aim of the proposed multilayer structure is to improve the sensitivity of the SPRB through monitoring of the biomolecular interactions of DNA hybridization. Significant sensitivity improvement is obtained for the GG-SPRB compared with the conventional SPRB. The result of the numerical investigation of the GG-SPRB is presented and compared with a theoretically developed multilayer matrix formalism, and a good agreement has been observed. In addition, an optimization of the grating dimensions including volume factor, grating depth, grating angle, grating period, and grating geometry (e.g., rectangular, sinusoidal and triangular) is presented. The outcome of the investigation presented in this paper identifies desired functioning conditions corresponding to the best design parameters for the GG-SPRB.

A mathematical model of neuromuscular adaptation to resistance training and its application in a computer simulation of accommodating loads

A large corpus of data obtained by means of empirical study of neuromuscular adaptation is currently of limited use to athletes and their coaches. One of the reasons lies in the unclear direct practical utility of many individual trials. This paper introduces a mathematical model of adaptation to resistance training, which derives its elements from physiological fundamentals on the one side, and empirical findings on the other. The key element of the proposed model is what is here termed the athlete’s capability profile. This is a generalization of length and velocity dependent force production characteristics of individual muscles, to an exercise with arbitrary biomechanics. The capability profile, a two-dimensional function over the capability plane, plays the central role in the proposed model of the training-adaptation feedback loop. Together with a dynamic model of resistance the capability profile is used in the model’s predictive stage when exercise performance is simulated using a numerical approximation of differential equations of motion. Simulation results are used to infer the adaptational stimulus, which manifests itself through a fed back modification of the capability profile. It is shown how empirical evidence of exercise specificity can be formulated mathematically and integrated in this framework. A detailed description of the proposed model is followed by examples of its application—new insights into the effects of accommodating loading for powerlifting are demonstrated. This is followed by a discussion of the limitations of the proposed model and an overview of avenues for future work.

Fracture modelling of DP780 sheets using a hybrid experimental-numerical method and two-dimensional digital image correlation

This paper is concerned with the construction of fracture envelopes of DP780 sheets using two methods: a hybrid experimental-numerical method; two-dimensional digital image correlation (2D-DIC). For the hybrid method, four types of ductile fracture tests were carried out covering a wide range of stress states on specimens: with a central hole; two symmetric circular notches; flat grooved; and diagonally double-notched. Based on the fracture strain and loading paths identified with finite element simulation, a fracture envelope was obtained by employing the three-parameter modified Mohr-Coulomb fracture model. In addition, the fracture surface strain was directly measured using 2D-DIC. Loading histories of each test were extracted from a surface element of a three dimensional finite element model. The comparison of fracture envelopes constructed by the two methods reveals that there is little difference. Thus, it can be concluded that 2D-DIC is applicable to fracture modelling of DP780 sheets despite the assumption of the plane stress condition even after necking

A non-associated plasticity model with anisotropic and nonlinear kinematic hardening for simulation of sheet metal forming

A material model for more thorough analysis of plastic deformation of sheet materials is presented in this paper. This model considers the following aspects of plastic deformation behavior of sheet materials: (1) the anisotropy in yield stresses and in work hardening by using Hill's 1948 quadratic yield function and non-constant stress ratios which leads to different flow stress hardening in different directions, (2) the anisotropy in plastic strains by using a quadratic plastic potential function and non-associated flow rule, also based on Hill's 1948 model and r-values, and (3) the cyclic hardening phenomena such as the Bauschinger effect, permanent softening and transient behavior for reverse loading by using a coupled nonlinear kinematic hardening model. Plasticity fundamentals of the model were derived in a general framework and the model calibration procedure was presented for the plasticity formulations. Also, a generic numerical stress integration procedure was developed based on backward-Euler method, so-called multi-stage return mapping algorithm. The model was implemented in the framework of the finite element method to evaluate the simulation results of sheet metal forming processes. Different aspects of the model were verified for two sheet metals, namely DP600 steel and AA6022 aluminum alloy. Results show that the new model is able to accurately predict the sheet material behavior for both anisotropic hardening and cyclic hardening conditions. The drawing of channel sections and the subsequent springback were also simulated with this model for different drawbead configurations. Simulation results show that the current non-associated anisotropic hardening model is able to accurately predict the sidewall curl in the drawn channel sections.

Numerical modeling and analysis for forming process of dual-phase 980 steel exposed to infrared local heating

A recent experiment confirmed that the infrared (IR) local heating method drastically reduces springback of dual-phase (DP) 980 sheets. In the experiment, only the plastic deformation zone of the sheets was locally heated using condensed IR heating. The heated sheets were then deformed by V-bending or 2D-draw bending. Although the experimental observation proved the merit of using the IR local heating to reduce springback, numerical modeling has not been reported. Numerical modeling has been required to predict springback and improve the understanding of the forming process. This paper presents a numerical modeling for V-bending and 2D-draw bending of DP 980 sheets exposed to the IR local heating with the finite element method (FEM). For describing the thermo-mechanical behavior of the DP 980 sheet, a flow stress model which includes a function of temperature and effective plastic strain was newly implemented into Euler-backward stress integration method. The numerical analysis shows that the IR local heating reduces the level of stress in the deformation zone, although it heats only the limited areas, and then it reduces the springback. The simulation also provides a support that the local heating method has an advantage of shape accuracy over the method to heat the material as a whole in V-bending. The simulated results of the springback in both V-bending and 2D-draw bending also show good predictions.

Consensus building in a local context for the AHP-GDM with the individual numerical scale and prioritization method

The methodology for selecting the individual numerical scale and prioritization method has recently been presented and justified in the analytic hierarchy process (AHP). In this study, we further propose a novel AHP-group decision making (GDM) model in a local context (a unique criterion), based on the individual selection of the numerical scale and prioritization method. The resolution framework of the AHP-GDM with the individual numerical scale and prioritization method is first proposed. Then, based on linguistic Euclidean distance (LED) and linguistic minimum violations (LMV), the novel consensus measure is defined so that the consensus degree among decision makers who use different numerical scales and prioritization methods can be analyzed. Next, a consensus reaching model is proposed to help decision makers improve the consensus degree. In this consensus reaching model, the LED-based and LMV-based consensus rules are proposed and used. Finally, a new individual consistency index and its properties are proposed for the use of the individual numerical scale and prioritization method in the AHP-GDM. Simulation experiments and numerical examples are presented to demonstrate the validity of the proposed model.

Towards simulation of NASA35 axial compressor

Purpose – This paper aims to validate and analyse the NASA35 axial compressor performance based on a numerical approach. Design/methodology/approach – Knowledge about flow property change during compressor operation at high and relatively low speed is still limited. This work provides a numerical approach to address these problems. Validation of numerical methods is proposed to generate confidence the numerical approach adopted, and after that, analysis of compressor performance at different operation conditions is carried out. Findings – The numerical methods proposed are proved capable in predicting compressor performance. Changes of flow property during compressor operation are discussed and explained. Research limitations/implications – The current numerical work is carried out based on the first stage of the NASA35 axial compressor, where the interactive effects from adjacent stage are not counted in. Furthermore, the steady-state simulation enforces an averaging of flow at rotor-stator interface, where the transient rotor-stator interaction is removed. Practical implications – This work validates the numerical methods used in the prediction of NASA35 axial compressor performance, and a similar numerical approach can be used for other turbomachinery simulation cases. Originality/value – This work reinforces the understanding of axial compressor operation and provides reliable results for further investigation of a similar type of compressor. In addition, details of flow field within the NASA35 compressor during operation are given and explained which experiments still have difficult to achieve.

Measuring eco-roof mitigation on flash floods via GIS simulation

Purpose: The rapid and ongoing expansion of urbanised impervious areas could lead to more frequent flood inundation in urban flood-prone regions. Nowadays, urban flood inundation induced by rainstorm is an expensive natural disaster in many countries. In order to reduce the flooding risk, eco-roof systems (or green roof systems) could be considered as an effective mechanism of mitigating flooding disasters through their rainwater retention capability. However, there is still a lack of examining the stormwater management tool. The purpose of this paper is to evaluate the effects on flooding disaster from extensive green roofs. Design/methodology/approach: Based on geographical information system (GIS) simulation, this research presents a frame of assessing eco-roof impacts on urban flash floods. The approach addresses both urban rainfall-runoff and underground hydrologic models for traditional impervious and green roofs. Deakin University’s Geelong Waurn Ponds campus is chosen as a study case. GIS technologies are then utilised to visualise and analyse the effects on flood inundation from surface properties of building roofs. Findings: The results reveal that the eco-roof systems generate varying degrees of mitigation of urban flood inundation with different return period storms. Originality/value: Although the eco-roof technology is considered as an effective stormwater management tool, it is not commonly adopted and examined in urban floods. This study will bring benefits to urban planners for raising awareness of hazard impacts and to construction technicians for considering disaster mitigation via roof technologies. The approach proposed here could be used for the disaster mitigation in future urban planning.

Numerical investigation of the turbulent energy budget in the wake of freely oscillating elastically mounted cylinder at low reduced velocities

We present a numerical study of the turbulent kinetic energy budget in the wake of cylinders undergoing Vortex-Induced Vibration (VIV). We show three-dimensional Large Eddy Simulations (LES) of an elastically mounted circular cylinder in the synchronization regime at Reynolds number of Re=8000. The Immersed Boundary Method (IBM) is used to account for the presence of the cylinder. The flow field in the wake is decomposed using the triple decomposition splitting the flow variables in mean, coherent and stochastic components. The energy transfer between these scales of motions are then studied and the results of the free oscillation are compared to those of a forced oscillation. The turbulent kinetic energy budget shows that the maximum amplitude of VIV is defined by the ability of the mean flow to feed energy to the coherent structures in the wake. At amplitudes above this maximum amplitude, the energy of the coherent structures needs to be fed additionally by small scale, stochastic energy in form of backscatter to sustain its motion. Furthermore, we demonstrate that the maximum amplitude of the VIV is defined by the integral length scale of the turbulence in the wake