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 937 records · Page 52

Validation of Computational Fluid Dynamics Simulations for Biological Performance Assessment in Hydropower units (Final Report)

The biological performance assessment (BioPA) toolset developed by Pacific Northwest National Laboratory (PNNL) estimates the relative biological performance of fish passage at a hydroelectric power turbine unit. The tool is based on the use of computational fluid dynamics (CFD) and fish biological response relationships. The recent release, BioPA-v3, is based on directly computed trajectory and collision of material Lagrangian particles using CFD simulation codes rather than the prior version that relies on Tecplot to compute streamtrace trajectories. Before modifying the toolset, a series of validation tests were performed at the various steps of modification in the toolset. Validation is a critical step of any numerical investigation that reflects the accuracy and reliability of the predicted results. It raises the confidence level of the user to use the modified version of the BioPA toolset. Several test cases were simulated and compared, where available, to observed data. The trajectory and collision of the small spherical and cylindrical particles in a water flume were compared to in-house experiments. The CFD predicted collision rate and flow field compared well with experimental observation for vane array and large cylinder as target bodies. Next, the CFD-predicted flow field and hydraulic performance of a laboratory-scale model of a Francis turbine was also successfully validated. Note that the trajectory of the particles is significantly affected by the flow field in such extreme conditions. In addition to the particle trajectories and flow field, the collision detection method employed in the CFD simulations was also successfully validated. The CFD predicted impact velocity, collision time, velocity, and trajectory of a sphere excellently matched with analytical value for a bouncing ball in the elastic collision. A similar approach was also tested and successfully validated for a collision of sphere with a 45° inclined plane. After successfully validating different cases, the BioPA toolset was modified to use direct output of the CFD prediction and the new version can be used in evaluating biological performance at hydroelectric turbines.

13 HYDRO ENERGY↗

Establishing a relationship between maximum torque production of isolated joints to simulate EVA ratchet push-pull maneuver: A case study

As manned exploration of space continues, analytical evaluation of human strength characteristics is critical. These extraterrestrial environments will spawn issues of human performance which will impact the designs of tools, work spaces, and space vehicles. Computer modeling is an effective method of correlating human biomechanical and anthropometric data with models of space structures and human work spaces. The aim of this study is to provide biomechanical data from isolated joints to be utilized in a computer modeling system for calculating torque resulting from any upper extremity motions: in this study, the ratchet wrench push-pull operation (a typical extravehicular activity task). Established here are mathematical relationships used to calculate maximum torque production of isolated upper extremity joints. These relationships are a function of joint angle and joint velocity.

Pandya, Abhilash↗

Evaluation of Mixed-Phase Microphysics Within Winter Storms Using Field Data and In Situ Observations

Snow prediction within models is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33-km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Cold-season Precipitation Experiment (GCPEx) experiment over southern Ontario, as well as a few years (12 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. The GCPex presentation will focus on verification using aircraft spirals through warm frontal snow band event on 18 February 2012. All the BMPs realistically simulated the structure of the band and the vertical distribution of snow/ice aloft, except the SBU-YLIN overpredicted slightly and Thompson (THOM) underpredicted somewhat. The Morrison (MORR) scheme produced the best slope size distribution for snow, while the Stony Brook (SBU) underpredicted and the THOM slightly overpredicted. Those schemes that have the slope intercept a function of temperature (SBU and WSM6) tended to perform better for that parameter than others, especially the fixed intercept in Goddard. Overall, the spread among BMPs was smaller than in other studies, likely because there was limited riming with the band. For the 15 cases at SBNY, which include moderate and heavy riming events, the non-spherical snow assumption (THOM and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM is approx. 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM, and MORR schemes were approx. 0.25 m/s too slow, while the SBU was 0.25 to 0.5 m/s too fast because of some excessive cloud water issues.

Colle, Brian A.↗

Observations of wind farm wake recovery at an operating wind farm

Abstract. The interplay of momentum surrounding wind farms significantly influences wake recovery, affecting the speed at which wakes return to their freestream velocities. Under stable atmospheric conditions, wind farm wakes can extend over considerable distances, leading to sustained vertical momentum flux downstream, with variations observed throughout the diurnal cycle. Particularly in regions such as the US Great Plains, stable conditions can induce low-level jets (LLJs), impacting wind farm performance and power output. This study examines the implications of wake recovery using long-term observations of vertical momentum flux profiles across diverse atmospheric conditions. In these observations, several key findings were observed, such as (a) LLJ heights being altered downstream of a wind farm, especially when the LLJs are below 250 m above ground level; (b) a notable impact of LLJ height on wake recovery being observed using momentum flux profiles at upwind and downwind locations, wherein LLJs between 250 and 500 m above ground level resulted in larger momentum transfer within the wake (i.e., smaller velocity deficit) compared to LLJs below 250 m above ground level; (c) the largest momentum flux variability being observed during stable atmospheric conditions, with non-negligible variability observed during neutral and unstable atmospheric conditions; (d) detection of wake effects almost always being observed throughout the atmospheric boundary layer height; and finally (e) enhancement of wake recovery being observed in the presence of propagating gravity waves. These insights deepen our understanding of the intricate dynamics governing wake recovery in wind farms, advancing efforts to model and predict their behavior across varying atmospheric contexts. In addition, the performance of large-eddy-simulation-based semi-empirical internal boundary layer height model estimates incorporating real-world atmospheric and turbine inputs was evaluated using observations during LLJ conditions.

17 WIND ENERGY↗

Wave propagation in anisotropic medium due to an oscillatory point source with application to unidirectional composites

The far-field displacements in an infinite transversely isotropic elastic medium subjected to an oscillatory concentrated force are derived. The concepts of velocity surface, slowness surface and wave surface are used to describe the geometry of the wave propagation process. It is shown that the decay of the wave amplitudes depends not only on the distance from the source (as in isotropic media) but also depends on the direction of the point of interest from the source. As an example, the displacement field is computed for a laboratory fabricated unidirectional fiberglass epoxy composite. The solution for the displacements is expressed as an amplitude distribution and is presented in polar diagrams. This analysis has potential usefulness in the acoustic emission (AE) and ultrasonic nondestructive evaluation of composite materials. For example, the transient localized disturbances which are generally associated with AE sources can be modeled via this analysis. In which case, knowledge of the displacement field which arrives at a receiving transducer allows inferences regarding the strength and orientation of the source, and consequently perhaps the degree of damage within the composite.

Williams, J. H., Jr.↗

New echocardiographic windows for quantitative determination of aortic regurgitation volume using color Doppler flow convergence and vena contracta

Color Doppler images of aortic regurgitation (AR) flow acceleration, flow convergence (FC), and the vena contracta (VC) have been reported to be useful for evaluating severity of AR. However, clinical application of these methods has been limited because of the difficulty in clearly imaging the FC and VC. This study aimed to explore new windows for imaging the FC and VC to evaluate AR volumes in patients and to validate this in animals with chronic AR. Forty patients with AR and 17 hemodynamic states in 4 sheep with strictly quantified AR volumes were evaluated. A Toshiba SSH 380A with a 3.75-MHz transducer was used to image the FC and VC. After routine echo Doppler imaging, patients were repositioned in the right lateral decubitus position, and the FC and VC were imaged from high right parasternal windows. In only 15 of the 40 patients was it possible to image clearly and measure accurately the FC and VC from conventional (left decubitus) apical or parasternal views. In contrast, 31 of 40 patients had clearly imaged FC regions and VCs using the new windows. In patients, AR volumes derived from the FC and VC methods combined with continuous velocity agreed well with each other (r = 0.97, mean difference = -7.9 ml +/- 9.9 ml/beat). In chronic animal model studies, AR volumes derived from both the VC and the FC agreed well with the electromagnetically derived AR volumes (r = 0.92, mean difference = -1.3 +/- 4.0 ml/beat). By imaging from high right parasternal windows in the right decubitus position, complementary use of the FC and VC methods can provide clinically valuable information about AR volumes.

NASA Discipline Cardiopulmonary↗

Airfoil self noise - Effect of scale

Key data from a comprehensive airfoil broadband self-noise study are reported. Attention here is restricted to two-dimensional sharp trailing-edge models. The models include seven NACA 0012 airfoil sections and five 'flat plate' sections with chordlengths ranging from 2.54 to 60.96 cm. Testing parameters include flow velocity, angle of attack to the flow, and boundary layer turbulence through natural transition and through tripping. Detailed aerodynamic measurements, of pertinence to the scaling problem of airfoil self-noise, were conducted in the near-wake of the trailing edges. Presented are mean and rms turbulent velocity profiles as well as boundary layer thicknesses and integral thickness parameters for a large range of conditions. The noise spectra of the self-noise sources were determined by the use of a cross-spectral technique. The spectra were normalized using the measured aerodynamic parameters in order to evaluate the most common scaling law now in use. An examination of the Reynolds number dependence of the overall self-noise levels has revealed a new and useful scaling result. This result appears to quantify the transition between turbulent-boundary-layer trailing edge noise and laminar-boundary-layer vortex shedding noise.

Brooks, T. F.↗

Prediction of Unsteady Transitional Layers in Turbomachinery Using Navier Stokes Equations

The objective of the research reported in this presentation is to develop computational techniques for the prediction of unsteady transitional flows associated with the rotor stator interaction in turbomachinery. Three low-Reynolds number turbulence models are incorporated in two unsteady Navier-Stokes codes (one is pressure based and the other is time marching with Runge-Kutta time stepping) and evaluated for accuracy in predicting the onset and the end of unsteady transitional patches due to wake passing. The best model is then used for modification and improvement for the leading edge effect. An existing steady Navier-Stokes code was modified to include pseudo-time stepping, which provided acceleration from 5 to 25 times that of the original code. A systematic validation procedure was implemented to assess the effects of the grid, artificial dissipation, physical, and the pseudo-time step for an accurate prediction of transitional flows resulting from the rotor-stator interaction. The ability of the Navier-Stokes code to predict the unsteady transitional flow on a turbomachinery blade is demonstrated. The unsteady pressure and velocity fields are in good agreement with the experimental data and the prediction from the Euler/boundary layer approach. The numerical solver was able to capture all zones (wake induced transitional strip, wake induced turbulent strip, calmed region, etc.) associated with wake induced transition in a compressor cascade. Another significant step is the assessment of k-epsilon turbulence models, including the leading edge modifications. Best results were obtained from the FLB model. The LB model predicted earlier inception of the transition and shorter transition length. Modification of the k-epsilon model was found to be essential for an accurate prediction of the unsteady transitional flow in a compressor cascade. The CH model failed to predict the unsteady transitional flow. Predicted boundary layer was turbulent from the leading edge, even with the modification of the k-epsilon model near the stagnation point. A comparison of the instantaneous shape factor, the skin friction coefficient, and the momentum thickness indicates that the Navier-Stokes predictions are reasonably good. Interaction between the upstream wake and stator wake results in shedding of unsteady vortices from the trailing edge and increased dissipation in the stator wake and, as a consequence, increased rate of decay of the stator wake. The procedure developed at Penn State should aid the designers of turbomachinery in allowing for the unsteady transitional flows due to rotor-stator interaction.

Lakshminarayana, B.↗

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

VanOverbeke, Thomas J.↗

Expansion-driven Weibel instability in magnetohydrodynamics: Linear theory for static systems

The Weibel dispersion relation is obtained in the magnetohydrodynamic (MHD) approximation by including a tensor expression for temperature. MHD gives an upper cutoff wavenumber for Weibel growth identical to kinetic theory if the electron drift velocity is included in the evolution of the electron temperature. MHD overestimates growth rates compared to kinetic theory and gives maximum growth at larger wavenumbers, but it only leads to divergent results as temperature anisotropy tends to infinity. Thermal conduction in MHD lowers the growth rates and shifts maximum growth to smaller wavenumbers, but flux-limited or nonlocal thermal conduction is found to significantly limit these effects. The results are used as a first step in evaluating the potential of MHD to simulate expansion-driven Weibel instability, which occurs due to the temperature tensor only cooling in the direction of expansion. Finally, the small temperature anisotropy found in expansion-driven Weibel instability means that MHD could be an adequate model because growth rates are much less than the electron plasma frequency, and wavelengths are much greater than the Debye length.

Davies, J. R. [Univ. of Rochester, NY (United Stat↗

The effects of Venus' thermal structure on buoyant magma ascent

The recent Magellan images have revealed a broad spatial distribution of surface volcanism on Venus. Previous work in modeling the ascent of magma on both Venus and Earth has indicated that the planetary thermal structure significantly influences the magmatic cooling rates and thus the amount of magma that can be transported to the surface before solidification. In order to understand which aspects of the thermal structure have the greatest influence on the cooling of ascending magma, we have constructed magma cooling curves for both plutonic and crack buoyant ascent mechanisms, and evaluated the curves for variations in the planetary mantle temperature, thermal gradient curvature with depth, surface temperature gradient, and surface temperature. The planetary thermal structure is modeled as T/T(sub 0) = 1-tau(1-Z/Z(sub 0)(exp n), where T is the temperature, T(sub 0) is the source depth temperature, tau = 1-(T(sub s)/T(sub 0)) where T(sub s) is the planetary surface temperature, Z is the depth, Z(sub 0) is the source depth, and n is a constant that controls thermal gradient curvature with depth. The equation is used both for mathematical convenience and flexibility, as well as its fit to the thermal gradients predicted by the cooling half-space models. We assume a constant velocity buoyant ascent, body-averaged magma temperatures and properties, an initially crystal-free magma, and the same liquidus and solidus for both Venus and Earth.

Sakimoto, S. E. H.↗

FUN3D Mesh Adaptation Simulations of 4th AIAA High-Lift Prediction Workshop Case 1a

This paper describes in detail the FUN3D Mesh adaptation simulation results for Case 1a of the 4th AIAA High-Lift Prediction Workshop (HLPW-4). The Stabilized Finite Element library with FUN3D is used along with the Negative Spalart-Allmaras One-Equation Model to simulate fully turbulent free-air flow around the High Lift version of the Common Research Model at three flap settings. Meshes are adapted using multi-scale and goal-oriented metrics. The results are presented and the discussion focuses on the predicted differences between the three flap settings. Coefficients of lift, drag, and pitching moment are predicted within 7% of values measured in the QinetiQ 5m Wind Tunnel at each flap setting. Predicted surface pressures and skin friction distributions show inconsistent discrepancies with measurements across the flap settings. Consequently, flap increments for these figures of merit are not well predicted. These predictions are consistent with submissions to the HLPW-4 from other RANS-SA methods. Visualizations of Liutex vortex cores colored by angular velocity are presented as possible new quantitative and qualitative perspective for solution evaluation.

High-Lift Prediction Workshop↗

Effects of geometry and jet velocity on noise associated with an upper-surface-blowing model

The noise characteristics associated with various upper surface blowing configurations were investigated using a small model consisting of a plate and flap assembly (simulated wing with flap) attached to a rectangular nozzle. Nozzle aspect ratio, flow-run length, and flap-deflection angle were the experimental parameters studied. Three nozzle-exit velocities were used. The normalized noise spectra obtained for different nozzle aspect ratios proved to be similar in terms of Strouhal number based on jet velocity and flow-run length. Consequently, the need for knowing local flow velocity and length scales (for example, at the flap trailing edge) as required in some of the existing noise prediction schemes is eliminated. Data are compared with results computed from three different noise prediction schemes, and the validity of each scheme is assessed. A simple method is proposed to evaluate the frequency dependence of acoustic shielding obtained with the simulated wing flap.

Clark, L. R.↗

Mars Observer Propulsion and Pyrotechnics Corrective Actions Test Program Status-1999

An extensive propulsion and pyrotechnic test program has been in progress at the NASA White Sands Test Facility since 1995. This program created the capabilities to: accurately measure and characterize pyrovalve combustion product blow-by into propellant systems; characterize valve operation using a Velocity Interferometer System for Any Reflector (VISAR); and evaluate hydrazine and monomethylhydrazine thermal decomposition initiated by blow-by. These capabilities were further utilized and refined this year. Low blow-by pyrovalves manufactured by Conax Florida Corporation continued to be evaluated as a potential corrective measure for blow-by induced propellant explosions. Development and testing of various advanced pyrovalves and investigation of explosion mechanisms also continued. Current and near-term testing includes: evaluation of 3/8 in. Conax pyrovalves and other commercially available valves; development and testing of advanced pyrovalve subcomponent technologies including a zero blow-by pyrovalve ram, composite overwrapped ram cylinder, and a zero particulate generating poppet; investigation of non-destructive evaluation techniques to evaluate pyrovalve ram seals; and testing and modeling of pyrotechnically induced explosive hydrazine decomposition. Evaluation of 3/8 in. Conax valves will include operational margin testing to be accomplished at NASA Langley Research Center. The test program also seeks to compile and format significant amounts of data from this and other pyrovalve test programs to generate a pyrovalve applications handbook. The handbook will facilitate formation of standards that ensure safe spacecraft applications. Current data and future plans are discussed, and community interaction is encouraged.

Saulsberry, Regor↗

Theoretical evaluation of a V/STOL fighter model utilizing the PAN AIR code

The PAN AIR computer code was investigated as a tool for predicting closely coupled aerodynamic and propulsive flowfields of arbitrary configurations. The NASA/Ames V/STOL fighter model, a configuration of complex geometry, was analyzed with the PAN AIR code. A successful solution for this configuration was obtained when the nozzle exit was treated as an impermeable surface and no wakes were included around the nozzle exit. When separated flow was simulated from the end of the nacelle, requiring the use of wake networks emanating from the nozzle exit, a number of problems were encountered. A circular body nacelle model was used to investigate various techniques for simulating the exhaust plume in PAN AIR. Several approaches were tested and eliminated because they could not correctly simulate the interference effects. Only one plume modeling technique gave good results. A PAN AIR computation that used a plume shape and inflow velocities obtained from the Navier-Stokes solution for the plume produced results for the effects of power that compared well with experimental data.

Howell, G. A.↗

Solving the Orszag–Tang vortex magnetohydrodynamics problem with physics-constrained convolutional neural networks

We study the 2D Orszag–Tang vortex magnetohydrodynamics (MHD) problem through the use of physics-constrained convolutional neural networks (PCNNs) for forecasting the density, ρ, and the magnetic field, B, as well as the prediction of B given the velocity field v of the fluid. In addition to translation equivariance from the convolutional architecture, other physics constraints were embedded: absence of magnetic monopoles, non-negativity of ρ, use of only relevant variables, and the periodic boundary conditions of the problem. The use of only relevant variables and the hard constraint of non-negative ρ were found to facilitate learning greatly. The divergenceless condition ∇·B=0 was implemented as a hard constraint up to machine precision through the use of a magnetic potential to define B=∇×A. Residual networks and data augmentation were also used to improve performance. This allowed for some of the residual models to function as surrogate models and provide reasonably accurate simulations. For the prediction task, the PCNNs were evaluated against a physics-informed neural network, which had the ideal MHD induction equation as a soft constraint. Several models were able to generate highly accurate fields, which are visually almost indistinguishable and have low mean squared error. Only methods with built-in hard constraints produced physical fields with ∇·B=0. The use of PCNNs for MHD has the potential to produce physically consistent real-time simulations to serve as virtual diagnostics in cases where inferences must be made with limited observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗

Mesh-Based Entry Vehicle and Explosive Debris Re-Contact Probability Modeling

The risk to a crewed vehicle arising from potential re-contact with fragments from an explosive breakup of any jettisoned spacecraft segments during entry has long sought to be quantified. However, great difficulty lies in efficiently capturing the potential locations of each fragment and their collective threat to the vehicle. The method presented in this paper addresses this problem by using a stochastic approach that discretizes simulated debris pieces into volumetric cells, and then assesses strike probabilities accordingly. Combining spatial debris density and relative velocity between the debris and the entry vehicle, the strike probability can be calculated from the integral of the debris flux inside each cell over time. Using this technique it is possible to assess the risk to an entry vehicle along an entire trajectory as it separates from the jettisoned segment. By decoupling the fragment trajectories from that of the entry vehicle, multiple potential separation maneuvers can then be evaluated rapidly to provide an assessment of the best strategy to mitigate the re-contact risk.

McPherson, Mark A.↗