Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computation time”

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

Check out of rebuilt optical cross-beam units

Preliminary evaluations of limited cross beam data taken with the rebuilt optical cross beam units at the Gun Barrel Hill Test Site of ESSA are presented. The evaluations were made using the analogue Princeton Time Correlator. A study of these evaluations indicates that the rebuilt units are performing satisfactorily. However, there are indications that, in order to obtain useful information from the cross-correlations, it is essential to make judicious choices of computation time and frequency band pass filtering for the analysis of the data.

Huff, L. L.↗

Flap-lag-torsional dynamics of extensional and inextensional rotor blades in hover and in forward flight

The reduction of the O(cu epsilon) integro differential equations to ordinary differential equations using a set of orthogonal functions is described. Attention was focused on the hover flight condition. The set of Galerkin integrals that appear in the reduced equations was evaluated by making use of nonrotating beam modes. Although a large amount of computer time was needed to accomplish this task, the Galerkin integrals so evaluated were stored on tape on a permanent basis. Several of the coefficients were also obtained in closed form in order to check the accuracy of the numerical computations. The equilibrium solution to the set of 3n equations obtained was determined as the solution to a minimization problem.

Dasilva, C.↗

The efficient numerical calculation of condensational cloud drop growth

A modified Runge-Kutta integration technique is applied to condensational cloud droplet growth calculations, using the device of limiting dr, the growth of a droplet of radius r over the time step dt to some function of the absolute value of r(eq)-r, where r(eq) denotes the equilibrium value. The use of variable time steps to obtain an efficient choice of dt is also discussed. Results using these techniques are compared with those obtained with the standard fourth-order Runge-Kutta method, showing that the modified technique reduces the requisite computer time by a factor of about one hundred.

Robinson, N. F.↗

Three parallel computation methods for structural vibration analysis

The Lanczos (1950), multisectioning, and subspace iteration sequential methods for vibration analysis presently used as bases for three parallel algorithms are noted, in the aftermath of three example problems, to maintain reasonable accuracy in the computation of vibration frequencies. Significant computation time reductions are obtained as the number of processors increases. An analysis is made of the performance of each method, in order to characterize relative strengths and weaknesses as well as to identify those parameters that most strongly affect computation efficiency.

Storaasli, Olaf↗

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗

RApid Blade and Blade-Vortex InTeraction (RABBIT) Reimagined: Algorithm Improvement and User-Friendly Interface Implementation

Rotorcraft noise source identification and reduction is crucial to the emergence of the Urban Air Mobility (UAM) market. One key rotor noise source is Blade Vortex Interaction (BVI), caused by the rotor blades interacting with the rotor wake. The low-fidelity RApid Blade and Blade-Vortex InTeraction (RABBIT) tool was created to predict the location and characteristics of BVI noise. RABBIT enables engineers to quickly design and simulate UAM configurations and understand BVI occurrences for a desired flight condition. Furthermore, RABBIT can also identify the time and location of blade overlap for coaxial rotors. RABBIT is not a general acoustic prediction tool, but instead utilizes vortex and wake parameters to visualize and predict BVI aerodynamics. RABBIT algorithm has been further developed including reduction of computational time, accounting for BVI source time arrival, improved vortex strength modeling, and BVI impulse factor. A user-friendly app is presented with the intent of bringing RABBIT to the masses to allow engineers from various educational backgrounds to design and begin to understand BVI. Furthermore, RABBIT capabilities are highlighted by a parameter study highlighting key parameters that impact BVIalong with highlighting RABBIT’s ability to identify BVI occurrences from other rotors.

RApid↗

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

Comparing quantile regression forest and mixture density long short-term memory models for probabilistic post-processing of satellite precipitation-driven streamflow simulations

Abstract. Deep learning (DL) and machine learning (ML) are widely used in hydrological modelling, which plays a critical role in improving the accuracy of hydrological predictions. However, the trade-off between model performance and computational cost has always been a challenge for hydrologists when selecting a suitable model, particularly for probabilistic post-processing with large ensemble members. This study aims to systematically compare the quantile regression forest (QRF) model and countable mixtures of asymmetric Laplacians long short-term memory (CMAL-LSTM) model as hydrological probabilistic post-processors. Specifically, we evaluate their ability in dealing with biased streamflow simulations driven by three satellite precipitation products across 522 nested sub-basins of the Yalong River basin in China. Model performance is comprehensively assessed using a series of scoring metrics from both probabilistic and deterministic perspectives. Our results show that the QRF model and the CMAL-LSTM model are comparable in terms of probabilistic prediction, and their performances are closely related to the flow accumulation area (FAA) of the sub-basin. The QRF model outperforms the CMAL-LSTM model in most sub-basins with smaller FAA, while the CMAL-LSTM model has an undebatable advantage in sub-basins with FAA larger than 60 000 km2 in the Yalong River basin. In terms of deterministic predictions, the CMAL-LSTM model is preferred, especially when the raw streamflow is poorly simulated and used as input. However, setting aside the differences in model performance, the QRF model with 100-member quantiles demonstrates a noteworthy advantage by exhibiting a 50 % reduction in computation time compared to the CMAL-LSTM model with the same ensemble members in all experiments. As a result, this study provides insights into model selection in hydrological post-processing and the trade-offs between model performance and computational efficiency. The findings highlight the importance of considering the specific application scenario, such as the catchment size and the required accuracy level, when selecting a suitable model for hydrological post-processing.

Geology↗

Flow in nonrotating passages of radial inflow turbines

The analysis of irrotational incompressible flow field in the stator unit of a radial inflow turbine is presented. The solution in the combined scroll-nozzle assembly is complicated by the domain geometry and by its multiconnectivity. This model is necessary, however, in order to provide a better understanding of the mutual interaction effects of these two components on the flow field. The finite element method is used in the solution which is limited to the two dimensional case. A substructuring technique is adopted in the computational procedure and results in considerable savings in both computer time and core storage requirements. The results are presented for the flow velocity magnitude and direction in the scroll and through the various nozzles, for two nozzle blade geometries. In addition, the mass flow rates in the different nozzles are computed and their deviations from the mean value determined.

Baskharone, E.↗

Maximum likelihood identification of aircraft parameters with unsteady aerodynamic modelling

A simplified aerodynamic force model based on the physical principle of Prandtl's lifting line theory and trailing vortex concept has been developed to account for unsteady aerodynamic effects in aircraft dynamics. Longitudinal equations of motion have been modified to include these effects. The presence of convolution integrals in the modified equations of motion led to a frequency domain analysis utilizing Fourier transforms. This reduces the integro-differential equations to relatively simple algebraic equations, thereby reducing computation time significantly. A parameter extraction program based on the maximum likelihood estimation technique is developed in the frequency domain. The extraction algorithm contains a new scheme for obtaining sensitivity functions by using numerical differentiation. The paper concludes with examples using computer generated and real flight data

Keskar, D. A.↗

On-Line Robust Modal Stability Prediction using Wavelet Processing

Wavelet analysis for filtering and system identification has been used to improve the estimation of aeroservoelastic stability margins. The conservatism of the robust stability margins is reduced with parametric and nonparametric time- frequency analysis of flight data in the model validation process. Nonparametric wavelet processing of data is used to reduce the effects of external disturbances and unmodeled dynamics. Parametric estimates of modal stability are also extracted using the wavelet transform. Computation of robust stability margins for stability boundary prediction depends on uncertainty descriptions derived from the data for model validation. The F-18 High Alpha Research Vehicle aeroservoelastic flight test data demonstrates improved robust stability prediction by extension of the stability boundary beyond the flight regime. Guidelines and computation times are presented to show the efficiency and practical aspects of these procedures for on-line implementation. Feasibility of the method is shown for processing flight data from time- varying nonstationary test points.

Brenner, Martin J.↗

Numerical Simulation of Inlet Bleed with Circular Holes on Plate Under Shock-Wave/Boundary-Layer Interactions

A numerical study was performed to investigate the shock-wave/boundary-layer interactions on a flat plate with bleed through one or more circular holes that vent into a plenum. The bleed-hole patterns considered for the study include in-line multiple holes and staggered multiple-row holes that are configured to simulate the patterns used in inlet bleed systems of high performance aircraft. The focus of the study was to examine how the bleed through multiple holes affect bleed rate and the pressure and Mach number distributions. Since the bleed performance was found sensitive to the change in bleed conditions, a computational procedure was developed to give a good turnaround computational time for parametric studies involving changes in bleed hole geometry and the structure of shock-wave/boundary-layer flowfield. The procedure includes the grid-generation methodology and the flow simulation with solutions from the Navier-Stokes equations. The computational techniques permit analysis of complex bleed systems and make possible the investigation of a broader range of design variables associated with inlet bleed operation.

Chyu, Wei J.↗

An experimental and computational investigation of the steady and unsteady transonic flow fields about an airfoil in a solid-wall test channel

An experimental and computational investigation of the steady and unsteady transonic flow field about a thick airfoil is described. An operational computer code for solving the two-dimensional, compressible Navier-Stokes equations for flow over airfoils was modified to include solid-wall, slipflow boundary conditions to properly assess the code and help guide the development of improved turbulence models. Steady and unsteady flow fields about an 18% thick circular arc airfoil at Mach numbers of 0.720, 0.754, and 0.783 and a chord Reynolds number of 11 x 10 to the 6th are predicted and compared with experiment. For the first time, computed results for unsteady turbulent flows with separation caused by a shock wave were obtained which qualitatively reproduce the time-dependent aspects of experiments. Features such as the intensity and reduced frequency of airfoil surface-pressure fluctuations, oscillatory regions of trailing-edge and shock-induced separation, and the Mach number range for unsteady flows were all qualitatively reproduced.

Levy, L. L., Jr.↗

Toward prediction of turbulent atmospheric flows over propagating oceanic waves via machine-learning augmented large-eddy simulation

Wind-wave interactions have important effects on the energy harvesting of offshore wind farms. High-fidelity large-eddy simulation (LES) is a powerful approach for investigating wind-wave interactions in turbulent oceanic environments. Due to the large scale of the flow domain and the high grid resolution required to resolve multi-scale flow motions, however, brute-force LES of wind-wave interactions is computationally very expensive. We propose augmenting brute-force LES via machine-learning data-driven modeling (ML-LES) to dramatically reduce the computational time required to obtain converged turbulence statistics when the brute-force approach is employed. Namely, we employ a convolutional neural network (CNN) autoencoder trained and validated with LES data sets to develop a highly efficient ML-LES approach for computing turbulence statistics from just a few snapshots of instantaneous LES flow fields. Further, our results demonstrate the accuracy and efficiency of ML-LES in predicting the mean velocity, velocity fluctuations, and turbulence kinetic energy in highly stretched computational grid systems required to carry out simulations in real-life oceanic environments.

42 ENGINEERING↗

In-flight thrust determination on a real-time basis

A real time computer program was implemented on a F-15 jet fighter to monitor in-flight engine performance of a Digital Electronic Engine Controlled (DEES) F-100 engine. The application of two gas generator methods to calculate in-flight thrust real time is described. A comparison was made between the actual results and those predicted by an engine model simulation. The percent difference between the two methods was compared to the predicted uncertainty based on instrumentation and model uncertainty and agreed closely with the results found during altitude facility testing. Data was obtained from acceleration runs of various altitudes at maximum power settings with and without afterburner. Real time in-flight thrust measurement was a major advancement to flight test productivity and was accomplished with no loss in accuracy over previous post flight methods.

Ray, R. J.↗

A comparison between numerically modelled and experimentally measured loss mechanisms in wave rotors

A numerical model has been developed which is capable of predicting the performance of a wave rotor (pressure exchanger) of specified geometry over a wide range of operating conditions. The model can account for the major loss mechanisms of leakage from the tube ends, fluid viscosity, heat transfer to the tube wails, finite tube opening time, shock waves, and non-uniform port flows. It is a one dimensional flow model which follows a single tube as it rotates past the various stationary ports. Since the model is relatively simple (i.e., one dimensional) it uses little computer time. This makes it suitable for design as well as analytical purposes. This paper will present a brief description of the model then discuss a comparison between the model predictions and several wave rotor experiments.

Paxson, Daniel E.↗

A Comparison Between Numerically Modelled and Experimentally Measured Loss Mechanisms in Wave Rotors

A numerical model has been developed which is capable of predicting the performance of a wave rotor (pressure exchanger) of specified geometry over a wide range of operating conditions. The model can account for the major loss mechanisms of leakage from the tube ends, fluid viscosity, heat transfer to the tube walls, finite tube opening time, shock waves, and non-uniform port flows. It is a one dimensional flow model which follows a single tube as it rotates past the various stationary ports. Since the model is relatively simple (i.e. one dimensional) it uses little computer time. This makes it suitable for design as well as analytical purposes. This paper will present a brief description of the model then discuss a comparison between the model predictions and several wave rotor experiments.

Paxson, Daniel E.↗