SEARCH · Engineering Papers

Results for “Numerical approximations”

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 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Protostellar formation in rotating interstellar clouds. IV Nonisothermal collapse

Radiative transfer in the Eddington approximation is included in a multidimensional, self-gravitational, hydrodynamical computer code. Details of the numerical solution and thermodynamic relations are given. Comparison calculations with previous spherically symmetrical models of protostellar collapse are used to validate the basic approach and the artifices which allow the explicit hydrodynamics code to follow the accretion of gas onto a quasi-equilibrium core. A series of axisymmetric models is used to investigate the importance of rotation in collapsing clouds, as the initial amount of angular momentum is lowered, with an emphasis on the possible formation of rings. Rings readily form even in the nonisothermal regime except for very low initial angular momenta; even these clouds may experience ring formation prior to reaching stellar densities. The models imply that other effects (such as gravitational torques or turbulent viscosity) may be necesary to avoid binary formation and thus result in a presolar nebula consistent with the assumptions of either Safronov or Cameron.

Boss, A. P.↗

NASA NTRS1985

A comparative study of the nonuniqueness problem of the potential equation

The nonuniqueness problem occurring at transonic speeds with the conservative potential equation is investigated numerically. The study indicates that the problem is not an inviscid phenomenon, but results from approximate treatment of shock waves inherent in the conservative potential model. A new bound on the limit of validity of the conservative potential model is proposed.

Salas, M. D.↗

NASA NTRS1986

Phase propagation of the solar modulation of galactic cosmic rays

The phase of the 11-yr galactic cosmic ray variation, due to a varying rate of emission of long-lived, propagating regions of enhanced scattering, travels faster than the scattering regions themselves. The radial speed of the 11-yr phase in the quasi-steady force field approximation is exactly twice the speed of the individual episodic decreases. A time-dependent numerical solution for 1-GeV protons at 1 and 30 AU gives a phase speed which is 1.85 times the propagation speed of the individual decreases.

Forman, M. A.↗

NASA NTRS1987

Space Station lubrication considerations

Future activities in space will require the use of large structures and high power availability in order to fully exploit opportunities in Earth and stellar observations, space manufacturing and the development of optimum space transportation vehicles. Although these large systems will have increased capabilities, the associated development costs will be high, and will dictate long life with minimum maintenance. The Space Station provides a concrete example of such a system; it is approximately one hundred meters in major dimensions and has a life requirement of thirty years. Numerous mechanical components will be associated with these systems, a portion of which will be exposed to the space environment. If the long life and low maintenance goals are to be satisfied, lubricants and lubrication concepts will have to be carefully selected. Current lubrication practices are reviewed with the intent of determining acceptability for the long life requirements. The effects of exposure of lubricants and lubricant binders to the space environment are generally discussed. Potential interaction of MoS2 with atomic oxygen, a component of the low Earth orbit environment, appears to be significant.

Leger, Lubert J.↗

NASA NTRS1987

A search for fast planetary-scale instabilities in the equatorial stratosphere

A method proposed by Salby for computing space-time Fourier transforms of asynoptic satellite data has been applied to stratospheric brightness temperature measurements from January and February 1979 in a search for two fast-moving low-wavenumber instabilities predicted in a numerical study by Lynch. The method is capable of producing spectra down to periods of approximately one day with no distortion due to the asynopticity of the data. No spectral peaks corresponding to the predicted modes were found; taking the noise level as the upper limit for their amplitudes at the altitudes observed, an upper limit of approximately 0.6 to 1.2 degrees Kelvin may be set for their peak amplitudes near the mesopause.

Lait, Leslie R.↗

NASA NTRS1988

Estimation of time- and state-dependent delays and other parameters in functional differential equations

A parameter estimation algorithm is developed which can be used to estimate unknown time- or state-dependent delays and other parameters (e.g., initial condition) appearing within a nonlinear nonautonomous functional differential equation. The original infinite dimensional differential equation is approximated using linear splines, which are allowed to move with the variable delay. The variable delays are approximated using linear splines as well. The approximation scheme produces a system of ordinary differential equations with nice computational properties. The unknown parameters are estimated within the approximating systems by minimizing a least-squares fit-to-data criterion. Convergence theorems are proved for time-dependent delays and state-dependent delays within two classes, which say essentially that fitting the data by using approximations will, in the limit, provide a fit to the data using the original system. Numerical test examples are presented which illustrate the method for all types of delay.

Murphy, K. A.↗

NASA NTRS1989

Guidance and Control strategies for aerospace vehicles

A neighboring optimal guidance scheme was devised for a nonlinear dynamic system with stochastic inputs and perfect measurements as applicable to fuel optimal control of an aeroassisted orbital transfer vehicle. For the deterministic nonlinear dynamic system describing the atmospheric maneuver, a nominal trajectory was determined. Then, a neighboring, optimal guidance scheme was obtained for open loop and closed loop control configurations. Taking modelling uncertainties into account, a linear, stochastic, neighboring optimal guidance scheme was devised. Finally, the optimal trajectory was approximated as the sum of the deterministic nominal trajectory and the stochastic neighboring optimal solution. Numerical results are presented for a typical vehicle. A fuel-optimal control problem in aeroassisted noncoplanar orbital transfer is also addressed. The equations of motion for the atmospheric maneuver are nonlinear and the optimal (nominal) trajectory and control are obtained. In order to follow the nominal trajectory under actual conditions, a neighboring optimum guidance scheme is designed using linear quadratic regulator theory for onboard real-time implementation. One of the state variables is used as the independent variable in reference to the time. The weighting matrices in the performance index are chosen by a combination of a heuristic method and an optimal modal approach. The necessary feedback control law is obtained in order to minimize the deviations from the nominal conditions.

Hibey, J. L.↗

NASA NTRS1990

Estimation of time- and state-dependent delays and other parameters in functional differential equations

Murphy, K. A.↗

NASA NTRS1992

Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method

A lattice Boltzmann model is presented which gives the complete Navier-Stokes equation and may provide an efficient parallel numerical method for solving various fluid problems. The model uses the single-time relaxation approximation and a particular Maxwell-type distribution. The model eliminates exactly (1) the non-Galilean invariance caused by a density-dependent coefficient in the convection term and (2) a velocity-dependent equation of state.

Chen, Hudong↗

NASA NTRS1993

The dynamics of magnetic flux rings

The evolution of magnetic fields in the presence of turbulent convection is examined using results of numerical simulations of closed magnetic flux tubes embedded in a steady 'ABC' flow field, which approximate some of the important characteristics of a turbulent convecting flow field. Three different evolutionary scenarios were found: expansion to a steady deformed ring; collapse to a compact fat flux ring, separated from the expansion type of behavior by a critical length scale; and, occasionally, evolution toward an advecting, oscillatory state. The work suggests that small-scale flows will not have a strong effect on large-scale, strong fields.

Deluca, E. E.↗

NASA NTRS1950

The effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings

An analysis is presented of the effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings with partial-span constant-percent-chord ailerons extending inboard from the wing tip. The geometric variables considered are aspect ratio, taper ratio, aileron span, and aileron chord. The shape of the wing-torsional-stiffness curve is assumed and the twisting moment is considered to result solely from the pressure distribution caused by aileron deflection, so that the necessity of using a successive-approximation method is avoided. Because of the complexity of the equations resulting from the analysis, numerical calculations from the equations are presented in a series of figures. A computational form is provided to be used in conjunction with these figures so that calculations can be made without reference to the analysis.

Tucker, Warren A↗

NASA NTRS1995

Computational strategies for tire monitoring and analysis

Computational strategies are presented for the modeling and analysis of tires in contact with pavement. A procedure is introduced for simple and accurate determination of tire cross-sectional geometric characteristics from a digitally scanned image. Three new strategies for reducing the computational effort in the finite element solution of tire-pavement contact are also presented. These strategies take advantage of the observation that footprint loads do not usually stimulate a significant tire response away from the pavement contact region. The finite element strategies differ in their level of approximation and required amount of computer resources. The effectiveness of the strategies is demonstrated by numerical examples of frictionless and frictional contact of the space shuttle Orbiter nose-gear tire. Both an in-house research code and a commercial finite element code are used in the numerical studies.

Danielson, Kent T.↗

NASA NTRS1998

Application of the ENS3DAE Euler/Navier-Stokes Aeroelastic method

This paper discusses recent applications of the ENS3DAE computational aeroelasticity method. In particular, it describes aeroelastic and unsteady aerodynamics calculations performed on wings with trailing edge control surfaces. These simulations include the investigation of control reversal for a structurally flexible wing with a deflected control surface, and a static and dynamic analysis of a rigid wing with an oscillating control surface. The two sets of calculations were performed independently on different wings using different grid topologies. The control reversal simulation represents an inviscid Euler static aeroelastic analysis of a thin wing with a rectangular planform. The geometry of this wing makes it suitable for computations using more approximate, inviscid aerodynamics methods. Thus, the results of the present Euler computations are compared with numerical data generated by a validated computational aeroelasticity code which uses a simpler aerodynamic formulation. The second illustrated case involves the simulation of a significantly more complex flowfield and the static and dynamic analyses of this geometry were performed using the viscous Navier-Stokes equation option in ENS3DAE. Results of both the steady and unsteady calculations on this wing are compared with existing experimental data.

Schuster, David M.↗

NASA NTRS1998

The Geometric Spreading of Coronal Plumes and Coronal Holes

The geometric spreading in plumes and in the interplume region in coronal holes is calculated, using analytic and numerical theoretical models, between 1.0 and 5.0 solar radius. We apply a two-scale approximation that permits the rapid local spreading at the base of plumes (f(sub t)) to be evaluated separately from the global spreading (f(sub g)) imposed by coronal hole geometry. We show that f(sub t) can be computed from a potential-field model and f(sub g) can be computed from global magnetohydrodynamic simulations of coronal structure. The approximations are valid when the plasma beta is mail with respect to unity and for a plume separation small with respect to a solar radius.

Suess, S. T.↗

NASA NTRS1999

Solidification Dynamics of Silver Drops in a Free Fall Environment

Silver drops (99.9%, 4, 5, 7, and 9 mm diameter) were levitated, melted, and released to fall through Marshall Space Flight Center's 105m drop tube in helium - 6% hydrogen and pure argon atmospheres. By systematically varying the initial superheat condition of the drop the extent of solidification prior to impact ranged from complete to none during the approximately 4.6s of free fall time. Comparison of the experimental observations is made with numerical solutions to a model of the heat transfer and solidification kinetics associated with cooling of the drop during free fall, particularly with regard to the fraction of liquid transformed. Analysis reveals the relative importance of the initial parameters affecting the cooling and solidification rates within the drop. A discussion of the conditions under which the actual observations deviate from the assumptions used in the model is presented.

Grugel, Richard N.↗

NASA NTRS2005

A First and Second Order Moment Approach to Probabilistic Control Synthesis

This paper presents a robust control design methodology based on the estimation of the first two order moments of the random variables and processes that describe the controlled response. Synthesis is performed by solving an multi-objective optimization problem where stability and performance requirements in time- and frequency domains are integrated. The use of the first two order moments allows for the efficient estimation of the cost function thus for a faster synthesis algorithm. While reliability requirements are taken into account by using bounds to failure probabilities, requirements related to undesirable variability are implemented by quantifying the concentration of the random outcome about a deterministic target. The Hammersley Sequence Sampling and the First- and Second-Moment- Second-Order approximations are used to estimate the moments, whose accuracy and associated computational complexity are compared numerically. Examples using output-feedback and full-state feedback with state estimation are used to demonstrate the ideas proposed.

Crespo, Luis G.↗

NASA NTRS2011

3D Reconstruction of a Rotating Erupting Prominence

A bright prominence associated with a coronal mass ejection (CME) was seen erupting from the Sun on 9 April 2008. This prominence was tracked by both the Solar Terrestrial Relations Observatory (STEREO) EUVI and COR1 telescopes, and was seen to rotate about the line of sight a it erupted; therefore, the event has been nicknamed the "Cartwheel CME." The threads of the prominence in the core of the CME quite clearly indicate the structure of a weakly to moderately twisted flux rope throughout the field of view, up to heliocentric heights of 4 solar radii. Although the STEREO separation was 48 deg, it was possible to match some sharp features in the later part of the eruption as seen in the 304 A line in EUVI and in the H-alpha-sensitive bandpass of COR I by both STEREO Ahead and Behind. These features could then be traced out in three-dimensional space, and reprojected into a view in which the eruption is directed toward the observer. The reconstructed view shows that the alignment of the prominence to the vertical axis rotates as it rises up to a leading-edge height of approximately equal to 2.5 solar radii, and then remains approximately constant. The alignment at 2.5 solar radii differs by about 115 deg from the original filament orientation inferred from H-alpha and EUV data, and the height profile of the rotation, obtained here for the first time, shows that two thirds of the total rotation are reached within approximately equal to 0.5 solar radii above the photosphere. These features are well reproduced by numerical simulations of an unstable moderately twisted flux rope embedded in external flux with a relatively strong shear field component.

Thompson, W. T.↗

NASA NTRS2022

Hybrid data-driven and model-informed online tool wear detection in milling machines

Precision machining tool wear is responsible for low product throughput and quality. Monitoring the tool wear online is vital to prevent degradation in machining quality. However, direct real-time tool wear measurement is not practical. This paper presents residual-based anomaly detection models, combining a hybrid model comprised of a physics-based model and a data-driven model (a decision tree or a neural network) to predict signals of interest (e.g., power or forces) under nominal conditions, followed by Page’s cumulative sum test for detecting tool wear on-line using the computer numerical control machine measurements. The most informative features are ranked using dynamic programming and its approximation variants from real-time measurements and machine settings, such as the width of cut, depth of cut, feed rate and spindle speed, that serve as inputs to the predictive models. The baseline nominal model is incrementally updated with experimental data via a gradient boosted adaptation model to generate the residuals that account for discrepancies between the actual machine data under normal conditions and the baseline nominal model predictions. The hybrid model is validated against 20 Mazak milling machine experimental tests and one Haas run-to-failure experiment. The proposed anomaly detector is applied to synthetic data from simulations of the physics-based model at different operating conditions, measurement noise levels, and tool wear levels, and the methods were able to achieve an overall 92% accuracy in data with 1% noise. The anomaly detection methods based on hybrid model reduced the false alarms of either the data-driven or physical-based models alone, and are found to be capable of good online detection of tool wear.

Online anomaly detection↗