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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Optimal feedback control of turbulent channel flow

Feedback control equations were developed and tested for computing wall normal control velocities to control turbulent flow in a channel with the objective of reducing drag. The technique used is the minimization of a 'cost functional' which is constructed to represent some balance of the drag integrated over the wall and the net control effort. A distribution of wall velocities is found which minimizes this cost functional some time shortly in the future based on current observations of the flow near the wall. Preliminary direct numerical simulations of the scheme applied to turbulent channel flow indicates it provides approximately 17 percent drag reduction. The mechanism apparent when the scheme is applied to a simplified flow situation is also discussed.

Bewley, Thomas↗

NASA NTRS1993

Antarctic Meteorite Location and Mapping Project (AMLAMP): Antarctic meteorite location map series explanatory text and user's guide to AMLAMP data

This technical report is an update to LPI Technical Report 89-02, which contained data and information that was current to May 1987. Since that time approximately 4000 new meteorites have been collected, mapped, and characterized, mainly from the numerous ice fields in the Allan Hills-David Glacier region, from the Pecora Escarpment and Moulton Escarpment in the Thiel Mountains-Patuxent region, the Wisconsin Range region, and from the Beardmore region. Meteorite location maps for ice fields from these regions have been produced and are available. This report includes explanatory texts for the maps of new areas and provides information on updates of maps of the areas covered in LPI Technical Report 89-02. Sketch maps and description of locales that have been searched and have yielded single or few meteorites are also included. The meteorite listings for all the ice fields have been updated to include any classification changes and new meteorites recovered from ice fields in the Allan Hills-David Glacier region since 1987. The text has been reorganized and minor errors in the original report have been corrected. Computing capabilities have improved immensely since the early days of this project. Current software and hardware allow easy access to data over computer networks. With various commercial software packages, the data can be used many different ways, including database creation, statistics, and mapping. The databases, explanatory texts, and the plotter files used to produce the meteorite location maps are available through a computer network. Information on how to access AMLAMP data, its formats, and ways it can be used are given in the User's Guide to AMLAMP Data section. Meteorite location maps and thematic maps may be ordered from the Lunar and Planetary Institute. Ordering information is given in Appendix A.

Schutt, J.↗

NASA NTRS1994

Analysis and control of hourglass instabilities in underintegrated linear and nonlinear elasticity

Methods are described to identify and correct a bad finite element approximation of the governing operator obtained when under-integration is used in numerical code for several model problems: the Poisson problem, the linear elasticity problem, and for problems in the nonlinear theory of elasticity. For each of these problems, the reason for the occurrence of instabilities is given, a way to control or eliminate them is presented, and theorems of existence, uniqueness, and convergence for the given methods are established. Finally, numerical results are included which illustrate the theory.

Jacquotte, Olivier P.↗

NASA NTRS1994

Recent deformation rates on Venus

Constraints on the recent geological evolution of Venus may be provided by quantitative estimates of the rates of the principal resurfacing processes, volcanism and tectonism. This paper focuses on the latter, using impact craters as strain indicators. The total postimpact tectonic strain lies in the range 0.5-6.5%, which defines a recent mean strain rate of 10(exp -18)-10(exp -17)/s when divided by the mean surface age. Interpretation of the cratering record as one of pure production requires a decline in resurfacing rates at about 500 Ma (catastrophic resurfacing model). If distributed tectonic resurfacing contributed strongly before that time, as suggested by the widespread occurrence of tessera as inliers, the mean global strain rate must have been at least approximately 10(exp -15)/s, which is also typical of terrestrial active margins. Numerical calculations of the response of the lithosphere to inferred mantle convective forces were performed to test the hypothesis that a decrease in surface strain rate by at least two orders of magnitude could be caused by a steady decline in heat flow over the last billion years. Parameterized convection models predict that the mean global thermal gradient decreases by only about 5 K/km over this time; even with the exponential dependence of viscosity upon temperature, the surface strain rate drops by little more than one order of magnitude. Strongly unsteady cooling and very low thermal gradients today are necessary to satisfy the catastrophic model. An alternative, uniformitarian resurfacing hypothesis holds that Venus is resurfaced in quasi-random 'patches' several hundred kilometers in size that occur in response to changing mantle convection patterns.

Grimm, Robert E.↗

NASA NTRS1994

Equivalent plate structural modeling for wing shape optimization including transverse shear

A new technique for structural modeling of airplanes wings is presented taking transverse shear effects into account. The kinematic assumptions of first-order shear deformation plate theory in combination with numerical analysis, where simple polynomials are used to define geometry, construction, and displacement approximations, lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps, and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of wing modeling techniques based on classical plate theory are discussed, and the improved accuracy of the new equivalent plate technique is demonstrated through comparison with finite element analysis and test results. Expressions for analytical derivatives of stiffness, mass, and load terms with respect to wing shape are given. Based on these, it is possible to obtain analytic sensitivities of displacements, stresses, and natural frequencies with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli↗

NASA NTRS1996

Heat Transfer Experiments in the Internal Cooling Passages of a Cooled Radial Turbine Rotor

An experimental study was conducted (1) to experimentally measure, assess and analyze the heat transfer within the internal cooling configuration of a radial turbine rotor blade and (2) to obtain heat transfer data to evaluate and improve computational fluid dynamics (CFD) procedures and turbulent transport models of internal coolant flows. A 1.15 times scale model of the coolant passages within the NASA LERC High Temperature Radial Turbine was designed, fabricated of Lucite and instrumented for transient beat transfer tests using thin film surface thermocouples and liquid crystals to indicate temperatures. Transient heat transfer tests were conducted for Reynolds numbers of one-fourth, one-half, and equal to the operating Reynolds number for the NASA Turbine. Tests were conducted for stationary and rotating conditions with rotation numbers in the range occurring in the NASA Turbine. Results from the experiments showed the heat transfer characteristics within the coolant passage were affected by rotation. In general, the heat transfer increased and decreased on the sides of the straight radial passages with rotation as previously reported from NASA-HOST-sponsored experiments. The heat transfer in the tri-passage axial flow region adjacent to the blade exit was relatively unaffected by rotation. However, the heat transfer on one surface, in the transitional region between the radial inflow passage and axial, constant radius passages, decreased to approximately 20 percent of the values without rotation. Comparisons with previous 3-D numerical studies indicated regions where the heat transfer characteristics agreed and disagreed with the present experiment.

Johnson, B. V.↗

NASA NTRS1996

Simplified Phased-Mission System Analysis for Systems with Independent Component Repairs

Accurate analysis of reliability of system requires that it accounts for all major variations in system's operation. Most reliability analyses assume that the system configuration, success criteria, and component behavior remain the same. However, multiple phases are natural. We present a new computationally efficient technique for analysis of phased-mission systems where the operational states of a system can be described by combinations of components states (such as fault trees or assertions). Moreover, individual components may be repaired, if failed, as part of system operation but repairs are independent of the system state. For repairable systems Markov analysis techniques are used but they suffer from state space explosion. That limits the size of system that can be analyzed and it is expensive in computation. We avoid the state space explosion. The phase algebra is used to account for the effects of variable configurations, repairs, and success criteria from phase to phase. Our technique yields exact (as opposed to approximate) results. We demonstrate our technique by means of several examples and present numerical results to show the effects of phases and repairs on the system reliability/availability.

Somani, Arun K.↗

NASA NTRS1994

Modeling Ionospheric Convection During a Major Geomagnetic Storm on October 22-23, 1981

Following the passage of an interplanetary shock at approximately 0500 UT, a major geomagnetic storm developed on October 22-23, 1981. Numerous auroral substorms occurred during this storm leading to an AE index greater than 1000 nT. We have used the expanding/contracting polar cap (ECPC) model (Moses et al., 1989) and data from the Dynamics Explorer 2 spacecraft to study the ionospheric electric fields for 12 consecutive traversals of the polar regions. The ECPC model can determine the voltage drops across the dayside merging and nightside reconnection gaps. We determined the relationship of the AL index (i.e., the intensity of the westward electrojet) to the nightside reconnection potential drop. An excellent linear correlation was found between the nightside reconnection gap voltage drop and the AL index. These results show that the solar wind strongly drives the magnetosphere-ionosphere system throughout the geomagnetic storm. A substantial level of dayside merging seems to occur throughout the event. Nightside reconnection varies from satellite pass to satellite pass and within the substorm recovery phase. We find that tail reconnection is an important feature of the recovery phase of substorms.

Moses, J. J.↗

NASA NTRS1997

The Geometric Spreading of Coronal Plumes and Coronal Holes

The geometric spreading in plumes and in the interplume region in coronal holes is calculated, using analytic and numerical theoretical models, between 1.0 and 5.0 solar radius. We apply two scale approximation that permits the rapid local spreading at the base of plumes to be evaluated separately from the global spreading imposed by coronal hole geometry. We show that fl can be computed from a potential field model and fg can be computed from global magnetohydrodynamic simulations of coronal structure. The approximations are valid when the plasma beta is small with respect to unity and for a plume separation small with respect to a solar radius.

Suess, S. T.↗

NASA NTRS1994

Equivalent Plate Structural Modeling for Wing Shape Optimization Including Transverse Shear

A new technique for structural modeling of airplane wings is presented taking transverse shear effects into account. The kinematic assumptions of first-order shear deformation plate theory In combination with numerical analysis, where simple polynomials are used to define geometry, construction, and displacement approximations, lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps, and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of wing modeling techniques based on classical plate theory are discussed, and the Improved accuracy of the new equivalent plate technique is demonstrated through comparison with finite element analysis and test results. Expressions for analytical derivatives of stiffness, mass, and load terms with respect to wing shape are given. Based on these, it is possible to obtain analytic sensitivities of displacements, stresses, and natural frequencies with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli↗

NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗

NASA NTRS2005

Variability of the ocean-induced magnetic field predicted at sea surface and at satellite altitudes

Spatial and temporal variability of the magnetic field component induced by ocean circulation is investigated on the basis of a standard thin-shell approximation of electro- and magneto-static equations. Well-known difficulties of numerical solution of the governing equations are resolved by reducing the problem to an equation for the electric field potential,(ef) as opposed to a more conventional approach focused on the vertical jump, (psi), of the magnetic field potential across a combined ocean/ marine-sediment-layer spherical shell. The present formulation permits using more realistic input data on ocean currents and ultimately yields much greater (by at least an order of magnitude) values of the magnetic field at sea surface than predicted in earlier studies. Such large values are comparable to, and in some cases exceed, magnetic field variations caused by lithospheric and ionospheric sources on monthly to interannual timescales. At the 400-km altitude (of CHAMP satellite), the field attains 6 nT. The model predictions show favorable comparisons with some in situ measurements as well as with Challenging Minisatellite Payload (CHAMP) satellite magnetometer data.

ocean induced magnetic field↗

NASA NTRS2011

Designing the Space Shuttle Propulsion System

The major elements of the Space Shuttle Main Propulsion System include two reusable solid rocket motors integrated into recoverable solid rocket boosters, an expendable external fuel and oxidizer tank, and three reusable Space Shuttle Main Engines. Both the solid rocket motors and space shuttle main engines ignite prior to liftoff, with the solid rocket boosters separating about two minutes into flight. The external tank separates after main engine shutdown and is safely expended in the ocean. The SSME's, integrated into the Space Shuttle Orbiter aft structure, are reused after post landing inspections. Both the solid rocket motors and the space shuttle main engine throttle during early ascent flight to limit aerodynamic loads on the structure. The configuration is called a stage and a half as all the propulsion elements are active during the boost phase, and the SSME's continue operation to achieve orbital velocity approximately eight and a half minutes after liftoff. Design and performance challenges were numerous, beginning with development work in the 1970 s. The solid rocket motors were large, and this technology had never been used for human space flight. The SSME s were both reusable and very high performance staged combustion cycle engines, also unique to the Space Shuttle. The multi body side mount configuration was unique and posed numerous integration and interface challenges across the elements. Operation of the system was complex and time consuming. This paper discusses a number of the system level technical challenges including development and operations.

Owen, James↗

NASA NTRS2022

Extension and Validation of the Pendulum Model for Longitudinal Solar Prominence Oscillations

Context. Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semicircular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not semicircular. Aims. Our aim is to understand the effects of including the nonuniform solar gravity on longitudinal oscillations and explore the validity of the pendulum model with different flux-tube geometries. Methods. We first derived the equation describing the motion of the plasma along the flux tube including the effects of nonuniform gravity, yielding corrections to the original pendulum model. We also computed the full numerical solutions for the normal modes and compared them with the new pendulum approximation. Results. We find that the nonuniform gravity introduces a significant modification in the pendulum model. We also found a cut-off period; i.e., the longitudinal oscillations cannot have a period longer than 167 min. In addition, considering different tube geometries, the period depends almost exclusively on the radius of curvature at the bottom of the dip. Conclusions. We conclude that nonuniform gravity significantly modifies the pendulum model. These corrections are important for prominence seismology, because the inferred values of the radius of curvature and minimum magnetic-field strength differ substantially from those of the old model. However, we find that the corrected pendulum model is quite robust and is still valid for noncircular dips.

Sun: corona↗

DOE OSTI2022

June 2022 Permeability Measurements in P-Tunnel, Nevada National Security Site

This document describes gas permeability measurements made in June 2022 in P-Tunnel at the Nevada National Security Site (NNSS). Four boreholes (AC-1, GI-5, GI-6, and HF-1) were tested with a dual packer system. The pressure versus flow rate data that were collected are analyzed using a numerical simulation method that captures compressible gas flow in partially saturated porous media using a close approximation to the geometry of the test setup. The first part of this document begins with data collection and a description of the packer assembly followed by details on data acquisition, measurement instrumentation, and equipment followed by a description of pre-test system verification. Next, the implementation of testing in the field is described, including how specific locations were chosen for testing. The second part of this document describes the numerical simulations used to back-out permeability from the pressure versus flow rate data. First, the numerical model is described, including assumptions. Next the details of the mesh that represents the experiment are described. A method for generating type curves in pressure versus flow rate space is presented and finally the data are mapped to permeability. We conclude with a discussion of the results.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Performance Evaluations of Noisy Approximate Quantum Fourier Arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route towards a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted sum optimization problems in data processing and machine learning and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance, and identify the most effective approximation depths for quantum add and quantum multiply within the given context. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges - relevant to signed, unsigned, modular and non-modular versions - that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze performance trends in our results and speculate on possible future development within this computational paradigm.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Performance evaluations of signed and unsigned noisy approximate quantum Fourier arithmetic

The Quantum Fourier Transform (QFT) grants competitive advantages, especially in resource usage and circuit approximation, for performing arithmetic operations on quantum computers, and offers a potential route toward a numerical quantum-computational paradigm. In this paper, we utilize efficient techniques to implement QFT-based integer addition and multiplications. These operations are fundamental to various quantum applications including Shor’s algorithm, weighted-sum optimization problems in data processing and machine learning, and quantum algorithms requiring inner products. We carry out performance evaluations of these implementations based on IBM’s superconducting-qubit architecture using different compatible noise models. We isolate the sensitivity of the component quantum circuits on both one-/two-qubit gate error rates, and the number of the arithmetic operands’ superposed integer states. We analyze performance and identify the most effective approximation depths for unsigned quantum addition and quantum multiplication within the given context. We then perform a similar analysis of signed addition and compare to the unsigned results. We observe significant dependency of the optimal approximation depth on the degree of machine noise and the number of superposed states in certain performance regimes. Finally, we elaborate on the algorithmic challenges—relevant to signed, unsigned, modular and non-modular versions—that could also be applied to current implementations of QFT-based subtraction, division, exponentiation, and their potential tensor extensions. Here, we analyze the performance trends in our results and speculate on possible future developments within this computational paradigm.

Computational models↗

NASA NTRS1974

Static, stability, and dynamic analysis of shells of revolution by numerical integration - A comparison

Recent innovations in digital computer technology have enabled engineers to analyze shell structures of complex configurations without unduly restrictive approximations. An attempt is made to compare the various programs now generally available from the point of view of the advantages of the relative technique utilized, as well as the programmed state of the art. Many of the comparisons are based on the sample problems solved by the STARS-2 system of programs. These examples indicate both the structural detail which can be analyzed by, and the analytical capabilities available in, the numerical shell-of-revolution programs. All advantages and differences are demonstrated by use of solutions for realistic shell problems in the areas of statics, stability, vibrations, and dynamic response of shells subjected to time-dependent loadings.

Svalbonas, V.↗