SEARCH · Engineering Papers
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.
A two-step Taylor-Galerkin finite element method with Lapidus' artificial viscosity scheme is applied to several test cases for internal compressible inviscid flow problems. Investigations for the effect of supersonic/subsonic inlet and outlet boundary conditions on computational results are particularly emphasized.
A program for the generation of unstructured grids over complex configurations, VGRID3D, is described. The grid elements (triangles on the surfaces and tetrahedra in the field) are generated starting from the surface boundaries towards the interior of the computational domain using the Advancing Front Method.
The feasibility of changing the bolt shape from circular to elliptical in order to increase the joint strength was investigated. An analytical method using a cosine series to represent the bearing stress was derived, and a boundary collocation method was used to determine the unknown coefficients of this cosine series. Stresses at the hole edge predicted by the analytical method agreed very well with finite element solutions. Failure analyses of joints in two different laminates were performed, and each laminate exhibited a different joint failure mode. Results demonstrate that the joint strength of both laminates can be improved substantially by changing the bolt shape to elliptical. The joint that failed in a bearing mode showed a greater strength increase compared to the joint that failed in a shearing mode.
The dynamics of a drop of a Newtonian liquid that is pendant from or sessile on a solid rod that is forced to undergo time-periodic oscillations along its axis is studied theoretically. The free boundary problem governing the time evolution of the shape of the drop and the flow field inside it is solved by a method of lines using a finite element algorithm incorporating an adaptive mesh. When the forcing amplitude is small, the drop approaches a limit cycle at large times and undergoes steady oscillations thereafter. However, drop breakup is the consequence if the forcing amplitude exceeds a critical value. Over a wide range of amplitudes above this critical value, drop ejection from the rod occurs during the second oscillation period from the commencement of rod motion. Remarkably, the shape of the interface at breakup and the volume of the primary drop formed are insensitive to changes in forcing amplitude. The interface shape at times close to and at breakup is a multi-valued function of distance measured along the rod axis and hence cannot be described by recently popularized one-dimensional approximations. The computations show that drop ejection occurs without the formation of a long neck. Therefore, this method of drop formation holds promise of preventing formation of undesirable satellite droplets.
Abstract Nb is an elemental superconductor with a critical temperature of 9.3 K and is widely used to fabricate superconducting radiofrequency (SRF) cavities for particle accelerators. However, microstructural defects in Nb, such as grain boundaries (GBs) and dislocations, can act as pinning centers for magnetic flux that can degrade SRF cavity performance. Hydrogen contamination is also detrimental to SRF cavity performance due to the formation of normal conducting hydrides during cool down. In this study, disc shaped Nb bi-crystals extracted from high-purity large-grain Nb slices were investigated to study the effects of GBs, hydrogen, and dislocations on superconducting properties. Grain orientation and GB misorientation were measured using Laue x-ray diffraction and electron backscattered diffraction (EBSD) analyses. Cryogenic magneto-optical imaging was used to directly observe magnetic flux penetration below T c = 9.3 K. Damage caused by low temperature precipitation of hydrides and their dissolution upon reheating after cryogenic cycles was examined using electron channeling contrast imaging, and EBSD. The relationships between hydride formation, dislocation content, GBs, cryo-cooling, heat treatment (HT), and flux penetration indicate that both GB character and hydrogen content affect magnetic flux penetration. Such flux penetration could be facilitated by dislocation structures and low angle GBs resulting from hydride precipitation and HT.
A velocity-pressure integrated, mixed interpolation, Galerkin finite element method for the Navier-Stokes equations is presented. In the method, the velocity variables were interpolated using complete quadratic shape functions and the pressure was interpolated using linear shape functions. For the two dimensional case, the pressure is defined on a triangular element which is contained inside the complete biquadratic element for velocity variables; and for the three dimensional case, the pressure is defined on a tetrahedral element which is again contained inside the complete tri-quadratic element. Thus the pressure is discontinuous across the element boundaries. Example problems considered include: a cavity flow for Reynolds number of 400 through 10,000; a laminar backward facing step flow; and a laminar flow in a square duct of strong curvature. The computational results compared favorable with those of the finite difference methods as well as experimental data available. A finite elememt computer program for incompressible, laminar flows is presented.
Viscous flow over discrete or distributed surface roughness has great implications for hypersonic flight due to aerothermodynamic considerations related to laminar-turbulent transition. Current prediction capability is greatly hampered by the limited knowledge base for such flows. To help fill that gap, numerical computations are used to investigate the intricate flow physics involved. An unstructured mesh, compressible Navier-Stokes code based on the space-time conservation element, solution element (CESE) method is used to perform time-accurate Navier-Stokes calculations for two roughness shapes investigated in wind tunnel experiments at NASA Langley Research Center. It was found through 2D parametric study that at subcritical Reynolds numbers of the boundary layers, absolute instability resulting in vortex shedding downstream, is likely to weaken at supersonic free-stream conditions. On the other hand, convective instability may be the dominant mechanism for supersonic boundary layers. Three-dimensional calculations for a rectangular or cylindrical roughness element at post-shock Mach numbers of 4.1 and 6.5 also confirm that no self-sustained vortex generation is present.
Fatigue cracks are a major problem in designing structures subjected to cyclic loading. Cracks frequently occur in structures such as aircraft and spacecraft. The inspection intervals of many aircraft structures are based on crack-propagation lives. Therefore, improved prediction of propagation lives under flight-load conditions (variable-amplitude loading) are needed to provide more realistic design criteria for these structures. The main thrust was to develop a three-dimensional, nonlinear, elastic-plastic, finite element program capable of extending a crack and changing boundary conditions for the model under consideration. The finite-element model is composed of 8-noded (linear-strain) isoparametric elements. In the analysis, the material is assumed to be elastic-perfectly plastic. The cycle stress-strain curve for the material is shown Zienkiewicz's initial-stress method, von Mises's yield criterion, and Drucker's normality condition under small-strain assumptions are used to account for plasticity. The three-dimensional analysis is capable of extending the crack and changing boundary conditions under cyclic loading.
The internal fluid mechanics research program in inlets, ducts, and nozzles is described. The program consists of a balanced effort between the development of computational tools and the conduct of experimental research. The experiments are designed to better understand the fluid flow physics, to develop new or improved flow models, and to provide benchmark quality data sets for validation of the computational methods. The inlet, duct, and nozzle research program is described according to three major classifications of flow phenomena: highly three-dimensional flow fields; shock-boundary layer interactions; and shear layer control. Specific examples of current and future elements of the research program are described for each of these phenomena. In particular, the highly three-dimensional flow field phenomena is highlighted by describing the computational and experiemental research program in transition ducts having a round-to-rectangular area variation. In the case of shock-boundary layer interactions, the specific details of research for normal shock-boundary layer interactions are described. For shear layer control research in vortex generators and the use of aerodynamic excitation for enhancement of the jet mixing process are described. Future research in inlets, ducts, and nozzles will include more emphasis on three-dimensional full Navier-Stokes methods and corresponding experiments designed to concentrate on the appropriate three-dimensional fluid flow physics.
The elasticity problem for a thick plate consisting of two bonded dissimilar homogeneous layers is considered. It is assumed that the plate is infinite, contains a circular hole, and is subjected to axisymmetric external loads. A technique suitable to the formulation of three-dimensional boundary value problems is introduced and applied to the three-dimensional elasticity solution of an infinite plate with a hole. The effect of Poisson's ratio is studied in some detail, and the axisymmetric problem of a laminated plate is solved. The stresses at the interface are calculated and are observed to become unbounded at the hole boundary. An approximate treatment of the singularity is presented and the stress intensity factors are computed. Some of the results are compared with those found from the finite element method. The agreement between the two sets with regard to trend and magnitude is found to be very satisfactory.
An advanced algorithm for simulating the probabilistic distribution of structural responses due to statistical uncertainties in loads, geometry, material properties, and boundary conditions is reported. The method effectively combines an advanced algorithm for calculating probability levels for multivariate problems (fast probability integration) together with a general-purpose finite-element code for stress, vibration, and buckling analysis. Application is made to a space propulsion system turbine blade for which the geometry and material properties are treated as random variables.
An approach to treating initial-boundary value problems by finite volume methods is described, in which the parallel between differential and difference arguments is closely maintained. By using intrinsic geometrical properties of the volume elements, it is possible to describe discrete versions of the div, curl, and grad operators which lead, using summation-by-parts techniques, to familiar energy equations as well as the div curl = 0 and curl grad = 0 identities. For the diffusion equation, these operators describe compact schemes whose convergence is assured by the energy equations and which yield both the potential and the flux vector with second order accuracy. A simplified potential form is especially useful for obtaining numerical results by multigrid and alternating direction implicit (ADI) methods. The treatment of general curvilinear coordinates is shown to result from a specialization of these general results.
New developments in a posteriori error estimates, smart algorithms, and h- and h-p adaptive finite element methods are discussed in the context of two- and three-dimensional compressible and incompressible flow simulations. Applications to rotor-stator interaction, rotorcraft aerodynamics, shock and viscous boundary layer interaction and fluid-structure interaction problems are discussed.
The fatigue damage of a space propulsion system component is computed using probabilistic structural analysis methods. The analysis takes into account the variations in static and dynamic loads, the uncertainty in structural damping, and the scatter in material fatigue resistance. The key elements of the probabilistic approach include: (1) a numerical engine model for describing the global component interface loads consistent with engine balance, (2) models for computing the local mode boundary conditions on the component, (3) a static structural analysis model that captures the strains at a damage critical location as a function of engine performance variables, (4) a finite element model for assessment of the random amplitude stress response due to random base and pressure excitations with uncertain power and correlation length, and (5) advanced first order reliability methods for computing the probabilities associated with the fatigue damage.
The coupling of the combined source vortex distribution of Green's potential flow function with contemporary numerical techniques is shown to provide accurate, efficient, and stable solutions to subsonic inviscid analysis and design problems for multi-element airfoils. The analysis problem is solved by direct calculation of the surface singularity distribution required to satisfy the flow tangency boundary condition. The design or inverse problem is solved by an iteration process. In this process, the geometry and the associated pressure distribution are iterated until the pressure distribution most nearly corresponding to the prescribed design distribution is obtained. Typically, five iteration cycles are required for convergence. A description of the analysis and design method is presented, along with supporting examples.
Energy storage materials undergo significant stresses during charge/discharge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), image-based modeling techniques are used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electrochemical-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in capturing crack propagation due to mesh dependency. Additionally, commonly used damage models, such as the continuous damage model and the cohesive zone model, often have steep tradeoffs between discontinuous field accuracy and computational expense. In this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022.
A global-local finite element analysis procedure is developed based on the fast convergent nature of FEM in displacement. A special scheme is used to utilize the global displacement solution as boundary conditions for local regions of interest. In the local region, a refined mesh is used for further stress analysis. This global-local procedure can be easily programmed in parallel on MIMD multi-processor computers for significant time savings. A Sequent Balance 21000 system is used for demonstrating the parallel programming.
Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.