Effects of hibernation on bone in yellow-bellied marmots

Disuse osteoporosis is a problem for people with spinal cord injury or stroke, patients confined to bed rest, and astronauts exposed to microgravity. Unlike most mammals however, bears have been shown to prevent bone loss during hibernation, a seasonal period of disuse. Similarly, studies in ground squirrels indicate preservation of whole bone strength during hibernation, though evidence suggests there may be some increased osteocytic osteolysis. Uncovering the mechanism by which these animals prevent bone loss during hibernation could lead to an improved treatment for osteoporosis in humans. Marmots are a good animal model for these studies because they are small enough to easily house in an animal facility yet still utilize intracortical remodeling like humans and bears, and unlike smaller rodents like squirrels. Marmots preserve bone mechanical and microstructural properties during hibernation. Bone mechanical and geometrical properties are not diminished in post-hibernation samples compared to pre-hibernation samples. Mineral content, measured by ash fraction, was higher in post-hibernation samples (p = 0.0003). Haversian porosity as well as remodeling cavity density were not different (p > 0.38) between pre- and post-hibernation samples. Similarly, average lacunar area, lacunar density, and lacunar porosity were all lower (p < 0.0001) in post-hibernation samples. Trabecular thickness was larger in posthibernation samples (p = 0.0058). Bone volume fraction was not different between groups, but approached significance (p = 0.0725). Further studies in marmots and other hibernators could help uncover the mechanism that allows hibernators to prevent disuse osteoporosis during hibernation.

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.

Decentralized or centralized production : impacts to the environment, industry, and the economy

Since product take-back is mandated in Europe, and has effects for producers worldwide including the U.S., designing efficient forward and reverse supply chain networks is becoming essential for business viability. Centralizing production facilities may reduce costs but perhaps not environmental impacts. Decentralizing a supply chain may reduce transportation environmental impacts but increase capital costs. Facility location strategies of centralization or decentralization are tested for companies with supply chains that both take back and manufacture products. Decentralized and centralized production systems have different effects on the environment, industry and the economy. Decentralized production systems cluster suppliers within the geographical market region that the system serves. Centralized production systems have many suppliers spread out that meet all market demand. The point of this research is to help further the understanding of company decision-makers about impacts to the environment and costs when choosing a decentralized or centralized supply chain organizational strategy. This research explores; what degree of centralization for a supply chain makes the most financial and environmental sense for siting facilities; and which factories are in the best location to handle the financial and environmental impacts of particular processing steps needed for product manufacture. This research considered two examples of facility location for supply chains when products are taken back; the theoretical case involved shoe resoling and a real world case study considered the location of operations for a company that reclaims multiple products for use as material inputs. For the theoretical example a centralized strategy to facility location was optimal: whereas for the case study a decentralized strategy to facility location was best. In conclusion, it is not possible to say that a centralized or decentralized strategy to facility location is in general best for a company that takes back products. Each company’s specific concerns, needs, and supply chain details will determine which degree of centralization creates the optimal strategy for siting their facilities.

Analysis of a novel waste heat recovery mechanism for an I.C. engine

Typical internal combustion engines lose about 75% of the fuel energy through the engine coolant, exhaust and surface radiation. Most of the heat generated comes from converting the chemical energy in the fuel to mechanical energy and in turn thermal energy is produced. In general, the thermal energy is unutilized and thus wasted. This report describes the analysis of a novel waste heat recovery (WHR) system that operates on a Rankine cycle. This novel WHR system consists of a second piston within the existing piston to reduce losses associated with compression and exhaust strokes in a four-cycle engine. The wasted thermal energy recovered from the coolant and exhaust systems generate a high temperature and high pressure working fluid which is used to power the modified piston assembly. Cycle simulation shows that a large, stationary natural gas spark ignition engine produces enough waste heat to operate the novel WHR system. With the use of this system, the stationary gas compression ignition engine running at 900 RPM and full load had a net increase of 177.03 kW (240.7 HP). This increase in power improved the brake fuel conversion efficiency by 4.53%.

Emission cross sections for neutral xenon impacted by Xe+ and Xe2+

Ion impact emission cross sections for eleven transitions from the 5p56p configuration to the 5p56s configuration of neutral xenon occurring in the spectral region between 700 nm and 1000 nm have been measured experimentally. Collisions between both singly- and doublyionized xenon and neutral xenon have been studied. These cross sections are of primary use in the development of a spectrographic diagnostic for Hall effect thruster plasma. A detailed discussion of the experimental methods and the subsequent data reduction is included. The results are presented and the importance of these data for spectrographic emission models of Hall effect thruster plasmas is discussed.

Bislug flow in circular and noncircular channels and the role of interface stretching on energy dissipation

The area of microfluidics has increased in popularity with such fields as MEMS, microreactors, microscaleheat exchangers, etc. A comprehensive understanding of dissipation mechanisms for fluid flow in microchannels is required to accurately predict the behavior in these small systems. Tests were performed using a constant pressure potential created by two immiscible fluids juxtaposed in a microchannel. This study focused on the flow and dissipation mechanisms in round and square microchannels. There are four major dissipation mechanisms in slug flow; wall shear, dissipation at the contact line, menisci interaction and the stretching of the interface. A force balance between the internal driving potential, viscous drag and interface stretching was used to develop a model for the prediction of the velocity of a bislug in a microchannel. Interface stretching is a dissipation mechanism that has been included due to the unique system properties and becomes increasingly more important as the bislug decreases in length.

Synthesis of computations and experiments for obtaining pulsatile gas flow rates from dynamic pressure difference measurements across an orifice-plate meter

The use of conventional orifice-plate meter is typically restricted to measurements of steady flows. This study proposes a new and effective computational-experimental approach for measuring the time-varying (but steady-in-the-mean) nature of turbulent pulsatile gas flows. Low Mach number (effectively constant density) steady-in-the-mean gas flows with large amplitude fluctuations (whose highest significant frequency is characterized by the value fF) are termed pulsatile if the fluctuations have a direct correlation with the time-varying signature of the imposed dynamic pressure difference and, furthermore, they have fluctuation amplitudes that are significantly larger than those associated with turbulence or random acoustic wave signatures. The experimental aspect of the proposed calibration approach is based on use of Coriolis-meters (whose oscillating arm frequency fcoriolis >> fF) which are capable of effectively measuring the mean flow rate of the pulsatile flows. Together with the experimental measurements of the mean mass flow rate of these pulsatile flows, the computational approach presented here is shown to be effective in converting the dynamic pressure difference signal into the desired dynamic flow rate signal. The proposed approach is reliable because the time-varying flow rate predictions obtained for two different orifice-plate meters exhibit the approximately same qualitative, dominant features of the pulsatile flow.

Effect of segmented anodes on the performance and plume of a hall thruster

Development of alternative propellants for Hall thruster operation is an active area of research. Xenon is the current propellant of choice for Hall thrusters, but can be costly in large thrusters and for extended test periods. Condensible propellants may offer an alternative to xenon, as they will not require costly active pumping to remove from a test facility, and may be less expensive to purchase. A method has been developed which uses segmented electrodes in the discharge channel of a Hall thruster to divert discharge current to and from the main anode and thus control the anode temperature. By placing a propellant reservoir in the anode, the evaporation rate, and hence, mass flow of propellant can be controlled. Segmented electrodes for thermal control of a Hall thruster represent a unique strategy of thruster design, and thus the performance of the thruster must be measured to determine the effect the electrodes have on the thruster. Furthermore, the source of any changes in thruster performance due to the adjustment of discharge current between the shims and the main anode must be characterized. A Hall thruster was designed and constructed with segmented electrodes. It was then tested at anode voltages between 300 and 400 V and mass flows between 4 and 6 mg/s, as well as 100%, 75%, 50%, 25%, and <5% of the discharge current on the shim electrodes. The level of current on the shims was adjusted by changing the shim voltage. At each operating point, the thruster performance, plume divergence, ion energy, and multiply charged ion fraction were measured performance exhibited a small change with the level of discharge current on the shim electrodes. Thrust and specific impulse increased by as much as 6% and 7.7%, respectively, as discharge current was shifted from the main anode to the shims at constant anode voltage. Thruster efficiency did not change. Plume divergence was reduced by approximately 4 degrees of half-angle at high levels of current on the shims and at all combinations of mass flow and anode voltage. The fraction of singly charged xenon in the thruster plume varied between approximately 80% and 95% as the anode voltage and mass flow were changed, but did not show a significant change with shim current. Doubly and triply charged xenon made up the remainder of the ions detected. Ion energy exhibited a mixed behavior. The highest voltage present in the thruster largely dictated the most probable energy; either shim or anode voltage, depending on which was higher. The overall change in most probable ion energy was 20-30 eV, the majority of which took place while the shim voltage was higher than the anode voltage. The thrust, specific impulse, plume divergence, and ion energy all indicate that the thruster is capable of a higher performance output at high levels of discharge current on the shims. The lack of a change in efficiency and fraction of multiply charged ions indicate that the thruster can be operated at any level of current on the shims without detrimental effect, and thus a condensible propellant thruster can control the anode temperature without a decrease in efficiency or a change in the multiply charged ion fraction.

Implementation of a phenomenological evaporation model into a porous network simulation for water management in low temperature fuel cells

A phenomenological transition film evaporation model was introduced to a pore network model with the consideration of pore radius, contact angle, non-isothermal interface temperature, microscale fluid flows and heat and mass transfers. This was achieved by modeling the transition film region of the menisci in each pore throughout the porous transport layer of a half-cell polymer electrolyte membrane (PEM) fuel cell. The model presented in this research is compared with the standard diffusive fuel cell modeling approach to evaporation and shown to surpass the conventional modeling approach in terms of predicting the evaporation rates in porous media. The current diffusive evaporation models used in many fuel cell transport models assumes a constant evaporation rate across the entire liquid-air interface. The transition film model was implemented into the pore network model to address this issue and create a pore size dependency on the evaporation rates. This is accomplished by evaluating the transition film evaporation rates determined by the kinetic model for every pore containing liquid water in the porous transport layer (PTL). The comparison of a transition film and diffusive evaporation model shows an increase in predicted evaporation rates for smaller pore sizes with the transition film model. This is an important parameter when considering the micro-scaled pore sizes seen in the PTL and becomes even more substantial when considering transport in fuel cells containing an MPL, or a large variance in pore size. Experimentation was performed to validate the transition film model by monitoring evaporation rates from a non-zero contact angle water droplet on a heated substrate. The substrate was a glass plate with a hydrophobic coating to reduce wettability. The tests were performed at a constant substrate temperature and relative humidity. The transition film model was able to accurately predict the drop volume as time elapsed. By implementing the transition film model to a pore network model the evaporation rates present in the PTL can be more accurately modeled. This improves the ability of a pore network model to predict the distribution of liquid water and ultimately the level of flooding exhibited in a PTL for various operating conditions.

Micro combined heat and power laboratory development

Micro Combined Heat and Power (Micro-CHP) system produces both electricity and heat required for residential or small business applications. Use of Micro-CHP in a residential application not only creates energy and economic savings but also reduces the carbon foot print of the house or small business. Additionally, micro-CHP can subsidize its cost of operation by selling excess electricity produced back to the grid. Even though Micro-CHP remains attractive on paper, high initial cost and optimization issues in residential scale heat and electrical requirement has kept this technology from becoming a success. To understand and overcome all disadvantages posed my Micro-CHP system, a laboratory is developed to test different scenarios of Micro-CHP applications so that we can learn and improve the current technology. This report focuses on the development of this Micro-CHP laboratory including installation of Ecopower micro-CHP unit, developing fuel line and exhaust line for Ecopower unit, design of electrical and thermal loop, installing all the instrumentation required for data collection on the Ecopower unit and developing controls for heat load simulation using thermal loop. Also a simulation of Micro-CHP running on Syngas is done in Matlab. This work was supported through the donation of ‘Ecopower’ a Micro-CHP unit by Marathon Engine and through the support of Michigan Tech REF-IF grand.

Adaptive control of sinusoidal brushless DC motor actuators

Electrical Power Assisted Steering system (EPAS) will likely be used on future automotive power steering systems. The sinusoidal brushless DC (BLDC) motor has been identified as one of the most suitable actuators for the EPAS application. Motor characteristic variations, which can be indicated by variations of the motor parameters such as the coil resistance and the torque constant, directly impart inaccuracies in the control scheme based on the nominal values of parameters and thus the whole system performance suffers. The motor controller must address the time-varying motor characteristics problem and maintain the performance in its long service life. In this dissertation, four adaptive control algorithms for brushless DC (BLDC) motors are explored. The first algorithm engages a simplified inverse dq-coordinate dynamics controller and solves for the parameter errors with the q-axis current (iq) feedback from several past sampling steps. The controller parameter values are updated by slow integration of the parameter errors. Improvement such as dynamic approximation, speed approximation and Gram-Schmidt orthonormalization are discussed for better estimation performance. The second algorithm is proposed to use both the d-axis current (id) and the q-axis current (iq) feedback for parameter estimation since id always accompanies iq. Stochastic conditions for unbiased estimation are shown through Monte Carlo simulations. Study of the first two adaptive algorithms indicates that the parameter estimation performance can be achieved by using more history data. The Extended Kalman Filter (EKF), a representative recursive estimation algorithm, is then investigated for the BLDC motor application. Simulation results validated the superior estimation performance with the EKF. However, the computation complexity and stability may be barriers for practical implementation of the EKF. The fourth algorithm is a model reference adaptive control (MRAC) that utilizes the desired motor characteristics as a reference model. Its stability is guaranteed by Lyapunov’s direct method. Simulation shows superior performance in terms of the convergence speed and current tracking. These algorithms are compared in closed loop simulation with an EPAS model and a motor speed control application. The MRAC is identified as the most promising candidate controller because of its combination of superior performance and low computational complexity. A BLDC motor controller developed with the dq-coordinate model cannot be implemented without several supplemental functions such as the coordinate transformation and a DC-to-AC current encoding scheme. A quasi-physical BLDC motor model is developed to study the practical implementation issues of the dq-coordinate control strategy, such as the initialization and rotor angle transducer resolution. This model can also be beneficial during first stage development in automotive BLDC motor applications.

Integrated solar energy and absorption cooling model for HVAC (heating, ventilating, and air conditioning) applications in buildings

The demands in production and associate costs at power generation through non renewable resources are increasing at an alarming rate. Solar energy is one of the renewable resource that has the potential to minimize this increase. Utilization of solar energy have been concentrated mainly on heating application. The use of solar energy in cooling systems in building would benefit greatly achieving the goal of non-renewable energy minimization. The approaches of solar energy heating system research done by initiation such as University of Wisconsin at Madison and building heat flow model research conducted by Oklahoma State University can be used to develop and optimize solar cooling building system. The research uses two approaches to develop a Graphical User Interface (GUI) software for an integrated solar absorption cooling building model, which is capable of simulating and optimizing the absorption cooling system using solar energy as the main energy source to drive the cycle. The software was then put through a number of litmus test to verify its integrity. The litmus test was conducted on various building cooling system data sets of similar applications around the world. The output obtained from the software developed were identical with established experimental results from the data sets used. Software developed by other research are catered for advanced users. The software developed by this research is not only reliable in its code integrity but also through its integrated approach which is catered for new entry users. Hence, this dissertation aims to correctly model a complete building with the absorption cooling system in appropriate climate as a cost effective alternative to conventional vapor compression system.

Automated optimization of polymer extrusion dies using finite element analysis

An extrusion die is used to continuously produce parts with a constant cross section; such as sheets, pipes, tire components and more complex shapes such as window seals. The die is fed by a screw extruder when polymers are used. The extruder melts, mixes and pressures the material by the rotation of either a single or double screw. The polymer can then be continuously forced through the die producing a long part in the shape of the die outlet. The extruded section is then cut to the desired length. Generally, the primary target of a well designed die is to produce a uniform outlet velocity without excessively raising the pressure required to extrude the polymer through the die. Other properties such as temperature uniformity and residence time are also important but are not directly considered in this work. Designing dies for optimal outlet velocity variation using simple analytical equations are feasible for basic die geometries or simple channels. Due to the complexity of die geometry and of polymer material properties design of complex dies by analytical methods is difficult. For complex dies iterative methods must be used to optimize dies. An automated iterative method is desired for die optimization. To automate the design and optimization of an extrusion die two issues must be dealt with. The first is how to generate a new mesh for each iteration. In this work, this is approached by modifying a Parasolid file that describes a CAD part. This file is then used in a commercial meshing software. Skewing the initial mesh to produce a new geometry was also employed as a second option. The second issue is an optimization problem with the presence of noise stemming from variations in the mesh and cumulative truncation errors. In this work a simplex method and a modified trust region method were employed for automated optimization of die geometries. For the trust region a discreet derivative and a BFGS Hessian approximation were used. To deal with the noise in the function the trust region method was modified to automatically adjust the discreet derivative step size and the trust region based on changes in noise and function contour. Generally uniformity of velocity at exit of the extrusion die can be improved by increasing resistance across the die but this is limited by the pressure capabilities of the extruder. In optimization, a penalty factor that increases exponentially from the pressure limit is applied. This penalty can be applied in two different ways; the first only to the designs which exceed the pressure limit, the second to both designs above and below the pressure limit. Both of these methods were tested and compared in this work.

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.