SEARCH · Engineering Papers

Results for “Perturbative methods”

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,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Dynamic sensitivity analysis of frictional contact/impact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response to variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the nodal displacements, and the Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. Temporal integration is performed by using Newmark method. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Karaoglan, Levent↗

NASA NTRS1996

Effect of a hypergravity environment on cortical bone elasticity in rats

There is considerable interest in determining whether hypergravity can be used as a countermeasure for microgravity-induced bone loss. This study was conducted on 20 immature male rats in order to investigate possible elastic adaptations of cortical bone in rapidly growing rats exposed to chronic hypergravity. Ten rats were continuously centrifuged for 14 days at twice gravitational acceleration (2G) on a 12.75 foot radius centrifuge and 10 rats concurrently acted as stationary controls. The effect of hypergravity on the elastic characteristics of cortical bone was quantified via ultrasonic wave propagation. Propagation velocities of longitudinal and shear waves were measured through cubic cortical specimens from the posterior femoral diaphyses. Density was measured with an Archimedes' technique. The orthotropic elastic properties were calculated and used to compare the difference between groups. Results showed an average increase in both the Young's moduli (Eii, + 2.2%) and shear moduli (Gij, + 4.3%) with a statistically significant increase only in G12 (+15.7%, P = 0.046). The ratio of transverse to axial strain (Poisson's ratio, nuij) demonstrated statistically significant changes in nu12, nu21, nu13, and nu31 (P < 0.05). These findings suggest that although slight elastic changes were incurred via a hypergravity environment, the treatment level or duration in this study do not dramatically perturb the normal elastic behavior of cortical bone and that dramatic biomechanical differences noted in previous studies were due more to structural changes than material elasticity changes. Hypergravity applied post facto to a microgravity environment would offer further illucidation of this method as treatment for a degenerative spaceflight experience.

Non-NASA Center↗

NASA NTRS2017

Oblate-Earth Effects on the Calculation of Ec During Spacecraft Reentry

The bulge in the Earth at its equator has been shown to lead to a clustering of natural decays biased to occur towards the equator and away from the orbit's extreme latitudes. Such clustering must be considered when predicting the Expectation of Casualty (Ec) during a natural decay because of the clustering of the human population in the same lower latitudes. This study expands upon prior work, and formalizes the correction that must be made to the calculation of the average exposed population density as a result of this effect. Although a generic equation can be derived from this work to approximate the effects of gravitational and atmospheric perturbations on a final decay, such an equation averages certain important subtleties in achieving a best fit over all conditions. The authors recommend that direct simulation be used to calculate the true Ec for any specific entry as a more accurate method. A generic equation is provided, represented as a function of ballistic number and inclination of the entering spacecraft over the credible range of ballistic numbers.

Bacon, John B.↗

NASA NTRSDate not supplied

Computational Study of Low-Speed Boundary-Layer Laminar-Turbulent Transition over a 30º Swept Plate in a Channel

Reliable boundary-layer laminar-turbulent transition predictions are critical for the design of aircraft that use natural or controlled laminar flow to achieve low drag and better fuel efficiency. During cruise flight through low turbulence, the three-dimensional boundary layer over swept surfaces can transition prematurely due to stationary crossflow instabilities induced by surface roughness. Parabolized stability equations analyses have been successfully used to explain the different phases of the transition process under well controlled perturbations such as steps or discrete roughness elements. However, methodologies are still lacking to model the boundary layer receptivity to distributed random roughness, and the following non-linear instability stages leading to transition. For engineering applications, several semi-empirical methods based on linear amplification factors or local helicity have found growing adoption, but their reliability is limited by their level of simplification and the sparsity of data they are based on. To address this gap, a new study of the effect of distributed random surface roughness on crossflow transition was started, using a 30◦ swept plate as a test vehicle. The plate, which is mounted in the 2 Foot × 3 Foot Low-Speed Boundary-Layer Channel at the NASA Langley Research Center, has been tested with micro-scale roughness of varied geometry near the leading edge. In order to analyze the flow in the test section and prepare Direct Numerical Simulations with resolved roughness, Reynolds-Averaged Navier-Stokes (RANS) simulations of the flow over the plate installed in the wind tunnel test section were performed using the Langtry-Menter transition model with stationary crossflow extension. This paper presents the results of these simulations, with comparison to the experiment.

ARMD↗

NASA NTRS2024

An Extensible Perturbed Parameter Ensemble for the Community Atmosphere Model Version 6

This paper documents the methodology and preliminary results from a Perturbed Parameter Ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2-CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations are analyzed and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks. The use of various emulators is explored in the multi- dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that using PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2-CAM6 PPE is publicly available, and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g. coupling to the land surface).

Machine learning↗

NASA NTRS1975

Hybrid computer methods for direct functional optimization

Control and trajectory optimization involves the minimization of a performance index (PI) of integral form where some optimal control law exists in a dynamic system. In this paper, a hybrid minicomputer with an adaptive random-search algorithm implements an iterative search for the optimal control. The search assumes that some initial control is randomly perturbed and a fast analog computer generates respective PI from the analog response of the dynamic system. An improved PI informs the digital computer to utilize the perturbed control as a basis for the next iteration; otherwise a new perturbation replaces the old perturbation in the next iteration. The search terminates when no further improvements occur.

Andrews, M.↗

NASA NTRS1979

First results of a study on turbulent boundary layers in oscillating flow with a mean adverse pressure gradient

The development of a turbulent unsteady boundary layer with a mean pressure gradient strong enough to induce separation, in order to complete the extend results obtained for the flat plate configuration is presented. The longitudinal component of the velocity is measured using constant temperature hot wire anemometer. The region where negative velocities exist is investigated with a laser Doppler velocimeter system with BRAGG cells. The boundary layer responds by forced pulsation to the perturbation of potential flow. The unsteady effects observed are very important. The average location of the zero skin friction point moves periodically at the perturbation frequency. Average velocity profiles from different instants in the cycle are compared. The existence of a logarithmic region enables a simple calculation of the maximum phase shift of the velocity in the boundary layer. An attempt of calculation by an integral method of boundary layer development is presented, up to the point where reverse flow starts appearing.

Houdeville, R.↗

NASA NTRS1988

Multiple directed graph large-class multi-spectral processor

Numerical analysis techniques for the interpretation of high-resolution imaging-spectrometer data are described and demonstrated. The method proposed involves the use of (1) a hierarchical classifier with a tree structure generated automatically by a Fisher linear-discriminant-function algorithm and (2) a novel multiple-directed-graph scheme which reduces the local maxima and the number of perturbations required. Results for a 500-class test problem involving simulated imaging-spectrometer data are presented in tables and graphs; 100-percent-correct classification is achieved with an improvement factor of 5.

Casasent, David↗

NASA NTRS2011

Time-Distance Helioseismology Data-Analysis Pipeline for Helioseismic and Magnetic Imager Onboard Solar Dynamics Observatory (SDO-HMI) and Its Initial Results

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.

Sun: helioseismology↗

NASA NTRSDate not supplied

Steady and Unsteady Simulations of Transonic Truss-BracedWing Aircraft for Flight Dynamic Stability Analysis

This paper presents steady and unsteady simulations of the Mach 0.8 and Mach 0.745Transonic Truss-Braced Wing (TTBW) aircrafts using the high-fidelity CFD solver FUN3D for the flight dynamic stability analysis. The steady-state stability derivatives with respect to the angle of attack, angle of sideslip, and airspeed are calculated with perturbations in the angle of attack, angle of sideslip, and Mach number, respectively. A series of unsteady CFD simulations conducted for the TTBW aircrafts in pitch oscillation at various reduced frequencies. The dynamic stability derivatives are estimated using a frequency domain estimation method. The results are then compared to the results obtained from the VSPAERO stability analysis.

CFD↗

NASA NTRS2008

Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio

Gradient approximation methods commonly used in unstructured-grid finite-volume schemes intended for solutions of high Reynolds number flow equations are studied comprehensively. The accuracy of gradients within cells and within faces is evaluated systematically for both node-centered and cell-centered formulations. Computational and analytical evaluations are made on a series of high-aspect-ratio grids with different primal elements, including quadrilateral, triangular, and mixed element grids, with and without random perturbations to the mesh. Both rectangular and cylindrical geometries are considered; the latter serves to study the effects of geometric curvature. The study shows that the accuracy of gradient reconstruction on high-aspect-ratio grids is determined by a combination of the grid and the solution. The contributors to the error are identified and approaches to reduce errors are given, including the addition of higher-order terms in the direction of larger mesh spacing. A parameter GAMMA characterizing accuracy on curved high-aspect-ratio grids is discussed and an approximate-mapped-least-square method using a commonly-available distance function is presented; the method provides accurate gradient reconstruction on general grids. The study is intended to be a reference guide accompanying the construction of accurate and efficient methods for high Reynolds number applications

Thomas, James↗

NASA NTRS1994

Enhanced simulation software for rocket turbopump, turbulent, annular liquid seals

One of the main objectives of this work is to develop a new dynamic analysis for liquid annular seals with arbitrary profile and to analyze a general distorted interstage seal of the space shuttle main engine high pressure oxygen turbopump (SSME-ATD-HPOTP). The dynamic analysis developed is based on a method originally proposed by Nelson and Nguyen. A simpler scheme based on cubic splines is found to be computationally more efficient and has better convergence properties at higher eccentricities. The first order solution of the original analysis is modified by including a more exact solution that takes into account the variation of perturbed variables along the circumference. A new set of equations for dynamic analysis are derived based on this more general model. A unified solution procedure that is valid for both Moody's and Hirs' friction models is presented. Dynamic analysis is developed for three different models: constant properties, variable properties, and thermal effects with variable properties. Arbitrarily varying seal profiles in both axial and circumferential directions are considered. An example case of an elliptical seal with varying degrees of axial curvature is analyzed in detail. A case study based on predicted clearances of an interstage seal of the SSME-ATD-HPOTP is presented. Dynamic coefficients based on external specified load are introduced to analyze seals that support a preload. The other objective of this work is to study the effect of large rotor displacements of SSME-ATD-HPOTP on the dynamics of the annular seal and the resulting transient motion. One task is to identify the magnitude of motion of the rotor about the centered position and establish limits of effectiveness of using current linear models. This task is accomplished by solving the bulk flow model seal governing equations directly for transient seal forces for any given type of motion, including motion with large eccentricities. Based on the above study, an equivalence is established between linearized coefficients based transient motion and the same motion as predicted by the original governing equations. An innovative method is developed to model nonlinearities in an annular seal based on dynamic coefficients computed at various static eccentricities. This method is thoroughly tested for various types of transient motion using bulk flow model results as a benchmark.

Padavala, Satya↗

NASA NTRS1976

Filamentation instability of a large amplitude whistler wave

Experiments performed in a pulsed afterglow plasma column with specified parameters revealed a filamentation instability of a large-amplitude whistler wave launched from antennas which produce a diverging energy flow in the linear regime. The difficult problem of diagnosing local density perturbations in a magnetized plasma in the presence of large-amplitude RF signals and nonuniform anisotropic electron distributions is discussed. Since Langmuir probes are too unreliable under these conditions, a microwave probe is developed based on the principle of the cavity shift method. The temporal and spatial evolution of the duct formation and wave propagation are shown. The role of the observed electron heating in the filamentation process is examined. The observed self-focusing process appears to be the result of the reinforcing interaction between the wave-induced density depression and the density-induced wave refraction. Interesting applications are noted.

Stenzel, R. L.↗

NASA NTRS1979

The rotational motion of an earth orbiting gyroscope according to the Einstein theory of general relativity

The classical Poisson equations of rotational motion are used to study the attitude motions of an earth orbiting, rapidly spinning gyroscope perturbed by the effects of general relativity (Einstein theory). The center of mass of the gyroscope is assumed to move about a rotating oblate earth in an evolving elliptic orbit which includes all first-order oblateness effects produced by the earth. A method of averaging is used to obtain a transformation of variables, for the nonresonance case, which significantly simplifies the Poisson differential equations of motion of the gyroscope. Long-term solutions are obtained by an exact analytical integration of the simplified transformed equations. These solutions may be used to predict both the orientation of the gyroscope and the motion of its rotational angular momentum vector as viewed from its center of mass. The results are valid for all eccentricities and all inclinations not near the critical inclination.

Hoots, F. R.↗

NASA NTRS1988

Mapping electrodynamic features of the high-latitude ionosphere from localized observations - Technique

This paper describes a novel procedure for mapping high-latitude electric fields and currents and their associated magnetic variations, using sets of localized observational data derived from different types of measurements. The technique provides a formalism for incorporating simultaneously such different classes of data as electric fields from radars and satellites, electric currents from radars, and magnetic perturbations at the ground and at satellite heights; the technique also uses available statistical information on the averages and variances of electrodynamic fields. The technique provides a more rigorous way of quantitatively estimating high-latitude electric field and current patterns than other methods and has a capability to quantify the errors in the mapped fields, based on the distribution of available data, their errors, and the statistical variances of the fields. The technique is illustrated by an application to a substorm which was analyzed by Kamide et al. (1982) by an earlier technique.

Richmond, A. D.↗

NASA NTRS1992

A self-consistent field method for galactic dynamics

The present study describes an algorithm for evolving collisionless stellar systems in order to investigate the evolution of systems with density profiles like the R exp 1/4 law, using only a few terms in the expansions. A good fit is obtained for a truncated isothermal distribution, which renders the method appropriate for galaxies with flat rotation curves. Calculations employing N of about 10 exp 6-7 are straightforward on existing supercomputers, making possible simulations having significantly smoother fields than with direct methods such as tree-codes. Orbits are found in a given static or time-dependent gravitational field; the potential, phi(r, t) is revised from the resultant density, rho(r, t). Possible scientific uses of this technique are discussed, including tidal perturbations of dwarf galaxies, the adiabatic growth of central masses in spheroidal galaxies, instabilities in realistic galaxy models, and secular processes in galactic evolution.

Hernquist, Lars↗

NASA NTRS2001

Sound Emission of Rotor Induced Deformations of Generator Casings

The casing of large electrical generators can be deformed slightly by the rotor's magnetic field. The sound emission produced by these periodic deformations, which could possibly exceed guaranteed noise emission limits, is analysed analytically and numerically. From the deformation of the casing, the normal velocity of the generator's surface is computed. Taking into account the corresponding symmetry, an analytical solution for the acoustic pressure outside the generator is round in terms of the Hankel function of second order. The normal velocity or the generator surface provides the required boundary condition for the acoustic pressure and determines the magnitude of pressure oscillations. For the numerical simulation, the nonlinear 2D Euler equations are formulated In a perturbation form for low Mach number Computational Aeroacoustics (CAA). The spatial derivatives are discretized by the classical sixth-order central interior scheme and a third-order boundary scheme. Spurious high frequency oscillations are damped by a characteristic-based artificial compression method (ACM) filter. The time derivatives are approximated by the classical 4th-order Runge-Kutta method. The numerical results are In excellent agreement with the analytical solution.

Polifke, W.↗

NASA NTRS2004

High-Order Shock-Capturing Methods for Modeling Dynamics of the Solar Atmosphere

We use one-dimensional high-order central shock capturing numerical methods to study the response of various model solar atmospheres to forcing at the solar surface. The dynamics of the atmosphere is modeled with the Euler equations in a variable-sized flux tube in the presence of gravity. We study dynamics of the atmosphere suggestive of spicule formation and coronal oscillations. These studies are performed on observationally-derived model atmospheres above the quiet sun and above sunspots. To perform these simulations, we provide a new extension of existing second- and third- order shock-capturing methods to irregular grids. We also solve the problem of numerically maintaining initial hydrostatic balance via the introduction of new variables in the model equations and a careful initialization mechanism. We find several striking results: all model atmospheres respond to a single impulsive perturbation with several strong shock waves consistent with the rebound-shock model. These shock waves lift material and the transition region well into the initial corona, and the sensitivity of this lift to the initial impulse depends non-linearly on the details of the atmosphere model. We also reproduce an observed 3-minute coronal oscillation above sunspots compared to 5-minute oscillations above the quiet sun.

Bryson, Steve↗