Search NASA⌕ Search

SEARCH · Search NASA

Results for “prediction correction”

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 343 records · Page 19

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗

Turbulence Simulations of Transonic Flows over an NACA-0012 Airfoil

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012 airfoil at different angles of attack covering pre-and post-buffet-onset regimes. The freestream Mach number is 0.75,and the Reynolds number based on the chord length is 10million. These conditions are the same as the wind-tunnel experimental conditions of McDevitt and Okuno (1985). The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre-and post-buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the span wise direction, WMLES predictions show buffeting consistent with the experiment.

Transonic Buffet↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

A New Global Regression Analysis Method for the Prediction of Wind Tunnel Model Weight Corrections

A new global regression analysis method is discussed that predicts wind tunnel model weight corrections for strain-gage balance loads during a wind tunnel test. The method determines corrections by combining "wind-on" model attitude measurements with least squares estimates of the model weight and center of gravity coordinates that are obtained from "wind-off" data points. The method treats the least squares fit of the model weight separate from the fit of the center of gravity coordinates. Therefore, it performs two fits of "wind- off" data points and uses the least squares estimator of the model weight as an input for the fit of the center of gravity coordinates. Explicit equations for the least squares estimators of the weight and center of gravity coordinates are derived that simplify the implementation of the method in the data system software of a wind tunnel. In addition, recommendations for sets of "wind-off" data points are made that take typical model support system constraints into account. Explicit equations of the confidence intervals on the model weight and center of gravity coordinates and two different error analyses of the model weight prediction are also discussed in the appendices of the paper.

strain-gage balance↗

Kalman-like estimation for static distributed systems Antenna shape from radiation measurements

This paper advances an approach to the determination of shape of static distributed systems. It also illustrates the application of the approach to the problems of surface diagnosis of large parabolic reflectors. The estimation methods developed combine in an optimal sense the information from an elliptic model of the structure and from measurements of the structural deflection and of the far-field pattern changes due to the structural deformation. The estimators have a predictor-corrector structure, quite similar to that of a Kalman filter. The system model is first used to obtain a predicted estimate. A correction term is then added to the prediction to obtain the final state estimate. The relative weighting between prediction and correction terms is determined by an estimator gain. As in a Kalman filter, the estimator gain can be expressed in terms of the state estimation error covariance.

Rodriguez, G.↗

Assessment of Turbulent Shock-Boundary Layer Interaction Computations Using the OVERFLOW Code

The performance of two popular turbulence models, the Spalart-Allmaras model and Menter s SST model, and one relatively new model, Olsen & Coakley s Lag model, are evaluated using the OVERFLOWcode. Turbulent shock-boundary layer interaction predictions are evaluated with three different experimental datasets: a series of 2D compression ramps at Mach 2.87, a series of 2D compression ramps at Mach 2.94, and an axisymmetric coneflare at Mach 11. The experimental datasets include flows with no separation, moderate separation, and significant separation, and use several different experimental measurement techniques (including laser doppler velocimetry (LDV), pitot-probe measurement, inclined hot-wire probe measurement, preston tube skin friction measurement, and surface pressure measurement). Additionally, the OVERFLOW solutions are compared to the solutions of a second CFD code, DPLR. The predictions for weak shock-boundary layer interactions are in reasonable agreement with the experimental data. For strong shock-boundary layer interactions, all of the turbulence models overpredict the separation size and fail to predict the correct skin friction recovery distribution. In most cases, surface pressure predictions show too much upstream influence, however including the tunnel side-wall boundary layers in the computation improves the separation predictions.

Oliver, A. B.↗

Three-Dimensional, Space-Dependent Mesoscale Diffusivity: Derivation and Implications

Recently, we presented a parameterization of an arbitrary tracer 3D mesoscale flux that describes both diabatic and adiabatic regimes without using arbitrary tapering functions. However, we did not parameterize the mesoscale diffusivity, which is the subject of this work. A key difference between the present and previous diffusivity parameterizations is that in the latter, the two main ingredients, mesoscale drift velocity and eddy kinetic energy, were not parameterized but determined using present data, which deprives the models of predictive power. Since winds, stratification, etc., are predicted to change in the future, use of these parameterizations to study future climate scenarios becomes questionable. In this work, we parameterize drift velocity and eddy kinetic energy (vertical–horizontal components), which we first assess with data [WOCE, TOPEX/Poseidon (T/P), and North Atlantic Tracer Release Experiment (NATRE)] and then use in a coarse-resolution stand-alone ocean code under Coordinated Ocean-Ice Reference Experiment I (CORE-I) forcing. We present results for the global ocean temperature and salinity, Atlantic overturning circulation, meridional heat transport, and Drake Passage transport, which we compare with several previous studies. The temperature drift is less than that of five of seven previous OGCMs, and the salinity drift is among the smallest in those studies. The predicted winter Antarctic Circumpolar Current mixed layer depths (MLDs) are in good agreement with the data. Predicting the correct MLD is important in climate studies since models that predict very deep mixed layers transfer more of the radiative perturbation to the deep ocean, reducing surface warming (and vice versa).

Eddies; Mesoscale processes; Mixing; Ocean models;↗

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗

Correcting Distortions In Optical Correlators

Coordinate transformation maps object space to correlation space. Theory developed to predict and correct typical anamorphic coordinate-transformation errors in off-axis Vander Lugt optical correlators.

Chrien, Thomas G.↗

Detecting Controller Malfunctions in Electromagnetic Environments: Modeling and Estimation of Nominal System Function - Part 1

A strategy for detecting control law calculation errors in critical flight control computers during laboratory validation testing is presented. This paper addresses Part I of the detection strategy which involves the use of modeling of the aircraft control laws and the design of Kalman filters to predict the correct control commands. Part II of the strategy which involves the use of the predicted control commands to detect control command errors is presented in the companion paper.

Weinstein, Bernice↗

Transonic flutter study of a wind-tunnel model of an arrow-wing supersonic transport

A 1/20-size, low-speed flutter model of the SCAT-15F complete airplane was tested on cables to simulate a near free-flying condition. Only the model wing and fuselage were flexible. Flutter boundaries were measured for a nominal configuration and a configuration with wing fins removed at Mach numbers M from 0.76 to 1.2. For both configurations, the transonic dip in the wing flutter dynamic pressure q boundary was relatively small and the minimum flutter q occurred near M = 0.92. Removing the wing fins increased the flutter q about 14 percent and changed the flutter mode from symmetric to antisymmetric. Vibration and flutter analyses were made using a finite-element structural representation and subsonic kernel-function aerodynamics. For the nominal configuration, the analysis (using calculated modal data) predicted the experimental flutter q levels within 10 percent but did not predict the correct flutter mode at the higher M. For the configuration without wing fins, the analysis predicted 16 to 36 percent unconservative (higher than experimental) flutter q levels and showed extreme sensitivity to mass representation details that affected wing tip mode shapes. For high subsonic M, empennage aerodynamics had a significant effect on the predicted flutter boundaries of several symmetric modes.

Ruhlin, C. L.↗

Transonic flutter study of a wind-tunnel model of an arrow-wing supersonic transport

A 1/20-size, low-speed flutter model of the SCAT-15F complete airplane was tested in a wind-tunnel on cables to simulate a near free-flying condition. Only the model wing and fuselage were flexible. Flutter boundaries were measured for a nominal configuration and a configuration with wing fins removed at Mach numbers M from 0.76 to 1.2. For both configurations, the transonic dip in the wing flutter dynamic pressure q boundary was relatively small and the minimum flutter q occurred near M = 0.92. Removing the wing fins increased the flutter q about 14 percent and changed the flutter mode from symmetric to antisymmetric. Vibration and flutter analyses were made using a finite-element structural representation and subsonic kernel-function aerodynamics. For the nominal configuration, the analysis (using calculated modal data) predicted the experimental flutter q levels within 10 percent but did not predict the correct flutter mode at the higher M. For the configuration without wing fins, the analysis predicted 16 to 36 percent unconservative (higher than experimental) flutter q levels and showed extreme sensitivity to mass shapes.

Ruhlin, C. L.↗

An improved far field drag calculation method for nonlinear CFD codes

An improved method to calculate drag based on far field conditions is presented and demonstrated. The method is illustrated using two CFD codes, PARC and CFL3D, and examining their ability to preserve lift and drag coefficients along grid line defined integration paths. The flow fields were generated by solving the Euler equations for a NACA 0012 airfoil at a free stream Mach number of 0.8 and angle of attack of 1.25 degrees. In comparison to force coefficients obtained by surface pressure integration, neither code acceptably preserved both force coefficients throughout the near and far fields. An investigation into the relationship between numerical prediction error and calculated force coefficients revealed a direct connection between solution mass conservation error and force coefficients error. A method to correct the predicted force coefficients based on integrated mass conservation error is described and demonstrated. The corrected force coefficients for both codes are shown to be more accurate than the uncorrected values. The correction method is applicable to both two and three dimensions and is independent of the algorithm used to generate the flow field.

Sirbaugh, James R.↗

Prediction of transition on a flat plate under the influence of free-stream turbulence using low-Reynolds-number two-equation turbulence models

The prediction characteristics and capabilities of two popular two-equation low-Reynolds-number turbulence models (Jones, Launder 1972; Lam, Bremhorst 1981) have been evaluated with respect to the prediction of transition of a flat plate under the influence of free-stream turbulence. The sensitivity of the predictions to free-stream turbulence intensity, initial starting location of the calculation, and the assumed initial starting profiles for k and epsilon has been determined and presented. Although both models predict the correct qualitative characteristics of transition, they also exhibit significant quantitative deficiencies with regard to both the predicted location and the length of transition. A modification to the production term in the turbulent kinetic energy equation is proposed which is based on a simple stability criterion and correlated to the free-stream turbulence level. The modification becomes inactive in the fully turbulent regime, but is shown to improve both the qualitative and quantitative characteristics of the transition predictions.

Schmidt, Rodney C.↗

Dipole gravitational radiation in Rosen's theory of gravity - Observable effects in the binary system PSR 1913+16

It is shown that Rosen's (1973) bimetric theory of gravity predicts the emission of dipole gravitational radiation from binary systems containing neutron stars, such as the binary pulsar PSR 1913+16, which causes rapid changes in orbital period. The theory also predicts sizable corrections to masses inferred from orbital data and periastron-shift data. It is demonstrated that this prediction is inconsistent with the observed upper limit on period changes unless the system consists of two neutron stars whose masses differ by less than 0.3 solar mass, or a neutron star of mass less than 0.4 solar mass and a companion which must be a rapidly rotating white dwarf or a helium main-sequence star. Because Rosen's theory is in agreement with all solar-system experiments to date, this represents a feasible test of its viability.

Will, C. M.↗

Updates to the NuMI Flux Simulation at MicroBooNE (V.1.0)

The NuMI flux prediction being used by all experiments at Fermilab that are sensitive to the NuMI beam is based on Geant v4.9.2.03 (G4.9) and uses the default FTFP-BERT physics list with no custom cross sections. In order to correct these predictions to match world data, PPFX, the Package to Predict the FluX, is used to reweight the flux prediction. In kinematic regions of overlap between the prediction and data, mainly NA49 measurements are used to correct for various hadron production processes. Specifcally, PPFX looks at the entire ancestry chain of the neutrino, up to the incident 120 GeV proton beam, and applies corrections for various hadron interactions based on said data. These corrections are implemented as a function of the Feynman-x (x F ) and transverse momentum (p T ) of the incident hadron.

43 PARTICLE ACCELERATORS↗

Polarization in massive X-ray binaries. I - A low-inclination model for Cygnus X-1

The possibility that variable linear polarization in massive X-ray binaries is produced by electron scattering in an asymmetric stellar wind is investigated. The stellar wind is asymmetric because of the gravitational field of the secondary (X-ray source). The degree of asymmetry and the magnitude of the linear polarization are controlled by the degree to which the primary star fills its Roche lobe. For the well-observed X-ray binary Cyg X-1, the present model can produce the correct magnitude for the polarization. Provided that the inclination of the system is less than about 20 deg, the present model should also predict the correct phase dependence of the polarization. Modifications to the model are described which would enable it to apply to systems with higher inclination.

Friend, D. B.↗