This study aimed to develop a three-dimensional finite element model of a functionally graded femoral prosthesis. The model consisted of a femoral prosthesis created from functionally graded materials (FGMs), cement, and femur. The hip prosthesis was composed of FGMs made of titanium alloy, chrome-cobalt, and hydroxyapatite at volume fraction gradient exponents of 0, 1, and 5, respectively. The stress was measured on the femoral prosthesis, cement, and femur. Stress on the neck of the femoral prosthesis was not sensitive to the properties of the constituent material. However, stress on the stem and cement decreased proportionally as the volume fraction gradient exponent of the FGM increased. Meanwhile, stress became uniform on the cement mantle layer. In addition, stress on the femur in the proximal part increased and a high surface area of the femoral part was involved in absorbing the stress. As such, the stress-shielding area decreased. The results obtained in this study are significant in the design and longevity of new prosthetic devices because FGMs offer the potential to achieve stress distribution that more closely resembles that of the natural bone in the femur.
A kinetic model incorporating adsorption, desorption and biodegradation processes was developed to describe the bioregeneration of granular activated carbon (GAC) loaded with 4-chlorophenol (4-CP) and 2,4-dichlorophenol (2,4-DCP), respectively, in simultaneous adsorption and biodegradation processes. The model was numerically solved and the results showed that the kinetic model was well-fitted (R(2)>0.83) to the experimental data at different GAC dosages and at various initial 4-CP and 2,4-DCP concentrations. The rate of bioregeneration in simultaneous adsorption and biodegradation processes was influenced by the ratio of initial chlorophenol concentration to GAC dosage. Enhancement in the rate of bioregeneration was achieved by using the lowest ratio under either one of the following experimental conditions: (1) increasing initial chlorophenol concentration at constant GAC dosage and (2) increasing GAC dosage at constant initial chlorophenol concentration. It was found that the rate enhancement was more pronounced under the second experimental condition.
The wrist is the most complex joint for virtual three-dimensional simulations, and the complexity is even more pronounced when dealing with skeletal disorders of the joint such, as rheumatoid arthritis (RA). In order to analyse the biomechanical difference between healthy and diseased joints, three-dimensional models of these two wrist conditions were developed from computed tomography images. These images consist of eight carpal bones, five metacarpal bones, the distal radius and ulna. The cartilages were developed based on the shape of the available articulations and ligaments were simulated via mechanical links. The RA model was developed accurately by simulating all ten common criteria of the disease related to the wrist. Results from the finite element (FE) analyses showed that the RA model produced three times higher contact pressure at the articulations compared to the healthy model. Normal physiological load transfer also changed from predominantly through the radial side to an increased load transfer approximately 5% towards the ulnar. Based on an extensive literature search, this is the first ever reported work that simulates the pathological conditions of the rheumatoid arthritis of the wrist joint.
Anaerobic treatment processes to remove organic matter from palm oil mill effluent (POME) have been used widely in Malaysia. Still the amounts of total organic and total mineral released from POME that may cause degradation of the receiving environment need to be verified. This paper proposes the use of the hydrodynamic equations to estimate performance of the cascaded anaerobic ponds (CAP) and to calculate amounts of total organic matter and total mineral released from POME. The CAP efficiencies to remove biochemical oxygen demands, chemical oxygen demands, total solids and volatile solids (VS) as high as 94.5, 93.6, 96.3 and 98.2 %, respectively, are estimated. The amounts of total organic matter and total mineral as high as 538 kg VS/day and 895 kg FS/day, respectively, released from POME to the receiving water are calculated. The implication of the proposed hydrodynamic equations contributes to more versatile environmental assessment techniques, sometimes replacing laboratory analysis.
A femur phantom made of wax and a real human bone was used to study the dose during radiographical procedures. The depth dose inside the phantom was determined using DOSXYZnrc, a Monte Carlo simulation software. The results were verified with measurements using TLD-100H. It was found that for 2.5 mm aluminium filtered 84-kVp X-rays, the radiation dose in the bone reached 57 % higher than the surface dose, i.e. 3.23 mGy as opposed to 2.06 mGy at the surface. The use of real bone introduces variations in the bone density in the DOSXYZnrc model, resulting in a lower attenuation effect than expected from solid bone tissues.
Optimization of thermo-alkaline disintegration of sewage sludge for enhanced biogas yield was carried out using response surface methodology (RSM) and Box-Behnken design of experiment. The individual linear and quadratic effects as well as the interactive effects of temperature, NaOH concentration and time on the degree of disintegration were investigated. The optimum degree of disintegration achieved was 61.45% at 88.50 °C, 2.29 M NaOH (24.23%w/w total solids) and 21 min retention time. Linear and quadratic effects of temperature are most significant in affecting the degree of disintegration. The coefficient of determination (R(2)) of 99.5% confirms that the model used in predicting the degree of disintegration process has a very good fitness with the experimental variables. The disintegrated sludge increased the biogas yield by 36%v/v compared to non-disintegrated sludge. The RSM with Box-Behnken design is an effective tool in predicting the optimum degree of disintegration of sewage sludge for increased biogas yield.
One of the available tools for mapping the geographical distribution and potential suitable habitats is species distribution models. These techniques are very helpful for finding poorly known distributions of species in poorly sampled areas, such as the tropics. Maximum Entropy (MaxEnt) is a recently developed modeling method that can be successfully calibrated using a relatively small number of records. In this research, the MaxEnt model was applied to describe the distribution and identify the key factors shaping the potential distribution of the vulnerable Malayan Sun Bear (Helarctos malayanus) in one of the main remaining habitats in Peninsular Malaysia. MaxEnt results showed that even though Malaysian sun bear habitat is tied with tropical evergreen forests, it lives in a marginal threshold of bio-climatic variables. On the other hand, current protected area networks within Peninsular Malaysia do not cover most of the sun bears potential suitable habitats. Assuming that the predicted suitability map covers sun bears actual distribution, future climate change, forest degradation and illegal hunting could potentially severely affect the sun bear's population.
The spinal facet joints are known to be an important component in the kinematics and the load transmission of the spine. The articular cartilage in the facet joint is prone to degenerative changes which lead to back pain and treatments for the condition have had limited long term success. There is currently a lack of information on the basic biomechanical properties of the facet joint cartilage which is needed to develop tissue substitution or regenerative interventions. In the present study, the thickness and biphasic properties of ovine facet cartilage were determined using a combination of indentation tests and computational modelling. The equilibrium biphasic Young's modulus and permeability were derived to be 0.76±0.35 MPa and 1.61±1.10×10⁻¹⁵ m⁴/(Ns) respectively, which were within the range of cartilage properties characterised from the human synovial joints. The average thickness of the ovine facet cartilage was 0.52±0.10 mm, which was measured using a needle indentation test. These properties could potentially be used for the development of substitution or tissue engineering interventions and for computational modelling of the facet joint. Furthermore, the developed method to characterise the facet cartilage could be used for other animals or human donors.
This study was conducted to characterize the growth of and aflatoxin production by Aspergillus flavus on paddy and to develop kinetic models describing the growth rate as a function of water activity (a(w)) and temperature.
This study uses computer simulations to study the impact of a motorcycle with the conventional w-beam guardrail. A three-dimensional computer simulation of a scaled hybrid III 50th-percentile male dummy mounted on a motorcycle and colliding with a w-beam guardrail is carried out. A multi-body model of the motorcycle and finite element model of the guardrail are developed using commercially available software. The simulation model is validated with a physical crash test conducted with same initial impact configurations. Impacts at speeds of 32, 48, and 60 km/h at an impact angle at 45 degrees are considered. The predicted forces and accelerations are compared with the biomechanical limits for each body part and the risk of injury to the rider are evaluated. Speed was found to have a significant influence on the level of injury to the head, neck, chest, and femur. A significant reduction of the severity of injuries was found when the impact speed was reduced from 60 to 32km/h. The accelerations experienced by the head and chest are found to be higher than safe levels for impact speeds of 48 and 60 km/h. The biomechanical limit for the right femur is exceeded at all three considered impact speeds. Neck injuries are also a concern, with the predicted tension values and neck bending extent being higher than the biomechanical limit for the 60 km/h impact speed. In light of these results, it is suggested that the design of guardrails should be reviewed with a focus on the safety of motorcyclists.
Chemical inactivation of microorganisms is a common process widely employed in many fields such as in treatment of water, preservation in food industry and antimicrobial treatments in healthcare. For economy of applications and efficiency of treatment establishment the minimum dosage of breakpoint in the chemical application becomes essential. Even though experimental investigations have been extensive, theoretical understanding of such processes are demanding. Commonly employed theoretical analyses for the inactivation of microorganisms and depletion of chemicals include kinetics expressing the rates of depletion of chemical and microorganisms. The terms chemical demand (x) and specific disinfectant demand (alpha) are often used in theoretical modeling of inactivation. The value of specific disinfectant demand (alpha) has always been assumed to be a constant in these models. Intracellular concentration built up within the cells of the microorganisms during inactivation could lead to possible weakening effects of microorganisms thereby requiring lower doses as disinfection proceeds makes the assumption of constant alpha inaccurate. Model equations are formulated based on these observations co-relating the parameters alpha and x with a progressive inactivation (N/N(0)). The chemical concentration (C) is also presented in terms of the inactivation time (t) and the survival ratio (N/N(0)) for given pH and temperature conditions. The model is examined using experimentally verified Ct data of Giardia Cysts/chlorine system. The respective values of x for different survival ratios were evaluated from the data using MatLab software. Proposed model correlating for the disinfectant demand (x) with the survival ratio (N/N(0)) fits satisfactorily with those evaluated from data. The rate constants for different pH and temperature conditions are evaluated which showed compatibility with the Arrhenius model. The dependence of frequency factors with pH indicated compatibility with accepted models. The Ct values regenerated with the kinetic data shows a very accurate fit with published data.
Activated sludge models (ASMs) have been widely used as a basis for further model development in wastewater treatment processes. Values for parameters to be used are vital for the accuracy of the modeling approach. A continuous stirred tank reactor (CSTR), as open respirometer with continuous flow for 20 h is used in ASMs. The dissolved oxygen (DO) profile for 11 days was monitored. It was found the mass transfer coefficient K(La) is 0.3 h(-1) during lag and start feed phase and 0.01 h(-1) during stop feed phase, while the heterotrophic yield coefficient Y(H) is 0.44. Some of the chemical oxygen demand (COD) fractionations of palm oil mill effluent (POME) using respirometric test in ASM models are S(s) 50 mg/L, S(I) 16,600 mg/L, X(S) 25,550 mg/L, and X(I) 2,800 mg/L. The comparison of experimental and ASM1 from OUR concentration is found to fit well.
This study aimed to develop an optimal continuous procedure of lipase-catalyzes transesterification of waste cooking palm oil in a packed bed reactor to investigate the possibility of large scale production further. Response surface methodology (RSM) based on central composite rotatable design (CCRD) was used to optimize the two important reaction variables packed bed height (cm) and substrate flow rate(ml/min) for the transesterification of waste cooking palm oil in a continuous packed bed reactor. The optimum condition for the transesterification of waste cooking palm oil was as follows: 10.53 cm packed bed height and 0.57 ml/min substrate flow rate. The optimum predicted fatty acid methyl ester (FAME) yield was 80.3% and the actual value was 79%. The above results shows that the RSM study based on CCRD is adaptable for FAME yield studied for the current transesterification system. The effect of mass transfer in the packed bed reactor has also been studied. Models for FAME yield have been developed for cases of reaction control and mass transfer control. The results showed very good agreement compatibility between mass transfer model and the experimental results obtained from immobilized lipase packed bed reactor operation, showing that in this case the FAME yield was mass transfer controlled.
Chromatin morphologies in human breast cancer cells treated with an anti-cancer agent are analyzed at their early stage of programmed cell death or apoptosis. The gray-level images of nuclear chromatin are modelled as random fields. We used two-dimensional isotropic generalized Cauchy field to characterize local self-similarity and global long-range dependence behaviors in the image spatial data. Generalized Cauchy field allows the description of fractal behavior inferred from fractal dimension and the long-range dependence inferred from correlation exponent to be carried out independently. We demonstrated the usefulness of locally self-similar random fields with long-range dependence for modelling chromatin condensation.
A sequential optimization based on statistical design and one-factor-at-a-time (OFAT) method was employed to optimize the media constituents for the improvement of citric acid production from oil palm empty fruit bunches (EFB) through solid state bioconversion using Aspergillus niger IBO-103MNB. The results obtained from the Plackett-Burman design indicated that the co-substrate (sucrose), stimulator (methanol) and minerals (Zn, Cu, Mn and Mg) were found to be the major factors for further optimization. Based on the OFAT method, the selected medium constituents and inoculum concentration were optimized by the central composite design (CCD) under the response surface methodology (RSM). The statistical analysis showed that the optimum media containing 6.4% (w/w) of sucrose, 9% (v/w) of minerals and 15.5% (v/w) of inoculum gave the maximum production of citric acid (337.94 g/kg of dry EFB). The analysis showed that sucrose (p<0.0011) and mineral solution (p<0.0061) were more significant compared to inoculum concentration (p<0.0127) for the citric acid production.
In 2008, the IUCN threat status of the Asian tapir (Tapirus indicus) was reclassified from 'vulnerable' to 'endangered'. The latest distribution map from the IUCN Red List suggests that the tapirs' native range is becoming increasingly fragmented in Peninsular Malaysia, but distribution data collected by local researchers suggest a more extensive geographical range. Here, we compile a database of 1261 tapir occurrence records within Peninsular Malaysia, and demonstrate that this species, indeed, has a much broader geographical range than the IUCN range map suggests. However, extreme spatial and temporal bias in these records limits their utility for conservation planning. Therefore, we used maximum entropy (MaxEnt) modeling to elucidate the potential extent of the Asian tapir's occurrence in Peninsular Malaysia while accounting for bias in existing distribution data. Our MaxEnt model predicted that the Asian tapir has a wider geographic range than our fine-scale data and the IUCN range map both suggest. Approximately 37% of Peninsular Malaysia contains potentially suitable tapir habitats. Our results justify a revision to the Asian tapir's extent of occurrence in the IUCN Red List. Furthermore, our modeling demonstrated that selectively logged forests encompass 45% of potentially suitable tapir habitats, underscoring the importance of these habitats for the conservation of this species in Peninsular Malaysia.
The development of accurate computational models of biological processes is fundamental to computational systems biology. These models are usually represented by mathematical expressions that rely heavily on the system parameters. The measurement of these parameters is often difficult. Therefore, they are commonly estimated by fitting the predicted model to the experimental data using optimization methods. The complexity and nonlinearity of the biological processes pose a significant challenge, however, to the development of accurate and fast optimization methods. We introduce a new hybrid optimization method incorporating the Firefly Algorithm and the evolutionary operation of the Differential Evolution method. The proposed method improves solutions by neighbourhood search using evolutionary procedures. Testing our method on models for the arginine catabolism and the negative feedback loop of the p53 signalling pathway, we found that it estimated the parameters with high accuracy and within a reasonable computation time compared to well-known approaches, including Particle Swarm Optimization, Nelder-Mead, and Firefly Algorithm. We have also verified the reliability of the parameters estimated by the method using an a posteriori practical identifiability test.
The healing process of ruptured tendons is problematic due to scar tissue formation and deteriorated material properties, and in some cases, it may take nearly a year to complete. Mechanical loading has been shown to positively influence tendon healing; however, the mechanisms remain unclear. Computational mechanobiology methods employed extensively to model bone healing have achieved high fidelity. This study aimed to investigate whether an established hyperelastic fibre-reinforced continuum model introduced by Gasser, Ogden and Holzapfel (GOH) can be used to capture the mechanical behaviour of the Achilles tendon under loading during discrete timepoints of the healing process and to assess the model's sensitivity to its microstructural parameters. Curve fitting of the GOH model against experimental tensile testing data of rat Achilles tendons at four timepoints during the tendon repair was used and achieved excellent fits ([Formula: see text]). A parametric sensitivity study using a three-level central composite design, which is a fractional factorial design method, showed that the collagen-fibre-related parameters in the GOH model-[Formula: see text] and [Formula: see text]-had almost equal influence on the fitting. This study demonstrates that the GOH hyperelastic fibre-reinforced model is capable of describing the mechanical behaviour of healing tendons and that further experiments should focus on establishing the structural and material parameters of collagen fibres in the healing tissue.
The epidemiology of dengue fever is characterized by highly seasonal, multi-annual fluctuations, and the irregular circulation of its four serotypes. It is believed that this behaviour arises from the interplay between environmental drivers and serotype interactions. The exact mechanism, however, is uncertain. Constraining mathematical models to patterns characteristic to dengue epidemiology offers a means for detecting such mechanisms. Here, we used a pattern-oriented modelling (POM) strategy to fit and assess a range of dengue models, driven by combinations of temporary cross protective-immunity, cross-enhancement, and seasonal forcing, on their ability to capture the main characteristics of dengue dynamics. We show that all proposed models reproduce the observed dengue patterns across some part of the parameter space. Which model best supports the dengue dynamics is determined by the level of seasonal forcing. Further, when tertiary and quaternary infections are allowed, the inclusion of temporary cross-immunity alone is strongly supported, but the addition of cross-enhancement markedly reduces the parameter range at which dengue dynamics are produced, irrespective of the strength of seasonal forcing. The implication of these structural uncertainties on predicted vulnerability to control is also discussed. With ever expanding spread of dengue, greater understanding of dengue dynamics and control efforts (e.g. a near-future vaccine introduction) has become critically important. This study highlights the capacity of multi-level pattern-matching modelling approaches to offer an analytic tool for deeper insights into dengue epidemiology and control.