Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity steps”

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 685 records · Page 38

Microscopic dynamics of stress relaxation in a nanocolloidal soft glass

Here, we report x-ray photon correlation spectroscopy (XPCS) experiments with in situ rheometry performed on a soft glass composed of a concentrated suspension of charged silica nanoparticles subjected to step strains that induce yielding and flow. The XPCS measurements characterize the particle-scale and mesoscale motions within the glass that underlie the highly protracted decay of the macroscopic stress following the step strains. These dynamics are anisotropic, with slow, convective particle motion along the direction of the preceding shear that persists for surprisingly large times and that is accompanied by intermittent motion in the perpendicular (vorticity) direction. A close correspondence between the convective dynamics and stress relaxation is demonstrated by power-law scaling between the characteristic velocity of the collective particle motion and the rate of stress decay.

36 MATERIALS SCIENCE↗

Catalytic Reduction of Residual Oxygen for the Purification of Carbon Dioxide Streams from Pressurized Oxy-Combustors

Pressurized oxy-combustion (OC) is regarded as an energy-efficient alternative to atmospheric OC. To utilize the captured carbon dioxide (CO2) from pressurized OC for enhanced oil recovery or carbon storage in geological formations, the CO2 product stream must meet the required purity specifications. One major impurity in OC flue gases is residual O2 (e.g., 3 vol%). In this study, a catalytic approach has been investigated to remove the residual O2 through a reduction reaction with injected methane (CH4) in a reverse-flow fixed-bed (RFFB) reactor. The performance was tested for a supported copper (Cu)-based catalyst synthesized by incipient wetness impregnation and a cobalt-manganese (CoMn) oxide catalyst synthesized by a coprecipitation method. A commercial Cu-based catalyst was also studied as the reference. Experimental results revealed that under the reverse-flow operation, the temperature profiles in the reactor were symmetrical on both sides of the catalyst zone and were more homogenous than those under the conventional one-direction-flow operation. Among the tested catalysts, the CoMn catalyst was the most reactive for promoting the reduction of O2 with CH4, able to attain >95% O2 removal at 475 °C and a gas hourly space velocity (GHSV) of 18,000 h-1 with a feed gas containing 3 vol% O2. Under all experiments, the concentrations of carbon monoxide (CO) were lower than 30 ppmv in the effluent gas streams. The reaction between O2 and CH4 could be depicted by using the Mars–van Krevelen redox mechanism based on a two-step reduction–oxidation reaction pathway.

36 MATERIALS SCIENCE↗

Scanning thermal plumes

Over a three-year period 800 thermal line scans of power plant plumes were made by an airborne scanner, with ground truth measured concurrently at the plants. Computations using centered finite differences in the thermal scanning imagery show a lower bound in the horizontal temperature gradient in excess of 1.6 C/m. Gradients persist to 3 m below the surface. Vector plots of the velocity of thermal fronts are constructed by tracing the front motion in successive thermal images. A procedure is outlined for the two-point ground calibration of a thermal scanner from an equation describing the scanner signal and the voltage for two known temperatures. The modulation transfer function is then calculated by input of a thermal step function and application of digital time analysis techniques using Fast Fourier Transforms to the voltage output. Field calibration tests are discussed. Data accuracy is limited by the level of ground truth effort chosen.

Scarpace, F. L.↗

Mean Winds and Tides over Urbana, Illinois During November 16-22, 1981

The Urbana Meteor Radar uses a high power pulse transmitter to achieve high echo rates, and an interferometer to spatially locate the individual echoes. The irregularly spaced observations of range, radial velocity, and direction cosines are reduced to a rectangular grid, evenly spaced in time and height, of estimates of the zonal and meridional wind. This is accomplished using the method of GROVES whereby a constant horizontal wind is fitted, using least squares, to all echoes occurring within a fixed time/height interval. This interpolation region is then stepped in time and height. the size of the interpolation region is 2 hours by 5 km, the time step is approximately 11 minutes, and the height step is 2 km. The filtering of high temporal and spatial frequencies is somewhat reduced by weighting the echoes according to their distance from the center of the interpolation box.

Tetenbaum, D.↗

Prediction of time-dependent fuselage pressures in the wake of a helicopter rotor

A method is described for the analysis of the unsteady, incompressible potential flow associated with a helicopter. The fuselage, rotor, and rotor wake are considered together in forward flight. This method is particularly useful in low advance ratio flight due to the contribution of the deformed wake to the airloads imposed on the fuselage. The rotor geometry is prescribed and the unsteady wake geometry is computed from the local flow perturbation velocities. The fuselage is modeled as a non-lifting body of source panels and the rotor and its wake are modeled as a full vortex lattice. The rotor geometry is arbitrary and several rotor blades can be represented. The unsteady airloads are computed in the presence of the deformed rotor wake by a time-stepping technique. The wake is started impulsively from rest, allowing a natural convection of the wake with time.

Berry, John D.↗

An innovations approach to decoupling of multibody dynamics and control

The problem of hinged multibody dynamics is solved using an extension of the innovations approach of linear filtering and prediction theory to the problem of mechanical system modeling and control. This approach has been used quite effectively to diagonalize the equations for filtering and prediction for linear state space systems. It has similar advantages in the study of dynamics and control of multibody systems. The innovations approach advanced here consists of expressing the equations of motion in terms of two closely related processes: (1) the innovations process e, a sequence of moments, obtained from the applied moments T by means of a spatially recursive Kalman filter that goes from the tip of the manipulator to its base; (2) a residual process, a sequence of velocities, obtained from the joint-angle velocities by means of an outward smoothing operations. The innovations e and the applied moments T are related by means of the relationships e = (I - L)T and T = (I + K)e. The operation (I - L) is a causal lower triangular matrix which is generated by a spatially recursive Kalman filter and the corresponding discrete-step Riccati equation. Hence, the innovations and the applied moments can be obtained from each other by means of a causal operation which is itself casually invertible.

Rodriguez, G.↗

Effect of Reynolds number and turbulence on airfoil aerodynamics at -90-degree incidence

A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered-grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the noniterative solution of the flowfield and satisfies the continuity equation to machine zero at each time step. The method is evaluated in terms of its stability to predict two-dimensional flow about an airfoil at -90-deg incidence for varying Reynolds number and laminar/turbulent models. The variations of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisom indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Stremel, Paul M.↗

Pulse Detonation Engine for Advanced Oxy-Combustion of Coal-Based Fuel for Direct Power Extraction Applications

Pressure gain combustion and magnetohydrodynamic (MHD) systems have the potential to provide a step increase in the efficiency of combined-cycle power plants. Specific advantages include a net pressure increase to the system instead of a pressure drop, the high temperature of the detonation waves can increase the efficiency of power extraction compared to other processes, significant thermal energy can be released in a compact region, and the high velocities of the flow increase extraction of electrical power. In summary, a pressure gain combustor coupled with a MHD has the potential to be transformative. Despite the potential advantages, relatively little research has been conducted considering coupled pressure gain combustion systems with MHD systems. With this background and motivation, the overall goal of this effort was to advance the knowledge, technology, and computational tools associated with coupled detonation and MHD systems. A joint experimental and computational approach was used while seeking to accomplish the goals of this work. Specifically, two pulse detonation engines were developed and used for the experiments to produce detonations. Detonation speeds were measured for a variety of flow and fuel conditions (e.g., methane, propane, with coal particles). Preliminary electrical conductivity measurements were collected. An extensive amount of research was performed to identify sensitivities of detonation behaviors to the presence of combustion products. Computationally, a twofold approach was used in this work. First, a solver was developed for solving the governing equations for a reactive flow with coupled detonation and ionization chemistry. The solver was applied to study the impacts of seed material ionization on detonation. Second, a conservation element-solution element (CE-SE) based numerical solver for detonation studies with a reduced reaction mechanism for oxy-methane combustion was developed and verified on standard test cases. Key findings and contributions from this work are as follows. A system was developed for injecting powderized coal, or other seeding material, into a pulse-detonation system. The influence of a combustion product (i.e., CO 2 ) on detonation behavior was identified. Knowledge gained from this work is applicable to devices such as rotating detonation engines, where combustion products mix with fresh reactants. A system for measuring the electrical conductivity of the exhaust from a pulse-detonation engine was developed. The open-source solver, Clawpack, was extended to solve the reactive Euler equations for simulating detonations. A coupled combustion and ionization chemistry was developed in a single chemical kinetic model for methane oxidation. This model can be used to solve coupled MHD and detonation simulations. It was found that parasitic interactions from ionization chemistry with the magnetic field can reduce the detonation velocity by up to 8%, with a potential impact on power extraction of 15%. It is recommended that interactions between the detonation front and MHD field be considered. Higher gas temperatures and velocities were achieved owing to oxy-fuel detonations. Use of radical dissociation reactions in the reduced reaction mechanism, was found to be critical in predicting detonation temperature and velocity accurately.

01 COAL, LIGNITE, AND PEAT↗

Pioneer Venus Orbiter planar retarding potential analyzer plasma experiment

The retarding potential analyzer (RPA) on the Pioneer Venus Orbiter Mission measures most of the thermal plasma parameters within and near the Venusian ionosphere. Parameters include total ion concentration, concentrations of the more abundant ions, ion temperatures, ion drift velocity, electron temperature, and low-energy (0-50 eV) electron distribution function. Several functions not previously used in RPA's were developed and incorporated into this instrument to accomplish these measurements on a spinning spacecraft with a small bit rate. The more significant functions include automatic electrometer ranging with background current compensation; digital, quadratic retarding potential step generation for the ion and low-energy electron scans; a current sampling interval of 2 ms throughout all scans; digital logic inflection point detection and data selection; and automatic ram direction detection.

Knudsen, W. C.↗

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗

Sensitivity analysis for the dynamic response of viscoplastic shells of revolution

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic axisymmetric response of viscoplastic shells of revolution. The analytical formulation is based on Reissner's large deformation shell theory with the effects of transverse shear deformation, rotatory inertia and moments turning around the normal to the middle surface included. The material model is chosen to be isothermal viscoplasticity, and an associated flow rule is used with a von Mises effective stress. A mixed formulation is used with the fundamental unknowns consisting of six stress resultants, three generalized displacements and three velocity components. Spatial discretization is performed using finite elements, with discontinuous stress resultants across element interfaces. The temporal integration is performed by using an explicit central difference scheme (leap-frog method) with an implicit constitutive update. The sensitivity coefficients are evaluated using a direct differentiation approach. Numerical results are presented for a spherical cap subjected to step loading, and a circular plate subjected to impulsive loading. The sensitivity coefficients are generated by evaluating the derivatives of the response quantities with respect to thickness, mass density, Young's modulus, and two of the material parameters characterizing the viscoplastic response. Time histories of the response and sensitivity coefficients are presented, along with spatial distributions of these quantities at selected times.

Kulkarni, Makarand↗

Weightless bubble lattices: A case of froth wicking

In the absence of gravity drainage, froth wicking draws excess fluid onto a bubble lattice. Capillary forces only cause fluid transport; a moving front moves stably and without fluid fingering along a constant velocity bubble-fluid contact line. This percolation of fluid crawling up the lattice shows fluid coverage on lattice borders varies linearly with available surface area (proportional to lattice perimeter) and fluid accelerates through regions or nests of high bubble density (number of bubbles/sq cm). The development of nearly two-dimensional bubble lattices in variable gravity (step function between 0.01 and 1.8 times earthly gravity) are examined experimentally and a zeroth-order model for froth wetting is presented, which captures many of the principal observations. Possible applications for bubble lattices include adhesion casting of metals and separation of biological cells, bacteria, and particles.

Noever, David A.↗

Update on ASME AG-1 and Revision of the DOE Nuclear Air Cleaning Handbook - 20536

The American Society of Mechanical Engineers AG-1 Code on Nuclear Air and Gas Treatment is updated every two years. This paper provides an overview of issues being discussed within AG-1 committee meetings and revisions that will be included in the next version of the code. Additional information is provided on issues being discussed ranging from media velocity limits, control of airflow using actual versus and standard cubic feet per minute. Navigation and identification of applicable high-efficiency particulate air (HEPA) filter code sections can be confusing for those unfamiliar with AG-1. A general reorganization of filtration code sections clarifies geometry, flow, and performance envelopes for nuclear grade HEPA filters. The initial step consists of establishment of Section FN covering qualification of conventional, ceramic, metal, and high strength fibrous glass media. Future filter sections are currently expected to be characterized by geometry and flow direction (axial vs. radial). The Department of Energy's (DOE) Nuclear Air Cleaning Handbook (NACH) is also undergoing revision. An update on the status of that effort, along with issues and changes being considered. The original charter of the AG-1 Code was specific to commercial nuclear power plants, and to safety-related air cleaning systems. This has evolved over the years to include normal non-safety ventilation systems and process gas systems. Corollary guidance for DOE facilities has been available in several editions of the Nuclear Air Cleaning Handbook, starting with ORNL-NSIC 65. The latest edition of the Handbook, 2003, is undergoing revision and modifications. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The effect of shearing strain-rate on the ultimate shearing resistance of clay

An approach for investigating the shearing resistance of cohesive soils subjected to a high rate of shearing strain is described. A fast step-loading torque apparatus was used to induce a state of pure shear in a hollow cylindrical soil specimen. The relationship between shearing resistance and rate of shear deformation was established for various soil densities expressed in terms of initial void ratio or water content. For rate of shearing deformation studies, the shearing resistance increases initially with shearing velocity, but subsequently reaches a terminal value as the shearing velocity increases. The terminal shearing resistance is also found to increase as the density of the soil increases. The results of this investigation are useful in the rheological study of clay. It is particularly important for mobility problems of soil runways, since the soil resistance is found to be sensitive to the rate of shearing.

Cheng, R. Y. K.↗

An Empirical Method Permitting Rapid Determination of the Area, Rate and Distribution of Water-Drop Impingement on an Airfoil of Arbitrary Section at Subsonic Speeds

An empirical method for the determination of the area, rate, and distribution of water-drop impingement on airfoils of arbitrary section is presented. The procedure represents an initial step toward the development of a method which is generally applicable in the design of thermal ice-prevention equipment for airplane wing and tail surfaces. Results given by the proposed empirical method are expected to be sufficiently accurate for the purpose of heated-wing design, and can be obtained from a few numerical computations once the velocity distribution over the airfoil has been determined. The empirical method presented for incompressible flow is based on results of extensive water-drop. trajectory computations for five airfoil cases which consisted of 15-percent-thick airfoils encompassing a moderate lift-coefficient range. The differential equations pertaining to the paths of the drops were solved by a differential analyzer. The method developed for incompressible flow is extended to the calculation of area and rate of impingement on straight wings in subsonic compressible flow to indicate the probable effects of compressibility for airfoils at low subsonic Mach numbers.

Bergrun, N. R.↗

The effect of Reynolds number and turbulence on airfoil aerodynamics at -90 degrees incidence

A method has been developed for calculating the viscous flow about airfoils in with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the direct solution of the flow field and satisfies the continuity equation to machine zero at each time-step. The method is evaluated in terms of its ability to predict two-dimensional flow about an airfoil at -90 degrees incidence for varying Reynolds number and different boundary layer models. A laminar and a turbulent boundary layer models. A laminar and a turbulent boundary layer model are considered in the evaluation of the method. The variation of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisons indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Stremel, Paul M.↗

Real-Time Aerodynamic Flow and Data Visualization in an Interactive Virtual Environment

Significant advances have been made to non-intrusive flow field diagnostics in the past decade. Camera based techniques are now capable of determining physical qualities such as surface deformation, surface pressure and temperature, flow velocities, and molecular species concentration. In each case, extracting the pertinent information from the large volume of acquired data requires powerful and efficient data visualization tools. The additional requirement for real time visualization is fueled by an increased emphasis on minimizing test time in expensive facilities. This paper will address a capability titled LiveView3D, which is the first step in the development phase of an in depth, real time data visualization and analysis tool for use in aerospace testing facilities.

Schwartz, Richard J.↗

TOWARD A METHOD FOR SCALING HUMAN BODY MODELS IN AN IMU-BASED WORKFLOW

BACKGROUND Scaled biomechanical models can more accurately inform crew health decisions when tailored to the wide range of astronaut sizes. One component to improve scaling of existing models to better represent each unique astronaut’s size is the individual length scaling of limbs. Traditionally, these lengths are determined by motion capture or manual measurement. A new method is herein proposed for length scaling which can be done by measuring linear and angular accelerations at a desired point during isolated motion around a point of rotation, then calculating the distance between the desired point and point of rotation. When an Inertial Measurement Unit (IMU) device is placed at the distal point of a limb, the isolated motion is about that limb’s proximal joint. For example, to measure forearm length, an IMU is placed at the wrist, the point of rotation is at the elbow, and the isolated motion is forearm flexion and extension. These calculated lengths are then used to scale models to each unique astronaut’s size, thereby improving the applicability of the model. This method of scaling limb segments can be used for any limb that has an easily defined proximal joint for the limb to rotate around including hands, arms, legs, feet. Utilizing IMUs for data collection also provides the synergistic ability to record data without a dedicated space in a room with many cameras, therefore reducing the data collection footprint, or record data where optical motion capture is not possible, such as inside a spacesuit. METHODS AND RESULTS To test this method, upper body data collection was performed with 5 Xsens DOT IMUs on a single subject. IMUs consist of an accelerometer, a gyroscope, and a magnetometer which collect linear acceleration, angular velocity, and magnetic fluctuations, respectively. Before any ground-based laboratory collection, the magnetic fluctuations are used to correct the heading of the IMU in space relative to the Earth’s magnetic field. The direct measurement of angular velocity is integrated to calculate angular acceleration. Then the linear acceleration ( a ) and angular acceleration (α) are solved using r = at/α to calculate the radius, which in this case is the distance between the IMU and the point of rotation (i.e., segment length). The IMU must be placed at the most distal point of the limb being measured (i.e., ankle if measuring lower leg length) and the test plan must consist of an isolated motion about that limb’s proximal joint (i.e., knee flexion and extension if measuring lower leg length). The distances (radii) calculated at every time interval were filtered (bandpass filter keeping 5th-90th percentile data) to eliminate outliers and spurious data that occur when the isolated motion was stopped or nearly stopped. The remaining distances were averaged, resulting in the calculated limb length. Scaling factors were then computed by dividing the calculated limb length by the unscaled model’s length. These scale factors are plugged into the Scale Tool in OpenSim [1,2] to apply the scaling to the OpenSim Full Body Rajagopal Model [3,4]. Manual measurements of limb lengths were taken before data collection started and used for comparing against the calculated lengths. The Anthropometric Survey of US Army Personnel (ANSUR II) [5] was also used as a third source of reference for limb length measurements. The following measurements were retrieved from the subject before data collection: 34.5 cm from L1 to C7 (thorax), 25.7 cm from C7 Joint Center (JC) to head vertex (neck and head), 36.3 cm from shoulder JC to elbow JC (humerus), 29.5 cm from elbow JC to wrist JC (forearm), and 16.2 cm from clavicle to acromion (clavicle). Of those five, forearm and humerus lengths were calculated using this proposed method to obtain preliminary results. The forearm length after filtering and averaging was calculated to be 37.6 cm. This is a 28% overestimation from the measured forearm length (29.5 cm). The humerus length after filtering and averaging was calculated to be 47.8 cm. This is a 31% difference from the measured subject length (36.3 cm). Sources of error include imperfect isolated motion (method currently expects that motion should be perfectly circular in a 2D plane, include no rotation of the IMU, and be relatively smooth; a more secure IMU attachment method will help), unrefined filter techniques (removed highest and lowest values with 20% high and low pass filters and no smoothing filters), arbitrary removal of stopped or near stopped data (kept data for only a short range before and after the angular velocity peaking), and a more representative method for removing gravitational acceleration is needed (current method is to zero all accelerations against a baseline taken just before the isolated motion which does not account for the gravitational acceleration changed due to IMU rotation during movement). Addressing these error sources will improve the accuracy of the limb length calculation. Next steps include creating a method for whole-body scaling estimation using individual limb scale factors. Continued pursuit of these techniques is expected to enable acquiring anthropometric information using only IMUs in real-time.

E. K. Marecki↗