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 289 records · Page 16

Development, characterization, and modeling of a high-performance Ru/B2CA catalyst for ammonia synthesis

This paper documents the development and performance of a nano-phase Ru catalyst on a (BaO) x (CaO) y (Al 2 O 3 ) support. Extensive screening of the support’s ternary composition shows the best stoichiometry is (BaO) 2 (CaO)(Al 2 O 3 ), denoted B2CA. The paper first describes catalyst preparation and characterization. The paper reports a detailed 12-step reaction mechanism that represents ammonia synthesis over wide ranges of temperature, pressure, space velocity, and feed composition. Additionally, the mechanism is developed and validated using results of packed-bed experiments. The elementary reaction pathways consider surface adsorbates, including catalyst-poisoning behaviors. The rate expressions include important coverage-dependent activation barriers. Machine learning models assist interpretation of the catalyst-support interactions. The detailed chemistry is much more predictive than is possible with global representations (N 2 + 3H 2 ⇌ 2NH 3 ). The validated models can be applied to assist optimizing reactor design and operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formation of chromospheric resonance line profiles in supergiants

The formation of chromospheric resonance line profiles is examined for the case of the relatively low-density atmospheres appropriate to late-type supergiants. The effects of partial frequency redistribution control the emergent line profiles, to the extent that even an isothermal atmosphere can give rise to apparent emission features. Schematic model chromospheres are used to demonstrate the effects of different velocity fields and temperature-density structures on the line profiles, as a first step toward a clearer understanding of the chromospheric line profiles and why supergiants have broad emission lines (the Wilson-Bappu effect). Neither the Doppler nor the mass column density explanation alone can explain the Wilson-Bappu effect in supergiants, but both theories may play a role.

Basri, G. S.↗

Computer coordination of limb motion for locomotion of a multiple-armed robot for space assembly

Consideration is given to a possible robotic system for the construction of large space structures, which may be described as a multiple general purpose arm manipulator vehicle that can walk over the structure under construction to a given site for further work. A description is presented of the locomotion of such a vehicle, modeling its arms in terms of a currently available industrial manipulator. It is noted that for whatever maximum speed of operation is chosen, rapid changes in robot velocity create situations in which already-selected handholds are no longer practical. A step is added to the 'free gait' walking algorithm in order to solve this problem.

Klein, C. A.↗

Rapid acceleration of energetic particles in the 1982 February 8 solar flare

Hard X-ray and gamma-ray measurements of the February 8, 1982 (1250 UT) solar flare made with spectrometers aboard the ISEE 3 and SMM spacecraft show that bursts of photons from 40 keV to 40 MeV are coincident within + or - 1 s. This indicates that the acceleration of particles to relativistic velocities can occur promptly (within 1 s). As far as the energetic particles interacting at the sun are concerned, the 'two-phase' process, which requires several minutes for particle acceleration to relativistic velocities, cannot account for the present gamma-ray observations. Even the 'two-step' acceleration process, which is assumed to occur in a small magnetic loop, is severely constrained by the present observations.

Kane, S. R.↗

Universal limiter for transient interpolation modeling of the advective transport equations: The ULTIMATE conservative difference scheme

A fresh approach is taken to the embarrassingly difficult problem of adequately modeling simple pure advection. An explicit conservative control-volume formation makes use of a universal limiter for transient interpolation modeling of the advective transport equations. This ULTIMATE conservative difference scheme is applied to unsteady, one-dimensional scalar pure advection at constant velocity, using three critical test profiles: an isolated sine-squared wave, a discontinuous step, and a semi-ellipse. The goal, of course, is to devise a single robust scheme which achieves sharp monotonic resolution of the step without corrupting the other profiles. The semi-ellipse is particularly challenging because of its combination of sudden and gradual changes in gradient. The ULTIMATE strategy can be applied to explicit conservation schemes of any order of accuracy. Second-order schemes are unsatisfactory, showing steepening and clipping typical of currently popular so-called high resolution shock-capturing of TVD schemes. The ULTIMATE third-order upwind scheme is highly satisfactory for most flows of practical importance. Higher order methods give predictably better step resolution, although even-order schemes generate a (monotonic) waviness in the difficult semi-ellipse simulation. Little is to be gained above ULTIMATE fifth-order upwinding which gives results close to the ultimate for which one might hope.

Leonard, B. P.↗

Pressure algorithm for elliptic flow calculations with the PDF method

An algorithm to determine the mean pressure field for elliptic flow calculations with the probability density function (PDF) method is developed and applied. The PDF method is a most promising approach for the computation of turbulent reacting flows. Previous computations of elliptic flows with the method were in conjunction with conventional finite volume based calculations that provided the mean pressure field. The algorithm developed and described here permits the mean pressure field to be determined within the PDF calculations. The PDF method incorporating the pressure algorithm is applied to the flow past a backward-facing step. The results are in good agreement with data for the reattachment length, mean velocities, and turbulence quantities including triple correlations.

Anand, M. S.↗

An assessment of full potential and Euler solutions for self-generated rotor blade-vortex interactions

An assessment is made of the accuracy and suitability of the three-dimensional full potential and Euler equations in the modeling of subcritical and supercritical rotor blade-vortex interactions (BVI). The interaction velocity fields are obtained through a nonlinear superposition of the rotor and the vortex wake flow fields. Vortex effects are simulated using the velocity 'transpiration' approach. A modified surface condition is prescribed and enforced at each time step of the computations to satisfy the tangency boundary condition. Potential blade-vortex encounters are identified and tracked in time at equal increments of rotor azimuth using the lifting-line helicopter/rotor trim code CAMRAD. This information is then utilized in interpolation routines within the flow solvers to compute the instantaneous position(s) of the interaction vortex element(s) with respect to the blade for the time-accurate calculations. Results for subcritical self-generated BVI and for supercritical interactions resulting from an upstream-generated vortex are presented.

Hassan, A. A.↗

MSFC Stream Model Preliminary Results: Modeling Recent Leonid and Perseid Encounters

The cometary meteoroid ejection model of Jones and Brown (1996b) was used to simulate ejection from comets 55P/Tempel-Tuttle during the last 12 revolutions, and the last 9 apparitions of 109P/Swift-Tuttle. Using cometary ephemerides generated by the Jet Propulsion Laboratory s (JPL) HORIZONS Solar System Data and Ephemeris Computation Service, two independent ejection schemes were simulated. In the first case, ejection was simulated in 1 hour time steps along the comet s orbit while it was within 2.5 AU of the Sun. In the second case, ejection was simulated to occur at the hour the comet reached perihelion. A 4th order variable step-size Runge-Kutta integrator was then used to integrate meteoroid position and velocity forward in time, accounting for the effects of radiation pressure, Poynting-Robertson drag, and the gravitational forces of the planets, which were computed using JPL s DE406 planetary ephemerides. An impact parameter was computed for each particle approaching the Earth to create a flux profile, and the results compared to observations of the 1998 and 1999 Leonid showers, and the 1993 and 2004 Perseids.

Cooke, William J.↗

Weight Assessment for Fuselage Shielding on Aircraft With Open-Rotor Engines and Composite Blade Loss

The Federal Aviation Administration (FAA) has been engaged in discussions with airframe and engine manufacturers concerning regulations that would apply to new technology fuel efficient "openrotor" engines. Existing regulations for the engines and airframe did not envision features of these engines that include eliminating the fan blade containment systems and including two rows of counter-rotating blades. Damage to the airframe from a failed blade could potentially be catastrophic. Therefore the feasibility of using aircraft fuselage shielding was investigated. In order to establish the feasibility of this shielding, a study was conducted to provide an estimate for the fuselage shielding weight required to provide protection from an open-rotor blade loss. This estimate was generated using a two-step procedure. First, a trajectory analysis was performed to determine the blade orientation and velocity at the point of impact with the fuselage. The trajectory analysis also showed that a blade dispersion angle of 3deg bounded the probable dispersion pattern and so was used for the weight estimate. Next, a finite element impact analysis was performed to determine the required shielding thickness to prevent fuselage penetration. The impact analysis was conducted using an FAA-provided composite blade geometry. The fuselage geometry was based on a medium-sized passenger composite airframe. In the analysis, both the blade and fuselage were assumed to be constructed from a T700S/PR520 triaxially-braided composite architecture. Sufficient test data on T700S/PR520 is available to enable reliable analysis, and also demonstrate its good impact resistance properties. This system was also used in modeling the surrogate blade. The estimated additional weight required for fuselage shielding for a wing- mounted counterrotating open-rotor blade is 236 lb per aircraft. This estimate is based on the shielding material serving the dual use of shielding and fuselage structure. If the shielding material is not used for dual purpose, and is only used for shielding, then the additional weight per aircraft is estimated to be 428 lb. This weight estimate is based upon a number of assumptions that would need to be revised when applying this concept to an actual airplane design. For example, the weight savings that will result when there is no fan blade containment system, manufacturing limitations which may increase the weight where variable thicknesses was assumed, engine placement on the wing versus aft fuselage, etc.

aircraft safety↗

Transient boiling heat transfer in saturated liquid nitrogen and F113 at standard and zero gravity

Transient and steady state nucleate boiling in saturated LN2 and F113 at standard and near zero gravity conditions were investigated for the horizontal up, vertical and horizontal down orientations of the heating surface. Two distinct regimes of heat transfer mechanisms were observed during the interval from the step increase of power input to the onset of nucleate boiling: the conduction and convection dominated regimes. The time duration in each regime was considerably shorter with LN2 than with F113, and decreased as heat flux increased, as gravity was reduced, and as the orientation was changed from horizontal up to horizontal down. In transient boiling, boiling initiates at a single point following the step increase in power, and then spreads over the surface. The delay time for the inception of boiling at the first site, and the velocity of spread of boiling varies depending upon the heat flux, orientation, body force, surface roughness and liquid properties, and are a consequence of changes in boundary layer temperature levels associated with changes in natural convection. Following the step increase in power input, surface temperature overshoot and undershoot occur before the steady state boiling temperature level is established.

Oker, E.↗

Cavitating Acoustics Modeling in Sierra/SD

Motivated by the need to simulate the effects of underwater explosion on ship structures, we develop a new cavitating acoustics formulation. The proposed approach is consistent with existing methods where the cavitation phenomenon is captured with a bilinear constitutive law. However, the new formulation is in terms of velocity potential, as opposed to the existing displacement-potential and pressure formulations. Also unique to the proposed formulation is a new generalized time-stepping procedure specific to cavitating acoustics, which has the ability to introduce numerical damping to control frothing. Numerical examples of varying complexity are presented to illustrate the effectiveness of the proposed approach and the ability to use velocity potential as a primary field variable for cavitating acoustics simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗

Application of unsteady airfoil theory to rotary wings

A clarification is presented on recent work concerning the application of unsteady airfoil theory to rotary wings. The application of this theory may be seen as consisting of four steps: (1) the selection of an appropriate unsteady airfoil theory; (2) the resolution of that velocity which is the resultant of aerodynamic and dynamic velocities at a point on the elastic axis into radial, tangential and perpendicular components, and the angular velocity of a blade section about the deformed axis; (3) the expression of lift and pitching moments in terms of the three components; and (4) the derivation of explicit expressions for the components in terms of flight velocity, induced flow, rotor rotational speed, blade motion variables, etc.

Kaza, K. R. V.↗

Flowfield measurements in a separated and reattached flat plate turbulent boundary layer

The separation and reattachment of a large-scale, two-dimensional turbulent boundary layer at low subsonic speed on a flat plate has been studied experimentally. The separation bubble was 55 cm long and had a maximum bubble thickness, measured to the height of the mean dividing streamline, of 17 cm, which was twice the thickness of the inlet boundary layer. A combination of laser velocimetry, hot-wire anemometry, pneumatic probing techniques, and flow visualization were used as diagnostics. Principal findings were that an outer inviscid rotational flow was defined which essentially convected over the blockage associated with the inner, viscously dominated bubble recirculation region. A strong backflow region in which the flow moved upstream 100 percent of the time was measured near the test surface over the central 35 percent of the bubble. A laminar backflow boundary layer having pseudo-turbulent characteristics including a log-linear velocity profile was generated under the highly turbulent backflow. Velocity profile shapes in the reversed flow region matched a previously developed universal backflow profile at the upstream edge of the separation region but not in the steady backflow region downstream. A smoke flow visualization movie and hot-film measurements revealed low frequency nonperiodic flapping at reattachment. However, forward flow fraction data at reattachment and mean velocity profiles in the redeveloping boundary layer downstream of reattachment correlated with backward-facing step data when the axial dimension was scaled by the distance from the maximum bubble thickness to reattachment.

Patrick, William P.↗

A test of a vortex method for the computation of flap side edge noise

Upon approach to landing, a major source location of airframe noise occurs at the side edges of the part span, trailing edge flaps. In the vicinity of these flaps, a complex arrangement of spanwise flow with primary and secondary tip vortices may form. Each of these vortices is observed to become fully three-dimensional. In the present study, a numerical model is developed to investigate the noise radiated from the side edge of a flap. The inherent three-dimensionality of this flow forces us to carefully consider a numerical scheme which will be both accurate in its prediction of the flow acoustics and also computationally efficient. Vortex methods have offered a fast and efficient means of simulating many two and three-dimensional, vortex dominated flows. In vortex methods, the time development of the flow is tracked by following exclusively the vorticity containing regions. Through the Biot-Savart law, knowledge of the vorticity field enables one to obtain flow quantities at any desired location during the flow evolution. In the present study, a numerical procedure has been developed which incorporates the Lagrangian approach of vortex methods into a calculation for the noise radiated by a flow-surface interaction. In particular, the noise generated by a vortex in the presence of a flat half plane is considered. This problem serves as a basic model of flap edge flow. It also permits the direct comparison between our computed results and previous acoustic analyses performed for this problem. In our numerical simulations, the mean flow is represented by the complex potential W(z) = Aiz(exp l/2), which is obtained through conformal mapping techniques. The magnitude of the mean flow is controlled by the parameter A. This mean flow has been used in the acoustic analysis by Hardin and is considered a reasonable model of the flow field in the vicinity of the edge and away from the leading and trailing edges of the flap. To represent the primary vortex which occurs near the flap, a point vortex is introduced just below the flat half plane. Using a technique from panel methods, boundary conditions on the flap surface are satisfied by the introduction of a row of stationary point vortices along the extent of the flap. At each time step in the calculation, the strength of these vortices is chosen to eliminate the normal velocity at intermediary collocation points. The time development of the overall flow field is then tracked using standard techniques from vortex methods. Vortex trajectories obtained through this computation are in good agreement with those predicted by the analytical solution given by Hardin, thus verifying the viability of this procedure for more complex flow arrangements. For the flow acoustics, the Ffowcs Williams-Hawkings equation is numerically integrated. This equation supplies the far field acoustic pressure based upon pressures occurring along the flap surface. With our vortex method solution, surface pressures may be obtained with exceptional resolution. The Ffowcs Williams-Hawkings equation is integrated using a spatially fourth order accurate Simpson's rule. Rational function interpolation is used to obtain the surface pressures at the appropriate retarded times. Comparisons between our numerical results for the acoustic pressure and those predicted by the Hardin analysis have been made. Preliminary results indicate the need for an improved integration technique. In the future, the numerical procedure developed in this study will be applied to the case of a rectangular flap of finite thickness and ultimately modified for application to the fully three-dimensional problem.

Martin, James E.↗

The spiral-compact galaxy pair AM 2208-251: Computer simulations versus observations

The system AM2208-251 is a roughly edge-on spiral extending east-west with a smaller round compact E system about 60 arcsec east of the spiral nucleus along the major axis of the spiral. Bertola, Huchtmeier, and Zeilinger (1990) have presented optical spectroscopic as well as single dish 21 cm observations of this system. Their spectroscopic data show, via emission lines lambda lambda 3727-29A, a rising rotation curve near the nucleus. These spectroscopic observations may indicate a tidal interaction in the system. In order to learn more about such pairs, the authors simulated the interaction using the computer model developed by Miller (1976 a,b, 1978) and modified by the authors (Byrd 1986, 1987, 1988). To do the simulation they need an idea of the mutual orbits of the two galaxies. Their computer model is a two-dimensional polar N-body program. It consists of a self-gravitating disk of particles, within an inert axially symmetric stabilizing halo potential. The particles are distributed in a 24(radial) by 36(azimuthal) polar grid. Self consistent calculations can be done only within the grid area. The disk is modeled with a finite Mestel disk, where all the particles initially move in circular orbits with constant tangential velocities (Mestel 1963), resulting in a flat rotation curve. The gas particles in the spiral's disk, which make up 30 percent of its mass, collide in the following manner. The number of particles in each bin of the polar grid is counted every time step. If it is greater than a given critical density, all the particles in the bin collide, obtaining in the result the same velocities, equal to the average for the bin. This process produces clumps of gas particles-the star formation sites. The authors suppress the collision in the inner part of the disk (within the circle r = 6) to represent the hole seen in the gas in the nuclear bulge of spirals. They thus avoid spurious effects due to collisions in that region. They also varied the size of the collisional bins, which did not affect their conclusions.

Klaric, Mario↗

Transverse injection into Mach 2 flow behind a rearward-facing step - A 3-D, compressible flow test case for hypersonic combustor CFD validation

A spatially-complete data set of the important primitive flow variables is presented for the complex, nonreacting, 3D unit combustor flow field employing transverse injection into a Mach 2 flow behind a rearward-facing step. A unique wind tunnel facility providing the capability for iodine seeding was built specifically for these measurements. Two optical techniques based on laser-induced-iodine fluorescence were developed and utilized for nonintrusive, in situ flow field measurements. LDA provided both mean and fluctuating velocity component measurements. A thermographic phosphor wall temperature measurement technique was developed and employed. Data from the 2D flow over a rearward-facing step and the complex 3D mixing flow with injection are reported.

Mcdaniel, James C.↗

Acceleration and Velocity Sensing from Measured Strain

A simple approach for computing acceleration and velocity of a structure from the strain is proposed in this study. First, deflection and slope of the structure are computed from the strain using a two-step theory. Frequencies of the structure are computed from the time histories of strain using a parameter estimation technique together with an autoregressive moving average model. From deflection, slope, and frequencies of the structure, acceleration and velocity of the structure can be obtained using the proposed approach. Simple harmonic motion is assumed for the acceleration computations, and the central difference equation with a linear autoregressive model is used for the computations of velocity. A cantilevered rectangular wing model is used to validate the simple approach. Quality of the computed deflection, acceleration, and velocity values are independent of the number of fibers. The central difference equation with a linear autoregressive model proposed in this study follows the target response with reasonable accuracy. Therefore, the handicap of the backward difference equation, phase shift, is successfully overcome.

shape sensing↗