Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Load 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 937 records · Page 52

Development, Demonstration, and Analysis of an Integrated Iodine Hall Thruster Feed System

The design of an in-space iodine-vapor-fed Hall effect thruster propellant management system is described. The solid-iodine propellant tank has unique issues associated with the microgravity environment, requiring a solution where the iodine is maintained in intimate thermal contact with the heated tank walls. The flow control valves required alterations from earlier iterations to survive for extended periods of time in the corrosive iodine-vapor environment. Materials have been selected for the entire feed system that can chemically resist the iodine vapor, with the design now featuring Hastelloy or Inconel for almost all the wetted components. An integrated iodine feed system/Hall thruster demonstration unit was fabricated and tested, with all control being handled by an onboard electronics card specifically designed to operate the feed system. Structural analysis shows that the feed system can survive launch loads after the implementation of some minor reinforcement. Flow modeling, while still requiring significant additional validation, is presented to show its potential in capturing the behavior of components in this low-flow, low-pressure system.

Polzin, Kurt A.↗

Fluid shear stress stimulates prostaglandin and nitric oxide release in bone marrow-derived preosteoclast-like cells

Bone is a porous tissue that is continuously perfused by interstitial fluid. Fluid flow, driven by both vascular pressure and mechanical loading, may generate significant shear stresses through the canaliculi as well as along the bone lining at the endosteal surface. Both osteoblasts and osteocytes produce signaling factors such as prostaglandins and nitric in response to fluid shear stress (FSS); however, these humoral agents appear to have more profound affects on osteoclast activity at the endosteal surface. We hypothesized that osteoclasts and preosteoclasts may also be mechanosensitive and that osteoclast-mediated autocrine signaling may be important in bone remodeling. In this study, we investigated the effect of FSS on nitric oxide (NO), prostaglandin E(2) (PGE(2)), and prostacyclin (PGI(2)) release by neonatal rat bone marrow-derived preosteoclast-like cells. These cells were tartrate-resistant acid phosphatase (TRAP) positive, weakly nonspecific esterase (NSE) positive, and capable of fusing into calcitonin-responsive, bone-resorbing, multinucleated cells. Bone marrow-derived preosteoclast-like cells exposed for 6 h to a well-defined FSS of 16 dynes/cm(2) produced NO at a rate of 7.5 nmol/mg protein/h, which was 10-fold that of static controls. This response was completely abolished by 100 microM N(G)-amino-L-arginine (L-NAA). Flow also stimulated PGE(2) production (3.9 microg/mg protein/h) and PGI(2) production (220 pg/mg protein/h). L-NAA attenuated flow-induced PGE(2) production by 30%, suggesting that NO may partially modulate PGE(2) production. This is the first report demonstrating that marrow derived cells are sensitive to FSS and that autocrine signaling in these cells may play an important role in load-induced remodeling and signal transduction in bone. Copyright 2000 Academic Press.

NASA Discipline Cell Biology↗

Navier-Stokes analysis of an oxidizer turbine blade with tip clearance

The Gas Generator Oxidizer Turbine (GGOT) Blade is being analyzed by various investigators under the NASA MSFC sponsored Turbine Stage Technology Team design effort. The present work concentrates on the tip clearance region flow and associated losses; however, flow details for the passage region are also obtained in the simulations. The present calculations simulate the rotor blade row in a rotating reference frame with the appropriate coriolis and centrifugal acceleration terms included in the momentum equation. The upstream computational boundary is located about one axial chord from the blade leading edge. The boundary conditions at this location were determined by using a Euler analysis without the vanes to obtain approximately the same flow profiles at the rotor as were obtained with the Euler stage analysis including the vanes. Inflow boundary layer profiles are then constructed assuming the skin friction coefficient at both the hub and the casing. The downstream computational boundary is located about one axial chord from the blade trailing edge, and the circumferentially averaged static pressure at this location was also obtained from the Euler analysis. Results were obtained for the 3-D baseline GGOT geometry at the full scale design Reynolds number. Details of the clearance region flow behavior and blade pressure distributions were computed. The spanwise variation in blade loading distributions are shown, and circumferentially averaged spanwise distributions of total pressure, total temperature, Mach number, and flow angle are shown at several axial stations. The spanwise variation of relative total pressure loss shows a region of high loss in the region near the casing. Particle traces in the near tip region show vortical behavior of the fluid which passes through the clearance region and exits at the downstream edge of the gap.

Gibeling, Howard J.↗

The organized nature of a turbulent trailing vortex

The turbulence structure of a trailing vortex produced at the juncture of a flow aligned cylinder and a pair of oppositely loaded airfoils is analyzed. The freestream turbulence intensity in this study varies from 0.32 to 1.48 percent, the vortex Reynold number varies from 15000 to 25000, and the Rossby number varies from 0.65 to 0.81. Within this parameter range, it is shown that the screens, but not the freestream turbulence level, are able to produce significant variations in the turbulence structure of the vortex, and that the turbulent structure is determined by the Rossby number and not the vortex Reynolds number. It is noted that the core is dynamic and an organized exchange of momentum takes place between the outer flow and the core region of the vortex. The vortex structure in the trailing vortex having the lowest Rossby number is considered.

Bandyopadhyay, Promode R.↗

Strongly Heated Turbulent Flow in a Channel with Pin Fins

Large-eddy simulations (LES) were performed to study the turbulent flow in a channel of height H with a staggered array of pin fins with diameter D = H/2 as a function of heating loads that are relevant to the cooling of turbine blades and vanes. The following three heating loads were investigated—wall-to-coolant temperatures of Tw/Tc = 1.01, 2.0, and 4.0 - where the Reynolds number at the channel inlet was 10,000 and the back pressure at the channel outlet was 1 bar. For the LES, two different subgrid-scale models—the dynamic kinetic energy model (DKEM) and the wall-adapting local eddy-viscosity model (WALE)—were examined and compared. This study was validated by comparing with data from direct numerical simulation and experimental measurements. The results obtained show high heating loads to create wall jets next to all heated surfaces that significantly alter the structure of the turbulent flow. Results generated on effects of heat loads on the mean and fluctuating components of velocity and temperature, turbulent kinetic energy, the anisotropy of the Reynolds stresses, and velocity-temperature correlations can be used to improve existing RANS models.

42 ENGINEERING↗

On the Feasibility of Market Manipulation and Energy Storage Arbitrage via Load-Altering Attacks

Around the globe, electric power networks are transforming into complex cyber–physical energy systems (CPES) due to the accelerating integration of both information and communication technologies (ICT) and distributed energy resources. While this integration improves power grid operations, the growing number of Internet-of-Things (IoT) controllers and high-wattage appliances being connected to the electric grid is creating new attack vectors, largely inherited from the IoT ecosystem, that could lead to disruptions and potentially energy market manipulation via coordinated load-altering attacks (LAAs). In this article, we explore the feasibility and effects of a realistic LAA targeted at IoT high-wattage loads connected at the distribution system level, designed to manipulate local energy markets and perform energy storage (ES) arbitrage. Realistic integrated transmission and distribution (T&D) systems are used to demonstrate the effects that LAAs have on locational marginal prices at the transmission level and in distribution systems adjacent to the targeted network.

25 ENERGY STORAGE↗

Application of CFD Analysis to Design Support and Problem Resolution for ASRM and RSRM

The use of Navier-Stokes CFD codes to predict the internal flow field environment in a solid rocket motor is a very important analysis element during the design phase of a motor development program. These computational flow field solutions uncover a variety of potential problems associated with motor performance as well as suggesting solutions to these problems. CFD codes have also proven to be of great benefit in explaining problems associated with operational motors such as in the case of the pressure spike problem with the STS-54B flight motor. This paper presents results from analyses involving both motor design support and problem resolution. The issues discussed include the fluid dynamic/mechanical stress coupling at field joints relative to significant propellant deformations, the prediction of axial and radial pressure gradients in the motor associated with motor performance and propellant mechanical loading, the prediction of transition of the internal flow in the motor associated with erosive burning, the accumulation of slag at the field joints and in the submerged nozzle region, impingement of flow on the nozzle nose, and pressure gradients in the nozzle region of the motor. The analyses presented in this paper have been performed using a two-dimensional axisymmetric model. Fluent/BFC, a three dimensional Navier-Stokes flow field code, has been used to make the numerical calculations. This code utilizes a staggered grid formulation along with the SIMPLER numerical pressure-velocity coupling algorithm. Wall functions are used to represent the character of the viscous sub-layer flow, and an adjusted k-epsilon turbulence model especially configured for mass injection internal flows, is used to model the growth of turbulence in the motor port. Conclusions discussed in this paper consider flow field effects on the forward, center, and aft propellant grains except for the head end star grain region of the forward propellant segment. The field joints and the submerged nozzle are discussed as well. Conclusions relative to both the design evaluation of the ASRM and the RSRM scenarios explaining the pressure spikes were based on the flow field solutions presented in this paper.

Dill, Richard A.↗

Prediction of Undsteady Flows in Turbomachinery Using the Linearized Euler Equations on Deforming Grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. The model accounts for blade loading, blade geometry, shock motion, and wake motion. Assuming that the unsteadiness in the flow is small relative to the nonlinear mean solution, the unsteady Euler equations can be linearized about the mean flow. This yields a set of linear variable coefficient equations that describe the small amplitude harmonic motion of the fluid. These linear equations are then discretized on a computational grid and solved using standard numerical techniques. For transonic flows, however, one must use a linear discretization which is a conservative linearization of the non-linear discretized Euler equations to ensure that shock impulse loads are accurately captured. Other important features of this analysis include a continuously deforming grid which eliminates extrapolation errors and hence, increases accuracy, and a new numerically exact, nonreflecting far-field boundary condition treatment based on an eigenanalysis of the discretized equations. Computational results are presented which demonstrate the computational accuracy and efficiency of the method and demonstrate the effectiveness of the deforming grid, far-field nonreflecting boundary conditions, and shock capturing techniques. A comparison of the present unsteady flow predictions to other numerical, semi-analytical, and experimental methods shows excellent agreement. In addition, the linearized Euler method presented requires one or two orders-of-magnitude less computational time than traditional time marching techniques making the present method a viable design tool for aeroelastic analyses.

Clark, William S.↗

Sensitivity Analysis of Tunable Equation of State Material Model In Pulsed Mercury Target Simulation

A pulsed neutron spallation target is subjected to very short but intense loads from repeated proton pulses. Approximately 60% of the energy from each proton pulse is deposited into the mercury target material and the stainless-steel target structure, leading to a high-pressure region in both the stationary target structure and the flowing mercury. The high-pressure region propagates and leads to fluid-structure interaction. The resultant loading on the target structure containing liquid mercury is difficult to predict, although various simulation approaches and material models for the mercury have been tried. To date, the best match of simulation to experimental data is obtained by using an equation of state (EOS) material model with a specified tensile cutoff pressure, which simulates the cavitation threshold. The inclusion of a threshold to represent cavitation is key to the successful predictions of stress waves triggered by the high-energy pulse striking the mercury and vessel. However, recent measurements of target structure strain show that significant discrepancies remain between the measured and simulated strain values in the EOS mercury model. These differences grow when noncondensable helium gas is intentionally injected into the flowing mercury to reduce the loading on the structure. An EOS-based proportional–integral–derivative (PID) mercury model has been proposed to reduce the gap between the measured and simulated vessel strain responses for targets with gas injection. The conceptual and numerical description and initial investigation of the PID model are presented in previous work. Further studies of this PID model — including the sensitivity of the structure’s strain response to model parameters (the tensile cutoff, PID parameters Kp, Ki, and Kd) — are reported in this article. Results show the strain response is more sensitive to changes in the tensile cutoff value than to changes in the model parameters Kp, Ki, and Kd. These results will aid in future work where the model parameters will be optimized to match simulation data to strain measurements.

Lin, Lianshan↗

Comprehensive Forced Response Analysis of J2X Turbine Bladed-Discs with 36- Degree Variation in CFD Loading

At NASA/MSFC, Structural Dynamics personnel continue to perform advanced analysis for the turbomachinery in the J2X Rocket Engine, which is under consideration for the new Space Launch System. One of the most challenging analyses in the program is predicting turbine blade structural capability. Resonance was predicted by modal analysis, so comprehensive forced response analyses using high fidelity cyclic symmetric finite element models were initiated as required. Analysis methodologies up to this point have assumed the flow field could be fully described by a sector, so the loading on every blade would be identical as it travelled through it. However, in the J2X the CFD flow field varied over the 360 deg of a revolution because of the flow speeds and tortuous axial path. MSFC therefore developed a complex procedure using Nastran Dmap's and Matlab scripts to apply this circumferentially varying loading onto the cyclically symmetric structural models to produce accurate dynamic stresses for every blade on the disk. This procedure is coupled with static, spin, and thermal loading to produce high cycle fatigue safety factors resulting in much more accurate analytical assessments of the blades.

Elrod, David↗

Performance of a Compression-ignition Engine with a Precombustion Chamber Having High-Velocity Air Flow

Presented here are the results of performance tests made with a single-cylinder, four stroke cycle, compression-ignition engine. These tests were made on a precombustion chamber type of cylinder head designed to have air velocity and tangential air flow in both the chamber and cylinder. The performance was investigated for variable load and engine speed, type of fuel spray, valve opening pressure, injection period and, for the spherical chamber, position of the injection spray relative to the air flow. The pressure variations between the pear-shaped precombustion chamber and the cylinder for motoring and full load conditions were determined with a Farnboro electric indicator. The combustion chamber designs tested gave good mixing of a single compact fuel spray with the air, but did not control the ensuing combustion sufficiently. Relative to each other, the velocity of air flow was too high, the spray dispersion by injection too great, and the metering effect of the cylinder head passage insufficient. The correct relation of these factors is of the utmost importance for engine performance.

Spanogle, J A↗

Aerothermal Uncertainty Quantification of Deployable Entry Technologies Using Multi-Fidelity Modeling

The objective of this work was to investigate the use of a co-Kriging based multi-fidelity modeling approach with a Monte Carlo uncertainty quantification analyses of surface heating on hypersonic inflatable aerodynamic decelerator vehicles and adaptable, deployable entry placement technology vehicles in Mars entry. A previously developed co-Kriging based multi-fidelity modeling approach was used to model the laminar and turbulent convective and radiative heat fluxes along the vehicle. Monte Carlo analyses of surface heat load for nine different vehicle nose radii was performed for both vehicles, assuming the same thermal protection system is used for both vehicles. The maximum heat loads and heat load uncertainties were found to occur at either the stagnation point or the vehicle shoulder. The maximum 95% confidence intervals for the total heat load for the hypersonic inflatable aerodynamic decelerator vehicle were found to be [-9.27, 7.81] % of the nominal value and [-6.09, 13.1] % of the nominal value for laminar and turbulent flow, respectively. The maximum 95% confidence intervals for the total heat load for the adaptable, deployable entry placement technology vehicle were found to be [-8.85, 7.22] % of the nominal value and [-9.33, 9.35] % of the nominal value for laminar and turbulent flow, respectively. The difference in heating uncertainties based on vehicle geometry between the two vehicle technologies was found to be minimal.

Mario Santos↗

Unsteady airloads in supercritical transonic flows

Results obtained from a simplified theory of unsteady perturbations of supercritical two-dimensional transonic flows, introduced in an earlier paper, are presented. Unsteady loads generated by an oscillating flap and by airfoils oscillating as a whole are given with comparisons to experimental results and finite difference solutions. The theory, which was originally formulated for symmetric airfoils at zero mean angle of attack, is extended to treat asymmetric mean flows.

Williams, M. H.↗

Calculation of unsteady linearized Euler flows in cascades using harmonically deforming grids

A method for calculating unsteady, inviscid, compressible flows in cascades is presented. Using the linearized Euler technique, the flow is decomposed into a steady or mean flow plus a harmonically varying small disturbance flow. The equations that describe the small disturbance flow are linear variable coefficient equations, and are solved using a pseudo-time time marching Lax-Wendroff technique. Unlike previous linearized methods, however, the solution is computed on a harmonically deforming computational grid that conforms to the motion of the vibrating airfoils. The mean flow and perturbation flow solutions are defined in the deforming coordinate system rather than in a coordinate system fixed in space. Hence, no extrapolation terms are required to implement the upwash boundary conditions at the airfoil surfaces, significantly improving the accuracy of the method. For transonic flow calculations, unsteady shock motions are modelled using shock capturing. The unsteady loads due to the shock motion are then seen as pressure impulses. Representative computational results are presented for transonic channel flows and subsonic and transonic cascade flows.

Hall, Kenneth C.↗

On the Importance of Spatial Resolution for Flap Side Edge Noise Prediction

A spatial resolution study of flap tip flow and the effects on the farfield noise signature for an 18%-scale, semispan Gulfstream aircraft model are presented. The NASA FUN3D unstructured, compressible Navier-Stokes solver was used to perform the highly resolved, time-dependent, detached eddy simulations of the flow field associated with the flap for this high-fidelity aircraft model. Following our previous work on the same model, the latest computations were undertaken to determine the causes of deficiencies observed in our earlier predictions of the steady and unsteady surface pressures and off-surface flow field at the flap tip regions, in particular the outboard tip area, where the presence of a cavity at the side-edge produces very complex flow features and interactions. The present results show gradual improvement in steady loading at the outboard flap edge region with increasing spatial resolution, yielding more accurate fluctuating surface pressures, off-surface flow field, and farfield noise with improved high-frequency content when compared with wind tunnel measurements. The spatial resolution trends observed in the present study demonstrate that the deficiencies reported in our previous computations are mostly caused by inadequate spatial resolution and are not related to the turbulence model.

Mineck, Raymond E.↗

Aeroelastic Modeling and Full-Scale Loads Measurements for Investigation of Single-Axis PV Tracker Wind-Driven Dynamic Instabilities

While wind tunnel testing and proprietary industry modeling tools have been used for years to design and develop PV tracker systems, wind induced dynamic failures are becoming more prevalent and high visibility. Designers and manufactures have reacted by developing new systems or add-on products to improve system dynamics adding to overall system costs. To better understand the physics and aeroelastic behavior that leads to dynamic instabilities NREL researches have taken up a first of a kind study to both develop open source aeroelastic modeling tools and measure tracker loads on a single axis full-scale tracker system at a high wind site. These modeling tools can simulate the fluid-structure interaction driving torsional instabilities under a wide range of turbulent inflow conditions and stow angles. The results of these simulations reveal the flow features associated with panel rotation and the induced loads in the structure. The experimental load measurements can be used to validate models and provide a quantifiable understanding of tracker dynamics, critical loads paths, identify instability markers, inform resilient design, and suggest the most favorable stow approach.

14 SOLAR ENERGY↗

Flexure and the role of inplane force around coronae on Venus

Large coronae on Venus, such as Artemis and Latona, are rimmed by conspicuous trenches and associated outer rises. Sandwell and Schubert have observed that these systems resemble terrestrial subduction zones in planform and have succeeded in fitting an elastic plate bending equation to the inferred flexural topography. However, the first zero crossing bending moments required are -2.5 x 10(exp 17) N for Artemis and -5.0 x 10(exp 16) N for Latona. Since these moments are similar in magnitude to those of subducting slabs on Earth, a rollback subduction mechanism was proposed to explain the flexure around the largest coronae, although a differential thermal subsidence model is sufficient to account for the topography around some coronae. The purpose is to investigate the effect of inplane force as a possible alternative to large applied moments in producing flexure at Artemis and Latona. The close correlation of gravity to topography on Venus implies the absence of a low viscosity zone and the strong coupling of the lithosphere to mantle convection. If coronae are the surface manifestations of mantle plumes, they may be the sites of active convective stress coupling. As the upwelling reaches the lithosphere, it spreads radially outward, inducing shear tractions on the base of the plate. In addition, the hot, expanding corona may load the surrounding plate horizonally. Both the basal shear stresses and radial loading can be treated as an equivalent compressive inplane force in the mechanical lithosphere, which contributes to the bending of the outlying plate. Using a model that relates inplane force to the measured gravity anomalies, a rough value of the inplane force at Artemis was calculated. Recent Pioneer Venus spherical harmonic gravity models indicate a geoid anomaly of about 75 m over Artemis, which corresponds to an estimated inplane force on the order of -1x10(exp 13) N/m. The gravity model is unable to resolve Latona, but an inplane force of similar dimensions is assumed. The maximum possible inplane force based on the expected rheology can be constrained by using the approximate 5 K/km thermal gradient inferred from the best fit 30 km elastic plate at Artemis and Latona. For a dry olivine flow law in the upper mantle, the compressional load limit of the 60 km thick mechanical lithosphere is -4 x 10(exp 13) N/m. This value is equivalent to a load of -8 x 10(exp 13) N/m on a 30 km thick elastic plate.

Brown, C. David↗