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 595 records · Page 33

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗

Time dependent wave envelope finite difference analysis of sound propagation

A transient finite difference wave envelope formulation is presented for sound propagation, without steady flow. Before the finite difference equations are formulated, the governing wave equation is first transformed to a form whose solution tends not to oscillate along the propagation direction. This transformation reduces the required number of grid points by an order of magnitude. Physically, the transformed pressure represents the amplitude of the conventional sound wave. The derivation for the wave envelope transient wave equation and appropriate boundary conditions are presented as well as the difference equations and stability requirements. To illustrate the method, example solutions are presented for sound propagation in a straight hard wall duct and in a two dimensional straight soft wall duct. The numerical results are in good agreement with exact analytical results.

Baumeister, K. J.↗

On the factorization of block-tridiagonals without storage constraints

In many programs solving difference equations, problem size is restricted by the number of available memory cells. A strategy has been developed to permit trade-offs between the number of floating point operations required and storage requirements for the solution of certain problems such as block tridiagonal systems of equations. This is done by recomputing some intermediate results instead of storing them. Reducing the storage to the square root of the current requirement will roughly double the number of computations. In theory, if m is the order of each sub-matrix in the block tridiagonal matrix, one can solve any linear system with only 5 sq m + 1 temporary storage cells. This method lends itself to efficient use on computers with parallel processing or vector processing architectures. On these computers the larger number of floating point operations is more than offset by the decrease in I/O and the increased percentage of vector operations made possible by this algorithm.

Merriam, M. L.↗

Active vibration control for flexible rotor by optimal direct-output feedback control

Experimental research tests were performed to actively control the rotor vibrations of a flexible rotor mounted on flexible bearing supports. The active control method used in the tests is called optimal direct-output feedback control. This method uses four electrodynamic actuators to apply control forces directly to the bearing housings in order to achieve effective vibration control of the rotor. The force actuators are controlled by an analog controller that accepts rotor displacement as input. The controller is programmed with experimentally determined feedback coefficients; the output is a control signal to the force actuators. The tests showed that this active control method reduced the rotor resonance peaks due to unbalance from approximately 250 micrometers down to approximately 25 micrometers (essentially runout level). The tests were conducted over a speed range from 0 to 10,000 rpm; the rotor system had nine critical speeds within this speed range. The method was effective in significantly reducing the rotor vibration for all of the vibration modes and critical speeds.

Nonami, Kenzou↗

Detection of potential space station control/structure interaction with CO-ST-IN

The NASA Lewis Research Center is concerned with the potential of interaction between space station controllers and the solar PV array structures. The models required to handle this problem are very large, and automated methods were developed for the transfer of data between structural dynamic and control system analysis software. These methods emphasize the need to achieve accurate coupled analysis results while using as small a model as possible. Specific tools which help the analyst in this regard include modal order techniques, the use of mode acceleration to calculate internal loads and stresses and the transfer of Craig-Bampton components to reduce problems associated with modal sufficiency. These techniques were applied to a space station model with 366 modes below 1 Hz. Attitude control, and alpha and beta joint control were simulated. The inclusion of alpha and beta joint controllers is important when examining overall space station dynamics. An initial choice of control parameters does indicate a potential for control/structure interaction during reboost. As expected this is exacerbated by increasing the rate gain and decreasing the hysteresis of the reaction control system (RCS) in order to improve rigid body performance.

Carney, Kelly↗

Active vibration control for flexible rotor by optimal direct-output feedback control

Experimental research tests were performed to actively control the rotor vibrations of a flexible rotor mounted on flexible bearing supports. The active control method used in the tests is called optimal direct-output feedback control. This method uses four electrodynamic actuators to apply control forces directly to the bearing housings in order to achieve effective vibration control of the rotor. The force actuators are controlled by an analog controller that accepts rotor displacement as input. The controller is programmed with experimentally determined feedback coefficients; the output is a control signal to the force actuators. The tests showed that this active control method reduced the rotor resonance peaks due to unbalance from approximately 250 microns down to approximately 25 microns (essentially runout level). The tests were conducted over a speed range from 0 to 10,000 rpm; the rotor system had nine critical speeds within this speed range. The method was effective in significantly reducing the rotor vibration for all of the vibration modes and critical speeds.

Nonami, K.↗

A compressed marine data set for geomagnetic field modeling

Some 13 million scalar magnetic field data points that have been collected from the world's ocean areas reside in the collection of the National Geophysical Data Center. In order to derive a suitable data set for modeling the geomagnetic field of the earth, each ship track is divided into 220 km segments. The distribution of the reduced data in position, time and local time is discussed. The along-track filtering process described has proved to be an effective method of condensing large numbers of shipborne magnetic data into a manageable and meaningful data set for field modeling. This process also provides the benefits of smoothing short-wavelength crystal anomalies, discarding data recorded during magnetically noisy periods, and assigning reasonable error estimates to be utilized in the least squares modeling.

Langel, R. A.↗

Measuring Massive Black Hole Binaries with LISA

The coalescence of two massive black holes produces gravitational waves (GWs) which can be detected by the space-based detector LISA. By measuring these waves, LISA can determine the various parameters which characterize the source. Measurements of the black hole masses and spins will provide information about the growth of black holes and their host galaxies over time. Measurements of a source's sky position and distance may help astronomers identify an electromagnetic counterpart to the GW event. The counterpart's redshift, combined with the GW-measured luminosity distance, can then be used to measure the Hubble constant and the dark energy parameter $w$. Because the potential science output is so high, it is useful to know in advance how well LISA can measure source parameters for a wide range of binaries. We calculate expected parameter estimation errors using the well-known Fisher matrix method. Our waveform model includes the physics of spin precession, as well as subleading harmonics. When these higher-order effects are not included, strong degeneracies between some parameters cause them to be poorly determined by a GW measurement. When precession and subleading harmonics are properly included, the degeneracies are broken, reducing parameter errors by one to several orders of magnitude.

Lang, Ryan N.↗

Localized Triple Modular Redundancy vs. Distributed Triple Modular Redundancy on a ProASIC3E Reprogrammable FPGA

Field programmable gate arrays (FPGA) are used in every space application. Currently, most space flight applications use radiation hardened (RH) FPGAs, which are very expensive. There is a desire to use cheaper, commercial off the shelf reprogrammable FPGAs, which are more susceptible to radiation effects known as single-event effects (SEE). The RH parts have SEE and total ionizing dose (TID) hardened elements pre-integrated into the part. This means that the designer does not need to implement any hardening techniques while configuring the device. The COTS parts on the other hand must be mitigated by design in order to insure any form of mitigation. The design techniques this project examines concern the use of localized triple modular redundancy (LTMR) and distributed triple modular redundancy (DTMR). LTMR triples every flip flop in the device architecture while DTMR triples everything except for the global routes (clocks, resets, and enables). The testing was performed on a ProASIC3E FPGA at the Texas A&M cyclotron facility. Two design architectures were used: shift registers and counters, both with LTMR and DTMR mitigation techniques. The test results prove that DTMR is more effective at reducing SEE than LTMR. We also determined that there was not a significant difference between the use of shift registers and counters for test purposes. More testing is required to obtain additional linear energy transfer values for each architecture and mitigation technique in order to determine the most cost-effective method of SEE mitigation.

McGuffey, Alex↗

(abstract) Optimal Low Thrust Trajectories Using Differential Inclusion Concepts

Low thrust propulsion systems typically have their greatest benefit for high energy missions or missions with large post-launch maneuver requirements. Missions which have been examined include main belt asteroid rendezvous, comet rendezvous, outer planet and Mercury orbiters, Pluto flyby, and solar probe missions. Low thrust mission design software used to determine these trajectories is based on two distinct formulations of the optimal control problem: the indirect and direct methods. The traditional approach (indirect) is to use the calculus of variations to obtain first order necessary conditions on the states and costates. In contrast, direct methods are conceptually different in that no explicit integration takes place. A direct method based on differential inclusion concepts has been developed and successfully used to compute low thrust trajectories. This new approach removes explicit control dependence from the problem thereby reducing the dimension of the parameter space for the nonlinear programming problem. Also when compared to other direct methods, fewer nonlinear constraints are required to represent the dynamics of the problem.

low↗

Intrastep, Stage-Value Predictors for Diagonally-Implicit Runge-Kutta Methods

To better identify the necessary attributes of good stage-value predictors (SVPs), numerous SVPs are designed for an existing: ESDIRK4(3)7L[2]SA [26] and a new: ESDIRK4(3)8L[2]SA scheme.1 Both are stifflyaccurate, stage-order two, explicit, singly-diagonally implicit Runge–Kutta (ESDIRK) schemes. Tradeoffs are studied in the parameter spaces enforcing the constraints on accuracy, linear stability, nonlinear stability and coefficient size to determine which objectives correlate with effective predictors. The SVPs are tested on three challenging external aerodynamics problems [107 − 108 degrees of freedom (DoFs)], each with a different level of stiffness. The problems include two 3D airfoils simulations and one canonical turbulence simulation. All simulations use the compressible Navier-Stokes equations (CNSE). An entropy stable spectral collocation formulation is used for discretizing the spatial terms in the equations. Simulations are performed at a wide variety of temporal error tolerances. Problems that are sufficiently stiff (e.g., lax temporal error tolerances) benefit from SVPs designed with second-order accuracy and stability properties: A-stability, and L-stability, rather than high accuracy constraints. Simulations with modest stiffness (e.g., strict error tolerances) are better suited for SVPs designed using high accuracy constraints. Designing SVPs with enhanced stability properties is tedious but worthwhile. Simulation times are reduced with optimal SVPs by as much as 100% on some stages, with combined stepwise improvements of between 50 − 100% for both methods. A comparative study is performed with the two aforementioned methods as well as four other ESDIRKs. The newly designed ESDIRK4(3)8L[2]SA with γ ≈ 1/10, proves to be the most efficient of the six tested ESDIRK schemes simulating the CNSE.

Diagonally-Implicit Runge-Kutta↗

Method and apparatus for rapid thrust increases in a turbofan engine

Upon a landing approach, the normal compressor stator schedule of a fan speed controlled turbofan engine is temporarily varied to substantially close the stators to thereby increase the fuel flow and compressor speed in order to maintain fan speed and thrust. This running of the compressor at an off-design speed substantially reduces the time required to subsequently advance the engine speed to the takeoff thrust level by advancing the throttle and opening the compressor stators.

Cornett, J. E.↗

Two-dimensional mesh embedding for Galerkin B-spline methods

A number of advantages result from using B-splines as basis functions in a Galerkin method for solving partial differential equations. Among them are arbitrary order of accuracy and high resolution similar to that of compact schemes but without the aliasing error. This work develops another property, namely, the ability to treat semi-structured embedded or zonal meshes for two-dimensional geometries. This can drastically reduce the number of grid points in many applications. Both integer and non-integer refinement ratios are allowed. The report begins by developing an algorithm for choosing basis functions that yield the desired mesh resolution. These functions are suitable products of one-dimensional B-splines. Finally, test cases for linear scalar equations such as the Poisson and advection equation are presented. The scheme is conservative and has uniformly high order of accuracy throughout the domain.

Shariff, Karim↗

Complementary Reliability-Based Decodings of Binary Linear Block Codes

This correspondence presents a hybrid reliability-based decoding algorithm which combines the reprocessing method based on the most reliable basis and a generalized Chase-type algebraic decoder based on the least reliable positions. It is shown that reprocessing with a simple additional algebraic decoding effort achieves significant coding gain. For long codes, the order of reprocessing required to achieve asymptotic optimum error performance is reduced by approximately 1/3. This significantly reduces the computational complexity, especially for long codes. Also, a more efficient criterion for stopping the decoding process is derived based on the knowledge of the algebraic decoding solution.

Fossorier, Marc P. C.↗

Application of the Finite Element Method in Atomic and Molecular Physics

The finite element method (FEM) is a numerical algorithm for solving second order differential equations. It has been successfully used to solve many problems in atomic and molecular physics, including bound state and scattering calculations. To illustrate the diversity of the method, we present here details of two applications. First, we calculate the non-adiabatic dipole polarizability of Hi by directly solving the first and second order equations of perturbation theory with FEM. In the second application, we calculate the scattering amplitude for e-H scattering (without partial wave analysis) by reducing the Schrodinger equation to set of integro-differential equations, which are then solved with FEM.

Shertzer, Janine↗

Accurate Adaptive Level Set Method and Sharpening Technique for Three Dimensional Deforming Interfaces

In this paper, we demonstrate improved accuracy of the level set method for resolving deforming interfaces by proposing two key elements: (1) accurate level set solutions on adapted Cartesian grids by judiciously choosing interpolation polynomials in regions of different grid levels and (2) enhanced reinitialization by an interface sharpening procedure. The level set equation is solved using a fifth order WENO scheme or a second order central differencing scheme depending on availability of uniform stencils at each grid point. Grid adaptation criteria are determined so that the Hamiltonian functions at nodes adjacent to interfaces are always calculated by the fifth order WENO scheme. This selective usage between the fifth order WENO and second order central differencing schemes is confirmed to give more accurate results compared to those in literature for standard test problems. In order to further improve accuracy especially near thin filaments, we suggest an artificial sharpening method, which is in a similar form with the conventional re-initialization method but utilizes sign of curvature instead of sign of the level set function. Consequently, volume loss due to numerical dissipation on thin filaments is remarkably reduced for the test problems

Kim, Hyoungin↗

Establishing and Monitoring an Aseptic Workspace for Building the MOMA Mass Spectrometer

Mars Organic Molecule Analyzer (MOMA) is an instrument suite on the European Space Agency (ESA) ExoMars 2020 Rover, and the Mass Spectrometer (MOMA-MS) is being built at Goddard Space Flight Center (GSFC). MOMA-MS is a life-detection instrument and thus falls in the most stringent category of Planetary Protection (PP) biological cleanliness requirements. Less than 0.03 spore/m2 are allowed in the instrument sample path. In order to meet these PP requirements, MOMA-MS must be built and maintained in a low bioburden environment. The MOMA-MS project at GSFC maintains three clean rooms with varying levels of bioburden control. The Aseptic Assembly Clean room has the highest level of control, applying three different bioburden reducing methods: 70% Isopropyl Alcohol (IPA), 7.5% Hydrogen Peroxide, and Ultra-Violet C (UVC) light. The three methods are used in rotation and each kills microorganisms by a different mechanism, reducing the likelihood of microorganisms developing resistance to all three. The Integration and Mars Chamber Clean rooms use less biocidal cleaning, with the option to deploy extra techniques as necessary. To support the monitoring of clean rooms and verification that MOMA-MS hardware meets PP requirements, a new Planetary Protection lab was established that currently has the capabilities of standard growth assays for spore or vegetative bacteria, rapid bioburden analysis that detects Adenosine Triphosphate (ATP), plus autoclave and Dry Heat microbial Reduction (DHMR) verification. The clean rooms are monitored for vegetative microorganisms and by rapid ATP assay, and a clear difference in bioburden is observed between the aseptic and other clean room.

Contamination Control Engineering↗

Efficacy of FTIR Analysis in Determining CO2 Loading on Diglycolamine

In support of advanced air revitalization technologies to enable human spaceflight beyond low earth orbit, performance studies have been conducted using a liquid amine, Diglycolamine (DGA) between teams at NASA's Johnson Spaceflight Center (JSC) and Ames Research Center (ARC). Liquid amines have been used in regenerable earth-based systems to remove CO2 from industrial systems as well as for closed-environment air revitalization because they can be regenerated at lower temperatures than solid sorbent systems. As an additional advantage to solid sorbent-based systems, liquid sorbents can be cycled between an adsorbing contactor and degassing chamber, thereby reducing system complexity by operation in a continuous loop. In an effort to inform a regeneration system design for micro-gravity applications, ARC has performed a number of tests to characterize the degas mechanics of DGA. In order to accurately measure the amount of CO2 captured or released by the amine, methods such as gravimetric weighing and chemical desorption are reasonable, however the first iteration test setup for a scaled down degas system required analysis on small sample sizes. Fourier-transform infrared spectroscopy (FTIR) analysis was experimentally evaluated to analyze CO2 concentration because it can produce measurements with sample sizes on the order of 100's of μL. Calibration against chemical desorption showed relatively good correlation and test data showed reasonable adherence to expected trends, however more extensive testing should be conducted to fully validate the usage of FTIR to determine CO2 loading on DGA.

Diglycolamine↗