Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computing”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Phase-space sensitivity (weight functions) of 3 MeV proton diagnostics

Abstract The weight function W ( X ) is a diagnostic sensitivity to phase-space variables X that relates the measured signal C to the distribution function F ( X ) through the equation C = ∫ W ( X ) F ( X ) d X . In the present work, an algorithm to calculate W for a diagnostic that measures 3 MeV protons produced in d(d,p)t fusion reactions between a fast ion and a thermal deuteron is developed. The emitted protons escape the tokamak on curved orbits and are detected. These curved orbits constitute effective diagnostic ‘sightlines’. The presented algorithm accounts for the complications associated with these curved sightlines. An initial calculation of time-reversed proton orbits computes effective solid angles and sightlines for the relevant range of incident proton velocity vectors. These precomputed orbits are then used within the framework of FIDASIM (2020 Plasma Phys. Cont. Fusion 62 105008) to calculate the reactivity averaged over the thermal distribution of the ‘target’ deuterons and the probability that a fast ion of specified energy and pitch has a gyroangle that is consistent with the kinematic equations along each of the sightlines. Comparisons with analytic formulas and with independent calculations for the Mega Amp Spherical Tokamak 3 MeV proton diagnostic verify the algorithm.

Physics↗
NASA NTRS1978

Ionospheric and satellite observations for studying the dynamic behavior of typhoons and the detection of severe storms and tsunamis

Atmospheric acoustic-gravity waves associated with severe thunderstorms, tornadoes, typhoons (hurricanes) and tsunamis can be studied through the coupling between the ionosphere and the troposphere. Reverse ray tracing computations of acoustic-gravity waves observed by an ionospheric Doppler sounder array show that wave sources are in the nearby storm systems and that the waves are excited prior to the storms. Results show that ionospheric observations, together with satellite observations, can contribute to the understanding of the dynamical behavior of typhoons, severe storms and tsunamis.

Hung, R. J.↗
NASA NTRS1979

Case studies of gravity waves associated with isolated tornadic storms on January 13, 1976

Penetrative convection, thunderstorms, squall lines, etc., all generate atmospheric gravity waves which can be observed by a ground-based ionospheric Doppler sounder array. Sources of these waves can be determined from reverse ray tracing computations. Case studies of gravity waves associated with isolated tornadic storms on January 13, 1976 were summarized to establish the minimum data sampling time required for correct spectral analysis and ray tracing computations. It was concluded that the data sampling time can be reduced to two to three times the wave period while still obtaining a reasonably good power spectral density. It was also demonstrated that the data sampling time can be reduced to two to three times the time delay of the wave arrival between two station pairs while still obtaining a justifiably good cross-spectral analysis. Computed source locations of the observed gravity waves are compared with conventional and satellite meteorological data.

Hung, R. J.↗
NASA NTRS1995

Off-design computer code for calculating the aerodynamic performance of axial-flow fans and compressors

An off-design axial-flow compressor code is presented and is available from COSMIC for predicting the aerodynamic performance maps of fans and compressors. Steady axisymmetric flow is assumed and the aerodynamic solution reduces to solving the two-dimensional flow field in the meridional plane. A streamline curvature method is used for calculating this flow-field outside the blade rows. This code allows for bleed flows and the first five stators can be reset for each rotational speed, capabilities which are necessary for large multistage compressors. The accuracy of the off-design performance predictions depend upon the validity of the flow loss and deviation correlation models. These empirical correlations for the flow loss and deviation are used to model the real flow effects and the off-design code will compute through small reverse flow regions. The input to this off-design code is fully described and a user's example case for a two-stage fan is included with complete input and output data sets. Also, a comparison of the off-design code predictions with experimental data is included which generally shows good agreement.

Schmidt, James F.↗
NASA NTRS1981

Computation of unsteady turbulent boundary layers with flow reversal and evaluation of two separate turbulence models

A procedure which solves the governing boundary layer equations within Keller's box method was developed for calculating unsteady laminar flows with flow reversal. This method is extended to turbulent boundary layers with flow reversal. Test cases are used to investigate the proposition that unsteady turbulent boundary layers also remain free of singularities. Turbulent flow calculations are performed. The governing equations for both models are solved. As in laminar flows, the unsteady turbulent boundary layers are free from singularities, but there is a clear indication of rapid thickening of the boundary layer with increasing flow reversal. Predictions of both turbulence models are the same for all practical purposes.

Cebeci, T.↗
DOE OSTI2021

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗
NASA NTRS2019

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↗
NASA NTRS2013

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: 1) a k-epsilon turbulence model computation on a 6.8 million point grid using wall functions, 2) a k-epsilon turbulence model computation on a 14 million point grid integrating to the wall, and 3) a k-omega turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-epsilon and k-omega computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-epsilon computations which extended from 70% to 100% span at the exit rating plane, whereas the k-omega computation had reversed flow from 95% to 100% span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-epsilon case resulted in an under-prediction in adiabatic efficiency of 8.3 points, whereas the k-omega case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗
NASA NTRS2013

Computational Study of the CC3 Impeller and Vaneless Diffuser Experiment

Centrifugal compressors are compatible with the low exit corrected flows found in the high pressure compressor of turboshaft engines and may play an increasing role in turbofan engines as engine overall pressure ratios increase. Centrifugal compressor stages are difficult to model accurately with RANS CFD solvers. A computational study of the CC3 centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort to understand potential causes of RANS CFD mis-prediction in these types of geometries. Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO Parallel Version 4 code: (1) a k-ε turbulence model computation on a 6.8 million point grid using wall functions, (2) a k-ε turbulence model computation on a 14 million point grid integrating to the wall, and (3) a k-ω turbulence model computation on the 14 million point grid integrating to the wall. It was found that all three cases compared favorably to data from inlet to impeller trailing edge, but the k-ε and k-ω computations had disparate results beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was observed in the k-ε computations which extended from 70 to 100 percent span at the exit rating plane, whereas the k-ω computation had reversed flow from 95 to 100 percent span. Compared to experimental data at near-peak-efficiency, the reversed flow region in the k-ε case resulted in an underprediction in adiabatic efficiency of 8.3 points, whereas the k-ω case was 1.2 points lower in efficiency.

Kulkarni, Sameer↗
DOE OSTI2022

Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems

Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2004

Advanced Applications of Adifor 3.0 for Efficient Calculation of First-and Second-Order CFD Sensitivity Derivatives

This final report will document the accomplishments of the work of this project. 1) The incremental-iterative (II) form of the reverse-mode (adjoint) method for computing first-order (FO) aerodynamic sensitivity derivatives (SDs) has been successfully implemented and tested in a 2D CFD code (called ANSERS) using the reverse-mode capability of ADIFOR 3.0. These preceding results compared very well with similar SDS computed via a black-box (BB) application of the reverse-mode capability of ADIFOR 3.0, and also with similar SDs calculated via the method of finite differences. 2) Second-order (SO) SDs have been implemented in the 2D ASNWERS code using the very efficient strategy that was originally proposed (but not previously tested) of Reference 3, Appendix A. Furthermore, these SO SOs have been validated for accuracy and computational efficiency. 3) Studies were conducted in Quasi-1D and 2D concerning the smoothness (or lack of smoothness) of the FO and SO SD's for flows with shock waves. The phenomenon is documented in the publications of this study (listed subsequently), however, the specific numerical mechanism which is responsible for this unsmoothness phenomenon was not discovered. 4) The FO and SO derivatives for Quasi-1D and 2D flows were applied to predict aerodynamic design uncertainties, and were also applied in robust design optimization studies.

Taylor, Arthur C., III↗
NASA NTRS1992

Automatic differentiation as a tool in engineering design

Automatic Differentiation (AD) is a tool that systematically implements the chain rule of differentiation to obtain the derivatives of functions calculated by computer programs. AD is assessed as a tool for engineering design. The forward and reverse modes of AD, their computing requirements, as well as approaches to implementing AD are discussed. The application of two different tools to two medium-size structural analysis problems to generate sensitivity information typically necessary in an optimization or design situation is also discussed. The observation is made that AD is to be preferred to finite differencing in most cases, as long as sufficient computer storage is available; in some instances, AD may be the alternative to consider in lieu of analytical sensitivity analysis.

Barthelemy, Jean-Francois↗
DOE OSTI2021

An empirical investigation of organic software product lines

Abstract Software product line engineering is a best practice for managing reuse in families of software systems that is increasingly being applied to novel and emerging domains. In this work we investigate the use of software product line engineering in one of these new domains, synthetic biology. In synthetic biology living organisms are programmed to perform new functions or improve existing functions. These programs are designed and constructed using small building blocks made out of DNA. We conjecture that there are families of products that consist of common and variable DNA parts, and we can leverage product line engineering to help synthetic biologists build, evolve, and reuse DNA parts. In this paper we perform an investigation of domain engineering that leverages an open-source repository of more than 45,000 reusable DNA parts. We show the feasibility of these new types of product line models by identifying features and related artifacts in up to 93.5% of products, and that there is indeed both commonality and variability. We then construct feature models for four commonly engineered functions leading to product lines ranging from 10 to 7.5 × 10 20 products. In a case study we demonstrate how we can use the feature models to help guide new experimentation in aspects of application engineering. Finally, in an empirical study we demonstrate the effectiveness and efficiency of automated reverse engineering on both complete and incomplete sets of products. In the process of these studies, we highlight key challenges and uncovered limitations of existing SPL techniques and tools which provide a roadmap for making SPL engineering applicable to new and emerging domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Computationally Guided Design of Multiple Impurities Tolerant Electrode (Final Report)

The current project was based on a combined experimental and computational approach, which can help recommend better cathode materials under multiple impurities conditions. The PI will mainly take in charge of experimental and computational thermodynamics of the selected cathode materials for the SOFC applications under multiple impurities. While the co-PI will run the electrochemical tests of the cathodes recommended and eventually the long-term degradation tests. At the end of the project, a multiple tolerant cathode material based on the combined experimental and computational approach will be recommended and the reliability of the commonly used accelerated testing will be evaluated. It will address multiple impurities poisoning effect of SO2, CO2, Cr and H2O on the LSM, LSCF and LNO cathodes by identifying the formation of the detrimental secondary phases by XRD, SEM and TEM techniques. And further recommended cathode material will be subjected to electrochemical testing and the most promising ones will be applied to long-term tests. The hybrid approach the PI proposed will not only be applied to the design of multiple impurities tolerant cathodes but will also be considered in the future oxygen electrode applications in SOECs or reversible SOCs. This hybrid computational and experimental approach includes four sections: 1) Investigation of single impurity poisoning on LSM, LSCF and LNO cathodes in the presence of SO2, CO2, Cr. In this section, LSM, LSCF cathodes from FuelCellMaterials and LNO cathodes from Sol-Gel synthesis will be heat-treated in the above single impurity. And the formation of the secondary phases as well as the corresponding simulations will be cross compared, which shows good agreement between each other. Meanwhile, the accelerated testing approach will be evaluated in these systems compared with the previous published work to understand the reliability of the approach. 2) Investigation of multiple impurities poisoning on LSM, LSCF and LNO cathodes in the presence of Cr+H2O, SO2+Cr and SO2+Cr+H2O conditions. We have also applied these 3 candidate cathodes under these multiple impurities’ conditions and the long-term degradation mechanism of the multiple impurities will be understood with the help of the combined experimental and computational approach. Meanwhile, the synergistic effect of those impurities will be compared with the individual ones in the same cathode system to further reveal the actual operating conditions. 3) Electrochemical testing and polarization of the recommended cathode. We have demonstrated very low polarization resistance in LSCF (core)-LSM (shell) electrodes using MSD process and analyzed impedance spectra using DRT analysis and confirming that the low polarization resistance in LSM infiltrated MSD cells is due to reduction in polarization resistance associated with O2-adsorption process. 4). Long-term degradation testing. We have developed a versatile MSD-based process to deposit various continuous coatings onto porous scaffolds and established baseline for longer term Cr-impurity testing in future projects.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Rapid estimation of frequency response functions by close-range photogrammetry

The accuracy of a rapid method which estimates the frequency response function from stereoscopic dynamic data is computed. It is shown that reversal of the order of the operations of coordinate transformation and Fourier transformation, which provides a significant increase in computational speed, introduces error. A portion of the error, proportional to the perturbation components normal to the camera focal planes, cannot be eliminated. The remaining error may be eliminated by proper scaling of frequency data prior to coordinate transformation. Methods are developed for least squares estimation of the full 3x3 frequency response matrix for a three dimensional structure.

Tripp, J. S.↗
DOE OSTI2021

Solar Energetic Proton Access to the Inner Magnetosphere During the September 7–8, 2017 Event

Here, the access of solar energetic protons into the inner magnetosphere on September 7–8, 2017 is investigated by following reversed proton trajectories to compute the proton cutoff energy using the Dartmouth geomagnetic cutoff code (Kress et al., 2010, https://doi.org/10.1029/2009sw000488). The cutoff energies for protons coming from the west and east direction, the minimum and maximum cutoff energy respectively, are calculated every 5 min along the orbit of Van Allen Probes using TS07 and the Lyon-Fedder-Mobarry (LFM) MHD magnetic field model. The result shows that the cutoff energy increases significantly as the radial distance decreases, and that the cutoff energy decreases with the building up of the ring current during magnetic storms. Solar wind dynamic pressure also affects cutoff suppression (Kress et al., 2004, https://doi.org/10.1029/2003gl018599). The LFM-RCM model shows stronger suppression of cutoff energy than TS07 during strong solar wind driving conditions. The simulation result is compared with proton flux measurements, showing consistent variation of the cutoff location during the September 7–8, 2017 geomagnetic storm.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1977

Detection of severe storms through a tropospheric-ionospheric coupling mechanism

Acoustic-gravity waves were detected by a ground-based ionospheric sounding array, and the location of the wave generation source was determined by a reverse group ray path computation. Computed sources of these waves were located near locations where tornadoes touched down from 2 to 4 hours later. It is suggested that the overshooting and ensuing collapse of convective turrets may be responsible for generating the acoustic-gravity waves observed.

Hung, R. J.↗