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.
Computer program has been written to analyze group of 100-node areas and then provide for summation of any number of 100-node areas to obtain temperature profile. SINDA program options offer user variety of methods for solution of thermal analog modes presented in network format.
Numerical procedures for accurate determination of elastic stress intensity factors for the general two-dimensional crack problem are reviewed. The elastic perfectly plastic state of crack tip deformation is studied by a finite element procedure. Elastic-plastic fields in the immediate vicinity of a crack tip are determined numerically by finite element procedures based on asymptotic studies of crack tip singularities in plastic materials. The small-scale yielding problem is modeled, and expressions for crack tip opening displacement, shear singularity amplitude, and plastic zone extent are derived. A finite element solution to the large-scale yielding of a circumferentially cracked round tension bar is obtained. The three-dimensional aspects of flawed structures and numerical methods of treating them are studied. Ductile fracture mechanisms, in particular crack tip fracture on the microscale, are discussed.
A computer program designed to optimize propeller characteristics was developed by combining two main programs: the first is the optimization program based on the gradient algorithm; the second is based on a propeller blade element theory and uses an aerodynamics subprogram to approximate the lift and drag characteristics of the NACA 16-series airfoil section. To evaluate the propeller program alone (with its aerodynamics subprogram), propeller characteristics were computed and compared to those from wind tunnel investigations conducted on three different NACA propellers. Although the thrust and power coefficients which were computed using the blade element theory were generally higher than the experimental results for two of the three propellers, the corresponding efficiencies showed good agreement for all three propellers. The propeller optimization program was then used to study the NACA 4-(5)(08)-03 propeller at various Mach numbers from 0.175 to 0.60. Improvements in propeller efficiency and thrust were obtained through the use of the propeller optimization program.
A numerical method for computing incompressible turbulent flows is presented. The method is tested by calculating laminar recirculating flows and is applied in conjunction with a modified Kappa-epsilon model to compute the flow over a backward-facing step. In the laminar regime, the computational results are in good agreement with the experimental data. The turbulent flow study shows that the reattachment length is underpredicted by the standard Kappa-epsilon model. The addition of a term to the standard model that accounts for the effects of rotation on turbulent flow improves the results in the recirculation region and increases the computed reattachment length.
A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.
The numerical algorithms used to simulate the advection, diffusion, sedimentation, coagulation, and condensational growth of atmospheric aerosols are described. The model can be used in one, two, or three spatial dimensions. The continuity equation in a generalized horizontal and vertical coordinate system is developed, which allows the model to be quickly adapted to a wide variety of dynamical models of global or regional scale. Algorithms are developed to treat the various physical processes, and the results of simulations are presented, which show the strengths and weaknesses of these algorithms. Although the emphasis is on the modeling of aerosols, the work is also applicable to the simulations of the transport of gases.
The Research Institute for Advanced Computer Science (RIACS) was established at the NASA Ames Research Center in June of 1983. RIACS is privately operated by the Universities Space Research Association (USRA), a consortium of 64 universities with graduate programs in the aerospace sciences, under several Cooperative Agreements with NASA. RIACS's goal is to provide preeminent leadership in basic and applied computer science research as partners in support of NASA's goals and missions. In pursuit of this goal, RIACS contributes to several of the grand challenges in science and engineering facing NASA: flying an airplane inside a computer; determining the chemical properties of materials under hostile conditions in the atmospheres of earth and the planets; sending intelligent machines on unmanned space missions; creating a one-world network that makes all scientific resources, including those in space, accessible to all the world's scientists; providing intelligent computational support to all stages of the process of scientific investigation from problem formulation to results dissemination; and developing accurate global models for climatic behavior throughout the world. In working with these challenges, we seek novel architectures, and novel ways to use them, that exploit the potential of parallel and distributed computation and make possible new functions that are beyond the current reach of computing machines. The investigation includes pattern computers as well as the more familiar numeric and symbolic computers, and it includes networked systems of resources distributed around the world. We believe that successful computer science research is interdisciplinary: it is driven by (and drives) important problems in other disciplines. We believe that research should be guided by a clear long-term vision with planned milestones. And we believe that our environment must foster and exploit innovation. Our activities and accomplishments for the calendar year 1987 and our plans for 1988 are reported.
Accurate numerical non-reflecting boundary conditions are important in all the proposed benchmark problems of the Second Workshop. Recently, a new absorbing boundary condition has been developed using Perfectly Matched Layer (PML) equations for the Euler equations. In this approach, a region with a width of a few grid points is introduced adjacent to the non-reflecting boundaries. In the added region, Perfectly Matched Layer equations are constructed and applied so that the out-going waves are absorbed inside the layer with little reflection to the interior domain. It will be demonstrated in the present paper that the proposed absorbing boundary condition is quite general and versatile, applicable to radiation boundaries as well as inflow and outflow boundaries. It is also easy to implement. The emphasis of the paper will be on the application of the PML absorbing boundary condition to problems in Categories 1, 2, and 3. In Category 1, solutions of problems 1 and 2 are presented. Both problems are solved using a multi-domain polar grid system. Perfectly Matched Layer equations for a circular boundary are constructed and their effectiveness assessed. In Category 2, solutions of problem 2 are presented. Here, in addition to the radiation boundary conditions at the far field in the axisymmetric coordinate system, the inflow boundary condition at the duct inlet is also dealt with using the proposed Perfectly Match Layer equations. At the inlet, a PML domain is introduced in which the incident duct mode is simulated while the waves reflected from the open end of the duct are absorbed at the same time. In Category 3, solutions of all three problems are presented. Again, the PML absorbing boundary condition is used at the inflow boundary so that the incoming vorticity wave is simulated while the outgoing acoustic waves are absorbed with very little numerical reflection. All the problems are solved using central difference schemes for spatial discretizations and the optimized Low-Dissipation and Low-Dispersion Runge-Kutta scheme for the time integration. Issues of numerical accuracy and efficiency are also addressed.
Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation to the exact problem and uncertainties exist. There is a method to approximate the errors in CFD via Richardson's Extrapolation. This method is based off of progressive grid refinement. To estimate the errors, the analyst must interpolate between at least three grids. This paper describes a study to find an appropriate interpolation scheme that can be used in Richardson's extrapolation or other uncertainty method to approximate errors.
Numerical Modeling of Proppant Transport and Coverage in Rock Fractures
Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.
Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.
A numerical scheme suitable for the computation of both the near field acoustic sources and the far field sound produced by turbulent free shear flows utilizing the Navier-Stokes equations is presented. To produce stable numerical schemes in the presence of shear, damping terms must be added to the boundary conditions. The numerical technique and boundary conditions are found to give stable results for computations of spatially evolving mixing layers.
A numerical method for computing the subsonic flow around multi-element airfoil systems was developed, allowing for flow separation at one or more elements. Besides multiple rear separation also sort bubbles on the upper surface and cove bubbles can approximately be taken into account. Also, compressibility effects for pure subsonic flow are approximately accounted for. After presentation the method is applied to several examples and improved in some details. Finally, the present limitations and desirable extensions are discussed.
Capabilities of several computational schemes combined and enhanced. Advanced Rocket Injector-Combustor Code is computer program for numerical simulation of two-phase flows of chemically reacting or nonreacting components. Applied to sprays of fuel in injectors and combustors, flows in heat exchangers, evaporators, and mass-diffusion systems; and two-dimensional or axisymmetric flows containing liquid or solid suspensions. Used in research for development and verification of various physical models.
A numerical method to compute a two-dimensional hypersonic flowfield that is ionized and in thermochemical nonequilibrium has been developed. Such a flowfield is described by coupled time-dependent partial differential equations for the conservation of species mass, mass-average momentum, vibrational energy of each diatomic species, electron energy, and total mass-averaged energy. The steady-state solution to these fully coupled equations is obtained using an implicit Gauss-Seidel line relaxation technique. The computed electron densities in the flowfield compare well with experimental results.
Implicit finite difference schemes are often the preferred numerical schemes in computational fluid dynamics, requiring less stringent stability bounds than the explicit schemes. Each iteration in an implicit scheme involves global data dependencies in the form of second and higher order recurrences. Efficient parallel implementations of such iterative methods are considerably more difficult and non-intuitive. The parallelization of the implicit schemes that are used for solving the Euler and the thin layer Navier-Stokes equations and that require inversions of large linear systems in the form of block tri-diagonal and/or block penta-diagonal matrices is discussed. Three-dimensional cases are emphasized and schemes that minimize the total execution time are presented. Partitioning and scheduling schemes for alleviating the effects of the global data dependencies are described. An analysis of the communication and the computation aspects of these methods is presented. The effect of the boundary conditions on the parallel schemes is also discussed.
A time dependent numerical procedure that calculates the viscous flow in the shock layer around the shuttle nose has been developed. It has been applied to the free stream conditions that may be encountered in the shuttle flight trajectory and found to be efficient and in satisfactory agreement with other steady and unsteady techniques. Two features associated with the present formulation are responsible for its successful applications. The first is that the flow field computation is made within the shock layer; the second is that the mapping of the computational region places proportionately more mesh points in the vicinity of the wall than in the remaining computational region. The numerical results may be used for low altitude flight to determine the edge conditions for a thin boundary layer, to provide initial boundary layer profiles for downstream boundary layer calculations, and to start supersonic computations. For high altitude flight, it may be used to determine the complete flow field in the nose region and also provide a starter solution for subsequent supersonic flow computations.