Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interior point 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.

At least 145 records · Page 8

Radiation/convection coupling in rocket motors and plumes

The three commonly used propellant systems - H2/O2, RP-1/O2, and solid propellants - primarily radiate as molecular emitters, non-scattering small particles, and scattering larger particles, respectively. Present technology has accepted the uncoupling of the radiation analysis from that of the flowfield. This approximation becomes increasingly inaccurate as one considers plumes, interior rocket chambers, and nuclear rocket propulsion devices. This study will develop a hierarchy of methods which will address radiation/convection coupling in all of the aforementioned propulsion systems. The nature of the radiation/convection coupled problem is that the divergence of the radiative heat flux must be included in the energy equation and that the local, volume-averaged intensity of the radiation must be determined by a solution of the radiative transfer equation (RTE). The intensity is approximated by solving the RTE along several lines of sight (LOS) for each point in the flowfield. Such a procedure is extremely costly; therefore, further approximations are needed. Modified differential approximations are being developed for this purpose. It is not obvious which order of approximations are required for a given rocket motor analysis. Therefore, LOS calculations have been made for typical rocket motor operating conditions in order to select the type approximations required. The results of these radiation calculations, and the interpretation of these intensity predictions are presented herein.

Farmer, R. C.↗

Multiple point least squares equalization in a room

Equalization filters designed to minimize the mean square error between a delayed version of the original electrical signal and the equalized response at a point in a room have previously been investigated. In general, such a strategy degrades the response at positions in a room away from the equalization point. A method is presented for designing an equalization filter by adjusting the filter coefficients to minimize the sum of the squares of the errors between the equalized responses at multiple points in the room and delayed versions of the original, electrical signal. Such an equalization filter can give a more uniform frequency response over a greater volume of the enclosure than can the single point equalizer above. Computer simulation results are presented of equalizing the frequency responses from a loudspeaker to various typical ear positions, in a room with dimensions and acoustic damping typical of a car interior, using the two approaches outlined above. Adaptive filter algorithms, which can automatically adjust the coefficients of a digital equalization filter to achieve this minimization, will also be discussed.

Elliott, S. J.↗

Asymptotic modal analysis and statistical energy analysis

The sound field of a structural-acoustic enclosure was subject to experimental analysis and theoretical description in order to develop an efficient and accurate method for predicting sound pressure levels in enclosures such as aircraft fuselages. Asymptotic Modal Analysis (AMA) is the method under investigation. AMA is derived from classical modal analysis (CMA) by considering the asymptotic limit of the sound pressure level as the number of acoustic and/or structural modes approaches infinity. Using AMA, results identical to those of Statistical Energy Analysis (SEA) were obtained for the spatially-averaged sound pressure levels in the interior. AMA is systematically derived from CMA and therefore the degree of generality of the end result can be adjusted through the choice of appropriate simplifying assumptions. For example, AMA can be used to obtain local sound pressure levels at particular points inside the enclosure, or to include the effects of varying the size and/or location of the sound source. AMA theoretical results were compared with CMA theory and also with experiment for the case where the structural-acoustic enclosure is a rectangular cavity with part of one wall flexible and vibrating, while the rest of the cavity is rigid.

Dowell, Earl H.↗

How Do You Go From a Concept Idea to a NASA Selected Mission? Formulating the Psyche Discovery Mission with JPL's Concurrent Engineering Teams

JPL’s Office of Formulation provides continuity of support and access to domain subject matter experts, as Principal Investigators mature their mission concepts from “cocktail napkin” ideas to Preliminary Design Reviews [1]. Using NASA’s Psyche mission as a case study, we describe JPL’s concurrent engineering A-Team and Team X support to the Psyche competed concept study team in the areas of 1) Initial Feasibility, 2) Trade Space Exploration, 3) Spacecraft Point Design and Cost Estimate, 4) Science, Technical, Management, and Cost Review, and 5) Strategy and Communication Development. NASA’s Psyche Discovery-class mission started as a grassroots idea from Principal Investigator L.T. ElkinsTanton. Is there a compelling Discovery mission to visit the interior of a body for the first time, by sending a mission to an iron metal asteroid? In less than five years the Psyche concept was selected as a mission under NASA’s Discovery Program. While Psyche had a dedicated concept development team [2], they utilized JPL’s concurrent engineering teams, methods, analysis tools, and subject matter experts throughout their mission concept formulation lifecycle

Ziemer, John↗

Application of two-point implicit central-difference methods to hyperbolic systems

This paper presents a general solution algorithm for the set of difference equations that arise when two-point central differences are used to approximate the flux difference terms in systems of hyperbolic differential equations. The general algorithm eliminates the weak points associated with the nonstandard algorithm reported by Wornom and Hafez (1986). The disadvantages of their algorithm relate to its implementation. It consists of separate algorithms for subsonic, supersonic, sonic and shock cells, applied individually, which presents a major bookkeeping problem when multiple sonic and shock cells are present. The general algorithm eliminates this problem and introduces an improved shock treatment which produces shocks with at most one interior shock point.

Wornom, Stephen F.↗

Space station interior noise analysis program

Documentation is provided for a microcomputer program which was developed to evaluate the effect of the vibroacoustic environment on speech communication inside a space station. The program, entitled Space Station Interior Noise Analysis Program (SSINAP), combines a Statistical Energy Analysis (SEA) prediction of sound and vibration levels within the space station with a speech intelligibility model based on the Modulation Transfer Function and the Speech Transmission Index (MTF/STI). The SEA model provides an effective analysis tool for predicting the acoustic environment based on proposed space station design. The MTF/STI model provides a method for evaluating speech communication in the relatively reverberant and potentially noisy environments that are likely to occur in space stations. The combinations of these two models provides a powerful analysis tool for optimizing the acoustic design of space stations from the point of view of speech communications. The mathematical algorithms used in SSINAP are presented to implement the SEA and MTF/STI models. An appendix provides an explanation of the operation of the program along with details of the program structure and code.

Stusnick, E.↗

In-situ synchrotron X-ray micro-diffraction investigation of ultra-low-strain deformation microstructure in laminated Ti-Al composites

An ultra-low-strain deformation microstructure was revealed for the first time non-destructively in the bulk interior of an annealed laminated Ti-Al composite-in the "fully recrystallized" Al layer-by a synchrotron-based micro-diffraction technique, namely differential aperture X-ray microscopy (DAXM), through real space mapping with a very high angular resolution (0.01 degrees). This ultra-low-strain deformation microstructure was found to result from the thermal stress, induced during cooling after annealing, due to the different coefficients of thermal expansion for the Ti and Al layers. The annealed sample was further tensile deformed to a strain of 1.66% and followed by in situ DAXM and analyzed by various misorientation mapping methods. The results pointed to the important effects of the initial microstructure and the interface constraint, as well as the grain size and crystal orientation, on the plastic deformation. A gradient in dislocation density from the layer interface to the center was found in the Al layer of the annealed sample, and this gradient increased slightly during tensile deformation. The variation of the dislocation density was further discussed based on the activation and interaction of dislocations in grains of different sizes and orientations during plastic deformation. The findings of this study provided valuable insights in understanding the constraint effect of the laminated metal composite and the design of novel composite materials. (c) 2020 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Crystallographic orientation↗

Vibrations and structureborne noise in space station

Analytical models were developed to predict vibrations and structureborne noise generation of cylindrical and rectangular acoustic enclosures. These models are then used to determine structural vibration levels and interior noise to random point input forces. The guidelines developed could provide preliminary information on acoustical and vibrational environments in space station habitability modules under orbital operations. The structural models include single wall monocoque shell, double wall shell, stiffened orthotropic shell, descretely stiffened flat panels, and a coupled system composed of a cantilever beam structure and a stiffened sidewall. Aluminum and fiber reinforced composite materials are considered for single and double wall shells. The end caps of the cylindrical enclosures are modeled either as single or double wall circular plates. Sound generation in the interior space is calculated by coupling the structural vibrations to the acoustic field in the enclosure. Modal methods and transfer matrix techniques are used to obtain structural vibrations. Parametric studies are performed to determine the sensitivity of interior noise environment to changes in input, geometric and structural conditions.

Vaicaitis, R.↗

Impact-Seismic Investigations of the InSight Mission

Impact investigations will be an important aspect of the InSight mission. One of the scientific goals of the mission is a measurement of the current impact rate at Mars. Impacts will additionally inform the major goal of investigating the interior structure of Mars. In this paper, we review the current state of knowledge about seismic signals from impacts on the Earth, Moon, and laboratory experiments. We describe the generalized physical models that can be used to explain these signals. A discussion of the appropriate source time function for impacts is presented, along with spectral characteristics including the cutoff frequency and its dependence on impact momentum. Estimates of the seismic efficiency (ratio between seismic and impact energies) vary widely. Our preferred value for the seismic efficiency at Mars is 5 × 10−4, which we recommend using until we can measure it during the InSight mission, when seismic moments are not used directly. Effects of the material properties at the impact point and at the seismometer location are considered. We also discuss the processes by which airbursts and acoustic waves emanate from bolides, and the feasibility of detecting such signals. We then consider the case of impacts on Mars. A review is given of the current knowledge of present-day cratering on Mars: the current impact rate, characteristics of those impactors such as velocity and directions, and the morphologies of the craters those impactors create. Several methods of scaling crater size to impact energy are presented. The Martian atmosphere, although thin, will cause fragmentation of impactors, with implications for the resulting seismic signals. We also benchmark several different seismic modeling codes to be used in analysis of impact detections, and those codes are used to explore the seismic amplitude of impactinduced signals as a function of distance from the impact site. We predict a measurement of the current impact flux will be possible within the timeframe of the prime mission (one Mars year) with the detection of ∼ a few to several tens of impacts. However, the error bars on these predictions are large. Specific to the InSight mission, we list discriminators of seismic signals from impacts that will be used to distinguish them from marsquakes. We describe the role of the InSight Impacts Science Theme Group during mission operations, including a plan for possible night-time meteor imaging. The impacts detected by these methods during the InSight mission will be used to improve interior structure models, measure the seismic efficiency, and calculate the size frequency distribution of current impacts.

Ingrid Daubar↗

The helium abundance of Jupiter from Voyager

The helium abundance in the Jovian atmosphere is derived from Voyager 1 data by two methods. The first method uses only infrared spectra from selected locations on the planet while the second method uses a thermal profile independently derived from radio occultation measurements and infrared spectra recorded near the occultation point. A hydrogen mole fraction of 0.897 plus or minus 0.030 is obtained from the first method, while the second method gives 0.880 plus or minus 0.036, corresponding to helium mass fractions of 0.19 plus or minus 0.05 and 0.21 plus or minus 0.06, respectively. The estimated errors for the first method are primarily due to systematic uncertainties in the H2 and He absorption coefficients, while those for the second method result mainly from errors in the radio occultation profile and are less well known. Random errors in the measured infrared spectra are found to be negligible in both cases. The results are consistent with a uniform mix of hydrogen and helium within Jupiter's interior, but a modest amount of helium depletion (Delta Y equal to or less than 0.05) cannot be excluded.

Gautier, D.↗

Solar Seismology from Space. a Conference at Snowmass, Colorado

The quality of the ground based observing environment suffers from several degrading factors: diurnal interruptions and thermal variations, atmospheric seeing and transparency fluctuations and adverse weather interruptions are among the chief difficulties. The limited fraction of the solar surface observable from only one vantage point is also a potential limitation to the quality of the data available without going to space. Primary conference goals were to discuss in depth the scientific return from current observations and analyses of solar oscillations, to discuss the instrumental and site requirements for realizing the full potential of the seismic analysis method, and to help bring new workers into the field by collecting and summarizing the key background theory. At the conclusion of the conference there was a clear consensus that ground based observation would not be able to provide data of the quality required to permit a substantial analysis of the solar convection zone dynamics or to permit a full deduction of the solar interior structure.

Ulrich, R. K.↗

Matrix operator theory of radiative transfer. I - Rayleigh scattering.

An entirely rigorous method for the solution of the equations for radiative transfer based on the matrix operator theory is reviewed. The advantages of the present method are: (1) all orders of the reflection and transmission matrices are calculated at once; (2) layers of any thickness may be combined, so that a realistic model of the atmosphere can be developed from any arbitrary number of layers, each with different properties and thicknesses; (3) calculations can readily be made for large optical depths and with highly anisotropic phase functions; (4) results are obtained for any desired value of the surface albedo including the value unity and for a large number of polar and azimuthal angles; (5) all fundamental equations can be interpreted immediately in terms of the physical interactions appropriate to the problem; and (6) both upward and downward radiance can be calculated at interior points from relatively simple expressions.

Plass, G. N.↗

Analysis of liquid-metal-jet impingement cooling in a corner region and for a row of jets

A conformal mapping method was used to analyze liquid-metal-jet impingement heat transfer. The jet flow region and energy equation are transformed to correspond to uniform flow in a parallel plate channel with nonuniform heat addition along a portion of one wall. The exact solution for the wall-temperature distribution was obtained in the transformed channel, and the results are mapped back into the physical plane. Two geometries are analyzed. One is for a single slot jet directed either into an interior corner formed by two flat plates, or over the external sides of the corner; the flat plates are uniformly heated, and the corner can have various included angles. The heat-transfer coefficient at the stagnation point at the apex of the plates is obtained as a function of the corner angle, and temperature distributions are calculated along the heated walls. The second geometry is an infinite row of uniformly spaced parallel slot jets impinging normally against a uniformly heated plate. The heat-transfer behavior is obtained as a function of the spacing between the jets. Results are given for several jet Peclet numbers from 5 to 50.

Siegel, R.↗

Transformation of two and three-dimensional regions by elliptic systems

The research during this period continued to expand the class of numerical algorithms that can be accurately and efficiently implemented on overlapping grids. Whereas previous calculations have been used to solve elliptic equations and to find the steady-state solution of parabolic equations, the present work is aimed towards developing time-accurate solution techniques for parabolic and hyperbolic equations. The primary difficulty here is in the correct treatment of the interior boundary nodes that must be updated at each iteration. The implementation of explicit methods is straightforward. However, the common practice of lagging these values when using an implicit methods leads to inconsistencies in the difference equation. One way to avoid this problem is to alternately calculate with an implicit and an explicit method on each subgrid. With this procedure, the explicit method generates boundary values at the next time level which are then used by the implicit step. It can be shown that when a backward implicit method is combined with a forward explicit method, the composite method is second order accurate and unconditionally stable for linear problems. A second area in which progress can be reported is in the distribution of grid points on curves and surfaces.

Mastin, E. W.↗

Optical Polarization of Light from a Sorghum Canopy Measured Under Both a Clear and an Overcast Sky

Introduction: We tested the hypothesis that the optical polarization of the light reflected by a sorghum canopy is due to a Fresnel-type redirection, by sorghum leaf surfaces, of light from an unpolarized light source, the sun or overcast sky, toward the measuring sensor. If it can be shown that the source of the polarization of the light scattered by the sorghum canopy is a first surface, Fresnel-type reflection, then removing this surface reflected light from measurements of canopy reflectance presumably would allow better insight into the biochemical processes such as photosynthesis and metabolism that occur in the interiors of sorghum canopy leaves. Methods: We constructed a tower 5.9m tall in the center of a homogenous sorghum field. We equipped two Barnes MMR radiometers with polarization analyzers on the number 1, 3 and 7 Landsat TM wavelength bands. Positioning the radiometers atop the tower, we collected radiance data in 44 view directions on two days, one day with an overcast sky and the other, clear and sunlit. From the radiance data we calculated the linear polarization of the reflected light for each radiometer wavelength channel and view direction. Results and Discussion: Our experimental results support our hypothesis, showing that the amplitude of the linearly polarized portion of the light reflected by the sorghum canopy varied dramatically with view azimuth direction under a point source, the sun, but the amplitude varied little with view azimuth direction under the hemispherical source, the overcast sky. Under the clear sky, the angle of polarization depended upon the angle of incidence of the sunlight on the leaf, while under the overcast sky the angle of polarization depended upon the zenith view angle. These results support a polarized radiation transport model of the canopy that is based upon a first surface, Fresnel reflection from leaves in the sorghum canopy.

cloud cover↗

Diffraction by a symmetric material junction simulated with generalized sheet transistion conditions

When diffraction by a material discontinuity in a thick dielectric/ferrite slab is treated for the dual integral equation approach (a variation of the Weiner-Hopf method), the resulting solution displays unknown constants, characteristic of the incompleteness associated with the generalized transition conditions employed in modeling a thick slab. The constants can be shown to depend on the geometry and properties of the discontinuity. Their explicit determination by enforcing field continuity across the slab junction is treated. The field internal to the slab is determined by first expressing the field as an expansion of discrete and continuous eigenmodes valid in the exterior and interior slab region. The expansion or mode coefficients are then identified by recasting the Weiner-Hopf solution in a suitable form, enabling the solution constants to be numerically computed by matching the field expansions to the left and right of the slab discontinuity at a discrete set of points.

Ricoy, Mark A.↗

Applying Dynamical Systems Theory to Optimize Libration Point Orbit Stationkeeping Maneuvers for WIND

NASA's WIND mission has been operating in a large amplitude Lissajous orbit in the vicinity of the interior libration point of the Sun-Earth/Moon system since 2004. Regular stationkeeping maneuvers are required to maintain the orbit due to the instability around the collinear libration points. Historically these stationkeeping maneuvers have been performed by applying an incremental change in velocity, or (delta)v along the spacecraft-Sun vector as projected into the ecliptic plane. Previous studies have shown that the magnitude of libration point stationkeeping maneuvers can be minimized by applying the (delta)v in the direction of the local stable manifold found using dynamical systems theory. This paper presents the analysis of this new maneuver strategy which shows that the magnitude of stationkeeping maneuvers can be decreased by 5 to 25 percent, depending on the location in the orbit where the maneuver is performed. The implementation of the optimized maneuver method into operations is discussed and results are presented for the first two optimized stationkeeping maneuvers executed by WIND.

WIND↗

A group finder algorithm optimised for the study of local galaxy environments

Context. The majority of galaxy group catalogues available in the literature use the popular friends-of-friends algorithm which links galaxies using a linking length. One potential drawback to this approach is that clusters of points can be linked with thin bridges which may not be desirable. In order to study galaxy groups, it is important to obtain realistic group structures. Aim. Here, in this study, we present a new simple group finder algorithm, TD-ENCLOSER, that finds the group that encloses a target galaxy of interest. Methods. TD-ENCLOSER is based on the kernel density estimation method which treats each galaxy, represented by a zero-dimensional particle, as a two-dimensional circular Gaussian. The algorithm assigns galaxies to peaks in the density field in order of density in descending order (‘top down’) so that galaxy groups ‘grow’ around the density peaks. Outliers in under-dense regions are prevented from joining groups by a specified hard threshold, while outliers at the group edges are clipped below a soft (blurred) interior density level. Results. The group assignments are largely insensitive to all free parameter variations apart from the hard density threshold and the kernel standard deviation, although this is a known feature of density-based group finder algorithms and it operates with a computing speed that increases linearly with the size of the input sample. In preparation for a companion paper, we also present a simple algorithm to select unique representative groups when duplicates occur. Conclusions. TD-ENCLOSER is tested on a mock galaxy catalogue using a smoothing scale of 0.3 Mpc and is found to be able to recover the input group distribution with sufficient accuracy to be applied to observed galaxy distributions.

79 ASTRONOMY AND ASTROPHYSICS↗