Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced order 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.

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

At least 343 records · Page 19

Automated procedure for design of wing structures to satisfy strength and flutter requirements

A pilot computer program was developed for the design of minimum mass wing structures under flutter, strength, and minimum gage constraints. The wing structure is idealized by finite elements, and second-order piston theory aerodynamics is used in the flutter calculation. Mathematical programing methods are used for the optimization. Computation times during the design process are reduced by three techniques. First, iterative analysis methods used to reduce significantly reanalysis times. Second, the number of design variables is kept small by not using a one-to-one correspondence between finite elements and design variables. Third, a technique for using approximate second derivatives with Newton's method for the optimization is incorporated. The program output is compared witH previous published results. It is found that some flutter characteristics, such as the flutter speed, can display discontinous dependence on the design variables (which are the thicknesses of the structural elements). It is concluded that it is undesirable to use such quantities in the formulation of the flutter constraint.

Haftka, R. T.↗

Decomposition-aggregation stability analysis of the spinning Skylab

Stability of an 11-th order linear model of the spinning Skylab is determined by the decomposition-aggregation method based upon the comparison principle and vector Liapunov functions. To reduce the inherent conservativeness of the method an optimization problem is formulated and resolved producing the optimum comparison system. The system provides the best estimate of the stability region of the important structural parameter - asymmetry in the boom settings.

Cuk, S. M.↗

Analysis of transonic flow about lifting wing-body configurations

An analytical solution was obtained for the perturbation velocity potential for transonic flow about lifting wing-body configurations with order-one span-length ratios and small reduced-span-length ratios and equivalent-thickness-length ratios. The analysis is performed with the method of matched asymptotic expansions. The angles of attack which are considered are small but are large enough to insure that the effects of lift in the region far from the configuration are either dominant or comparable with the effects of thickness. The modification to the equivalence rule which accounts for these lift effects is determined. An analysis of transonic flow about lifting wings with large aspect ratios is also presented.

Barnwell, R. W.↗

OH measurement by laser light absorption

Since the first attempt to measure atmospheric hydroxyl radicals by optical absorption in 1975 (Perner et al., 1976) this method has been continuously developed further and its major obstacles and limitations are known today. The laser beam needs to be expanded in order to reduce the beam divergence. At the same time the energy density of the laser beam which produces OH via ozone photolysis is reduced to such an extent that the self-produced OH concentration ranges well below the atmospheric value. Atmospheric absorptions should be observed over a wide spectral range so that not only the OH radicals are properly identified by several rotational lines but their absorption can be corrected for interfering absorptions from other air constituents as SO2, CH2O, CS2, etc., which can be identified in a wide spectral range with more confidence. Air turbulence demands fast spectral scanning or probing on and off the absorption line. Energy requirements should be kept small in field operations. In the experiment frequency doubled dye laser pulses at 308 nm are produced. The picosecond light pulses are expected to show a smooth profile (light intensity against wavelength) which will be broadened to the required spectral width according to the uncertainty principle. The pump laser will be an optoacoustically modulated Nd:YAG laser.

Perner, D.↗

Theoretical and experimental study of flow-control devices for inlets of indraft wind tunnels

The design of closed circuit wind tunnels has historically been performed using rule of thumb which have evolved over the years into a body of useful guidelines. The development of indraft wind tunnels, however, has not been as well documented. The design of indraft wind tunnels is therefore generally performed using a more intuitive approach, often resulting in a facility with disappointing flow quality. The primary problem is a lack of understanding of the flow in the inlet as it passes through the required antiturbulence treatment. For wind tunnels which employ large contraction ratio inlets, this lack of understanding is not serious since the relatively low velocity of the flow through the inlet treatment reduces the sensitivity to improper inlet design. When designing a small contraction ratio inlet, much more careful design is needed in order to reduce the flow distortions generated by the inlet treatment. As part of the National Full Scale Aerodynamics Complex Modification Project, 2-D computational methods were developed which account for the effect of both inlet screens and guide vanes on the test section velocity distribution. Comparisons with experimental data are presented which indicate that the methods accurately compute the flow distortions generated by a screen in a nonuniform velocity field. The use of inlet guide vanes to eliminate the screen induced distortion is also demonstrated both computationally and experimentally. Extensions of the results to 3-D is demonstrated and a successful wind tunnel design is presented.

Ross, James C.↗

Practical model reduction for flexible structures

The most accepted methods for developing models of flexible structures for the analysis, design and simulation of control systems are finite element methods and extraction from experimental data. Using either of these methods, models are often produced with modes that do not significantly add to the fidelity of the model. This paper presents techniques for eliminating these modes. Particular attention is paid to multiple-input, multiple-output systems. First, control system models developed for flexible structures, using finite element methods and experimental data, are briefly discussed. the shortcomings of using models with unnecessary and/or residual modes are delineated. Then, two techniques for reducing the order of models are presented; the first is applicable to single-input, single-output systems and the second for multiple-input, multiple-output systems. Finally, both are illustrated using model data from the proposed NASA Shuttle-C.

Mitchell, Jerrel R.↗

Robust Assignment Of Eigensystems For Flexible Structures

Improved method for placement of eigenvalues and eigenvectors of closed-loop control system by use of either state or output feedback. Applied to reduced-order finite-element mathematical model of NASA's MAST truss beam structure. Model represents deployer/retractor assembly, inertial properties of Space Shuttle, and rigid platforms for allocation of sensors and actuators. Algorithm formulated in real arithmetic for efficient implementation. Choice of open-loop eigenvector matrix and its closest unitary matrix believed suitable for generating well-conditioned eigensystem with small control gains. Implication of this approach is that element of iterative search for "optimal" unitary matrix appears unnecessary in practice for many test problems.

Juang, Jer-Nan↗

Exercise Countermeasures for Bone Loss During Space Flight: A Method for the Study of Ground Reaction Forces and their Implications for Bone Strain

Effective countermeasures to prevent loss of bone mineral during long duration space flight remain elusive. Despite an exercise program on MIR flights, the data from LeBlanc et al. (1996) indicated that there was still a mean rate of loss of bone mineral density in the proximal femur of 1.58% per month (n=18, flight duration 4 - 14.4 months). The specific mechanisms regulating bone mass are not known, but most investigators agree that bone maintenance is largely dependent upon mechanical demand and the resultant local bone strains. A plausible hypothesis is that bone loss during space flight, such as that reported by LeBlanc et al. (1996), may result from failure to effectively load the skeleton in order to generate localized bone strains of sufficient magnitude to prevent disuse osteoporosis. A variety of methods have been proposed to simulate locomotor exercise in reduced gravity. In such simulations, and in an actual microgravity environment, a gravity replacement load (GRL) must always be added to return the exercising subject to the support surface and the resulting skeletal load is critically dependent upon the magnitude of the GRL. To our knowledge, GRLs during orbital flight have only been measured once (on STS 81) and it is likely that most or all prior treadmill exercise in space has used GRLs that were less than one body weight. McCrory (1997) has shown that subjects walking and running in simulated zero-G can tolerate GRLs of 1 if an appropriate harness is used. Several investigators have attempted to measure in vivo strains and forces in the bones of humans, but have faced ethical and technical limitations. The anteromedial aspect of the tibial midshaft has been a common site for the placement of strain gauges; one reason to measure strains in the anterior tibia is that this region is surgically accessible. Aamodt et al. (1997) were able to measure strains on the lateral surface of the proximal femur only because their experimental subjects were already scheduled for hip surgery. Lu et al. (1997) used an instrumented massive proximal femoral prosthesis along with electromyographic measurements to demonstrate that femoral forces depend on muscular activity. These analyses of in vivo bone mechanics are valuable. The invasive nature of the procedures involved, however, limits both the number of subjects and the number of strain gauge locations. Further, the results of these studies may be confounded by the inclusion of subjects with pathological conditions. Gross et al. (1992) measured strain at three locations on the equine third metacarpal and used those data to construct a computer model of the internal strain environment of the bone. An analogous placement of multiple gauges in living humans would be difficult and potentially hazardous because of the depth of soft tissue overlying the tibia and femur.

Peterman, M.↗

Direct Detection Doppler Lidar for Spaceborne Wind Measurement

The theory of double edge lidar techniques for measuring the atmospheric wind using aerosol and molecular backscatter is described. Two high spectral resolution filters with opposite slopes are located about the laser frequency for the aerosol based measurement or in the wings of the Rayleigh - Brillouin profile for the molecular measurement. This doubles the signal change per unit Doppler shift and improves the measurement accuracy by nearly a factor of 2 relative to the single edge technique. For the aerosol based measurement, the use of two high resolution edge filters reduces the effects of background, Rayleigh scattering, by as much as an order of magnitude and substantially improves the measurement accuracy. Also, we describe a method that allows the Rayleigh and aerosol components of the signal to be independently determined. A measurement accuracy of 1.2 m/s can be obtained for a signal level of 1000 detected photons which corresponds to signal levels in the boundary layer. For the molecular based measurement, we describe the use of a crossover region where the sensitivity of a molecular and aerosol-based measurement are equal. This desensitizes the molecular measurement to the effects of aerosol scattering and greatly simplifies the measurement. Simulations using a conical scanning spaceborne lidar at 355 nm give an accuracy of 2-3 m/s for altitudes of 2-15 km for a 1 km vertical resolution, a satellite altitude of 400 km, and a 200 km x 200 km spatial.

Korb, C. Laurence↗

Optimizing Vetoes for Gravitational-wave Transient Searches

Interferometric gravitational-wave detectors like LIGO, GEO600 and Virgo record a surplus of information above and beyond possible gravitational-wave events. These auxiliary channels capture information about the state of the detector and its surroundings which can be used to infer potential terrestrial noise sources of some gravitational-wave-like events. We present an algorithm addressing the ordering (or equivalently optimizing) of such information from auxiliary systems in gravitational-wave detectors to establish veto conditions in searches for gravitational-wave transients. The procedure was used to identify vetoes for searches for unmodelled transients by the LIGO and Virgo collaborations during their science runs from 2005 through 2007. In this work we present the details of the algorithm; we also use a limited amount of data from LIGO's past runs in order to examine the method, compare it with other methods, and identify its potential to characterize the instruments themselves. We examine the dependence of Receiver Operating Characteristic curves on the various parameters of the veto method and the implementation on real data. We find that the method robustly determines important auxiliary channels, ordering them by the apparent strength of their correlations to the gravitational-wave channel. This list can substantially reduce the background of noise events in the gravitational-wave data. In this way it can identify the source of glitches in the detector as well as assist in establishing confidence in the detection of gravitational-wave transients.

Transient↗

High Order Numerical Methods for LES of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten's subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.

High Order↗

High Order Numerical Methods for LES of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten’s subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.

LES↗

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗

An interactive editor for definition of touch-sensitive zones for a graphic display

In the continuing effort to develop more efficient man-machine communications methods, touch displays have shown potential as straightforward input systems. The development of software necessary to handle such systems, however, can become tedious. In order to reduce the need for redundant programming, a touch editor has been developed which allows a programmer to interactively define touch-sensitive areas for a graphic display. The information produced during the editing process is written to a data file, which can be accessed easily when needed by an application program. This paper outlines the structure, logic, and use of the editor, as well as the hardware with which it is presently compatible.

Monroe, Burt L., III↗

Parameter identification methods for improving structural dynamic models

There is an increasing need to develop Parameter Identification methods for improving structural dynamic models, based on the inability of engineers to produce mathematical models which correlate with experimental data. This research explores the efficiency of combining Component Mode Synthesis (substructuring) methods with Parameter Identification procedures in order to improve analytical modeling of structural components and their connections. Improvements are computed in terms of physical stiffness and damping parameters in order that the physical characteristics of the model can be better understood. Connections involving both viscous and friction damping are investigated. Substructuring methods are utilized to reduce the complexity of the identification problem. Component and inter-component structural connection properties are evaluated and identified independently, thus simplifying the identification problem. It is shown that modal test data is effective for identifying modeling problems associated with structural components, and for determining the stiffness and damping properties of intercomponent connections. In general, Parameter Identification is improved when greater quantities of experimental data are available.

Lawrence, Charles↗

New literal approximations for the longitudinal dynamic characteristics of flexible flight vehicles

The goal of the literal approximation method is to obtain simple literal (analytical) approximations for key dynamic characteristics of flexible flight vehicles. A basic question regarding the method is its usefulness as an additional design tool for existing design and simulation procedures. Two aspects of this question are: (1) ease of derivation and use of the literal approximations, and (2) the suitability of one set of literal approximations to describe the dynamics of a large set of significantly different vehicles. These issues are addressed by incorporating symbolic manipulation software into the literal approximation method for the analysis of a fifth order model of the longitudinal dynamics of a flexible flight vehicle. The automated literal approximation generated in this fashion reduces the manual derivation time by an approximate factor of four. A single set of literal approximations is shown to provide adequate approximations for the dynamics of significantly different flight vehicles configurations, such as an aircraft, a missile, and a hypersonic vehicle.

Livneh, Rafael↗

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Yu, Sheng-Tao↗

Nutrient Recovery of Plant Leachates Under Thermal, Biological, and Photocatalytic Pretreatments

Nutrient recovery has always been a problem for long distance and long-term space missions. To allow humans to man these missions, a steady source of oxygen, water, and food are necessary for survival beyond Earth's atmosphere. Plants are currently an area of interest since they are capable of providing all three resources for life sustainability. We are currently interested in nutrient recovery for future plant growth and simple aqueous leachate extractions can recover some of the nutrients. However, leaching plants also removes water-soluble organic plant wastes, which inhibits plant growth if not separated properly. To combat the issues with waste and maximize nutrient recovery, we are attempting to pre-treat the plant matter using biological, thermal, and photocatalytic methods before subjecting the solution with variable-strength acid digestion. For the biological method, the inoculums: mixed heterotrophic/nitrifying bioreactor effluent and Trichoderma vessei are used in an attempt to liberate more nutrients from the plant matter. For the thermal method, plants are subjected to varying temperatures at different retention times to determine nutrient recovery. Lastly, the photocatalytic method utilizes TiO (sub 2)'s oxidizing abilities under specific pHs and retention times to reduce organic wastes and improve nutrient gains. A final acid digestion serves to liberate nutrients even further in order to maximize recovery. So far, we have tested ideal acid digestion variables for practicality and performance in our experiments. We found that a low retention time of 10 minutes and a high acid concentration of 0.1 and 1 mole HCl were the most effective at nutrient recovery. For space travel purposes, 0.1 mole currently looks like a viable acid digestion to use since it is relatively effective and sustainable from a mass and energy balance if acid recovery can be performed on waste brines. Biological pretreatments do not look to be too effective and the thermal and photocatalytic methods may be preferred since they show a potential to recover more than 70 percent of the nutrients.

Photocatalytic↗