Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “1-dimensional systems”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

$2+1$ dimensional Floquet systems and lattice fermions: Exact bulk spectral equivalence

A connection has recently been proposed between periodically driven systems known as Floquet insulators in continuous time and static fermion theories in discrete time. This connection has been established in a (1+1) ( 1 + 1 ) -dimensional free theory, where an explicit mapping between the spectra of a Floquet insulator and a discrete-time Dirac fermion theory has been formulated. Here we investigate the potential of static discrete-time theories to capture Floquet physics in higher dimensions, where so-called anomalous Floquet topological insulators can emerge that feature chiral edge states despite having bulk bands with zero Chern number. Starting from a particular model of an anomalous Floquet system, we provide an example of a static discrete-time theory whose bulk spectrum is an exact analytic match for the Floquet spectrum. The spectra with open boundary conditions in a particular strip geometry also match up to finite-size corrections. However, the models differ in several important respects. The discrete-time theory is spatially anisotropic, so that the spectra do not agree for all lattice terminations, e.g. other strip geometries or on half spaces. This difference can be attributed to the fact that the static discrete-time model is quasi-one-dimensional in nature and therefore has a different bulk-boundary correspondence than the Floquet model.

Iadecola, Thomas (ORCID:0000000251456441)↗
NASA NTRS1994

Design of a variable area diffuser for a 15-inch Mach 6 open-jet tunnel

The Langley 15-inch Mach 6 High Temperature Tunnel was recently converted from a Mach 10 Hypersonic Flow Apparatus. This conversion was effected to improve the capability of testing in Mach 6 air at relatively high reservoir temperatures not previously possible at Langley. Elevated temperatures allow the matching of the Mach numbers, Reynolds numbers, and ratio of wall-to-adiabatic-wall temperatures (TW/Taw) between this and the Langley 20-inch Mach 6 CF4 Tunnel. This ratio is also matched for Langley's 31-inch Mach 10 Tunnel and is an important parameter useful in the simulation of slender bodies such as National Aerospace Plane (NASP) configurations currently being studied. Having established the nozzle's operating characteristics, the decision was made to install another test section to provide model injection capability. This test section is an open-jet type, with an injection system capable of injecting a model from retracted position to nozzle centerline between 0.5 and 2 seconds. Preliminary calibrations with the new test section resulted in Tunnel blockage. This blockage phenomenon was eliminated when the conical center body in the diffuser was replaced. The issue then, is to provide a new and more efficient variable area diffuser configuration with the capability to withstand testing of larger models without sending the Tunnel into an unstart condition. Use of the 1-dimensional steady flow equation with due regard to friction and heat transfer was employed to estimate the required area ratios (exit area / throat area) in a variable area diffuser. Correlations between diffuser exit Mach number and area ratios, relative to the stagnation pressure ratios and diffuser inlet Mach number were derived. From these correlations, one can set upper and lower operating pressures and temperatures for a given diffuser throat area. In addition, they will provide appropriate input conditions for the full 3-dimensional computational fluid dynamics (CFD) code for further simulation studies.

Loney, Norman W.↗
NASA NTRS1992

Experiments with recursive estimation in astronomical image processing

Recursive estimation concepts were applied to image enhancement problems since the 70's. However, very few applications in the particular area of astronomical image processing are known. These concepts were derived, for 2-dimensional images, from the well-known theory of Kalman filtering in one dimension. The historic reasons for application of these techniques to digital images are related to the images' scanned nature, in which the temporal output of a scanner device can be processed on-line by techniques borrowed directly from 1-dimensional recursive signal analysis. However, recursive estimation has particular properties that make it attractive even in modern days, when big computer memories make the full scanned image available to the processor at any given time. One particularly important aspect is the ability of recursive techniques to deal with non-stationary phenomena, that is, phenomena which have their statistical properties variable in time (or position in a 2-D image). Many image processing methods make underlying stationary assumptions either for the stochastic field being imaged, for the imaging system properties, or both. They will underperform, or even fail, when applied to images that deviate significantly from stationarity. Recursive methods, on the contrary, make it feasible to perform adaptive processing, that is, to process the image by a processor with properties tuned to the image's local statistical properties. Recursive estimation can be used to build estimates of images degraded by such phenomena as noise and blur. We show examples of recursive adaptive processing of astronomical images, using several local statistical properties to drive the adaptive processor, as average signal intensity, signal-to-noise and autocorrelation function. Software was developed under IRAF, and as such will be made available to interested users.

Busko, I.↗
DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

Role of Gravity Waves in Determining Cirrus Cloud Properties

Cirrus clouds are important in the Earth's radiation budget. They typically exhibit variable physical properties within a given cloud system and from system to system. Ambient vertical motion is a key factor in determining the cloud properties in most cases. The obvious exception is convectively generated cirrus (anvils), but even in this case, the subsequent cloud evolution is strongly influenced by the ambient vertical motion field. It is well know that gravity waves are ubiquitous in the atmosphere and occur over a wide range of scales and amplitudes. Moreover, researchers have found that inclusion of statistical account of gravity wave effects can markedly improve the realism of simulations of persisting large-scale cirrus cloud features. Here, we use a 1 -dimensional (z) cirrus cloud model, to systematically examine the effects of gravity waves on cirrus cloud properties. The model includes a detailed representation of cloud microphysical processes (bin microphysics and aerosols) and is run at relatively fine vertical resolution so as to adequately resolve nucleation events, and over an extended time span so as to incorporate the passage of multiple gravity waves. The prescribed gravity waves "propagate" at 15 m s (sup -1), with wavelengths from 5 to 100 km, amplitudes range up to 1 m s (sup -1)'. Despite the fact that the net gravity wave vertical motion forcing is zero, it will be shown that the bulk cloud properties, e.g., vertically-integrated ice water path, can differ quite significantly from simulations without gravity waves and that the effects do depend on the wave characteristics. We conclude that account of gravity wave effects is important if large-scale models are to generate realistic cirrus cloud property climatology (statistics).

OCStarr, David↗
DOE OSTI2026

Spontaneously Broken Noninvertible Symmetries in Transverse-Field Ising Qudit Chains

Recent developments have revealed that symmetries need not form a group, but instead can be noninvertible. Here we use analytical arguments and numerical evidence to illuminate how spontaneous symmetry breaking of a noninvertible symmetry is similar yet distinct from ordinary, invertible, symmetry breaking. We consider one-dimensional chains of group-valued qudits, whose local Hilbert space is spanned by elements of a finite group 𝐺 (reducing to ordinary qubits when 𝐺=ℤ 2 ). We construct Ising-type transverse-field Hamiltonians with Rep⁡(𝐺) symmetry whose generators multiply according to the tensor product of irreducible representations (irreps) of the group 𝐺 . For non-Abelian 𝐺 , the symmetry is noninvertible. In the symmetry broken phase there is one ground state per irrep on a closed chain. The symmetry breaking can be detected by local order parameters but, unlike the invertible case, different ground states have distinct entanglement patterns. We show that for each irrep of dimension greater than one the corresponding ground state exhibits string order, entanglement spectrum degeneracies, and has gapless edge modes on an open chain—features usually associated with symmetry-protected topological order. Consequently, domain wall excitations behave as one-dimensional non-Abelian anyons with nontrivial internal Hilbert spaces and fusion rules. Our Letter identifies properties of noninvertible symmetry breaking that existing quantum hardware can probe.

1-dimensional spin chains↗
NASA NTRS2007

GRMHD and GRPIC Simulations

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous code. The simulation results show the jet formations from a geometrically thin accretion disk near a nonrotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field configuration including issues for future research. A General Relativistic Particle-in-Cell Code (GRPIC) has been developed using the Kerr-Schild metric. The code includes kinetic effects, and is in accordance with GRMHD code. Since the gravitational force acting on particles is extreme near black holes, there are some difficulties in numerically describing these processes. The preliminary code consists of an accretion disk and free-falling corona. Results indicate that particles are ejected from the black hole. These results are consistent with other GRMHD simulations. The GRPIC simulation results will be presented, along with some remarks and future improvements. The emission is calculated from relativistic flows in black hole systems using a fully general relativistic radiative transfer formulation, with flow structures obtained by GRMHD simulations considering thermal free-free emission and thermal synchrotron emission. Bright filament-like features protrude (visually) from the accretion disk surface, which are enhancements of synchrotron emission where the magnetic field roughly aligns with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We would like to extend this research using GRPIC simulations and examine a possible new mechanism for certain X-ray quasi-periodic oscillations (QPOs) observed in blackhole X-ray binaries.

Nishikawa, K.-I.↗
DOE OSTI2024

Slow thermalization and subdiffusion in $U$(1) conserving Floquet random circuits

Random quantum circuits are paradigmatic models of minimally structured and analytically tractable chaotic dynamics. We study a family of Floquet unitary circuits with Haar random U(1) charge conserving dynamics; the minimal such model has nearest-neighbor gates acting on spin-1/2 qubits, and a single layer of even/odd gates repeated periodically in time. We find that this minimal model is not robustly thermalizing at numerically accessible system sizes, and displays slow subdiffusive dynamics for long times. We map out the thermalization dynamics in a broader parameter space of charge conserving circuits, and understand the origin of the slow dynamics in terms of proximate localized and integrable regimes in parameter space. In contrast, we find that small extensions to the minimal model are sufficient to achieve robust thermalization; these include (i) increasing the interaction range to three-site gates, (ii) increasing the local Hilbert space dimension by appending an additional unconstrained qubit to the conserved charge on each site, or (iii) using a larger Floquet period comprised of two independent layers of gates. Furthermore, our results should inform future studies of charge conserving circuits, which are relevant for a wide range of topical theoretical questions.

1-dimensional spin chains↗
DOE OSTI2023

On the relation between the subadditivity cone and the quantum entropy cone

Given a multipartite quantum system, what are the possible ways to impose mutual independence among some subsystems, and the presence of correlations among others, such that there exists a quantum state which satisfies these demands? This question and the related notion of a pattern of marginal independence (PMI) were introduced in [1], and then argued in [2] to be central in the derivation of the holographic entropy cone. Here we continue the general information theoretic analysis of the PMIs allowed by strong subadditivity (SSA) initiated in [1]. We show how the computation of these PMIs simplifies when SSA is replaced by a weaker constraint, dubbed Klein’s condition (KC), which follows from the necessary condition for the saturation of subadditivity (SA). Formulating KC in the language of partially ordered sets, we show that the set of PMIs compatible with KC forms a lattice, and we investigate several of its structural properties. One of our main results is the identification of a specific lower dimensional face of the SA cone that contains on its boundary all the extreme rays (beyond Bell pairs) that can possibly be realized by quantum states. We verify that for four or more parties, KC is strictly weaker than SSA, but nonetheless the PMIs compatible with SSA can easily be derived from the KC-compatible ones. For the special case of 1-dimensional PMIs, we conjecture that KC and SSA are in fact equivalent. To make the presentation self-contained, we review the key ingredients from lattice theory as needed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Challenges in Pulsed-Field Gradient Nuclear Magnetic Resonance on Magnetically Heterogeneous Interfaces: Sequence and Field-Dependent Apparent Diffusion Coefficients

It is well known that the internal gradient (gi) that exists within pores haunts the diffusion coefficient (D) as measured by the pulsed-field gradient (PFG) nuclear magnetic resonance (NMR). Several PFG-NMR methods developed to determine an accurate D were not successful. Then, the steady-state diffusion coefficient (Dapp,8) for the cation [C4mim]+ of [C4mim][Tf2N]; [1-butyl-3-methylimidazolium][bis(trifluoromethylsulfonyl)imde] ionic liquid confined in ordered mesoporous carbon (OMC) were determined by comparing Dapp,8 obtained from 1H PFG-NMR performed with three different stimulated echo sequences: STE, APFG, and MPFG under the two external magnetic field strength, B0 = 9.4 and 14.1 Tesla. The measured Dapp,8 which is an order of magnitude smaller than D of bulk [C4mim][Tf2N], is in good agreement between APFG and MPFG both in B0 = 9.4 and 14.1 Tesla. However, the strong gi artifact, which caused apparent diffusion coefficient (Dapp) depending strongly and weakly on B0 and temperature, respectively, in diffusion-time dependent Dapp, Dapp(?) obtained from a sequence with monopolar gradients (STE) was suppressed by using sequences employing bipolar gradients (APFG and MPFG) in the region of steady-state diffusion. But incompletely suppressed gi artifact resulting in the different behaviors of the early part of Dapp(?) between the sequences leads a ˜ 0.6 and 0.9 in MPFG and APFG, respectively, in the relationship between mean squared displacement and diffusion time: = 2Dta, where a = 0.5 and 1 for 1-dimensional single file diffusion and 3-dimensional bulk diffusion, respectively. The above observations clearly show that the diffusion behavior of ions/molecules within the pores and pore structure, such as the surface-to-volume ratio? (D?_app (?)=D_0 [1-4/(9vp) S/V v(D_0 ?)]) and tortuosity (T = D0/Dapp,8), are possible to be misunderstood, especially in the systems with a non-negligible gi. This work demonstrates that it may be necessary to test several PFG sequences under multiple external magnetic fields for the correct determination of the diffusion behavior of ions/molecules in the pores with a larger internal gradient, gi.

Han, Kee Sung↗
DOE OSTI2025

Single-Walled Zeolitic Nanotube–Poly(oxazoline) Nanocomposites as Heterogeneous Catalysts for Acid–Base Cascade Reactions

Zeolites with a unique, 1-dimensional form factor were recently discovered − zeolite nanotubes (ZNTs). Here we describe the synthesis and characterization of NaH-ZNTpoly( oxazoline) composites targeting liquid-phase acid−base cascade catalysis. NaH-ZNT, a one-dimensional zeolite analogue with mesoporosity (3−4 nm) associated with nanotubes and inherent Brønsted acid sites associated with the microporous zeolite domains, is functionalized with poly(oxazoline)-based triblock copolymers with varying molecular weights (3−17 kDa). The composites are characterized using N2 sorption, STEM, FTIR, and elemental analysis, confirming successful grafting and preservation of the zeolite nanotube structure. The composites’ catalytic performance is evaluated through separate acid and base reactions, followed by a combined cascade of a deacetalization−Knoevenagel condensation for the synthesis of chalcone compounds. High initial reaction rates are demonstrated, but modest overall cascade product formation rates are observed, attributed to interactions between Brønsted acid sites and base amine groups that occur in the polymer-grafted systems. Physical mixtures of NaH-ZNT-SH and poly(oxazoline)s, lacking covalent linkages between ZNT and the polymer, support this supposition. This work demonstrates the potential of NaH-ZNT-poly(oxazoline) composites for liquid-phase cascade catalysis for synthesizing compounds of potential medicinal interest, highlighting the benefits of the grafting-to approach as well as the need for further optimization of the catalytic performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Higher Dimensionality in the Mg–Co–B System: Synthesis and Structure of Incommensurate Composite Mg 1+ε Co 4 B 4

Guided by high-temperature in situ X-ray diffraction, the discovery and synthesis of Mg 1+ε Co 4 B 4 (ε ≈ 0.272) using a MgH 2 hydride precursor is reported, along with a detailed crystal structure description and measurement of magnetic properties. The mismatch in lattice periodicities between Mg and Co–B substructures places Mg 1+ε Co 4 B 4 in the family of incommensurate composite crystals and prompted structural refinement in a (3 + 1)-dimensional model. The structure of Mg 1+ε Co 4 B 4 (P4 2 /ncm(00γ)s00s, a = 6.75847(7) Å, c = 3.94007(8) Å, q = (0, 0, 1.2721(3))) was refined from neutron powder diffraction and high-resolution powder X-ray diffraction data and confirmed by scanning transmission electron microscopy and electron diffraction. Mg 1+ε Co 4 B 4 is isostructural to Nd 1+ε Fe 4 B 4 and several related ternary borides with 0.07 ≤ ε ≤ 0.17, with Mg occupying the rare-earth site. Satellite reflections in the electron diffraction patterns hinted at positional modulation of the transition metal–boron substructure by Mg atoms, but this could not be refined from the neutron or X-ray diffraction data. Low-temperature magnetic measurements show no indications of long-range magnetic ordering or superconductivity down to 5 K. DFT calculations confirmed the absence of a magnetically ordered ground state and the stability of a 5:4 supercell (ε = 0.25) relative to the fully commensurate structure. Neutron diffraction and synthesis from elemental Mg demonstrated that Mg 1+ε Co 4 B 4 is not a hydrogen-stabilized phase. Mg 1+ε Co 4 B 4 represents the second compound reported in the Mg–Co–B system and the first superspace symmetry model of a Nd 1+ε Fe 4 B 4 -type incommensurate composite compound refined from powder diffraction data.

chemical structure↗
DOE OSTI2023

Confinement of many-body Bethe strings

Based on the Bethe ansatz approach and inelastic neutron scattering experiments, we reveal the evolution of confinement of many-body Bethe strings in ordered regions of the quasi-one-dimensional antiferromagnet YbAlO 3 . In the antiferromagnetic phase, the spin dynamics is dominated by confined length-1 Bethe strings, whose dominancy in the high-energy branch of the excitation spectrum yields to confined length-2 Bethe strings when the material is tuned to the spin-density-wave phase. In the thermal-induced disordered region, the confinement effect disappears, and the system restores the conventional quantum integrable physics of the one-dimensional Heisenberg model. Finally, our results establish a unified picture based on a Bethe string for the spin dynamics in different magnetic phases of YbAlO 3 , and thus provide profound insight into many-body quantum magnetism.

1-dimensional spin chains↗
DOE OSTI2025

Holographic View of Mixed-State Symmetry-Protected Topological Phases in Open Quantum Systems

We establish a holographic duality between d -dimensional mixed-state symmetry-protected topological (mSPT) phases and ( d + 1 ) -dimensional subsystem symmetry-protected topological (SSPT) states. Specifically, we show that the reduced density matrix of the boundary layer of a ( d + 1 ) -dimensional SSPT state with subsystem symmetry S and global symmetry G corresponds to a d -dimensional mSPT phase with strong S and weak G symmetries. Conversely, we demonstrate that the wave function of an SSPT state can be constructed by replicating the density matrix of the corresponding lower-dimensional mSPT phase. This mapping links the density matrix in lower dimensions to the entanglement properties of higher-dimensional wave functions, providing an approach for analyzing nonlinear quantities and quantum information metrics in mixed-state systems. Our duality offers a new perspective for studying intrinsic mSPT phases that are unique to open quantum systems, without pure-state analogues. We show that strange correlators and twisted Rényi- N correlators can diagnose these nontrivial phases and we explore their connection to strange correlators in pure-state SSPT phases. Furthermore, we discuss several implications of this holographic duality, including a method for preparing intrinsic mSPT phases through the duality. Published by the American Physical Society 2025

Sun, Shijun (ORCID:0000000168880030)↗
DOE OSTI2023

Signal Processing of Multiplexed Optical PWM Signals for Sensor Arrays in Nuclear Environments

Safe and effective generation of terrestrial nuclear power greatly benefits from the actionable data provided by the array of sensors located throughout a plant to provide a holistic online indication of reactor operation. This array includes sensors for monitoring coolant flow and pressure, temperature and heat transfer, radiation levels, structure health monitoring, and other critical parameters for reactor operation. While sensors capable of measuring these parameters have been developed, the electronics used in the pre-amplification and analog-to-digital conversion of the small signals they produce are extremely sensitive and susceptible to damage by high-temperatures and radiation environments nuclear reactors encounter while generating power. The small signals from sensors in nuclear power plants (NPPs) are transmitted over long cable runs which introduce dispersion artifacts into the signals of interest as well as electromagnetic interference (EMI) from lighting fixtures, pumps, mains electricity, and other equipment. To overcome these challenges, a front-end digitization (FREND) platform has been developed to use radiation-tolerant electronics to multiplex, amplify, and optically encode signals from an array of sensors for transmission over an optical fiber to mitigate dispersion and EMI artifacts from long runs of electrical cabling. To recover the optically transmitted data, a signal processing scheme based on 1-dimensional template matching to an indexing channel is described herein and demonstrated to have an effective bit-depth of 9.2 bits (1%). This scheme has been validated in proof-of-concept, non-nuclear testing and preliminary experimental results show good agreement between the measured optical output and and sensor input signals. The FREND system represents a low-loss data link between sensors in nuclear environments and data acquisition hardware which is aimed at improving the signal-to-noise ratio of the data acquired from these sensors to provide better information to operators and researchers.

Sweeney, Dan↗
NASA NTRS1973

Multispectral combination and display of ERTS-1 data

A significant problem in the use of ERTS-1 data is the extraction of information pertinent to each application and the presentation of that information in a form most suitable to users. When the information is to be displayed for visual study by an observer, the problem can be reduced to two steps: (1) Dimensionality reduction, an objective procedure which attempts to preserve most of the ERTS-1 information in a smaller number of components. (2) Display of the reduced number of components for optimum visibility by an observer. A specific dimensionality reduction technique has been applied to ERTS-1 data for several geographical areas in California and distinct types of earth resources. In the display of the reduced number of components, consideration has to be given to properties of the human visual system and the statistics of the data to be displayed. Previous work on digital image enhancement was applied to this problem to generate color composites which contain and display most of the information provided by the ERTS-1 sensors. Results of this approach were interesting, both in terms of the small mean-square caused by the dimensionality reduction, as well as for the examples of enhanced images that have been obtained.

Algazi, V. R.↗
NASA NTRS1985

A cluster analysis method for identification of subpopulations of cells in flow cytometric list-mode arrays

A specialized program was developed for flow cytometric list-mode data using an heirarchical tree method for identifying and enumerating individual subpopulations, the method of principal components for a two-dimensional display of 6-parameter data array, and a standard sorting algorithm for characterizing subpopulations. The program was tested against a published data set subjected to cluster analysis and experimental data sets from controlled flow cytometry experiments using a Coulter Electronics EPICS V Cell Sorter. A version of the program in compiled BASIC is usable on a 16-bit microcomputer with the MS-DOS operating system. It is specialized for 6 parameters and up to 20,000 cells. Its two-dimensional display of Euclidean distances reveals clusters clearly, as does its 1-dimensional display. The identified subpopulations can, in suitable experiments, be related to functional subpopulations of cells.

Li, Z. K.↗
NASA NTRS1993

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.↗