SEARCH · Engineering Papers
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.
Methods and lessons learned in developing GUI-enabled radiative modeling suite for particle-based simulations including accelerated MCRT, distribution factor surrogate-models, and volumetric participating media
Highly accelerated boundary layer turbulence structure, measuring intermittency factor distribution with hot-wire probe
Laboratory and field tests with an experimental jack and an NX-borehole jack are reported. The following conclusions were made: Under borehole jack loading, a circular opening in a brittle solid fails by tensile fracturing when the bearing plate width is not too small. Two proposed contact stress distributions can explain the mechanism of tensile fracturing. The contact stress distribution factor is a material property which can be determined experimentally. The borehole tensile strength is larger than the rupture flexural strength. Knowing the magnitude and orientation of the in situ stress field, borehole jack test results can be used to determine the borehole tensile strength. Knowing the orientation of the in situ stress field and the flexural strength of the rock substance, the magnitude of the in situ stress components can be calculated. The detection of very small cracks is essential for the accurate determination of the failure loads which are used in the calculation of strengths and stress components.
The design of a real-time rem-rad dosimeter with sufficient generality for inclusion of dose distribution factors for space applications is discussed. This generalized dosimetric system is only slightly more complex than dosimeters in current use.
From 1974 to 1977 the application of remote sensing methods in coastal areas and tidal bays and estuaries was investigated on the German coast of the North Sea. Aerial photographs were taken using different films; (1) color, (2) color infrared, and (3) black and white films. Scanner recordings were taken by an 11 channel scanner. Ground truth measurements of radiation and measurements of meteorological elements were carried out. For mapping the morphology in mudflat areas a digital texture analysis was developed, by which measurement of the change of image structures cased by distributing factors, such as changing illumination, is possible.
The ride quality experienced by passengers is a function of airframe rigid-body, elastic dynamic responses, autopilot, and stability augmentation system control inputs. A frequency response method has been developed to select sinusoidal elevator input time histories yielding vertical load factor distributions, within a given limit, as a function of fuselage station. The numerical technique is illustrated by applying two-degree-of-freedom short-period and first symmetric mode equations of motion to a B-1 aircraft at Mach 0.85 during sea level flight conditions.
A simple dosimeter design has been developed by NASA to monitor the space proton dose (in rads) to a distributed body organ as a linear combination of ion chambers with varying wall thickness. Estimated doses are given for ion chambers of thickness 2, 3, 4, and 5 g/sq cm. The analytical equation used to calculate the dose distribution factor is also given.
Relativistic beaming of radiation is generally invoked to explain the extreme properties of blazars. The simplest synchrotron-self-Compton (SSC) model is applied to all blazars that have been observed in the X-ray region and have known VLBI sizes, and the minimum kinematic Doppler factor for each source is derived. All high polarization quasars (HPQs), but only half of the 30 BL Lac objects in the sample require relativistic enhancement of radiation due to beaming effects. The different Doppler factor distributions for the two classes of sources argue against the hypothesis that BL Lacs could be low-luminosity HPQs whose strong emission lines are swamped out by the boosted nonthermal continuum. Thus HPQs and BL Lacs could be intrinsically different in their physical structure.
A computer simulation for 'running wet' and evaporative aircraft anti-icing systems is developed. The model is based on the analysis of the liquid water film which forms in the regions of direct impingement and, then, breaks up near the impingement limits into rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted using a stability analysis theory and the laws of mass energy conservation. The solid structure is modeled as a multiple layer wall. The anti-icing system modeled is of the thermal type utilizing hot air and/or electrical heating elements embedded within the wall layers. Experimental observations revealing some of the basic physics of the water flow on the surface are presented. Detailed qualitative documentation of the tests are given. Several numerical examples are considered, and the effect of some of the involved parameters on the system performance are investigated.
A numerical simulation for 'running wet' aircraft anti-icing systems is developed. The model includes breakup of the water film, which exists in regions of direct impingement, into individual rivulets. The wetness factor distribution resulting from the film breakup and the rivulet configuration on the surface are predicted in the numerical solution procedure. The solid wall is modeled as a multilayer structure and the anti-icing system used is of the thermal type utilizing hot air and/or electrical heating elements embedded with the layers. Details of the calculation procedure and the methods used are presented.
Bone densitometry has previously been used to obtain cross-sectional properties of bone in a single scan plane. Using three non-coplanar scans, we have extended the method to obtain the principal area Moments of inertia and orientations of the principal axes at each cross-section along the length of the scan. Various 5 aluminum phantoms were used to examine scanner characteristics to develop the highest accuracy possible for in vitro non-invasive analysis of mass distribution. Factors considered included X-ray photon energy, initial scan orientation, the included angle of the 3 scans, and Imin/Imax ratios. Principal moments of inertia were accurate to within 3.1% and principal angles were within 1 deg. of the expected value for phantoms scanned with included angles of 60 deg. and 90 deg. at the higher X-ray photon energy. Low standard deviations in error also 10 indicate high precision of calculated measurements with these included angles. Accuracy and precision decreased slightly when the included angle was reduced to 30 deg. The method was then successfully applied to a pair of excised cadaveric tibiae. The accuracy and insensitivity of the algorithms to cross-sectional shape and changing isotropy (Imin/Imax) values when various included angles are used make this technique viable for future in vivo studies.
The adjoint parabold stability equations (PSE) formulation is used to calculate the boundary layer receptivity to localized surface roughness and suction for compressible boundary layers. Receptivity efficiency functions predicted by the adjoint PSE approach agree well with results based on other nonparallel methods including linearized Navier-Stokes equations for both Tollmien-Schlichting waves and crossflow instability in swept wing boundary layers. The receptivity efficiency function can be regarded as the Green's function to the disturbance amplitude evolution in a nonparallel (growing) boundary layer. Given the Fourier transformed geometry factor distribution along the chordwise direction, the linear disturbance amplitude evolution for a finite size, distributed nonuniformity can be computed by evaluating the integral effects of both disturbance generation and linear amplification. The synergistic approach via the linear adjoint PSE for receptivity and nonlinear PSE for disturbance evolution downstream of the leading edge forms the basis for an integrated transition prediction tool. Eventually, such physics-based, high fidelity prediction methods could simulate the transition process from the disturbance generation through the nonlinear breakdown in a holistic manner.
Using a 3-D relativistic particle-in-cell (RPIC) code, we have investigated particle acceleration associated with a relativistic electron-positron (cold) jet propagating into ambient electron-positron and electron-ion plasmas without initial magnetic fields in order to investigate the nonlinear stage of the Weibel instability. We have also performed simulations with broad Lorentz factor distribution of jet electrons and positrons, which are assumed to be created by the photon annihilation. The growth time and nonlinear saturation levels depend on the initial jet parallel velocity distributions and ambient plasma. Simulations show that the Weibel instability created in the collisionless shocks accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The nonlinear fluctuation amplitude of densities, currents, electric, and magnetic fields in the electron-ion ambient plasma are larger than those in the electron-positron ambient plasma. We have shown that plasma instabilities driven by these streaming electron-positron pairs are responsible for the excitation of near-equipartition, turbulent magnetic fields. These fields maintain a strong saturated level on timescales much longer than the electron skin depth at least for the duration of the simulations. Our results reveal the importance of the electromagnetic filamentation instability in ensuring an effective coupling between electron-positron pairs and ions, and may help explain the origin of large upstream fields in GRB shock.
Hypersonic flows over circular cones constitute one of the most important generic configurations for fundamental aerodynamic and aerothermodynamic studies. In this paper, numerical computations are carried out for Mach 6 flows over a 7-degree half-angle cone with two different flow incidence angles and a compression cone with a large concave curvature. Instability wave and transition-related flow physics are investigated using a series of advanced stability methods ranging from conventional linear stability theory (LST) and a higher-fidelity linear and nonlinear parabolized stability equations (PSE), to the 2D eigenvalue analysis based on partial differential equations. Computed N factor distribution pertinent to various instability mechanisms over the cone surface provides initial assessments of possible transition fronts and a guide to corresponding disturbance characteristics such as frequency and azimuthal wave numbers. It is also shown that strong secondary instability that eventually leads to transition to turbulence can be simulated very efficiently using a combination of advanced stability methods described above.
We propose a study of the Dense Warm Interstellar Medium (D-WIM) via the two fine structure lines of ionized nitrogen (\NII). This new phase of the interstellar medium (ISM) has density 10 $\cmv \le n(e) \le$ 100 $\cmv$, and has effectively been identified throughout much of the Milky Way and in nearby galaxies. :The unexpected conditiosn in the D-WIM immediately raise important questions both about this phase of the ISM and its relationship to the ISM as a while. These include: what is the structure of the D-WIM -- e.g. filling factor, distribution, and morphology?; what is the origin of the D-WIM?; what supports its ionization and what is its evolution given huge overpressure relative to the $\ge$ 100 times more diffuse Warm Ionized Medium (WIM)?; what is its possible relationship to stars - is it the the remnants of very old HII regions that were produced by clusters of masive stars? The PRIMA survey requires a coverage of a large fraction of the plane, plus selected regions at high galactic latitudes. A sensitivity of 10$^{-9}$ $Wm^{-2}sr^{-1}$ would be highly desirable. A spectral resolution $R$ of 10$^4$ will isolate individual spectral features. $R$ = 10$^3$ would be acceptable to measure the total intensity. This PRIMA survey could be part of a broader survey of the ISM, and include a number of other important ISM tracers such as [CII] 158 \um, [OI] 63 \um\ and [OI] 146 \um, depending on wavelength coverage of the instrument. The sensitivity requirement for study of the D-WIM is likely more stringent than for the other lines, and may require special consideration for successful implementation as well as for determining the area that can be observed.
The motivation for the development of this tool and the underlying algorithms and methods was to enable the development of high-temporal resolution, realistic time-series data for quasi-static time-series (QSTS) analysis of distribution systems. Often, aggregated load profile data for a distribution circuit is available (e.g. feeder loading data collected via SCADA at the utility substation) and, while this data is typically accurate it masks the considerable variability of the 100’s or 1000’s of individual loads connected on the circuit. This tool was developed to model both the increased variability expected for these individual loads (e.g. the load of a single distribution transformer connected to 8-12 houses) and the expected diversity between loads on the circuit. It is important to note that the difference in variability and diversity, in the context of this tool, is that variability modeling only adds representative variability due to disaggregated load characteristics (e.g. the presence in the load profile of loads turning off and on like an air conditioner/oven) while the average energy profile remains the same as the user supplied power profile. Diversity modeling generates multiple individual load profiles which, in aggregate, sum to the user supplied power profile. Diversity is effectively variability in the energy usage over longer periods of time than seen in the variability model. Put another way, variability modeling supplies the expected variability due to the operation of various end-use loads and diversity modeling supplies the usage differences due to human behavior, schedules, etc. This load modeling tool was developed for use in generating data for distribution systems. Modeling is summarized by two major functions: 1) taking low resolution load profiles and adding intra-seconds variability onto the profiles, and 2) taking a user supplied load profile and distribution factors and adding both diversity and variability to the user supplied profile.
While wind energy production loss due to unavailability, environmental impacts, curtailment, and other causes has been studied and characterized at the utility-scale wind farm level, observation-based characterization of project loss is lacking for distributed wind energy, particularly for projects involving small wind turbines. Contemporary tools and research that support pre-construction distributed wind energy characterization present a wide range of default loss factors to convert gross energy estimates to net: 7-18%. We hypothesize that we can use generation observations from operational distributed wind projects to develop more accurate representations of loss. Using a density-based filtering technique on distributed wind power generation timeseries, we determine periods of typical performance and use them with regression algorithms in a measure-correlate-predict fashion to simulate what the generation would have been during periods of atypical or unreported performance. From there, the actual versus predicted generation leads to the establishment of observation-informed loss factors (median = 17%) for small, single turbine installation distributed wind projects.
This report gives the results of pressure-distribution tests made to study the effects on lateral stability of changing the span-load distribution on a rectangular monoplane wing model of fairly thick section. Three methods of changing the distribution were employed: variation in profile along the span to a thin symmetrical section at the tip, twist from +5 degrees to -15 degrees at the tip, and sweepback from +20 degrees to -20 degrees. The tests were conducted in a 5-foot closed-throat atmospheric wind tunnel. The investigation shows the following results: (1) change in profile along the span from the NACA-84 at the root to the NACA-M2 at the tip considerably reduces lateral instability, but also reduces the general effectiveness of the wing. (2) washout up to 11 degrees progressively reduces maximum lateral instability. (3) transition from sweepforward to sweepback gradually reduces the useful angle-of-attack range, but has no clearly defined effect on maximum lateral instability.