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 685 records · Page 38

Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation

Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.

Swanson, R. C.↗

Direct simulations of wall-bounded compressible turbulence

A study has been initiated on the effects of compressibility (Mach number) on turbulent boundary layers. An understanding of both qualitative (turbulence structures, physical processes) and quantitative (turbulence statistics) effects are desired. This understanding should lead to better turbulence models for applications involving supersonic wall-bounded flows. Direct numerical simulations of an idealized problem will be used to accomplish these objectives. Among several possibilities we chose plane Couette flow with constant-temperature walls as the first problem to be studied. The lack of mean streamwise pressure gradient plus isothermal walls implies that both horizontal directions can be assumed to be homogeneous and that the flow can reach a statistically steady state. Together, these features greatly simplify the calculations and analyses of the results. To date, an algorithm has been developed and implemented (but not yet fully tested) for the accurate solution of the Navier-Stokes equations with the assumptions noted above.

Buell, J. C.↗

Simulation of charged particle transport on the MPP

Computations of cosmic ray transport based upon finite difference methods are afflicted by instabilities, inaccuracies, and artifacts. To avoid these problems, a Monte Carlo formulation was developed which is closely related not only to the finite difference formulation, but also to the underlying physics of transport phenomena. Implementations of this approach are currently running on the Massively Parallel Processor at Goddard, whose enormous computing power overwhelms the poor statistical accuracy that usually limits the use of stochastic methods. In a Monte Carlo simulation of rectilinear transport, the coherent and diffusive effects that appeared are in good quantitative agreement with both finite difference and analytical calculations.

Earl, James A.↗

Statistical Mechanics and Dynamics of the Outer Solar System.I. The Jupiter/Saturn Zone

We report on numerical simulations designed to understand how the solar system evolved through a winnowing of planetesimals accreeted from the early solar nebula. This sorting process is driven by the energy and angular momentum and continues to the present day. We reconsider the existence and importance of stable niches in the Jupiter/Saturn Zone using greatly improved numerical techniques based on high-order optimized multi-step integration schemes coupled to roundoff error minimizing methods.

numerical simulations planetesimals planet formati↗

Software for Simulating Remote Sensing Systems

The Application Research Toolbox (ART) is a collection of computer programs that implement algorithms and mathematical models for simulating remote sensing systems. The ART is intended to be especially useful for performing design-tradeoff studies and statistical analyses to support the rational development of design requirements for multispectral imaging systems. Among other things, the ART affords a capability to synthesize coarser-spatial-resolution image-data sets from finer-spatial-resolution data sets and multispectral-image-data products from hyperspectral-image-data products. The ART also provides for synthesis of image-degradation effects, including point-spread functions, misregistration of spectral images, and noise. The ART can utilize real or synthetic data sets, along with sensor specifications, to create simulated data sets. In one example of a typical application, simulated data pertaining to an existing multispectral sensor system are used to verify the data collected by the system in operation. In the case of a proposed sensor system, the simulated data can be used to conduct trade studies and statistical analyses to ensure that the sensor system will satisfy the requirements of potential scientific, academic, and commercial user communities.

Zanoni, Vicki↗

Software for Simulating Remote Sensing Systems

The Application Research Toolbox (ART) is a collection of computer programs that implement algorithms and mathematical models for simulating remote sensing systems. The ART is intended to be especially useful for performing design-tradeoff studies and statistical analyses to support the rational development of design requirements for multispectral imaging systems. Among other things, the ART affords a capability to synthesize coarser-spatial-resolution image-data sets from finer-spatial-resolution data sets and multispectral-image-data products from hyperspectral-image-data products. The ART also provides for synthesis of image-degradation effects, including point-spread functions, misregistration of spectral images, and noise. The ART can utilize real or synthetic data sets, along with sensor specifications, to create simulated data sets. In one example of a typical application, simulated data pertaining to an existing multispectral sensor system are used to verify the data collected by the system in operation. In the case of a proposed sensor system, the simulated data can be used to conduct trade studies and statistical analyses to ensure that the sensor system will satisfy the requirements of potential scientific, academic, and commercial user communities.

Vicki Zanoni↗

Application of a Monte Carlo accuracy assessment tool to TDRS and GPS

In support of a NASA study on the application of radio interferometry to satellite orbit determination, MITRE developed a simulation tool for assessing interferometric tracking accuracy. Initially, the tool was applied to the problem of determining optimal interferometric station siting for orbit determination of the Tracking and Data Relay Satellite (TDRS). Subsequently, the Orbit Determination Accuracy Estimator (ODAE) was expanded to model the general batch maximum likelihood orbit determination algorithms of the Goddard Trajectory Determination System (GTDS) with measurement types including not only group and phase delay from radio interferometry, but also range, range rate, angular measurements, and satellite-to-satellite measurements. The user of ODAE specifies the statistical properties of error sources, including inherent observable imprecision, atmospheric delays, station location uncertainty, and measurement biases. Upon Monte Carlo simulation of the orbit determination process, ODAE calculates the statistical properties of the error in the satellite state vector and any other parameters for which a solution was obtained in the orbit determination. This paper presents results from ODAE application to two different problems: (1)determination of optimal geometry for interferometirc tracking of TDRS, and (2) expected orbit determination accuracy for Global Positioning System (GPS) tracking of low-earth orbit (LEO) satellites. Conclusions about optimal ground station locations for TDRS orbit determination by radio interferometry are presented, and the feasibility of GPS-based tracking for IRIDIUM, a LEO mobile satellite communications (MOBILSATCOM) system, is demonstrated.

Pavloff, Michael S.↗

Excercise Within LBNP as an Artificial Gravity Countermeasure

Previous exercise in space has lacked sufficient loads to maintain preflight cardiovascular and musculoskeletal mass and function. Lower body negative pressure (LBNP) produces a static force equivalent to one Earth body weight by each 52 mm Hg of LBNP during supine posture. LBNP also provides transmural blood pressures simulating upright exercise. Thus, this artificial-gravity concept may help maintain cardiovascular and musculoskeletal systems of crewmembers during prolonged exposure to microgravity. Currently available, bungee cord assisted, treadmill exercise is limited by harness discomfort, lower than normal loads, abnormal post-flight gait, and the absence of gravitational blood pressures within the vascular system. PURPOSE: This project evaluates a method to create artificial gravity using supine LBNP treadmill exercise to prevent loss of physiologic function in microgravity simulated by 30 days of bed rest. Identical twins were used as volunteers so that statistical power could be maximized. This countermeasure is being transitioned to space flight. CURRENT STATUS OF RESEARCH Methods: Six sets of identical twins (6 females and 14 males, 21-36 years) remained in 6 head-down tilt (HDT) bed rest for 30 days to simulate prolonged microgravity. Six subjects were randomly selected to exercise supine in an LBNP chamber for 40 minutes six days per week (EX group), while their twin brothers served as non-exercise controls (CON). Pressure within the exercise LBNP chamber was adjusted to increase load, hence increasing exercise intensity. During supine treadmill exercise, LBNP (52-63 mmHg) was applied to produce foot ward forces equivalent to those for upright running on Earth at 1.0-1.2 times body weight (BW) and subjects performed an interval exercise protocol (40-80% peak exercise capacity [VO2pk]). Five minutes of resting LBNP immediately followed each exercise session. Results: Orthostatic tolerance time decreased significantly after 30 days bed rest in the CON group, but was relatively maintained in the EX group. VO2pk was maintained in EX males, but not in CON males. Isokinetic knee strength (extension, peak torque) decreased significantly in CON males, but was preserved in EX males. The EX group had significantly higher spine muscle strength after bed rest than the CON group. The cross-sectional area of spinal muscle at L4/5 level decreased significantly in the CON group but not in the EX group. Urinary n-telopeptide excretion, an index of bone resorption, was increased during bed rest in CON, but not in EX subjects. This suggests protection by LBNP exercise against the increase in bone resorption typically seen in simulated and actual microgravity. Significant changes in bone mineral density (BMD) in the spine and ribs were observed in CON subjects, but not in EX subjects. Conclusions: Our treadmill exercise protocol within LBNP plus a short period of post-exercise LBNP maintains orthostatic responses, upright exercise capacity and other important physiologic parameters during bed rest. These results document the efficacy of our apparatus and exercise protocol for maintaining physiologic structure and function during long-duration microgravity as simulated by 30 days of HDT bed rest. FUTURE PLANS: More sets of female identical twins are needed to reach significance. The LBNP exercise chamber will be redesigned for flight.

Hargens, A. R.↗

Juneau Airport Doppler Lidar Deployment: Extraction of Accurate Turbulent Wind Statistics

A 2 micrometer pulsed Doppler lidar was deployed to the Juneau Airport in 1998 to measure turbulence and wind shear in and around the departure and arrival corridors. The primary objective of the measurement program was to demonstrate and evaluate the capability of a pulsed coherent lidar to remotely and unambiguously measure wind turbulence. Lidar measurements were coordinated with flights of an instrumented research aircraft operated by representatives of the University of North Dakota (UND) under the direction of the National Center for Atmospheric Research (NCAR). The data collected is expected to aid both turbulence characterization as well as airborne turbulence detection algorithm development activities within NASA and the FAA. This paper presents a summary of the deployment and results of analysis and simulation which address important issues regarding the measurement requirements for accurate turbulent wind statistics extraction.

Hannon, Stephen M.↗

Piloted simulation evaluation of pitch control designs for highly augmented STOVL aircraft

Analyses of design variations on a pitch axis stabilization and command augmentation system (SCAS) for a STOVL fighter aircraft are performed in a moving base simulation experiment. The primary goal of this study is to determine if turbulence-induced control activity could be reduced by modifying SCAS parameters while keeping the response-to-command characteristics of the baseline system that provide Level 1 flying qualities. Pilot ratings and control utilization statistics for the baseline system are in agreement with similar data gathered in a prior simulation test involving the same aircraft and control system.

Engelland, S. A.↗

Tensoral: A system for post-processing turbulence simulation data

Many computer simulations in engineering and science -- and especially in computational fluid dynamics (CFD) -- produce huge quantities of numerical data. These data are often so large as to make even relatively simple post-processing of this data unwieldy. The data, once computed and quality-assured, is most likely analyzed by only a few people. As a result, much useful numerical data is under-utilized. Since future state-of-the-art simulations will produce even larger datasets, will use more complex flow geometries, and will be performed on more complex supercomputers, data management issues will become increasingly cumbersome. My goal is to provide software which will automate the present and future task of managing and post-processing large turbulence datasets. My research has focused on the development of these software tools -- specifically, through the development of a very high-level language called 'Tensoral'. The ultimate goal of Tensoral is to convert high-level mathematical expressions (tensor algebra, calculus, and statistics) into efficient low-level programs which numerically calculate these expressions given simulation datasets. This approach to the database and post-processing problem has several advantages. Using Tensoral the numerical and data management details of a simulation are shielded from the concerns of the end user. This shielding is carried out without sacrificing post-processor efficiency and robustness. Another advantage of Tensoral is that its very high-level nature lends itself to portability across a wide variety of computing (and supercomputing) platforms. This is especially important considering the rapidity of changes in supercomputing hardware.

Dresselhaus, Eliot↗

Flying the TRMM Satellite in a GCM

By incorporating the Tropical Rainfall Measurement Mission (TRMM) satellite orbital information into the Colorado State University General Circulation Model (CSU GCM), we are able to 'fly' a satellite in the GCM, and sample the simulated atmosphere in the same way as the TRMM sensors sample the real atmosphere. The TRMM-sampled statistics for precipitation and radiative fluxes at annual, intraseasonal, monthly-mean and seasonal-mean diurnal time scales are evaluated by comparing the satellite-sampled against fully-sampled simulated atmospheres. The sampling rates of the TRMM sensors are significantly affected by the sensors' swath widths. The TRMM Microwave Imager (TMI) and the Visible Infrared Scanner (VIRS) sample each 2.25 x 2.25 degree grid box in the Tropics and subtropics about once per day, but at a different local time every day, while the Precipitation Radar (PR) and the Clouds and the Earth's Radiant Energy System (CERES) sensor visit each grid box about once every three days and twice per day, respectively. Besides inadequate samplings resulting from sensors' swath widths, there is a large, systematic diurnal undersampling associated with TRMM's orbital geometry for grid boxes away from the Equator. When only one month of TRMM data are used, this diurnal undersampling can lead to more daytime samples relative to nighttime samples in one hemisphere, and more nighttime samples relative to daytime samples in the other hemisphere. The resulting sampling biases (3-6 W m(exp-2)) are very pronounced in outgoing longwave radiation (OLR) over the subtropical land masses. The sampling errors in OLR monthly- and seasonal-means are less than 8 W m(exp-2) (5%) for each 2.25 x 2.25 degree grid box. The OLR monthly- and seasonal-means are not sensitive to diurnal undersamplings associated with the TRMM orbits and sensors' swath widths. However, this is not the case for total precipitation. Diurnal undersampling could produce errors as large as 20% in the Tropics and 40% in the subtropics, for the zonally averaged monthly-mean rain rates. The TRMM orbits sample each 2.25 x 2.25 degree grid box in the Tropics and subtropics 1-6 times for each hour of the day within a single season. The seasonal-mean diurnal cycles of precipitation and OLR are not well sampled for any one grid box. By either accumulating the satellite data for a long enough period, or averaging the data over a large area with a relatively uniform diurnal signal, the diurnal cycles of precipitation and OLR can be satisfactorily sampled. The effects TRMM sampling errors on the inferred tropical-mean hydrologic cycle and radiative fluxes are also evaluated. There are strong spurious oscillations associated with TRMM's orbital geometry, with periods of 23 days and 3-4 months, in tropical-mean daily and monthly precipitation. While the relative fluctuations of the sampled-OLR are negligible, the relative fluctuations of the sampled precipitation have magnitudes similar to those of the observed climate variability. Caution must therefore be used when applying TRMM observations of tropical-mean precipitation to interpret climate variations at intraseasonal and interannual scales.

Lin, Xin↗

The Implementation of Maximum Likelihood Estimation in Space Launch System Vehicle Design

The Space Launch System uses a Maximum Likelihood Estimation process in conjunction with Design of Experiments to develop statistically representative vehicles for the Block 1 configuration. These vehicles are then used to estimate maximum load conditions for simulating stressing cases in other simulations. This paper discusses the modeling process and how SLS captures manufacturing uncertainty in the launch vehicle design. It also provides an overview of the differences between Block 1 statistical representations. This paper also discusses proper grid choice as well as which uncertainties drive the vehicle design.

Stein, W. B.↗

Large-Eddy Simulation of Turbulent Wall-Pressure Fluctuations

Large-eddy simulations of a turbulent boundary layer with Reynolds number based on displacement thickness equal to 3500 were performed with two grid resolutions. The computations were continued for sufficient time to obtain frequency spectra with resolved frequencies that correspond to the most important structural frequencies on an aircraft fuselage. The turbulent stresses were adequately resolved with both resolutions. Detailed quantitative analysis of a variety of statistical quantities associated with the wall-pressure fluctuations revealed similar behavior for both simulations. The primary differences were associated with the lack of resolution of the high-frequency data in the coarse-grid calculation and the increased jitter (due to the lack of multiple realizations for averaging purposes) in the fine-grid calculation. A new curve fit was introduced to represent the spanwise coherence of the cross-spectral density.

Singer, Bart A.↗

Investigation of radiative interaction in laminar flows using Monte Carlo simulation

The Monte Carlo method (MCM) is employed to study the radiative interactions in fully developed laminar flow between two parallel plates. Taking advantage of the characteristics of easy mathematical treatment of the MCM, a general numerical procedure is developed for nongray radiative interaction. The nongray model is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. To validate the Monte Carlo simulation for nongray radiation problems, the results of radiative dissipation from the MCM are compared with two available solutions for a given temperature profile between two plates. After this validation, the MCM is employed to solve the present physical problem and results for the bulk temperature are compared with available solutions. In general, good agreement is noted and reasons for some discrepancies in certain ranges of parameters are explained.

Liu, Jiwen↗

Characterization Study of TestBed Infrastructure Performance in a Distributed Simulation Environment: Baseline Analysis

Characterization of the performance of Air Traffic Management Exploration (ATM-X) TestBed integration environment has been investigated and documented for one system configuration for progressively increasing traffic. Several statistical parameters were used to assess the performance of the TestBed distributed system such as mean, standard deviation, skewness, and kurtosis of latency, and update rate for aircraft state messages that are transmitted through the simulated system under investigation. It is necessary to assess the performance characteristics of distributed systems in terms of the indicated statistical parameters mentioned above. It is critical to verify the system performance with respect to a researcher’s required system performance. Computer host specifications are documented in terms of Central Processing Unit (CPU) clock speed and core count. Transmission Control Protocol/ Internet Protocol (TCP/IP) message protocol was used for data transmission. The system network topology also contributes to the latency and update rate variations from the one imposed by the data source. The motivation for selecting the TestBed infrastructure as the focus of this study can be attributed to the number of services and capabilities it provides that help simplify the process of preparing and conducting a simulation. These capabilities include an easy to use GUI for simulation configuration, access to TestBed library by the end-user of other simulation software components, a modular adapter paradigm that allows simple connectivity of external software to TestBed, connectivity with other simulation laboratories, and a Software Development Kit (SDK) for quicker development. Two types of traffic generators, Air Traffic Generator (ATG) and Multi Aircraft Control System (MACS) were used to generate messages that were injected into the TestBed distributed environment. Eight different air traffic scenarios with progressively increasing loads were generated for each air traffic simulator. The corresponding air traffic loads between the two simulators had an identical number of aircraft per scenario, but different flight plans. It was observed that the performance of MACS degraded for air traffic scenarios containing more than 200 aircraft (37.5 KB/s nominal throughput). However, ATG performed adequately under all tested air traffic loads up to 1200 aircraft (225. KB/s nominal throughput). The tests show that MACS exhibits better latency performance with smaller aircraft loads when compared to ATG. The tests also show that the TestBed infrastructure successfully transmits 1200 aircraft without significant degradation of its performance. From the latency trends for both MACS and ATG, it is clear that as aircraft load increases, the latency in the system increases as well as its standard deviation. Likewise, the trends for the update data rate for both MACS and ATG show that as the aircraft load increases, so does the standard deviation and mean of the update rates which can be attributed to the performance of MACS and ATG applications. The analysis of the results of this study have proven that the overall system performance is dependent on the individual performance of each system component that is connected to TestBed, which subsequently propagates into the system. All TestBed characterization tests were conducted in SimLabs at NASA Ames Research Center in November 2019. This study addresses the need for a baseline TestBed characterization, and the results will serve as a reference for more complex simulation systems.

Air Traffic Management simulations↗

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗

Numerical Study of Pressure Fluctuations due to High-Speed Turbulent Boundary Layers

Direct numerical simulations (DNS) are used to examine the pressure fluctuations generated by fully developed turbulence in supersonic turbulent boundary layers with an emphasis on both pressure fluctuations at the wall and the acoustic fluctuations radiated into the freestream. The wall and freestream pressure fields are first analyzed for a zero pressure gradient boundary layer with Mach 2.5 and Reynolds number based on momentum thickness of approximately 2835. The single and multi-point statistics reported include the wall pressure fluctuation intensities, frequency spectra, space-time correlations, and convection velocities. Single and multi-point statistics of surface pressure fluctuations show good agreement with measured data and previously published simulations of turbulent boundary layers under similar flow conditions. Spectral analysis shows that the acoustic fluctuations outside the boundary layer region have much lower energy content within the high-frequency region. The space-time correlations reflect the convective nature of the pressure field both at the wall and in the freestream, which is characterized by the downstream propagation of pressure-carrying eddies. Relative to those at the wall, the pressure-carrying eddies associated with the freestream signal are larger and convect at a significantly lower speed. The preliminary DNS results of a Mach 6 boundary layer show that the pressure rms in the freestream region is significantly higher than that of the lower Mach number case.

Duan, Lian↗