Abstract Primary blast injury (PBI), which relates gross blast-related trauma or traces of injury in air-filled tissues or those tissues adjacent to air-filled regions (rupture/lesions, contusions, hemorrhaging), has been documented in a number of marine mammal species after blast exposure [1, 2, 3]. However, very little is known about marine mammal susceptibility to PBI except in rare cases of opportunistic studies. As a result, traditional techniques rely on analyses using small-scale terrestrial mammals as surrogates for large-scale marine mammals. For an In-house Laboratory Independent Research (ILIR) project sponsored by the Office of Naval Research (ONR), researchers at the Naval Undersea Warfare Center, Division Newport (NUWCDIVNPT), have undertaken a broad 3-year effort to integrate computational fluid-structure interaction techniques with marine mammal anatomical structure. The intent is to numerically simulate the dynamic response of a marine mammal thoracic cavity and air-filled lungs to shock loading, to enhance understanding of marine mammal lungs to shock loading in the underwater environment. In the absence of appropriate test data from live marine mammals, a crucial part of this work involves code validation to test data for a suitable surrogate test problem. This research employs a surrogate of an air-filled spherical membrane structure subjected to shock loading as a first order approximation to understanding marine mammal lung response to underwater explosions (UNDEX). This approach incrementally improves upon the currently used one-dimensional spherical air bubble approximation to marine mammal lung response by providing an encapsulating boundary for the air. The encapsulating structure is membranous, with minimal simplified representation not accounting for marine mammal species-specific and individual animal differences in tissue composition, rib mechanics, and mechanical properties of interior lung tissue. NUWCDIVNPT partnered with the Naval Submarine Medical Research Laboratory (NSMRL) to design and execute a set of experiments to investigate the shock response of an air-filled rubber dodgeball in a shallow underwater environment. These tests took place in the 2.13 m (7-ft) diameter pressure tank at the University of Rhode Island, with test measurements including pressure data and digital image correlation (DIC) data captured with high-speed cameras in a stereo setup. The authors developed 3-dimensional computational models of the dodgeball experiments using Dynamic System Mechanics Advanced Simulation (DYSMAS), a Navy fluid-structure interaction code. DYSMAS models of a variety of different problems involving submerged pressure vessel structures responding to hydrostatic and/or UNDEX loading have been validated against test data [4]. Proper validation of fluid structure interaction simulations is quite challenging, requiring measurements in both the fluid and structure domains. This paper details the development of metrics for comparison between test measurements and simulation results, with a discussion of potential sources of uncertainty.
Underwater implosion experiments were conducted with thin-wall glass spheres to determine the influence that structural failure has on the pressure pulse. Four experiments were conducted with glass spheres having an outside diameter of 7.62 cm, thickness of 0.762 mm, and an estimated buckling pressure of 7.57 MPa. The experiments were performed in a pressure vessel at a hydrostatic pressure of 6.996 MPa. The average peak pressure of the implosion pressure pulse was 26.1 MPa, measured at a radial distance of 10.16 cm from the sphere center. A computational fluid structure interaction model was developed to assess how the failure rate of the glass structure influences the pressure time history. The model employed a specified glass failure sequence that is uniform in time and space. It was found that for the conditions of the test, a glass failure rate of 275 m/s provided a reasonable representation of the test data. The test data and the model results show that the failure time history of the structure has a significant influence on an implosion pressure pulse. Computational prediction of an implosion pressure pulse needs to include the failure time history of the structure; otherwise it will overpredict the pressure time history.
This paper presents an experimental investigation of convective heat transfer for laminar gas flow through a microchannel. A test stand was set up to impose thermal boundary conditions of constant temperature gradient along the microchannel length. Additionally, thin film temperature sensors were developed and used to directly measure the microchannel surface temperature. Heat transfer experiments were conducted with laminar nitrogen gas flow, in which the outlet Ma was between 0.10 and 0.42. The experimental measurements of inlet and outlet gas temperature and the microchannel wall temperature were used to validate a two-dimensional numerical model for gaseous flow in microchannel. The model was then used to determine local values of Ma, Re, and Nu. The numerical results show that after the entrance region, Nu approaches 8.23, the fully developed value of Nu for incompressible flow for constant wall heat flux if Nu is defined based on (Tw−Tref) and plotted as a function of the new dimensionless axial length, X*=(x∕2H)(Ma2)∕(RePr).
Poiseuille number, the product of friction factor and Reynolds number (fRe) for quasi-fully-developed gas microchannel flow in the slip flow regime, was obtained numerically based on the arbitrary-Lagrangian-Eulerian method. Two-dimensional compressible momentum and energy equations were solved for a wide range of Reynolds and Mach numbers for constant wall temperatures that are lower or higher than the inlet temperature. The channel height ranges from 2 μm to 10 μm and the channel aspect ratio is 200. The stagnation pressure p stg is chosen such that the exit Mach number ranges from 0.1 to 1.0. The outlet pressure is fixed at atmospheric conditon. Mach and Knudsen numbers are systematically varied to determine their effects on fRe. The correlation for fRe for the slip flow is obtained from that of fRe of no-slip flow and incompressible theory as a function of Mach and Knudsen numbers. The results are in excellent agreement with the available experimental measurements. It was found that fRe is a function of Mach and Knudsen numbers and is different from the values by 96/(1+12Kn) obtained from the incompressible flow theory.
In this paper an experimental investigation of convective heat transfer is presented for laminar gas flow through a microchannel. A test stand was setup to impose thermal boundary conditions of constant temperature gradient along the microchannel length. Additionally, thin film temperature sensors were developed and used to directly measure the microchannel surface temperature. Heat transfer experiments were conducted in the laminar flow regime with the outlet Ma between 0.10 and 0.42. The experimental measurements of inlet and outlet gas temperature and the microchannel wall temperature were used to validate a 2D numerical model for gaseous flow in microchannel. The model was then used to determine bcal values of Ma, Re, and Nu. The numerical results show that after the entrance region, Nu approaches 8.23, the fully developed value of Nu for incompressible flow for constant wall heat flux if Nu is defined based on (Tw-Taw) and plotted as a function of the new dimensionless axial length, X* = (x/2H)(Ma2)/(RePr).
This paper presents an experimental investigation of laminar gas flow through microchannels. The independent variables: relative surface roughness, Knudsen number and Mach number were systematically varied to determine their influence on the friction factor. The microchannels were etched into silicon wafers, capped with glass, and have hydraulic diameters between 5 and 96 μm. The pressure was measured at seven locations along the channel length to determine local values of Knudsen number, Mach number and friction factor. All measurements were made in the laminar flow regime with Reynolds numbers ranging from 0.1 to 1000. The results show close agreement for the friction factor in the limiting case of low Ma and low Kn with the incompressible continuum flow theory. The effect of compressibility is observed to have a mild (8 percent) increase in the friction factor as the Mach number approaches 0.35. A 50 percent decrease in the friction factor was seen as the Knudsen number was increased to 0.15. Finally, the influence of surface roughness on the friction factor was shown to be insignificant for both continuum and slip flow regimes.
Two-dimensional compressible momentum and energy equations are solved in a parallel plate channel to obtain the effect of compressibility on gaseous flow characteristics in micro-channels. The numerical methodology is based on arbitrary-Lagrangian–Eulerian method. The computations were performed for a wide range of Reynolds number and Mach number and for both ‘no heat conduction’ and isothermal flow conditions. The results are presented in form of the product of friction factor and Reynolds number (f·Re) expressed as a function of both Reynolds number and Mach number and was compared with available experimental results. It was found that f·Re is mainly a function of Mach number and is different from the incompressible value of 96 for parallel-plate channels.
This paper presents an experimental investigation on compressible gas flow through microchannels with a constant area, rectangular cross-section. The microchannels are etched into silicon wafers, capped with smooth glass, and have hydraulic diameters between 4 and 100 μm. All measurements were made in the laminar flow regime with Reynolds numbers ranging from 0.02 to 1000. Smooth channels were obtained by etching (100) silicon wafers with potassium hydroxide (KOH) solution. Rough channel surfaces were obtained by etching (110) silicon wafers with KOH. The investigation shows that the friction factor for both smooth and rough microchannels compares closely with continuum theory.
This paper presents an experimental investigation on nitrogen and helium flow through microchannels etched in silicon with hydraulic diameters between 10 and 40 microns, and Reynolds numbers ranging from 0.3 to 600. The objectives of this research are (1) to fabricate microchannels with uniform surface roughness and local pressure measurement; (2) to determine the friction factor within the locally fully developed region of the microchannel; and (3) to evaluate the effect of surface roughness on momentum transfer by comparison with smooth microchannels. The friction factor results are presented as the product of friction factor and Reynolds number plotted against Reynolds number. The following conclusions have been reached in the present investigation: (1) microchannels with uniform corrugated surfaces can be fabricated using standard photolithographic processes; and (2) surface features with low aspect ratios of height to width have little effect on the friction factor for laminar flow in microchannels.
This paper presents an experimental investigation on nitrogen and helium flow in microchannels etched in silicon with hydraulic diameters of 9.7, 19.6, and 46.6 μm, and Reynolds numbers ranging from 0. 2 to 1000. The objectives of this research are (1) to measure the pressure distribution along the length of a microchannel; and (2) to determine the friction factor within the fully developed region of the microchannel. The pressure distribution is presented as absolute local pressure plotted against the distance from the microchannel inlet. The friction factor results are presented as the product of friction factor and Reynolds number plotted against Reynolds number with the outlet Knudsen number, Kn, as a curve parameter. The following conclusions have been reached in the present investigation: (1) Pressure losses at the microchannel entrance can be significant; (2) the product, f*Re, when measured sufficiently far away from the entrance and exit is a constant in the laminar flow region; and (3) the friction factor decreases as the Knudsen number increases.