Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximations”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Navier-Stokes Flowfield Simulation of Boeing 747-200 as Platform for SOFIA

Steady and unsteady viscous, three-dimensional flowfields are calculated using a thin layer approximation of Navier-Stokes equations in conjunction with Chimera overset grids. The finite-difference numerical scheme uses structured grids and a pentadiagonal flow solver called "OVERFLOW". The configuration of Boeing 747-200 has been chosen as one of configurations to be used as a platform for the SOFIA (Stratospheric Observatory For Infrared Astronomy). Initially, the steady flowfield of the full aircraft is calculated for the clean configuration (without a cavity to house telescope). This solution is then used to start the unsteady flowfield of a configuration containing cavity housing the observation telescope and its peripheral units. Analysis of unsteady flowfield in the cavity and its influence on the tail empennage, as well as the noise due to turbulence and optical quality of the flow are the main focus of this study. For the configuration considered here, the telescope housing cavity is located slightly downstream of the portwing. The entire flow-field is carefully constructed using 45 overset grids and consists of nearly 4 million grid points. All the computations axe done at one freestream flow condition of M(sub infinity) = 0.85, alpha = 2.5deg, and a Reynolds of Re = 1.85x10deg

Srinivasan, G.R.↗
NASA NTRS2006

Sulfur during the Transition from Anoxic to Oxic Atmospheres

The invention of oxygenic photosynthesis was likely accompanied by the introduction of large amounts of O2 and complementary reduced gases (chiefly CH4) into the atmosphere. To first approximation the venting of O2 and CH4 are stochiometrically linked. We therefore present a suite of numerical photochemical models that address the anoxic-oxic transition in an atmosphere driven by large linked inputs of biogenic 02 and CH4. We find in general that, in steady state, there are two solutions, one oxic and the other anoxic. The anoxic solution appears to be linearly stable. If volcanic SO2 fluxes are large, S disproportionates into oxidized (H2S04) and reduced (S8) exit channels. As elemental sulfur is insoluble it provides a means of preserving photochemical mass-independent fractionation (MIF). On the other hand, if the source of volcanic SO2 is smaller than today, all S can leave the atmosphere as S8. Under these conditions there would be no MIF signal. The oxic solution appears to be linearly unstable. In the oxic solutions S is invariably oxidized to sulfate, and the MIF signal would be absent. The transitional atmosphere is relatively unstable and is also the most photochemically active. Consequently it is the transitional atmosphere, not the oxic or anoxic atmospheres, that has the lowest CH4 levels and weakest greenhouse warming. As a practical matter we expect the transitional atmospheres to vary strongly in response to diurnal and seasonal biological forcing.

Zahnle, Kevin↗
NASA NTRS2017

Online Visualization and Value Added Services of MERRA-2 Data at GES DISC

NASA climate reanalysis datasets from MERRA-2, distributed at the Goddard Earth Sciences Data and Information Services Center (GES DISC), have been used in broad research areas, such as climate variations, extreme weather, agriculture, renewable energy, and air quality, etc. The datasets contain numerous variables for atmosphere, land, and ocean, grouped into 95 products. The total archived volume is approximately 337 TB ( approximately 562K files) at the end of October 2017. Due to the large number of products and files, and large data volumes, it may be a challenge for a user to find and download the data of interest. The support team at GES DISC, working closely with the MERRA-2 science team, has created and is continuing to work on value added data services to best meet the needs of a broad user community. This presentation, using aerosol over Asia Monsoon as an example, provides an overview of the MERRA-2 data services at GES DISC, including: How to find the data? How many data access methods are provided? What are the best data access methods for me? How do download the subsetted (parameter, spatial, temporal) data and save in preferred spatial resolution and data format? How to visualize and explore the data online? In addition, we introduce a future online analytic tool designed for supporting application research, focusing on long-term hourly time-series data access and analysis.

Shen, Suhung↗
NASA NTRS2018

Conceptual Design of an Experiment to Study Dust Destruction by Astrophysical Shock Waves

A novel laboratory experimental design is described that will investigate the processing of dust grains in astrophysical shocks. Dust is a ubiquitous ingredient in the interstellar medium (ISM) of galaxies; however, its evolutionary cycle is still poorly understood. Especially shrouded in mystery is the efficiency of grain destruction by astrophysical shocks generated by expanding supernova remnants. While the evolution of these remnants is fairly well understood, the grain destruction efficiency in these shocks is largely unknown. The experiments described herein will fill this knowledge gap by studying the dust destruction efficiencies for shock velocities in the range of approximately 10-30 kilometers per second (microns per nanosecond), at which most of the grain destruction and processing in the ISM takes place. The experiments focus on the study of grain-grain collisions by accelerating small (1 millimeter) dust particles into a large (approximately 5-10 millimeter diameter) population; this simulates the astrophysical system well in that the more numerous, small grains impact and collide with the large population. Facilities that combine the versatility of high-power optical lasers with the diagnostic capabilities of X-ray free-electron lasers, e.g., the Matter in Extreme Conditions instrument at the SLAC (originally named Stanford Linear Accelerator Center) National Accelerator Laboratory, provide an ideal laboratory environment to create and diagnose dust destruction by astrophysically relevant shocks at the micron scale.

Manuel, M. J.-E.↗
NASA NTRS1993

Mapping the stability region of the 3:2 Neptune-Pluto resonance

Pluto and Charon are most likely the remnants of a large number of objects that existed in the Uranus-Neptune region at early epochs of the solar system. Numerical integrations have shown that, in general, such objects were ejected from the planetary region on timescales of approximately 10(exp 7) years after Neptune and Uranus reached their current masses. It is thought that the Pluto-Charon system survived to current times without being dynamically removed in this way because it is trapped in a set of secular and mean motion resonances with Neptune. The best-known Pluto-Neptune orbit coupling is the 3:2 mean motion resonance discovered almost 30 years ago by C. Cohen and E. Hubbard. These workers showed that the resonance angle, delta is equivalent to 3(lambda(sub P)) - 2(lambda(sub N)) - omega-bar(sub P) where omega-bar(sub P) is the longitude of perihelion of the Pluto-Charon system, and lambda(sub N) and lambda(sub P) are the mean longitude of Neptune and Pluto-Charon respectively, librates about 180 deg with an amplitude, A(sub delta), of 76 deg. A numerical simulation project to map out the stability region of the 3:2 resonance is reported. The results of these simulations are important to understanding whether Pluto's long-term heliocentric stability requires only the 3:2 resonance, or whether it instead requires one or more of the other Pluto-Neptune resonances. Our study also has another important application. By investigating stability timescales as a function of orbital elements, we gain insight into the fraction of orbital phase space which the stable 3:2 resonance occupies. This fraction is directly related to the probability that the Pluto-Charon system (and possibly other small bodies) could have been captured into this resonance.

Levison, H. F.↗
NASA NTRS1995

Numerical Solution of Multi-Dimensional Hyperbolic Conservation Laws on Unstructured Meshes

The lecture material will discuss the application of one-dimensional approximate Riemann solutions and high order accurate data reconstruction as building blocks for solving multi-dimensional hyperbolic equations. This building block procedure is well-documented in the nationally available literature. The relevant stability and convergence theory using positive operator analysis will also be presented. All participants in the minisymposium will be asked to solve one or more generic test problems so that a critical comparison of accuracy can be made among differing approaches.

Barth, Timothy J.↗
DOE OSTI2022

CFD-trained ANN Model for Approximating Near-occupant Condition in Real-time Simulations

The main drawback of Computational Fluid Dynamics (CFD) simulations has been the time and resource consuming nature which is not suitable for real-time applications. In this work, we first generated numerous CFD models of a given indoor space to obtain airspeed, temperature, and mean radiant temperature near an occupant as training data. Several artificial neural networks (ANN) models were trained using this CFD simulated data to approximate near real-time environmental conditions for a given occupant. This trained ANN model approach is a part of a real-time simulation of building operations using a combination of software and real hardware (HVAC equipment) approaches. The preliminary results suggest that the CFD- generated training data and the trained ANN model can accurately approximate such conditions in a real-time application, a method that has great potential in building simulation and building digital twin areas of research.

Zhang, Yun↗
NASA NTRS1982

Stability analysis of flexible wind turbine blades using finite element method

Static vibration and flutter analysis of a straight elastic axis blade was performed based on a finite element method solution. The total potential energy functional was formulated according to linear beam theory. The inertia and aerodynamic loads were formulated according to the blade absolute acceleration and absolute velocity vectors. In vibration analysis, the direction of motion of the blade during the first out-of-lane and first in-plane modes was examined; numerical results involve NASA/DOE Mod-0, McCauley propeller, north wind turbine and flat plate behavior. In flutter analysis, comparison cases were examined involving several references. Vibration analysis of a nonstraight elastic axis blade based on a finite element method solution was performed in a similar manner with the straight elastic axis blade, since it was recognized that a curved blade can be approximated by an assembly of a sufficient number of straight blade elements at different inclinations with respect to common system of axes. Numerical results involve comparison between the behavior of a straight and a curved cantilever beam during the lowest two in-plane and out-of-plane modes.

Kamoulakos, A.↗
NASA NTRS2005

Tools for Designing and Analyzing Structures

Structural Design and Analysis Toolset is a collection of approximately 26 Microsoft Excel spreadsheet programs, each of which performs calculations within a different subdiscipline of structural design and analysis. These programs present input and output data in user-friendly, menu-driven formats. Although these programs cannot solve complex cases like those treated by larger finite element codes, these programs do yield quick solutions to numerous common problems more rapidly than the finite element codes, thereby making it possible to quickly perform multiple preliminary analyses - e.g., to establish approximate limits prior to detailed analyses by the larger finite element codes. These programs perform different types of calculations, as follows: 1. determination of geometric properties for a variety of standard structural components; 2. analysis of static, vibrational, and thermal- gradient loads and deflections in certain structures (mostly beams and, in the case of thermal-gradients, mirrors); 3. kinetic energies of fans; 4. detailed analysis of stress and buckling in beams, plates, columns, and a variety of shell structures; and 5. temperature dependent properties of materials, including figures of merit that characterize strength, stiffness, and deformation response to thermal gradients

Luz, Paul L.↗
NASA NTRS1996

A wavelet-optimized, very high order adaptive grid and order numerical method

Differencing operators of arbitrarily high order can be constructed by interpolating a polynomial through a set of data followed by differentiation of this polynomial and finally evaluation of the polynomial at the point where a derivative approximation is desired. Furthermore, the interpolating polynomial can be constructed from algebraic, trigonometric, or, perhaps exponential polynomials. This paper begins with a comparison of such differencing operator construction. Next, the issue of proper grids for high order polynomials is addressed. Finally, an adaptive numerical method is introduced which adapts the numerical grid and the order of the differencing operator depending on the data. The numerical grid adaptation is performed on a Chebyshev grid. That is, at each level of refinement the grid is a Chebvshev grid and this grid is refined locally based on wavelet analysis.

Jameson, Leland↗
NASA NTRS1972

Current drainage to a high voltage probe in a dilute plasma.

The current drainage from a plasma through approximately 0.05 cm diameter holes in Kapton H polymide film, FEP type C, and quartz (7940 Corning glass) placed on a probe at voltages up to 2000 volts dc have been determined both analytically and experimentally. A Laplace field was used to numerically predict an upper limit for the drainage current. The measured current was less than the calculated current for the FEP and quartz, and the two were approximately in agreement for Kapton H. The measured currents were on the order of 5 microamperes with the calculated currents on the order of 10 microamperes.

Grier, N. T.↗
NASA NTRS1991

An accelerated lambda iteration method for multilevel radiative transfer. I - Non-overlapping lines with background continuum

A method is presented for solving multilevel transfer problems when nonoverlapping lines and background continuum are present and active continuum transfer is absent. An approximate lambda operator is employed to derive linear, 'preconditioned', statistical-equilibrium equations. A method is described for finding the diagonal elements of the 'true' numerical lambda operator, and therefore for obtaining the coefficients of the equations. Iterations of the preconditioned equations, in conjunction with the transfer equation's formal solution, are used to solve linear equations. Some multilevel problems are considered, including an eleven-level neutral helium atom. Diagonal and tridiagonal approximate lambda operators are utilized in the problems to examine the convergence properties of the method, and it is found to be effective for the line transfer problems.

Rybicki, G. B.↗
NASA NTRS1983

Uniform asymptotic approximations for duct eigenfunctions in a thin boundary layer flow

Analytical approximations for the acoustic modes in a duct carrying a uniform core flow with a thin shear layer at the walls are developed using the Method of Matched Asymptotic Expansions. Both two-dimensional and cylindrical duct propagation are considered. Numerical results for eigenvalues calculated using the theory are presented for the two dimensional problem and compared with results from earlier analyses. It is found that the new approximations yield a significant increase in accuracy.

Myers, M. K.↗
DOE OSTI2020

Approximate inverse-based block preconditioners in poroelasticity

We focus on the fully implicit solution of the linear systems arising from a three-field mixed finite element approximation of Biot’s poroleasticity equations. The objective is to develop algebraic block preconditioners for the efficient solution of such systems by Krylov subspace methods. In this work, we investigate the use of approximate inverse-based techniques to decouple the native system of equations and obtain explicit sparse approximations of the Schur complements related to the physics-based partitioning of the unknowns by field type. Here, the proposed methods are tested in various numerical experiments including real-world applications dealing with petroleum and geotechnical engineering.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Applications of computer algebra to distributed parameter systems

In the analysis of vibrations of continuous elastic systems, one often encounters complicated transcendental equations with roots directly related to the system's natural frequencies. Typically, these equations contain system parameters whose values must be specified before a numerical solution can be obtained. The present paper presents a method whereby the fundamental frequency can be obtained in analytical form to any desired degree of accuracy. The method is based upon truncation of rapidly converging series involving inverse powers of the system natural frequencies. A straightforward method to developing these series and summing them in closed form is presented. It is demonstrated how Computer Algebra can be exploited to perform the intricate analytical procedures which otherwise would render the technique difficult to apply in practice. We illustrate the method by developing two analytical approximations to the fundamental frequency of a vibrating cantilever carrying a rigid tip body. The results are compared to the numerical solution of the exact (transcendental) frequency equation over a range of system parameters.

Storch, Joel A.↗
NASA NTRS1991

An Evaluation of the Iosipescu Specimen for Composite Materials Shear Property Measurement

A detailed evaluation of the suitability of the Iosipescu specimen tested in the modified Wyoming fixture is presented. An experimental investigation using conventional strain gage instrumentation and moire interferometry is performed. A finite element analysis of the Iosipescu shear test for unidirectional and cross-ply composites is used to assess the uniformity of the shear stress field in the vicinity of the notch, and demonstrate the effect of the nonuniform stress field upon the strain gage measurements used for the determination of composite shear moduli. From the test results for graphite-epoxy laminates, it is shown that the proximity of the load introduction point to the test section greatly influences the individual gage readings for certain fiber orientations but the effect upon shear modulus measurement is relatively unimportant. A numerical study of the load contact effect shows the sensitivity of some fiber configurations to the specimen/fixture contact mechanism and may account for the variations in the measured shear moduli. A comparison of the strain gage readings from one surface of a specimen with corresponding data from moire interferometry on the opposite face documented an extreme sensitivity of some fiber orientations to eccentric loading which induced twisting and yielded spurious shear stress-strain curves. In the numerical analysis, it is shown that the Iosipescu specimens for different fiber orientations have to be modeled differently in order to closely approximate the true loading conditions. Correction factors are needed to allow for the nonuniformity of the strain field and the use of the average shear stress in the shear modulus evaluation. The correction factors, which are determined for the region occupied by the strain gage rosette, are found to be dependent upon the material orthotropic ratio and the finite element models. Based upon the experimental and numerical results, recommendations for improving the reliability and accuracy of the shear modulus values are made, and the implications for shear strength measurement discussed. Further application of the Iosipescu shear test to woven fabric composites is presented. The limitations of the traditional strain gage instrumentation on the satin weave and high tow plain weave fabrics is discussed. Test results of a epoxy based aluminum particulate composite is also presented. A modification of the Iosipescu specimen is proposed and investigated experimentally and numerically. It is shown that the proposed new specimen design provides a more uniform shear stress field in the test section and greatly reduces the normal and shear stress concentrations in the vicinity of the notches. While the fabrication and the material cost of the proposed specimen is tremendously reduced, it is shown the accuracy of the shear modulus measurement is not sacrificed.

Ho, Henjen↗
DOE OSTI2021

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive “full-wave” simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. Furthermore, the numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Non-numeric computation for high eccentricity orbits

Geocentric orbits of large eccentricity (e = 0.9 to 0.95) are significantly perturbed in cislunar space by the sun and moon. The time-history of the height of perigee, subsequent to launch, is particularly critical. The determination of 'launch windows' is mostly concerned with preventing the height of perigee from falling below its low initial value before the mission lifetime has elapsed. Between the extremes of high accuracy digital integration of the equations of motion and of using an approximate, but very fast, stability criteria method, this paper is concerned with the developement of a method of intermediate complexity using non-numeric computation. The computer is used as the theory generator to generalize Lidov's theory using six osculating elements. Symbolic integration is completely automatized and the output is a set of condensed formulae well suited for repeated applications in launch window analysis. Examples of applications are given.

Sridharan, R.↗