Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Numerical evaluation of the radiation from unbaffled, finite plates using the FFT

An iteration technique is described which numerically evaluates the acoustic pressure and velocity on and near unbaffled, finite, thin plates vibrating in air. The technique is based on Rayleigh's integral formula and its inverse. These formulas are written in their angular spectrum form so that the fast Fourier transform (FFT) algorithm may be used to evaluate them. As an example of the technique the pressure on the surface of a vibrating, unbaffled disk is computed and shown to be in excellent agreement with the exact solution using oblate spheroidal functions. Furthermore, the computed velocity field outside the disk shows the well-known singularity at the rim of the disk. The radiated fields from unbaffled flat sources of any geometry with prescribed surface velocity may be evaluated using this technique. The use of the FFT to perform the integrations in Rayleigh's formulas provides a great savings in computation time compared with standard integration algorithms, especially when an array processor can be used to implement the FFT.

Williams, E. G.↗
NASA NTRS1990

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems, task 1: Ducted propfan analysis

The time-dependent three-dimensional Euler equations of gas dynamics were solved numerically to study the steady compressible transonic flow about ducted propfan propulsion systems. Aerodynamic calculations were based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. An implicit residual smoothing operator was used to aid convergence. Two calculation grids were employed in this study. The first grid utilized an H-type mesh network with a branch cut opening to represent the axisymmetric cowl. The second grid utilized a multiple-block mesh system with a C-type grid about the cowl. The individual blocks were numerically coupled in the Euler solver. Grid systems were generated by a combined algebraic/elliptic algortihm developed specifically for ducted propfans. Numerical calculations were initially performed for unducted propfans to verify the accuracy of the three-dimensional Euler formulation. The Euler analyses were then applied for the calculation of ducted propfan flows, and predicted results were compared with experimental data for two cases. The three-dimensional Euler analyses displayed exceptional accuracy, although certain parameters were observed to be very sensitive to geometric deflections. Both solution schemes were found to be very robust and demonstrated nearly equal efficiency and accuracy, although it was observed that the multi-block C-grid formulation provided somewhat better resolution of the cowl leading edge region.

Hall, Edward J.↗
NASA NTRS2002

Computation of High Speed Jet Noise

The objective of this work was to use direct numerical simulation (DNS) techniques to study the physics of noise generation by a high-speed turbulent jet. A Mach 0.9, Reynolds number 3,600 jet was selected because of available experimental data. New numerical methods for generating disturbances at the nozzle and computing far-field sound were developed and reported in the course of this work. Over 25 million mesh points were used in the simulations which ran for over 50,000 timesteps and required over 50,000 processor hours on state-of-the-art parallel computer systems to complete. Figures show a visualization of the jet and sound field, a comparison of the mean flow development with the experiment, a directivity comparison with the experiment, and time spectrum comparison with the experiment. Agreement is seen to be excellent. These are fully document in the attached references. Full details of the work, detailed achievements and conclusions are discussed in appendices, which are copies of publications that resulted from this work. We have studied noise mechanisms in supersonic jets, the refraction of sound by turbulence in subsonic jets, and noise sources in conjunction with a DNS of a Mach 0.9 jet.

Freund, Jonathan B.↗
NASA NTRS2005

Controlling Reflections from Mesh Refinement Interfaces in Numerical Relativity

A leading approach to improving the accuracy on numerical relativity simulations of black hole systems is through fixed or adaptive mesh refinement techniques. We describe a generic numerical error which manifests as slowly converging, artificial reflections from refinement boundaries in a broad class of mesh-refinement implementations, potentially limiting the effectiveness of mesh- refinement techniques for some numerical relativity applications. We elucidate this numerical effect by presenting a model problem which exhibits the phenomenon, but which is simple enough that its numerical error can be understood analytically. Our analysis shows that the effect is caused by variations in finite differencing error generated across low and high resolution regions, and that its slow convergence is caused by the presence of dramatic speed differences among propagation modes typical of 3+1 relativity. Lastly, we resolve the problem, presenting a class of finite-differencing stencil modifications which eliminate this pathology in both our model problem and in numerical relativity examples.

Baker, John G.↗
NASA NTRS1991

Numerical and literal aeroelastic-vehicle-model reduction for feedback control synthesis

The simplification of a high-order, literal model for large flexible aircraft is discussed. Areas of model fidelity that are critical if the model is to be used for control law synthesis are presented. Several simplification techniques, some new and some widely available, that can deliver the necessary model fidelity are presented and applied to a model from the literature. The techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory, is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle vibrational characteristics and rigid-body and aeroelastic stability derivatives, thus providing insight into the underlying causes for critical dynamic characteristics.

Newman, Brett↗
DOE OSTI2023

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗
NASA NTRS2002

Land Surface Data Assimilation and the Gulf Coast Sea-Breeze

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The focus of this paper is to examine how the satellite assimilation technique impacts simulations of near-surface meteorology on the 4- to 48-hour time scale. The technique is implemented within the PSU/NCAR MM5 V34 and applied on a 36-km CONUS domain and a 12-km nest centered over Oklahoma. Two Friday periods during 2000 were selected for study, one during the warm season (May-June) and the other in the cold season (Nov-Dec). Bulk verification statistics (BIAS and RMSE) of surface air temperature and dewpoint show that the assimilation technique can improve numerical simulations on both regional and continental scales. Comparison of the simulated surface energy and water fluxes with observations from Energy Balance Bowen Ratio measurements taken at CART/ARM sites in the Southern Great Plains showed the assimilation technique produces a realistic argument of the Bowen ratio in both seasons. Additional runs with the Oregon State University land surface scheme as implemented within MM5 will be made and comparisons with assimilation results will be presented at the meeting.

Lapenta, William M.↗
NASA NTRS1972

Analysis and design of three dimensional supersonic nozzles. Volume 1: Nozzle-exhaust flow field analysis by a reference plane characteristics technique

A second order numerical method employing reference plane characteristics has been developed for the calculation of geometrically complex three dimensional nozzle-exhaust flow fields, heretofore uncalculable by existing methods. The nozzles may have irregular cross sections with swept throats and may be stacked in modules using the vehicle undersurface for additional expansion. The nozzles may have highly nonuniform entrance conditions, the medium considered being an equilibrium hydrogen-air mixture. The program calculates and carries along the underexpansion shock and contact as discrete discontinuity surfaces, for a nonuniform vehicle external flow.

Dash, S.↗
NASA NTRS1993

Techniques and resources for storm-scale numerical weather prediction

The topics discussed include the following: multiscale application of the 5th-generation PSU/NCAR mesoscale model, the coupling of nonhydrostatic atmospheric and hydrostatic ocean models for air-sea interaction studies; a numerical simulation of cloud formation over complex topography; adaptive grid simulations of convection; an unstructured grid, nonhydrostatic meso/cloud scale model; efficient mesoscale modeling for multiple scales using variable resolution; initialization of cloud-scale models with Doppler radar data; and making effective use of future computing architectures, networks, and visualization software.

Droegemeier, Kelvin↗
NASA NTRS1995

Nonperiodic variations in astrophysical systems: Investigating frequency evolution

We present a method related to the wavelet transform, the Gabor transform, for investigating astronomical time series containing nonconsistent frequencies. Instances in which such data sets may arise include variable star light curves, numerical studies of the gravitational three-body problem, X-ray binaries, and signals from more exotic objects such as planets around pulsars and mass infall from accretion structures onto compact objects. As an illustration of its power, we apply the technique to a numerical data set of a gravitational three-body interaction and to photometry of the rapidly oscillating peculiar A star HD 60435. In the three-body example, the method provides an insightful shorthand that allows for the determination of episodes where the system behaves as two nearly Keplerian orbits. For HD 60435, the power in the main frequency exhibits unusual evolution over the duration of the observation.

Boyd, Patricia T.↗
NASA NTRS2001

Assimilation of GOES Land Surface Data: Benefits to Numerical Weather Prediction

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The technique has been employed on a semi-operational basis at the GHCC within the PSU/CAR MM5 since 1 November 1998. A one-way nested grid configuration was employed with a 75 kin CONUS domain and a 25 km grid over the southeastern United States. Initial conditions were obtained from the 12 UTC Early Eta Data Assimilation System analyses and lateral boundary conditions from the Early Eta forecast available at 3 hour intervals. The satellite-derived land surface temperature tendencies and insolation were assimilated between two and five hours (1400 and 1650 UTC) of the forecast. We performed the assimilation on the Southeastern domain only. In addition, a control run without assimilation was performed to provide insight into the performance of the assimilation technique.

Lapenta, William M.↗
NASA NTRS1977

A balanced diagnostic system compatible with a barotropic prognostic model

A system of diagnostic equations for the velocity field, or wind laws, for a barotropic primitive-equation model of large-scale atmospheric flow is derived. Attention is given to the classical balance equation and its ellipticity condition. Numerical solutions of the diagnostic system are presented, including examples of cases of the mixed elliptic-hyperbolic type and cases with non-zero divergence. Procedures for implementing such a system are outlined, along with a review of factors in using the technique for operational numerical weather prediction.

Ghil, M.↗
NASA NTRSDate not supplied

Dense Feature Tracking of Atmospheric Winds with Deep Optical Flow

Atmospheric winds are a key physical phenomenon impacting natural hazards, energy transport, ocean currents, large-scale circulation, and ecosystem fluxes. Observing winds is a complex process and presents a large gap in NASA’s Earth Observation System. Atmospheric motion vectors (AMVs) aim to fill this gap by making numerical estimates of cloud movement between sequences of multi-spectral satellite images, tracking clouds and water vapor. Recent imaging hardware and software advancements have enabled the use of numerical optical flow techniques to produce accurate and dense vector fields outperforming traditional methods. This work presents WindFlow as the first machine learning based system for feature tracking atmospheric motion using optical flow. Due to the lack of large-scale satellite-based observations, we leverage high-resolution numerical simulations from NASA's GEOS-5 Nature Run to perform supervised learning and transfer to satellite images. We demonstrate that our approach using deep learning based optical flow scales to ultra-high-resolution images of size 2881x5760 with less than 1 m/s bias and 2.5 m/s average error. Four network and learning architectures are compared and it is found that recurrent all-pairs field transforms (RAFT) produces the lowest errors on all metrics for wind speed and direction. Results on held out numerical outputs shows RAFT's good performance in each of the spatial, temporal, and physical dimensions. A comparison between WindFlow and an operational AMV product against rawinsonde observations show that RAFT transfers across simulations and thermal infrared satellite observations. This work shows that machine learning based optical flow is an efficient approach to generating robust feature tracking for AMVs consistently over large regions.

Atmospheric winds↗
NASA NTRS1966

An orthonormalization procedure for multivariable function approximation

Where a function of several variables is given numerically in tabular form, an orthonormalization technique allows an approximation of the numerical data to be determined in a convenient functional form. In this technique, the speed and accuracy of coefficient computation are much improved.

Ingram, H. L.↗