Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 1,459 records · Page 81

Realistic Covariance Prediction for the Earth Science Constellation

Routine satellite operations for the Earth Science Constellation (ESC) include collision risk assessment between members of the constellation and other orbiting space objects. One component of the risk assessment process is computing the collision probability between two space objects. The collision probability is computed using Monte Carlo techniques as well as by numerically integrating relative state probability density functions. Each algorithm takes as inputs state vector and state vector uncertainty information for both objects. The state vector uncertainty information is expressed in terms of a covariance matrix. The collision probability computation is only as good as the inputs. Therefore, to obtain a collision calculation that is a useful decision-making metric, realistic covariance matrices must be used as inputs to the calculation. This paper describes the process used by the NASA/Goddard Space Flight Center's Earth Science Mission Operations Project to generate realistic covariance predictions for three of the Earth Science Constellation satellites: Aqua, Aura and Terra.

Duncan, Matthew↗

Highly Scalable Matching Pursuit Signal Decomposition Algorithm

Matching Pursuit Decomposition (MPD) is a powerful iterative algorithm for signal decomposition and feature extraction. MPD decomposes any signal into linear combinations of its dictionary elements or atoms . A best fit atom from an arbitrarily defined dictionary is determined through cross-correlation. The selected atom is subtracted from the signal and this procedure is repeated on the residual in the subsequent iterations until a stopping criterion is met. The reconstructed signal reveals the waveform structure of the original signal. However, a sufficiently large dictionary is required for an accurate reconstruction; this in return increases the computational burden of the algorithm, thus limiting its applicability and level of adoption. The purpose of this research is to improve the scalability and performance of the classical MPD algorithm. Correlation thresholds were defined to prune insignificant atoms from the dictionary. The Coarse-Fine Grids and Multiple Atom Extraction techniques were proposed to decrease the computational burden of the algorithm. The Coarse-Fine Grids method enabled the approximation and refinement of the parameters for the best fit atom. The ability to extract multiple atoms within a single iteration enhanced the effectiveness and efficiency of each iteration. These improvements were implemented to produce an improved Matching Pursuit Decomposition algorithm entitled MPD++. Disparate signal decomposition applications may require a particular emphasis of accuracy or computational efficiency. The prominence of the key signal features required for the proper signal classification dictates the level of accuracy necessary in the decomposition. The MPD++ algorithm may be easily adapted to accommodate the imposed requirements. Certain feature extraction applications may require rapid signal decomposition. The full potential of MPD++ may be utilized to produce incredible performance gains while extracting only slightly less energy than the standard algorithm. When the utmost accuracy must be achieved, the modified algorithm extracts atoms more conservatively but still exhibits computational gains over classical MPD. The MPD++ algorithm was demonstrated using an over-complete dictionary on real life data. Computational times were reduced by factors of 1.9 and 44 for the emphases of accuracy and performance, respectively. The modified algorithm extracted similar amounts of energy compared to classical MPD. The degree of the improvement in computational time depends on the complexity of the data, the initialization parameters, and the breadth of the dictionary. The results of the research confirm that the three modifications successfully improved the scalability and computational efficiency of the MPD algorithm. Correlation Thresholding decreased the time complexity by reducing the dictionary size. Multiple Atom Extraction also reduced the time complexity by decreasing the number of iterations required for a stopping criterion to be reached. The Course-Fine Grids technique enabled complicated atoms with numerous variable parameters to be effectively represented in the dictionary. Due to the nature of the three proposed modifications, they are capable of being stacked and have cumulative effects on the reduction of the time complexity.

Christensen, Daniel↗

Dynamic Analysis of Sounding Rocket Pneumatic System Revision

The recent fusion of decades of advancements in mathematical models, numerical algorithms and curve fitting techniques marked the beginning of a new era in the science of simulation. It is becoming indispensable to the study of rockets and aerospace analysis. In pneumatic system, which is the main focus of this paper, particular emphasis will be placed on the efforts of compressible flow in Attitude Control System of sounding rocket.

Armen, Jerald↗

Use of MODIS Cloud Top Pressure to Improve Assimilation Yields of AIRS Radiances in GSI

Improvements to global and regional numerical weather prediction have been demonstrated through assimilation of data from NASA's Atmospheric Infrared Sounder (AIRS). Current operational data assimilation systems use AIRS radiances, but impact on regional forecasts has been much smaller than for global forecasts. Previously, it has been shown that cloud top designation associated with quality control procedures within the Gridpoint Statistical Interpolation (GSI) system used operationally by a number of Joint Center for Satellite Data Assimilation (JCSDA) partners may not provide the best representation of cloud top pressure (CTP). Because this designated CTP determines which channels are cloud-free and, thus, available for assimilation, ensuring the most accurate representation of this value is imperative to obtaining the greatest impact from satellite radiances. This paper examines the assimilation of hyperspectral sounder data used in operational numerical weather prediction by comparing analysis increments and numerical forecasts generated using operational techniques with a research technique that swaps CTP from the Moderate-resolution Imaging Spectroradiometer (MODIS) for the value of CTP calculated from the radiances within GSI.

Zavodsky, Bradley↗

Thermionic Properties of Carbon Based Nanomaterials Produced by Microhollow Cathode PECVD

Thermionic emission is the process in which materials at sufficiently high temperature spontaneously emit electrons. This process occurs when electrons in a material gain sufficient thermal energy from heating to overcome the material's potential barrier, referred to as the work function. For most bulk materials very high temperatures (greater than 1500 K) are needed to produce appreciable emission. Carbon‐based nanomaterials have shown significant promise as emission materials because of their low work functions, nanoscale geometry, and negative electron affinity. One method of producing these materials is through the process known as microhollow cathode PECVD. In a microhollow cathode plasma, high energy electrons oscillate at very high energies through the Pendel effect. These high energy electrons create numerous radical species and the technique has been shown to be an effective method of growing carbon based nanomaterials. In this work, we explore the thermionic emission properties of carbon based nanomaterials produced by microhollow cathode PECVD under a variety of synthesis conditions. Initial studies demonstrate measureable current at low temperatures (approximately 800 K) and work functions (approximately 3.3 eV) for these materials.

Haase, John R.↗

Using Graph Coloring to Compute Total Derivatives More Efficiently in OpenMDAO

When they are applicable, gradient based optimization algorithms are the most efficient way to solve design optimization problems. Although gradient based methods are generally efficient, they can be made significantly more so through the usage of analytic techniques to compute the necessary total derivatives. The traditional forward (direct) and reverse (adjoint) analytic techniques have computational costs that scale linearly with the number of design variables and the number of constraints, respectively. In this work, we present an application of a graph coloring algorithm to the analytic techniques for computing total derivative Jacobians in order to achieve much better computational scaling than the pure analytic methods can provide alone. A detailed theoretical explanation of how coloring algorithms interact with analytic derivative methods is presented that illustrates specific types of sparsity patterns that must be present in total derivative Jacobians in order for this coloring technique to be effective. The new technique has been implemented as a feature in the OpenMDAO framework and the implementation is demonstrated on two example problems. The performance on the example problems up to 50% reduction in compute cost for optimizations with bi-directional coloring compared to traditional constraint aggregation. Additionally, the results show how coloring technique alleviates some of the numerical difficulties that constraint aggregation can cause, leading to the ability to solve larger problems. It is expected that the new method will have wide applicability to multidisciplinary optimization problems, and that its availability in OpenMDAO will offer significant computational savings for users without the need for them to implement the coloring algorithm themselves.

Trajectory↗

Pterodactyl: Guidance and Control of a Symmetric Deployable Entry Vehicle using an Aerodynamic Control System

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles (DEVs). This paper presents the guidance and control design and analysis for a mechanically-deployed DEV with a symmetric aeroshell and aerodynamic control surfaces. We present guidance tuning using both bank angle modulation with the Fully-Numerical Predictor Corrector Entry Guidance technique (FNPEG) and angle of attack and sideslip modulation with FNPEG with uncoupled range control (FNPEG URC). We show the control system design and simulation results for tracking both types of guidance commands.

Benjamin W L Margolis↗

Pterodactyl: Guidance and Control of a Symmetric Deployable Entry Vehicle using an Aerodynamic Control System

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles (DEVs). This paper presents the guidance and control design and analysis for a mechanically-deployed DEV with a symmetric aeroshell and aerodynamic control surfaces. We present guidance tuning using both bank angle modulation with the Fully-Numerical Predictor Corrector Entry Guidance technique (FNPEG) and angle of attack and sideslip modulation with FNPEG with uncoupled range control (FNPEG URC). We show the control system design and simulation results for tracking both types of guidance commands.

Benjamin Margolis↗

Fault Detection and Performance Monitoring of Propellers in Electric UAV

Unmanned aerial vehicles (UAVs) are used in various industries such as agriculture and logistics, to name but a few, where their applications are beyond basic mapping, surveillance, and photography. In near future, UAVs are expected to be used in package delivery service and larger electric vertical takeoff and landing vehicles will be employed for urban air mobility applications (air taxi). Thus, several electric propulsion systems will enter the low-altitude airspace with frequent take offs and landings. To achieve state-of-the-art safety standards under such high traffic density, UAVs will require in-time fault detection and performance monitoring of critical powertrain components. This work focuses on propeller blade performance and damage detection in electric UAVs. Propellers are the fastest moving component in an UAV; even a minor defect in the propeller blades could cause performance deterioration, with consequent challenges in flying through the planned trajectory or adhere to the safety requirements of the operation. Monitoring and updating aerodynamic efficiency of each rotor would therefore enable the detection of off-nominal propeller conditions thus magnifying the state-awareness of powertrain monitoring systems based on the acquired electrical signals. We use an extended Kalman filter-based parameter estimation algorithm that incorporates time history responses from UAV powertrain in conjunction with a full powertrain system model. Propeller fault detection is achieved by incorporating the aerodynamic parameters of propeller into the powertrain model. The proposed technique is successfully validated with numerical simulations.

Unmanned aerial vehicles↗

Pterodactyl: Coupled 6-Dof Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles (DEVs). This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance technique (FNPEG) and angle of attack and sideslip modulation with FNPEG uncoupled range control (FNPEG URC). These, with a linear quadratic regulator (LQR) controller, had previously been presented in uncoupled 3-DOF trajectories. This work will show results from fully coupled 6-DOF simulations, leveraging recent advancements in trajectory simulation software, namely the Julia-based Genesis package. Preliminary results show good performance for angle of attack and sideslip modulation when using a controller designed at high dynamic pressure conditions.

range targeting↗

Computationally Efficient Frequency Domain Method for Dynamic Control Derivative Estimation with Application to Transonic Truss-Braced Wing

This paper presents a computationally efficient method for dynamic control derivative estimation technique via control surface oscillation numerical experiments. Unsteady RANS CFD simulations in FUN3D are performed to simulate the control surface oscillations with a prescribed truncated square wave containing sufficient frequencies of interest. The truncated square wave oscillation offers the computational efficiency which reduces the computational cost by almost an order of magnitude compared to a sine wave oscillation. The nonlinear effect of large control surface oscillation amplitudes creates a spillover effect whereby the frequency response at the same input frequencies contains not only the linear aerodynamic response but also nonlinear aerodynamic response. A correction procedure is developed to remove the spillover effect from the linear aerodynamic response. A frequency-domain regression is performed to estimate the dynamic control derivatives after the correction. The results generally agree with the previous results obtained from the sine wave oscillation.

Stability and Control↗

Atomic Layer Deposition for Materials-Based H2 Storage: Opportunities and Limitations

The transportation demands in our growing hydrogen (H2) economy requires robust storage systems. The current commercially implemented technology in fuel cell cars relies on the well-established technology of hydrogen gas compressed to 350-700 bar, depending on the application. The compressed gas tanks in use today are bulky and cost intensive. To address this challenge, material-based storage is one of the long-term alternatives considered and constitutes the focus of this talk. Material-based storage is broadly defined as hydrogen bound to solid materials, with its binding strength varying from physisorption to porous materials, such as zeolites and metal organic frameworks, to chemisorption in (complex) metal hydrides. The ultimate targets set by the U.S. Department of Energy for this technology include a system gravimetric capacity of 6.5 wt% and volumetric capacity of 40 g/L at 100 bar, operating temperatures ranging between -40 C and +40 C and adsorption/desorption timescales of < 5 min. Storage in the form of physisorbed or chemisorbed hydrogen has guided the materials research, which metal- organic framework and (complex) metal hydrides being the most promising materials classes. A variety of these materials have met one or more of the targets, but it has remained elusive for a single material system to meet all these stringent requirements. Nano-encapsulation and low-concentration chemical additives have previously been employed separately to overcome such challenges. Functionalization via atomic layer deposition (ALD), however, offers unique characteristics that make it suitable for both, nano-encapsulation and "doping" with low-concentration additives. This deposition technique has sub-monolayer thickness control, is highly conformal in high-surface area materials and is self-limiting, i.e., once the gas-phase precursor reacts with the available surface sites, the surface reactions stop. In this presentation, we will showcase examples where ALD, more generally vapor-phase functionalization, on (complex) metal hydrides and organic frameworks has improved the material properties for H2 storage. In our first study, Al2O3 was deposited on magnesium borohydride, Mg(BH4)2, at room temperature using trimethylaluminum (TMA) and water. From our findings, encouraging initial results were obtained: the H2 desorption temperature was lowered by 60-120 degrees C, the desorbed gravimetric H2 capacity at temperatures < 250 degrees C was doubled, and the desorption kinetics increased by a factor of ~6 compared to uncoated Mg(BH4)2. However, hydrolysis reactions caused by residual surface -OH groups from the ALD water-pulse degraded the sample substantially. Through this study, the use of TMA was observed to be highly reactive with the Mg(BH4)2 surface species, leading to the strategy of exposing the sample only to TMA. The relative mole fraction of the vapor-phase additive can be precisely controlled with the number of TMA pulses and pulse duration. By tuning these parameters, we show in our second study that 10 pulses of TMA at ambient conditions were able to decrease the H2 desorption temperature by ~100 degrees C while retaining <95 % of its H2 capacity. We applied this approach with other additives such as BBr3, TiCl4 and tetrahydrofuran, and demonstrate that this unique approach opens the door to a new class of molecular additives and catalysts which cannot easily be introduced with conventional mechano-chemical or solvent-based techniques for (complex) metal-hydrides. In the case of metal organic frameworks, ALD is a promising technique to functionalize the pores with metal atoms or functional groups able to tune the gas selectivity and binding energy, for which ~15 kJ/mol has been established as the optimal value for H2 storage in sorbent materials. Vapor-phase techniques such as ALD have numerous benefits over other functionalization tools for H2 storage materials opening innumerable opportunities to the field. To exploit these opportunities, the current limitations on room temperature and water-less ALD processes needs to be overcome which will greatly expand the possibilities of encapsulation and incorporation of additives for organic frameworks and (complex) metal hydrides.

atomic layer deposition↗

Binary Oscillatory Crossflow Electrophoresis

We present preliminary results of our implementation of a novel electrophoresis separation technique: Binary Oscillatory Cross flow Electrophoresis (BOCE). The technique utilizes the interaction of two driving forces, an oscillatory electric field and an oscillatory shear flow, to create an active binary filter for the separation of charged species. Analytical and numerical studies have indicated that this technique is capable of separating proteins with electrophoretic mobilities differing by less than 10%. With an experimental device containing a separation chamber 20 cm long, 5 cm wide, and 1 mm thick, an order of magnitude increase in throughput over commercially available electrophoresis devices is theoretically possible.

Molloy, Richard F.↗

Automated monitor and control for deep space network subsystems

The problem of automating monitor and control loops for Deep Space Network (DSN) subsystems is considered and an overview of currently available automation techniques is given. The use of standard numerical models, knowledge-based systems, and neural networks is considered. It is argued that none of these techniques alone possess sufficient generality to deal with the demands imposed by the DSN environment. However, it is shown that schemes that integrate the better aspects of each approach and are referenced to a formal system model show considerable promise, although such an integrated technology is not yet available for implementation. Frequent reference is made to the receiver subsystem since this work was largely motivated by experience in developing an automated monitor and control loop for the advanced receiver.

Smyth, P.↗

Flow field analysis of aircraft configurations using a numerical solution to the three-dimensional unified supersonic/hypersonic small disturbance equations, part 1

The unified small disturbance equations are numerically solved using the well-known Lax-Wendroff finite difference technique. The method allows complete determination of the inviscid flow field and surface properties as long as the flow remains supersonic. Shock waves and other discontinuities are accounted for implicity in the numerical method. This technique was programed for general application to the three-dimensional case. The validity of the method is demonstrated by calculations on cones, axisymmetric bodies, lifting bodies, delta wings, and a conical wing/body combination. Part 1 contains the discussion of problem development and results of the study. Part 2 contains flow charts, subroutine descriptions, and a listing of the computer program.

Gunness, R. C., Jr.↗

Preliminary results of polarization signatures for glacial moraines in the Mono Basin, Eastern Sierra Nevada

The valleys of the Mono Basin contain several sets of lateral and terminal moraines representing multiple stages of glaciation. The semi-arid climate with slow weathering rates preserved sequences of nested younger moraines within older ones. There is a well established relative chronology and recently exposure dating provided a new set of numerical dates. The moraines span the late Wisconsin (11-25 ka) to the Illinoian (130-190 ka) glaciations. The Mono Basin area was used as a 'calibration site' to establish remote dating techniques for eventual transfer to the more inaccessible but geomorphically and climatically similar moraines of the South American Andes Mountains. Planned polarimetric synthetic aperture radar (SAR) imagery acquired by JPL AIRSAR (South American Campaign) and SIR-C (Andes super-site) are analyzed to establish chronologies of previously undated moraine sequences in a study of Pleistocene climatic change in the Southern Hemisphere. The dry climate and sparse vegetation is also favorable for correlation of ground surface roughness with radar polarization signature. The slow weathering processes acting over thousands of years reduce the size, frequency, and angularity of surface boulders while increasing soil development on the moraines. Field observations based on this hypothesis result in relative ages consistent with those inferred from nested position within the valley. Younger moraines, therefore, will appear rougher than the older smoother moraines at scales measurable at AIRSAR wavelengths. Previously documented effects of ground surface roughness on polarization signatures suggest that analysis of moraine polarization signatures can be useful for relative dating. The technique may be extended to predict numerical ages. The data set reported were acquired on 8 Sep. 1989 with the JPL Airborne SAR (AIRSAR) collecting polarimetric imagery at C- (5.6 cm), L- (24 cm), and P-band (68 cm) with a flight-line parallel to the strike of the mountains. Phase calibration was performed on the analyzed scene by setting the co-phase of a smooth lake to zero as described. Absolute amplitude calibration was not possible because corner reflectors were not deployed.

Forster, Richard R.↗

Hybridization effect on the x-ray absorption spectra for actinide materials: Application to PuB 4

Studying the local moment and 5 f -electron occupations sheds insight into the electronic behavior in actinide materials. X-ray absorption spectroscopy (XAS) has been a powerful tool to reveal the valence electronic structure when assisted with theoretical calculations. However, the analysis currently taken in the community on the branching ratio of the XAS spectra generally does not account for the hybridization effects between local f orbitals and conduction states. In this paper, we discuss an approach which employs the density functional theory plus Gutzwiller rotationally invariant slave boson method to obtain a local Hamiltonian for the single-impurity Anderson model, and calculates the XAS spectra by the exact diagonalization (ED) method. A customized numerical routine was implemented for the ED XAS part of the calculation. By applying this technique to the recently discovered 5 f -electron topological Kondo insulator Pu B 4 , we determined the signature of 5 f -electronic correlation effects in the theoretical x-ray spectra. Furthermore, we found that the Pu 5 f - 6 d hybridization effect provides an extra channel to mix the j = 5 / 2 and 7 / 2 orbitals in the 5 f valence. As a consequence, the resulting electron occupation number and spin-orbit coupling strength deviate from the intermediate-coupling regime.

36 MATERIALS SCIENCE↗

A multistage time-stepping scheme for the Navier-Stokes equations

A class of explicit multistage time-stepping schemes is used to construct an algorithm for solving the compressible Navier-Stokes equations. Flexibility in treating arbitrary geometries is obtained with a finite-volume formulation. Numerical efficiency is achieved by employing techniques for accelerating convergence to steady state. Computer processing is enhanced through vectorization of the algorithm. The scheme is evaluated by solving laminar and turbulent flows over a flat plate and an NACA 0012 airfoil. Numerical results are compared with theoretical solutions or other numerical solutions and/or experimental data.

Swanson, R. C.↗