SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Very flexible computer-assisted numerical analysis is used to solve dynamic fluid-flow equations characterizing computer-controlled heat dissipation system developed for Space lab. Losses caused by bends, ties, fittings, valves, and like are easily included, and analysis can solve both steady-state and transient cases. It can also interact with parallel thermal analysis.
OBJECTIVE: The goal of this study was to analyze left ventricular outflow tract systolic acceleration (LVOT(Acc)) during alterations in left ventricular (LV) contractility and LV filling. BACKGROUND: Most indexes described to quantify LV systolic function, such as LV ejection fraction and cardiac output, are dependent on loading conditions. METHODS: In 18 sheep (4 normal, 6 with aortic regurgitation, and 8 with old myocardial infarction), blood flow velocities through the LVOT were recorded using conventional pulsed Doppler. The LVOT(Acc) was calculated as the aortic peak velocity divided by the time to peak flow; LVOT(Acc) was compared with LV maximal elastance (E(m)) acquired by conductance catheter under different loading conditions, including volume and pressure overload during an acute coronary occlusion (n = 10). In addition, a clinically validated lumped-parameter numerical model of the cardiovascular system was used to support our findings. RESULTS: Left ventricular E(m) and LVOT(Acc) decreased during ischemia (1.67 +/- 0.67 mm Hg.ml(-1) before vs. 0.93 +/- 0.41 mm Hg.ml(-1) during acute coronary occlusion [p < 0.05] and 7.9 +/- 3.1 m.s(-2) before vs. 4.4 +/- 1.0 m.s(-2) during coronary occlusion [p < 0.05], respectively). Left ventricular outflow tract systolic acceleration showed a strong linear correlation with LV E(m) (y = 3.84x + 1.87, r = 0.85, p < 0.001). Similar findings were obtained with the numerical modeling, which demonstrated a strong correlation between predicted and actual LV E(m) (predicted = 0.98 [actual] -0.01, r = 0.86). By analysis of variance, there was no statistically significant difference in LVOT(Acc) under different loading conditions. CONCLUSIONS: For a variety of hemodynamic conditions, LVOT(Acc) was linearly related to the LV contractility index LV E(m) and was independent of loading conditions. These findings were consistent with numerical modeling. Thus, this Doppler index may serve as a good noninvasive index of LV contractility.
Tremendous progress has been made at the level of sequential computation in phylogenetics. However, little attention has been paid to parallel computation. Parallel computing is particularly suited to phylogenetics because of the many ways large computational problems can be broken into parts that can be analyzed concurrently. In this paper, we investigate the scaling factors and efficiency of random addition and tree refinement strategies using the direct optimization software, POY, on a small (10 slave processors) and a large (256 slave processors) cluster of networked PCs running LINUX. These algorithms were tested on several data sets composed of DNA and morphology ranging from 40 to 500 taxa. Various algorithms in POY show fundamentally different properties within and between clusters. All algorithms are efficient on the small cluster for the 40-taxon data set. On the large cluster, multibuilding exhibits excellent parallel efficiency, whereas parallel building is inefficient. These results are independent of data set size. Branch swapping in parallel shows excellent speed-up for 16 slave processors on the large cluster. However, there is no appreciable speed-up for branch swapping with the further addition of slave processors (>16). This result is independent of data set size. Ratcheting in parallel is efficient with the addition of up to 32 processors in the large cluster. This result is independent of data set size. c2001 The Willi Hennig Society.
We have examined the feasibility of using massively-parallel and vector-processing supercomputers to solve large-scale optimization problems for human movement. Specifically, we compared the computational expense of determining the optimal controls for the single support phase of gait using a conventional serial machine (SGI Iris 4D25), a MIMD parallel machine (Intel iPSC/860), and a parallel-vector-processing machine (Cray Y-MP 8/864). With the human body modeled as a 14 degree-of-freedom linkage actuated by 46 musculotendinous units, computation of the optimal controls for gait could take up to 3 months of CPU time on the Iris. Both the Cray and the Intel are able to reduce this time to practical levels. The optimal solution for gait can be found with about 77 hours of CPU on the Cray and with about 88 hours of CPU on the Intel. Although the overall speeds of the Cray and the Intel were found to be similar, the unique capabilities of each machine are better suited to different portions of the computational algorithm used. The Intel was best suited to computing the derivatives of the performance criterion and the constraints whereas the Cray was best suited to parameter optimization of the controls. These results suggest that the ideal computer architecture for solving very large-scale optimal control problems is a hybrid system in which a vector-processing machine is integrated into the communication network of a MIMD parallel machine.
Transient mixing driven by buoyancy occurs through the birth of a symmetric Rayleigh-Taylor morphology (RTM) structure for large length scales. Beyond its critical bifurcation the RTM structure exhibits self-similarity and occurs on smaller and smaller length scales. The dynamics of the RTM structure, its nonlinear growth and internal collision, show that its genesis occurs from an explosive bifurcation which leads to the overlap of resonance regions in phase space. This event shows the coexistence of regular and chaotic regions in phase space which is corroborated with the existence of horseshoe maps. A measure of local chaos given by the topological entropy indicates that as the system evolves there is growth of uncertainty. Breakdown of the dissipative RTM structure occurs during the transition from explosive to catastrophic bifurcation; this event gives rise to annihilation of the separatrices which drives overlap of resonance regions. The global bifurcation of explosive and catastrophic events in phase space for the large length scale of the RTM structure serves as a template for which mixing occurs on smaller and smaller length scales. Copyright 2004 American Institute of Physics.
The use of spectral techniques to quantify short term heart rate fluctuations on the order of seconds to minutes has helped define the autonomic contributions to beat-to-beat control of heart rate. We used similar techniques to quantify the entire spectrum (0.00003-1.0 Hz) of heart rate variability during 24 hour ambulatory ECG monitoring. The ECG from standard Holter monitor recordings from normal subjects was sampled with the use of a phase locked loop, and a heart rate time series was constructed at 3 Hz. Frequency analysis of the heart rate signal was performed after a nonlinear filtering algorithm was used to eliminate artifacts. A power spectrum of the entire 24 hour record revealed power that was inversely proportional to frequency, 1/f, over 4 decades from 0.00003 to 0.1 Hz (period approximately 10 hours to 10 seconds). Displaying consecutive spectra calculated at 5 minute intervals revealed marked variability in the peaks at all frequencies throughout the 24 hours, probably accounting for the lack of distinct peaks in the spectra of the entire records.
This study aims to the validation of the left ventricular (LV) volume estimates obtained by processing volumetric data utilizing a segmentation model based on level set technique. The validation has been performed by comparing real-time volumetric echo data (RT3DE) and magnetic resonance (MRI) data. A validation protocol has been defined. The validation protocol was applied to twenty-four estimates (range 61-467 ml) obtained from normal and pathologic subjects, which underwent both RT3DE and MRI. A statistical analysis was performed on each estimate and on clinical parameters as stroke volume (SV) and ejection fraction (EF). Assuming MRI estimates (x) as a reference, an excellent correlation was found with volume measured by utilizing the segmentation procedure (y) (y=0.89x + 13.78, r=0.98). The mean error on SV was 8 ml and the mean error on EF was 2%. This study demonstrated that the segmentation technique is reliably applicable on human hearts in clinical practice.
The present report describes a desktop computer-based method for the quantitative assessment of the area occupied by immunoreactive terminals in close apposition to nerve cells in relation to the perimeter of the cell soma. This method is based on Fast Fourier Transform (FFT) routines incorporated in NIH-Image public domain software. Pyramidal cells of layer V of the somatosensory cortex outlined by GABA immunolabeled terminals were chosen for our analysis. A Leitz Diaplan light microscope was employed for the visualization of the sections. A Sierra Scientific Model 4030 CCD camera was used to capture the images into a Macintosh Centris 650 computer. After preprocessing, filtering was performed on the power spectrum in the frequency domain produced by the FFT operation. An inverse FFT with filter procedure was employed to restore the images to the spatial domain. Pasting of the original image to the transformed one using a Boolean logic operation called 'AND'ing produced an image with the terminals enhanced. This procedure allowed the creation of a binary image using a well-defined threshold of 128. Thus, the terminal area appears in black against a white background. This methodology provides an objective means of measurement of area by counting the total number of pixels occupied by immunoreactive terminals in light microscopic sections in which the difficulties of labeling intensity, size, shape and numerical density of terminals are avoided.
The human anterior cingulate cortex is distinguished by the presence of an unusual cell type, a large spindle neuron in layer Vb. This cell has been noted numerous times in the historical literature but has not been studied with modern neuroanatomic techniques. For instance, details regarding the neuronal class to which these cells belong and regarding their precise distribution along both ventrodorsal and anteroposterior axes of the cingulate gyrus are still lacking. In the present study, morphological features and the anatomic distribution of this cell type were studied using computer-assisted mapping and immunocytochemical techniques. Spindle neurons are restricted to the subfields of the anterior cingulate cortex (Brodmann's area 24), exhibiting a greater density in anterior portions of this area than in posterior portions, and tapering off in the transition zone between anterior and posterior cingulate cortex. Furthermore, a majority of the spindle cells at any level is located in subarea 24b on the gyral surface. Immunocytochemical analysis revealed that the neurofilament protein triple was present in a large percentage of these neurons and that they did not contain calcium-binding proteins. Injections of the carbocyanine dye DiI into the cingulum bundle revealed that these cells are projection neurons. Finally, spindle cells were consistently affected in Alzheimer's disease cases, with an overall loss of about 60%. Taken together, these observations indicate that the spindle cells of the human cingulate cortex represent a morphological subpopulation of pyramidal neurons whose restricted distribution may be associated with functionally distinct areas.
We present a numerical method which provides the ability to analyze digitized microscope images of retinal explants and quantify neurite outgrowth. Few parameters are required as input and limited user interaction is necessary to process an entire experiment of images. This eliminates fatigue related errors and user-related bias common to manual analysis. The method does not rely on stained images and handles images of variable quality. The algorithm is used to determine time and dose dependent, in vitro, neurotoxic effects of 1 GeV per nucleon iron particles in retinal explants. No neurotoxic effects are detected until 72 h after exposure; at 72 h, significant reductions of neurite outgrowth occurred at doses higher than 10 cGy.
OBJECTIVE: A new computerized semi-automatic method for left ventricular (LV) chamber segmentation is presented. METHODS: The LV is imaged by real-time three-dimensional echocardiography (RT3DE). The surface detection model, based on level set techniques, is applied to RT3DE data for image analysis. The modified level set partial differential equation we use is solved by applying numerical methods for conservation laws. The initial conditions are manually established on some slices of the entire volume. The solution obtained for each slice is a contour line corresponding with the boundary between LV cavity and LV endocardium. RESULTS: The mathematical model has been applied to sequences of frames of human hearts (volume range: 34-109 ml) imaged by 2D and reconstructed off-line and RT3DE data. Volume estimation obtained by this new semi-automatic method shows an excellent correlation with those obtained by manual tracing (r = 0.992). Dynamic change of LV volume during the cardiac cycle is also obtained. CONCLUSION: The volume estimation method is accurate; edge based segmentation, image completion and volume reconstruction can be accomplished. The visualization technique also allows to navigate into the reconstructed volume and to display any section of the volume.
Development of a temporal query framework for time-oriented medical databases has hitherto been a challenging problem. We describe a novel method for the detection of hemodynamic events in multiparameter trends utilizing wavelet coefficients in a MySQL relational database. Storage of the wavelet coefficients allowed for a compact representation of the trends, and provided robust descriptors for the dynamics of the parameter time series. A data model was developed to allow for simplified queries along several dimensions and time scales. Of particular importance, the data model and wavelet framework allowed for queries to be processed with minimal table-join operations. A web-based search engine was developed to allow for user-defined queries. Typical queries required between 0.01 and 0.02 seconds, with at least two orders of magnitude improvement in speed over conventional queries. This powerful and innovative structure will facilitate research on large-scale time-oriented medical databases.
The chinchilla crista ampullaris was studied in 10 samples, each containing 32 consecutive ultrathin sections of the entire neuroepithelium. Dissector methods were used to estimate the incidence of various synaptic features, and results were confirmed in completely reconstructed hair cells. There are large regional variations in cellular and synaptic architecture. Type I and type II hair cells are shorter, broader, and less densely packed in the central zone than in the intermediate and peripheral zones. Complex calyx endings are most common centrally. On average, there are 15-20 ribbon synapses and 25-30 calyceal invaginations in each type I hair cell. Synapses and invaginations are most numerous centrally. Central type II hair cells receive considerably fewer afferent boutons than do peripheral type II hair cells, but have similar numbers of ribbon synapses. The numbers are similar because central type II hair cells make more synapses with the outer faces of calyx endings and with individual afferent boutons. Most afferent boutons get one ribbon synapse. Boutons without ribbon synapses were only found peripherally, and boutons getting multiple synapses were most frequent centrally. Throughout the neuroepithelium, there is an average of three to four efferent boutons on each type II hair cell and calyx ending. Reciprocal synapses are rare. Most synaptic ribbons in type I hair cells are spherules; those in type II hair cells can be spherical or elongated and are particularly heterogeneous centrally. Consistent with the proposal that the crista is concentrically organized, the intermediate and peripheral zones are each similar in their cellular and synaptic architecture near the base and near the planum. An especially differentiated subzone may exist in the middle of the central zone.
Variations in respiratory pattern influence the heart rate spectrum. It has been suggested, hence, that metronomic respiration should be used to correctly assess vagal modulation of heart rate by using spectral analysis. On the other hand, breathing to a metronome has been reported to increase heart rate spectral power in the high- or respiratory frequency region; this finding has led to the suggestion that metronomic respiration enhances vagal tone or alters vagal modulation of heart rate. To investigate whether metronomic breathing complicates the interpretation of heart rate spectra by altering vagal modulation, we recorded the electrocardiogram and respiration from eight volunteers during three breathing trials of 10 min each: 1) spontaneous breathing (mean rate of 14.4 breaths/min); 2) breathing to a metronome at the rate of 15, 18, and 21 breaths/min for 2, 6, and 2 min, respectively; and 3) breathing to a metronome at the rate of 18 breaths/min for 10 min. Data were also collected from eight volunteers who breathed spontaneously for 20 min and breathed metronomically at each subject's mean spontaneous breathing frequency for 20 min. Results from the three 10-min breathing trials showed that heart rate power in the respiratory frequency region was smaller during metronomic breathing than during spontaneous breathing. This decrease could be explained fully by the higher breathing frequencies used during trials 2 and 3 of metronomic breathing. When the subjects breathed metronomically at each subject's mean breathing frequency, the heart rate powers during metronomic breathing were similar to those during spontaneous breathing. Our results suggest that vagal modulation of heart rate is not altered and vagal tone is not enhanced during metronomic breathing.
The widespread availability of methods for noninvasive assessment of brain structure has enabled researchers to investigate neuroimaging correlates of normal aging, cerebrovascular disease, and other processes; we designate such studies as image-based clinical trials (IBCTs). We propose an architecture for a brain-image database, which integrates image processing and statistical operators, and thus supports the implementation and analysis of IBCTs. The implementation of this architecture is described and results from the analysis of image and clinical data from two IBCTs are presented. We expect that systems such as this will play a central role in the management and analysis of complex research data sets.