SEARCH · Engineering Papers

Results for “First Order Second Moment”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Part A: Nonprincipal-plane scattering from flat plates: Second-order and corner diffractions

Two models of a flat plate for nonprincipal-plane scattering are explored. The first is a revised version of the Physical Optics/Physical Theory of Diffraction (PO/PTD) model with second-order PTD equivalent currents included to account for second-order interactions among the plate edges. The second model uses a heurisitcally derived corner diffraction coefficient to account for the corner scattering mechanism. The patterns obtained using the newer models were compared to the data of previously reported models, the Moment Method (MM), and experimental results. Near normal incidence, all the models agreed; however, near grazing incidence a need for higher-order and corner diffraction mechanisms was noted. In many instances the second-order and corner-scattered fields which were formulated improved the results.

Balanis, Constantine A.↗

NASA NTRS1993

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗

NASA NTRS2007

NASA Langley's Approach to the Sandia's Structural Dynamics Challenge Problem

The objective of this challenge is to develop a data-based probabilistic model of uncertainty to predict the behavior of subsystems (payloads) by themselves and while coupled to a primary (target) system. Although this type of analysis is routinely performed and representative of issues faced in real-world system design and integration, there are still several key technical challenges that must be addressed when analyzing uncertain interconnected systems. For example, one key technical challenge is related to the fact that there is limited data on target configurations. Moreover, it is typical to have multiple data sets from experiments conducted at the subsystem level, but often samples sizes are not sufficient to compute high confidence statistics. In this challenge problem additional constraints are placed as ground rules for the participants. One such rule is that mathematical models of the subsystem are limited to linear approximations of the nonlinear physics of the problem at hand. Also, participants are constrained to use these models and the multiple data sets to make predictions about the target system response under completely different input conditions. Our approach involved initially the screening of several different methods. Three of the ones considered are presented herein. The first one is based on the transformation of the modal data to an orthogonal space where the mean and covariance of the data are matched by the model. The other two approaches worked solutions in physical space where the uncertain parameter set is made of masses, stiffnesses and damping coefficients; one matches confidence intervals of low order moments of the statistics via optimization while the second one uses a Kernel density estimation approach. The paper will touch on all the approaches, lessons learned, validation 1 metrics and their comparison, data quantity restriction, and assumptions/limitations of each approach. Keywords: Probabilistic modeling, model validation, uncertainty quantification, kernel density

Horta, Lucas G.↗

NASA NTRS1987

A higher-order theory for geometrically nonlinear analysis of composite laminates

A third-order shear deformation theory of laminated composite plates and shells is developed, the Navier solutions are derived, and its finite element models are developed. The theory allows parabolic description of the transverse shear stresses, and therefore the shear correction factors of the usual shear deformation theory are not required in the present theory. The theory also accounts for the von Karman nonlinear strains. Closed-form solutions of the theory for rectangular cross-ply and angle-ply plates and cross-ply shells are developed. The finite element model is based on independent approximations of the displacements and bending moments (i.e., mixed finite element model), and therefore, only C sup o -approximation is required. The finite element model is used to analyze cross-ply and angle-ply laminated plates and shells for bending and natural vibration. Many of the numerical results presented here should serve as references for future investigations. Three major conclusions resulted from the research: First, for thick laminates, shear deformation theories predict deflections, stresses and vibration frequencies significantly different from those predicted by classical theories. Second, even for thin laminates, shear deformation effects are significant in dynamic and geometrically nonlinear analyses. Third, the present third-order theory is more accurate compared to the classical and firt-order theories in predicting static and dynamic response of laminated plates and shells made of high-modulus composite materials.

Reddy, J. N.↗

NASA NTRS1980

An analysis of short pulse and dual frequency radar techniques for measuring ocean wave spectra from satellites

Scanning beam microwave radars were used to measure ocean wave directional spectra from satellites. In principle, surface wave spectral resolution in wave number can be obtained using either short pulse (SP) or dual frequency (DF) techniques; in either case, directional resolution obtains naturally as a consequence of a Bragg-like wave front matching. A four frequency moment characterization of backscatter from the near vertical using physical optics in the high frequency limit was applied to an analysis of the SP and DF measurement techniques. The intrinsic electromagnetic modulation spectrum was to the first order in wave steepness proportional to the large wave directional slope spectrum. Harmonic distortion was small and was a minimum near 10 deg incidence. NonGaussian wave statistics can have an effect comparable to that in the second order of scattering from a normally distributed sea surface. The SP technique is superior to the DF technique in terms of measurement signal to noise ratio and contrast ratio.

Jackson, F. C.↗

NASA NTRS1995

Modeling of nonequilibrium space plasma flows

Godunov-type numerical solution of the 20 moment plasma transport equations. One of the centerpieces of our proposal was the development of a higher order Godunov-type numerical scheme to solve the gyration dominated 20 moment transport equations. In the first step we explored some fundamental analytic properties of the 20 moment transport equations for a low b plasma, including the eigenvectors and eigenvalues of propagating disturbances. The eigenvalues correspond to wave speeds, while the eigenvectors characterize the transported physical quantities. In this paper we also explored the physically meaningful parameter range of the normalized heat flow components. In the second step a new Godunov scheme type numerical method was developed to solve the coupled set of 20 moment transport equations for a quasineutral single-ion plasma. The numerical method and the first results were presented at several national and international meetings and a paper describing the method has been published in the Journal of Computational Physics. To our knowledge this is the first numerical method which is capable of producing stable time-dependent solutions to the full 20 (or 16) moment set of transport equations, including the full heat flow equation. Previous attempts resulted in unstable (oscillating) solutions of the heat flow equations. Our group invested over two man-years into the development and implementation of the new method. The present model solves the 20 moment transport equations for an ion species and thermal electrons in 8 domain extending from a collision dominated to a collisionless region (200 km to 12,000 km). This model has been applied to study O+ acceleration due to Joule heating in the lower ionosphere.

Gombosi, Tamas↗

NASA NTRS2016

Nano-ADEPT Aeroloads Wind Tunnel Test

A wind tunnel test of the Adaptable Deployable Entry and Placement Technology (ADEPT) was conducted in April 2015 at the US Army's 7 by10 Foot Wind Tunnel located at NASA Ames Research Center. Key geometric features of the fabric test article were a 0.7 meter deployed base diameter, a 70 degree half-angle forebody cone angle, eight ribs, and a nose-to-base radius ratio of 0.7. The primary objective of this wind tunnel test was to obtain static deflected shape and pressure distributions while varying pretension at dynamic pressures and angles of attack relevant to entry conditions at Earth, Mars, and Venus. Other objectives included obtaining aerodynamic force and moment data and determining the presence and magnitude of any dynamic aeroelastic behavior (buzz/flutter) in the fabric trailing edge. All instrumentation systems worked as planned and a rich data set was obtained. This paper describes the test articles, instrumentation systems, data products, and test results. Four notable conclusions are drawn. First, test data support adopting a pre-tension lower bound of 10 foot pounds per inch for Nano-ADEPT mission applications in order to minimize the impact of static deflection. Second, test results indicate that the fabric conditioning process needs to be reevaluated. Third, no flutter/buzz of the fabric was observed for any test condition and should also not occur at hypersonic speeds. Fourth, translating one of the gores caused ADEPT to generate lift without the need for a center of gravity offset. At hypersonic speeds, the lift generated by actuating ADEPT gores could be used for vehicle control.

system engineering↗

NASA NTRS1993

Porosity and the ecology of icy satellites

The case for a significant role for porosity in the structure and evolution of icy bodies in the Solar System has been difficult to establish. We present a relevant new data set and a series of structure models including a mechanical compression, not thermal creep, model for porosity that accounts satisfactorily for observed densities, moments of inertia, geologic activity, and sizes of tectonic features on icy satellites. Several types of observational data sets have been used to infer significant porosity, but until recently, alternative explanations have been preferred. Our first area of concern is the occurrence of cryovolcanism as a function of satellite radius; simple radiogenic heating models of icy satellites suggest minimum radii for melting and surface cryovolcanism to be 400 to 500 km, yet inferred melt deposits are seen on satellites half that size. One possible explanation is a deep, low conductivity regolith which lowers conductivity and raises internal temperatures, but other possibilities include tidal heating or crustal compositions of low conductivity. Our second area of concern is the occurrence and magnitude of tectonic strain; tectonic structures have been seen on icy satellites as small as Mimas and Proteus. The structures are almost exclusively extensional, with only a few possible compression Al features, and inferred global strains are on the order of 1 percent expansion. Expansions of this order in small bodies like Mimas and prevention of late compressional tectonics due to formation of ice mantles in larger bodies like Rhea are attained only in structure models including low-conductivity, and thus possibly high porosity, crusts. Thirdly, inferred moments of inertia less than 0.4 in Mimas and Tethys can be explained by high-porosity crusts, but also by differentiation of a high density core. Finally, the relatively low densities of smaller satellites like Mimas and Miranda relative to larger neighbors can be explained by deep porosity, but also by bulk compositional differences. Recent work has strengthened the case for significant porosity. Halley's nucleus was found to have a density near 0.6 g/cu cm, Janus and Epimethus were proposed to have densities near 0.7 g/cu cm, densities almost certainly due to high porosity. The irregular-spherical shape transition of icy satellites was quantitatively explained by low conductivity regoliths. A creative structure/thermal history model for Mimas simultaneously accounts quantitatively for Mimas' low density and moment of inertia by invoking initial high-porosity and subsequent compaction in the deep interior by thermal creep. The main problem with this promising model is that approximately 7 percent predicts a reduction in Mimas' radius, implying significant compressional failure and prevention of extensional tectonics, in contradiction to the observed extensional features and inferred 1 percent expansion in radius.

Croft, Steven K.↗

NASA NTRS1992

A unified perturbation expansion for surface scattering

Starting with the extinction theorem, a perturbation expansion which, to first and second orders, converges over a wider domain than the small perturbation expansion and the momentum transfer expansion is presented. It is shown that, in the appropriate limits, both of these theories, as well as the two-scale expansion, are recovered. There is no adjustable parameter, such as a spectral split, in the theory. This theory is applied to random rough surfaces and derive analytic expressions for the coherent field and the bistatic cross section. Finally, a numerical test of the theory against method of moments results for Gaussian random rough surfaces with a power law spectrum is given. These results show that the expansion is ramarkably accurate over a large range of surface heights and slopes for both horizontal and vertical polarization.

Rodriguez, Ernesto↗

NASA NTRS1993

Plasma heating by collisionless magnetic reconnection - Analysis and computation

Analytic and numerical results on particle acceleration in 2D collisionless magnetic reconnection are presented. The particles are followed until they reach an outgoing flux surface at the same distance from the origin as the starting surface. The magnetic moment is not conserved for particles passing through the unmagnetized region around the X line at the origin. Other particles cross the separatrix without passing near the X line. The magnetic moment of the first class of outgoing particles is randomized, whereas it can be considered for the second class. The analytic model is based upon the observation of the final kinetic energy as a function of the initial conditions. Analytic results are shown to predict a Maxwellian tail for the distribution function in the perpendicular kinetic energy, with this energy much greater than the parallel kinetic energy. Numerical results showing that the predicted tail temperature agrees with the numerically computed temperature to within 10 percent over 4 orders of magnitude in the electric field are presented.

Moses, R. W.↗

NASA NTRS2007

Multi-Layer Arctic Mixed-Phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity Experiments

A cloud-resolving model (CRM) is used to simulate the multiple-layer mixed-phase stratiform (MPS) clouds that occurred during a three-and-a-half day subperiod of the Department of Energy-Atmospheric Radiation Measurement Program s Mixed-Phase Arctic Cloud Experiment (M-PACE). The CRM is implemented with an advanced two-moment microphysics scheme, a state-of-the-art radiative transfer scheme, and a complicated third-order turbulence closure. Concurrent meteorological, aerosol, and ice nucleus measurements are used to initialize the CRM. The CRM is prescribed by time-varying large-scale advective tendencies of temperature and moisture and surface turbulent fluxes of sensible and latent heat. The CRM reproduces the occurrences of the single- and double-layer MPS clouds as revealed by the M-PACE observations. However, the simulated first cloud layer is lower and the second cloud layer thicker compared to observations. The magnitude of the simulated liquid water path agrees with that observed, but its temporal variation is more pronounced than that observed. As in an earlier study of single-layer cloud, the CRM also captures the major characteristics in the vertical distributions and temporal variations of liquid water content (LWC), total ice water content (IWC), droplet number concentration and ice crystal number concentration (nis) as suggested by the aircraft observations. However, the simulated mean values differ significantly from the observed. The magnitude of nis is especially underestimated by one order of magnitude. Sensitivity experiments suggest that the lower cloud layer is closely related to the surface fluxes of sensible and latent heat; the upper cloud layer is probably initialized by the large-scale advective cooling/moistening and maintained through the strong longwave (LW) radiative cooling near the cloud top which enhances the dynamical circulation; artificially turning off all ice-phase microphysical processes results in an increase in LWP by a factor of 3 due to interactions between the excessive LW radiative cooling and extra cloud water; heating caused by phase change of hydrometeors could affect the LWC and cloud top height by partially canceling out the LW radiative cooling. It is further shown that the resolved dynamical circulation appears to contribute more greatly to the evolution of the MPS cloud layers than the parameterized subgrid-scale circulation.

Luo, Yali↗

NASA NTRS1997

The Benchmark Active Controls Technology Model Aerodynamic Data

The Benchmark Active Controls Technology (BACT) model is a part of the Benchmark Models Program (BMP). The BMP is a NASA Langley Research Center program that includes a series of models which were used to study different aeroelastic phenomena and to validate computational fluid dynamics codes. The primary objective of BACT testing was to obtain steady and unsteady loads, accelerations, and aerodynamic pressures due to control surface activity in order to calibrate unsteady CFD codes and active control design tools. Three wind-tunnel tests in the Transonic Dynamics Tunnel (TDT) have been completed. The first and parts of the second and third tests focused on collecting open-loop data to define the model's aeroservoelastic characteristics, including the flutter boundary across the Mach range. It is this data that is being presented in this paper. An extensive database of over 3000 data sets was obtained. This database includes steady and unsteady control surface effectiveness data, including pressure distributions, control surface hinge moments, and overall model loads due to deflections of a trailing edge control surface and upper and lower surface

Scott, Robert C.↗

NASA NTRS1938

Compressible flow about symmetrical Joukowski profiles

The method of Poggi is employed for the determination of the effects of compressibility upon the flow past an obstacle. A general expression for the velocity increment due to compressibility is obtained. The general result holds whatever the shape of the obstacle; but, in order to obtain the complete solution, it is necessary to know a certain Fourier expansion of the square of the velocity of flow past the obstacle. An application is made to the case flow of a symmetrical Joukowski profile with a sharp trailing edge, fixed in a stream of an arbitrary angle of attack and with the circulation determined by the Kutta condition. The results are obtained in a closed form and are exact insofar as the second approximation to the compressible flow is concerned, the first approximation being the result for the corresponding incompressible flow. Formulas for lift and moment analogous to the Blasius formulas in incompressible flow are developed and are applied to thin symmetrical Joukowski profiles for small angles of attack.

Kaplan, Carl↗

NASA NTRS1997

Compton Scattering by Static and Moving Media: The Transfer Equation and its Moments - Part 1

Compton scattering of photons by nonrelativistic particles is thought to play an important role in forming the radiation spectrum of many astrophysical systems. Here we derive the time-dependent photon kinetic equation that describes spontaneous and induced Compton scattering, as well as absorption and emission by static and moving media, the corresponding radiative transfer equation, and their zeroth and first angular moments, both in the system frame and in the frame comoving with the medium. We show that it is necessary to use the correct relativistic differential scattering cross section in order to obtain a photon kinetic equation that is correct to first order in Epsilon/m(sub e), T(sub e)/m(sub e), and V, where Epsilon is the photon energy, T(sub e) and m(sub e) are the electron temperature and rest mass, and V is the electron bulk velocity in units of the speed of light. We also demonstrate that the terms in the radiative transfer equation that are second order in V should usually be retained, because if the radiation energy density is sufficiently large, compared to the radiation flux, the effects of bulk Comptonization described by the terms that are second order in V can be as important as the effects described by the terms that are first order in V, even when V is small. The system- and fluid-frame equations that we derive are correct to first order in Epsilon/m(sub e). Our system-frame equations, which are correct to second order in V, may be used when V is not too large. Our fluid-frame equations, which are exact in V, may be used when V approaches 1. Both sets of equations are valid for systems of arbitrary optical depth and can therefore be used in both the free-streaming and diffusion regimes. We demonstrate that Comptonization by the electron bulk motion occurs whether or not the radiation field is isotropic or the bulk flow converges and that it is more important than thermal Comptonization if V(sup 2) is greater than 3T(sub e)/m(sub e).

Psaltis, Dimitrios↗

NASA NTRS1982

Computer program for post-flight evaluation of the control surface response for an attitude controlled missile

A FORTRAN IV coded computer program is presented for post-flight analysis of a missile's control surface response. It includes preprocessing of digitized telemetry data for time lags, biases, non-linear calibration changes and filtering. Measurements include autopilot attitude rate and displacement gyro output and four control surface deflections. Simple first order lags are assumed for the pitch, yaw and roll axes of control. Each actuator is also assumed to be represented by a first order lag. Mixing of pitch, yaw and roll commands to four control surfaces is assumed. A pseudo-inverse technique is used to obtain the pitch, yaw and roll components from the four measured deflections. This program has been used for over 10 years on the NASA/SCOUT launch vehicle for post-flight analysis and was helpful in detecting incipient actuator stall due to excessive hinge moments. The program is currently set up for a CDC CYBER 175 computer system. It requires 34K words of memory and contains 675 cards. A sample problem presented herein including the optional plotting requires eleven (11) seconds of central processor time.

Knauber, R. N.↗

NASA NTRS1988

Numerical study of the vortex burst phenomenon for delta wings

A flux-difference splitting scheme is employed to compute low-speed flows over a delta wing for angles of attack from 0 to 40 deg as steady-state solutions to the three-dimensional, Reynolds-averaged Navier-Stokes equations in their thin-layer approximation. The finite-difference scheme is made spatially second-order accurate by applying a total variation diminishing-like discretization to the inviscid fluxes and central differencing to the viscous shear fluxes. Using first-order accurate Euler backward-time differencing, an efficient implicit algorithm is constructed, which combines approximate factorization in cross planes with a symmetric planar Gauss-Seidel relaxation in the remaining third spatial direction. The geometry of the thin (maximum thickness is 0.021), slender (aspect ratio is unity), sharp-edged delta wing is taken from Hummel's (1967, 1978) wind tunnel model. Over the entire angle-of-attack range, the computed values of lift and pitching moment are in good agreement with the experimental data. Also details of the flow-fieldlike spanwise surface pressure distributions compare well with the experiment. Computed flow-field results with a bubble-type vortex burst are analyzed in detail.

Hartwich, PETER-M.↗

NASA NTRS1996

A Fundamental Study of Nucleate Pool Boiling Under Microgravity

An experimental study of incipient boiling in short-term microgravity and with a/g = +/- 1 for pool boiling was performed. Calibrated thin gold films sputtered on a smoothly polished quartz surface were used simultaneously for thermal-resistance measurements and heating of the boiling surface. The gold films were used for both transient and quasi-steady heating surface temperature measurements. Two test vessels were constructed for precise measurement and control of fluid temperature and pressure: a laboratory pool boiling vessel for the a/g = +/- 1 experiments and a pool boiling vessel designed for the 131 m free-fall in the NASA Lewis Research Center Microgravity Research Facility for the microgravity tests. Measurements included the heater surface temperature, the pressure near the heating surface, the bulk liquid temperatures. High speed photography (up to 1,000 frames per second) was used in the experiments. With high quality microgravity and the measured initial temperature of the quiescent test fluid, R113, the temperature distribution in the liquid at the moment of boiling inception resulting from an imposed step in heat flux is known with a certainty not possible previously. The types of boiling propagation across the large flat heating surface, some observed here for the first time, are categorized; the conditions necessary for their occurrence are described. Explosive boiling propagation with a striking pattern of small scale protuberances over the entire vapor mass periphery not observed previously at low heat flux levels (on the order of 5 W/cm(exp 2)) is described. For the heater surface with a/g = -1, a step in the heater surface temperature of short duration was imposed. The resulting liquid temperature distribution at the moment of boiling inception was different from that obtained with a step in heat flux.

Ervin, Jamie S.↗

NASA NTRS1994

Modeling of aircraft unsteady aerodynamic characteristics. Part 1: Postulated models

A short theoretical study of aircraft aerodynamic model equations with unsteady effects is presented. The aerodynamic forces and moments are expressed in terms of indicial functions or internal state variables. The first representation leads to aircraft integro-differential equations of motion; the second preserves the state-space form of the model equations. The formulations of unsteady aerodynamics is applied in two examples. The first example deals with a one-degree-of-freedom harmonic motion about one of the aircraft body axes. In the second example, the equations for longitudinal short-period motion are developed. In these examples, only linear aerodynamic terms are considered. The indicial functions are postulated as simple exponentials and the internal state variables are governed by linear, time-invariant, first-order differential equations. It is shown that both approaches to the modeling of unsteady aerodynamics lead to identical models.

Klein, Vladislav↗