SEARCH · Search NASA
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.
The widely transonic swept-wing code, FL022, was found to have an error in the transformed flow equation in the computational domain. The revised version of the code correctly accounted for the non-straight leading edge geometry and its effect on the pressure distribution.
The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.
The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.
ANISAP is a 3-D finite element FORTRAN 77 computer code for linear elastic, small strain, analysis of laminated composites with arbitrary geometry including free edges and holes. Individual layers may be isotropic or transversely isotropic in material principal coordinates; individual layers may be rotated off-axis about a global z-axis. The laminate may be a hybrid. Three different isoparametric elements, variable order of gaussian integration, calculation of stresses at element boundaries, and loading by either nodal displacement of forces are included in the program capability. Post processing capability includes failure analysis using the tensor polynominal failure criterion.
Crack propagation in a rotating inner raceway of a high speed roller bearing is analyzed using the boundary integral equation method. The method consists of an edge crack in a plate under tension, upon which varying Hertzian stress fields are superimposed. A computer program for the boundary integral equation method was written using quadratic elements to determine the stress and displacement fields for discrete roller positions. Mode I and Mode II stress intensity factors and crack extension forces G sub 00 (energy release rate due to tensile opening mode) and G sub r0 (energy release rate due to shear displacement mode) were computed. These calculations permit determination of that crack growth angle for which the change in the crack extension forces is maximum. The crack driving force was found to be the alternating mixed-mode loading that occurs with each passage of the most heavily loaded roller. The crack is predicted to propagate in a step-like fashion alternating between radial and inclined segments, and this pattern was observed experimentally. The maximum changes DeltaG sub 00 and DeltaG sub r0 of the crack extension forces are found to be good measures of the crack propagation rate and direction.
A vectorized 3D Navier-Stokes code has been implemented on CYBER 205 for solving the supersonic laminar flow over a swept fin/flat plate junction. The code extends MacCormack's predictor-corrector finite volume scheme to a generalized coordinate system in a locally one dimensional time split fashion. A systematic parametric study is conducted to examine the effect of fin sweep on the computed flow field. Calculated results for the pressure distribution on the flat plate and fin leading edge are compared with the experimental measurements of a right angle blunt fin/flat plate junction. The decrease in the extent of the separated flow region and peak pressure on the fin leading edge, and weakening of the two reversed supersonic zones with increase in fin sweep have been clearly observed in the numerical simulation.
Aluminum reinforced with carbon fibers superior to copper in some respects. Lightweight composite material has high thermal conductivity. Consists of aluminum matrix containing graphite fibers, all oriented in same direction. Available as sheets, tubes, and bars. Thermal conductivity of composite along fibers rises above that of pure copper over substantial range of temperatures. Graphite/aluminum composite useful in variety of heat-transfer applications in which reduction of weight critical. Used to conduct heat in high-density, high-speed integrated-circuit packages for computers and in base plates for electronic equipment. Also used to carry heat away from leading edges of wings in high-speed airplanes.
An engineering prediction method and associated computer code VTXCHN to predict nose vortex shedding from circular and noncircular forebodies with sharp chine edges in subsonic flow at angles of attack and roll are presented. Axisymmetric bodies are represented by point sources and doublets, and noncircular cross sections are transformed to a circle by either analytical or numerical conformal transformations. The lee side vortex wake is modeled by discrete vortices in crossflow planes along the body; thus the three-dimensional steady flow problem is reduced to a two-dimensional, unsteady, separated flow problem for solution. Comparison of measured and predicted surface pressure distributions, flow field surveys, and aerodynamic characteristics are presented for noncircular bodies alone and forebodies with sharp chines.
An engineering prediction method and associated computer code VTXCHN to predict nose vortex shedding from circular and noncircular forebodies with sharp chines edges in subsonic flow at angles of attack and roll are presented. Axisymmetric bodies are represented by point sources and doublets, and noncircular cross sections are transformed to a circle by either analytical or numerical conformal transformations. The lee side vortex wake is modeled by discrete vortices in crossflow planes along the body; thus the three-dimensional steady flow problem is reduced to a two-dimensional, unsteady, separated flow problem for solution. Comparison of measured and predicted pressure distributions, flow field surveys, and aerodynamic characteristics are presented for noncircular bodies alone and forebodies with sharp chines.
The past several decades have seen the rise of large, highly interactive systems that are on the forward edge of technology. As a result of this growth and the increased usage of digital systems (computers and software), the concept of systems engineering has gained increasing attention. Some of this attention is no doubt due to large program failures which possibly could have been avoided, or at least mitigated, through the use of systems engineering principles. The complexity of modern day weapon systems requires conscious application of systems engineering concepts to ensure producible, operable and supportable systems that satisfy mission requirements. Although many authors have traced the roots of systems engineering to earlier dates, the initial formalization of the systems engineering process for military development began to surface in the mid-1950s on the ballistic missile programs. These early ballistic missile development programs marked the emergence of engineering discipline 'specialists' which has since continued to grow. Each of these specialties not only has a need to take data from the overall development process, but also to supply data, in the form of requirements and analysis results, to the process. A number of technical instructions, military standards and specifications, and manuals were developed as a result of these development programs. In particular, MILSTD-499 was issued in 1969 to assist both government and contractor personnel in defining the systems engineering effort in support of defense acquisition programs. This standard was updated to MIL-STD499A in 1974, and formed the foundation for current application of systems engineering principles to military development programs.
The boundary layer receptivity process due to the interaction of three-dimensional slow and fast acoustic disturbances with a blunted flat plate is numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 10(exp 6)/inch. The computations are performed with and without two-dimensional isolated roughness element located near the leading edge. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The simulations showed that the linear instability waves are generated very close to the leading edge. The wavelength of the disturbances inside the boundary layer first increases gradually and becomes longer than the wavelength for the instability waves within a short distance from the leading edge. The wavelength then decreases gradually and merges with the wavelength for the Tollmien-Schlichting wave. The initial amplitudes of the instability waves near the neutral points, the receptivity coefficients, are about 1.20 and 0.07 times the amplitude of the free-stream disturbances for the slow and the fast waves respectively. It was also revealed that small isolated roughness element does not enhance the receptivity process for the given nose bluntness.
The boundary layer receptivity process due to the interaction of three-dimensional slow and fast acoustic disturbances with a blunted flat plate is numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 106/inch. The computations are performed with and without two-dimensional isolated roughness element located near the leading edge. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The simulations showed that the linear instability waves are generated very close to the leading edge. The wavelength of the disturbances inside the boundary layer first increases gradually and becomes longer than the wavelength for the instability waves within a short distance from the leading edge. The wavelength then decreases gradually and merges with the wavelength for the Tollmien_Schlichting wave. The initial amplitudes of the instability waves near the neutral points, the receptivity coefficients, are about 1.20 and 0.07 times the amplitude of the free-stream disturbances for the slow and the fast waves respectively. It was also revealed that small isolated roughness element does not enhance the receptivity process for the given nose bluntness.
Arc jet Computer Vision (arcjetCV) is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.
Arc jet Computer Vision (arcjetCV) is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.
For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.
Two research programs are described which directly relate to the application of natural laminar flow (NLF) technology to transonic transport-type wind planforms. Each involved using state-of-the-art computational methods to design three-dimensional wing contours which generate significant runs of favorable pressure gradients. The first program supported the Variable Sweep Transition Flight Experiment and involves design of a full-span glove which extends from the leading edge to the spoiler hinge line on the upper surface of an F-14 outer wing panel. Boundary-layer and static-pressure data will be measured on this design during the supporting wind-tunnel and flight tests. These data will then be analyzed and used to infer the relationship between crossflow and Tollmein-Schlichting disturbances on laminar boundary-layer transition. A wing was designed computationally for a corporate transport aircraft in the second program. The resulting wing design generated favorable pressure gradients from the leading edge aft to the mid-chord on both upper and lower surfaces at the cruise design point. Detailed descriptions of the computational design approach are presented along with the various constraints imposed on each of the designs. Wing surface pressure distributions, which support the design objective and were derived from transonic three-dimensional analyses codes, are also presented. Current status of each of the research programs is included in the summary.
Since its inception 30 years ago, LAMMPS has grown to be a world-class molecular dynamics code and a cornerstone of computational materials science research. This project aimed to keep LAMMPS at the forefront of molecular dynamics simulations by adapting LAMMPS to the latest developments in machine learning technology and hardware. Initially, the project set out to provide a unified implementation of active learning for efficient training data generation in LAMMPS, but the research trajectory pivoted to address more immediate and impactful opportunities. On the hardware side, recent record-breaking molecular dynamics simulations were developed on the Cerebras wafer-scale AI chip, and this project has developed an interface between LAMMPS and the hardware-specific molecular dynamics code to accelerate and simplify development and user adoption. On the software side, PyTorch’s Ahead-of-Time (AOT) compilation features promised increased performance for state-of-the-art equivariant neural network potentials, and this project laid the groundwork for their adoption in LAMMPS, resulting in a nearly 20x acceleration in extreme cases. Combined with a comprehensive benchmark study of LAMMPS across all current exascale systems, this project has reinforced LAMMPS’s role as a versatile, high-performance tool for current and future materials science applications.
To better understand and address the challenges faced in computing the acoustics of flow fields, test problems must be considered. In the present study, the sound radiated by the interaction of a flat plate with an oncoming gust containing a two component, mean velocity is computed. The gust has a uniform mean flow in x with Mach number M(infinity) equal to 0.5. The gust's mean velocity in y is of smaller amplitude and is given by: v = 0.1 sin(pi/8(x/M(sub infinity) - t)). This problem has been posed for an upcoming ICASE/LaRC workshop on benchmark problems in computational aeroacoustics. A plate with a length of 30 units in x is used. The plate is assumed to be infinitesimally thin and is centered at the origin. All variables are made dimensionless using the scales specified. Acoustic quantities are obtained by numerically integrating the linearized Euler equations. Integration is performed on the computational domain -100.0 less than or equal to x less than or equal to 100.0, -100.0 less than or equal to y less than or equal to 100.0, using unit length grid spacing in x and in y. An integration scheme is sought which will provide accurate solution to the small quantities of interest at a minimal computational expense. Results indicate that with the given discretization a scheme of minimal fourth order accuracy might be adequate to approximate the waves within the given flow. Thus, a variation of the MacCormack scheme with fourth order accuracy in space and second order accuracy in time was chosen. A scheme with sixth order accuracy in space has also been implemented and results compared with those of the fourth order accurate scheme. To ensure no mass flux, zero normal velocity is assigned at the plate. This condition will induce a discontinuity in the pressure across the plate location. Values for the perturbation pressure p' along the surface of the plate are obtained using a one-sided, third order Taylor expansion, such that p'(sub y) = O. In accordance with the Kutta condition, perturbation pressure at the trailing edge is assigned to zero. In the far field, radiation boundary conditions have been implemented. The effectiveness of the far field conditions are validated by computing in a larger computational domain and comparing the results. Early in time, sound waves begin to radiate from the plate. A Doppler effect is observed. After the initial transients disappear, the strongest waves leave the trailing edge at an approximate 45 degree angle. The intensity pattern of pressure fluctuations shows five lobes (of increasing magnitude with increasing downstream direction) emerging from the plate. Undesirable short wave contaminants are observed in the computed pressure distribution along the plate surface. For a more accurate solution at small scales, a more refined discretization will be required.