Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear extrapolation method”

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.

88 records · Page 5

Linear magnetization feature associated with Rima Sirsalis

Measurements of the magnetic fields by the electron reflection method in the neighborhood of the long structural rille Rima Sirsalis show that a magnetic field with a strength of at least 100 nT (100 gammas) is present over a region of the order of 10 km in width and at least 300 km long. The center of the magnetized region closely parallels and is centered on the rille. The linear magnetization feature extends at least to latitude 8 deg S, 60 km beyond the place where the rille disappears at the edge of Oceanus Procellarum. This extension is coincident with the extrapolation of the rille based on photographs. However, the magnetization is much weaker or entirely absent at 5 deg S and has vanished at 0 deg latitude. These results suggest that the rille is indeed a structural feature and has associated with it magnetization, either in the form of intrusive magnetized rock or of a gap in an otherwise more or less uniformly magnetized layer of rock of large extent in two dimensions. Furthermore, the rille structure evidently is present for some distance beneath the lava flows of the Oceanus Procellarum basin.

Anderson, K. A.↗

Numerical Study of Base Pressure Characteristic Curve for a Four-Engine Clustered Nozzle Configuration

Excessive base heating has been a problem for many launch vehicles. For certain designs such as the direct dump of turbine exhaust in the nozzle section and at the nozzle lip of the Space Transportation Systems Engine (STME), the potential burning of the turbine exhaust in the base region has caused tremendous concern. Two conventional approaches have been considered for predicting the base environment: (1) empirical approach, and (2) experimental approach. The empirical approach uses a combination of data correlations and semi-theoretical calculations. It works best for linear problems, simple physics and geometry. However, it is highly suspicious when complex geometry and flow physics are involved, especially when the subject is out of historical database. The experimental approach is often used to establish database for engineering analysis. However, it is qualitative at best for base flow problems. Other criticisms include the inability to simulate forebody boundary layer correctly, the interference effect from tunnel walls, and the inability to scale all pertinent parameters. Furthermore, there is a contention that the information extrapolated from subscale tests with combustion is not conservative. One potential alternative to the conventional methods is computational fluid dynamics (CFD), which has none of the above restrictions and is becoming more feasible due to maturing algorithms and advancing computer technology. It provides more details of the flowfield and is only limited by computer resources. However, it has its share of criticisms as a predictive tool for base environment. One major concern is that CFD has not been extensively tested for base flow problems. It is therefore imperative that CFD be assessed and benchmarked satisfactorily for base flows. In this study, the turbulent base flowfield of a experimental investigation for a four-engine clustered nozzle is numerically benchmarked using a pressure based CFD method. Since the cold air was the medium, accurate prediction of the base pressure distributions at high altitudes is the primary goal. Other factors which may influence the numerical results such as the effects of grid density, turbulence model, differencing scheme, and boundary conditions are also being addressed. Preliminary results of the computed base pressure agreed reasonably well with that of the measurement. Basic base flow features such as the reverse jet, wall jet, recompression shock, and static pressure field in plane of impingement have been captured.

Wang, Ten-See↗

3-D inelastic analysis methods for hot section components (base program)

A 3-D inelastic analysis methods program consists of a series of computer codes embodying a progression of mathematical models (mechanics of materials, special finite element, boundary element) for streamlined analysis of combustor liners, turbine blades, and turbine vanes. These models address the effects of high temperatures and thermal/mechanical loadings on the local (stress/strain) and global (dynamics, buckling) structural behavior of the three selected components. These models are used to solve 3-D inelastic problems using linear approximations in the sense that stresses/strains and temperatures in generic modeling regions are linear functions of the spatial coordinates, and solution increments for load, temperature and/or time are extrapolated linearly from previous information. Three linear formulation computer codes, referred to as MOMM (Mechanics of Materials Model), MHOST (MARC-Hot Section Technology), and BEST (Boundary Element Stress Technology), were developed and are described.

Wilson, R. B.↗

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗

Development of MCAERO wing design panel method with interactive graphics module

A reliable and efficient iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical pressure distribution. The design process is initialized by using MCAERO (MCAIR 3-D Subsonic Potential Flow Analysis Code) to analyze a baseline configuration. A second program DMCAERO is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter by applying a first-order expansion to the baseline equations in MCAERO. This matrix is calculated only once but is used in each iteration cycle to calculate the geometry perturbation and to analyze the perturbed geometry. The potential on the new geometry is calculated by linear extrapolation from the baseline solution. This extrapolated potential is converted to velocity by numerical differentiation, and velocity is converted to pressure by using Bernoulli's equation. There is an interactive graphics option which allows the user to graphically display the results of the design process and to interactively change either the geometry or the prescribed pressure distribution.

Hawk, J. D.↗

Determining The Nanoflare Heating Frequency of an X-Ray Bright Point Observed by MaGIXS

Nanoflares are thought to be one of the prime candidates that can heat the solar corona to its multi-million kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing simulated nanoflare-heated plasma emissions with the observed emission. Using HYDRAD coronal loop simulations, we model the emission from an X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with nearest-available observations from the Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO) and X-Ray Telescope (XRT) onboard Hinode observatory. The length and magnetic field strength of the coronal loops are derived from the linear-force-free extrapolation of the observed photospheric magnetogram by Helioseismic and Magnetic Imager (HMI) onboard SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the loop length and magnetic field strength. The simulation results are then compared and matched against the measured intensity from AIA, XRT, and MaGIXS. Our model results indicate the observed emissions from the XBP under study could be well matched by a distribution of nanoflares with average delay times 1500 s to 3000 s, which suggest that the heating is dominated by high-frequency events. Further, we demonstrate the high sensitivity of MaGIXS and XRT to diagnose the heating frequency using this method, while AIA passbands are found to be the least sensitive.

coronal heating↗

Exploring the Relationship Between Reflectance Red Edge and Chlorophyll Content in Slash Pine

Chlorophyll is a key indicator of the physiological status of a forest canopy. However, its distribution may vary greatly in time and space, so that the estimation of chlorophyll content of canopies or branches by extrapolation from leaf values obtained by destructive sampling is labor intensive and potentially inaccurate. Chlorophy11 content is related positively to the point of maximum slope in vegetation reflectance spectra which occurs at wavelengths between 690-740 nm and is known as the "red edge." The red edge of needles on individual slash pine (Piniis elliottii Engelm.) branches and in whole forest canopies was measured with a spectroradiometer. Branches were measured on the ground against a spectrally flat reflectance target and canopies were measured from observation towers against a spectrally variable understory and forest floor. There was a linear relationship between red edge and chlorophyll content of branches (R(exp 2) = 0.91). Measurements of the red edge and this relationship were used to estimate the chlorophyll content of other branches with an error that was lower than that associated with the colorimetric (laboratory) method. There was no relationship between the red edge and the chlorophyll content of whole canopies. This can be explained by the overriding influence of the understory and forest floor, an influence that was illustrated by spectral mixture modeling. The results suggest that the red edge could be used to estimate the chlorophyll content in branches but it is unlikely to be of value for the estimation of chlorophyll content in canopies unless the canopy cover is high.

Curran, Paul J.↗

An investigation of new methods for estimating parameter sensitivities

Parameter sensitivity is defined as the estimation of changes in the modeling functions and the design variables due to small changes in the fixed parameters of the formulation. There are currently several methods for estimating parameter sensitivities requiring either difficult to obtain second order information, or do not return reliable estimates for the derivatives. Additionally, all the methods assume that the set of active constraints does not change in a neighborhood of the estimation point. If the active set does in fact change, than any extrapolations based on these derivatives may be in error. The objective here is to investigate more efficient new methods for estimating parameter sensitivities when the active set changes. The new method is based on the recursive quadratic programming (RQP) method and in conjunction a differencing formula to produce estimates of the sensitivities. This is compared to existing methods and is shown to be very competitive in terms of the number of function evaluations required. In terms of accuracy, the method is shown to be equivalent to a modified version of the Kuhn-Tucker method, where the Hessian of the Lagrangian is estimated using the BFS method employed by the RPQ algorithm. Inital testing on a test set with known sensitivities demonstrates that the method can accurately calculate the parameter sensitivity. To handle changes in the active set, a deflection algorithm is proposed for those cases where the new set of active constraints remains linearly independent. For those cases where dependencies occur, a directional derivative is proposed. A few simple examples are included for the algorithm, but extensive testing has not yet been performed.

Beltracchi, Todd J.↗

Thermal Conductance of Pressed Bimetal Contact Pairs at Liquid Nitrogen Temperatures

Large Dewars often use aluminum radiation shields and stainless steel vent lines. A simple, low cost method of making thermal contact between the shield and the line is to deform the shield around the line. A knowledge of the thermal conductance of such a joint is needed to thermally analyze the system. The thermal conductance of pressed metal contacts consisting of one aluminum and one stainless steel contact has been measured at 77 K, with applied forces from 8.9 N to 267 N. Both 5052 or 5083 aluminum were used as the upper contact. The lower contact was 304L stainless steel. The thermal conductance was found to be linear in temperature over the narrow temperature range of measurement. As the force was increased, the thermal conductance ranged from roughly 9 to 21 mW/K within a range of errors from 3% to 8%. Within the range of error no difference could be found between the using either of the aluminum alloys as the upper contact. Extrapolating the data to zero applied force does not result in zero thermal conductance. Possible causes of this anomalous effect are discussed.

Kittle, Peter↗

Consistent Pl Analysis of Aqueous Uranium-235 Critical Assemblies

The lethargy-dependent equations of the consistent Pl approximation to the Boltzmann transport equation for slowing down neutrons have been used as the basis of an IBM 704 computer program. Some of the effects included are (1) linearly anisotropic center of mass elastic scattering, (2) heavy element inelastic scattering based on the evaporation model of the nucleus, and (3) optional variation of the buckling with lethargy. The microscopic cross-section data developed for this program covered 473 lethargy points from lethargy u = 0 (10 Mev) to u = 19.8 (0.025 ev). The value of the fission neutron age in water calculated here is 26.5 square centimeters; this value is to be compared with the recent experimental value given as 27.86 square centimeters. The Fourier transform of the slowing-down kernel for water to indium resonance energy calculated here compared well with the Fourier transform of the kernel for water as measured by Hill, Roberts, and Fitch. This method of calculation has been applied to uranyl fluoride - water solution critical assemblies. Theoretical results established for both unreflected and fully reflected critical assemblies have been compared with available experimental data. The theoretical buckling curve derived as a function of the hydrogen to uranium-235 atom concentration for an energy-independent extrapolation distance was successful in predicting the critical heights of various unreflected cylindrical assemblies. The critical dimensions of fully water-reflected cylindrical assemblies were reasonably well predicted using the theoretical buckling curve and reflector savings for equivalent spherical assemblies.

Fieno, Daniel↗

Design Process for High Speed Civil Transport Aircraft Improved by Neural Network and Regression Methods

A key challenge in designing the new High Speed Civil Transport (HSCT) aircraft is determining a good match between the airframe and engine. Multidisciplinary design optimization can be used to solve the problem by adjusting parameters of both the engine and the airframe. Earlier, an example problem was presented of an HSCT aircraft with four mixed-flow turbofan engines and a baseline mission to carry 305 passengers 5000 nautical miles at a cruise speed of Mach 2.4. The problem was solved by coupling NASA Lewis Research Center's design optimization testbed (COMETBOARDS) with NASA Langley Research Center's Flight Optimization System (FLOPS). The computing time expended in solving the problem was substantial, and the instability of the FLOPS analyzer at certain design points caused difficulties. In an attempt to alleviate both of these limitations, we explored the use of two approximation concepts in the design optimization process. The two concepts, which are based on neural network and linear regression approximation, provide the reanalysis capability and design sensitivity analysis information required for the optimization process. The HSCT aircraft optimization problem was solved by using three alternate approaches; that is, the original FLOPS analyzer and two approximate (derived) analyzers. The approximate analyzers were calibrated and used in three different ranges of the design variables; narrow (interpolated), standard, and wide (extrapolated).

Hopkins, Dale A.↗

Human Tracking Performance in Uncoupled and Coupled Two-Axis Systems

This report presents tile results of an experimental and analytical study of human performance in uncoupled and coupled control systems. Human pilot performance in single and two-axis systems was mathematically modeled by linear second-order describing functions. Model parameters were determined using model matching techniques. Analysis of the models showed that the amplitude ratio and phase lead of the describing function increased with training indicating an increase in open loop bandwidth. The phase margin also decreased with training. Increasing the plant lag time constant resulted in an increase in the model lead time constant and a decrease in the zero frequency gain. No significant difference was found to exist in the normalized tracking error per axis between the two-axis tasks and the single-axis tasks. However tile model lead time constant was significantly greater in two-axis tracking. Manual tracking of two-axis systems with cross-coupling was studied experimentally and analytically. Approximate methods for modeling two-axis performance were developed and checked using a precise spectral analysis approach. Coupled and uncoupled, symmetrical and asymmetrical two-axis performance was compared. The results show that modeling of cross-coupled systems is feasible and that trained subjects are capable of decoupling the axes of some systems. A methodology study compared the identification performance of continuous, iterative, and extrapolation model matching techniques. An iterative technique employing sensitivity equations for the generation of influence coefficients was found to be the best technique due to its excellent identification accuracy and ease of implementation. Convergence in iterative techniques can be improved substantially by equalizing the parameter adjustment rates and limiting the maximum correction per iteration.

Todosiev, E. P.↗

Elevated temperature viscous remanent magnetization of natural and synthetic multidomain magnetite

The time-temperature relationship of magnetization is a subject of much interest and debate by paleomagnetists, rock magnetists, and magnetic anomaly modellers. We have investigated this relationship by studying the viscous remanent magnetization (VRM) of coarse-grained multidomain (MD) magnetite. Our experiments covered the temperature range from 22 to 400 C, times from minutes to months, and included both Australian granulite samples and multidomain magnetite samples synthesized by the glass ceramic method. VRM acquisition was found to generally increase with temperature but not always at the rate predicted from classical thermal fluctuation theories. Thermal cycling between room temperature (at which the measurements were made) and the VRM acquisition temperature sharply decreased the temperature dependence of the VRM. Room temperature VRM acquisition accelerates with time when plotted on a semilog plot, whereas at elevated temperature the curves are quasilinear against log(time) for both the natural and synthetic samples. This change in behavior may suggest a variation in the VRM acquisition mechanism as a function of temperature for MD magnetite. The granulites have a nearly linear increase in VRM acquisition rate with temperature whereas the glass ceramics display little change in the acquisition rate between 22 and 200 C, but increase by nearly a factor of 3 by 400 C. The increase in VRM of the glass ceramics between 200 and 400 C is in general qualitative agreement with thermal fluctuation theory. There was no systematic change in the rate of VRM acquisition with grain size for the multidomain magnetites used in this study. Elevated temperature (e.g., 400 C) VRM acquisition by the deep crustal granulites, if extrapolated over the Brunhes chron, would produce a magnetization of several A/m which, if true, is of the order required by models for the source of long-wavelength magnetic anomalies.

Kelso, Paul R.↗

Investigation of HZETRN 2010 as a Tool for Single Event Effect Qualification of Avionics Systems

NASA's future missions are focused on deep space for human exploration that do not provide a simple emergency return to Earth. In addition, the deep space environment contains a constant background Galactic Cosmic Ray (GCR) radiation exposure, as well as periodic Solar Particle Events (SPEs) that can produce intense amounts of radiation in a short amount of time. Given these conditions, it is important that the avionics systems for deep space human missions are not susceptible to Single Event Effects (SEE) that can occur from radiation interactions with electronic components. The typical process to minimizing SEE effects is through using heritage hardware and extensive testing programs that are very costly. Previous work by Koontz, et al. [1] utilized an analysis-based method for investigating electronic component susceptibility. In their paper, FLUKA, a Monte Carlo transport code, was used to calculate SEE and single event upset (SEU) rates. This code was then validated against in-flight data. In addition, CREME-96, a deterministic code, was also compared with FLUKA and in-flight data. However, FLUKA has a long run-time (on the order of days), and CREME-96 has not been updated in several years. This paper will investigate the use of HZETRN 2010, a deterministic transport code developed at NASA Langley Research Center, as another tool that can be used to analyze SEE and SEU rates. The benefits to using HZETRN over FLUKA and CREME-96 are that it has a very fast run time (on the order of minutes) and has been shown to be of similar accuracy as other deterministic and Monte Carlo codes when considering dose [2, 3, 4]. The 2010 version of HZETRN has updated its treatment of secondary neutrons and thus has improved its accuracy over previous versions. In this paper, the Linear Energy Transfer (LET) spectra are of interest rather than the total ionizing dose. Therefore, the LET spectra output from HZETRN 2010 will be compared with the FLUKA and in-flight data to validate HZETRN 2010 as a computational tool for SEE qualification by analysis. Furthermore, extrapolation of these data to interplanetary environments at 1 AU will be investigated to determine whether HZETRN 2010 can be used successfully and confidently for deep space mission analyses.

Rojdev, Kristina↗

Scaling effects in the static and dynamic response of graphite-epoxy beam-columns

Scale model technology represents one method of investigating the behavior of advanced, weight-efficient composite structures under a variety of loading conditions. It is necessary, however, to understand the limitations involved in testing scale model structures before the technique can be fully utilized. These limitations, or scaling effects, are characterized. in the large deflection response and failure of composite beams. Scale model beams were loaded with an eccentric axial compressive load designed to produce large bending deflections and global failure. A dimensional analysis was performed on the composite beam-column loading configuration to determine a model law governing the system response. An experimental program was developed to validate the model law under both static and dynamic loading conditions. Laminate stacking sequences including unidirectional, angle ply, cross ply, and quasi-isotropic were tested to examine a diversity of composite response and failure modes. The model beams were loaded under scaled test conditions until catastrophic failure. A large deflection beam solution was developed to compare with the static experimental results and to analyze beam failure. Also, the finite element code DYCAST (DYnamic Crash Analysis of STructure) was used to model both the static and impulsive beam response. Static test results indicate that the unidirectional and cross ply beam responses scale as predicted by the model law, even under severe deformations. In general, failure modes were consistent between scale models within a laminate family; however, a significant scale effect was observed in strength. The scale effect in strength which was evident in the static tests was also observed in the dynamic tests. Scaling of load and strain time histories between the scale model beams and the prototypes was excellent for the unidirectional beams, but inconsistent results were obtained for the angle ply, cross ply, and quasi-isotropic beams. Results show that valuable information can be obtained from testing on scale model composite structures, especially in the linear elastic response region. However, due to scaling effects in the strength behavior of composite laminates, caution must be used in extrapolating data taken from a scale model test when that test involves failure of the structure.

Jackson, Karen E.↗

Metrology of Large Parts

As discussed in the first chapter of this book, there are many different methods to measure a part using optical technology. Chapter 2 discussed the use of machine vision to measure macroscopic features such as length and position, which was extended to the use of interferometry as a linear measurement tool in chapter 3, and laser or other trackers to find the relation of key points on large parts in chapter 4. This chapter looks at measuring large parts to optical tolerances in the sub-micron range using interferometry, ranging, and optical tools discussed in the previous chapters. The purpose of this chapter is not to discuss specific metrology tools (such as interferometers or gauges), but to describe a systems engineering approach to testing large parts. Issues such as material warpage and temperature drifts that may be insignificant when measuring a part to micron levels under a microscope, as will be discussed in later chapters, can prove to be very important when making the same measurement over a larger part. In this chapter, we will define a set of guiding principles for successfully overcoming these challenges and illustrate the application of these principles with real world examples. While these examples are drawn from specific large optical testing applications, they inform the problems associated with testing any large part to optical tolerances. Manufacturing today relies on micrometer level part performance. Fields such as energy and transportation are demanding higher tolerances to provide increased efficiencies and fuel savings. By looking at how the optics industry approaches sub-micrometer metrology, one can gain a better understanding of the metrology challenges for any larger part specified to micrometer tolerances. Testing large parts, whether optical components or precision structures, to optical tolerances is just like testing small parts, only harder. Identical with what one does for small parts, a metrologist tests large parts and optics in particular to quantify their mechanical properties (such as dimensions, mass, etc); their optical prescription or design (i.e. radius of curvature, conic constant, vertex location, size); and their full part shape. And, just as with small parts, a metrologist accomplishes these tests using distance measuring instruments such as tape measures, inside micrometers, coordinate measuring machines, distance measuring interferometers; angle measuring instruments such as theodolites, autocollimators; and surface measuring instruments including interferometers, stylus profilers, interference microscopes, photogrammetric cameras, or other tools. However, while the methodology may be similar, it is more difficult to test a large object for the simple reason that most metrologists do not have the necessary intuition. The skills used to test small parts or optics in a laboratory do not extrapolate to testing large parts in an industrial setting any more than a backyard gardener might successfully operate a farm. But first, what is a large part? A simple definition might be the part's size or diameter. For optics and diffuse surface parts alike, the driving constraint is ability to illuminate the part's surface. For reflective convex mirrors, large is typically anything greater than 1 meter. But, for refractive optics, flats or convex mirrors, large is typically greater than 0.5 meter. While a size definition is simple, it may be less than universal. A more nuanced definition might be that a large part is any component which cannot be easily tested in a standard laboratory environment, on a standard vibration isolated table using standard laboratory infrastructure. A micro-switch or a precision lens might be easily measured to nanometer levels under a microscope in a lab, but a power turbine spline or a larger telescope mirror will not fit under that microscope and may not even fit on the table.

Stahl, H. Philip↗