Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape Matching”

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 91 records · Page 5

Calibration improvements to electronically scanned pressure systems and preliminary statistical assessment

Orifice-to-orifice inconsistencies in data acquired with an electronically-scanned pressure system at the beginning of a wind tunnel experiment forced modifications to the standard, instrument calibration procedures. These modifications included a large increase in the number of calibration points which would allow a critical examination of the calibration curve-fit process, and a subsequent post-test reduction of the pressure data. Evaluation of these data has resulted in an improved functional representation of the pressure-voltage signature for electronically-scanned pressures sensors, which can reduce the errors due to calibration curve fit to under 0.10 percent of reading compared to the manufacturer specified 0.10 percent of full scale. Application of the improved calibration function allows a more rational selection of the calibration set-point pressures. These pressures should be adjusted to achieve a voltage output which matches the physical shape of the pressure-voltage signature of the sensor. This process is conducted in lieu of the more traditional approach where a calibration pressure is specified and the resulting sensor voltage is recorded. The fifteen calibrations acquired over the two-week duration of the wind tunnel test were further used to perform a preliminary, statistical assessment of the variation in the calibration process. The results allowed the estimation of the bias uncertainty for a single instrument calibration; and, they form the precursor for more extensive and more controlled studies in the laboratory.

Everhart, Joel L.↗

Stereo-Based Region-Growing using String Matching

We present a novel stereo algorithm based on a coarse texture segmentation preprocessing phase. Matching is performed using a string comparison. Matching sub-strings correspond to matching sequences of textures. Inter-scanline clustering of matching sub-strings yields regions of matching texture. The shape of these regions yield information concerning object's height, width and azimuthal position relative to the camera pair. Hence, rather than the standard dense depth map, the output of this algorithm is a segmentation of objects in the scene. Such a format is useful for the integration of stereo with other sensor modalities on a mobile robotic platform. It is also useful for localization; the height and width of a detected object may be used for landmark recognition, while depth and relative azimuthal location determine pose. The algorithm does not rely on the monotonicity of order of image primitives. Occlusions, exposures, and foreshortening effects are not problematic. The algorithm can deal with certain types of transparencies. It is computationally efficient, and very amenable to parallel implementation. Further, the epipolar constraints may be relaxed to some small but significant degree. A version of the algorithm has been implemented and tested on various types of images. It performs best on random dot stereograms, on images with easily filtered backgrounds (as in synthetic images), and on real scenes with uncontrived backgrounds.

Mandelbaum, Robert↗

Light-Hydrocarbon Bearing Solids on Planetesimal 5145 Pholus

Object 5145 Pholus (=1992 AD) is a planetesimal in an orbit that crosses those of Saturn, Uranus, and Neptune (period 92.7 years). It is particularly notable because of its extreme red color, corroborated by several observing teams. A spectrum of Pholus obtained in 1992 shows a strong absorption band with a characteristic shape at 2.27 micrometers, plus a weaker band at 1.7 micrometers. A better spectrum of the 2.0-2.5 micrometer region in 1993 confirms the position and shape of the 2.27-micrometer band. The color and spectral bands are identified with the aliphatic-rich and high H/C organic solid called asphaltite, which in a terrestrial setting originates from thermal processing of products of biological activity. In Pholus, this material is attributed to formation from radiation processing of ices on grains in the interstellar medium. Laboratory spectra of asphaltite and related materials have been published by Moroz et al., while Cloutis showed similar bands in comparable materials and identified them as the overtone and combination bands of C-H stretching and bending modes in CH2 and CH3 groups. Asphaltites, kerites, and anthraxolites are solid non-graphite members of a sequence ranging from oil to graphite; diffuse reflectance spectra of suites of these intermediate materials show color characteristics similar to those of the low-albedo asteroids (C,P,D), although specific identifications have not been made because of the lack of distinct absorption bands in the spectra of most low-albedo solar system bodies. In the case of Pholus, however, the primary band is strong; its wavelength and its shape, plus the match of the extremely red color, leads us to the identification of aliphatic-rich, asphaltite-like organic solid. The C, P, and D-type asteroids vary in degree of "redness", but are all less red than Pholus. Pholus and the Ctype asteroids are the end members of a sequence that represents the radiation processing of hydrocarbons, with Pholus being the least processed. Solar irradiation processes and heating reduce the H/C and aliphatic content of hydrocarbons preserved from the interstellar medium, and in the end produce opaque solids of neutral reflectance, including the kerogens (similar to anthraxolites) found in profusion in the carbonaceous meteorites.

Cruikshank, Dale P.↗

Light-Hydrocarbon Bearing Solids on Planetesimal 5145 Pholus

Object 5145 Pholus (=1992 AD) is a planetesimal in an orbit that crosses those of Saturn, Uranus, and Neptune (period 92.7 years). It is particularly notable because of its extreme red color, corroborated by several observing teams. A spectrum of Pholus obtained in 1992 shows a strong absorption band with a characteristic shape at 2.27 micron, plus a weaker band at 1.7 microns. A better spectrum of the 2.0-2.5 micron region in 1993 confirms the position and shape of the 2.27 micron band. The color and spectral bands are identified with the aliphatic-rich and high H/C organic solid called asphaltite, which in a terrestrial setting originates from thermal processing of products of biological activity. In Pholus, this material is attributed to formation from radiation processing of ices on grains in the interstellar medium. Laboratory spectra of asphaltite and related materials have been published by Moroz et al., while Cloutis showed similar bands in comparable materials and identified them as the overtone and combination bands of C-H stretching and bending modes in CH2 and CH3 groups. Asphaltites, kerites, and anthraxolites are solid non-graphite members of a sequence ranging from oil to graphite; diffuse reflectance spectra of suites of these intermediate materials show color characteristics similar to those of the low-albedo asteroids (C,P,D), although specific identifications have not been made because of the lack of distinct absorption bands in the spectra of most low-albedo solar system bodies. In the case of Pholus, however, the primary band is strong; its wavelength and its shape, plus the match of the extremely red color, leads us to the identification of aliphatic-rich, asphaltite-like organic solid. The C, P, and D-type asteroids vary in degree of 'redness', but are all less red than Pholus. Pholus and the C-type asteroids are the end members of a sequence that represents the radiation processing of hydrocarbons, with Pholus being the least processed. Solar irradiation processes and heating reduce the H/C and aliphatic content of hydrocarbons preserved from the interstellar medium, and in the end produce opaque solids of neutral reflectance, including the kerogens (similar to anthraxolites) found in profusion in the carbonaceous meteorites.

Cruikshank, Dale P.↗

Investigation of Atmospheric Boundary-Layer Effects on Launch-Vehicle Ground Wind Loads

A launch vehicle ground-wind-loads program is underway at the NASA Langley Transonic Dynamics Tunnel. The objectives are to quantify key aerodynamic and structural characteristics that impact the occurrence of large wind-induced oscillations of a launch vehicle when exposed to ground winds prior to launch. Of particular interest is the dynamic response of a launch vehicle when a von Kármán vortex street forms in the wake of the vehicle resulting in quasiperiodic lift and drag forces. Vehicle response to these quasiperiodic forces can become quite large when the frequency of vortex shedding nears that of a lowly-damped structural mode thereby exciting a resonant response. Wind approaching the vehicle can be characterized by a varying speed with height and turbulence content. The combination of both the varying speed and turbulence content is referred to herein as the atmospheric boundary-layer. The importance of the atmospheric boundary-layer upon launch vehicle wind-induced oscillation response has long been questioned, and its effects are not well understood. Although there are several facilities around the world dedicated to replicating atmospheric boundary layers, the development of such a boundary layer in a wind tunnel capable of producing flight-representative Reynolds numbers for aeroelastically-scaled launch vehicle models has only recently been accomplished. The NASA Langley Transonic Dynamics Tunnel is capable of simulating flight-representative Reynolds numbers of launch vehicles on the pad and is uniquely capable of replicating many fluid-structure scaling parameters typical of aeroelastic tests. Recent test efforts successfully developed representative atmospheric boundary-layers for three launch sites in the Transonic Dynamics Tunnel, thereby allowing all known aerodynamic and fluid-structure coupling parameters to be simultaneously simulated for those sites. Dynamic aeroelastically-scaled models representative of typical large launch vehicles were constructed for testing. Aeroelastic scaling includes matching geometry, mode shapes, reduced frequencies, damping, running mass ratios, and running stiffness ratios. The models were tested in smooth uniform flow and then immersed in the atmospheric boundary-layer for comparison of these effects. Dynamic data were acquired measuring unsteady pressure, acceleration, and base bending moment. It was discovered that peak dynamic loads resulting from resonant wind-induced oscillation response are similar when acquired in either smooth uniform flow or an atmospheric boundary-layer. This indicates that resonant lock-in events are minimally impacted by representative turbulence and/or wind profile. Alternately, nonresonant wind-induced oscillation response events are stronger when acquired in an atmospheric boundary-layer. This indicates that a lowly-damped structural response will increase when exposed to an increased magnitude of random excitation, which is consistent with historical comparisons. Loads created by the resonant response events were substantially stronger than those from the nonresonant response events. Therefore, if testing is done to simply identify worst-case conditions and load magnitude, then smooth uniform flow is likely an adequate test technique. However, if nonresonant response loads are of primary interest, then atmospheric boundary-layer simulation is required.

Thomas G Ivanco↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Integrated freeform optical couplers

Reflecting light beams off of microscale three-dimensional (3D) freeform surfaces can yield highly efficient coupling into and out of optical waveguides, optical fibers, and photonic chips. The structure of the 3D freeform reflective surface determines the shape of the reflected beam. This allows freeform reflectors to control the mode profile, rotation angle, and divergence angle of light beams. Control of beam shape enables mode matching between source output mode and target input mode, which results in low-loss optical coupling. An inventive freeform reflective surface can direct light beams in plane or out of plane via specular reflection or total internal reflection. A photonic integrated circuit with this type of freeform optical coupler can operate with a bandwidth range of at least 400 nm, potentially encompassing all visible or telecommunications wavelengths, and can be volume manufactured in photonic chips.

Yu, Shaoliang↗

A Data-Driven Pivot-Point-Based Time-Series Feeder Load Disaggregation Method

The load profile at a feeder-head is usually known to utility engineers while the nodal load profiles are not. However, the nodal load profiles are increasingly important for conducting time-series analysis in distribution systems. Therefore, in this paper, we present a pivot-point based, two-stage feeder load disaggregation algorithm using smart meter data. The two stages are load profile selection (LPS) and load profile allocation (LPA). In the LPS stage, a random load profile selection process is first executed to meet the load diversity requirement. Then, a few pairs of pivot points are selected as the matching targets. After that, a matching algorithm will run repetitively to select one load profile at a time for matching the reference load profile at the pivot points. In the LPA stage, the LPS selected load profiles are allocated to each load node on the feeder considering distribution transformer loading limits, load composition, and square-footage. The proposed method is validated using actual data collected in a North Carolina service area. Finally, simulation results show that the proposed method can generate a unique load shape for each load node while match the shape of their aggregated profile with the actual feeder head load profile.

42 ENGINEERING↗

Imaging the effect of dipole tilt on magetotail boundaries

Eight years of Interplanetary Monitoring Platform (IMP) 8, four years of International Sun Earth Explorer (ISEE) 2, and one year of IMP 7 magnetometer data have been combined to produce an 'image' of the average magnetic field for a YZ cross section (aberrated GSM) of the magnetotail at a downtail distance of 25 R(sub E). The shape of the neutral sheet and magnetopause boundaries can be observed directly from the images. A fitting function that qualitatively matches the observed boundary shape can then be chosen. This approach improves on previous fits to possible unsuitable functional forms specified independently of the data. In addition, as a refinement of previous studies, we have corrected for varying solar wind dynamic pressure and the effects of tail flaring. We find the magnetopause is displaced above the XY plane with increasing dipole tilt. The neutral sheet is found to curve slightly more than the model of Fairlfield (1980) during times of large dipole tilt and near the flanks appears to differ substantially from the neutral sheet shape given by the analytic model of Voigt (1984), the more recent neutral sheet model of Dandouras (1988) based on the Voigt model, and the semi-empirical model of Tsyganenko (1989).

Hammond, C. Max↗

Build-Up Approach to Updating the Mock Quiet Spike Beam Model

When a new aircraft is designed or a modification is done to an existing aircraft, the aeroelastic properties of the aircraft should be examined to ensure the aircraft is flight worthy. Evaluating the aeroelastic properties of a new or modified aircraft can include performing a variety of analyses, such as modal and flutter analyses. In order to produce accurate results from these analyses, it is imperative to work with finite element models (FEM) that have been validated by or correlated to ground vibration test (GVT) data, Updating an analytical model using measured data is a challenge in the area of structural dynamics. The analytical model update process encompasses a series of optimizations that match analytical frequencies and mode shapes to the measured modal characteristics of structure. In the past, the method used to update a model to test data was "trial and error." This is an inefficient method - running a modal analysis, comparing the analytical results to the GVT data, manually modifying one or more structural parameters (mass, CG, inertia, area, etc.), rerunning the analysis, and comparing the new analytical modal characteristics to the GVT modal data. If the match is close enough (close enough defined by analyst's updating requirements), then the updating process is completed. If the match does not meet updating-requirements, then the parameters are changed again and the process is repeated. Clearly, this manual optimization process is highly inefficient for large FEM's and/or a large number of structural parameters. NASA Dryden Flight Research Center (DFRC) has developed, in-house, a Mode Matching Code that automates the above-mentioned optimization process, DFRC's in-house Mode Matching Code reads mode shapes and frequencies acquired from GVT to create the target model. It also reads the current analytical model, as we11 as the design variables and their upper and lower limits. It performs a modal analysis on this model and modifies it to create an updated model that has similar mode shapes and frequencies as those of the target model. The Mode Matching Code output frequencies and modal assurance criteria (MAC) values that allow for the quantified comparison of the updated model versus the target model. A recent application of this code is the F453 supersonic flight testing platform, NASA DFRC possesses a modified F-15B that is used as a test bed aircraft for supersonic flight experiments. Traditionally, the finite element model of the test article is generated. A GVT is done on the test article ta validate and update its FEM. This FEM is then mated to the F-15B model, which was correlated to GVT data in fall of 2004, A GVT is conducted with the test article mated to the aircraft, and this mated F-15B/ test article FEM is correlated to this final GVT.

Herrera, Claudia Y.↗

A Three-Dimensional Coupled Internal/External Simulation of a Film-Cooled Turbine Vane

A three-dimensional Navier-Stokes simulation has been performed for a realistic film-cooled turbine vane using the LeRC-HT code. The simulation includes the flow regions inside the coolant plena and film cooling holes in addition to the external flow. The vane is the subject of an upcoming NASA Glenn Research Center experiment and has both circular cross-section and shaped film cooling holes. This complex geometry is modeled using a multi-block grid which accurately discretizes the actual vane geometry including shaped holes. The simulation matches operating conditions for the planned experiment and assumes periodicity in the spanwise direction on the scale of one pitch of the film cooling hole pattern. Two computations were performed for different isothermal wall temperatures, allowing independent determination of heat transfer coefficients and film effectiveness values. The results indicate separate localized regions of high heat transfer coefficient values, while the shaped holes provide a reduction in heat flux through both parameters. Hole exit data indicate rather simple skewed profiles for the round holes, but complex profiles for the shaped holes with mass fluxes skewed strongly toward their leading edges.

Heidmann, James D.↗

Increased Efficiency LED

In a semiconductor light emitting diode (LED) a large proportion of the light produced will be internally reflected when it strikes the end face of the LED at a semiconductor-air interface. This is due to the high index of refraction of most LED semiconductor materials. This problem may be partially overcome by modification of the shape of the LED by a reverse taper of the sides of the LED. Light is redirected by the taper to strike the interface at an angle closer to normal. This allows light to exit the LED that would be totally internally reflected in an untapered LED. The novelty of the present invention lies in the tapering of the sides of the LED for increased through transmission of light. Prior art devices have made use surface modifications such as roughening, addition of guiding layers, use of index matching materials and hemispherical shaping of the emitting end of the LED. Current technology LEDs have a transmission efficiency of only a few percent. An increase in this efficiency would be of great value to any technology that makes use of LEDs for light generation. The present invention is related to LAR 1 5050-SB, Collection of Light from Optical Fibers of NA greater than 1.

Egalon, Claudio Oliviera↗

Tuning Perovskites’ Hydration-Induced Chemical Expansion with Octahedral Tilt Angles

Hydration-induced strains in proton-conducting oxides compromise chemo-mechanical stability when these materials are applied in protonic ceramic electrochemical cells. Here, to develop design principles for zero-strain materials, we systematically studied the hydration coefficients of chemical expansion (CCEs) in perovskite (Sr, Ba)(Ce, Zr, Y)O 3–x solid solutions with in situ dilatometry and thermogravimetric analysis in the range of 430–630 °C. By including and decoupling a wide range of tolerance factors and lattice parameters, we were able to identify a minimum in hydration CCEs (0–0.02) at intermediate tolerance factor values (t ≈ 0.95). Conversely, despite expectations of lower CCEs in larger unit cells, no general trend in CCE versus lattice parameter was found, and opposite trends could be seen for Sr(Ce, Zr, Y)O 3–x versus Ba(Ce, Zr, Y)O 3–x separately. In situ neutron diffraction (ND) enabled atomistic insight. Upon decreasing t, chemical strain anisotropy increased, but this trend did not match the U-shaped dependence of macroscopic CCEs on t. Instead, perovskites with intermediate t, hosting intermediate octahedral tilt angles in the nominally dry state, underwent the largest change in the B–O–B angles during hydration. Accommodating hydration through decreasing B–O–B angles is beneficial because it does not result in large lattice parameter changes. We propose an intermediate tolerance factor as a simple structural descriptor to enable near-zero hydration strains in proton-conducting perovskites.

36 MATERIALS SCIENCE↗

Refractive axicon for X-ray microscopy applications: design, optimization, and experiment

In a full-field transmission X-ray microscopy (TXM) setup, a condenser X-ray optical element is used to illuminate the sample by condensing the X-ray beam delivered by the synchrotron storage ring. On-going and future upgrades of synchrotron facilities to diffraction-limited storage rings will pose new challenges to these TXM setups, such as much smaller X-ray beams on the condenser. Here, we demonstrate that a refractive axicon can be used as an X-ray beam shaper to match the ring-shaped aperture of the condenser. Aiming at more efficient use of the incoming X-ray intensity, we explore several axicon designs both analytically and with numerical simulations. The axicons were produced by two-photon polymerization 3D printing on thin silicon nitride membrane substrates. The first characterization of the axicon was carried out at the TOMCAT beamline of the Swiss Light Source (Switzerland).

47 OTHER INSTRUMENTATION↗

Dimming the Lights: 2D Simulations of Deflagrations of Hybrid C/O/Ne White Dwarfs Using FLASH

The dimmest and most numerous outlier of the Type Ia supernova population, the Type Iax event, is increasingly being found in the results of observational campaigns. There is currently no single accepted model to describe these events. This 2D study explores the viability of modeling Type Iax events as a hybrid C/O/Ne white dwarf progenitor undergoing a deflagration using the multiphysics software FLASH. This hybrid was created using the stellar evolution code MESA, and its C-depleted core and mixed structure have demonstrated lower yields than traditional C/O progenitors in previous deflagration-to-detonation studies. To generate a sample, 30 “realizations” of this simulation were performed, the only difference being the shape of the initial match head used to start the deflagration. Consistent with earlier work, these realizations produce the familiar hot dense bound remnant surrounded by sparse ejecta. Our results indicate that the majority of the star remains unburned (∼70%) and bound (>90%). Our realizations produce total ejecta yields on the order of 10 −2 –10 −1 M ☉ , ejected 56 Ni yields on the order of 10 −4 –10 −2 M ☉ , and ejecta kinetic energies on the order of 10 48 –10 49 erg. Compared to yields inferred from recent observations of the dimmest Type Iax events—SN 2007qd, SN 2008ha, SN 2010ae, SN 2019gsc, SN 2019muj, SN 2020kyg, and SN 2021fcg—our simulation produces comparable 56 Ni yields but too-small total yields and kinetic energies. Reignition of the remnant is also seen in some realizations.

Hydrodynamical simulations↗

Collisionless tearing in a field-reversed sheet pinch assuming nonparallel propagation

The problem of collisionless linear tearing is examined assuming a wave vector with a component normal to the equilibrium field. The geometry is defined and the general form of the linear dispersion equation is calculated. The linear theory results when k is parallel to B are reviewed, and Ampere's law is calculated for the external adiabatic region when k times B does not equal zero, using two-fluid theory. A solution is obtained for the approximate form of the perturbed currents and vector potential assuming quasi-parallel k. The resonant current contributions within the singular layer are calculated, obtaining an estimate of the dispersion equation. The form of the adiabatic currents within the singular layer is calculated, showing that an x-z current system persists even in the limit k perpendicular to B goes to zero. Finally, the perturbed vector potential solutions across the singular layer are matched to obtain the shape of the complete eigenfunction.

Quest, K. B.↗

The optical structure of the central core in the peculiar supernova remnant CTB 80

Deep and high resolution H-alpha + forbidden N II, O III, and red continuum interference-filter images of the bright central radio core of the supernova remnant CTB 80 are presented. These photographs reveal a surprisingly filamentary structure not visible from previous imaging and indicate considerable morphological differences between forbidden O III and H-alpha + forbidden N II emission regions. The H-alpha + forbidden N II filaments appear elongated in the E-W direction like that of the radio emission plateau surrounding the central core. The forbidden O III emission is much smaller in angular size but matches well the position, shape and size of the core's 6 cm radio emission. Also detected are the two possible optical counterparts reported by Blair and Schild (1984) to the remnants' X-ray point source and it is noted that both may have coincident H-alpha + forbidden N II emission.

Fesen, R. A.↗

Generic hypersonic vehicle performance model

An integrated computational model of a generic hypersonic vehicle was developed for the purpose of determining the vehicle's performance characteristics, which include the lift, drag, thrust, and moment acting on the vehicle at specified altitude, flight condition, and vehicular configuration. The lift, drag, thrust, and moment are developed for the body fixed coordinate system. These forces and moments arise from both aerodynamic and propulsive sources. SCRAMjet engine performance characteristics, such as fuel flow rate, can also be determined. The vehicle is assumed to be a lifting body with a single aerodynamic control surface. The body shape and control surface location are arbitrary and must be defined. The aerodynamics are calculated using either 2-dimensional Newtonian or modified Newtonian theory and approximate high-Mach-number Prandtl-Meyer expansion theory. Skin-friction drag was also accounted for. The skin-friction drag coefficient is a function of the freestream Mach number. The data for the skin-friction drag coefficient values were taken from NASA Technical Memorandum 102610. The modeling of the vehicle's SCRAMjet engine is based on quasi 1-dimensional gas dynamics for the engine diffuser, nozzle, and the combustor with heat addition. The engine has three variable inputs for control: the engine inlet diffuser area ratio, the total temperature rise through the combustor due to combustion of the fuel, and the engine internal expansion nozzle area ratio. The pressure distribution over the vehicle's lower aft body surface, which acts as an external nozzle, is calculated using a combination of quasi 1-dimensional gas dynamic theory and Newtonian or modified Newtonian theory. The exhaust plume shape is determined by matching the pressure inside the plume, calculated from the gas dynamic equations, with the freestream pressure, calculated from Newtonian or Modified Newtonian theory. In this manner, the pressure distribution along the vehicle after body expansion surface is then determined. The aerodynamic modeling, the engine modeling, and the exhaust plume analysis are described in more detail. A description of the computer code used to perform the above calculations is given and an input/output example is then given. The computer code is available on a Macintosh floppy disk.

Chavez, Frank R.↗