SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗

NASA NTRS1990

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗

NASA NTRS2016

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗

NASA NTRS1993

Three-dimensional flow analysis inside turbomachinery stages with steady and unsteady Navier-Stokes method

This study presents a numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed in the present study. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field of a subsonic turbine stage and the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid.

Copenhaver, W. W.↗

NASA NTRS1972

Optimal trajectories between Earth and Mars in their true planetary orbits.

The optimal transfers from Earth to Mars and from Mars to Earth, considering the actual planetary orbits, are presented as functions of the corresponding idealized Hohmann transfers. The numerically exact two-impulse optimal trajectories are given in graphical form for all possible Hohmann windows. The two-impulse transfers which are absolute optimals and those for which a third impulse provides the absolute optimal are delineated. These data are designed to provide all the information necessary for quick orbit calculations for preliminary Martian mission analysis. In this form, they are as easy to use as the standard Hohmann transfer approximations and provide much greater accuracy.

Gravier, J. P.↗

NASA NTRS1975

Jupiter orbiter lifetime: The hazard of Galilean satellite collision

The four Galilean satellites of Jupiter present a long-term collision hazard to an uncontrolled orbiting spacecraft that repeatedly enters the spatial region occupied by the satellites. Satellite close encounters and the likelihood of collision over a wide range of initial orbit conditions were analyzed. The effect of orbit inclination was of key interest. The scope of the analysis was restricted to orbital dynamic considerations alone, i.e. the question of biological contamination given the event of collision was not considered. A quarantine or orbiter lifetime of 50 years was assumed. This time period begins at spacecraft shutdown following completion of the mission objectives. A numerical approach was adopted wherein each initial orbit is propagated for 50 years, and satellite closest encounter distances recorded on every revolution. The computer program includes approximations of the three major perturbation effects on the long-term motion of the orbiter: (1) Jupiter oblateness, (2) solar gravity, and (3) satellite gravity.

Friedlander, A. L.↗

NASA NTRS1977

Digital color analysis of color-ratio composite LANDSAT scenes

A method is presented that can be used to calculate approximate Munsell coordinates of the colors produced by making a color composite from three registered images. Applied to the LANDSAT MSS data of the Goldfield, Nevada, area, this method permits precise and quantitative definition of the limonitic areas originally observed in a LANDSAT color ratio composite. In addition, areas of transported limonite can be discriminated from the limonite in the hydrothermally altered areas of the Goldfield mining district. From the analysis, the numerical distinction between limonitic and nonlimonitic ground is generally less than 3% using the LANDSAT bands and as much as 8% in ratios of LANDSAT MSS bands.

Raines, G. L.↗

NASA NTRS1985

Electrophoretic mobility of cells in a vertical Ficoll gradiant

The upward migration of living cells and test particles under the influence of a constant electric field in a low conductivity Ficoll gradient occurs at nearly constant velocity. Viscosity and neutral polymer concentration affect migration rate. Decreasing viscosity speeds up the particle migration, decreasing neutral polymer (Ficoll) concentration, slows particle migration, since electrophoretic mobility increases approximately linearly with neutral polymer concentration. Neutral polymers interact with the cell surface to effectively raise its zeta potential. An analytic function was developed from the known dependence of these physical variables on migration distance; the analysis expresses migration velocity as an explicit function of position in the density gradient. It predicts an almost linear increase in velocity of about 12 to 16% over the working region of the gradient. It was numerically integrated and correctly predicts cell migration distance vs time curves without the use of any fitted parameters. The resulting migration curves follow the expected slowly varying exponential form that closely resembles a straight line. The ability to determine standard electrophoretic mobilities from such curves depends on knowledge of the effect of Ficoll on the zeta potential of the cell type that is separated.

Plank, L. D.↗

NASA NTRS1992

A general numerical model for wave rotor analysis

Wave rotors represent one of the promising technologies for achieving very high core temperatures and pressures in future gas turbine engines. Their operation depends upon unsteady gas dynamics and as such, their analysis is quite difficult. This report describes a numerical model which has been developed to perform such an analysis. Following a brief introduction, a summary of the wave rotor concept is given. The governing equations are then presented, along with a summary of the assumptions used to obtain them. Next, the numerical integration technique is described. This is an explicit finite volume technique based on the method of Roe. The discussion then focuses on the implementation of appropriate boundary conditions. Following this, some results are presented which first compare the numerical approximation to the governing differential equations and then compare the overall model to an actual wave rotor experiment. Finally, some concluding remarks are presented concerning the limitations of the simplifying assumptions and areas where the model may be improved.

Paxson, Daniel W.↗

NASA NTRS1993

A model for the beam-filling effect associated with the microwave retrieval of rain

Estimating rain rate from environmental microwave emissions is hampered by several difficulties. One of these difficulties is known as the beam-filling effect. Beam filling is the systematic error introduced when the microwave radiometer's field of view is not filled with uniform rain. Beam filling can have dramatic effects on rain-rate estimation, causing rain rates to be underestimated by as much as a factor of 2. The present study derives an approximate expression for beam filling that provides, in principle, a way to estimate this effect. In addition, this study deals only with single-channel microwave rain estimation over the ocean. The final results reveal that beam filling is essentially determined by the freezing level, the mean fraction of the footprint raining, and the footprint-averaged rain rate. Also, the numerical results appear to agree with other empirical studies. Furthermore, the analysis brings to light an interesting connection with rain threshold techniques for estimating area-averaged rain rates.

Graves, Charles E.↗

NASA NTRS1987

A linear MHD instability analysis of solar mass ejections with gravitation

The linear MHD instability of a cylindrical plasma is used to investigate the origin of solar mass ejections, and the dispersion relation is solved numerically. The initial plasma-flow velocity is found to have a significant effect on the instability criteria and growth rate, and the instability growth-rate is shown to be larger in cases where plasma flow exists, relative to the static case. Results suggest that the plasma column may break into small pieces. Assuming a thin-tube approximation, gravity is found to have little effect on the instability of quasi-horizontal ejection, but to have considerable effect on the vertical ejection. In considering the gravitational force, an exact analytical solution is found for the vertical case, while asymptotic solutions are given for the horizontal and oblique cases.

Song, M. T.↗

NASA NTRS1974

Parameter estimation for terrain modeling from gradient data

A method is developed for modeling terrain surfaces for use on an unmanned Martian roving vehicle. The modeling procedure employs a two-step process which uses gradient as well as height data in order to improve the accuracy of the model's gradient. Least square approximation is used in order to stochastically determine the parameters which describe the modeled surface. A complete error analysis of the modeling procedure is included which determines the effect of instrumental measurement errors on the model's accuracy. Computer simulation is used as a means of testing the entire modeling process which includes the acquisition of data points, the two-step modeling process and the error analysis. Finally, to illustrate the procedure, a numerical example is included.

Dangelo, K. R.↗

NASA NTRS1985

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 3rd, California State University, Long Beach, CA, January 21-24, 1985, Proceedings

Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.

Source record↗

NASA NTRS1990

Domain decomposition methods in aerodynamics

Compressible Euler equations are solved for two-dimensional problems by a preconditioned conjugate gradient-like technique. An approximate Riemann solver is used to compute the numerical fluxes to second order accuracy in space. Two ways to achieve parallelism are tested, one which makes use of parallelism inherent in triangular solves and the other which employs domain decomposition techniques. The vectorization/parallelism in triangular solves is realized by the use of a recording technique called wavefront ordering. This process involves the interpretation of the triangular matrix as a directed graph and the analysis of the data dependencies. It is noted that the factorization can also be done in parallel with the wave front ordering. The performances of two ways of partitioning the domain, strips and slabs, are compared. Results on Cray YMP are reported for an inviscid transonic test case. The performances of linear algebra kernels are also reported.

Venkatakrishnan, V.↗

NASA NTRS1992

Numerical simulation and analysis of a prefrontal squall line. I - Observations and basic simulation results. II - Propagation of the squall line as an internal gravity wave

Results of an observational and numerical study of the squall line that occurred on June 17-18, 1978, are presented. This squall line was initially triggered by the strong surface convergence along a cold front and stretched from Illinois to the Texas Panhandle. The squall line was aligned with the surface front during its initial development, but then propagated faster than the front, resulting in a separation of approximately 200 km by 0300 UTC and 300-400 km by 0600 UTC. The Colorado State University Regional Atmospheric Modeling System is used to model the squall-line development and propagation. Results are described from several experiments that tested the sensitivity to the use of the Kuo-type cumulus parameterization scheme and grid-scale microphysical processes. The movement of the squall line in the model is shown to be due to the propagation of a deep tropospheric internal gravity wave in a CISK-like process. The thermal and dynamic perturbations associated with the hypothesized wave are shown to be consistent with internal gravity wave theory.

Cram, Jennifer M.↗

NASA NTRS1997

Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars

In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.

Mathews, Grant J.↗

NASA NTRSDate not supplied

History and Status of the CDISC Aerodynamic Design Method

This paper will review the development and application of the CDISC aerodynamic design method that began in the early 1980s and continues in active use today. The method uses an iterative, knowledge-based approach that has been coupled with numerous flow solvers utilizing a variety of grid types and levels of flow physics. It has been applied to configurations varying in complexity from 2-D airfoils to full 3-D aircraft at flow conditions ranging from low-speed, high-lift to supersonic cruise. The knowledge-based approach provides a very rapid design process, with designs typically completed in approximately 1-2x the time required for an initial converged analysis of the baseline. The evolution of the method as it tracked the expansive growth of grid and flow solver technology, along with some major applications over 5 decades, are given, followed by a description of the philosophy and process used in design. Finally, an example of its use in a preliminary design environment is given.

CDISC↗

NASA NTRSDate not supplied

History and Status of the CDISC Aerodynamic Design Method

CDISC↗