Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “experimental algorithms”

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,333 records · Page 74

Climatespark: an In-Memory Distributed Computing Framework for Big Climate Data Analytics

The unprecedented growth of climate data creates new opportunities for climate studies, and yet big climate data pose a grand challenge to climatologists to efficiently manage and analyze big data. The complexity of climate data content and analytical algorithms increases the difficulty of implementing algorithms on high performance computing systems. This paper proposes an in-memory, distributed computing framework, ClimateSpark, to facilitate complex big data analytics and time-consuming computational tasks. Chunking data structure improves parallel I/O efficiency, while a spatiotemporal index is built for the chunks to avoid unnecessary data reading and preprocessing. An integrated, multi-dimensional, array-based data model (ClimateRDD) and ETL operations are developed to address big climate data variety by integrating the processing components of the climate data lifecycle. ClimateSpark utilizes Spark SQL and Apache Zeppelin to develop a web portal to facilitate the interaction among climatologists, climate data, analytic operations and computing resources (e.g., using SQL query and Scala/Python notebook). Experimental results show that ClimateSpark conducts different spatiotemporal data queries/analytics with high efficiency and data locality. ClimateSpark is easily adaptable to other big multiple- dimensional, array-based datasets in various geoscience domains.

Hu, Fei↗

Steady and unsteady solutions of the incompressible Navier-Stokes equations

An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The equations are solved with a line-relaxation scheme that allows the use of very large pseudotime steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. The steady-state solution of flow through a square duct with a 90-deg bend is computed, and the results are compared with experimental data. Good agreement is observed. Computations of unsteady flow over a circular cylinder are presented and compared to other experimental and computational results. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Rogers, Stuart E.↗

Point-implicit relaxation strategies for viscous, hypersonic flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's symmetric total variation diminishing scheme. Newton relaxation of the fully coupled equation set is employed on a cell-to-cell basis. Under-relaxation of the inviscid and over-relaxation of the viscous contributions to the residual are implemented. Computational work is easily partitioned among many processors in an asynchronous, dynamic mode for convergence acceleration. An overview of the physical models employed herein for thermochemical nonequilibrium is included. Several test cases and comparisons with experimental data are presented involving hypersonic flow over blunt bodies which illustrate the qualitative and quantitative capabilities of this approach.

Gnoffo, P. A.↗

Prospects of the ICESat-2 Laser Altimetry Mission for Savanna Ecosystem Structural Studies Based on Airborne Simulation Data

The next planned spaceborne lidar mission is the Ice, Cloud and land Elevation Satellite 2 (ICESat-2), which will use the Advanced Topographic Laser Altimeter System (ATLAS) sensor, a photon counting technique. To pre-validate the capability of this mission for studying three dimensional vegetation structure in savannas, we assessed the potential of the measurement approach to estimate canopy height in an oak savanna landscape. We used data from the Multiple Altimeter Beam Experimental Lidar (MABEL), an airborne photon counting lidar sensor developed by NASA's Goddard Space Flight Center. ATLAS-like data was generated using the MATLAS simulator, which adjusts MABEL data's detected number of signal and noise photons to that expected from the ATLAS instrument. Transects flown over the Tejon ranch conservancy in Kern County, California, USA were used for this work. For each transect we chose to use data from the near infrared channel that had the highest number of photons. We segmented each transect into 50 m, 25 m and 14 m long blocks and aggregated the photons in each block into a histogram based on their elevation values. We then used an automated algorithm to identify cut off points where the cumulative density of photons from the highest elevation indicates the presence of the canopy top and likewise where such cumulative density from the lowest elevation indicates the mean terrain elevation. MABEL derived height metrics were moderately correlated to discrete return lidar (DRL) derived height metrics r(sub 2) and RMSE values ranging from 0.60 to 0.73 and 2.9 m to 4.4 m respectively) but MATLAS simulation resulted in more modest correlations with DRL indices r(sub 2) ranging from 0.5 to 0.64 and RMSE from 3.6 m to 4.6 m). Simulations also indicated that the expected number of signal photons from ATLAS will be substantially lower, a situation that reduces canopy height estimation precision especially in areas of low density vegetation cover. On the basis of the simulated data, there is reason to believe that the ability of ICESat-2 to estimate height in savannas will be comparable to the original ICESat mission although the respective sensors have different measurement principles.

Gwenzi, David↗

Development of a Sunspot Tracking System

Large solar flares produce a significant amount of energetic particles which pose a hazard for human activity in space. In the hope of understanding flare mechanisms and thus better predicting solar flares, NASA's Marshall Space Flight Center (MSFC) developed an experimental vector magnetograph (EXVM) polarimeter to measure the Sun's magnetic field. The EXVM will be used to perform ground-based solar observations and will provide a proof of concept for the design of a similar instrument for the Japanese Solar-B space mission. The EXVM typically operates for a period of several minutes. During this time there is image motion due to atmospheric fluctuation and telescope wind loading. To optimize the EXVM performance an image motion compensation device (sunspot tracker) is needed. The sunspot tracker consists of two parts, an image motion determination system and an image deflection system. For image motion determination a CCD or CID camera is used to digitize an image, than an algorithm is applied to determine the motion. This motion or error signal is sent to the image deflection system which moves the image back to its original location. Both of these systems are under development. Two algorithms are available for sunspot tracking which require the use of only one row and one column of image data. To implement these algorithms, two identical independent systems are being developed, one system for each axis of motion. Two CID cameras have been purchased; the data from each camera will be used to determine image motion for each direction. The error signal generated by the tracking algorithm will be sent to an image deflection system consisting of an actuator and a mirror constrained to move about one axis. Magnetostrictive actuators were chosen to move the mirror over piezoelectrics due to their larger driving force and larger range of motion. The actuator and mirror mounts are currently under development.

Taylor, Jaime R.↗

Numerical modeling of turbulent flow in a channel

Two-dimensional incompressible turbulent flow in a channel with a backward-facing step was studied numerically by Chorin's Random Vortex Method (RVM), an algorithm capable of tracing the action of elementary turbulent eddies and their cumulative effects without imposing any restrictions upon their motions. The step occurs in one side of a channel with otherwise flat, parallel walls; its height equals 1/3, 1/4 or 1/5 the width of the channel downstream. The main objective was to investigate the behavior of the large-scale turbulent eddies in a flow and the flow characteristics in the separated shear layer, the reattached zone, and the rebuilding boundary layer after reattachment. The unsteady vorticity field and the distribution of time-averaged turbulent statistics were obtained. The effects of expansion step height and initial boundary layer state were also studied. Comparisons were made with the available experimental results. The agreement is satisfactory in the velocity profiles and in the reattachment length, and fairly good in the turbulence profiles. Also a mechanism of the development of the reattaching turbulent flow was suggested by the numerical results.

Dai, Y. W.↗

Do nonlinearities play a significant role in short term, beat-to-beat variability?

Numerous studies of short-term beat-to-beat variability in cardiovascular signals have not resolved the debate about the completeness of linear analysis techniques. This aim of this paper is to evaluate further the role of nonlinearities in short-term, beat-to-beat variability. We compared linear autoregressive moving average (ARMA) and nonlinear neural network (NN) models for predicting instantaneous heart rate (HR) and mean arterial blood pressure (BP) from past HR and BP. To evaluate these models, we used HR and BP time series from the MIMIC database. Experimental results indicate that NN-based nonlinearities do not play a significant role and suggest that ARMA linear analysis techniques provide adequate characterization of the system dynamics responsible for generating short-term, beat-to-beat variability.

NASA Discipline Cardiopulmonary↗

Surface texture characterization by angular distributions of scattered light

Work at the National Bureau of Standards to develop an on-line optical measurement device and attendant algorithms for automated optical scattering measurements of machined metal surfaces are described. The surfaces could be milled, ground or lapped, and the system is intended to categorize the resulting surface characteristics. The optical device consists of a He-Ne laser which is shone on a surface. The scattered light is captured by a semicircular array of 87 detector elements rotated over the scanned area. The light signals are processed through a digital voltmeter and then an A/D converter. The signals are then stored for later comparisons with optical scattering data obtained by optical instruments used with stylus instruments for characterizing surface topographies. A theoretical model has been defined which relates light scattering and the characteristics of the surface roughness. Initial experimental results with a sinusoidal surface have indicated that although the system can follow the trend of the roughness, the roughness amplitude is as yet uncertainly defined and the computations require excessive time.

Gilsinn, D. E.↗

Turbomachinery noise

Summarized here are key advances in experimental techniques and theoretical applications which point the way to a broad understanding and control of turbomachinery noise. On the experimental side, the development of effective inflow control techniques makes it possible to conduct, in ground based facilities, definitive experiments in internally controlled blade row interactions. Results can now be valid indicators of flight behavior and can provide a firm base for comparison with analytical results. Inflow control coupled with detailed diagnostic tools such as blade pressure measurements can be used to uncover the more subtle mechanisms such as rotor strut interaction, which can set tone levels for some engine configurations. Initial mappings of rotor wake-vortex flow fields have provided a data base for a first generation semiempirical flow disturbance model. Laser velocimetry offers a nonintrusive method for validating and improving the model. Digital data systems and signal processing algorithms are bringing mode measurement closer to a working tool that can be frequently applied to a real machine such as a turbofan engine. On the analytical side, models of most of the links in the chain from turbomachine blade source to far field observation point have been formulated. Three dimensional lifting surface theory for blade rows, including source noncompactness and cascade effects, blade row transmission models incorporating mode and frequency scattering, and modal radiation calculations, including hybrid numerical-analytical approaches, are tools which await further application.

Groeneweg, John F.↗

Efficient simulation of incompressible viscous flow over single and multi-element airfoils

Incompressible viscous turbulent flows over single- and multiple-element airfoils are numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The solution algorithm uses the method of pseudocompressibility with an upwind-differencing scheme for the convective fluxes and an implicit line-relaxation scheme to study high-lift take-off and landing configurations and to compute lift and drag at various angles of attack up to stall. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil. The approach used for multiple-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared: a patched system of grids and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface-pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time on a CRAY YMP per element in the airfoil configuration.

Rogers, Stuart E.↗

Efficient simulation of incompressible viscous flow over multi-element airfoils

The incompressible, viscous, turbulent flow over single and multi-element airfoils is numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The computer code uses the method of pseudo-compressibility with an upwind-differencing scheme for the convective fluxes and an implicit line-relaxation solution algorithm. The motivation for this work includes interest in studying the high-lift take-off and landing configurations of various aircraft. In particular, accurate computation of lift and drag at various angles of attack, up to stall, is desired. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil; an accurate prediction of stall is obtained. The approach used for multi-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared: a patched system of grids, and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time (on a CRAY YMP) per element in the airfoil configuration.

Rogers, Stuart E.↗

Computation of separated flows by a vortex-tracing algorithm

Numerical solutions for two-dimensional, time-dependent, separated flows around bodies are obtained, using a new version of the vortex method. This method provides an efficient representation of flows involving large regions of separation. The modifications incorporated in the new version, which improve its accuracy, versatility, and computing speed, are described. The computer cost is only of the order of the 3/2 power on N, instead of N-squared, for each step with N vortices. Arbitrary shapes can be treated; a conformal mapping is not required. Special attention is paid to the viscous character of the solution and to the accurate computation of the pressure distribution at the body surface. The vortex solution for the outer flow is coupled to an inner solution for the attached part of the boundary layer. Numerical results are presented for several bluff bodies exhibiting dependence on Reynolds number, for stationary airfoils under steady or transient conditions and for oscillating airfoils, including dynamic stalls. These results are compared with other available results, analytical or experimental, and demonstrate the enhanced reliability and accuracy of the improved method.

Spalart, P. R.↗

Experimental testing using a remotely operated flap with a 2D high lift model

The problem of efficient wind tunnel testing for multi-element airfoils was first addressed by the author during a previous ASEE fellowship. A modern three element model with internal actuators to position a flap in two degrees of freedom was designed and later built. Some preliminary testing proved that the approach was viable. The purpose of this summer's work was to fully develop experimental methods including efficient data acquisition. The final goal is to develop dense data sets for both lift and drag measurements as a function of flap position for both take-off and landing configurations. The model has a span of 36 in. and chord of 18 in. and is currently being fitted for a 3 ft. x 4 ft. low speed wind tunnel. The flap was reworked to allow all pressure taps to function after initial tests showed two blocked ports. The serial method of obtaining pressures from the surface taps was found to be exceedingly slow so a new method using 12 pressure transducers and a 12 port parallel scanning valve were developed. A new automated data acquisition and control algorithm was developed using LabView software and a PC platform. Flow two-dimensionality is currently under investigation with boundary layer control by blowing; this was previously omitted for initial testing. By the end of the summer a detailed data set (uncorrected) consisting of lift coefficient versus flap position for the landing configuration should be available.

Landman, Drew↗

Solar Thermal Concept Evaluation

Concentrated solar thermal energy can be utilized in a variety of high temperature applications for both terrestrial and space environments. In each application, knowledge of the collector and absorber's heat exchange interaction is required. To understand this coupled mechanism, various concentrator types and geometries, as well as, their relationship to the physical absorber mechanics were investigated. To conduct experimental tests various parts of a 5,000 watt, thermal concentrator, facility were made and evaluated. This was in anticipation at a larger NASA facility proposed for construction. Although much of the work centered on solar thermal propulsion for an upper stage (less than one pound thrust range), the information generated and the facility's capabilities are applicable to material processing, power generation and similar uses. The numerical calculations used to design the laboratory mirror and the procedure for evaluating other solar collectors are presented here. The mirror design is based on a hexagonal faceted system, which uses a spherical approximation to the parabolic surface. The work began with a few two dimensional estimates and continued with a full, three dimensional, numerical algorithm written in FORTRAN code. This was compared to a full geometry, ray trace program, BEAM 4, which optimizes the curvatures, based on purely optical considerations. Founded on numerical results, the characteristics of a faceted concentrator were construed. The numerical methodologies themselves were evaluated and categorized. As a result, the three-dimensional FORTRAN code was the method chosen to construct the mirrors, due to its overall accuracy and superior results to the ray trace program. This information is being used to fabricate and subsequently, laser map the actual mirror surfaces. Evaluation of concentrator mirrors, thermal applications and scaling the results of the 10 foot diameter mirror to a much larger concentrator, were studied. Evaluations, recommendations and pit falls regarding the structure, materials and facility design are presented.

Hawk, Clark W.↗

Simulation of Conformal Spiral Slot Antennas on Composite Platforms

During the course of the grant, we wrote and distributed about 12 reports and an equal number of journal papers supported fully or in part by this grant. The list of reports (title & abstract) and papers are given in Appendices A and B. This grant has indeed been instrumental in developing a robust hybrid finite element method for the analysis of complex broadband antennas on doubly curved platforms. Previous to the grant, our capability was limited to simple printed patch antennas on mostly planar platforms. More specifically: (1) mixed element formulations were developed and new edge-based prisms were introduced; (2) these elements were important in permitting flexibility in geometry gridding for most antennas of interest; (3) new perfectly matched absorbers were introduced for mesh truncations associated with highly curved surfaces; (4) fast integral algorithms were introduced for boundary integral truncations reducing CPU time from O(N-2) down to O(N-1.5) or less; (5) frequency extrapolation schemes were developed for efficient broadband performance evaluations. This activity has been successfully continued by NASA researchers; (6) computer codes were developed and extensively tested for several broadband configurations. These include FEMA-CYL, FEMA-PRISM and FEMA-TETRA written by L. Kempel, T. Ozdemir and J. Gong, respectively; (7) a new infinite balun feed was designed nearly constant impedance over the 800-3000 MHz operational band; (8) a complete slot spiral antenna was developed, fabricated and tested at NASA Langley. This new design is a culmination of the projects goals and integrates the computational and experimental efforts. this antenna design resulted in a U.S. patent and was revised three times to achieve the desired bandwidth and gain requirements from 800-3000 MHz.

Volakis, J. L.↗

Lidar measurements of cloud extinction coefficient distribution and its forward scattering phase function according to multiply scattered lidar returns

The study of the inversion algorithm for the single scatter lidar equation, for quantitative determination of cloud (or aerosol) optical properties, has received much attention over the last thirty years. Some of the difficulties associated with the solution of this equation are not yet solved. One problem is that a single scatter lidar equation has two unknowns. Because of this, the determination of the far-end boundary value, in the case of Klett's algorithm, is a problem if the atmosphere is optically inhomogeneous. Another difficulty concerns multiple scattering. There is a large error in the extinction distribution solution, in many cases, if only the single scattering component is considered, while neglecting the multiple scattering component. However, the use of multiple scattering in the remote sensing of aerosol or cloud optical properties is promising. In our early study, an inversion method for simultaneous determination of the cloud (or aerosol) Extinction Coefficient Distribution (ECD) and its Forward Scattering Phase Function (FSPF) was proposed according to multiply scattered lidar returns with two fields of view for the receiver. The method is based on a parameterized multiple scatter lidar equation. This paper is devoted to further numerical tests and an experimental study of lidar measurements of cloud ECD and FSPF using this method.

Qiu, Jinhuan↗

A globally efficient means of distributing UTC time and frequency through GPS

Time and frequency outputs comparable in quality to the best laboratories have been demonstrated on an integrated system suitable for field application on a global basis. The system measures the time difference between 1 pulse-per-second (pps) signals derived from local primary frequency standards and from a multi-channel GPS C/A receiver. The measured data is processed through optimal SA Filter algorithms that enhance both the stability and accuracy of GPS timing signals. Experiments were run simultaneously at four different sites. Even with large distances between sites, the overall results show a high degree of cross-correlation of the SA noise. With sufficiently long simultaneous measurement sequences, the data shows that determination of the difference in local frequency from an accepted remote standard to better than 1 x 10(exp -14) is possible. This method yields frequency accuracy, stability, and timing stability comparable to that obtained with more conventional common-view experiments. In addition, this approach provides UTC(USNO MC) in real time to an accuracy better than 20 ns without the problems normally associated with conventional common-view techniques. An experimental tracking loop was also set up to demonstrate the use of enhanced GPS for dissemination of UTC(USNO MC) over a wide geographic area. Properly disciplining a cesium standard with a multi-channel GPS receiver, with additional input from USNO, has been found to permit maintaining a timing precision of better than 10 ns between Palo Alto, CA and Washington, DC.

Kusters, John A.↗

Network Analysis of Rodent Transcriptomes in Spaceflight

Network analysis methods leverage prior knowledge of cellular systems and the statistical and conceptual relationships between analyte measurements to determine gene connectivity. Correlation and conditional metrics are used to infer a network topology and provide a systems-level context for cellular responses. Integration across multiple experimental conditions and omics domains can reveal the regulatory mechanisms that underlie gene expression. GeneLab has assembled rich multi-omic (transcriptomics, proteomics, epigenomics, and epitranscriptomics) datasets for multiple murine tissues from the Rodent Research 1 (RR-1) experiment. RR-1 assesses the impact of 37 days of spaceflight on gene expression across a variety of tissue types, such as adrenal glands, quadriceps, gastrocnemius, tibalius anterior, extensor digitorum longus, soleus, eye, and kidney. Network analysis is particularly useful for RR-1 -omics datasets because it reinforces subtle relationships that may be overlooked in isolated analyses and subdues confounding factors. Our objective is to use network analysis to determine potential target nodes for therapeutic intervention and identify similarities with existing disease models. Multiple network algorithms are used for a higher confidence consensus.

genomics↗