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 901 records · Page 50

Computational Analysis of Water Film Thickness During Rain Events for Assessing Hydroplaning Risk Part 2: Rough Road Surfaces

Hydroplaning occurs when a water film forms on a road during a rain event and vehicles are traveling at a speed that does not provide sufficient time for the tires to push the water film out of the tire path. Under these conditions, the tire loses contact with the road and the driver may lose control over the vehicle. The present study is a part of a multiyear research effort aimed at gaining a better understanding of the factors that contribute to hydroplaning risk in order to minimize the occurrence of hydroplaning accidents. Three-dimensional computational fluid dynamics (CFD) analysis was used to investigate the water film thickness on multilane roads, one of the key parameters in evaluating hydroplaning risk. In the first part of the study, water films forming on wide roadways, i.e. multilane highways with 2, 4, and 6 lanes per side, were analyzed, with varying cross slope, longitudinal slope, and rainfall rate. Roadways with and without curbs and drainage were included in the analysis. The analysis in the first part was limited to nearly smooth asphalt or concrete surfaces because the maximum roughness height that can be specified in a CFD model using wall functions to determine shear stress at the road surface boundary is one half the thickness of computational cell layer adjacent to the boundary and the need to resolve the vertical velocity distribution in a thin water film requires cell layers that are a fraction of a millimeter thick. In this part of the study, the analysis was extended to consider the influence of pavement roughness over a range extending up 3.6 mm on the development of water film thickness. Various methods of modeling roughness were investigated and are presented: the roughness height model, a meshed-out geometry of the macrotexture of the pavement, and a porous region model. The computational results were analyzed and compared with experimental measurements performed by Gallaway et al.

42 ENGINEERING↗

Simulation evaluation of transition and hover flying qualities of a mixed-flow, remote-lift STOVL aircraft

Using a generalized simulation model developed for piloted evaluations of STOVL aircraft, an initial fixed-base simulation of a mixed-flow, remote-lift configuration has been completed. Objectives were to evaluate the integration of the aircraft's flight and propulsion controls to achieve good flying qualities throughout the low-speed flight envelope; to determine control power used during transition, hover, and vertical landing; and to evaluate the transition flight envelope considering the influence of thrust deflection of the remote-lift component. Pilots' evaluations indicated that Level 1 flying qualities could be achieved for deceleration to hover in instrument conditions, for airfield landings, and for recovery to a small ship when attitude and velocity stabilization and command augmentation control modes were provided. Level 2 flying qualities were obtained for these same tasks when only the attitude command mode was used, leaving the pilot to perform the task of thrust management required to control the flight-path and speed in transition and the horizontal and vertical translational velocities in hover. Thrust margins were defined for vertical landing as a function of ground effect and hot-gas ingestion.

Franklin, James A.↗

Investigation at Mach Numbers of 0.60 to 3.50 of Blended Wing-Body Combinations with Cambered and Twisted Wings with Diamond, Delta and Arrow Plan Forms

This investigation is a continuation of the experimental and theoretical evaluation of blended wing-body combinations. The basic diamond, delta, and arrow plan forms which had an aspect ratio of 2 with leading-edge sweeps of 45.00 deg., 59.04 deg., and 70.82 deg. and trailing edge of -45.00 deg., -18.43 deg., and 41.19 deg., respectively, are used herein as standards for evaluating the effects of camber and warp. The wing thickness distributions were computed by varying the section shape along with the body radii (blending process) to match the prescribed area distribution and wing plan form. The wing camber and warp were computed to try to obtain nearly elliptical spanwise and chordwise load distributions for each plan form and thus to obtain low drag due to lift for a range of Mach numbers for which the velocities normal to the wing leading edge are subsonic. Elliptical chordwise load distributions were not possible for the plan forms and design conditions selected, so these distributions were somewhat different for each plan form. The models were tested with transition fixed at Mach numbers from 0.60 to 3.50 and at Reynolds numbers, based on the mean aerodynamic chord of the wing, of roughly 4,000,000 to 9,000,000. At speeds where the velocities normal to the wing leading edges were supersonic, an increase in the experimental wave-drag coefficients due to camber and twist was evident, but this penalty decreased with increased sweep. Thus the minimum wave-drag coefficients for the cambered arrow model were almost identical with the zero-lift wave- drag coefficients for the uncambered arrow model at all test Mach numbers.

Holdaway, George H.↗

Isolating Curvature Effects in Computing Wall-Bounded Turbulent Flows

The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

Rumsey, Christopher L.↗

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗

Analysis on Evaluations of Monterey Bay Aquarium Research Institute’s Wave Energy Converter’s Field Data Using WEC-Sim and Gazebo: A Simulation Tool Comparison

Although many studies have validated wave energy converter (WEC) numerical models against scaled prototype experimental data, there remains a notable lack of validation using data from full-scale deployed WECs. This paper compares two numerical models of Monterey Bay Aquarium Research Institute’s Wave Energy Converter (MBARI-WEC), a two-body point absorber with an electro-hydraulic power take-off system (PTO). The models are implemented in WEC-Sim/Simscape and Gazebo Simulator. A statistical analysis of the models was performed, and field results were obtained to compare the models’ accuracy in predicting the RMS piston velocity, RMS motor speed, and mean electric power compared to field data for 56 observations across varying sea states. The Gazebo model demonstrated a closer agreement across all three parameters for a majority of the observations. When compared to the field data, the Gazebo and WEC-Sim models exhibited average mean electric power overestimations of 13% and 22%, respectively.

16 TIDAL AND WAVE POWER↗

Sensitivity analysis of numerical modeling input parameters on floating offshore wind turbine loads

Abstract. Floating wind turbines must withstand a unique and challenging set of loads from the wind and ocean environment. To de-risk development, accurate predictions of these loads are necessary. Uncertainty in modeling predictions leads to larger required safety factors, increasing project costs and the levelized cost of energy. Complex aero-hydro-elastic modeling tools use many input parameters to represent the wind, waves, current, aerodynamic loads, hydrodynamic loads, and structural properties. It is helpful to understand which of these parameters ultimately drives a design. In this work, an ultimate and fatigue-proxy load sensitivity analysis was performed with 35 different input parameters, using an elementary effects approach to identify the most influential parameters for a case study involving the National Renewable Energy Laboratory (NREL) 5 MW baseline wind turbine atop the OC4-DeepCwind semisubmersible during normal operation. The importance of each parameter was evaluated using 14 response quantities of interest across three operational wind speed conditions. The study concludes that turbulent wind velocity standard deviation is the parameter with the strongest sensitivity; this value is important not just for turbine loads, but also for the global system response. The system center of mass in the wind direction is found to have the highest impact on the system rotation and tower loads. The current velocity is found to be the most dominating parameter for the system global motion and consequently the mooring loads. All tested wind turbulence parameters in addition to the standard deviation are also found to be influential. Wave characteristics are influential for some fatigue-proxy loading but do not significantly impact the extreme ultimate loads in these operational load cases. The required number of random seeds for stochastic environmental conditions is considered to ensure that the sensitivities are due to the input parameters and not due to the seed. The required number of analysis points in the parameter space is identified so that the conclusions represent a global sensitivity. The results are specific to the platform, turbine, and choice of parameter ranges, but the demonstrated approach can be applied widely to guide focus in parameter uncertainty.

17 WIND ENERGY↗

Space shuttle post-entry and landing analysis. Volume 2: Appendices

Four candidate navigation systems for the space shuttle orbiter approach and landing phase are evaluated in detail. These include three conventional navaid systems and a single-station one-way Doppler system. In each case, a Kalman filter is assumed to be mechanized in the onboard computer, blending the navaid data with IMU and altimeter data. Filter state dimensions ranging from 6 to 24 are involved in the candidate systems. Comprehensive truth models with state dimensions ranging from 63 to 82 are formulated and used to generate detailed error budgets and sensitivity curves illustrating the effect of variations in the size of individual error sources on touchdown accuracy. The projected overall performance of each system is shown in the form of time histories of position and velocity error components.

Crawford, B. S.↗

Acoustic resonances and sound scattering by a shear layer

The energy reflection coefficient is evaluated numerically for plane waves incident on a plane shear layer having a linear velocity profile. The shear layer is found to exhibit no resonances and no Brewster angles. The behavior of the reflection coefficient depends crucially on the parameter tau, a nondimensional measure of the disturbance Strouhal number with respect to the disturbance Mach number in the mean flow direction. For moderate values of tau, the amplified reflection regime degenerates into the total reflection one, whereas in the ordinary reflection regime the variation of the reflection coefficient with tau depends on whether or not the corresponding vortex sheet has a Brewster angle. The results indicate that caution should be exercised in uncritically modeling a finite thickness shear layer by a corresponding vortex sheet.

Koutsoyannis, S. P.↗

Two-Dimensional Dynamic Simulation of a Continuous Foil Bearing

In this paper, the two dimensional(radial and circumferential) transient Navier-Stokes equations are used to solve the hydrodynamic problem in conjunction with the time dependent motion of the journal, and the deformable, spring supported foil. The elastic deformation of the foil and its supports are simulated by a finite element model. The time-dependent Navier-Stokes formulation is used to solve for the interaction between the fluid lubricant, the motion of the journal and the deformable foil boundary. The steady state, the quasi-transient and the full transient dynamic simulation of the foil-fluid journal interaction are examined on a comparative basis. For the steady state simulation, the fluid lubricant pressures are evaluated for a particular journal position, by means of an iterative scheme until convergence is achieved in both the fluid pressures and the corresponding foil deformation. For the quasi-transient case, the transient motion of the journal is calculated using a numerical integration scheme for the velocity and displacement of the journal. The deformation of the foil is evaluated through numerical iteration in feedback mode with the fluid film pressure generated by the journal motion until convergence at every time step is achieved. For the full transient simulation, a parallel real-time integration scheme is used to evaluate simultaneously the new journal position and the new deformed shape of the foil at each time step. The pressure of the fluid lubricant is iterated jointly with the corresponding journal position and the deformed foil geometry until convergence is achieved. A variable time-stepping Newmark-Beta integration procedure is used to evaluate the transient dynamics at each time step of the bearing.

Braun, M. Jack↗

SiC and Si3N4 Recession Due to SiO2 Scale Volatility Under Combustor Conditions

SiC and Si3N4 materials were tested under various turbine engine combustion environments, chosen to represent either conventional fuel-lean or fuel-rich mixtures proposed for high speed aircraft. Representative CVD, sintered, and composite materials were evaluated in both furnace and high pressure burner rig exposure. While protective SiO2 scales form in all cases, evidence is presented to support paralinear growth kinetics, i.e. parabolic growth moderated simultaneously by linear volatilization. The volatility rate is dependent on temperature, moisture content, system pressure, and gas velocity. The burner tests were used to map SiO2 volatility (and SiC recession) over a range of temperature, pressure, and velocity. The functional dependency of material recession (volatility) that emerged followed the form: exp(-QIRT) * P(exp x) * v(exp y). These empirical relations were compared to rates predicted from the thermodynamics of volatile SiO and SiO(sub x)H(sub Y) reaction products and a kinetic model of diffusion through a moving, boundary layer. For typical combustion conditions, recession of 0.2 to 2 micron/h is predicted at 1200- 1400C, far in excess of acceptable long term limits.

Smialek, James L.↗

Icing Formation Strains on Airfoil Skin

Aircraft icing has long presented a safety hazard to the aerospace industry. One problem inherent to aircraft icing is predicting the adhesive properties of the ice-aircraft interface; values commonly reported span multiple orders of magnitude, making shedding prediction and surface evaluation difficult. An effort has been undertaken at the NASA Glenn Research Center to understand the cause of the scatter in the data. Geometry dependencies, variation in testing procedures, and residual stresses in ice samples have previously been identified as significant contributors to the problem. Residual stresses were singled out for this study, which investigated the strains at the leading edge of a Common Research Model (CRM) wind tunnel model in the NASA Glenn Icing Research Tunnel. Sprays were conducted at a particle size median volumetric diameter (MVD) of 20 µm, with liquid water contents (LWC) spanning from 0.4 to 3.0 g/m3. Total temperatures were varied between 0 and -25 °C, and test section velocity between 51 and 118 m/s. The leading edge was instrumented with custom strain gages, resistance temperature detectors (RTDs), optical fibers, and a pressure mat. Results indicate Δ strains, which were the difference between recorded strains and those predicted from isotropic thermal expansion of aluminum, below -400 µε with consistent patterns of strain across a wide variety of tunnel conditions. Isotropic, linear elastic finite element analysis (FEA) was conducted by varying the temperature of simplified airfoil and ice geometries, producing similar strains to those observed in testing. This allowed induced stresses to be predicted. Predicted stresses were above the expected tensile strength for ice, and suggest a source for previously observed cracking in tunnel ice.

Ice Adhesion↗

Preliminary Result of Neoclassical Impurity Transport Affected by ICRH with a 4-D Fokker-Planck Code

To understand the role of ion cyclotron resonance heating (ICRH) in mitigating tungsten accumulation, we evaluate the impact of ICRH-induced poloidal asymmetries on neoclassical impurity transport using a newly developed 4-D Fokker–Planck code (FP4D). Unlike conventional approaches that separate RF heating and transport modelling, FP4D solves both fast-ion distribution formation and impurity transport self-consistently without bounce-averaging. Benchmarking against the NEO code confirms the accuracy of the neoclassical physics in FP4D. Using FP4D, we analyze anisotropic and non-Maxwellian minority distributions under ICRH and their effects on the equilibrium potential and tungsten transport. Results show that anisotropy induced poloidal potential variations can reduce inward tungsten flux, depending on the modelling approach. Our findings highlight the importance of resolving the full velocity and poloidal structure of minority species in predicting impurity transport under strong RF heating conditions.

Jeong, Yunho [Hanyang Univ., Seoul (Korea, Republi↗

The structure of evaporating and combusting sprays: Measurements and predictions

The structure of particle-laden jets and nonevaporating and evaporating sprays was measured in order to evaluate models of these processes. Three models are being evaluated: (1) a locally homogeneous flow model, where slip between the phases is neglected and the flow is assumed to be in local thermodynamic equilibrium; (2) a deterministic separated flow model, where slip and finite interphase transport rates are considered but effects of particle/drop dispersion by turbulence and effects of turbulence on interphase transport rates are ignored; and (3) a stochastic separated flow model, where effects of interphase slip, turbulent dispersion and turbulent fluctuations are considered using random sampling for turbulence properties in conjunction with random-walk computations for particle motion. All three models use a k-e-g turbulence model. All testing and data reduction are completed for the particle laden jets. Mean and fluctuating velocities of the continuous phase and mean mixture fraction were measured in the evaporating sprays.

Shuen, J. S.↗

Structure of nonevaporating sprays. I - Initial conditions and mean properties

Structure measurements were completed within the dilute portion of axisymmetric nonevaporating sprays. Measurements included: mean velocities, velocity fluctuations, and Reynolds stress of the gas phase; and mean velocities, fluctuating velocities, mass flux, and diameter distributions of the drop phase. The measurements were used to evaluate three typical methods of analyzing sprays: (1) a locally homogeneous flow (LHF) analysis, where slip between the phases is neglected; (2) a deterministic separated flow (DSF) analysis, where slip is considered but effects of drop interactions with turbulence are ignored; and (3) a stochastic separated flow (SSF) analysis, where effects of both slip and turbulence are considered using random-walk computations for drop motion. Measurements of initial conditions of both phases near the injector, mean gas-phase properties, and liquid flux distributions are described herein. Best agreement between predictions and measurements was obtained with the SSF model, which provided a reasonable representation of turbulent dispersion of drops. A companion paper presents additional measurements of drop continuous phase properties.

Solomon, A. S. P.↗

Dynamics of the nighttime thermosphere at Arecibo

Incoherent scatter radar observations of the nighttime F layer at Arecibo, Puerto Rico, are used to determine the O+ diffusion velocity at different altitudes. Further analysis allows evaluation of the neutral wind and the ion-drag force in the direction of the magnetic meridian. The local acceleration of the meridional neutral wind is also determined. The possibility of evaluating the viscous force using incoherent scatter radar data was also investigated. Preliminary results indicate that, at certain times and at altitudes above about 350 km, viscous drag may be a significant term in the neutral equation of motion. Combining these results allows height profiles of the meridional pressure gradient to be deduced. The pressure gradients thus derived is compared with that determined from measurements of the horizontal temperature gradient and that given by the MSIS model atmosphere.

Burnside, R. G.↗

Excitation of the low lying vibrational levels of H2O by O(3P) as measured on Spacelab 2

The data from the infrared telescope (IRT), which was flown on space shuttle Challenger Spacelab 2 mission (July 1985), were originally reported by Koch et al. (1987) as originating from near orbital emissions, primarily H2O. In this study, analysis of this data was extended to determine the collisional cross sections for the excitation of the low lying vibrational levels of H2O, present in the orbiter cloud, by atmospheric O(3P). The evaluation of the contribution to the measured signal from solar excitation and ram O excitation of outgassing H2O permits the determination of the H2O column density and the excitation cross section of the (101) level at an O(3P) velocity of approximately 7.75 km/s. Contributions to the radiation in the 1.7-3.0 micron band by transitions from the (100), (001), and multiquantum excited levels are discussed. The findings of the study are (1) the IRT data for the 4.5-9.5 micron and the nighttime data for the 1.7-3.0 micron sensors are consistent with being explained by collision excitation of H2O by O(3P), (2) diurnal variations of 4.5-9.5 micron intensities follow the model predicted O density for a full orbit, (3) daytime increases in the H2O cloud density were not evident, (4) the cross sections for the collisional excitation process are derived and compared to values computated by Johnson (1986) and Redmon et al. (1986), (5) theoretical investigation suggests greater than 60% of the radiation from H2O is a result of multiphoton emission resulting from collisional multiquanta excitation, and (6) the large daytime increase in the 1.7-3.0 micron intensity data suggests that O(+) may likely be instrumental in producing excited H2O(+) through charge exchange.

Meyerott, R. E.↗

Subsonic Analysis of 0.04-Scale F-16XL Models Using an Unstructured Euler Code

The subsonic flow field about an F-16XL airplane model configuration was investigated with an inviscid unstructured grid technique. The computed surface pressures were compared to wind-tunnel test results at Mach 0.148 for a range of angles of attack from 0 deg to 20 deg. To evaluate the effect of grid dependency on the solution, a grid study was performed in which fine, medium, and coarse grid meshes were generated. The off-surface vortical flow field was locally adapted and showed improved correlation to the wind-tunnel data when compared to the nonadapted flow field. Computational results are also compared to experimental five-hole pressure probe data. A detailed analysis of the off-body computed pressure contours, velocity vectors, and particle traces are presented and discussed.

Lessard, Wendy B.↗