Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity model evaluation”

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 163 records · Page 9

Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska

The Kodiak Islands are located approx.120 to 250 km from the Alaska-Aleutian Trench and are within the southern extent of the 1964 Prince William Sound (M(sub W) = 9.2) earthquake rupture and aftershock zone. Here we report new campaign GPS results (1993-2001) from northeastern Kodiak and reprocessed GPS results (1993-1997) from southwestern Kodiak. The rate and orientation of the horizontal velocities, relative to a fixed North America, range from 29.7 +/- 1.7 mm/yr at N30.3degW +/- 3.3deg, located approx.120 km from the deepest point of the trench, to 8.0 +/- 1.3 mm/yr at N62.4degW +/- 9.3deg, located approx.230 km from the trench. We evaluated alternate models of coseismic and interseismic slip to test the importance of the mechanisms that account for surface deformation rates. Near the Gulf of Alaska coastal region of Kodiak the horizontal velocity can be accounted for primarily by the viscoelastic response to plate motion and a locked main thrust zone (MTZ), down-dip creep, and to a lesser extent, slip in the 1964 earthquake. Further inland the dominant mechanisms that account for post-1964 uplift rates are time-dependent, down-dip creep and a locked MTZ; for the horizontal velocity component southwest translation of western Kodiak may be important as well. Based on the pre-1964 and post-1964 earthquake pattern of interseismic earthquakes, we suggest that between the occurrences of great earthquakes like the 1964 event, more moderate to large earthquakes occur in the southwestern Kodiak region than near northeastern Kodiak .

Sauber, Jeanne↗

A method for predicting the noise levels of coannular jets with inverted velocity profiles

A coannular jet was equated with a single stream equivalent jet with the same mass flow, energy, and thrust. The acoustic characteristics of the coannular jet were then related to the acoustic characteristics of the single jet. Forward flight effects were included by incorporating a forward exponent, a Doppler amplification factor, and a Strouhal frequency shift. Model test data, including 48 static cases and 22 wind tunnel cases, were used to evaluate the prediction method. For the static cases and the low forward velocity wind tunnel cases, the spectral mean square pressure correlation coefficients were generally greater than 90 percent, and the spectral sound pressure level standard deviation were generally less than 3 decibels. The correlation coefficient and the standard deviation were not affected by changes in equivalent jet velocity. Limitations of the prediction method are also presented.

Russell, J. W.↗

A study of Reynolds-stress closure model

A hybrid model of the Reynolds stress closure was developed. This model was tested for various sizes of step flow, and the computed Reynolds stress behavior was compared with experimental data. The third order closure model was reviewed. Transport equations for the triple velocity correlation were developed and implemented in a numerical code to evaluate the behavior of the triple velocity products in various regions of the flow field including recirculating, reattaching, and redeveloping flow domains.

Amano, R. S.↗

Wave propagation in composite media and material characterization

Characteristics of wave propagation in an undamaged composite medium are influenced by many factors, the most important of which are: microstructure, constituent properties, interfaces, residual stress fields, and ply lay-ups. Measurements of wave velocities, attenuation, and dispersion provide a powerful tool for nondestructive evaluation of these properties. Recent developments are reviewed for modeling ultrasonic wave propagation in fiber and particle-reinforced composite media. Additionally, some modeling studies are reviewed for the effects of interfaces and layering on attenuation and dispersion. These studies indicate possible ways of characterizing material properties by ultrasonic means.

Datta, Subhendu K.↗

Development and evaluation of a musculoskeletal model of the elbow joint complex

This paper describes the development and evaluation of a musculoskeletal model that represents human elbow flexion-extension and forearm pronation-supination. The length, velocity, and moment arm for each of the eight musculotendon actuators were based on skeletal anatomy and position. Musculotendon parameters were determined for each actuator and verified by comparing analytical torque-angle curves with experimental joint torque data. The parameters and skeletal geometry were also utilized in the musculoskeletal model for the analysis of ballistic elbow joint complex movements. The key objective was to develop a computational model, guided by parameterized optimal control, to investigate the relationship among patterns of muscle excitation, individual muscle forces, and movement kinematics. The model was verified using experimental kinematic, torque, and electromyographic data from volunteer subjects performing ballistic elbow joint complex movements.

Gonzalez, Roger V.↗

Crashworthy Evaluation of a 1/5-Scale Model Composite Fuselage Concept

A 1/5-scale model composite fuselage concept for light aircraft and rotorcraft has been developed to satisfy structural and flight loads requirements and to satisfy design goals for improved crashworthiness. The 1/5-scale model fuselage consists of a relatively rigid upper section which forms the passenger cabin, a stiff structural floor, and an energy absorbing subfloor which is designed to limit impact forces during a crash event. The focus of the present paper is to describe the crashworthy evaluation of the fuselage concept through impact testing and finite element simulation using the nonlinear, explicit transient dynamic code, MSC/DYTRAN. The energy absorption behavior of two different subfloor configurations was determined through quasi-static crushing tests. For the dynamic evaluation, each subfloor configuration was incorporated into a 1/5-scale model fuselage section, which was impacted at 31 ft/s vertical velocity onto a rigid surface. The experimental data demonstrate that the fuselage section with a foam-filled subfloor configuration satisfied the impact design requirement. In addition, the fuselage section maintained excellent energy absorption behavior for a 31 ft/s vertical drop test with a 15 deg-roll impact attitude. Good correlation was obtained between the experimental data and analytical results for both impact conditions.

Jackson, Karen E.↗

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗

Numerical Simulations of STOVL Hot Gas Ingestion in Ground Proximity Using a Multigrid Solution Procedure

A multi grid solution procedure for the numerical simulation of turbulent flows in complex geometries has been developed. A Full Multigrid-Full Approximation Scheme (FMG-FAS) is incorporated into the continuity and momentum equations, while the scalars are decoupled from the multi grid V-cycle. A standard kappa-Epsilon turbulence model with wall functions has been used to close the governing equations. The numerical solution is accomplished by solving for the Cartesian velocity components either with a traditional grid staggering arrangement or with a multiple velocity grid staggering arrangement. The two solution methodologies are evaluated for relative computational efficiency. The solution procedure with traditional staggering arrangement is subsequently applied to calculate the flow and temperature fields around a model Short Take-off and Vertical Landing (STOVL) aircraft hovering in ground proximity.

Wang, Gang↗

GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design

GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Velocity selection at large undercooling in a two-dimensional nonlocal model of solidification

The formation of needle-crystal dendrites from an undercooled melt is investigated analytically, applying the method of Caroli et al. (1986) to Langer's (1980) symmetric two-dimensional nonlocal model of solidification with finite anisotropy in the limit of large undercooling. A solution based on the WKB approximation is obtained, and a saddle-point evaluation is performed. It is shown that needle-crystal solutions exist only if the capillary anisotropy is nonzero, in which case a particular value of the growth velocity can be selected. This finding and the expression for the dependence of the selected velocity on the singular perturbation parameter and the strength of the anisotropy are found to be in complete agreement with the results of a boundary-layer model (Langer and Hong, 1986).

Barbieri, Angelo↗

Convective Updraft and Downdraft Characteristics of Continental Mesoscale Convective Systems in the Model Gray Zone

Here, the "gray zone" of convective modeling is defined as the range of horizontal grid spacings (Δx) at which turbulent transport processes are only partially resolved by the dynamics of the numerical model. This zone typically covers Δx from a few kilometers to several hundred meters, wherein the realistic representation of convective cloud processes can be challenging. This study characterizes the convective draft behaviors at multiple Δx across the gray zone and determines the appropriate Δx that can reliably capture these salient convective properties. We perform an ensemble of idealized simulations of mesoscale convective systems (MCS) using the Weather Research and Forecasting model at various Δx from 4 km to 250 m over the central U.S. An evaluation of key MCS kinematic properties is constrained using unique, long-term vertical velocity estimates obtained by radar wind profilers deployed by the Department of Energy Atmospheric Radiation Measurement user facility. MCS simulations for all Δx tested overestimate (underestimate) the probabilities of convective updrafts (downdrafts) compared to the observations. In terms of the convective draft intensity, finer-Δx models overestimate the updraft intensity, while the opposite is found for downdrafts. Moving from Δx = 4 km to 250 m, downdrafts become stronger and more frequent especially at middle and upper levels, attributed to additional drag from an increasing graupel frequency and compensation of enhanced updraft velocity. Simulated draft characteristics, including core size, intensity, and probability of occurrence, exhibit pronounced changes at Δx greater than 500 m, but suggest less sensitivity when Δx is reduced below 500 m.

54 ENVIRONMENTAL SCIENCES↗

Toward integrating high-fidelity CFD approaches in the thermal-hydraulic analysis of turbulent dry cask systems

Nuclear power plants have been supplying resilient and reliable electricity for decades, contributing to energy independence of the U.S.. However, nuclear waste management remains one of the most significant challenges in the industry. The safety of dry cask storage systems relies heavily on their thermal-hydraulic performance. Computational Fluid Dynamics (CFD) simulations are often used to demonstrate this performance and ensure that the system design meets safety standards. This study presents reduced numerical models for various types of dry cask systems. These numerical models can produce efficient and fast results based on the employed modeling strategies. Additionally, the study uses a novel approach to high-fidelity simulations to evaluate modeling assumptions in dry cask modeling. Large Eddy Simulations (LES) are used for this purpose, particularly in regions where fluid velocity is relatively high and the turbulence characteristics become important. Furthermore, the results of these high-fidelity simulations will enhance the interpretation of outcomes produced from a lower-fidelity CFD model.

CFD↗

Orographic Cirrus and Its Radiative Forcing in NCAR CAM6

Cirrus clouds play an important role in the Earth's radiative energy budget, thereby affecting the climate state and climate change. Orographic gravity wave (OGW)-induced sub-grid scale vertical velocity (i.e., cooling rate) is not resolved by large-scale models and its impact on ice formation in cirrus clouds is not well quantified. In this study, one sub-grid scale OGW scheme (e.g., McFarlane) is used in the Community Atmosphere Model version 6 (CAM6) to generate vertical velocity variance (σw) for cirrus formation. Results from the default model and simulations with the OGW-induced σw are evaluated against the DOE ARM Small Particles in Cirrus (SPARTICUS) campaign observations. The OGW based on the McFarlane scheme increases the sub-grid scale σw over mountains compared to the default model and improves the model agreement with the SPARTICUS observations. Larger σ w due to OGWs can trigger more frequent homogeneous nucleation in orographic cirrus and generates a higher number concentration of ice crystals observed during the SPARTICUS campaign. Moreover, our evaluation of the model simulations against satellite observations indicates that the McFarlane scheme generates high in-cloud ice number concentrations (>200 L -1 ) in the upper troposphere over mountains and high plateaus at mid- and high-latitudes of the winter hemisphere as shown in the observations. More ice crystals with smaller sizes absorb more infrared radiation (+0.523 ± 0.125 W m -2 ). The net radiative cloud forcing change at the top of the atmosphere is +0.330 W m -2 due to the orographic cirrus.

54 ENVIRONMENTAL SCIENCES↗

Prediction of overall and blade-element performance for axial-flow pump configurations

A method and a digital computer program for prediction of the distributions of fluid velocity and properties in axial flow pump configurations are described and evaluated. The method uses the blade-element flow model and an iterative numerical solution of the radial equilbrium and continuity conditions. Correlated experimental results are used to generate alternative methods for estimating blade-element turning and loss characteristics. Detailed descriptions of the computer program are included, with example input and typical computed results.

Serovy, G. K.↗

Radiation transfer through a model sunspot

A cylindrically symmetric model for a sunspot atmosphere using the similarity principle of Schlueter and Temesvary for the magnetic field configuration is presented. The equations of magnetostatic equilibrium are used, augmented by a radial Evershed flow. The LTE radiative transfer equations for the Stokes vector were solved under a variety of conditions for a ray emerging from a typical penumbral point. The contribution from isolated lines to the broadband circular polarization in sunspot penumbrae is evaluated using a more realistic model sunspot atmosphere than has hitherto been considered. Results indicate that the inclusion of a velocity field along the magnetic field vector is unable to give a net circular polarization of sufficient magnitude, although the variation with the angle between the line-of-sight and the magnetic field vector is in qualitative agreement with observations. The corresponding results for the net linear polarization are satisfactory.

Landman, D. A.↗

A data base for three-dimensional all-interference code evaluation

A validation of a measured boundary condition technique was carried out to demonstrate the feasibility of a wall interference assessment/correction (WIAC) system. An experimental evaluation was also carried out to compare performances of various techniques, to define the number of necessary boundary measurements for accurate assessment/corrections and to define the envelope of test conditions for which accurate assessment/corrections are achieved. The relative merits of a WIAC system and an adaptive wall tunnel are compared. The measurement surface boundary data is performed with a system of two rotating pipes. These pipes sweep out a cylindrical measurement surface near the tunnel walls, approximately one inch from the wall at the closest point. The experimental model was specially designed and fabricated for the adaptive wall experiments. The model is a wing/tail/body configuration with swept lifting surface. The boundary data taken in Tunnel 1T with the rotating pipe system has been shown to offer several attractive features for WIAC code evaluation. Good spatial resolution of measurements is achieved and measurements are made upstream and downstream of the model. Also, two velocity components are determined.

Sickles, W. L.↗

Large-eddy simulation of a backward facing step flow using a least-squares spectral element method

We report preliminary results obtained from the large eddy simulation of a backward facing step at a Reynolds number of 5100. The numerical platform is based on a high order Legendre spectral element spatial discretization and a least squares time integration scheme. A non-reflective outflow boundary condition is in place to minimize the effect of downstream influence. Smagorinsky model with Van Driest near wall damping is used for sub-grid scale modeling. Comparisons of mean velocity profiles and wall pressure show good agreement with benchmark data. More studies are needed to evaluate the sensitivity of this method on numerical parameters before it is applied to complex engineering problems.

Chan, Daniel C.↗

DBD Plasma Actuators for Flow Control in Air Vehicles and Jet Engines - Simulation of Flight Conditions in Test Chambers by Density Matching

Dielectric Barrier Discharge (DBD) Plasma actuators for active flow control in aircraft and jet engines need to be tested in the laboratory to characterize their performance at flight operating conditions. DBD plasma actuators generate a wall-jet electronically by creating weakly ionized plasma, therefore their performance is affected by gas discharge properties, which, in turn, depend on the pressure and temperature at the actuator placement location. Characterization of actuators is initially performed in a laboratory chamber without external flow. The pressure and temperature at the actuator flight operation conditions need to be simultaneously set in the chamber. A simplified approach is desired. It is assumed that the plasma discharge depends only on the gas density, while other temperature effects are assumed to be negligible. Therefore, tests can be performed at room temperature with chamber pressure set to yield the same density as in operating flight conditions. The needed chamber pressures are shown for altitude flight of an air vehicle and for jet engines at sea-level takeoff and altitude cruise conditions. Atmospheric flight conditions are calculated from standard atmosphere with and without shock waves. The engine data was obtained from four generic engine models; 300-, 150-, and 50-passenger (PAX) aircraft engines, and a military jet-fighter engine. The static and total pressure, temperature, and density distributions along the engine were calculated for sea-level takeoff and for altitude cruise conditions. The corresponding chamber pressures needed to test the actuators were calculated. The results show that, to simulate engine component flows at in-flight conditions, plasma actuator should be tested over a wide range of pressures. For the four model engines the range is from 12.4 to 0.03 atm, depending on the placement of the actuator in the engine. For example, if a DBD plasma actuator is to be placed at the compressor exit of a 300 PAX engine, it has to be tested at 12.4 atm for takeoff, and 6 atm for cruise conditions. If it is to be placed at the low-pressure turbine, it has to be tested at 0.5 and 0.2 atm, respectively. These results have implications for the feasibility and design of DBD plasma actuators for jet engine flow control applications. In addition, the distributions of unit Reynolds number, Mach number, and velocity along the engine are provided. The engine models are non-proprietary and this information can be used for evaluation of other types of actuators and for other purposes.

Ashpis, David E.↗