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 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

Constituent properties, along with volume fraction, have a first order effect on the microscale fields within a composite material and influence the macroscopic response. Therefore, there is a need to assess the significance of stochastic variation in the constituent properties of composites at the higher scales. The effect of variability in the parameters controlling the time-dependent behavior, in a unidirectional SCS-6 SiC fiber-reinforced RBSN matrix composite lamina, on the residual stresses induced during processing is investigated numerically. The generalized method of cells micromechanics theory is utilized to model the ceramic matrix composite lamina using a repeating unit cell. The primary creep phases of the constituents are approximated using a Norton-Bailey, steady state, power law creep model. The effect of residual stresses on the proportional limit stress and strain to failure of the composite is demonstrated. Monte Carlo simulations were conducted using a normal distribution for the power law parameters and the resulting residual stress distributions were predicted.

processing↗

NASA NTRS1993

High-resolution seismic reflection survey at the Manson crater, Iowa

Approximately 17.4 km of high-resolution reflection data were acquired along an east-west radius of the Manson Impact Structure (MIS) to delineate the shallow (upper 300 m) subsurface structural configuration. The geometry of the shallow structure is poorly known due to a 30-90 m thick Pleistocene till cover. The resolution of the new seismic data is roughly 5-10 times that of existing Vibroseis data. Data quality varies rapidly along the line from exceptional to poor, due primarily to velocity variations associated with the geological complexity of the area. Preliminary results indicate subsurface structural blocks previously envisioned to be several hundreds of meters in size are actually an order of magnitude smaller and more complex. A seismogram-by-seismogram analysis is necessary to confidently identify intricate stratigraphic and structural relationships seen on preliminary CDP sections, as numerous faults, diffractions, and complicated reflection patterns create potential pitfalls.

Keiswetter, D. A.↗

NASA NTRS2012

Comparison of Aircraft Models and Integration Schemes for Interval Management in the TRACON

Reusable models of common elements for communication, computation, decision and control in air traffic management are necessary in order to enable simulation, analysis and assurance of emergent properties, such as safety and stability, for a given operational concept. Uncertainties due to faults, such as dropped messages, along with non-linearities and sensor noise are an integral part of these models, and impact emergent system behavior. Flight control algorithms designed using a linearized version of the flight mechanics will exhibit error due to model uncertainty, and may not be stable outside a neighborhood of the given point of linearization. Moreover, the communication mechanism by which the sensed state of an aircraft is fed back to a flight control system (such as an ADS-B message) impacts the overall system behavior; both due to sensor noise as well as dropped messages (vacant samples). Additionally simulation of the flight controller system can exhibit further numerical instability, due to selection of the integration scheme and approximations made in the flight dynamics. We examine the theoretical and numerical stability of a speed controller under the Euler and Runge-Kutta schemes of integration, for the Maintain phase for a Mid-Term (2035-2045) Interval Management (IM) Operational Concept for descent and landing operations. We model uncertainties in communication due to missed ADS-B messages by vacant samples in the integration schemes, and compare the emergent behavior of the system, in terms of stability, via the boundedness of the final system state. Any bound on the errors incurred by these uncertainties will play an essential part in a composable assurance argument required for real-time, flight-deck guidance and control systems,. Thus, we believe that the creation of reusable models, which possess property guarantees, such as safety and stability, is an innovative and essential requirement to assessing the emergent properties of novel airspace concepts of operation.

Neogi, Natasha↗

NASA NTRS2009

On the Importance of Small Ice Crystals in Tropical Anvil Cirrus

In situ measurements of ice crystal concentrations and sizes made with aircraft instrumentation over the past two decades have often indicated the presence of numerous relatively small (< 50 m diameter) crystals in cirrus clouds. Further, these measurements frequently indicate that small crystals account for a large fraction of the extinction in cirrus clouds. The fact that the instruments used to make these measurements, such as the Forward Scattering Spectrometer Probe (FSSP) and the Cloud Aerosol Spectrometer (CAS), ingest ice crystals into the sample volume through inlets has led to suspicion that the indications of numerous small ]crystals could be artifacts of large ]crystal shattering on the instrument inlets. We present new aircraft measurements in anvil cirrus sampled during the Tropical Composition, Cloud, and Climate Coupling (TC4) campaign with the 2 ] Dimensional Stereo (2D ]S) probe, which detects particles as small as 10 m. The 2D ]S has detector "arms" instead of an inlet tube. Since the 2D ]S probe surfaces are much further from the sample volume than is the case for the instruments with inlets, it is expected that 2D ]S will be less susceptible to shattering artifacts. In addition, particle inter ]arrival times are used to identify and remove shattering artifacts that occur even with the 2D ]S probe. The number of shattering artifacts identified by the 2D ]S interarrival time analysis ranges from a negligible contribution to an order of magnitude or more enhancement in apparent ice concentration over the natural ice concentration, depending on the abundance of large crystals and the natural small ]crystal concentration. The 2D ]S measurements in tropical anvil cirrus suggest that natural small ]crystal concentrations are typically one to two orders of magnitude lower than those inferred from CAS. The strong correlation between the CAS/2D ]S ratio of small ]crystal concentrations and large ]crystal concentration suggests that the discrepancy is likely caused by shattering of large crystals on the CAS inlet. We argue that past measurements with CAS in cirrus with large crystals present may contain errors due to crystal shattering, and past conclusions derived from these measurements may need to be revisited. Further, we present correlations between CAS spurious concentration and 2D ]S large ]crystal mass from spatially uniform anvil cirrus sampling periods as an approximate guide for estimating quantitative impact of large ]crystal shattering on CAS concentrations in previous datasets. We use radiative transfer calculations to demonstrate that in the maritime anvil cirrus sampled during TC4, small crystals indicated by 2D ]S contribute relatively little cloud extinction, radiative forcing, or radiative heating in the anvils, regardless of anvil age or vertical location in the clouds. While 2D ]S ice concentrations in fresh anvil cirrus may often exceed 1 cm.3, and are observed to exceed 10 cm.3 in turrets, they are typically ~0.1 cm.3 and rarely exceed 1 cm.3 (<1.4% of the time) in aged anvil cirrus. We hypothesize that isolated occurrences of higher ice concentrations in aged anvil cirrus may be caused by ice nucleation driven by either small ]scale convection or gravity waves. It appears that the numerous small crystals detrained from convective updrafts do not persist in the anvil cirrus sampled during TC ]4.

Jensen, E. J.↗

NASA NTRS2013

Analysis of the Effects of Streamwise Lift Distribution on Sonic Boom Signature

Investigation of sonic boom has been one of the major areas of study in aeronautics due to the benefits a low-boom aircraft has in both civilian and military applications. This work conducts a numerical analysis of the effects of streamwise lift distribution on the shock coalescence characteristics. A simple wing-canard-stabilator body model is used in the numerical simulation. The streamwise lift distribution is varied by fixing the canard at a deflection angle while trimming the aircraft with the wing and the stabilator at the desired lift coefficient. The lift and the pitching moment coefficients are computed using the Missile DATCOM v. 707. The flow field around the wing-canard- stabilator body model is resolved using the OVERFLOW-2 flow solver. Overset/ chimera grid topology is used to simplify the grid generation of various configurations representing different streamwise lift distributions. The numerical simulations are performed without viscosity unless it is required for numerical stability. All configurations are simulated at Mach 1.4, angle-of-attack of 1.50, lift coefficient of 0.05, and pitching moment coefficient of approximately 0. Four streamwise lift distribution configurations were tested.

Yoo, Paul↗

NASA NTRS2015

The Effect of Stochastically Varying Creep Parameters on Residual Stresses in Ceramic Matrix Composites

residual stress↗

NASA NTRS1995

An Analytic Approximation to Very High Specific Impulse and Specific Power Interplanetary Space Mission Analysis

A simple, analytic approximation is derived to calculate trip time and performance for propulsion systems of very high specific impulse (50,000 to 200,000 seconds) and very high specific power (10 to 1000 kW/kg) for human interplanetary space missions. The approach assumed field-free space, constant thrust/constant specific power, and near straight line (radial) trajectories between the planets. Closed form, one dimensional equations of motion for two-burn rendezvous and four-burn round trip missions are derived as a function of specific impulse, specific power, and propellant mass ratio. The equations are coupled to an optimizing parameter that maximizes performance and minimizes trip time. Data generated for hypothetical one-way and round trip human missions to Jupiter were found to be within 1% and 6% accuracy of integrated solutions respectively, verifying that for these systems, credible analysis does not require computationally intensive numerical techniques.

Williams, Craig Hamilton↗

NASA NTRS1998

Design of Neural Networks for Fast Convergence and Accuracy

A novel procedure for the design and training of artificial neural networks, used for rapid and efficient controls and dynamics design and analysis for flexible space systems, has been developed. Artificial neural networks are employed to provide a means of evaluating the impact of design changes rapidly. Specifically, two-layer feedforward neural networks are designed to approximate the functional relationship between the component spacecraft design changes and measures of its performance. A training algorithm, based on statistical sampling theory, is presented, which guarantees that the trained networks provide a designer-specified degree of accuracy in mapping the functional relationship. Within each iteration of this statistical-based algorithm, a sequential design algorithm is used for the design and training of the feedforward network to provide rapid convergence to the network goals. Here, at each sequence a new network is trained to minimize the error of previous network. The design algorithm attempts to avoid the local minima phenomenon that hampers the traditional network training. A numerical example is performed on a spacecraft application in order to demonstrate the feasibility of the proposed approach.

Maghami, Peiman G.↗

NASA NTRSDate not supplied

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗

NASA NTRS2023

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Chris Ostrom↗

NASA NTRS1998

Development of Advanced Methods of Structural and Trajectory Analysis for Transport Aircraft

This paper develops a near-optimal guidance law for generating minimum fuel, time, or cost fixed-range trajectories for supersonic transport aircraft. The approach uses a choice of new state variables along with singular perturbation techniques to time-scale decouple the dynamic equations into multiple equations of single order (second order for the fast dynamics). Application of the maximum principle to each of the decoupled equations, as opposed to application to the original coupled equations, avoids the two point boundary value problem and transforms the problem from one of a functional optimization to one of multiple function optimizations. It is shown that such an approach produces well known aircraft performance results such as minimizing the Brequet factor for minimum fuel consumption and the energy climb path. Furthermore, the new state variables produce a consistent calculation of flight path angle along the trajectory, eliminating one of the deficiencies in the traditional energy state approximation. In addition, jumps in the energy climb path are smoothed out by integration of the original dynamic equations at constant load factor. Numerical results performed for a supersonic transport design show that a pushover dive followed by a pullout at nominal load factors are sufficient maneuvers to smooth the jump.

Ardema, Mark D.↗

NASA NTRS1983

Tibial changes in experimental disuse osteoporosis in the monkey

The mechanical properties and structural changes in the monkey tibia with disuse osteoporosis and during subsequent recovery are investigated. Bone mending stiffness is evaluated in relation to microscopic changes in cortical bone and Norland bone mineral analysis. Restraint in the semireclined position is found to produce regional losses of bone most obviously in the anterior-proximal tibiae. After six months of restraint, the greatest losses of bone mineral in the proximal tibiae range from 23 percent to 31 percent; the largest changes in bone stiffness range from 36 percent to 40 percent. Approximately eight and one-half months of recovery are required to restore the normal bending properties. Even after 15 months of recovery, however, the bone mineral content does not necessarily return to normal levels. Histologically, resorption cavities in cortical bone are seen within one month of restraint; by two and one-half months of restraint there are large resorption cavities subperiosteally, endosteally, and intracortically. After 15 months of recovery, the cortex consists mainly of first-generation haversian systems. After 40 months, the cortex appears normal, with numerous secondary and tertiary generations of haversian systems.

Young, D. R.↗

NASA NTRS2012

Error Estimates of the Ares I Computed Turbulent Ascent Longitudinal Aerodynamic Analysis

Numerical predictions of the longitudinal aerodynamic characteristics for the Ares I class of vehicles, along with the associated error estimate derived from an iterative convergence grid refinement, are presented. Computational results are based on an unstructured grid, Reynolds-averaged Navier-Stokes analysis. The validity of the approach to compute the associated error estimates, derived from a base grid to an extrapolated infinite-size grid, was first demonstrated on a sub-scaled wind tunnel model at representative ascent flow conditions for which the experimental data existed. Such analysis at the transonic flow conditions revealed a maximum deviation of about 23% between the computed longitudinal aerodynamic coefficients with the base grid and the measured data across the entire roll angles. This maximum deviation from the wind tunnel data was associated with the computed normal force coefficient at the transonic flow condition and was reduced to approximately 16% based on the infinite-size grid. However, all the computed aerodynamic coefficients with the base grid at the supersonic flow conditions showed a maximum deviation of only about 8% with that level being improved to approximately 5% for the infinite-size grid. The results and the error estimates based on the established procedure are also presented for the flight flow conditions.

Abdol-Hamid, Khaled S.↗

NASA NTRS1998

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗

NASA NTRS1972

The motion and stability of a dual spin satellite during the momentum wheel spin-up maneuver

The stability of a dual-spin satellite system during the momentum wheel spin-up maneuver is treated both analytically and numerically. The dual-spin system consists of: a slowly rotating or despun main-body; a momentum wheel (or rotor) which is accelerated by a torque motor to change its initial angular velocity relative to the main part to some high terminal value; and a nutation damper. A closed form solution for the case of a symmetrical satellite indicates that when the nutation damper is physically constrained for movement (i.e. by use of a mechanical clamp) the magnitude of the vector sum of the transverse angular velocity components remains bounded during the wheel spin-up under the influence of a constant motor torque. The analysis is extended to consider such effects as: the motion of the nutation damper during spin-up; a non-uniform motor torque; and the effect of a non-symmetrical mass distribution in the main spacecraft and the rotor. An approximate analytical solution using perturbation techniques is developed for the case of a slightly asymmetric main spacecraft.

Bainum, P. M.↗

NASA NTRS1991

Hierarchic models for laminated plates

Structural plates and shells are three-dimensional bodies, one dimension of which happens to be much smaller than the other two. Thus, the quality of a plate or shell model must be judged on the basis of how well its exact solution approximates the corresponding three-dimensional problem. Of course, the exact solution depends not only on the choice of the model but also on the topology, material properties, loading and constraints. The desired degree of approximation depends on the analyst's goals in performing the analysis. For these reasons models have to be chosen adaptively. Hierarchic sequences of models make adaptive selection of the model which is best suited for the purposes of a particular analysis possible. The principles governing the formulation of hierarchic models for laminated plates are presented. The essential features of the hierarchic models described models are: (1) the exact solutions corresponding to the hierarchic sequence of models converge to the exact solution of the corresponding problem of elasticity for a fixed laminate thickness; and (2) the exact solution of each model converges to the same limit as the exact solution of the corresponding problem of elasticity with respect to the laminate thickness approaching zero. The formulation is based on one parameter (beta) which characterizes the hierarchic sequence of models, and a set of constants whose influence was assessed by a numerical sensitivity study. The recommended selection of these constants results in the number of fields increasing by three for each increment in the power of beta. Numerical examples analyzed with the proposed sequence of models are included and good correlation with the reference solutions was found. Results were obtained for laminated strips (plates in cylindrical bending) and for square and rectangular plates with uniform loading and with homogeneous boundary conditions. Cross-ply and angle-ply laminates were evaluated and the results compared with those of MSC/PROBE. Hierarchic models make the computation of any engineering data possible to an arbitrary level of precision within the framework of the theory of elasticity.

Actis, Ricardo Luis↗

NASA NTRS1994

Airfoil optimization by the one-shot method

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (Governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Language multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Kuruvila, G.↗

NASA NTRS1995

Airfoil Design and Optimization by the One-Shot Method

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Lagrange multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Kuruvila, G.↗