Plasmonic effect of annealed gold islands for improving efficiency of organic solar cells

Embedding metallic nanoparticles in organic solar cells can enhance the photoabsorption through light trapping processes. This paper investigates how gold islands obtained by annealing 1–5 nm thick Au layers affect the photoabsorption. Using finite-difference time-domain simulations, the cell efficiency for various island geometries and thicknesses are analyzed and the properties of the islands for maximal photocurrent are discussed. It is shown that a careful choice of size and concentration of gold islands could contribute to enhance the power conversion efficiencies when compared to standard organic solar cell devices. The conclusions are then compared to experimental data for thermally annealed gold islands in bulk heterojunction solar cells. The results of this paper will contribute to the optimization of plasmonic organic solar cell systems and will pave the way for the development of highly efficient organic solar cell devices.

Numerical investigation of case hardening of plant tissue during dying and its influence on the cellular level shrinkage

Dried plant food materials are one of the major contributors to the global food industry. Widening the fundamental understanding on different mechanisms of food material alterations during drying assists the development of novel dried food products and processing techniques. In this regard, case hardening is an important phenomenon, commonly observed during the drying processes of plant food materials, which significantly influences the product quality and process performance. In this work, a recent meshfree-based numerical model of the authors is further improved and used to simulate the influence of case hardening on shrinkage characteristics of plant tissues during drying. In order to model fluid and wall mechanisms in each cell, Smoothed Particle Hydrodynamics (SPH) and the Discrete Element Method (DEM) are used. The model is fundamentally more capable of simulating large deformation of multiphase materials, when compared with conventional grid-based modelling techniques such as Finite Element Methods (FEM) or Finite Difference Methods (FDM). Case hardening is implemented by maintaining distinct moisture levels in the different cell layers of a given tissue. In order to compare and investigate different factors influencing tissue deformations under case hardening, four different plant tissue varieties (apple, potato, carrot and grape) are studied. The simulation results indicate that the inner cells of any given tissue undergo limited shrinkage and cell wall wrinkling compared to the case hardened outer cell layers of the tissues. When comparing unique deformation characteristics of the different tissues, irrespective of the normalised moisture content, the cell size, cell fluid turgor pressure and cell wall characteristics influence the tissue response to case hardening.

The effect of thermal radiation on the natural convection boundary layer flow over a wavy horizontal surface

In this article, natural convection boundary layer flow is investigated over a semi-infinite horizontal wavy surface. Such an irregular (wavy) surface is used to exchange heat with an external radiating fluid which obeys Rosseland diffusion approximation. The boundary layer equations are cast into dimensionless form by introducing appropriate scaling. Primitive variable formulations (PVF) and stream function formulations (SFF) are independently used to transform the boundary layer equations into convenient form. The equations obtained from the former formulations are integrated numerically via implicit finite difference iterative scheme whereas equations obtained from lateral formulations are simulated through Keller-box scheme. To validate the results, solutions produced by above two methods are compared graphically. The main parameters: thermal radiation parameter and amplitude of the wavy surface are discussed categorically in terms of shear stress and rate of heat transfer. It is found that wavy surface increases heat transfer rate compared to the smooth wall. Thus optimum heat transfer is accomplished when irregular surface is considered. It is also established that high amplitude of the wavy surface in the boundary layer leads to separation of fluid from the plate.

Non-Newtonian mixed convection flow from a horizontal circular cylinder with uniform surface heat flux

Mixed convection laminar two-dimensional boundary-layer flow of non-Newtonian pseudo-plastic fluids is investigated from a horizontal circular cylinder with uniform surface heat flux using a modified power-law viscosity model, that contains no unrealistic limits of zero or infinite viscosity; consequently, no irremovable singularities are introduced into boundary-layer formulations for such fluids. The governing boundary layer equations are transformed into a non-dimensional form and the resulting nonlinear systems of partial differential equations are solved numerically applying marching order implicit finite difference method with double sweep technique. Numerical results are presented for the case of shear-thinning fluids in terms of the fluid temperature distributions, rate of heat transfer in terms of the local Nusselt number.

A computer code for enclosed natural convection in the boundary layer regime

A computer code is developed for the numerical prediction of natural convection in rectangular two-dimensional cavities at high Rayleigh numbers. The governing equations are retained in the primitive variable form. The numerical method is based on finite differences and an ADI scheme. Convective terms may be approximated with either central or hybrid differencing for greater stability. A non-uniform grid distribution is possible for greater efficiency. The pressure is dealt with via a SIMPLE type algorithm and the use of a fast elliptic solver for the solenoidal velocity correction field significantly reduces computing times. Preliminary results indicate that the code is reasonably accurate, robust and fast compared with existing benchmarks and finite difference based codes, particularly at high Rayleigh numbers. Extension to three-dimensional problems and turbulence studies in similar geometries is readily possible and indicated.

Implicit difference approximation for the time fractional diffusion equation

In this paper, we consider a time fractional diffusion equation on a finite domain. The equation is obtained from the standard diffusion equation by replacing the first-order time derivative by a fractional derivative (of order $0<\alpha<1$ ). We propose a computationally effective implicit difference approximation to solve the time fractional diffusion equation. Stability and convergence of the method are discussed. We prove that the implicit difference approximation (IDA) is unconditionally stable, and the IDA is convergent with $O(\tau+h^2)$, where $\tau$ and $h$ are time and space steps, respectively. Some numerical examples are presented to show the application of the present technique.

Creation of a validated 3D finite element model of an ovine tibia

Introduction Ovine models are widely used in orthopaedic research. To better understand the impact of orthopaedic procedures computer simulations are necessary. 3D finite element (FE) models of bones allow implant designs to be investigated mechanically, thereby reducing mechanical testing. Hypothesis We present the development and validation of an ovine tibia FE model for use in the analysis of tibia fracture fixation plates. Material & Methods Mechanical testing of the tibia consisted of an offset 3-pt bend test with three repetitions of loading to 350N and return to 50N. Tri-axial stacked strain gauges were applied to the anterior and posterior surfaces of the bone and two rigid bodies – consisting of eight infrared active markers, were attached to the ends of the tibia. Positional measurements were taken with a FARO arm 3D digitiser. The FE model was constructed with both geometry and material properties derived from CT images of the bone. The elasticity-density relationship used for material property determination was validated separately using mechanical testing. This model was then transformed to the same coordinate system as the in vitro mechanical test and loads applied. Results Comparison between the mechanical testing and the FE model showed good correlation in surface strains (difference: anterior 2.3%, posterior 3.2%). Discussion & Conclusion This method of model creation provides a simple method for generating subject specific FE models from CT scans. The use of the CT data set for both the geometry and the material properties ensures a more accurate representation of the specific bone. This is reflected in the similarity of the surface strain results.

Finite element shoulder models

Shoulder joint is a complex integration of soft and hard tissues. It plays an important role in performing daily activities and can be considered as a perfect compromise between mobility and stability. However, shoulder is vulnerable to complications such as dislocations and osteoarthritis. Finite element (FE) models have been developed to understand shoulder injury mechanisms, implications of disease on shoulder complex and in assessing the quality of shoulder implants. Further, although few, Finite element shoulder models have also been utilized to answer important clinical questions such as the difference between a normal and osteoarthritic shoulder joint. However, due to the absence of experimental validation, it is questionable whether the constitutive models applied in these FE models are adequate to represent mechanical behaviors of shoulder elements (Cartilages, Ligaments, Muscles etc), therefore the confidence of using current models in answering clinically relevant question. The main objective of this review is to critically evaluate the existing FE shoulder models that have been used to investigate clinical problems. Due concern is given to check the adequacy of representative constitutive models of shoulder elements in drawing clinically relevant conclusion. Suggestions have been given to improve the existing shoulder models by inclusion of adequate constitutive models for shoulder elements to confidently answer clinically relevant questions.

Finite element analysis investigation of response of pumpkin tissue under uniaxial loading

Finite element analysis (FEA) models of uniaxial loading of pumpkin peel and flesh tissues were developed and validated using experimental results. The tensile model was developed for both linear elastic and plastic material models, the compression model was develop d only with the plastic material model. The outcomes of force versus time curves obtained from FEA models followed similar pattern to the experimental curves however the curve resulted with linear elastic material properties had a higher difference with the experimental curves. The values of predicted forces were determined and compared with the experimental curve. An error indicator was introduced and computed for each case and compared. Additionally Root Mean Square Error (RMSE) values were also calculated for each model and compared. The results of modelling were used to develop material model for peel and flesh tissues in FEA modelling of mechanical peeling of tough skinned vegetables.

Arterial luminal curvature and fibrous-cap thickness affect critical stress conditions within atherosclerotic plaque: An in vivo MRI-based 2D finite-element study

High mechanical stress in atherosclerotic plaques at vulnerable sites, called critical stress, contributes to plaque rupture. The site of minimum fibrous cap (FC) thickness (FCMIN) and plaque shoulder are well-documented vulnerable sites. The inherent weakness of the FC material at the thinnest point increases the stress, making it vulnerable, and it is the big curvature of the lumen contour over FC which may result in increased plaque stress. We aimed to assess critical stresses at FCMIN and the maximum lumen curvature over FC (LCMAX) and quantify the difference to see which vulnerable site had the highest critical stress and was, therefore, at highest risk of rupture. One hundred patients underwent high resolution carotid magnetic resonance (MR) imaging. We used 352 MR slices with delineated atherosclerotic components for the simulation study. Stresses at all the integral nodes along the lumen surface were calculated using the finite-element method. FCMIN and LCMAX were identified, and critical stresses at these sites were assessed and compared. Critical stress at FC MIN was significantly lower than that at LCMAX (median: 121.55 kPa; inter quartile range (IQR) = [60.70-180.32] kPa vs. 150.80 kPa; IQR = [91.39-235.75] kPa, p < 0.0001). If critical stress at FCMIN was only used, then the stress condition of 238 of 352 MR slices would be underestimated, while if the critical stress at LCMAX only was used, then 112 out of 352 would be underestimated. Stress analysis at FCMIN and LCMAX should be used for a refined mechanical risk assessment of atherosclerotic plaques, since material failure at either site may result in rupture.

Identifying vulnerable carotid plaques in vivo using high resolution magnetic resonance imaging-based finite element analysis

Object. Individuals with carotid atherosclerosis develop symptoms following rupture of vulnerable plaques. Biomechanical stresses within this plaque may increase vulnerability to rupture. In this report the authors describe the use of in vivo carotid plaque imaging and computational mechanics to document the magnitude and distribution of intrinsic plaque stresses. Methods. Ten (five symptomatic and five asymptomatic) individuals underwent plaque characterization magnetic resonance (MR) imaging. Plaque geometry and composition were determined by multisequence review. Intrinsic plaque stress profiles were generated from 3D meshes by using finite element computational analysis. Differences in principal (shear) stress between normal and diseased sections of the carotid artery and between symptomatic and asymptomatic plaques were noted. Results. There was a significant difference in peak principal stress between diseased and nondiseased segments of the artery (mean difference 537.65 kPa, p < 0.05). Symptomatic plaques had higher mean stresses than asymptomatic plaques (627.6 kPa compared with 370.2 kPa, p = 0.05), which were independent of luminal stenosis and plaque composition. Conclusions. Significant differences in plaque stress exist between plaques from symptomatic individuals and those from asymptomatic individuals. The MR imaging-based computational analysis may therefore be a useful aid to identification of vulnerable plaques in vivo.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.

A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.