Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Analysis of unsteady compressible viscous layers

The development of an analysis to predict the unsteady compressible flows in blade boundary layers and wakes is presented. The equations that govern the flows in these regions are transformed using an unsteady turbulent generalization of the Levy-Lees transformation. The transformed equations are solved using a finite difference technique in which the solution proceeds by marching in time and in the streamwise direction. Both laminar and turbulent flows are studied, the latter using algebraic turbulence and transition models. Laminar solutions for a flat plate are shown to approach classical asymptotic results for both high and low frequency unsteady motions. Turbulent flat-plate results are in qualitative agreement with previous predictions and measurements. Finally, the numerical technique is also applied to the stator and rotor of a low-speed turbine stage to determine unsteady effects on surface heating. The results compare reasonably well with measured heat transfer data and indicate that nonlinear effects have minimal impact on the mean and unsteady components of the flow.

Power, G. D.↗
NASA NTRS2004

Subsurface Stress Fields in FCC Single Crystal Anisotropic Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent high cycle fatigue (HCF) failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and non-crystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is presented for evaluating the subsurface stresses in the elastic half-space, based on the adaptation of a stress function method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis (FEA). Effects of crystal orientation on stress response and fatigue life are examined. Obtaining accurate subsurface stress results for anisotropic single crystal contact problems require extremely refined three-dimensional (3-D) finite element grids, especially in the edge of contact region. Obtaining resolved shear stresses (RSS) on the principal slip planes also involves considerable post-processing work. For these reasons it is very advantageous to develop analytical solution schemes for subsurface stresses, whenever possible.

Arakere, Nagaraj K.↗
NASA NTRS1974

Evaluation of analytical procedures for prediction of turbulent boundary layers on a porous wall

An analytical study has been made to determine how well current boundary layer prediction techniques work when there is mass transfer normal to the wall. The data that were considered in this investigation were for two-dimensional, incompressible, turbulent boundary layers with suction and blowing. Some of the bleed data were taken in an adverse pressure gradient. An integral prediction method was used three different porous wall skin friction relations, in addition to a solid-surface relation for the suction cases. A numerical prediction method was also used. Comparisons were made between theoretical and experimental skin friction coefficients, displacement and momentum thicknesses, and velocity profiles. The integral method with one of the porous wall skin friction laws gave very good agreement with data for most of the cases considered. The use of the solid-surface skin friction law caused the integral to overpredict the effectiveness of the bleed. The numerical techniques also worked well for most of the cases.

Towne, C. E.↗
NASA NTRS1972

A discrete variable approximation to minimum weight panel designs subject to a supersonic flutter-speed constraint.

A numerical technique is presented that determines the minimum-weight thickness distribution for a one-dimensional, simply supported panel subjected to a parallel supersonic flow on one side. An aeroelastic eigenvalue, which characterizes flutter speed, is held constant. The governing differential equations are approximated by sets of difference equations adjoined to the weight function via a penalty function. A conjugate gradient method then solves the resulting sequence of unconstrained minimization problems. Numerical results are obtained without imposing a minimum thickness constraint; the effect of constant inplane tensile and compressive stresses is also investigated.

Pierson, B. L.↗
NASA NTRS1975

A nonlinear model for gas chromatograph systems

Fundamental engineering design techniques and concepts were studied for the optimization of a gas chromatograph-mass spectrometer chemical analysis system suitable for use on an unmanned, Martian roving vehicle. Previously developed mathematical models of the gas chromatograph are found to be inadequate for predicting peak heights and spreading for some experimental conditions and chemical systems. A modification to the existing equilibrium adsorption model is required; the Langmuir isotherm replaces the linear isotherm. The numerical technique of Crank-Nicolson was studied for use with the linear isotherm to determine the utility of the method. Modifications are made to the method eliminate unnecessary calculations which result in an overall reduction of the computation time of about 42 percent. The Langmuir isotherm is considered which takes into account the composition-dependent effects on the thermodynamic parameter, mRo.

Feinberg, M. P.↗
NASA NTRS1991

Electromagnetic scattering from arbitrary shaped conducting bodies coated with lossy materials of arbitrary thickness

A simple and efficient numerical technique is presented to solve the electromagnetic scattering problem of coated conducting bodies of arbitrary shape. The surface equivalence principle is used to formulate the problem in terms of a set of coupled integral equations involving equivalent electric and magnetic surface currents which represent boundary fields. The conducting structures and the dielectric materials are modeled by planar triangular patches, and the method of moments is used to solve the integral equations. Numerical results for scattering cross sections are given for various structures and compared with other available data. These results are proved accurate by a number of representative examples.

Rao, Sadasiva M.↗
NASA NTRS1984

Analysis and calculation of macrosegregation in a casting ingot, exhibits C and E

A computer model which describes the solidification of a binary metal alloy in an insulated rectangular mold with a temperature gradient is presented. A numerical technique, applicable to a broad class of moving boundary problems, was implemented therein. The solidification model described is used to calculate the macrosegregation within the solidified casting by coupling the equations for liquid flow in the solid/liquid or mushy zone with the energy equation for heat flow throughout the ingot and thermal convection in the bulk liquid portion. The rate of development of the solid can be automatically calculated by the model. Numerical analysis of such solidification parameters as enthalpy and boundary layer flow is displayed. On-line user interface and software documentation are presented.

Poirier, D. R.↗
NASA NTRS1990

3D Euler analysis of ducted propfan flowfields

A numerical method is presented for predicting the steady inviscid flow through a ducted propfan based on a time-marching solution of the three-dimensional Euler equations. A four-stage multiple-block Runge-Kutta finite volume numerical technique, utilizing implicit residual smoothing and a blended second and fourth difference dissipation, is applied to predict the transonic flowfield about both single-rotation and counter-rotations ducted rotors. Counter-rotation predictions are based on an average-passage system of equations approach. Calculations are performed for both a single sheared H-type grid system and a multiple-block grid system incorporating a C-type grid about the cowl. Numerical results are compared with experimental data for two cases: a low speed ducted propeller and a 1.15 pressure ratio fan stage.

Hall, Edward J.↗
NASA NTRS1974

Performance of finned thermal capacitors

The performance of typical thermal capacitors, both in earth and orbital environments, was investigated. Techniques which were used to make predictions of thermal behavior in a one-g earth environment are outlined. Orbital performance parameters are qualitatively discussed, and those effects expected to be important under zero-g conditions are outlined. A summary of thermal capacitor applications are documentated, along with significant problem areas and current configurations. An experimental program was conducted to determine typical one-g performance, and the physical significance of these data is discussed in detail. Numerical techniques were employed to allow comparison between analytical and experimental data.

Humphries, W. R.↗
NASA NTRS1993

Antenna pattern control using impedance surfaces

A hybrid numerical technique is developed for electrically large pyramidal horn antennas radiating in free space. A stepped-waveguide method is used to analyze the interior surfaces of the horn transition. The Electric Field Integral Equation (EFIE) is employed on the outer surfaces of the pyramidal horn including the radiating aperture. Meanwhile, the Magnetic Field Integral Equation (MFIE) is used on the aperture to relate the aperture fields and those in the horn transition. The resultant hybrid field integral equation (HFIE) is solved numerically by the method of moments. This formulation is both accurate and numerically stable so that high-gain microwave pyramidal horns can be analyzed rigorously. Far-field radiation patterns, both computed and measured, are presented for three electrically-large x-band horn antennas. The comparisons demonstrate that this method is accurate enough to predict the fine pattern structure at wide angles and in the back region. Computed far-field patterns and aperture field distribution of two smaller x-band horns are also presented along with a discussion on the validity of the approximate aperture field distributions routinely used in the analysis and design of pyramidal horns.

Balanis, Constantine A.↗
NASA NTRSDate not supplied

Convolution-Based Numerical Solutions of Transient Temperature Fields during Powder Bed Fusion Additive Manufacturing: Theory, Accuracy, and Computational Cost

Powder bed fusion (PBF) additive manufacturing has found numerous applications in the aerospace domain. However, components fabricated via PBF have a complex time-temperature history that significantly impacts subsequent mechanical performance. This study examines convolution-based numerical solutions of transient temperature fields that support simulations involving arbitrary beam shapes and paths during PBF. The convolutional approach is verified through comparisons with analytical solutions of the temperature field. The computational speed and accuracy of the method are assessed through comparisons with other explicit and implicit numerical techniques. In addition, the straightforward translation of the approach from a CPU to a GPU implementation and the resultant performance improvement are presented. The role of the technique in predicting microstructure evolution during PBF (for a greater process-structure-property-performance framework) is also demonstrated. This work supports the development of computational materials methods for understanding and controlling the time-temperature history during PBF.

powder bed fusion↗
NASA NTRS1974

An improved solution of the simplified two-phase plume

An exact solution of the two-phase, gas-particle, solid propellant rocket exhaust plume flowfield is presented. The good agreement between the exact solution and the experimental data proves that the accuracy of the present technique is acceptable for, at least, preliminary design calculations. A good comparison between the plume impingement effects predicted by the exact solution and those by the best computer code available has also been found for various propellant formulations and rocket sizes. This shows that the present exact solution does describe the basic characteristics of the two-phase plume. In case of a parametric analysis for obtaining an optimum design, the present exact solution will be a more convenient tool than any one of the numerical techniques especially when many parameters are involved in the optimization study. Several typical results are shown herein and their physical significance is discussed.

Wang, S. Y.↗
NASA NTRS1999

Monte Carlo Simulation of Alloy Design Techniques: Fracture and Welding Studied Using the BFS Method for Alloys

Large-scale simulations of dynamic processes at the atomic level have developed into one of the main areas of work in computational materials science. Until recently, severe computational restrictions, as well as the lack of accurate methods for calculating the energetics, resulted in slower growth in the area than that required by current alloy design programs. The Computational Materials Group at the NASA Lewis Research Center is devoted to the development of powerful, accurate, economical tools to aid in alloy design. These include the BFS (Bozzolo, Ferrante, and Smith) method for alloys (ref. 1) and the development of dedicated software for large-scale simulations based on Monte Carlo- Metropolis numerical techniques, as well as state-of-the-art visualization methods. Our previous effort linking theoretical and computational modeling resulted in the successful prediction of the microstructure of a five-element intermetallic alloy, in excellent agreement with experimental results (refs. 2 and 3). This effort also produced a complete description of the role of alloying additions in intermetallic binary, ternary, and higher order alloys (ref. 4).

Bozzolo, Guillermo H.↗
DOE OSTI2021

Intertwined States at Finite Temperatures in the Hubbard Model

We report significant advances in numerical techniques have enabled recent breakthroughs in the study of various properties of the Hubbard model - a seemingly simple, yet complex model of correlated electrons that has been a focus of study for more than half a century. In particular, it captures the essence of strong correlations, and is believed to possess various emergent, low energy states and collective excitations characteristic of cuprate high-temperature superconducting materials. While a thorough review of all activity is not possible here, we have focused the discussion on our recent work using unbiased, numerically exact, “brute force”, finite temperature quantum Monte Carlo methods. Our various studies reveal a rich variety of quantum liquid crystal phases, and complementary transport properties, which answer some questions, but certainly raise others concerning “strange metal” behavior and the ultimate fate of quasiparticles in the Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Learning the ground state of a non-stoquastic quantum Hamiltonian in a rugged neural network landscape

Strongly interacting quantum systems described by non-stoquastic Hamiltonians exhibit rich low-temperature physics. Yet, their study poses a formidable challenge, even for state-of-the-art numerical techniques. Here, we investigate systematically the performance of a class of universal variational wave-functions based on artificial neural networks, by considering the frustrated spin- 1/2 1 / 2 J_1-J_2 J 1 − J 2 Heisenberg model on the square lattice. Focusing on neural network architectures without physics-informed input, we argue in favor of using an ansatz consisting of two decoupled real-valued networks, one for the amplitude and the other for the phase of the variational wavefunction. By introducing concrete mitigation strategies against inherent numerical instabilities in the stochastic reconfiguration algorithm we obtain a variational energy comparable to that reported recently with neural networks that incorporate knowledge about the physical system. Through a detailed analysis of the individual components of the algorithm, we conclude that the rugged nature of the energy landscape constitutes the major obstacle in finding a satisfactory approximation to the ground state wavefunction, and prevents learning the correct sign structure. In particular, we show that in the present setup the neural network expressivity and Monte Carlo sampling are not primary limiting factors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

A numerical method for solving the three-dimensional parabolized Navier-Stokes equations

A numerical technique that solves the parabolized form of the Navier-Stokes equations is presented. Such a method makes it possible to obtain very detailed descriptions of the flowfield in a relatively modest CPU time. The present approach is based on a space-marching technique, uses a finite volume discretization and an upwind flux-difference splitting scheme for the evaluation of the inviscid fluxes. Second order accuracy is achieved following the guidelines of the the ENO schemes. The methodology is used to investigate three-dimensional supersonic viscous flows over symmetric corners. Primary and secondary streamwise vortical structures embedded in the boundary layer and originated by the interaction with shock waves are detected and studied. For purpose of validation, results are compared with experimental data extracted from literature. The agreement is found to be satisfactory. In conclusion, the numerical method proposed seems to be promising as it permits, at a reasonable computational expense, investigation of complex three-dimensional flowfields in great detail.

Dambrosio, Domenic↗
NASA NTRS1985

The method of complex characteristics for transonic airfoil design, with an application to compressors

The use of mathematical models to study physical problems of current interest to aeronautical engineers has been made possible by the development of numerical techniques to compute solutions of the differential equations of transonic aerodynamics. These advances have encouraged the improvement of supercritical wing technology. A method to determined steady, shockless flow of an inviscid, compressible fluid past a cascade of airfoils in the (x,y)-plane is considered, taking into account also the case of an isolated airfoil. The method of complex characteristics solves the equations in the hodograph plane by extending all variables into the complex domain, where the notion of type is no longer significant. Attention is given to the mathematical background, the method of complex characteristics, and numerical calculations.

Bledsoe, M.↗
NASA NTRS1990

Development of a linearized unsteady aerodynamic analysis for cascade gust response predictions

A method for predicting the unsteady aerodynamic response of a cascade of airfoils to entropic, vortical, and acoustic gust excitations is being developed. Here, the unsteady flow is regarded as a small perturbation of a nonuniform isentropic and irrotational steady background flow. A splitting technique is used to decompose the linearized unsteady velocity into rotational and irrotational parts leading to equations for the complex amplitudes of the linearized unsteady entropy, rotational velocity, and velocity potential that are coupled only sequentially. The entropic and rotational velocity fluctuations are described by transport equations for which closed-form solutions in terms of the mean-flow drift and stream functions can be determined. The potential fluctuation is described by an inhomogeneous convected wave equation in which the source term depends on the rotational velocity field, and is determined using finite-difference procedures. The analytical and numerical techniques used to determine the linearized unsteady flow are outlined. Results are presented to indicate the status of the solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described herein leads to very efficient predictions of cascade unsteady aerodynamic response phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.

Verdon, Joseph M.↗