Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “H CODES”

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 289 records · Page 16

Prediction of drag at subsonic and transonic speeds using Euler methods

A technique for the evaluation of aerodynamic drag from flowfield solutions based on the Euler equations is discussed. The technique is limited to steady attached flows around three-dimensional configurations in the absence of active systems such as surface blowing/suction and propulsion. It allows the decomposition of the total drag into induced drag and wave drag and, consequently, it provides more information on the drag sources than the conventional surface-pressure integration technique. The induced drag is obtained from the integration of the kinetic energy (per unit distance) of the trailing vortex system on a wake plane and the wave drag is obtained from the integration of the entropy production on a plane just downstream of the shocks. The drag-evaluation technique is applied to three-dimensional flowfield solutions for the ONERA M6 wing as well as an aspect-ratio-7 wing with an elliptic spanwise chord distribution and an NACA-0012 section shape. Comparisons between the drag obtained with the present technique and the drag based on the integration of surface pressures are presented for two Euler codes.

Nikfetrat, K.↗

Gross-fragmentation of meteoroids and bulk density of Geminids from photographic fireball records

The explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid moving in any actual atmosphere was published several years ago. The solution yields the theoretical relation of l, the distance flown by the meteoroid in its trajectory, as a function of time, t, assuming that the height, h, is a known function of l. The photographic records of meteors and fireballs are coded by time marks, using a rotating shutter or a similar device to break the moving image. Time is, thus, the independent variable and for each time mark on a meteoroid trajector, the observed distance along the trajectory, l sub obs, as well as the double- or multiple- station photographs of the same meteoroid. Applying this solution to all available Prairie Network (PN) fireball-records, we recognized that the majority of them gave good solutions with standard deviations somewhat bigger than the intrinsic geometrical precision of the data. We also noticed that, on an average, previous methods of evaluation of the meteoroid velocities (interpolation polynomials, numerical differenciation of the observed l sub obs) used up to only several tens of percent of the intrinsic precision of the PN observational data. When residuals of these solutions, i.e. l sub obs - l sub com, were represented as a function of time for about 75 percent of solutions. The remaining 25 percent of residuals showed systematic changes with time exceeding one standard deviation. We tried to explain these systematic time course of residuals by using different meteoroids first computed theoretically and then analyzed by the same model as the natural PN fireballs were. The conclusion of these model computations is that systematic time changes of residuals in the nonfragmenting model exceeding one standard of deviation are caused by sudden gross fragmentation at one or more trajectory points. Thus, we generalized the explicit solution of the drag and ablation equations of a single nonfragmenting meteoroid by allowing for one or more points, where sudden gross fragmentation can occur. Using this generalized solution, the distances along the meteoroid trajectory can be computed for any choice of input parameters and compared with the observed distances flown by the meteoroid. For the most precise and long fireball trajectories, the least-squares solution can, thus, yield the initial velocities, the ablation coefficients, the dynamical masses, the positions of gross-fragmentation points, and the terminal mass. At a gross-fragmentation point, the ratio of the main mass to all the remaining fragments can be compared with the dynamic mass determined from our gross-fragmentation model and, thus, the meteoroid bulk density can be evaluated. We applied the gross-fragmentation model to sever PN fireballs showing time changes of residuals, and we recognized that, in all these cases, the new computed bulk densities of meteoroids resulted higher in comparison with the meteoroid densities determined with the non-gross-fragmentation model. Other aspects of the study are discussed.

Ceplecha, Zdenek↗

Explicit robust schemes for implementation of a class of principal value-based constitutive models: Symbolic and numeric implementation

The issue of developing effective and robust schemes to implement a class of the Ogden-type hyperelastic constitutive models is addressed. To this end, special purpose functions (running under MACSYMA) are developed for the symbolic derivation, evaluation, and automatic FORTRAN code generation of explicit expressions for the corresponding stress function and material tangent stiffness tensors. These explicit forms are valid over the entire deformation range, since the singularities resulting from repeated principal-stretch values have been theoretically removed. The required computational algorithms are outlined, and the resulting FORTRAN computer code is presented.

Arnold, S. M.↗

Optimality Of Variable-Length Codes

Report presents analysis of performances of conceptual Rice universal noiseless coders designed to provide efficient compression of data over wide range of source-data entropies. Includes predictive preprocessor that maps source data into sequence of nonnegative integers and variable-length-coding processor, which adapts to varying entropy of source data by selecting whichever one of number of optional codes yields shortest codeword.

Yeh, Pen-Shu↗

Stellar models with microscopic diffusion and rotational mixing. 1: Application to the Sun

The Yale stellar evolution code has been modified to include the combined effects of diffusion and rotational mixing on H-1, H-4, and the trace elements He, Li-6, Li-7, and Be-9. The interaction between rotational mixing and diffusion is studied by calculating a number of calibrated solar models. The rotational mixing inhibits the diffusion in the outer parts of the models, leading to a decrease in the envelope diffusion by 25%-50%. Conversely, diffusion leads to gradients in mean molecular weight which can inhibit the rotational mixing. The degree to which gradients in mean molecular weight inhibit the rotational mixing is somewhat uncertain. A comparison with the observed solar oblateness suggests that gradients in the mean molecular weight play a smaller role in inhibiting the rotational mixing than previously believed. This is reinforced by the fact that the model with the standard value for the inhibiting effect of mean molecular weight on the rotational mixing depletes no Li on the main sequence. This is clear in contrast to the observations. A reduction in the inhibiting effect of mean molecular weight gradients by a factor of 10 loads to noticeable main-sequence Li depletion.

Chaboyer, Brian↗

Ground-Testing of One-Way Ranging from a Lunar Beacon Demonstrator Payload (Lunar Node -1)

Similar to lighthouses providing navigation aids and awareness to ships sailing along the coast, in-situ navigation beacons can be used to provide redundancy, awareness, and improved onboard position knowledge for vehicles operating in and around the lunar sphere of influence. Lunar Node - 1 is a radio frequency navigation beacon that is manifested on Intuitive Machines' NOVA-C lunar lander flight in 2022 for operation during cruise and from a mid-latitude landing location. This payload utilizes commercial components developed for Low Earth Orbit applications and packages them into an integrated package to provide a navigation reference signal as part of a large lunar-centric navigation network. As part of its operations, Deep Space Network ground stations will operate as the notional "user" of this service, allowing for performance characterization, calibration, and operation of the beacon payload. This paper provides an overview of the payload, how its fits into a large concept of operations, and ground testing results. A specific focus is on ground testing of the payload with DSN ground receivers as part of Radio Frequency (RF) Compatibility Testing. A primary goal of this testing was to characterize a pseudo-noise ranging code sequence generated by the payload. The results in this paper define the testing sequence, unique operational constraints with two-way ranging ground equipment, and applications of the ground-derived error model to predicted in-flight performance of the ranging code.

Evan J Anzalone↗

Mshpy23: A User-Friendly, Parameterized Model of Magnetosheath Conditions

Lunar Environment heliospheric X-ray Imager (LEXI) and Solar wind – Magnetosphere - Ionosphere Link Explorer (SMILE) will observe magnetosheath and its boundary motion in soft X-rays for understanding magnetopause reconnection modes under varioussolar wind conditions after their respective launches in 2024 and 2025. Magnetosheath conditions, namely, plasma density, velocity, and temperature, are key parameters for predicting and analyzing soft X-ray images from the LEXI and SMILE missions. We developed a user-friendly model of magnetosheath that parameterizes number density, velocity, temperature, and magnetic field by utilizing the global Magnetohydrodynamics (MHD) model as well as the pre-existing gas-dynamic and analytic models. Using this parameterized magnetosheath model, scientists can easily reconstruct expected soft X-ray images and utilize them for analysis of observed images of LEXI and SMILE without simulating the complicated global magnetosphere models. First, we created an MHD based magnetosheath model by running a total of 14 OpenGGCM global MHD simulations under 7 solar wind densities (1, 5, 10, 15, 20, 25, and 30cm−3) and 2 interplanetary magnetic field BZ components (± 4nT), and then parameterizing the results in new magnetosheath conditions. We compared the magnetosheath model result with THEMIS statistical data and it showed good agreement with a weighted Pearson correlation coefficient greater than 0.77, especially for plasma density and plasma velocity. Second, we compiled a suite of magnetosheath models incorporating previous magnetosheath models (gas-dynamic, analytic), and did two case studies to test the performance. The MHD based model was comparable to or better than the previous models while providing self consistency among the magnetosheath parameters. Third, we constructed a tool to calculate a soft X-ray image from any given vantage point, which can support the planning and data analysis of the aforementioned LEXI and SMILE missions. A release of the code has been uploaded to a Github repository.

Magnetosheath↗

DUKSUP: A Computer Program for High Thrust Launch Vehicle Trajectory Design and Optimization

From the late 1960s through 1997, the leadership of NASAs Intermediate and Large class unmanned expendable launch vehicle projects resided at the NASA Lewis (now Glenn) Research Center (LeRC). One of LeRCs primary responsibilities --- trajectory design and performance analysis --- was accomplished by an internally-developed analytic three dimensional computer program called DUKSUP. Because of its Calculus of Variations-based optimization routine, this code was generally more capable of finding optimal solutions than its contemporaries. A derivation of optimal control using the Calculus of Variations is summarized including transversality, intermediate, and final conditions. The two point boundary value problem is explained. A brief summary of the codes operation is provided, including iteration via the Newton-Raphson scheme and integration of variational and motion equations via a 4th order Runge-Kutta scheme. Main subroutines are discussed. The history of the LeRC trajectory design efforts in the early 1960s is explained within the context of supporting the Centaur upper stage program. How the code was constructed based on the operation of the AtlasCentaur launch vehicle, the limits of the computers of that era, the limits of the computer programming languages, and the missions it supported are discussed. The vehicles DUKSUP supported (AtlasCentaur, TitanCentaur, and ShuttleCentaur) are briefly described. The types of missions, including Earth orbital and interplanetary, are described. The roles of flight constraints and their impact on launch operations are detailed (such as jettisoning hardware on heating, Range Safety, ground station tracking, and elliptical parking orbits). The computer main frames on which the code was hosted are described. The applications of the code are detailed, including independent check of contractor analysis, benchmarking, leading edge analysis, and vehicle performance improvement assessments. Several of DUKSUPs many major impacts on launches are discussed including Intelsat, Voyager, Pioneer Venus, HEAO, Galileo, and Cassini.

Launch vehicle performance optimization↗

Carbon and nitrogen in Type 2 supernova diamonds

Abundant diamonds found in meteorites seem either to have condensed within supernova interiors during their expansions and coolings or to have been present around those explosions. Either alternative allows implantation of Xe-HL prior to interstellar mixing. A puzzling feature is the near normalcy of the carbon isotopes, considering that the only C-rich matter, the He-burning shell, is pure C-12 in that region. That last fact has caused many to associate supernova carbon with C-12 carbon, so that its SUNOCONS have been anticipated as very C-12-rich. We show that this expectation is misleading because the C-13-rich regions of Type 2's have been largely overlooked in this thinking. We here follow the idea that the diamonds nucleated in the C-12-rich He shell, the only C-rich site for nucleation, but then attached C-13-rich carbon during turbulent encounters with overlying C-13-rich matter. That is, the initial diamonds continued to grow during the same collisional encounters that cause the Xe-HL implantation. Instead of interacting with the small carbon mass having 13/12 = 0.2 in the upper He zone, however, we have calculated the remnants of the initial H-burning core, which left behind C-13-rich matter as it receded during core hydrogen burning. Howard et al. described why the velocity mixing would be essential to understanding the implantation of both the Xe-H and Xe-L components. Velocity mixing is now known to occur from the X-ray and gamma-ray light curves of supernova 1987A. Using the stellar evolution code developed at Goettingen, we calculated at Clemson the evolution of a grid of massive stars up to the beginning of core He burning. We paid attention to all H-burning reactions throughout the star, to the treatment of both convection and semiconvection, and to the recession of the outer boundary of the convective H-burning core as the star expands toward a larger redder state. This program was to generate a careful map of the CNO isotope distribution as He burning begins. Our result for the 30 solar mass star is shown.

Clayton, Donald D.↗

The spatial correlation properties of dark galaxy halos in a CDM universe

We use the Hierarchical Particle Mesh (HPM) N-body code written by J. V. Villumsen (Villumsen, 1989) to investigate the two-point spatial correlation function, xi(r), of dark galaxy halos as a function of halo mass and local environment (i.e. high, low, or average mass density). We assume a standard cold dark matter (CDM) universe (omega = 1, delta = 0, H sub 0 = 50,km/sec/Mpc). Because of the large dynamic ranges in mass and length that can be obtained with the HPM code, it is well-suited to an investigation of this sort.

Brainerd, Tereasa G.↗

Finite volume calculation of three-dimensional potential flow around a propeller

The finite volume scheme of Jameson (1977) is used to calculate potential flow around a propeller rotating at high speed. An H-type mesh is generated and used successfully in the calculations. A test calculation with a thick blade cross section shows that the present code is capable of computing the propeller flow at the advance Mach number 0.8. The possible physical mechanisms which may play an important role in the propeller aerodynamics are discussed.

Jou, W.-H.↗

Comparative study of x ray and microwave emissions during solar flares

The work supported by the grant consisted of two projects. The first project involved making detailed case studies of two flares using SMM data in conjunction with ground based observations. The first flare occurred at 1454 UT on June 20, 1989 and involved the eruption of a prominence near the limb. In the study we used data from many wavelength regimes including the radio, H-alpha, hard X-rays, and soft X-rays. We used a full gyrosynchrotron code to model the apparent presence of a 1.4 GHz source early in the flare that was in the form of a large coronal loop. The model results lead us to conclude that the initial acceleration occurs in small, dense loops which also produced the flare's hard X-ray emission. We also found evidence that a source at 1.4 GHz later in the event was due to second harmonic plasma emission. This source was adjacent to a leg of the prominence and comes from a dense column of material in the magnetic structure supporting the prominence. Finally, we investigated a source of microwaves and soft X-rays, occurring approximately 10 min after the hard X-ray peak, and calculate a lower limit for the density of the source. The second flare that was studied occurred at 2156 UT on June 20, 1989 and was observed with the VLA and the Owens Valley Radio Observatory (OVRO) Frequency Agile Array. We have developed a gyrosynchrotron model of the sources at flare peak using a new gyrosynchrotron approximation which is valid at very low harmonics of the gyrofrequency. We found that the accelerated particle densities of the sources decreased much more with radius from the source center than had been supposed in previous work, while the magnetic field varied less. We also used the available data to analyze a highly polarized source which appeared late in the flare. The second project involved compiling a statistical base for the relative timing of the hard X-ray peak, the turbulent and blue-shift velocities inferred from soft X-ray line emissions observed by SMM and the microwave peak as determined from ground-based observations. This timing was then used to aid the testing of newly developed global models for flares that incorporate the global magnetic topology as well as the electron dynamics that are responsible for the hard X-rays and microwaves.

Winglee, Robert M.↗

A Study of Ablation-Flowfield Coupling Relevant to the Orion Heatshield

The coupled interaction between an ablating surface and the surrounding aerothermal environment is studied. An equilibrium ablation model is coupled to the LAURA flowfield solver, which allows the char ablation rate (m-dot(sub c)) to be computed as part of the flowfield solution. The wall temperature (T(sub w)) and pyrolysis ablation rate (m-dot(sub g)) may be specified by the user, obtained from the steady-state ablation approximation, or computed from a a material response code. A 32 species thermochemical nonequilibrium flowfield model is applied, which permits the treatment of C, H, O, N, and Si containing species. Coupled ablation cases relevant to the Orion heatshield are studied. These consist of diffusion-limited oxidation cases with Avcoat as the ablation material. The m-dot(sub c) values predicted from the developed coupled ablation analysis were compared with those obtained from a typical uncoupled ablation analysis. The coupled results were found to be as much as 50% greater than the uncoupled values. This is shown to be a result of the cumulative effect of the two fundamental approximations inherent in the uncoupled analysis.

Johnston, Christopher O.↗

Electromagnetic scattering by coated convex surfaces and wedges simulated by approximate boundary conditions

Asymptotic/high-frequency solutions are developed for analyzing the non-specular scattering mechanisms associated with coated convex surfaces and edges simulated by approximate boundary conditions. In particular, the standard impedance boundary conditions (SIBC's) and the second order generalized impedance boundary conditions (GIBC's) are employed for a characterization of the edge diffraction, creeping wave, and surface diffracted wave contributions. To study the creeping wave and surface diffracted wave mechanisms, rigorous UTD (uniform geometrical theory of diffraction) diffraction coefficients are developed for a convex coated cylinder simulated with SIBC's and GIBC's. The ray solutions obtained remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. The diffraction coefficient for a convex cylinder are obtained via a generalization of the corresponding ones of the circular cylinder. To validate the asymptotic/high-frequency solution, integral equations are derived for both E and H-polarization and solved numerically using the method of moments. Results are presented for a single and three layered coated convex cylinder. Some insights are also provided on the accuracy of the employed GIBC's versus SIBC's for application to curved surfaces. To characterize the scattering by impedance wedges illuminated at skew incidence, diffraction coefficients are derived from an approximate solution of the governing functional difference equations. This solution exactly recovers the known ones for an impedance half plane or an arbitrary wedge at normal incidence, and to validate it for other wedge angles, a moment method code was used. Finally, to test the usefulness of the approximate skew incidence impedance wedge diffraction coefficient for three dimensional structures, equivalent currents are derived in the context of the physical theory of diffraction (PTD) for a finite length impedance wedge of arbitrary internal angles. These are incorporated in a standard general purpose PTD code and results are presented for a number of different impedance structures.

Syed, H. H.↗

Boundary layer integral matrix procedure code modifications and verifications

A summary of modifications to Aerotherm's Boundary Layer Integral Matrix Procedure (BLIMP) code is presented. These modifications represent a preliminary effort to make BLIMP compatible with other JANNAF codes and to adjust the code for specific application to rocket nozzle flows. Results of the initial verification of the code for prediction of rocket nozzle type flows are discussed. For those cases in which measured free stream flow conditions were used as input to the code, the boundary layer predictions and measurements are in excellent agreement. In two cases, with free stream flow conditions calculated by another JANNAF code (TDK) for use as input to BLIMP, the predictions and the data were in fair agreement for one case and in poor agreement for the other case. The poor agreement is believed to result from failure of the turbulent model in BLIMP to account for laminarization of a turbulent flow. Recommendations for further code modifications and improvements are also presented.

Evans, R. M.↗