A Modeling Study of Notch Noise Responses of Type III Units in the Gerbil Dorsal Cochlear Nucleus A computational model of the neural circuitry of the gerbil dorsal cochlear nucleus -LRB- DCN -RRB- , based on the MacGregor 's neuromime model , was used to simulate type III unit -LRB- P-cell -RRB- responses to notch noise stimuli . The DCN patch model is based on a previous computational model of the cat DCN -LSB- Hancock , K. E. , and H. F. Voigt . Ann . Biomed . Eng . 27:73 -- 87 , 1999 -RSB- . According to the experimental study of Parsons et al. -LSB- Ann . Biomed . Eng . 29:887 -- 896 , 2001 -RSB- , the responses of gerbil DCN type III units to notch noise stimuli are similar to those of cat DCN type IV units , which are thought to be spectral notch detectors . This suggests that type III units in the gerbil DCN may serve as spectral notch detectors . In this modeling study , a simplified notch noise response plot -- spike discharge rate vs. notch cutoff frequency plot -- was used to compare model responses to the experimental results . Parameter estimation and sensitivity analysis of three connection parameters within the DCN patch have been studied and shows the model is robust , providing reasonable fits to the experimental data from 14 of 15 type III units examined -LSB- work supported by a grant from NIDCD , Boston University 's Biomedical Engineering department and Hearing Research Center -RSB- . INTRODUCTION The cochlear nucleus -LRB- CN -RRB- is the sole target of auditory nerve fibers and is the first brainstem nucleus in the ascending auditory pathway . The dorsal part of the cochlear nucleus -LRB- DCN -RRB- is one of its three sub-nuclei and has initial integration functions . The DCN 's role in audition appears , in part , to be involved in sound localization in the median plane equidistant from the two ears where the traditional sound source location cues , interaural time and level differences , are absent . For broadband signals , the pinna attenuates energy in a narrow band whose center frequency can be correlated with sound source elevation . The pinnae-induced attenuation is called a spectral notch . These notches are seen in the head-related transfer functions -LRB- HRTFs -RRB- , which are the Fourier transforms of the impulse responses from free space sound sources to the tympanic membrane . These spectral features are found in human , cat and gerbil HRTFs .12,13 Principal cells of the DCN , type III units in the gerbil15 and type IV units in the cat ,22 exhibit sharp sensitivity to notch noise signals .14,18 This suggests that the DCN is a spectral notch detector across different species . The previous study of Hancock and Voigt6 used a computational model of the DCN to model the responses of three type IV units to notch noises and presented close fits to the physiological data recorded by Spirou and Young .18 The gerbil DCN has fewer type IV units than the cat DCN -LRB- 11 % type IV units in gerbil vs. 31 % in cat -RRB- 2,17 and is comprised primarily of type III units -LRB- 62 % -RRB- .2 As the principal units in gerbil DCN , type III units undoubtedly play an important role in gerbil sound localization . In this study , we used the same computational model to fit the physiological data obtained from the gerbil DCN by Parsons et al. 15 Despite the fact that type III units have quite different responses to tones than type IV units , this study shows that by modifying the connectivities within the neuron groups we can simulate both type III and type IV units as notch detectors . Also in this work we have good fits for 14 out of 15 type III units in gerbil ; this is substantially more modeling data than previously presented and points to the robustness of the DCN model . RESPONSE MAP SCHEME AND RATE VS. CUTOFF FREQUENCY PLOTS A response map -LRB- RM -RRB- is a plot of excitatory and inhibitory responses to tonal stimuli shown in a stimulus frequency vs. sound pressure level plane . It is a popular method to categorize the neurons in the CN of unanethetized , decerebrate animals .4,17,22,23 Type III units are characterized by their RMs -LRB- Fig. 1A -RRB- , which show a center V-shape excitation region flanked by side-band inhibition . The tip of the excitation region points to the unit 's best frequency -LRB- BF -RRB- , the frequency to which the unit was most sensitive . In Fig. 1A , there are seven discharge-rate vs. frequency curves shown for seven sound pressure levels presented . RMs of type IV units typically show an excitatory region over BF and an inhibitory area above that ranging over a broad frequency range ; sometimes a narrow excitatory area appears off BF -LRB- see Fig. 1 of Ref . 6 -RRB- . FIGURE 1 . Physiology of a type III unit : -LRB- A -RRB- Response map composed of rate vs. frequency curves of unit L021199 1.01 in Parsons et al. 15 . Excitatory responses are shaded black while inhibitory responses left blank . -LRB- B -RRB- Rate-level curves when given wideband noise and notch noise stimuli . Only four notch noise responses with notch -LRB- nw -RRB- widths 1.0 kHz , 2.0 kHz , 4.0 kHz , and 8.0 kHz are shown . -LRB- C -RRB- Rate vs. Cutoff frequency plot . Here notch widths varies systematically from 0.5 to 8.0 kHz in step of 0.5 kHz . The average rates are calculated from the range of 46 -- 54 dB SPL as shown by black bar in -LRB- B -RRB- . Responses to broadband noise , in addition to tones , are used to classify units in the DCN . We also used the same notch noise -LRB- NN -RRB- stimuli in the model simulations as in the physiological experiments .15 These have band reject notches of 30 dB centered at the BF with varying widths . In Fig. 1B , discharge-rate vs. level curves are plotted for broadband noise and four NN stimuli with widths of 1 , 2 , 4 , and 8 kHz . These responses were taken from a series of responses to BF-centered NN stimuli whose widths systematically varied from 0.5 to 8.5 kHz with a 0.5 kHz step . For Unit L021199 1.01 in Fig. 1 there were 16 NN stimuli applied with widths from 0.5 to 8.0 kHz . Only four NN responses are plotted for clarity . To summarize the information in these rate-level curves , Spirou and Young18 further characterized NN responses with rate vs. cutoff frequency plots -LRB- Fig. 1C -RRB- , where the firing rates to each NN were averaged over a range of sound levels from 46 to 54 dB SPL , where the inhibition reached the maximum , as indicated with black bar on the abscissa in Fig. 1B . The cutoff frequencies included both the notch noise 's lowpass and the highpass cutoff frequencies , thus the average rate of each notch noise response is plotted twice . In such plots , 11 of 15 units in Parsons et al. 's experimental data showed a specific `` Mexican Hat '' shape , indicating that when increasing the notch width , the average response rate dropped sharply from excitatory to inhibitory and then showed some degree of recovery . The unit in the figure was excited by the broad-band noise and NN with narrow-notch width -LRB- 0.5 kHz , not shown in Fig. 1B but shown in Fig. 1C -RRB- , but was inhibited by NN with wider widths -LRB- from 1.0 to 6.0 kHz , also shown in Fig. 1C -RRB- . Also it showed some degree of recovery to a less inhibitory stage when NN with an even wider notch was applied -LRB- from 6.5 to 8.0 kHz , shown in Fig. 1C -RRB- . This recovery was seen in 11 of 15 units in the paper of Parsons et al. 15 Similar responses to NN have been found in type IV units in cat by Spirou and Young .18 This demonstrated the sensitivity of type III units in gerbils to the NNs . Thus , type III units in gerbils are capable of detecting spectral notches as do type IV units in cats . FIGURE 2 . Conceptual model of gerbil dorsal cochlear nucleus circuitry . Only inhibitory connections are shown . The size of the triangles indicates connection strength . The horizontal axis represents the input from the auditory nerve fibers with increasing frequency , which target all cells . Modified from Davis et al. 2 CONCEPTUAL MODEL OF DCN CIRCUITRY To illustrate the functional role of type III units in the gerbil DCN , the conceptual model from Davis et al. 2 was chosen -LRB- see Fig. 2 -RRB- . The gerbil DCN conceptual model is more complex than the cat model since the principal cells in gerbils can be either type III units or type IV units , including subtypes of type IV units : type IV-i and type IV-T units . In addition , besides type II units and the wide band inhibitor -LRB- WBI -RRB- units , type III-i units were included in the circuitry to provide necessary inhibitory input to type IV-i units . In this paper , however , we only discuss the properties and responses of type III units , thus only type II units and WBI units are taken as inhibitors . Auditory nerve -LRB- AN -RRB- fibers provide excitatory inputs to type II units , type III units and WBI units . However , DCN units show spontaneous activity even after destruction of the cochlea .9 Thus non-specific afferents -LRB- NSA -RRB- were added to the principal cells in our computational model to provide noncochlear inputs to the DCN -LRB- Fig. 3 -RRB- . In order not to confuse physiologically defined unit types and model unit types , we use different names to distinguish them -LRB- see Table 1 -RRB- . FIGURE 3 . A patch of DCN model . Five cell populations are organized in a tonotopic manner with frequency step equal to 0.005 octaves . The excitatory and inhibitory connections are shown in a single frequency slice . P cells , representing type III units in gerbil and type IV units in cat ; I2 cells , type II units ; W cells , providing wide band inhibition ; AN fibers , the auditory nerve fiber input . Modified From Hancock and Voigt , 1999.6 Fusiform cells in gerbil DCN are mostly type III units in intracellular recording and marking experiments by Hancock and Voigt7 whereas type II units arise from vertical cells in the DCN21 and WBI units are thought to project from the posteroventral CN .20 In addition , WBI units send inhibitory connections to both type IV units -LRB- principal units in cat -RRB- and type II units in cat DCN . We are assuming that these connections to the gerbil principal units remain intact.TABLE 1 . Physiological unit types with corresponding model unit types.Physiological unit typeModel unit typeType III unitsP-cellsType II unitsI2-cellsWideband inhibitorsW-cellsAuditory nerve fibersAN-fibers METHODS Organization of the DCN Model The gerbil model is based on a previous computational model of the cat DCN6 -LRB- see Fig. 3 -RRB- . The model contains five cell groups that are arranged into 800 isofrequency slices centered at 5 kHz with two octaves below and above in intervals of 0.005 octaves . Figure 3 shows the connections among the model cells within a single frequency slice , where P-cells represent the principal cells , I2-cells represent interneuron inhibitors with type II responses and W-cells represent wide band inhibitors . AN-fibers excite P-cells , I2-cells and W-cells . W-cells inhibit P-cells and I2-cells . I2-cells inhibit P-cells . P-cells also are excited by non-specific afferents . The relationship among the model cells can be discussed in terms of the connection parameters as shown in Fig. 4 . For a given population of model cells -LRB- e.g. , P-cells -RRB- , the target cell `` B '' will receive N inputs from the `` A '' source cells , which are chosen from a band of `` A '' source cells with bandwidth BW , offset from the BF of B by C . The subscript indicates the source unit , A to the target unit , B . In the simulations shown here , all the P-cells ' BFs will be 5 kHz.FIGURE 4 . Group connection parameters . B is the target cell and a group of source cells A connect to it . and are the step increase and time constant in response to the input spikes . Here B cell receives N inputs from A cells , which are randomly chosen from a band of bandwidth and center frequency offset specified cells.TABLE 2 . Intrinsic parameters.Populationτm -LRB- ms -RRB- θ -LRB- mV -RRB- bkτk -LRB- ms -RRB- Ek -LRB- mV -RRB- Eex -LRB- mV -RRB- Ein -LRB- mV -RRB- W-cells5.04 .02.001.0 − 10 +70 − 10I2-cells6.014 .51.751.0 − 10 +70 − 10P-cells10.07 .52.001.0 − 10 +70 − 10 Neuron Model The neuron model -LRB- Fig. 5 -RRB- is based on the MacGregor neuromime model ,10 which is a parallel circuit model containing the membrane capacitance , leakage conductance , a potassium channel branch and the excitatory/inhibitory connection branches . Each excitatory or inhibitory input will add a branch , with the variable conductance controlled by the parameters , step σ and time constant τ.FIGURE 5 . -LRB- A -RRB- Vm is the membrane potential relative to rest ; Cm is membrane capacitance ; G is resting conductance ; Gk and Ek are variable conductance and reversal potential of potassium ; Gex/in and Eex/in are excitatory/inhibitory synaptic conductance and reversal potential . -LRB- B -RRB- Expression of connection synapse conductance . Target cell will undergo a step increase and exponential decay , time constant in response to input spikes . Event time is recorded when the membrane potential exceeds its threshold , θ , and the potassium conductance is activated to induce the cell 's refractory period . Thus , the computation is reduced by leaving out the details of action potential generation . The following equations describe the neuronmime model : where gk = Gk/G , gex = Gex/G , and gin = Gin/G ; G , Gk , Gex , and Gin represent respectively the conductances of neuron membrane , potassium channel , excitatory synaptic inputs and inhibitory synaptic inputs . The latter three are normalized by G to simplify the equations . Vm is the membrane potential . Ek , Eex , and Ein represent the reversal potentials . τm and τk are the membrane time constant and the refractory time constant . Equation -LRB- 1 -RRB- states the current fluxes across the neuronal membrane . Equation -LRB- 2 -RRB- states the change of the potassium conductance , where bk represents the sensitivity to potassium conductance . S is the spiking variable with two alternative values : 1 indicates a cell is firing.TABLE 3 . Connection parameters.Connection -LRB- octave -RRB- -LRB- octave -RRB- -LRB- ms -RRB- 0.0 Varies1400 .06100.00.4480.55100.00.448 Varies100 .00.1151.40100 / 0.20.115 Varies10 − 0.2 / 00.6210.11 n/an/a150 .153 For target cell B , the variable conductance that represents synapses from source cells A is described by where SA are the input spikes . Parameter Values In the model , we used four intrinsic -LRB- τm , θ , bk , and τk -RRB- and five connection parameters -LRB- , , , , and -RRB- . The intrinsic parameter set describes the membrane properties of the cell , which determines the steady firing rate and the threshold sound pressure level of the single cell . The connection parameter set describes the connections between two groups of cells ; this is used to specify the connectivity within the neuronal circuitry of the model . We used the intrinsic values in the former work of Hancock and Voigt6 -LRB- Table 2 -RRB- , as they have little or no effect on the spectral response characteristics . FIGURE 6 . Best fit of Unit L041499 13.01 and parameter sensitivity analysis . -LRB- A -RRB- Best fit of the model to physiological data . The best frequency of physiological data and the center frequency of the model data have been lined up through simply translation of the physiological data rate vs. cutoff frequency plot . The circles represents the physiological data and the solid line represents the best fit of the model . The dashed line represents the average spontaneous rate of the physiological data . The minimum error Фmin and the values of the three important parameters of the best fit are shown in the upper left . -LRB- B -RRB- -- -LRB- D -RRB- Contours of equal Ф/Фmin -LRB- 1.5 , 2 , 5 , 10 , and 15 , from inside to outside -RRB- showing the sensitivity of the fit to the parameter values . , , and are held to be the value in the best fit corresponsively in plots -LRB- B -RRB- , -LRB- C -RRB- , and -LRB- D -RRB- . FIGURE 7 . Best fit of Unit L110199 2.01 and parameter analysis . Organized as in Fig. 6 except that in the contours plots -LRB- B -- D -RRB- equal Φ/Φmin of 1.1 was added to show the center.FIGURE 8 . Best fit of Unit L122299 4.01 and parameter analysis . Organized as in Fig. 6 . MacGregor11 pointed out that the steady state effect of one population on its target is proportional only to the product , which enabled us to set only one parameter from these three to be a free parameter ; here we chose , the strength from cell A to B. Thus , only three connection parameters , , , and were taken into account in the simulations . The values of the parameters are listed in Table 3 . As in Table 3 , three parameters were set to be main parameters that vary : , , and . Center offset and were set to two possible values 0 / − 0.2 octaves and 0/0 .2 octaves respectively according to the physiological data in cross-correlation studies . was set to 0.1 to obtain type III RM properties . Non-specific afferent firing rate was set to 80 , 90 , or 120 spikes/s according to the various spontaneous rates of type III units . Auditory Nerve Model AN-fibers were based on the model described by Carney .1 The input to the AN-fiber is a sound pressure signal and the output is a spike train . Each fiber consists of a gammatone filter that provides the frequency selectivity , a nonlinear process that generates an inner hair cell potential and a compartmental model of neurotransmitter release . This in turn creates an instantaneous firing probability . AN-fiber thresholds and spontaneous activity rates were randomly assigned from physiological distributions of these parameters .5 In order to save time , all stimuli were processed by the auditory nerve filter bank once and the spike time responses were saved for use in subsequent simulations of the DCN circuitry . Simulation Protocol and Data Analysis All simulations were carried out on an IBM pSeries 655 , which is a 48-processor system composed of six nodes . Each p655 node consists of eight Power4 processors running at 1.1 GHz and sharing 16 GB of memory . There are three levels of cache on this machine . Each processor has a 32 KB L1 cache and then each pair of processors share a 1.41 MB L2 cache , and each p655 node shares a 128 MB L3 cache . To create RMs , the model was stimulated by 50-ms tone bursts presented every 250 ms with a 32-ms delay while sound pressure levels varied from 0 to 60 dB SPL in 2 dB SPL steps and the frequency varied in 0.1 octave steps within a three octave band above and below 5 kHz . Thus , there are 31 frequency slices in one RM simulation . The spikes of the last 40 ms of each trial were used to compute the spontaneous rate and the spikes of the last 40 ms of each tone burst were used to compute the driven rate . In this paper , RMs composed of rate vs. frequency curves for seven levels were plotted -LRB- 0 -- 60 dB SPL in 10 dB SPL steps -RRB- . This is because the physiological data were taken only to 60 dB SPL due to increasing recording artifacts above these levels . At each level , the horizontal line represents the unit 's spontaneous discharge rate ; excitation and inhibition are shown as rates above and below this line , respectively . The excitatory region was filled black and the inhibitory region left blank . We used the same notch noises and broadband noise stimuli to obtain rate-level data as in Hancock and Voigt .6 In the simulation the sound level was varied from 0 to 90 dB SPL in 2 dB SPL steps and the noise bursts were presented for 200 ms in 1000-ms trials . To obtain rate vs. notch cutoff frequency plots , the firing rates from 46 to 54 dB SPL were averaged to capture the most inhibitory area , consistent with Hancock and Voigt .6 These plots are sometimes referred to as `` Mexican Hats . '' TABLE 4 . Data fits for all 15 type III units.In RefUnitBFSRΦminNSA1Fig . 3AL122299 3.014.3383.490003.25 a0 .250.902 Fig. 3BL021199 1.014.7773.81200.203.25 a0 .451.803 Fig. 3CL110199 2.012.6132.480003.000.120.254 Fig. 4AL111799 9.014.992.8800.2 − 0.22.000.120.255 Fig. 4BL122299 4.0111.0144.7800 − 0.21.250.240.906 Fig. 4CL110199 11.024.21211.080001.00 a0 .150.05 a7Fig . 4DL102599 7.0110.2265.190001.500.250.708 Fig. 4EL041499 16.013.8308.1800.2 − 0.22.500.120.209 Fig. 4FL021199 4.017.1152.880003.25 a0 .180.2510 Fig. 4GL102599 8.016.4225.680001.500.120.4011 Fig. 4HL041499 13.011.2352.480001.250.240.8012 Fig. 4IL041499 1.021.399.8800.203.25 a0 .060.1513 Fig. 4JL041499 1.011.2495.8800 − 0.21.00 a0 .240.9014 Fig. 4KL040700 3.0111.4201.3800.203.000.240.2515 Fig. 4LL040700 7.0115.1240.4120002.750.452.00 aNote . Ref . refers to Parsons et al. 15 ; BF in kHz , SR in spikes/s , NSA in spikes/s . aThe parameter value reached the parameter limit . RESULTS Fit of the Model to Physiological Data and Parameter Sensitivity Analysis Discharge rate versus cutoff frequency plots were used to evaluate the effectiveness of our DCN model quantitatively when fit to the physiological data obtained by Parsons et al. 15 The physiological discharge rate versus cutoff frequency plots were translated along the frequency axis to align the BF of the unit to that of the model at 5 kHz without shape change in logarithmic frequency plots . This operation was done mathematically by scaling the cutoff frequencies by CF/BF , where CF is the center frequency of model and BF is the best frequency of physiological data . This operation also made the after-translation physiological data look like that obtained by notch noise stimuli centered at 5 kHz while notch widths scaled by the factor CF/BF . Then linear interpolation was used to obtain `` new '' physiological data corresponding to model data cutoff frequencies and this made comparison of fits among different units possible . Here we used an objective function , Φ , to demonstrate quantitatively the fit of one unit , or for comparisons across the units as done by Hancock and Voigt6 : where yd and ym are the mean firing rates in the physiological data and model results respectively , and N is the total number of cutoff frequencies used in the comparison . The value of the objective function is scaled by 500 to make it more readable , consistent with Hancock and Voigt .6 The three parameters , , and , were systematically varied to get the `` best fit '' -LRB- minimum Φ , Φmin -RRB- and also to show contours of equal goodness or quality of fit . Center offsets and , were set to 0 or − 0.2 octaves and 0 or 0.2 octaves respectively , and the fitting results were only slightly different with similar Φmin when best fits were obtained . Non-specific afferent -LRB- NSA -RRB- firing rate was set to 80 , 90 , or 120 spikes/s to adjust to the various spontaneous rates of type III units and this improved the fits , especially for some low-spontaneous rate units . An alternative way to adjust spontaneous rate was to vary the connection strength of NSA to P-cells as in Hancock and Voigt .6 These two methods show no difference from each other since P-cells are the only targets of the NSA and the different firing rates of NSA or the strength of the connection between NSA and P-cells will only affect the total NSA excitatory inputs to P-cells . Figures 6 -- 8A show the best fits to three physiological units labeled Unit L041499 13.01 , Unit L110199 2.01 and Unit L122299 4.01 . Figures 6 -- 8B -- D are contours showing parameter sensitivity analyses of these parameters . As shown in the figures , varied from 1.0 to 3.25 octaves in steps of 0.25 octaves ; varied from 0.05 to 0.5 or from 0.03 to 0.3 ; and varied from 0.1 to 1.0 or from 0.05 to 0.5 Very good fits are obtained for these three units . Not only are the values of Φmin small , but the centers of the contours are shown clearly and the contours themselves resemble ellipses.FIGURE 9 . Response maps and rate-level plots for notch noise stimuli . -LRB- A , B -RRB- Response map composed of rate vs. frequency curves plot of Unit L041499 13.01 and the model using the set of best fit parameters in Fig. 6 . The best frequency of this unit is 1.2 kHz and has been moved to 5 kHz with the same method in `` Mexican Hat '' plots for easier comparison to model RM . -LRB- C , D -RRB- Rate-level curves of broad-band noise and notch noise responses of the unit and the model data.FIGURE 10 . Response maps and rate-level plots responding to notch noise stimuli . Unit L110199 2.01 . The best frequency of this unit is 2.6 kHz . Picture is organized as in Fig. 9 . Table 4 shows the results for all 15 units . Of the 15 units from Parsons et al. , 15 eight units had best fits within closed ellipses . Another six units also gave good fits , although the parameter in these cases was actually at its limit -LRB- four reaching the upper limit and one reaching the lower limit of the parameter range , one reaching the upper limit of -RRB- . These 14 units had values of Φmin < 9.8 , which shows that the model is robust . The final unit was difficult to fit because the rate vs. notch noise cutoff frequency plot showed additional features -LRB- see Parsons et al. , Fig. 4C -RRB- .15 Comparisons Between Physiological Data and Model Data : Response Maps and Notch Noise Rate-Level Curves Since the rate vs. notch noise cutoff frequency curves are the functions that were modeled , to what extent do the model RMs and notch noise rate-level curves resemble their physiological counterparts ? Figures 9 -- 11A and B show the physiological RMs together with the model RMs for the units of Figures 6 -- 8 . The model RMs were simulated using the model parameters set found when Φ = Φmin . For two sets of RMs we found that the model RMs show insufficient excitatory regions across the frequency axis when high-level tone bursts are given . This occurs in another five units where the physiological data show wide excitatory regions and indicates that our model parameters may need adjustment to get better tone responses for these units . All three model RMs , however , show nice center excitatory regions and sideband inhibitory regions , which is characteristic of type III units . Figures 9 -- 11C and D show the rate-level curves corresponding to the BBN and NN stimuli . For Unit L041499 13.01 and Unit L110199 2.01 -LRB- Figures 9 and 10 -RRB- , the maximum notch widths could not reach 8.0 kHz since the units had low BFs -LRB- 1.2 and 2.6 kHz -RRB- . Unit L122299 4.01 has a higher BF at 11.0 kHz and the notch width was varied from 0.5 to 8.5 kHz as in Parsons et al. 15 In general , the responses to BBN are greatest in both model and physiological data , but there are great differences in shape , maximum rate , saturation rates and thresholds . The model 's curves for the NN stimuli are qualitatively similar to the physiology , but again several differences are observed . The model seems to be providing more inhibition than that found in the physiology . DISCUSSION Spirou and Young18 showed that type IV units in the cat DCN -LRB- projection neurons -RRB- are notch detectors . Hancock and Voigt6 used a computational model of cat DCN circuitry , whose P-cells had type IV response properties to show that these also have notch noise sensitivity and that the variability seen in the units ' responses to notch noise are accounted for by parameter variations rather than different circuit arrangements . The gerbil DCN shows far fewer type IV units than the cat and the projection neurons in gerbil DCN are predominately type III units . Parsons et al. 15 showed that gerbil type III units are sensitive to notch noises in much the same way as cat DCN type IV units . Thus , it appears that DCN projection neurons from both species are doing the same analysis on notch noise and thus may play similar roles in sound source detection in the median plane . This study demonstrates that the same computational model used to model the notch noise behavior in cat type IV units can be used to model type III unit notch noise sensitivity in gerbil DCN . The fact that we can use the same neural circuit to model the DCN in both species suggests that perhaps this aspect of the DCN neural circuit is invariant across species . The Parsons et al. 15 data set of 15 type III units provided a rich database to test our model . In all , the model provided reasonable fits to 14 type III units . Prior to this , notch noise responses from only three type IV units from cat were available for modeling .6,18 The objective function used here is the same one used in Hancock and Voigt6 ; this allowed direct comparisons to that study . In this study , the RMs of the model P-cells with the best-fit parameters were compared to the physiological units ' RMs . The model RMs had excitatory center regions with inhibitory surrounds . These are similar to the physiological RMs and correspond to type III units . The model rate-level curves to broadband noise were similar to the physiological data , but the rate-level curves to notch noise showed many differences . Of course , these data are not part of the objective function used and so we would expect differences . It is possible to modify the objective function to include the rate-level curves , the RMs , or both , and this may be done in the future.FIGURE 11 . Response maps and rate-level plots responding to notch noise stimuli . Unit L122299 4.01 . The best frequency of this unit is 11.0 kHz . Picture is organized as in Fig. 9 . Our model also shows reduced sensitivity to in half of the units . Four units reached the upper limit of octave and two units reached the lower limit of octave , which indicates that better fits might be obtained by extending the range of . This is not possible in the present version of the model , which is limited to a four octave frequency range . Non-specific afferent -LRB- NSA -RRB- firing rate affected the Φmin in a significant way for some units . This was not true , however , for most of the units , where NSA firing rates of 80 , 90 , and 120 spikes/s produced similar Φmin . The best-fit parameter values and sensitivity plots , however , would change significantly when different rates applied . In general , the parameters of the best fits showed no specific relations to unit BF or spontaneous rates . The DCN computational models for cat and gerbil are identical in form and differ only in connectivity . The model is still incomplete , however , in that it lacks model cartwheel cells , granule cells and other cell types known to exist . The CN also receives input from somatosensory neurons that subserve tactile and kinesthetic sensations from the trigeminal ganglion16 the interpolar and caudal spinal trigeminal ganglion ,8,25 cuneate nucleus19 as well as from other locations . These indicate that in addition to auditory function , the CN is involved in sensory integration at a very early stage in the brain . These are not currently in the model , but provide rich directions to explore in the future .

O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O

hearing
sound localization
parameter estimation
sensitivity analysis


The Effects of Time Varying Curvature on Species Transport in Coronary Arteries Alterations in mass transport patterns of low-density lipoproteins -LRB- LDL -RRB- and oxygen are known to cause atherosclerosis in larger arteries . We hypothesise that the species transport processes in coronary arteries may be affected by their physiological motion , a factor which has not been considered widely in mass transfer studies . Hence , we numerically simulated the mass transport of LDL and oxygen in an idealized moving coronary artery model under both steady and pulsatile flow conditions . A physiological inlet velocity and a sinusoidal curvature waveform were specified as velocity and wall motion boundary conditions . The results predicted elevation of LDL flux , impaired oxygen flux and low wall shear stress -LRB- WSS -RRB- along the inner wall of curvature , a predilection site for atherosclerosis . The wall motion induced changes in the velocity and WSS patterns were only secondary to the pulsatile flow effects . The temporal variations in flow and WSS due to the flow pulsation and wall motion did not affect temporal changes in the species wall flux . However , the wall motion did alter the time-averaged oxygen and LDL flux in the order of 26 % and 12 % respectively . Taken together , these results suggest that the wall motion may play an important role in coronary arterial transport processes and emphasise the need for further investigation . Introduction The formation of an atherosclerotic plaque is dependent upon a range of factors including local variations in haemodynamics , the transport of chemical species and cells from the blood to the arterial wall and their accumulation within the arterial wall .3,21,38 The local susceptibility of coronary arteries to atherosclerosis has been putatively linked to the physiological motion arising out of their anatomical attachment to the epicardial surface of the beating heart .7,8 Coronary arteries exhibit highly complex wall motion patterns such as bending , stretching , twisting , vessel torsion and vessel displacement due to the cardiac and respiratory motion of the heart .7,22 Despite their seemingly important role , to date , only a handful of studies have looked at the effects of wall motion on atherogenesis . Early studies investigated the axial and lateral movement of a straight tube model of the coronary arteries and showed that the wall shear stress -LRB- WSS -RRB- patterns were significantly affected by the axial movement and not so by the lateral movement .4,23 Other studies examined the effects of time dependent curvature on the flow field and demonstrated that the flow patterns were affected by the wall motion due to instantaneous changes in curvature and that the WSS patterns were strongly dependent on the frequency of curvature variation .18,24,32,33 Another study considered the cyclic flexion of coronary arteries and reported that the plaque advancement was proportional to the flexion angle .37 More recently , studies have concluded that the dynamic curvature effects could play an important role in the flow and WSS patterns in curved and in bifurcating arteries .28,40,41 Nevertheless , these studies reported that the temporal WSS variations produced by the wall motion was only secondary to that caused by the pulsatility of blood flow .28,41 So far studies have only paid attention to the haemodynamic effects of wall motion . However , it has now come to light that the transport of atherogenic macromolecules such as LDL and solutes such as oxygen from the blood to the arterial wall play an important role in atherogenesis .38 Interestingly , with the exception of a single study , no other past study has investigated the effects of species transport to the moving walls of coronary arteries .17 Even in that study only the effects of diameter variation -LRB- < 6 % of mean diameter -RRB- on oxygen transport was considered . Hence , the aim of our work was to determine the effects of wall motion on mass transport patterns in coronary arteries with physiological flow . For our investigation , we chose a geometric model of the left anterior descending coronary artery -LRB- LAD -RRB- because of its high predilection to atherosclerosis among the major coronary arteries .12 Computer aided design -LRB- CAD -RRB- models of the LAD were employed in this study as opposed to realistic geometries to simplify the complexity of the problem definition , solution and analysis of results so as to gain a fundamental understanding of the haemodynamic effects of wall motion in coronary arteries . The species of interest in this study were low-density lipoprotein -LRB- LDL -RRB- and oxygen . Methods Computational Model A 3-D computational model of the LAD -LRB- Fig. 1 -RRB- was constructed based on the average LAD dimensions -LRB- diameter -LRB- d -RRB- = 0.36 cm and mean radius of curvature -LRB- Rmean -RRB- = 4 cm -RRB- ,9,15,16 using the CFD pre-processor Gambit -LRB- v2 .0.4 -- Fluent Inc. , Lebanon , NH -RRB- . The total length -LRB- l -RRB- of the LAD considered in this study was 4.68 cm which was equivalent to 13 tube diameters . In order to make the simulation as simple as possible , the model was assumed to be rigid with a circular cross section of constant radius and a constant uniplanar curvature over its entire length at any instant in time . An unstructured grid composed of hexahedral mesh elements with a thin boundary layer was generated using the Cooper meshing scheme2 with strict adherence to mesh quality metrics . Figure 1 . A schematic model of the left anterior descending coronary artery -LRB- LAD -RRB- . Blood Flow Boundary Conditions The flow boundary conditions controlling the model included no slip at the tube wall , velocity inlet at the inlet , symmetry condition at the plane of symmetry and a traction free outflow condition . The LAD inlet velocity profile was assumed to be parabolic . A pulsatile physiological LAD inlet velocity waveform -LRB- Fig. 2 -RRB- derived from a previous publication by Marcus et al. 22 was applied at the inlet boundary through a user defined function . The time-averaged mean velocity -LRB- ν -RRB- , mean flow rate -LRB- qo -RRB- Reynolds number -LRB- Re -RRB- and Dean number -LRB- κ -RRB- corresponding to this flow waveform were 9.715 cm/s , 0.99 ml/s , 105 and 22.3 . The Dean number is defined as below Figure 2 . The LAD inlet flow and curvature waveform . The Womersley parameter -LRB- α -RRB- of the flow was 2.78 and the period of a cardiac cycle -LRB- T -RRB- was 0.8 s. Blood was considered to be an incompressible , homogeneous , Newtonian fluid with a viscosity -LRB- μ -RRB- of 0.0035 Pa s and a density -LRB- ρ -RRB- of 1050 kg/m3 . Past arterial flow dynamic studies have shown that these assumptions are reasonable enough to use in haemodynamic studies .16,29 Wall Motion Boundary Conditions The wall motion boundary condition used in this study was based on the model presented by Santamarina et al. 32 The main components of coronary artery motion are solid body motion and deformation . The deformation parameter -LRB- ε -RRB- was calculated as below where R is the radius of curvature of the LAD . The present study used a deformation parameter of 0.6 which falls within the physiological range of LAD deformations .15 A sinusoidal function as given below was used to specify the time dependent variation of curvature of the LAD -LRB- Fig. 2 -RRB- where ω is the angular frequency of deformation . A user-defined function was written to apply this curvature waveform as a grid motion boundary condition . The tube inlet was fixed constant in the simulations at all times . This allowed the centre of curvature to vary sinusoidally in a cardiac cycle . Mass Transfer Boundary Conditions The mass transport of species in the arterial lumen is mathematically described using the three dimensional convection-diffusion equation where C is the concentration of species , u is the velocity vector and D is the diffusivity of species in the blood . In this study , the transport of LDL and oxygen from the blood to the endothelium -LRB- fluid phase -RRB- , and across the endothelium -LRB- transendothelial -RRB- were modelled . The fluid phase transport of both LDL and oxygen takes place via convection and diffusion . However , the transendothelial transport of LDL and oxygen are governed by different mechanisms due to their relative size and the endothelial function . The semi permeable endothelial membrane , acts as a barrier for the entry of macromolecules such as LDL into the arterial wall , but allows blood plasma and small solutes such as oxygen to pass through it .10 This is reflected in the permeability of LDL -LRB- K = 2 × 10 − 8cm/s -RRB- which is five orders of magnitude smaller than that of oxygen -LRB- K = 2 × 10 − 03cm/s -RRB- .36 This indicates the need for modelling LDL and oxygen transport differently through appropriate mass transport boundary conditions . LDL transport was modelled as a concentration polarisation effect , a phenomenon which has been demonstrated to occur in large and medium sized arteries .5,11,39 Concentration polarisation is the accumulation and elevation of endothelial surface concentration of LDL due to plasma filtration through the semi permeable endothelium . LDL mass transport boundary conditions applied in this study were uniform inlet concentration , Co -LRB- = 1.2 mg/ml of blood -RRB- ,36 a zero gradient in concentration at the outlet and mass conservation of LDL at the luminal surface of the endothelium39 described as where Vw is the plasma filtration velocity normal to the wall -LRB- = 4 × 10 − 6cm/s -RRB- ,36 Cw is the endothelial surface concentration of LDL , n is the unit vector normal to the vessel wall , D is the physiological LDL diffusivity -LRB- = 5 × 10 − 8cm2/s -RRB- 39 and K is the overall mass transfer coefficient of LDL which was considered equivalent to the physiological endothelial permeability of LDL . Equation -LRB- 5 -RRB- was rearranged as below to implement a wall reaction boundary condition in the CFD solver . In this study , oxygen transport was modelled by applying a passive transport law for oxygen flux at the wall .31 It is known from previous estimates that the diffusion velocity of oxygen is orders of magnitude greater than the filtration velocity .1 Hence , the flux of oxygen to the wall is essentially the diffusive flux . To satisfy flux continuity at the wall the diffusive flux of oxygen to the wall should be equal to the wall flux of oxygen which can be modelled with a constant permeability wall boundary condition .31 Oxygen mass transport boundary conditions employed in this study were uniform inlet oxygen mass fraction -LRB- = 0.005 -RRB- ,20 a zero gradient in outlet concentration and a constant permeability wall boundary condition31 which is given as where the oxygen diffusivity was 1 × 10 − 5cm2/s1 . Numerical Method The commercial CFD solver Fluent -LRB- v6 .1.22 -- Fluent Inc. , Lebanon , NH -RRB- was used to simulate species transport in the moving LAD model . Fluent uses a finite volume technique to solve the three-dimensional -LRB- 3-D -RRB- unsteady equations of momentum , mass and species conservation . QUICK scheme was employed to solve the flow and species transport equations . The time discretization and pressure -- velocity coupling were achieved through second order time implicit , and PISO schemes respectively . The dynamic mesh model feature of the Fluent was used to implement the deformation of the model according to Eq . -LRB- 3 -RRB- . The dynamic mesh model uses an arbitrary Lagrangian -- Eulerian approach to move the mesh nodes . At each time step , the volume mesh is updated automatically by the solver based on the new co-ordinates of the boundary zones that were calculated by the user-defined function . At the beginning of each simulation , the radius of curvature of the model was at its maximum -LRB- = 6.4 cm -RRB- . At t = 0.2 s and 0.6 s the radius of curvature of the model was equivalent to its mean value . Modelling Approach To separate the effects of blood flow pulsatility and wall motion on mass transport patterns in coronary arteries so as to determine their relative importance , we undertook the following four different numerical simulations : Case 1 : Steady flow inlet -LRB- qo = 0.99 ml/s -RRB- in the static geometry with the mean radius of curvature -LRB- R = 4 cm -RRB- . Case 2 : Pulsatile flow inlet in the static geometry with the mean radius of curvature -LRB- R = 4 cm -RRB- , to isolate the haemodynamics effects of flow pulsatility.Case 3 : Steady flow inlet -LRB- qo = 0.99 ml/s -RRB- with dynamic geometry , to isolate the haemodynamics effects of wall motion.Case 4 : Pulsatile flow inlet with dynamic geometry , to study the combined effects of pulsatility and wall motion . Model Validation The computational method was validated by simulating Santamarina et al. 's 32 coronary artery motion study in a curved model of exactly similar dimensions under identical flow and wall motion conditions . The flow and wall shear stress -LRB- WSS -RRB- patterns agreed well with that of Santamarina et al 's study and the maximum difference in the wall shear rate was less than 2 % . The mesh and time step independence of the solutions were investigated by successively refining the grid and the time step size by a factor of two . Computational solutions obtained with three different mesh densities -LRB- 80,310 , 156,861 and 310,464 nodes -RRB- and time step sizes -LRB- 0.025 s , 0.0125 s and 0.00625 s -RRB- were then compared against each other -LRB- Fig. 3 -RRB- . Mesh and time step independence were considered to have been achieved when the maximum difference between the time-averaged normalised outer and inner wall flux of LDL -LRB- Cw/Co -RRB- between the successive simulations was less than 2 % at any axial location . The outer and inner walls were defined as the wall regions farthest and shortest from the centre of curvature to represent epicardial and myocardial surface of the curved coronary arteries . This study was only carried out for LDL transport in the dynamic model with pulsatile flow because this was the computationally most challenging problem -LRB- Schmidt number , μ/ρD = 6.7 × 105 & Peclet number , dv/D = 7 × 107 -RRB- encountered in this study . Based on these criteria , a mesh containing 156,861 nodes with a time step size of 0.0125 s was chosen as the optimum case for simulation . Figure 3 . Comparison of the time-averaged normalised LDL wall flux predicted with three different meshes and time step sizes to demonstrate mesh -LRB- top panel -RRB- and time step -LRB- bottom panel -RRB- independence of the computational solutions . Results Numerical simulations of oxygen and LDL transport in a 3-D curved tube model of the LAD were carried out under physiological flow and mass transport conditions . The computational solutions concerned the prediction of flow , WSS and mass transport patterns . The results will be presented only for the first 10 diameters length as this length is more relevant to realistic coronary flow situation as discussed elsewhere .32 Velocity Patterns The temporal velocity magnitude -LRB- coloured contour plots -RRB- and secondary flow patterns -LRB- vector plots -RRB- in the static and dynamic models of the LAD are shown in Fig. 4 . The fundamental flow patterns in a curved tube flow situation with a characteristic velocity skewing towards the outer wall of bend and the formation of Dean vortices were clearly seen in all the models -LRB- Fig. 4 -RRB- . Although both the pulsatility -LRB- Fig. 4b -RRB- and wall motion -LRB- Fig. 4c -RRB- produced time dependent velocity skewing , the flow pulsation had a greater effect on velocity skewing than did wall motion . The secondary flow was also affected by both the pulsatility and wall motion -LRB- Fig. 4e and f -RRB- . However , the wall motion induced time dependent curvature -LRB- Fig. 4f -RRB- seemed to have a major influence on the secondary flow patterns . Adding the dynamic motion to the pulsatile flow reduced skewing over the entire cycle -LRB- Fig. 4d -RRB- . Moreover , this combination led to a reduction in secondary flow when the time dependent vessel curvature was small and vice versa during the peak curvature phases -LRB- Fig. 4g -RRB- . In general , both the velocity skewing and secondary flow increased with increase in curvature and the flow rate -LRB- Fig. 4 -RRB- . Figure 4 . Temporal axial and secondary flow patterns in a plane located at 9 diameters from the inlet in the -LRB- a -RRB- steady flow static model -LRB- R = 4 cm , qo = 0.99 ml/s -RRB- , -LRB- b and e -RRB- pulsatile flow static model -LRB- R = 4 cm -RRB- , -LRB- c and f -RRB- steady flow dynamic model -LRB- qo = 0.99 ml/s -RRB- , and -LRB- d and g -RRB- pulsatile flow dynamic model . The outer and inner walls are located at the top right corner and top left corner of individual figures respectively . WSS Patterns The effects of pulsatility and wall motion on the instantaneous WSS distribution in a cardiac cycle at two axial locations -LRB- 5D & 9D -RRB- are plotted in Fig. 5 . The effect of motion in steady flow at the outer wall was to increase the WSS when the curvature was maximum . Along the inner wall the opposite effects were seen with a drop in the WSS at the maximum curvature . The wall motion also produced asymmetric WSS distribution over time in the steady flow dynamic model despite the imposed symmetric wall motion curve . Figure 5 . Temporal axial inner and outer WSS distribution at 5 diameters -LRB- top panel -RRB- and 9 diameters -LRB- bottom panel -RRB- from the LAD inlet . The pulsatility of blood flow caused time varying WSS distribution along the outer and inner walls . The WSS increased with the pulsatile flow rate and overall this increase was greater along the outer wall . Hence peak WSS was noted along the outer wall at maximum flow rate . Together , wall motion and flow pulsatility produced instantaneous variations in the outer and inner WSS which were somewhat similar to the WSS distribution in the pulsatile flow static model . Although the gross distribution of WSS primarily followed the flow rate curve , the dynamic curvature effects were evident in the central section where the location and magnitude of peak WSS were different from that of the pulsatile flow static model . The time-averaged axial WSS distribution at the outer and inner walls in all the models is shown in Fig. 6 . In general , the spatial dependence of time averaged WSS was evident in all the models . In the steady flow static model the outer WSS increased to a peak and thereafter was found to decrease very gradually . In contrast , the inner WSS decreased to a minimum and was thereafter found to increase slowly . In comparison , the flow pulsatility caused a big increase in the outer WSS in the pulsatile flow static model . However , the wall motion only produced a little change in the outer WSS in the steady flow dynamic model . Figure 6.Time-averaged axial WSS distribution along the inner and outer walls of the LAD . The combined effects of wall motion and pulsatility caused an increase in the outer WSS up to its peak value and a subsequent decrease in the WSS in comparison with that of the pulsatile flow static model , and an overall increase in the outer WSS when compared to the steady flow models . Along the inner wall the flow pulsatility and the wall motion produced contrasting effects whereby the former decreased the WSS and the later increased the WSS . Due to this opposing effect of wall motion and flow pulsatility their combination had only caused smaller variations in the WSS distribution . Oxygen Transport The time dependent non-dimensional oxygen wall flux -LRB- KCw/KCo -RRB- distribution at 5D and 9D is presented in Fig. 7 . The oxygen wall flux did not vary appreciably over time in a cardiac cycle . At any instance the outer wall oxygen flux was higher than that at the inner wall . These patterns were similar in all the models studied . When compared to the oxygen wall flux in the steady flow static model , pulsatility caused a reduction and wall motion an increase in the oxygen wall flux along the inner wall . Along the outer wall the effects of pulsatility and wall motion were opposite to this . In total , the net effect of pulsatility and wall motion was an increase in oxygen flux at the inner wall and decrease at the outer wall when compared to the steady flow static model . This effect was more prominent in the distal region -LRB- Fig. 7c and d -RRB- . Figure 7 . Temporal normalised oxygen flux to the inner and outer walls at 5 diameters -LRB- top panel -RRB- and 9 diameters -LRB- bottom panel -RRB- from the LAD inlet . The time-averaged non-dimensional oxygen wall flux to the outer and inner walls of all the four models are shown in Fig. 8 . The overall effects of pulsatility and wall motion were similar to that predicted with the time dependent oxygen wall flux distribution . The effect of flow pulsation was to decrease the oxygen flux to the inner wall . In contrast , the wall motion contributed to an increase in oxygen flux in the distal region of the inner wall . Hence the combined effects of wall motion and flow pulsatility were to decrease oxygen flux in the proximal and medial regions of the inner wall and a very gradual increase distally . Figure 8.Time-averaged normalised oxygen wall flux distribution along the inner and outer walls of the LAD . Along the outer wall after an initial sharp drop the flow pulsation caused a gradual increase and then a gradual decrease in the oxygen flux . The wall motion produced a net reduction in the oxygen flux to the outer wall from the inlet . When both the wall motion and flow pulsatility were included the oxygen flux decreased sharply closer to the inlet followed by a gradual increase in the medial region and a gradual decrease thereafter . Nevertheless , in all the models the time-averaged outer wall oxygen flux was more than that to the inner wall . LDL Transport The time dependent non-dimensional LDL wall flux -LRB- KCw/KCo -RRB- distribution at 5D and 9D is plotted in Fig. 9 . Like oxygen wall flux the LDL wall flux did not exhibit notable variations over time . The flow pulsatility caused an increase in LDL flux to the inner wall and a decrease in LDL flux to the outer wall when compared to the LDL flux distribution in the steady flow static model . In contrast , the wall motion induced LDL wall flux was lower along the inner wall . However , along the outer wall it produced a differential distribution whereby the LDL wall flux was lower than that in the steady flow static model at 5D and vice versa at 9D -LRB- Fig. 9c and d -RRB- . At 9D the outer wall LDL flux was greater than that along the inner wall . Also the outer wall flux in the steady flow dynamic model was higher than that in the pulsatile flow static model . Hence the net effect of wall motion and pulsatility produced LDL wall flux distribution in the pulsatile flow dynamic model such that it was between the range of LDL wall fluxes caused by the pulsatility and wall motion separately . Figure 9 . Temporal normalised LDL flux to the inner and outer walls at 5 diameters -LRB- top panel -RRB- and 9 diameters -LRB- bottom panel -RRB- from the LAD inlet . The time-averaged non-dimensional LDL wall flux to the outer and inner walls of all the four models are shown in Fig. 10 . In all the models the time-averaged LDL wall flux distribution was seen to vary appreciably with respect to the axial distance from the inlet . The flow pulsatility caused an increase and decrease in the flux along the inner and outer walls respectively when compared to the steady flow static model . In contrary , the wall motion caused a reduction in the inner wall LDL flux when compared to that of the steady flow static model . Whereas the outer wall LDL flux was lower and higher than that in the steady flow static model in the medial and distal regions respectively . The combined effects of pulsatility and wall motion led to a substantial increase in the outer wall LDL flux in the distal regions . In the proximal and medial regions the outer wall LDL flux was lower than that to the inner wall . However , the inner wall LDL flux increased to a maximum and was then seen to drop off to a value lower than that along the outer wall in the distal part . Figure 10.Time-averaged normalised LDL flux distribution along the inner and outer walls of the LAD . Discussion To determine the effects of wall motion , a factor which has received less attention in the previous mass transport studies in coronary arteries , mass transport patterns were computed in an idealised moving model of the LAD . In this study , LDL and oxygen transport were simulated because of their important role in atherogenesis . The flow patterns reported in the current study -LRB- Fig. 4 -RRB- qualitatively agree with the findings of past experimental and computational studies in that they all demonstrate a general trend of velocity skewing towards the outer wall of curvature and the appearance of in-plane secondary flow .16,27,29,32 The centripetal acceleration of the fluid due to the arterial curvature was primarily responsible for these flow patterns -LRB- Fig. 4a -RRB- . This basic flow feature was further altered by the flow pulsatility -LRB- Fig. 4b and e -RRB- and the dynamic curvature -LRB- Fig. 4c and f -RRB- . Their relative role on secondary flow can be characterised using the non-dimensional Dean number -LRB- κ -RRB- . From the definition of the Dean number -LRB- Eq . 1 -RRB- , it is clear that higher the flow rate or smaller the radius of curvature the stronger the Dean vortices will be . This was evident in the figures as the magnitude of skewing and the strength of Dean vortices increased with increasing flow rate and curvature and were maximum at maximum flow rate in the pulsatile flow static model and at maximum curvature in the steady flow dynamic model . Therefore , the combination of flow pulsatility and wall motion were responsible for the appearance of time dependent flow patterns in the pulsatile flow dynamic model -LRB- Fig. 4d and g -RRB- . The flow patterns seen in this study are qualitatively similar to that of Santamarina et al 's work which considered the effects of dynamic curvature .32 The flow patterns at t = 0.2 s and 0.6 s in the steady flow dynamic model where the time dependent curvature was equal to the mean curvature appeared to be different from that in the steady flow static model -LRB- Fig. 4c and f -RRB- . This suggests that a quasi-static model can not predict the flow patterns accurately and that the dynamic wall motion could be as important as the instantaneous geometry . Similar findings were reported by Ethier et al. who compared the flow patterns in frozen geometries and in a moving model of the right coronary artery segment .41 The WSS patterns calculated in the pulsatile flow static models are similar to those reported by other studies on coronary artery flow .16,27,29 Since the WSS is directly proportional to the gradient of velocity the WSS patterns followed the trends of velocity distribution . The velocity skewing towards the outer wall subjected the outer and inner walls to higher and lower WSS respectively except up to a distance of 1D where the entrance effects were dominant -LRB- Fig. 6 -RRB- . The instantaneous and mean WSS distribution in the pulsatile flow dynamic model closely followed that of the pulsatile flow static model except in the central section where the curvature was high . These results suggest that the effects of wall motion on temporal WSS patterns are only secondary to the pulsatile flow effects and agree well with the findings of other dynamic model studies .28,41 Concentration polarisation of LDL has been demonstrated to occur in large arteries by several past experimental and numerical studies .5,11,39 The build-up of LDL concentration at the apical surface of the endothelium due to concentration polarisation could increase the influx of LDL into the arterial wall .5,35 This , in turn , could initiate a sequence of biological and biophysical events that are associated with atherogenesis .5,25 Therefore , LDL transport was modelled as a concentration polarisation effect and the results were presented as normalised LDL wall flux -LRB- KCw/KCo -RRB- . The present study also considered oxygen transport because hypoxia or low oxygen tension in the arterial wall has long been considered to be a factor responsible for atherogenesis .1,17 It induces endothelial cell apoptosis , causes breakdown of endothelial barrier and forms interendothelial gaps which may cause lipid infiltration into the arterial wall . Recent evidences suggest that hypoxia may also increase vascular permeability to macromolecules by up regulating vascular endothelial growth factor -LRB- VEGF -RRB- release by vascular cells .38 Hence oxygen transport was modelled with a passive transport law for oxygen flux at the wall .31 To describe oxygen transport patterns normalised oxygen wall flux -LRB- KCw/KCo -RRB- was calculated . In general , the mass transport patterns were similar in all the models . The establishment of low velocity , low WSS regions along the inner wall due to reasons discussed earlier led to an overall elevation of LDL flux and a reduction in oxygen flux to the inner wall . In contrast , high axial flow and high WSS along the outer wall due to axial velocity skewing caused a reduction in LDL flux and an increase in oxygen flux to the outer wall . The higher the WSS the higher the convection and the convective transport of oxygen . Hence the oxygen flux to the outer wall was more than that to the inner wall . The opposite effects were noted for LDL transport because the degree of concentration polarisation is reduced where high flow and high WSS prevail .39 As a result , the LDL accumulation was generally lower along the outer wall and higher along the inner wall . The diffusive transport of LDL and oxygen however , had minor influence on their transport . This was especially true because of their constant low diffusion coefficients . These results qualitatively agree with mass transport studies in curved arterial transport models .17,17,39 Despite the large time dependent variations in the blood flow and WSS patterns , the flow pulsation only had a minimal effect on oxygen -LRB- Fig. 7 -RRB- and LDL -LRB- Fig. 9 -RRB- transport to the arterial wall . Consequently , the time-averaged oxygen and LDL wall fluxes -LRB- Figs. 8 and 10 -RRB- in the static pulsatile flow model were only slightly different from their instantaneous values . This less significant pulsatile effects on time-averaged mass transport patterns have also been reported by other studies .11,19,34 Schneiderman et al. 34 simulated oxygen transfer in straight arterial segments and concluded that the pulsatile flow has negligible effect on oxygen transport to the wall . Ma et al. 19 studied the effect of pulsatility on mass transfer in a sudden expansion region and reported that the time-averaged pulsatile mass transfer patterns were similar to the steady flow mass transfer patterns . Fatouraee et al. 11 simulated LDL concentration polarisation in a straight tube model of a carotid artery under pulsatile flow conditions . The variation of LDL wall concentration in their study was less than 0.014 % . Time dependent variations in mass transport patterns were also not present in the dynamic models -LRB- Figs. 7 and 9 -RRB- . As a result , the magnitude of time-averaged oxygen and LDL wall fluxes were only slightly different from their corresponding time dependent flux values -LRB- Figs. 8 and 10 -RRB- . Although the wall motion did not significantly affect the temporal mass transfer patterns in the dynamic models , an overall change in the mass transfer patterns was evident -LRB- Figs. 7 and 9 -RRB- . This change was also reflected in the time-averaged mass transfer patterns in the dynamic models -LRB- Figs. 8 and 10 -RRB- . The percentage change in time-averaged oxygen wall flux and LDL wall flux along the inner and outer walls between the static and dynamic models under both pulsatile and steady flow conditions was in the order of 26 % and 12 % respectively . The fact that the pulsatility did not alter mass transfer patterns and that there were variations in mass transport patterns in the dynamic models suggest that the wall motion may indeed be important in determining the time-averaged mass transport patterns in the dynamic models . These variations were not significant in the proximal regions of the LAD model . However , considerable changes in the mass transport patterns were present in the medial and distal regions of the LAD -LRB- Figs. 8 and 10 -RRB- . This regional difference occurred due to the specification of fixed inlet boundary condition . Although this is an effect of the fixed inlet boundary condition these results once again indicate that strong motion of the coronary arteries could have a significant impact on the luminal mass transport patterns . The species transport however , does not entirely depend on the WSS distribution . Karino et al drew a scatter plot to establish a relationship between the WSS and LDL surface concentration .39 The scatter in their data suggests that the species transport is also influenced by global flow patterns , path line of fluid flow , spatial location and residence time of LDL at a site . This was also found in the present study which noticed different values of LDL and oxygen wall flux for the same WSS value along the outer and inner walls -LRB- Figs. 6 , 8 and 10 -RRB- . The inner walls of a curved artery where low WSS prevails is known to be a predilection site for atherosclerosis .14,21 LDL and oxygen transport studies have reported impaired oxygen flux and enhanced LDL accumulation in the low shear regions .5,17,17,39 Similar to these findings the present study also predicted co-localisation of low WSS regions , high LDL flux regions and hypoxic regions along the inner wall of curvature -LRB- Figs. 6 and 10 -RRB- in both the static and dynamic models . Hence our results suggest that even when the wall motion is considered the susceptibility of inner wall regions of the curved arteries to atherosclerosis is maintained . Back et al. 1 investigated the relationship between oxygen transport and lipoprotein accumulation in various arterial configurations . They noticed a correspondence between oxygen deprivation -LRB- hypoxia -RRB- and enhanced lipoprotein accumulation in decelerated flow regions and separation locations and concluded that the oxygen transport behaves in an opposite way to lipoprotein accumulation . This contrasting behaviour of oxygen and LDL transport was also noticed in our present study . It would be interesting to study the indirect effects of hypoxia on LDL transport through pathways involving enhancement of endothelial permeability for macromolecules .38 The present study assumed a constant endothelial permeability of oxygen and LDL which , in reality , is spatially variable and is dependent on the concentration of species at the arterial wall surface and the WSS .26,36 Pulsatile variations in the WSS has been shown to alter the morphology of endothelial cells , increase macromolecular permeability to LDL , promote transendothelial diffusion and intimal entrapment of atherogenic particles and enhance oxygen wall flux .13,21 Therefore , it is possible that the temporal variations in WSS caused both by the pulsatility of blood flow and wall motion may have an effect on LDL and oxygen transport if shear dependent variations in physiological parameters including the endothelial permeability are considered in mass transport studies . So far their effects on species transport have only been considered in static models30 ,36 and it remains to be investigated how these factors would influence mass transport in moving arteries . This will be the subject of our future study . The limitations of this study centre on the assumptions that were made in the construction of computational models which include idealized geometry , wall motion , and biophysical factors governing LDL and oxygen transport . An idealised , uniformly curved model without complex geometric features such as tortuosity , twisting and non-uniform cross sections was used . Inclusion of these geometric features would be expected to create local regions of flow acceleration , deceleration , separation , recirculation and secondary motion1 ,27 which may exacerbate or reduce the effects of wall motion on haemodynamics and species transport . While we accept that the neglect of complex geometry may introduce some simplifications in the results it should be recognised that the geometric simplification enabled us to isolate the effects of geometry and to concentrate on the role of wall motion on species transport . The results predicted in this study are specific to the flow and wall motion used in our model and therefore would be expected to change if a different flow or wall motion were imposed . For instance , if the flow and wall motion were more out of phase it may be expected that this would reduce the instantaneous Dean number and therefore one might expect the impact of the dynamic wall motion to be decreased from what was predicted . Although the wall motion effects were secondary to the pulsatile flow effects this does not mean that the wall motion is unimportant . Given the considerable variability in coronary artery motion , geometry and flow patterns in vivo there may be instances where wall motion could play a more significant role and therefore its effects should be considered on a case by case basis . Our study assumed that the non-Newtonian effects caused by the presence of red blood cells are negligible and therefore used constant diffusivities for LDL and oxygen . However , the effective diffusivity of a species changes due to the shear induced microscopic motion of red blood cells .6 In regions of high WSS , the effective diffusivities of LDL and oxygen would increase . This would in turn enhance oxygen wall flux but reduce LDL wall flux due to LDL concentration polarisation . Therefore , the magnitude of LDL and oxygen wall fluxes predicted in our study would be expected to change but the general patterns of flux will remain the same . Conclusion In this study , we have incorporated the effects of idealised motion into a simulation of oxygen and LDL transport in the LAD . The LDL and oxygen transport were modelled differently due to their differences in in-vivo transport mechanisms . We carried out simulations under both the steady and pulsatile flow conditions in the static and dynamic models to quantify the relative importance of pulsatility and wall motion on mass transport patterns . This study showed that the effects of wall motion on mass transport in the coronary arteries could be modelled in the commercial CFD code Fluent . Our results predicted elevation of LDL flux , impaired oxygen flux and low wall shear stress -LRB- WSS -RRB- along the inner wall of curvature , a region known for its predilection to atherosclerosis . The temporal variations in velocity and WSS patterns in the dynamic model did not influence time dependent mass transport and were only secondary to the pulsatile flow effects . However , wall motion may influence mass transport if shear dependent variations in species permeability and other factors are considered in the calculations . Nevertheless , the wall motion did alter the time-averaged mass transfer in the medial and distal regions in the order of 26 % and 12 % for oxygen and LDL transport when compared to the corresponding static models . Taken together , these results suggest that wall motion may play an important role in coronary arterial transport processes . However , future studies on more physiologic models are warranted to gain a fundamental understanding of the role of wall motion on mass transport and atherosclerosis .

O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O B I I O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O B I I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O B I O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O B I O O O O O O O O O O O O O O O O O O O O O B I I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O B I O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O B I I O O O O O O O O O O O O O B I O O O O O B I I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O B I O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O B I O O O O O B I O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O B I O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O B I O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O B I O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O B I O O O O O O O O O O O O O O B I O O O O O O O O B I I O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O B I O O O O O O O O O B I O B I O O O O O O O O O O B O O O O O O O O O O O O O O O O B I I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O B I O O B I O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O B I O O O O O O O O O O O O O O O O O O O O O O O O O O O O O O B I O B I O O O

wall motion
coronary artery motion
mass transport
cfd
wall shear stress
ldl transport
oxygen transport