Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time reversal flow”

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.

At least 145 records · Page 8

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗

Integrated CMOS photodetectors and signal processing for very low-level chemical sensing with the bioluminescent bioreporter integrated circuit

We report an integrated CMOS microluminometer optimized for the detection of low-level bioluminescence as part of the bioluminescent bioreporter integrated circuit (BBIC). This microluminometer improves on previous devices through careful management of the sub-femtoampere currents, both signal and leakage, that flow in the front-end processing circuitry. In particular, the photodiode is operated with a reverse bias of only a few mV, requiring special attention to the reset circuitry of the current-to-frequency converter (CFC) that forms the front-end circuit. We report a sub-femtoampere leakage current and a minimum detectable signal (MDS) of 0.15 fA (1510 s integration time) using a room temperature 1.47 mm2 CMOS photodiode. This microluminometer can detect luminescence from as few as 5000 fully induced Pseudomonas fluorescens 5RL bacterial cells. c2002 Elsevier Science B.V. All rights reserved.

Non-NASA Center↗

Integrated Bosch Process System Models for In-Situ Oxygen and Carbon Production

In-Situ Resource Utilization (ISRU) technology is a vital component to NASA’s mission of a sustainable presence on the Moon and Mars. Local resources can be leveraged to reduce resupply frequency and mass. Elements of the Bosch process, combined with the carbothermal reduction process, can produce oxygen on the lunar surface with minimal consumables. The Bosch process can also produce oxygen on the Martian surface by using the CO 2 -rich environment. Between both systems, adsorption pump, solar thermal energy, carbon formation reactor, and water recovery subsystems are modeled and integrated to create a functional model in MATLAB software. The model is used to simulate performance of the system and reduce mass, power, and volume requirements. This integrated system model provides a tool to scale ISRU technologies for oxygen and carbon production. The MATLAB model is created by developing a system of independent subsystem models that are solved for their quasi-steady state values which can be integrated with respect to time to determine the change in current states. A flexible time stepping method is used to ensure a high level of accuracy during periods of rapid change while still making use of a simple explicit integration method. The flexible time step is calculated for each independent subsystem and the minimum value from those is used as the overall time step. A flexible time step is calculated by dividing a resolution value, or the maximum change per time step, by the variables current rate of change. The maximum value from all points in space is used for subsystem models that contain multiple values. The process is done for every variable that is being monitored in each subsystem and the global minimum is used as that iteration’s timestep. Several assumptions used in the MATLAB model for fluid flow dynamics, such as 1-D gas flow through the sorption pump, are supported by modeling in Ansys Fluent software. The Lunar oxygen production system is outlined in Fig. 1. The carbothermal reduction subsystem uses solar energy to heat a mixture of lunar regolith and carbon powder to produce carbon monoxide. To begin, the carbon monoxide feeds to the modified Bosch subsystem along with hydrogen gas. The reactants then enter the carbon formation reactor where water and carbon powder are produced. Solar thermal energy is used to add energy to the reactor, but waste heat from the carbothermal process is another potential heat source. The water is collected and electrolyzed to produce hydrogen which reenters the Bosch subsystem, and the oxygen is stored for downstream use. The carbon powder is collected and feeds back into the carbothermal subsystem. The Martian oxygen production system uses the full Bosch process and is outlined in Fig 2. A CO 2 adsorption pump thermally cycles to scrub and pressurize CO 2 from the environment. Along with an initial supply of hydrogen, the reactants enter the Reverse Water Gas Shift Reactor (RWGSR) which produces carbon monoxide and water. Carbon monoxide and unreacted hydrogen enter the carbon formation reactor to produce water and carbon powder. The water is collected from both reactors and electrolyzed to reintroduce hydrogen and store oxygen for propellant production or life support. Carbon is removed from the carbon formation reactor and stored. The adsorption pump utilizes rapid cycle temperature swings within a stack of zeolite coated surfaces. The subsystem model solves 1-D quasi-steady conservation laws of the quasi-steady form, shown in Eq. 1, for the gas stream and heat exchange liquid to predict performance parameters such as breakthrough capacity and optimum cycle time. The source term S is used to capture interactions between the fluid flows and the sorbent. A quasi-steady-state scheme is used where no time derivatives appear in the governing equations, except for those in the source terms. This results in an autonomous system, where ∂F/∂x = ƒ(F). The fluxes F are provided at the inlet, and an explicit method is used to solve for the spatial distribution of F. The heat and mass flows to the sorbent are then extracted from the source terms. These flows are numerically integrated to produce a 1-D solution for the system’s state as a function of both time and space. The body of the adsorption pump is separated into two semi-independent models: the heat exchanger fluid flow and gas flow through the zeolite coated surfaces. Both models are solved using the above-described method to find a 1-D solution as a function of space and interact only once a timestep is taken. The interaction point is the sorbent through which all heat transfer between the two models must occur. Sorbent mass adsorption is calculated using the Lagergren model, shown in Eq. 2, where the transfer coefficient, λ D , is found by solving a system of nondimensionalized equations derived by using the heat and mass transfer analogy for transport phenomena. Using Grade 544 Type 13X zeolite as the sorbent material, the equilibrium concentration, θ eq , is calculated using the k-site Langmuir isotherm and fit parameters. Additionally, the enthalpy of adsorption used in the model is computed by interpolation of available data [1]. The subsystem model was validated using the Rapid Cycle Temperature Swing Adsorption (RC-TSA) pump. The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC). The solar thermal energy subsystem focuses on a solar concentrator concept with a heat exchanger to heat the reactants before entering the carbon formation reactor. The subsystem model assumes a fixed solar flux and reflector efficiency to calculate the reactant temperature given the incoming temperature, pressure, and exchanger geometry. The receiver is a custom manufactured series of copper blocks with serpentine channels to increase its surface area and the residence time of the reactants to heat up to 550 °C. The subsystem model was validated using a heat exchanger developed at NASA Glenn Research Center (GRC).

In situ Resource Utilization↗

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗

The Effects of Thermal History on Nucleation of Tetragonal Lysozyme Crystals, or Hot Protein and Cold Nucleation

Chicken egg white lysozyme has a well characterized thermally driven phase transition. Between pH 4.2 and 5.2, the transition temperature, as defined by the point where the tetragonal and orthorhombic solubilities are equal, is a function of the pH, salt (precipitant) type and concentration, and most likely of the buffer concentration as well. This phase transition can be carried out with protein solution alone, prior to addition of precipitant solution. Warming a lysozyme solution above the phase transition point, then cooling it back below this point, has been shown to affect the subsequent nucleation rate, as determined by the numbers and size of crystals formed, but not the growth rate for the tetragonal crystal form . We have now measured the kinetics of this process and investigated its reversibility. The transition effects are progressive with temperature, having a half time of about 1 hour at 37C at pH 4.8. After holding a lysozyme solution at 37C (prior to addition of precipitant) for 16 hours, then cooling it back to 4C no return to the pre-warmed nucleation kinetics are observed after at least 4 weeks. Orthorhombic lysozyme crystals apparently do not undergo the flow-induced growth cessation of tetragonal lysozyme crystals. Putting the protein in the orthorhombic form does not affect the averaged face growth kinetics, only nucleation, for tetragonal crystals. This differential behaviour may be exploited to elucidate how and where flow affects the lysozyme crystal growth process. The presentation will focus on the results of these and ongoing studies in this area.

Burke, Michael↗

Torque characteristics of a 122-centimeter butterfly valve with a hydro/pneumatic actuator

Actuating torque data from field testing of a 122-centimeter (48 in.) butterfly valve with a hydro/pneumatic actuator is presented. The hydraulic cylinder functions as either a forward or a reverse brake. Its resistance torque increases when the valve speeds up and decreases when the valve slows down. A reduction of flow resistance in the hydraulic flow path from one end of the hydraulic cylinder to the other will effectively reduce the hydraulic resistance torque and hence increase the actuating torque. The sum of hydrodynamic and friction torques (combined resistance torque) of a butterfly valve is a function of valve opening time. An increase in the pneumatic actuating pressure will result in a decrease in both the combined resistance torque and the actuator opening torque; however, it does shorten the valve opening time. As the pneumatic pressure increases, the valve opening time for a given configuration approaches an asymptotical value.

Lin, F. N.↗

An Active Broad Area Cooling Model of a Cryogenic Propellant Tank with a Single Stage Reverse Turbo-Brayton Cycle Cryocooler

As focus shifts towards long-duration space exploration missions, an increased interest in active thermal control of cryogenic propellants to achieve zero boil-off of cryogens has emerged. An active thermal control concept of considerable merit is the integration of a broad area cooling system for a cryogenic propellant tank with a combined cryocooler and circulator system that can be used to reduce or even eliminate liquid cryogen boil-off. One prospective cryocooler and circulator combination is the reverse turbo-Brayton cycle cryocooler. This system is unique in that it has the ability to both cool and circulate the coolant gas efficiently in the same loop as the broad area cooling lines, allowing for a single cooling gas loop, with the primary heat rejection occurring by way of a radiator and/or aftercooler. Currently few modeling tools exist that can size and characterize an integrated reverse turbo-Brayton cycle cryocooler in combination with a broad area cooling design. This paper addresses efforts to create such a tool to assist in gaining a broader understanding of these systems, and investigate their performance in potential space missions. The model uses conventional engineering and thermodynamic relationships to predict the preliminary design parameters, including input power requirements, pressure drops, flow rate, cycle performance, cooling lift, broad area cooler line sizing, and component operating temperatures and pressures given the cooling load operating temperature, heat rejection temperature, compressor inlet pressure, compressor rotational speed, and cryogenic tank geometry. In addition, the model allows for the preliminary design analysis of the broad area cooling tubing, to determine the effect of tube sizing on the reverse turbo-Brayton cycle system performance. At the time this paper was written, the model was verified to match existing theoretical documentation within a reasonable margin. While further experimental data is needed for full validation, this tool has already made significant steps towards giving a clearer understanding of the performance of a reverse turbo-Brayton cycle cryocooler integrated with broad area cooling technology for zero boil-off active thermal control.

Guzik, Monica C.↗

Jet Noise Diagnostics Supporting Statistical Noise Prediction Methods

The primary focus of my presentation is the development of the jet noise prediction code JeNo with most examples coming from the experimental work that drove the theoretical development and validation. JeNo is a statistical jet noise prediction code, based upon the Lilley acoustic analogy. Our approach uses time-average 2-D or 3-D mean and turbulent statistics of the flow as input. The output is source distributions and spectral directivity. NASA has been investing in development of statistical jet noise prediction tools because these seem to fit the middle ground that allows enough flexibility and fidelity for jet noise source diagnostics while having reasonable computational requirements. These tools rely on Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) solutions as input for computing far-field spectral directivity using an acoustic analogy. There are many ways acoustic analogies can be created, each with a series of assumptions and models, many often taken unknowingly. And the resulting prediction can be easily reverse-engineered by altering the models contained within. However, only an approach which is mathematically sound, with assumptions validated and modeled quantities checked against direct measurement will give consistently correct answers. Many quantities are modeled in acoustic analogies precisely because they have been impossible to measure or calculate, making this requirement a difficult task. The NASA team has spent considerable effort identifying all the assumptions and models used to take the Navier-Stokes equations to the point of a statistical calculation via an acoustic analogy very similar to that proposed by Lilley. Assumptions have been identified and experiments have been developed to test these assumptions. In some cases this has resulted in assumptions being changed. Beginning with the CFD used as input to the acoustic analogy, models for turbulence closure used in RANS CFD codes have been explored and compared against measurements of mean and rms velocity statistics over a range of jet speeds and temperatures. Models for flow parameters used in the acoustic analogy, most notably the space-time correlations of velocity, have been compared against direct measurements, and modified to better fit the observed data. These measurements have been extremely challenging for hot, high speed jets, and represent a sizeable investment in instrumentation development. As an intermediate check that the analysis is predicting the physics intended, phased arrays have been employed to measure source distributions for a wide range of jet cases. And finally, careful far-field spectral directivity measurements have been taken for final validation of the prediction code. Examples of each of these experimental efforts will be presented. The main result of these efforts is a noise prediction code, named JeNo, which is in middevelopment. JeNo is able to consistently predict spectral directivity, including aft angle directivity, for subsonic cold jets of most geometries. Current development on JeNo is focused on extending its capability to hot jets, requiring inclusion of a previously neglected second source associated with thermal fluctuations. A secondary result of the intensive experimentation is the archiving of various flow statistics applicable to other acoustic analogies and to development of time-resolved prediction methods. These will be of lasting value as we look ahead at future challenges to the aeroacoustic experimentalist.

Bridges, James E.↗

DC and Structured Electric Fields Observed on the C/NOFS Satellite and Their Association with Longitude, Plasma Density, and Solar Activity

Observations of DC electric fields and associated E x B plasma drifts gathered by the Vector Electric Field Investigation (VEFI) on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite are presented. We show statistical averages of the vector fields and resulting E x B plasma flows for the first three years of operations as a function of season, longitude, local time, and Fl 0.7 conditions. Magnetic field data from the VEFI science magnetometer are used to compute the plasma flows. Although typically displaying eastward and outward-directed fields during the day and westward and downward-directed fields at night, the data from DC electric field detector often reveal variations from this pattern that depend on longitude, solar activity, and plasma density. Clear "wave-4" tidal effects in both electric field components have been detected and will be presented. Zonal plasma drifts show a marked variation with solar activity and may be used as a proxy for neutral winds at night. Evidence for pre-reversal enhancements in the meridional drifts that depend on solar activity is present for some longitudes, and are corroborated by clear evidence in the plasma density data that the spacecraft journeyed below the F-peak during evenings when the rise in the ionosphere is most pronounced. In addition to DC electric fields, the data reveal considerable electric field structures at large scales (approx 100's of km) that are usually confined to the nightside. Although such electric field structures are typically associated with plasma density depletions and structures, what is surprising is the number of cases in which large amplitude, structured DC electric fields are observed without a significant plasma density counterpart structure, including their appearance at times when the ambient plasma density appears relatively quiescent. We investigate the mapping of structured electric fields along magnetic field lines from distant locations and consider tropospheric thunderstorm sources to explain some of the observations. The dependence of the structuring on season and solar activity will be reported and discussed.

Pfaff, Robert↗

Performance of a three-dimensional Navier-Stokes code on CYBER 205 for high-speed juncture flows

A vectorized 3D Navier-Stokes code has been implemented on CYBER 205 for solving the supersonic laminar flow over a swept fin/flat plate junction. The code extends MacCormack's predictor-corrector finite volume scheme to a generalized coordinate system in a locally one dimensional time split fashion. A systematic parametric study is conducted to examine the effect of fin sweep on the computed flow field. Calculated results for the pressure distribution on the flat plate and fin leading edge are compared with the experimental measurements of a right angle blunt fin/flat plate junction. The decrease in the extent of the separated flow region and peak pressure on the fin leading edge, and weakening of the two reversed supersonic zones with increase in fin sweep have been clearly observed in the numerical simulation.

Lakshmanan, B.↗

Forebody Aerodynamics of the F-18 High Alpha Research Vehicle with Actuated Forebody Strakes

Extensive pressure measurements and off-surface flow visualization were obtained on the forebody and strakes of the NASA F-18 High Alpha Research Vehicle (HARV) equipped with actuated forebody strakes. Forebody yawing moments were obtained by integrating the circumferential pressures on the forebody and strakes. Results show that large yawing moments can be generated with forebody strakes. At a 50 -angle-of-attack, deflecting one strake at a time resulted in a forebody yawing moment control reversal for small strake deflection angles. However, deflecting the strakes differentially about a 20 symmetric strake deployment eliminated the control reversal and produced a near linear variation of forebody yawing moment with differential strake deflection. At an angle of attack of 50 and for 0 and 20 symmetric strake deployments, a larger forebody yawing moment was generated by the forward fuselage (between the radome and the apex of the leading-edge extensions) than on the radome where the actuated forebody strakes were located. Cutouts on the flight vehicle strakes that were not on the wind tunnel models are believed to be responsible for deficits in the suction peaks on the flight radome pressure distributions and differences in the forebody yawing moments.

Fisher, David F.↗

Satellite observations of zonal electric fields near sunrise in the equatorial ionosphere

We report here on a number of examples of anomalous enhancements of eastward electric fields near sunrise in the equatorial ionospheric F-region. These examples were selected from the data base of the equatorial satellite, San Marco D (1988), which measured ionospheric electric fields during a period of solar minimum. The eastward electric fields reported correspond to vertical plasma drifts. The examples studied here are similar in signature and polarity to the pre-reversal electric field enhancements seen near sunset from ground-based radar systems. The morphology of these sunrise events, which are observed on about 14% of the morning-side satellite passes, are studied as a function of local zonal velocity, magnetic activity, geographic longitude and altitude. The nine events studied occur at locations where the zonal plasma flow is generally measured to be eastward, but reducing as a function of local time and at satellite longitudes where the magnetic declination has the opposite polarity as the declination of the sunrise terminator.

Aggson, T. L.↗

Effect of Actuated Forebody Strakes on the Forebody Aerodynamics of the NASA F-18 HARV

Extensive pressure measurements and off-surface flow visualization were obtained on the forebody and strakes of the NASA F-18 High Alpha Research Vehicle (HARV) equipped with actuated forebody strakes. Forebody yawing moments were obtained by integrating the circumferential pressures on the forebody and strakes. Results show that large yawing moments can be generated with forebody strakes. At angles of attack greater than 40 deg., deflecting one strake at a time resulted in a forebody yawing moment control reversal for small strake deflection angles. At alpha = 40 deg. and 50 deg., deflecting the strakes differentially about a 20 deg. symmetric strake deployment eliminated the control reversal and produced a near linear variation of forebody yawing moment with differential strake deflection. At alpha = 50 deg. and for 0 deg. and 20 deg. symmetric strake deployments, a larger forebody yawing moment was generated by the forward fuselage (between the radome and the apex of the leading-edge extensions), than on the radome where the actuated forebody strakes were located. Cutouts on the flight vehicle strakes that were not on the wind tunnel models are believed to be responsible for deficits in the suction peaks on the flight radome pressure distributions and differences in the forebody yawing moments.

Fisher, David F.↗

Forebody Aerodynamics of the F-18 High Alpha Research Vehicle with Actuated Forebody Strakes

Extensive pressure measurements and off-surface flow visualization were obtained on the forebody and strakes of the NASA F-18 High Alpha Research Vehicle (HARV) equipped with actuated forebody strakes. Forebody yawing moments were obtained by integrating the circumferential pressures on the forebody and strakes. Results show that large yawing moments can be generated with forebody strakes. At a 50 deg-angle-of-attack, deflecting one strake at a time resulted in a forebody yawing moment control reversal for small strake deflection angles. However, deflecting the strakes differentially about a 20 deg symmetric strake deployment eliminated the control reversal and produced a near linear variation of forebody yawing moment with differential strake deflection. At an angle of attack of 50 deg and for 0 deg and 20 deg symmetric strake deployments, a larger forebody yawing moment was generated by the forward fuselage (between the radome and the apex of the leading-edge extensions) than on the radome where the actuated forebody strakes were located. Cutouts on the flight vehicle strakes that were not on the wind tunnel models are believed to be responsible for deficits in the suction peaks on the flight radome pressure distributions and differences in the forebody yawing moments.

Fisher, David F.↗

Ion convection velocity reversals in the dayside cleft

Nearly simultaneous measurements of energetic particle precipitation and vector ion velocity have been made for the first time on Atmosphere Explorer C. These data have revealed the existence of both shear and rotational reversals in the ion convection velocity within the region of dayside cleft particle precipitation. The ion velocity data have shown that substantial regions of the dayside polar cap boundary are nearly electric equipotentials at which cleft particle precipitation is also observed but which cannot be identified with a merging region producing open field lines. Asymmetries in the polar cap flow, first implied by Heppner, may be determined by the extent and location of the equipotential boundaries.

Heelis, R. A.↗

Experimental investigations of the time and flow-direction responses of shear-stress-sensitive liquid crystal coatings

Time and flow-direction responses of shearstress-sensitive liquid crystal coatings were explored experimentally. For the time-response experiments, coatings were exposed to transient, compressible flows created during the startup and off-design operation of an injector-driven supersonic wind tunnel. Flow transients were visualized with a focusing Schlieren system and recorded with a 1000 frame/sec color video camera. Liquid crystal responses to these changing-shear environments were then recorded with the same video system, documenting color-play response times equal to, or faster than, the time interval between sequential frames (i.e., 1 millisecond). For the flow-direction experiments, a planar test surface was exposed to equal-magnitude and known-direction surface shear stresses generated by both normal and tangential subsonic jet-impingement flows. Under shear, the sense of the angular displacement of the liquid crystal dispersed (reflected) spectrum was found to be a function of the instantaneous direction of the applied shear. This technique thus renders dynamic flow reversals or flow divergences visible over entire test surfaces at image recording rates up to 1 KHz. Extensions of the technique to visualize relatively small changes in surface shear stress direction appear feasible.

Reda, Daniel C.↗

A hydrodynamic model of an outer hair cell

On the model it is possible to measure the force and the force direction for each individual hair as a function of the flow direction and velocity. Measurements were made at the man flow velocity .01 m/s, which is equivalent to a flow velocity in the real ear of about 1 micrometer/s. The kinematic viscosity of the liquid used in the model was 10,000 times higher than the viscosity of perilymph to attain hydrodynamic equality. Two different geometries for the sterocilia pattern were tested. First the force distribution for a W-shaped sterocilia pattern was recorded. This is the sterocilia pattern found in all real ears. It is found that the forces acting on the hairs are very regular and perpendicular to the legs of the W when the flow is directed from the outside of the W. When the flow is reversed, the forces are not reversed, but are much more irregular. This can eventually explain the half wave rectification of the nerve signals. As a second experiment, the force distribution for a V-shaped sterocilia pattern was recorded. Here the forces were irregular both when the flow was directed into the V and when it was directed against the edge of the V.

Jacobson, B. O.↗

Rapid Structured Volume Grid Smoothing and Adaption Technique

A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reduction in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.

Alter, Stephen J.↗