Development of one-dimensional and two-dimensional computational tools that simulate steady internal condensing flows in terrestrial and zero-gravity environments

This dissertation presents an effective quasi one-dimensional (1-D) computational simulation tool and a full two-dimensional (2-D) computational simulation methodology for steady annular/stratified internal condensing flows of pure vapor. These simulation tools are used to investigate internal condensing flows in both gravity as well as shear driven environments. Through accurate numerical simulations of the full two dimensional governing equations, results for laminar/laminar condensing flows inside mm-scale ducts are presented. The methodology has been developed using MATLAB/COMSOL platform and is currently capable of simulating film-wise condensation for steady (and unsteady flows). Moreover, a novel 1-D solution technique, capable of simulating condensing flows inside rectangular and circular ducts with different thermal boundary conditions is also presented. The results obtained from the 2-D scientific tool and 1-D engineering tool, are validated and synthesized with experimental results for gravity dominated flows inside vertical tube and inclined channel; and, also, for shear/pressure driven flows inside horizontal channels. Furthermore, these simulation tools are employed to demonstrate key differences of physics between gravity dominated and shear/pressure driven flows. A transition map that distinguishes shear driven, gravity driven, and “mixed” driven flow zones within the non-dimensional parameter space that govern these duct flows is presented along with the film thickness and heat transfer correlations that are valid in these zones. It has also been shown that internal condensing flows in a micro-meter scale duct experiences shear driven flow, even in different gravitational environments. The full 2-D steady computational tool has been employed to investigate the length of annularity. The result for a shear driven flow in a horizontal channel shows that in absence of any noise or pressure fluctuation at the inlet, the onset of non-annularity is partly due to insufficient shear at the liquid-vapor interface. This result is being further corroborated/investigated by R. R. Naik with the help of the unsteady simulation tool. The condensing flow results and flow physics understanding developed through these simulation tools will be instrumental in reliable design of modern micro-scale and spacebased thermal systems.

Analytical investigation of squeeze film dampers

Squeeze film damping effects naturally occur if structures are subjected to loading situations such that a very thin film of fluid is trapped within structural joints, interfaces, etc. An accurate estimate of squeeze film effects is important to predict the performance of dynamic structures. Starting from linear Reynolds equation which governs the fluid behavior coupled with structure domain which is modeled by Kirchhoff plate equation, the effects of nondimensional parameters on the damped natural frequencies are presented using boundary characteristic orthogonal functions. For this purpose, the nondimensional coupled partial differential equations are obtained using Rayleigh-Ritz method and the weak formulation, are solved using polynomial and sinusoidal boundary characteristic orthogonal functions for structure and fluid domain respectively. In order to implement present approach to the complex geometries, a two dimensional isoparametric coupled finite element is developed based on Reissner-Mindlin plate theory and linearized Reynolds equation. The coupling between fluid and structure is handled by considering the pressure forces and structural surface velocities on the boundaries. The effects of the driving parameters on the frequency response functions are investigated. As the next logical step, an analytical method for solution of squeeze film damping based upon Green’s function to the nonlinear Reynolds equation considering elastic plate is studied. This allows calculating modal damping and stiffness force rapidly for various boundary conditions. The nonlinear Reynolds equation is divided into multiple linear non-homogeneous Helmholtz equations, which then can be solvable using the presented approach. Approximate mode shapes of a rectangular elastic plate are used, enabling calculation of damping ratio and frequency shift as well as complex resistant pressure. Moreover, the theoretical results are correlated and compared with experimental results both in the literature and in-house experimental procedures including comparison against viscoelastic dampers.

Post processing of multiple GPS receivers to enhance baseline accuracy

Inexpensive, commercial available off-the-shelf (COTS) Global Positioning Receivers (GPS) have typical accuracy of ±3 meters when augmented by the Wide Areas Augmentation System (WAAS). There exist applications that require position measurements between two moving targets. The focus of this work is to explore the viability of using clusters of COTS GPS receivers for relative position measurements to improve their accuracy. An experimental study was performed using two clusters, each with five GPS receivers, with a fixed distance of 4.5 m between the clusters. Although the relative position was fixed, the entire system of ten GPS receivers was on a mobile platform. Data was recorded while moving the system over a rectangular track with a perimeter distance of 7564 m. The data was post processed and yielded approximately 1 meter accuracy for the relative position vector between the two clusters.

Length of the annular regime for condensing flows inside a horizontal channel - the experimental determination of its values and its trends

The experiments observe and measure the length of the annular regime in fully condensing quasi-steady (steady-in-the-mean) flows of pure FC-72 vapor in a horizontal condenser (rectangular cross-section of 2 mm height, 15 mm width, and 1 m length). The sides and top of the duct are made of clear plastic that allows flow visualization. The experimental system in which this condenser is used is able to control and achieve different quasi-steady mass flow rates, inlet pressures, and wall cooling conditions (by adjustment of the temperature and flow rate of the cooling water flowing underneath the condensing-plate). The reported correlations and measurements for the annular length are also vital information for determining the length of the annular regime and proposing extended correlation (covering many vapors and a larger parameter set than the experimentally reported version here) by ongoing independent modeling and computational simulation approach.

DEVELOPMENT OF CONFORMABLE CNG TANK FOR AUTOMOTIVE APPLICATION

The purpose of this study is to design, develop and integrate a Compressed Natural Gas (CNG) tank that will have a conformable shape for efficient storage in a light-duty pick-up truck. The CNG tank will be a simple rectangular box geometry to demonstrate capability of non-cylindrical shapes. Using CAD drawings of the truck, a conformable tank will be designed to fit under the pick-up bed. The intent of the non-cylindrical CNG tank is to demonstrate improvement in size over the current solution, which is a large cylinder in the box of a pick-up truck. The geometry of the tank’s features is critical to its size and strength. The optimized tank design will be simulated with Finite Element Analysis (FEA) to determine critical stress regions, and appropriate design changes will be made to reduce stress concentration. Following the American National Standard Institute (ANSI) guide, different aluminum alloys will be optimized to obtain the best possible result for the CNG tank.

Differential optical densities of intraretinal spaces

PURPOSE: To test the hypothesis that hyporeflective spaces in the neuroretina found on optical coherence tomography (OCT) examination have different optical reflectivities according to whether they are associated with exudation or degeneration. METHODS: Retrospective analysis of eyes with idiopathic perifoveal telangiectasia (IPT), diabetic macular edema (DME), idiopathic central serous chorioretinopathy (CSC), retinitis pigmentosa (RP), or cone dystrophy (CD) and eyes of healthy control subjects. OCT scans were performed. Raw scan data were exported and used to calculate light reflectivity profiles. Reflectivity data were acquired by projecting three rectangular boxes, each 50 pixels long and 5 pixels wide, into the intraretinal cystoid spaces, centrally onto unaffected peripheral RPE, and onto the prefoveolar vitreous. Light reflectivity in the retinal pigment epithelium (RPE), vitreous, and intraretinal spaces for the different retinal conditions and control subjects were compared. RESULTS: Reflectivities of the vitreous and the RPE were similar among the groups. Hyporeflective spaces in eyes with exudation (DME, RP, and CSC) had higher reflectivity compared with the mean reflectivity of the vitreous, whereas the cystoid spaces in the maculae of the eyes without exudation (CD and IPT) had a lower reflectivity than did the normal vitreous. CONCLUSIONS: Analysis of the light reflectivity profiles may be a tool to determine whether the density of hyporeflective spaces in the macula is greater or less than that of the vitreous, and may be a way to differentiate degenerative from exudative macular disease.

Experimental investigation of regular fluids and nanofluids during flow boiling in a single microchannel at different heat fluxes and mass fluxes

The dissipation of high heat flux from integrated circuit chips and the maintenance of acceptable junction temperatures in high powered electronics require advanced cooling technologies. One such technology is two-phase cooling in microchannels under confined flow boiling conditions. In macroscale flow boiling bubbles will nucleate on the channel walls, grow, and depart from the surface. In microscale flow boiling bubbles can fill the channel diameter before the liquid drag force has a chance to sweep them off the channel wall. As a confined bubble elongates in a microchannel, it traps thin liquid films between the heated wall and the vapor core that are subject to large temperature gradients. The thin films evaporate rapidly, sometimes faster than the incoming mass flux can replenish bulk fluid in the microchannel. When the local vapor pressure spike exceeds the inlet pressure, it forces the upstream interface to travel back into the inlet plenum and create flow boiling instabilities. Flow boiling instabilities reduce the temperature at which critical heat flux occurs and create channel dryout. Dryout causes high surface temperatures that can destroy the electronic circuits that use two-phase micro heat exchangers for cooling. Flow boiling instability is characterized by periodic oscillation of flow regimes which induce oscillations in fluid temperature, wall temperatures, pressure drop, and mass flux. When nanofluids are used in flow boiling, the nanoparticles become deposited on the heated surface and change its thermal conductivity, roughness, capillarity, wettability, and nucleation site density. It also affects heat transfer by changing bubble departure diameter, bubble departure frequency, and the evaporation of the micro and macrolayer beneath the growing bubbles. Flow boiling was investigated in this study using degassed, deionized water, and 0.001 vol% aluminum oxide nanofluids in a single rectangular brass microchannel with a hydraulic diameter of 229 µm for one inlet fluid temperature of 63°C and two constant flow rates of 0.41 ml/min and 0.82 ml/min. The power input was adjusted for two average surface temperatures of 103°C and 119°C at each flow rate. High speed images were taken periodically for water and nanofluid flow boiling after durations of 25, 75, and 125 minutes from the start of flow. The change in regime timing revealed the effect of nanoparticle suspension and deposition on the Onset of Nucelate Boiling (ONB) and the Onset of Bubble Elongation (OBE). Cycle duration and bubble frequencies are reported for different nanofluid flow boiling durations. The addition of nanoparticles was found to stabilize bubble nucleation and growth and limit the recession rate of the upstream and downstream interfaces, mitigating the spreading of dry spots and elongating the thin film regions to increase thin film evaporation.

Predator-prey: interactions between the spiny waterflea (Bythotrephes longimanus) and pumpkinseed sunfish (Lepomis gibbosus)

Small pumpkinseed sunfish ( Lepomis gibbosus), were found to be capable of removing the spine of Bythotrephes longimanus, an invasive cladoceran. Because fish consumption may be important in the dispersal or control of Bythotrephes, aquarium feeding experiments were conducted to 1) establish if the spine removal behavior of the pumpkinseeds was locally unique; 2) quantify how frequently pumpkinseeds exhibit the behavior; 3) determine if pumpkinseed handle Bythotrephes more quickly than other species of fish; and 4) verify if Bythotrephes' resting eggs pass through the digestive systems of pumpkinseeds in viable condition. The experiments revealed that pumpkinseeds (45-70 mm TL) from two geographic regions were more successful (100%) at removing Bythotrephes' spine, and handled Bythotrephes more quickly than yellow perch (Perca flavescens) (49-57 mm TL) and smallmouth bass (Micropterus dolomieu) (50-57mm TL) used in the study. Of 244 live Bythotrephes' resting eggs fed to the pumpkinseeds, 104 (42.6%) passed through their digestive systems. From those eggs, only 10 successfully hatched. Preliminary enclosure experiments were carried out and indicated that pumpkinseeds will consume Bythotrephes in natural settings. These findings provide new evidence that certain fish, with specialized morphology for prey manipulation, have the ability to influence the distribution and establishment of Bythotrephes.

Variations of the FEAST Eigenvalue Algorithm

FEAST is a recently developed eigenvalue algorithm which computes selected interior eigenvalues of real symmetric matrices. It uses contour integral resolvent based projections. A weakness is that the existing algorithm relies on accurate reasoned estimates of the number of eigenvalues within the contour. Examining the singular values of the projections on moderately-sized, randomly-generated test problems motivates orthogonalization-based improvements to the algorithm. The singular value distributions provide experimentally robust estimates of the number of eigenvalues within the contour. The algorithm is modified to handle both Hermitian and general complex matrices. The original algorithm (based on circular contours and Gauss-Legendre quadrature) is extended to contours and quadrature schemes that are recursively subdividable. A general complex recursive algorithm is implemented on rectangular and diamond contours. The accuracy of different quadrature schemes for various contours is investigated.

Compressive strength of interbody cages in the lumbar spine: the effect of cage shape, posterior instrumentation and bone density

One goal of interbody fusion is to increase the height of the degenerated disc space. Interbody cages in particular have been promoted with the claim that they can maintain the disc space better than other methods. There are many factors that can affect the disc height maintenance, including graft or cage design, the quality of the surrounding bone and the presence of supplementary posterior fixation. The present study is an in vitro biomechanical investigation of the compressive behaviour of three different interbody cage designs in a human cadaveric model. The effect of bone density and posterior instrumentation were assessed. Thirty-six lumbar functional spinal units were instrumented with one of three interbody cages: (1) a porous titanium implant with endplate fit (Stratec), (2) a porous, rectangular carbon-fibre implant (Brantigan) and (3) a porous, cylindrical threaded implant (Ray). Posterior instrumentation (USS) was applied to half of the specimens. All specimens were subjected to axial compression displacement until failure. Correlations between both the failure load and the load at 3 mm displacement with the bone density measurements were observed. Neither the cage design nor the presence of posterior instrumentation had a significant effect on the failure load. The loads at 3 mm were slightly less for the Stratec cage, implying lower axial stiffness, but were not different with posterior instrumentation. The large range of observed failure loads overlaps the potential in vivo compressive loads, implying that failure of the bone-implant interface may occur clinically. Preoperative measurements of bone density may be an effective tool to predict settling around interbody cages.

Cell shape-dependent early responses of fibroblasts to cyclic strain

Randomly spread fibroblasts on fibronectin-coated elastomeric membranes respond to cyclic strain by a varying degree of focal adhesion assembly and actin reorganization. We speculated that the individual shape of the cells, which is linked to cytoskeletal structure and pre-stress, might tune these integrin-dependent mechanotransduction events. To this aim, fibronectin circles, squares and rectangles of identical surface area (2000μm(2)) were micro-contact printed onto elastomeric substrates. Fibroblasts plated on these patterns occupied the corresponding shapes. Cyclic 10% equibiaxial strain was applied to patterned cells for 30min, and changes in cytoskeleton and cell-matrix adhesions were quantified after fluorescence staining. After strain, megakaryocytic leukemia-1 protein translocated to the nucleus in most cells, indicating efficient RhoA activation independently of cell shape. However, circular and square cells (with radial symmetry) showed a significantly greater increase in the number of actin stress fibers and vinculin-positive focal adhesions after cyclic strain than rectangular (bipolar) cells of identical size. Conversely, cyclic strain induced larger changes in pY397-FAK positive focal complexes and zyxin relocation from focal adhesions to stress fibers in bipolar compared to symmetric cells. Thus, radially symmetric cells responded to cyclic strain with a larger increase in assembly, whereas bipolar cells reacted with more pronounced reorganization of actin stress fibers and matrix contacts. We conclude that integrin-mediated responses to external mechanical strain are differentially modulated in cells that have the same spreading area but different geometries, and do not only depend on mere cell size.

Reply to the editor

We appreciate the comments and concerns expressed by Arakawa and colleagues regarding our article, titled “Pulsatile control of rotary blood pumps: Does the modulation waveform matter?”1 Unfortunately, we have to disagree with Arakawa and colleagues. As is obvious from the title of our article, it investigates the effect of different waveforms on the heart–device interaction. In contrast to the authors' claim, this is the first article in the literature that uses basic waveforms (sine, triangle, saw tooth, and rectangular) with different phase shifts to examines their impact on left ventricular unloading. The previous publications2, 3 and 4 just varied the pump speed during systole and diastole, which was first reported by Bearnson and associates5 in 1996, and studied its effect on aortic pressure, coronary flow, and end-diastolic volume. We should mention that dp/dtmax is a load-sensitive parameter of contractility and not representative for the degree of unloading. Moreover, none of the aforementioned reports has studied mechanical unloading and in particular the stroke work of the left ventricle. Our method is unique because we do not just alternate between high and low speed but have accurate control of the waveform because of the direct drive system of Levitronix Technologies LLC (Waltham, Mass) and a custom-developed pump controller. Without referring, Arakawa and associates state “several previous studies have already reported the coronary flow diminishes as the left ventricular assist device support increases.” It should be noted that all the waveforms used in our study have 2000 rpm average value with 1000 rpm amplitude, which is not an excessive speed for the CentriMag rotary pump (Levitronix) to collapse the ventricle and diminish the coronary flow. We agree with Arakawa and coworkers that there is a need for a heart failure model to come to more relevant results with respect to clinical expectations. However, we have explored many existing models, including species and breeds that have a native proneness to cardiomyopathy, but all of them differ from the genetic presentation in humans. We certainly do not believe that the use of microembolization, in which the coronary circulation is impaired by the injection of microspheres, would form a good model from which to draw conclusions about coronary flow change under different loading conditions. A model would be needed in which either an infarct is created to mimic ischemic heart failure or the coronary circulation remains untouched to simulate, for instance, dilated cardiomyopathy. Furthermore, in discussion we clearly mention that “lack of heart failure is a major limitation of our study.” We also believe that unloading is not the only factor of the cardiac functional recovery, and an excessive unloading of the left ventricle might lead to cardiac tissue atrophy. Therefore, in our article we mention that control of the level of cardiac unloading by assist devices has been suggested as a mechanical tool to promote recovery, and more studies are required to find better strategies for the speed modulation of rotary pumps and to achieve an optimal heart load control to enhance myocardial recovery. Finally, there are many publications about pulsing rotary blood pumps and it was impossible to include them all. We preferred to reference some of the earlier basic works such as an original research by Bearnson and coworkers5 and another article published by our group,6 which is more relevant.

The origin of oriented lakes: evidence from the Bolivian Amazon

The presence of hundreds of rectangular and oriented lakes is one of the most striking characteristics of the Llanos de Moxos (LM) landscape in the Bolivian Amazon. Oriented lakes also occur in the Arctic coastal plains of Russia, Alaska and Canada and along the Atlantic Coastal Plain from northeast Florida to southeast New Jersey and along the coast of northeast Brazil. Many different mechanisms have been proposed for their formation. In the LM, Plafker's (1964) tectonic model, in which subsidence results from the propagation of bedrock faults through the foreland sediments, is the most accepted. However, this model has not been verified. Here, we present new results from stratigraphic transects across the borders of three rectangular and oriented lakes in the LM. A paleosol buried under mid-Holocene sediments is used as a stratigraphic marker to assess the vertical displacement of sediments on both sides of the alleged faults. Our results show that there is no vertical displacement and, therefore, that Plafker's model can be ruled out. We suggest that, among all the proposed mechanisms behind lake formation, the combined action of wind and waves is the most likely. The evidence from the LM provides new hints for the formation of oriented lakes worldwide.

The timing of significant arch dimensional changes with fixed orthodontic appliances: data from a multicenter randomised controlled trial.

OBJECTIVES To identify the timing of significant arch dimensional increases during orthodontic alignment involving round and rectangular nickel-titanium (NiTi) wires and rectangular stainless steel (SS). A secondary aim was to compare the timing of changes occurring with conventional and self-ligating fixed appliance systems. METHODS In this non-primary publication, additional data from a multicenter randomised trial initially involving 96 patients, aged 16 years and above, were analysed. The main pre-specified outcome measures were the magnitude and timing of maxillary intercanine, interpremolar, and intermolar dimensions. Each participant underwent alignment with a standard Damon (Ormco, Orange, CA) wire sequence for a minimum of 34 weeks. Blinding of clinicians and patients was not possible; however, outcome assessors and data analysts were kept blind to the appliance type during data analysis. RESULTS Complete data were obtained from 71 subjects. Significant arch dimensional changes were observed relatively early in treatment. In particular, changes in maxillary inter-first and second premolar dimensions occurred after alignment with an 0.014in. NiTi wire (P<0.05). No statistical differences in transverse dimensions were found between rectangular NiTi and working SS wires for each transverse dimension (P>0.05). Bracket type had no significant effect on the timing of the transverse dimensional changes. CONCLUSIONS Arch dimensional changes were found to occur relatively early in treatment, irrespective of the appliance type. Nickel-titanium wires may have a more profound effect on transverse dimensions than previously believed. CLINICAL SIGNIFICANCE On the basis of this research orthodontic expansion may occur relatively early in treatment. Nickel-titanium wires may have a more profound effect on transverse dimensions than previously believed.

Insight into the dynamic ligand exchange process involved in bipyridyl linked arene ruthenium metalla-rectangles

The dynamic ligand exchange behavior of cationic arene ruthenium metalla-rectangles of the type [(pcymene) 4Ru4(OOXOO)2(NXN)2]4+ (OOXOO ¼ oxalato, 2,5-dioxydo-1,4-benzoquinonato, 5,8-dioxydo-1,4-naphthoquinonato; NXN ¼ 4,40-bipyridine-H8, 4,40-bipyridine-D8) has been studied in solution. The robustness of the rectangular architecture has been evidenced by NMR and ESI mass spectrometry. Thermodynamic and kinetic aspects of the ligand exchange process have been explored using 1H/2D isotope labeling of the 4,40-bipyridine connectors. This study shows that ligand exchange does not proceed spontaneously for these metalla-assemblies, even at high temperature, unless an external stimulus is applied.