scispace - formally typeset
Search or ask a question

Showing papers by "YuanTong Gu published in 2014"


Journal ArticleDOI
TL;DR: In this article, the authors present preliminary 3D numerical simulations of a high-density radial-inflow ORC turbine in sensible geothermal conditions, using the refrigerant R143a as the high density working fluid.

96 citations


Journal ArticleDOI
TL;DR: A finite difference method to discretize the time variable and obtain a semi-discrete scheme and the convergence order is verified from a numerical example which is presented to describe a fractal model of mobile/immobile transport process with different problem domains.

80 citations


Journal ArticleDOI
TL;DR: In this paper, microscopic changes that occur in plant food materials during drying significantly influence the macroscopic properties and quality factors of the dried food materials, and it is critical to study the...
Abstract: Microscopic changes that occur in plant food materials during drying significantly influence the macroscopic properties and quality factors of the dried food materials. It is critical to study the ...

67 citations


Journal ArticleDOI
TL;DR: In this paper, a new one-dimensional space-fractional Boussinesq equation is proposed, and two novel numerical methods with a nonlocal operator (using nodal basis functions) for the space-of-partitioning Bousseinq equation are derived.

63 citations


Journal ArticleDOI
TL;DR: In this paper, a Striped Point-In-Box (S-PIB) algorithm is proposed to overcome the shortcoming of PIB algorithm that caused by points distribution, and the two timeconsuming steps of neighbor searching and contact searching are integrated into one searching step.

61 citations


Journal ArticleDOI
TL;DR: In this article, the mechanical properties of graphynes (α-, β-, γ- and 6,6,12-graphynes) under uniaxial tension deformation at different temperatures and strain rates are studied using molecular dynamics simulations.
Abstract: Graphyne is an allotrope of graphene. The mechanical properties of graphynes (α-, β-, γ- and 6,6,12-graphynes) under uniaxial tension deformation at different temperatures and strain rates are studied using molecular dynamics simulations. It is found that graphynes are more sensitive to temperature changes than graphene in terms of fracture strength and Young's modulus. The temperature sensitivity of the different graphynes is proportionally related to the percentage of acetylenic linkages in their structures, with the α-graphyne (having 100% of acetylenic linkages) being most sensitive to temperature. For the same graphyne, temperature exerts a more pronounced effect on the Young's modulus than fracture strength, which is different from that of graphene. The mechanical properties of graphynes are also sensitive to strain rate, in particular at higher temperatures.

48 citations


Journal ArticleDOI
TL;DR: In this paper, a single plant cell was modeled with smoothed particle hydrodynamics (SPH) and a discrete element method (DEM) to study the basic micromechanics that govern the cellular structural deformations during drying.

44 citations


Journal ArticleDOI
TL;DR: In this paper, a mesh-free particle based model was developed which is fundamentally capable of handling extreme deformations of plant tissues during drying, which can be used to simulate tissues under excessive moisture content reductions incorporating cell wall wrinkling.
Abstract: Fundamental understanding on microscopic physical changes of plant materials is vital to optimize product quality and processing techniques, particularly in food engineering. Although grid-based numerical modelling can assist in this regard, it becomes quite challenging to overcome the inherited complexities of these biological materials especially when such materials undergo critical processing conditions such as drying, where the cellular structure undergoes extreme deformations. In this context, a meshfree particle based model was developed which is fundamentally capable of handling extreme deformations of plant tissues during drying. The model is built by coupling a particle based meshfree technique: Smoothed Particle Hydrodynamics (SPH) and a Discrete Element Method (DEM). Plant cells were initiated as hexagons and aggregated to form a tissue which also accounts for the characteristics of the middle lamella. In each cell, SPH was used to model cell protoplasm and DEM was used to model the cell wall. Drying was incorporated by varying the moisture content, the turgor pressure, and cell wall contraction effects. Compared to the state of the art grid-based microscale plant tissue drying models, the proposed model can be used to simulate tissues under excessive moisture content reductions incorporating cell wall wrinkling. Also, compared to the state of the art SPH-DEM tissue models, the proposed model better replicates real tissues and the cell-cell interactions used ensure efficient computations. Model predictions showed good agreement both qualitatively and quantitatively with experimental findings on dried plant tissues. The proposed modelling approach is fundamentally flexible to study different cellular structures for their microscale morphological changes at dehydration.

43 citations


Journal ArticleDOI
01 Oct 2014-Carbon
TL;DR: In this paper, the thermal transport management of monolayer graphene allotrope nanoribbons (size ∼20 × 4 nm2) by the modulation of their structures via molecular dynamics simulations was studied.

38 citations


Journal ArticleDOI
TL;DR: In this paper, the authors employ the reverse non-equilibrium molecular dynamics simulation to elucidate the effect of various faulted stacking layers on the thermal conductivity (TC) of silicon nanowires.
Abstract: Faulted stacking layers are ubiquitously observed during the crystal growth of semiconducting nanowires (NWs). In this paper, we employ the reverse non-equilibrium molecular dynamics simulation to elucidate the effect of various faulted stacking layers on the thermal conductivity (TC) of silicon (Si) NWs. We find that the stacking faults can greatly reduce the TC of the Si NW. Among the different stacking faults that are parallel to the NW's axis, the 9R polytype structure, the intrinsic and extrinsic stacking faults (iSFs and eSFs) exert more pronounced effects in the reduction of TC than the twin boundary (TB). However, for the perpendicularly aligned faulted stacking layers, the eSFs and 9R polytype structures are observed to induce a larger reduction to the TC of the NW than the TB and iSFs. For all considered NWs, the TC does not show a strong relation with the increasing number of faulted stacking layers. Our studies suggest the possibility of tuning the thermal properties of Si NWs by altering the crystal structure via the different faulted stacking layers.

36 citations


Journal ArticleDOI
TL;DR: In this paper, the thermal conductivity of the two-dimensional carbon nanotube (CNT)-based architecture, which can be constructed through welding of single-wall CNTs by electron beam, was reported.
Abstract: We reported the thermal conductivity of the two-dimensional carbon nanotube (CNT)-based architecture, which can be constructed through welding of single-wall CNTs by electron beam. Using large-scale nonequilibrium molecular dynamics simulations, the thermal conductivity is found to vary with different junction types due to their different phonon scatterings at the junction. The strong length and strain dependence of the thermal conductivity suggests an effective avenue to tune the thermal transport properties of the CNT-based architecture, benefiting the design of nanoscale thermal rectifiers or phonon engineering.

Journal ArticleDOI
TL;DR: In this article, the lattice Boltzmann method was used to investigate the flow past four square cylinders in an in-line square configuration, and the gap spacing g = s/d was set at 1, 3 and 6 and Reynolds number ranging from Re = 60 to 175.
Abstract: In this paper two-dimensional (2-D) numerical investigation of flow past four square cylinders in an in-line square configuration are performed using the lattice Boltzmann method. The gap spacing g = s/d is set at 1, 3 and 6 and Reynolds number ranging from Re = 60 to 175. We observed four distinct wake patterns: (i) a steady wake pattern (Re = 60 and g = 1); (ii) a stable shielding wake pattern (80 ≤ Re ≤ 175 and g = 1); (iii) a wiggling shielding wake pattern (60 ≤ Re ≤ 175 and g = 3) and (iv) a vortex shedding wake pattern (60 ≤ Re ≤ 175 and g = 6). At g = 1, the Reynolds number is observed to have a strong effect on the wake patterns. It is also found that at g = 1, the secondary cylinder interaction frequency significantly contributes for drag and lift coefficients signal. It is found that the primary vortex shedding frequency dominates the flow and the role of secondary cylinder interaction frequency almost vanish at g = 6. It is observed that the jet between the gaps strongly influenced the wake interaction for different gap spacing and Reynolds number combination. To fully understand the wake transformations the details vorticity contour visualization, power spectra of lift coefficient signal and time signal analysis of drag and lift coefficients also presented in this paper.

Journal ArticleDOI
TL;DR: The results show that CO2 can form weak interactions with the BC59 cage in its neutral state and the interactions can be enhanced significantly by introducing an extra electron to the system.
Abstract: Recently, the capture and storage of CO2 have attracted research interest as a strategy to reduce the global emissions of greenhouse gases. It is crucial to find suitable materials to achieve an efficient CO2 capture. Here we report our study of CO2 adsorption on boron-doped C60 fullerene in the neutral state and in the 1e−-charged state. We use first principle density functional calculations to simulate the CO2 adsorption. The results show that CO2 can form weak interactions with the BC59 cage in its neutral state and the interactions can be enhanced significantly by introducing an extra electron to the system.

Journal ArticleDOI
01 Mar 2014-Carbon
TL;DR: In this paper, the formation of periodic ripples on the surface of adjacent graphitic grains has been observed for the first time on highly ordered pyrolitic graphite by UHV-STM.

Journal ArticleDOI
TL;DR: In this article, the authors carried out a first-time investigation on the tensile properties of the hybrid structures with different dopants and found that with the presence of dopants, the hybrid structure usually exhibit lower yield strength, Young's modulus, and earlier yielding compared to that of a pristine hybrid structure.
Abstract: Doping is an effective approach that allows for the intrinsic modification of the electrical and chemical properties of nanomaterials. Recently, a graphene and carbon nanotube hybrid structure (GNHS) has been reported, which extends the excellent properties of carbon-based materials to three dimensions. In this paper, we carried out a first-time investigation on the tensile properties of the hybrid structures with different dopants. It is found that with the presence of dopants, the hybrid structures usually exhibit lower yield strength, Young’s modulus, and earlier yielding compared to that of a pristine hybrid structure. For dopant concentrations below 2.5% no significant reduction of Young’s modulus or yield strength could be observed. For all considered samples, the failure is found to initiate at the region where the nanotubes and graphene sheets are connected. After failure, monatomic chains are normally observed around the failure region. Dangling graphene layers without the separation of a residual CNT wall are found to adhere to each other after failure with a distance of about 3.4 A. This study provides a fundamental understanding of the tensile properties of the doped graphene–nanotube hybrid structures, which will benefit the design and also the applications of graphene-based hybrid materials.

Journal Article
TL;DR: A meshfree particle based model was developed which is fundamentally capable of handling extreme deformations of plant tissues during drying and can be used to simulate tissues under excessive moisture content reductions incorporating cell wall wrinkling.
Abstract: Fundamental understanding on microscopic physical changes of plant materials is vital to optimize product quality and processing techniques, particularly in food engineering. Although grid-based numerical modelling can assist in this regard, it becomes quite challenging to overcome the inherited complexities of these biological materials especially when such materials undergo critical processing conditions such as drying, where the cellular structure undergoes extreme deformations. In this context, a meshfree particle based model was developed which is fundamentally capable of handling extreme deformations of plant tissues during drying. The model is built by coupling a particle based meshfree technique: Smoothed Particle Hydrodynamics (SPH) and a Discrete Element Method (DEM). Plant cells were initiated as hexagons and aggregated to form a tissue which also accounts for the characteristics of the middle lamella. In each cell, SPH was used to model cell protoplasm and DEM was used to model the cell wall. Drying was incorporated by varying the moisture content, the turgor pressure, and cell wall contraction effects. Compared to the state of the art grid-based microscale plant tissue drying models, the proposed model can be used to simulate tissues under excessive moisture content reductions incorporating cell wall wrinkling. Also, compared to the state of the art SPH-DEM tissue models, the proposed model better replicates real tissues and the cell-cell interactions used ensure efficient computations. Model predictions showed good agreement both qualitatively and quantitatively with experimental findings on dried plant tissues. The proposed modelling approach is fundamentally flexible to study different cellular structures for their microscale morphological changes at dehydration.

Journal ArticleDOI
TL;DR: In this article, the deformation properties of single chondrocytes have a dependency on the strain-rates, and they have been shown to be governed by both the cellular cytoskeleton and the intracellular fluid.
Abstract: Based on the characterization by Atomic Force Microscopy, we report that the mechanical property of single chondrocytes has dependency on the strain-rates. By comparing the mechanical deformation responses and the Young's moduli of living and fixed chondrocytes at four different strain-rates, we explore the deformation mechanisms underlying this dependency property. We found that the strain-rate-dependent mechanical property of living cells is governed by both of the cellular cytoskeleton and the intracellular fluid when the fixed chondrocytes are mainly governed by their intracellular fluid, which is called the consolidation-dependent deformation behavior. Finally, we report that the porohyperelastic constitutive material model which can capture the consolidation-dependent behavior of both living and fixed chondrocytes is a potential candidature to study living cell biomechanics.

Journal ArticleDOI
TL;DR: In this article, a focused ion beam (FIB) was used to manipulate the selected Na2Ti3O7 NW over a hole drilled in an indium tin oxide substrate.
Abstract: We report on the mechanical properties of sodium titanate nanowires (Na2Ti3O7 NW) through a combination of bending experiments and theoretical analysis. Na2Ti3O7 NWs with lateral dimensions ranging from 20–700 nm were synthesized by a hydrothermal approach. A focused ion beam (FIB) was used to manipulate the selected Na2Ti3O7 NW over a hole drilled in an indium tin oxide substrate. After welding the nanowire, a series of bending tests was performed. It was observed that the Na2Ti3O7 NW exhibits a brittle behavior, and a nonlinear elastic deformation was observed before failure. By using the modified Euler–Bernoulli beam theory, such nonlinear elastic deformation is found to originate from a combination of surface effects and axial elongation (arising from the bending deformation). The effective Young's modulus of the Na2Ti3O7 NW was found to be independent of the wire length, and ranges from 21.4 GPa to 45.5 GPa, with an average value of 33 ± 7 GPa. The yield strength of the Na2Ti3O7 NW is measured at 2.7 ± 0.7 GPa.

Journal ArticleDOI
TL;DR: In this article, a coupled mesh-free particle-based two-dimensional (2D) model was developed to simulate microscale deformations of plant cells during drying, which can be used as a building block for advanced tissue models which are highly applicable for product and process optimizations in Food Engineering.
Abstract: Plant based dried food products are popular commodities in global market where much research is focused to improve the products and processing techniques. In this regard, numerical modeling is highly applicable and in this work, a coupled meshfree particle-based two-dimensional (2-D) model was developed to simulate microscale deformations of plant cells during drying. Smoothed Particle Hydrodynamics (SPH) was used to model the viscous cell protoplasm (cell fluid) by approximating it to an incompressible Newtonian fluid. The visco-elastic characteristic of the cell wall was approximated to a Neo-Hookean solid material augmented with a viscous term and modeled with a Discrete Element Method (DEM). Compared to a previous work [44] , this study proposes three model improvements: linearly decreasing positive cell turgor pressure during drying, cell wall contraction forces and cell wall drying. The improvements made the model more comparable with experimental findings on dried cell morphology and geometric properties such as cell area, diameter, perimeter, roundness, elongation and compactness. This single cell model could be used as a building block for advanced tissue models which are highly applicable for product and process optimizations in Food Engineering.

Journal ArticleDOI
TL;DR: In this article, a numerical model of the composite portable water-filled barriers (PWFBs) consisting of a plastic shell, steel frame, water and foam was developed and validated against results from full scale experimental tests.

Journal ArticleDOI
TL;DR: Based on the large-scale molecular dynamics simulation and the Euler-Bernoulli beam theory, the influence from different pre-existing surface defects on the bending properties of Ag nanowires (NWs) is studied in this paper.

Journal ArticleDOI
TL;DR: 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, and the MSnHS is a good model for exploration of mechanical deformation responses of single osteocytes and chondROcytes.
Abstract: The aim of this paper is to determine the strain-rate-dependent mechanical behavior of living and fixed osteocytes and chondrocytes, in vitro. First, 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.

Journal ArticleDOI
TL;DR: In this paper, the fluid-structure interaction of portable water-filled barriers (PWFBs) under vehicular impact using several methods was investigated. And the authors found that water can be treated with the FEA-NSM model for velocities higher than 80 km/h −1.
Abstract: Portable water-filled barriers (PWFB) are roadside structures used to separate moving traffic from work-zones. Numerical PWFB modelling is preferred in the design stages prior to actual testing. This paper aims to study the fluid–structure interaction of PWFB under vehicular impact using several methods. The strategy to treat water as non-structural mass was proposed and the errors were investigated. It was found that water can be treated with the FEA–NSM model for velocities higher than 80 km h −1 . However, the full SPH/FEA model is still the best treatment for water and necessary for lower impact velocities. The findings in this paper can be used as guidelines for modelling and designing PWFB.

Journal ArticleDOI
TL;DR: In this paper, the binding ability of two-dimensional monolayer carbon nanomaterial on actin by molecular simulation to understand their adhesive characteristics on F-actin cytoskeleton was analyzed.
Abstract: The biosafety of carbon nanomaterial needs to be critically evaluated with both experimental and theoretical validations before extensive biomedical applications. In this Letter, we present an analysis of the binding ability of two-dimensional monolayer carbon nanomaterial on actin by molecular simulation to understand their adhesive characteristics on F-actin cytoskeleton. The modelling results indicate that the positively charged carbon nanomaterial has higher binding stability on actin. Compared to crystalline graphene, graphene oxide shows higher binding influence on actin when carrying positive surface charge. This theoretical investigation provides insights into the sensitivity of actin-related cellular activities on carbon nanomaterial.

Journal ArticleDOI
TL;DR: In this article, the Euler-Bernoulli beam model can only reliably be used when the single actin filament is of the order of hundreds of nanometre scale.
Abstract: The mechanical vibration properties of single actin filaments from 50 to 288 nm are investigated by the molecular dynamics simulation in this study. The natural frequencies obtained from the molecular simulations agree with those obtained from the analytical solution of the equivalent Euler–Bernoulli beam model. Through the convergence study of the mechanical properties with respect to the filament length, it was found that the Euler–Bernoulli beam model can only be reliably used when the single actin filament is of the order of hundreds of nanometre scale. This molecular investigation not only provides the evidence for the use of the continuum beam model in characterising the mechanical properties of single actin filaments, but also clarifies the criteria for the effective use of the Euler–Bernoulli beam model.

Journal ArticleDOI
TL;DR: In this article, the authors explore the possibility of tailoring the resonance properties of a bilayer graphene sheet (GS) with interlayer sp3 bonds and show that the sp3 bond can either degrade or enhance the resonance of the resonator depending on their density and location.
Abstract: Graphene-based resonators are envisioned to build the ultimate limit of two-dimensional nanoelectromechanical system due to their ultrasensitive detection of mass, force, pressure and charge. However, such application has been greatly impeded by their extremely low quality factor. In the present work, we explore, using the large-scale molecular dynamics simulation, the possibility of tailoring the resonance properties of a bilayer graphene sheet (GS) with interlayer sp3 bonds. For the bilayer GS resonator with interlayer sp3 bonds, we discovered that the sp3 bonds can either degrade or enhance the resonance properties of the resonator depending on their density and location. It is found that the distribution of sp3 bonds only along the edges of either pristine or hydrogenated bilayer GS, leads to a greatly enhanced quality factor. A quality factor of ~1.18×105 is observed for a 3.07×15.31 nm2 bilayer GS resonator with sp3 bonds, which is more than 30 times larger comparing with that of a pristine bilayer GS. The present study demonstrates that the resonance properties of a bilayer GS resonator can be tuned by introducing sp3 bonds. This finding provides a useful guideline for the synthesis of the bilayer GS for its application as a resonator component.

Journal Article
TL;DR: In this article, a mesh-free-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.
Abstract: 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.

01 Oct 2014
TL;DR: In this paper, a numerical model of a composite portable water-filled barrier (PWFB) consisting of a plastic shell, steel frame, water and foam was developed and validated against results from full scale experimental tests.
Abstract: Portable water-filled barriers (PWFBs) are roadside appurtenances that are used to prevent errant vehicles from penetrating into temporary construction zones on roadways A numerical model of the composite PWFB, consisting of a plastic shell, steel frame, water and foam was developed and validated against results from full scale experimental tests This model can be extended to larger scale impact cases, specifically ones that include actual vehicle models The cost-benefit of having a validated numerical model is significant and this allows the road barrier designer to conduct extensive tests via numerical simulations prior to standard impact tests Effects of foam cladding as additional energy absorption material in the PWFB was investigated Different types of foam were treated and it was found that XPS foam was the most suitable foam type Results from this study will aid PWFB designers in developing new generation of roadside structures which will provide enhanced road safety

Journal ArticleDOI
TL;DR: In this article, the effects of sudden differential heating and its resultant transient behavior on fluid flow and heat transfer characteristics have been shown for the range of governing parameters for the Rayleigh numbers ranging from 100 to 10,000.
Abstract: Numerically investigation of free convection within a porous cavity with differential heating has been performed using modified corrugated side walls. Sinusoidal hot left and cold right walls are assumed to receive sudden differentially heating where top and bottom walls are insulated. Air is considered as working fluid and is quiescent, initially. Numerical experiments reveal 3 distinct stages of developing pattern including initial stage, oscillatory intermediate, and finally steady-state condition. Implicit Finite Volume Method with TDMA solver is used to solve the governing equations. This study has been performed for the Rayleigh numbers ranging from 100 to 10,000. Outcomes have been reported in terms of isotherms, streamline, velocity and temperature plots and average Nusselt number for various Ra, corrugation frequency, and corrugation amplitude (CA). The effects of sudden differential heating and its resultant transient behavior on fluid flow and heat transfer characteristics have been shown for the range of governing parameters. The present results show that the transient phenomena are enormously influenced by the variation of the Rayleigh Number with CA and frequency.

Journal ArticleDOI
TL;DR: The aim of this study is to examine the strain-rate dependent behavior of a single chondrocyte using a computational approach based on finite element method (FEM).
Abstract: Various studies have been conducted to investigate the effects of impact loading on cartilage damage and chondrocyte death. These have shown that the rate and magnitude of the applied strain significantly influence chondrocyte death, and that cell death occurred mostly in the superficial zone of cartilage suggesting the need to further understand the fundamental mechanisms underlying the chondrocytes death induced at certain levels of strain-rate. To date there is no comprehensive study providing insight on this phenomenon. The aim of this study is to examine the strain-rate dependent behavior of a single chondrocyte using a computational approach based on finite element method (FEM). An FEM model was developed using various mechanical models, which were standard Neo-Hookean solid (SnHS), porohyperelastic (PHE) and poroviscohyperelastic (PVHE) to simulate atomic force microscopy (AFM) experiments of chondrocyte. The PVHE showed, it can capture both relaxation and loading rate dependent behaviors of chondrocytes, accurately compared to other models.