Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical simulation”

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 1,099 records · Page 61

Direct Numerical Simulation of Acoustic Noise Generation from the Nozzle Wall of a Hypersonic Wind Tunnel

The acoustic radiation from the turbulent boundary layer on the nozzle wall of a Mach 6 Ludwieg Tube is simulated using Direct Numerical Simulations (DNS), with the flow conditions falling within the operational range of the Mach 6 Hypersonic Ludwieg Tube, Braunschweig (HLB). The mean and turbulence statistics of the nozzle-wall boundary layer show good agreement with those predicted by Pate's correlation and Reynolds Averaged Navier-Stokes (RANS) computations. The rms pressure fluctuation P'(rms)/T(w) plateaus in the freestream core of the nozzle. The intensity of the freestream noise within the nozzle is approximately 20% higher than that radiated from a single at pate with a similar freestream Mach number, potentially because of the contributions to the acoustic radiation from multiple azimuthal segments of the nozzle wall.

Huang, Junji↗

Construction of Supersonic Flight Dynamics Test (SFDT-2) Vehicle Monte Carlo Splashdown Footprints for Use in Range Safety and Recovery Operations

The Low-Density Supersonic Dececlerator (LDSD) project performed the second test flight of a supersonic inflatable device (SIAD) and ring-sail parachute on June 8, 2015 over the Pacific Ocean west of Kauai. In order for recovery ships to quickly extract the test vehicle hardware from the ocean, and in the interest of safety to the population of the islands of Kauai and Nihau, statistical estimates of the splashdown location had to be performed. This paper describes the modeling assumptions used to generate the Monti Carlo splashdown trajectory simulation results and the use of the splashdown probability ellipses in support of satisfying range safety requirements.

Strauss, William D.↗

Photon Counting Detector Array Algorithms for Deep Space Optical Communications

For deep-space optical communications systems utilizing an uplink optical beacon, a single-photon-counting detector array on the flight terminal can be used to simultaneously perform uplink tracking and communications as well as accurate downlink pointing at photon-starved (pW/sq. m) power levels. In this paper, we discuss concepts and algorithms for uplink signal acquisition, tracking, and parameter estimation using a photon-counting camera. Statistical models of detector output data and signal processing algorithms are presented, incorporating realistic effects such as Earth background and detector/readout blocking. Analysis and simulation results are validated against measured laboratory data using state-of-the-art commercial photon-counting detector arrays, demonstrating sub-microradian tracking errors under channel conditions representative of deep space optical links.

beacon tracking↗

Practical Aspects of Multiple-Input Design for Aircraft System Identification Flight Tests

Practical aspects of multiple-input design methods for aircraft system identification flight tests are discussed. The orthogonal optimized multisine multiple-input design technique is explained and demonstrated, along with useful variations. Practical implementation and application issues associated with applying orthogonal optimized multisine inputs in flight are explained, and solutions are demonstrated using flight data from full-scale and subscale aircraft, as well as simulated flight data from an F-16 nonlinear simulation. Topics include truncated and extended excitation inputs, variations for control effectors with limited movement capabilities, aircraft dynamic response to excitation inputs, pilot implementation, statistical uncertainty quantification, multiple-input efficiency and effectiveness, flight testing with active feedback control, frequency-response estimation, control actuator modeling, and control interaction effects. Examples demonstrate that the orthogonal optimized multisine input design method is a practical, versatile, and effective approach for designing multiple excitation inputs for aircraft system identification flight tests.

Multiple-Input Design↗

Man overboard - What next?

A study was performed in the Virtual Interactive Environment Workstation (VIEW) at NASA Ames Research Center. Simulations were conducted to assess the feasibility and quantify the fuel and time reguirements for a stranded crewperson to return himself to a space station after an accidental separation. A hand-held thruster, similar to the Hand-Held Maneuvering Unit from the Gemini Program was used for propulsion. Thirty different separation scenarios were composed of three separation rates, five initial spin rates, and ability to use an attitude hold mode in a repeated measures design. Statistically significant results were produced by separation velocity. Fuel, maximum range, time to maximum range, maximum axial range, and final axial velocity increased with separation rate. VIEW was determined to be a useful device for simulating accidental separations, and a hand-held thruster is a viable alternative for accomplishing a self rescue.

Brody, Adam R.↗

The Vertical Error Characteristics of GOES-derived Winds: Description and Impact on Numerical Weather Prediction

Errors in the height assignment of some satellite-derived winds exist because the satellites sense radiation emitted from a finite layer of the atmosphere rather than a specific level. Potential problems in data assimilation may arise because the motion of a measured layer is often represented by a single-level value. In this research, cloud and water vapor motion winds that are derived from the Geostationary Operational Environmental Satellites (GOES winds) are compared to collocated rawinsonde observations (RAOBs). An important aspect of this work is that in addition to comparisons at each assigned height, the GOES winds are compared to the entire profile of the collocated RAOB data to determine the vertical error characteristics of the GOES winds. The impact of these results on numerical weather prediction is then investigated. The comparisons at individual vector height assignments indicate that the error of the GOES winds range from approx. 3 to 10 m/s and generally increase with height. However, if taken as a percentage of the total wind speed, accuracy is better at upper levels. As expected, comparisons with the entire profile of the collocated RAOBs indicate that clear-air water vapor winds represent deeper layers than do either infrared or water vapor cloud-tracked winds. This is because in cloud-free regions the signal from water vapor features may result from emittance over a thicker layer. To further investigate characteristics of the clear-air water vapor winds, they are stratified into two categories that are dependent on the depth of the layer represented by the vector. It is found that if the vertical gradient of moisture is smooth and uniform from near the height assignment upwards, the clear-air water vapor wind tends to represent a relatively deep layer. The information from the comparisons is then used in numerical model simulations of two separate events to determine the forecast impacts. Four simulations are performed for each case: 1) A control simulation that assimilates no satellite wind data, 2) assimilation of all GOES winds according to their assigned single level height, 3) assimilation of all GOES winds spread over multiple levels, and 4) assimilation of all GOES winds spread over multiple levels, but with variations in the vertical influence of clear-air water vapor winds based on the moisture profile in the model. In the first case, a strong mid-latitude cyclone is present and the use of the satellite data results in improved storm tracks during the initial approx. 36 h forecast period. This is because the satellite data improves the analysis of the environment into which the storm progresses. Statistics for mean wind vector and height differences show that, with the exception of the height field at later times in the first case, the use of GOES winds improves the simulation with time. The simulation results suggest that it is beneficial to spread the GOES wind information over multiple levels, particularly when the moisture profile is used to define the vertical influence.

Rao, P. Anil↗

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗

Continued Research into Characterizing the Preturbulence Environment for Sensor Development, New Hazard Algorithms and Experimental Flight Planning

The purpose of the research was to develop and test improved hazard algorithms that could result in the development of sensors that are better able to anticipate potentially severe atmospheric turbulence, which affects aircraft safety. The research focused on employing numerical simulation models to develop improved algorithms for the prediction of aviation turbulence. This involved producing both research simulations and real-time simulations of environments predisposed to moderate and severe aviation turbulence. The research resulted in the following fundamental advancements toward the aforementioned goal: 1) very high resolution simulations of turbulent environments indicated how predictive hazard indices could be improved resulting in a candidate hazard index that indicated the potential for improvement over existing operational indices, 2) a real-time turbulence hazard numerical modeling system was improved by correcting deficiencies in its simulation of moist convection and 3) the same real-time predictive system was tested by running the code twice daily and the hazard prediction indices updated and improved. Additionally, a simple validation study was undertaken to determine how well a real time hazard predictive index performed when compared to commercial pilot observations of aviation turbulence. Simple statistical analyses were performed in this validation study indicating potential skill in employing the hazard prediction index to predict regions of varying intensities of aviation turbulence. Data sets from a research numerical model where provided to NASA for use in a large eddy simulation numerical model. A NASA contractor report and several refereed journal articles where prepared and submitted for publication during the course of this research.

Kaplan, Michael L.↗

The Statistical Emulators of GGCMI Phase 2: Responses of Year-to-Year Variation of Crop Yield to CO 2 , Temperature, Water and Nitrogen Perturbations

Understanding the impact of climate change on year-to-year variation of crop yield is critical to global food stability and security. While crop model emulators are believed to be lightweight tools to replace the models, few emulators have been developed to capture such interannual variation of crop yield in response to climate variability. In this study, we developed a statistical emulator with a machine learning algorithm to reproduce the response of year-to-year variation of four crop yields to CO 2 (C), temperature (T), water (W), and nitrogen (N) perturbations defined in the Global Gridded Crop Model Intercomparison Project (GGCMI) phase 2. The emulators were able to explain more than 52 % of the variance of simulated yield and performed well in capturing the year-to-year variation of global average and gridded crop yield over current croplands in the baseline. With the changes in CO 2 –temperature–water–nitrogen (CTWN) perturbations, the emulators could reproduce the year-to-year variation of crop yield well over most current cropland. The variation of R and the mean absolute error was small under the single CTWN perturbations and dual-factor perturbations. These emulators thus provide statistical response surfaces of yield, including both its mean and interannual variability, to climate factors. They could facilitate spatiotemporal downscaling of crop model simulation, projecting the changes in crop yield variability in the future and serving as a lightweight tool for multi-model ensemble simulation. The emulators enhanced the flexibility of crop yield estimates and expanded the application of large-ensemble simulations of crop yield under climate change.

Climate change↗

Path Profiles of Cn2 Derived from Radiometer Temperature Measurements and Geometrical Ray Tracing

Atmospheric turbulence has significant impairments on the operation of Free-Space Optical (FSO) communication systems, in particular temporal and spatial intensity fluctuations at the receiving aperture resulting in power surges and fades, changes in angle of arrival, spatial coherence degradation, etc. The refractive index structure parameter Cn2 is a statistical measure of the strength of turbulence in the atmosphere and is highly dependent upon vertical height. Therefore to understand atmospheric turbulence effects on vertical FSO communication links such as space-to-ground links, it is necessary to specify Cn2 profiles along the atmospheric propagation path. To avoid the limitations on the applicability of classical approaches, propagation simulation through geometrical ray tracing is applied. This is achieved by considering the atmosphere along the optical propagation path as a spatial distribution of spherical bubbles with varying relative refractive index deviations representing turbulent eddies. The relative deviations of the refractive index are statistically determined from altitude-dependent and time-varying temperature fluctuations, as measured by a microwave profiling radiometer. For each representative atmosphere ray paths are analyzed using geometrical optics, which is particularly advantageous in situations of strong turbulence where there is severe wavefront distortion and discontinuity. The refractive index structure parameter is then determined as a function of height and time.

Atmospheric turbulence↗

Statistical Evaluation and Improvement of Methods for Combining Random and Harmonic Loads

Structures in many environments experience both random and harmonic excitation. A variety of closed-form techniques has been used in the aerospace industry to combine the loads resulting from the two sources. The resulting combined loads are then used to design for both yield/ultimate strength and high- cycle fatigue capability. This Technical Publication examines the cumulative distribution percentiles obtained using each method by integrating the joint probability density function of the sine and random components. A new Microsoft Excel spreadsheet macro that links with the software program Mathematica to calculate the combined value corresponding to any desired percentile is then presented along with a curve tit to this value. Another Excel macro that calculates the combination using Monte Carlo simulation is shown. Unlike the traditional techniques. these methods quantify the calculated load value with a consistent percentile. Using either of the presented methods can be extremely valuable in probabilistic design, which requires a statistical characterization of the loading. Additionally, since the CDF at high probability levels is very flat, the design value is extremely sensitive to the predetermined percentile; therefore, applying the new techniques can substantially lower the design loading without losing any of the identified structural reliability.

Brown, A. M.↗

Use of the tilt cue in a simulated heading tracking task

The task was performed with subjects using visual-only cues and combined visual and roll-axis motion cues. Half of the experimental trials were conducted with the simulator rotating about the horizontal axis; to suppress the tilt cue, the remaining trials were conducted with the simulator cab tilted 90 deg so that roll-axis motions were about earth vertical. The presence of the tilt cue allowed a substantial and statistically significant reduction in performance scores. When the tilt cue was suppressed, the availability of motion cues did not result in significant performance improvement. These effects were accounted for by the optimal-control pilot/vehicle model, wherein the presence or absence of various motion cues was represented by appropriate definition of the perceptual quantities assumed to be used by the human operator.

Levison, W. H.↗

Prediction of spacecraft potentials at geosynchronous orbit

Two relatively straightforward techniques are outlined for determining spacecraft potentials in the limit of a 'thick sheath' surrounding the spacecraft. A statistical model of the various features of the geosynchronous environment based on ATS-5 and ATS-6 data and an analytic model capable of detailed simulation of the low energy geosynchronous environment are also discussed. The results from these two environmental models are then combined with the charging models in order to provide estimates of the relationships between the geomagnetic index and spacecraft potential. The results are compared with actual potential measurements from ATS-5 and ATS-6.

Garrett, H. B.↗

Magnetic field directional discontinuities. I - Minimum variance errors

The paper deals with a statistical analysis of the errors associated with a minimum variance analysis of directional discontinuities by use of an idealized model of these discontinuities and various simulations, and also by an examination of actual Mariner 10 interplanetary magnetic field data. An empirical expression is derived for the magnitude of the error in an estimated discontinuity normal component, relative to the total field across the directional discontinuity. The analysis was performed primarily to aid in differentiating between interplanetary tangential and rotational discontinuities observed by Mariner 10.

Lepping, R. P.↗

Results of a joint NOAA/NASA sounder simulation study

This paper presents the results of a joint NOAA and NASA sounder simulation study in which the accuracies of atmospheric temperature profiles and surface skin temperature measuremnents retrieved from two sounders were compared: (1) the currently used IR temperature sounder HIRS2 (High-resolution Infrared Radiation Sounder 2); and (2) the recently proposed high-spectral-resolution IR sounder AMTS (Advanced Moisture and Temperature Sounder). Simulations were conducted for both clear and partial cloud conditions. Data were analyzed at NASA using a physical inversion technique and at NOAA using a statistical technique. Results show significant improvement of AMTS compared to HIRS2 for both clear and cloudy conditions. The improvements are indicated by both methods of data analysis, but the physical retrievals outperform the statistical retrievals.

Phillips, N.↗

Estimation of UH-60 blade aerodynamic loads and rotor impedance using generalized strain pattern/Kalman filter approach

The purpose of the study is to develop and verify a methodology capable of predicting the vibration levels and estimating the aerodynamic loads and rotor impedance of a rotorcraft blade. Simulated flight test data is generated, blade airloads and elastic hub motions are estimated from the simulated data through the use of the Kalman filter/smoother, simulation upgrading and parameter identification are performed, and the ability to identify rotor impedance from a simulation by isolating the rotor model and providing a prescribed motion for the hub as rotor excitation is demonstrated. It is pointed out that the statistical estimation procedure utilized in the proposed methodology minimizes the impact of sensor noise, truncation error, and instrumentation bias on the results.

Bruhis, Ofer↗

Statistics for laminar flamelet modeling

Statistical information required to support modeling of turbulent premixed combustion by laminar flamelet methods is extracted from a database of the results of Direct Numerical Simulation of turbulent flames. The simulations were carried out previously by Rutland (1989) using a pseudo-spectral code on a three dimensional mesh of 128 points in each direction. One-step Arrhenius chemistry was employed together with small heat release. A framework for the interpretation of the data is provided by the Bray-Moss-Libby model for the mean turbulent reaction rate. Probability density functions are obtained over surfaces of the constant reaction progress variable for the tangential strain rate and the principal curvature. New insights are gained which will greatly aid the development of modeling approaches.

Cant, R. S.↗

Narrowband and wideband characterisation of satellite mobile/PCN channel

This paper presents models characterizing satellite mobile channel. Statistical narrowband models based on the CSER high elevation angle channel measurement campaign are reported. Such models are understood to be useful for communication system simulations. It has been shown from the modelling results that for the mobile satellite links at high elevation angles line-of-sight (LOS) signal is available most of the time, even under the heavy shadowing conditions. Wideband measurement campaign which CSER is about to undertake, and subsequently the modelling approach to be adopted is also discussed. It is noted that a wideband channel model is expected to provide a useful tool in investigating CDMA applications.

Butt, G.↗