SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 1,405 records · Page 78

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Joint transform correlation with the synthetic estimation filter

The synthetic estimation filter (SEF) combines a set of input images into a much smaller set of filters. Those input images correspond to views of a known object taken with a variety of values of the pose parameter to be estimated. The response of the filter to variations in the pose parameter is tailored so that estimation can be done to a precision considerably finer than the intervals between filters in the pose space. The technique is extended to the joint transform correlator; as a reference, the composite image from which the SEF is conventionally calculated is used. Precautions are necessary for practical reasons. For example, the SEF algorithm may call for negative values in the image, an unrealizable condition with amplitude encoding, whereas with phase-mostly filtering in the SEF, negative image values are easily accommodated.

Monroe, Stanley E., Jr.↗

NASA NTRS1992

Dynamic interactions between hypersonic vehicle aerodynamics and propulsion system performance

Described here is the development of a flexible simulation model for scramjet hypersonic propulsion systems. The primary goal is determination of sensitivity of the thrust vector and other system parameters to angle of attack changes of the vehicle. Such information is crucial in design and analysis of control system performance for hypersonic vehicles. The code is also intended to be a key element in carrying out dynamic interaction studies involving the influence of vehicle vibrations on propulsion system/control system coupling and flight stability. Simple models are employed to represent the various processes comprising the propulsion system. A method of characteristics (MOC) approach is used to solve the forebody and external nozzle flow fields. This results in a very fast computational algorithm capable of carrying out the vast number of simulation computations needed in guidance, stability, and control studies. The three-dimensional fore- and aft body (nozzle) geometry is characterized by the centerline profiles as represented by a series of coordinate points and body cross-section curvature. The engine module geometry is represented by an adjustable vertical grid to accommodate variations of the field parameters throughout the inlet and combustor. The scramjet inlet is modeled as a two-dimensional supersonic flow containing adjustable sidewall wedges and multiple fuel injection struts. The inlet geometry including the sidewall wedge angles, the number of injection struts, their sweepback relative to the vehicle reference line, and strut cross-section are user selectable. Combustion is currently represented by a Rayleigh line calculation including corrections for variable gas properties; improved models are being developed for this important element of the propulsion flow field. The program generates (1) variation of thrust magnitude and direction with angle of attack, (2) pitching moment and line of action of the thrust vector, (3) pressure and temperature distributions throughout the system, and (4) performance parameters such as thrust coefficient, specific impulse, mass flow rates, and equivalence ratio. Preliminary results are in good agreement with available performance data for systems resembling the NASP vehicle configuration.

Flandro, G. A.↗

NASA NTRS1988

An adaptive identification and control scheme for large space structures

A unified identification and control scheme capable of achieving space at form performance objectives under nominal or failure conditions is described. Preliminary results are also presented, showing that the methodology offers much promise for effective robust control of large space structures. The control method is a multivariable, adaptive, output predictive controller called Model Predictive Control (MPC). MPC uses a state space model and input reference trajectories of set or tracking points to adaptively generate optimum commands. For a fixed model, MPC processes commands with great efficiency, and is also highly robust. A key feature of MPC is its ability to control either nonminimum phase or open loop unstable systems. As an output controller, MPC does not explicitly require full state feedback, as do most multivariable (e.g., Linear Quadratic) methods. Its features are very useful in LSS operations, as they allow non-collocated actuators and sensors. The identification scheme is based on canonical variate analysis (CVA) of input and output data. The CVA technique is particularly suited for the measurement and identification of structural dynamic processes - that is, unsteady transient or dynamically interacting processes such as between aerodynamics and structural deformation - from short, noisy data. CVA is structured so that the identification can be done in real or near real time, using computationally stable algorithms. Modeling LSS dynamics in 1-g laboratories has always been a major impediment not only to understanding their behavior in orbit, but also to controlling it. In cases where the theoretical model is not confirmed, current methods provide few clues concerning additional dynamical relationships that are not included in the theoretical models. CVA needs no a priori model data, or structure; all statistically significant dynamical states are determined using natural, entropy-based methods. Heretofore, a major limitation in applying adaptive controllers to LSS problems has been the requirement for accurate models. We feel that MPC coupled with the CVA identification scheme will resolve this limitation. The results are based on research performed under an SBIR grant, administered by the Air Force Office of Scientific Research.

Carroll, J. V.↗

NASA NTRS1995

Investigation of Convection and Pressure Treatment with Splitting Techniques

Treatment of convective and pressure fluxes in the Euler and Navier-Stokes equations using splitting formulas for convective velocity and pressure is investigated. Two schemes - controlled variation scheme (CVS) and advection upstream splitting method (AUSM) - are explored for their accuracy in resolving sharp gradients in flows involving moving or reflecting shock waves as well as a one-dimensional combusting flow with a strong heat release source term. For two-dimensional compressible flow computations, these two schemes are implemented in one of the pressure-based algorithms, whose very basis is the separate treatment of convective and pressure fluxes. For the convective fluxes in the momentum equations as well as the estimation of mass fluxes in the pressure correction equation (which is derived from the momentum and continuity equations) of the present algorithm, both first- and second-order (with minmod limiter) flux estimations are employed. Some issues resulting from the conventional use in pressure-based methods of a staggered grid, for the location of velocity components and pressure, are also addressed. Using the second-order fluxes, both CVS and AUSM type schemes exhibit sharp resolution. Overall, the combination of upwinding and splitting for the convective and pressure fluxes separately exhibits robust performance for a variety of flows and is particularly amenable for adoption in pressure-based methods.

Thakur, Siddharth↗

NASA NTRS2010

Investigating the Performance of the Jason-2/OSTM Radar Altimeter Over Lakes and Reservoirs

Many inland water investigations utilize archival and near-real time radar altimetry data to enable observation of the variation in surface water level. A multi-altimeter approach allows a more global outlook with improved spatial resolution, and combined long-term observations improve statistical analyses. Central to all programs is a performance assessment of each instrument. Here, we focus on data quantity and quality pertaining to the Poseidon-3 radar altimeter onboard the Jason-2/OSTM satellite.Utilizing an interim data set (IGDR), studies show that the new on-board DIODE/median and DIODE/DEM tracking modes are performing well, acquiring and maintaining the majority of lake and reservoir surfaces in varying terrains. The 20-Hz along-track resolution of the data, and particularly the availability of the range output from the ice-retracker algorithm, also improves the number of valid height measurements. Based on test-case lakes and reservoirs, output from the ice-retracker algorithm is also seen to have a clear advantage over the ocean-retracker having better height stability across calm and icy surfaces, a greater ability to gain coastline waters, and less sensitivity to loss of water surface when there is island contamination in the radar echo. Such on-board tracking and postprocessing retracking enables the lake waters to be quickly gained after coastline crossing. Values can range from <0.1 s to 2.5 s, but the majority of measurements are obtained in less than 0.4 s or <2.3 km from the coast. Validation exercises reveal that targets of 150 km2 surface area and 0.8 km width are able to be monitored offering greater potential to acquire lakes in the 100 C300 km2 size-category. Time series of height variations are also found to be accurate to 3 to 33 cm rms depending on target size and the presence of winter ice. These findings are an improvement over the IGDR/GDR results from the predecessor Jason-1 and TOPEX/Poseidon missions and can satisfy the accuracy requirements of both the science-related and operational lake study programs.

Birkett, C. M.↗

DOE OSTI2022

Beyond carbon flux partitioning: Carbon allocation and nonstructural carbon dynamics inferred from continuous fluxes

Carbon (C) allocation and nonstructural carbon (NSC) dynamics play essential roles in plant growth and survival under stress and disturbance. However, quantitative understanding of these processes remains limited. Here, in this work, we propose a framework where we connect commonly measured carbon cycle components (eddy covariance fluxes of canopy CO 2 exchange, soil CO 2 efflux, and allometry-based biomass and net primary production) by a simple mass balance model to derive ecosystem-level NSC dynamics (NSC i ), C translocation (dC i ), and the biomass production efficiency (BPE i ) in above- and belowground plant (i = agp and bgp) compartments. We applied this framework to two long-term monitored loblolly pine (Pinus taeda) plantations of different ages in North Carolina and characterized the variations of NSC and allocation in years under normal and drought conditions. The results indicated that the young stand did not have net NSC flux at the annual scale, whereas the mature stand stored a near-constant proportion of new assimilates as NSC every year under normal conditions, which was comparable in magnitude to new structural growth. Roots consumed NSC in drought and stored a significant amount of NSC post drought. The above- and belowground dC i and BPE i varied more from year to year in the young stand and approached a relatively stable pattern in the mature stand. The belowground BPE bgp differed the most between the young and mature stands and was most responsive to drought. With the internal C dynamics quantified, this framework may also improve biomass production estimation, which reveals the variations resulting from droughts. Overall, these quantified ecosystem-scale dynamics were consistent with existing evidence from tree-based manipulative experiments and measurements and demonstrated that combining the continuous fluxes as proposed here can provide additional information about plant internal C dynamics. Given that it is based on broadly available flux data, the proposed framework is promising to improve the allocation algorithms in ecosystem C cycle models and offers new insights into observed variability in soil–plant–climate interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Toward implementing autonomous adaptive data acquisition for scanning hyperspectral imaging of biological systems

Autonomous experimentation is an emerging area of research, primarily related to autonomous vehicles, scientific combinatorial discovery approaches in materials science and drug discovery, and iterative research loops of planning, experimentation, and analysis. However, autonomous approaches developed in these contexts are difficult to apply to high-dimensional mapping technologies, such as scanning hyperspectral imaging of biological systems, due to sample complexity and heterogeneity. We briefly cover the history of adaptive sampling algorithms and surrogate modeling in order to define autonomous adaptive data acquisition as an objective-based, flexible building block for future biological imaging experimentation driven by intelligent infrastructure. We subsequently summarize the recent implementations of autonomous adaptive data acquisition (AADA) for scanning hyperspectral imaging, assess how these address the difficulties of autonomous approaches in hyperspectral imaging, and highlight the AADA design variation from a goal-oriented perspective. Finally, we present a modular AADA architecture that embeds AADA-driven flexible building blocks to address the challenge of time resolution for high-dimensional scanning hyperspectral imaging of nonequilibrium dynamical systems. In our example research-driven experimental design case, we propose an AADA infrastructure for time-resolved, noninvasive, and label-free scanning hyperspectral imaging of living biological systems. This AADA infrastructure can accurately target the correct state of the system for experimental workflows that utilize subsequent expensive, high-information-content analytical techniques.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1977

A first-principles derivation of Doppler noise expected from solar wind density fluctuations

The level of Doppler noise (DN) expected from solar wind (SW) density fluctuations (DF) is derived beginning with the expression for refractive index variations. The calculation takes account of up and downlink paths and of the method actually used to produce the DN values. The usual assumptions that the DF are frozen in, that the large-scale radial variation can be separated from the DF, that the DF power spectrum is a power law with outer scale k sub 0 and that the DF are homogeneous on scales less than 2c delta t, delta t = sample time, are made. The result agrees quite well with previous observations of DN. Corrections for the finite number of points used in the actual algorithm are discussed.

Callahan, P. S.↗

NASA NTRS1981

An algorithm for minimum-cost set-point ordering in a cryogenic wind tunnel

An algorithm for minimum cost ordering of set points in a cryogenic wind tunnel is developed. The procedure generates a matrix of dynamic state transition costs, which is evaluated by means of a single-volume lumped model of the cryogenic wind tunnel and the use of some idealized minimum-costs, which is evaluated by means of a single-volume lumped model of the cryogenic wind tunnel and the use of some idealized minimum-cost state-transition control strategies. A branch and bound algorithm is employed to determine the least costly sequence of state transitions from the transition-cost matrix. Some numerical results based on data for the National Transonic Facility are presented which show a strong preference for state transitions that consume to coolant. Results also show that the choice of the terminal set point in an open odering can produce a wide variation in total cost.

Tripp, J. S.↗

NASA NTRS1986

On symmetric and upwind TVD schemes

The performance of the upwind and symmetric total variation diminishing (TVD) schemes in viscous and inviscid airfoil steady-state calculations is considered, and the extension of the implicit second-order-accurate TVD scheme for hyperbolic systems of conservative laws in curvilinear coordinates is discussed. For two-dimensional steady-state applications, schemes are implemented in a conservative noniterative alternating direction implicit form, and results illustrate that the algorithm produces a fairly good solution for an RAE2822 airfoil calculation. The study demonstrates that the symmetric TVD scheme is as accurate as the upwind TVD scheme, while requiring less computational effort than it.

Yee, H. C.↗

NASA NTRS1989

Determination of design and operation parameters for upper atmospheric research instrumentation to yield optimum resolution with deconvolution, appendix 4

The power spectrum for a stationary random process can be defined with the Wiener-Khintchine Theorem, which says that the power spectrum and the auto correlation function are a Fourier transform pair. To implement this theorem for signals that are discrete and of finite length we can use the Blackman-Tukey method. Blackman and Tukey (1958) show that a function w(tau), called a lag window, can be applied to the auto correlation estimates to obtain power spectrum estimates that are statistically stable. The Fourier transform of w(r) is called a spectral window. Typical choices for spectral windows show a distinct trade-off between the main lobe width and side lobe strength. A new idea for designing windows by taking linear combinations of the standard windows to produce hybrid windows was introduced by Smith (1985). We implement Smith's idea to obtain spectral windows with narrow main lobes and smaller (compared with typical windows) near side lobes. One of the main contributions of this thesis is that we show that Smith's problem is equivalent to a Quadratic Programming (QP) problem with linear equality and inequality constraints. A computer program was written to produce hybrid windows by setting up and solving the QP problem. We also developed and solved two variations of the original problem. The two variations involved changing the inequality constraints in both cases from non negativity on the combination coefficients to non negativity on the hybrid lag window itself. For the second variation, the window functions used to construct the hybrid window were changed to a frequency-variable set of truncated cosinusoids. A series of tests was run with the three computer programs to investigate the behavior of the hybrid spectral and lag windows. Emphasis was put on obtaining spectral windows with both relatively narrow main lobes and the lowest possible (for these algorithms) near side lobes. Some success was achieved for this goal. A 10 dB peak side lobe reduction over the rectangular spectral window without significant main lobe broadening was achieved. Also, average side lobe levels of -117 dB were reached at a cost of doubling the main lobe width (at the -3 dB point).

Ioup, George E.↗

NASA NTRS2010

Detection of Multi-Layer and Vertically-Extended Clouds Using A-Train Sensors

The detection of mUltiple cloud layers using satellite observations is important for retrieval algorithms as well as climate applications. In this paper, we describe a relatively simple algorithm to detect multiple cloud layers and distinguish them from vertically-extended clouds. The algorithm can be applied to coincident passive sensors that derive both cloud-top pressure from the thermal infrared observations and an estimate of solar photon pathlength from UV, visible, or near-IR measurements. Here, we use data from the A-train afternoon constellation of satellites: cloud-top pressure, cloud optical thickness, the multi-layer flag from the Aqua MODerate-resolution Imaging Spectroradiometer (MODIS) and the optical centroid cloud pressure from the Aura Ozone Monitoring Instrument (OMI). For the first time, we use data from the CloudSat radar to evaluate the results of a multi-layer cloud detection scheme. The cloud classification algorithms applied with different passive sensor configurations compare well with each other as well as with data from CloudSat. We compute monthly mean fractions of pixels containing multi-layer and vertically-extended clouds for January and July 2007 at the OMI spatial resolution (l2kmx24km at nadir) and at the 5kmx5km MODIS resolution used for infrared cloud retrievals. There are seasonal variations in the spatial distribution of the different cloud types. The fraction of cloudy pixels containing distinct multi-layer cloud is a strong function of the pixel size. Globally averaged, these fractions are approximately 20% and 10% for OMI and MODIS, respectively. These fractions may be significantly higher or lower depending upon location. There is a much smaller resolution dependence for fractions of pixels containing vertically-extended clouds (approx.20% for OMI and slightly less for MODIS globally), suggesting larger spatial scales for these clouds. We also find higher fractions of vertically-extended clouds over land as compared with ocean, particularly in the tropics and summer hemisphere.

Joiner, J.↗

NASA NTRS2015

Composite Characterization Using Laser Doppler Vibrometry and Multi-Frequency Wavenumber Analysis

NASA has recognized the need for better characterization of composite materials to support advances in aeronautics and the next generation of space exploration vehicles. An area of related research is the evaluation of impact induced delaminations. Presented is a non-contact method of measuring the ply depth of impact delamination damage in a composite through use of a Scanning Laser Doppler Vibrometer (SLDV), multi-frequency wavenumber analysis, and a wavenumber-ply correlation algorithm. A single acquisition of a chirp excited lamb wavefield in an impacted composite is post-processed into a numerous single frequency excitation wavefields through a deconvolution process. A spatially windowed wavenumber analysis then extracts local wavenumbers from the wavefield, which are then correlated to theoretical dispersion curves for ply depth determination. SLDV based methods to characterize as-manufactured composite variation using wavefield analysis will also be discussed.

Juarez, Peter↗

DOE OSTI2026

Statorless mixed flow turbine for transonic pulsating inflow

Effective harvesting of power from high speed highly transient inflows such as the outflow of rotating detonation combustors (RDCs) is key to achieving their promised cycle efficiency step jump. To increase power density and efficiency simultaneously, a concept that can directly ingest transonic outflow, addressing the choking is needed, i.e. without additional transition elements. A new statorless design was assessed using a comprehensive approach to quantify all contributions to loss generation in transient flows, locally and globally. The turbine rotor was designed under steady flow conditions with a genetic algorithm. The comparison of power and loss generation between steady and unsteady flow results shows that a design methodology under steady-state conditions is suitable to characterize the performance of different designs. Three design families were further assessed under highly transient transonic conditions including traveling shock waves with a relative total pressure amplitude of 149.6% around the mean value and inflow angle variations from -26.5 deg to +51.8 deg. The oblique shock impinging on the pressure side (PS) of the turbine augments the shaft power extraction. When the oblique shock reflects between the pressure side and suction side (SS), its strength diminishes. This paper provides design guidelines on efficient turbine work extraction from the shocks emanating from detonation combustors.

rotating detonation engines↗

DOE OSTI2023

Online data-driven changepoint detection for high-dimensional dynamical systems

In this study, the detection of anomalies or transitions in complex dynamical systems is of critical importance to various applications. In this study, we propose the use of machine learning to detect changepoints for high-dimensional dynamical systems. Here, changepoints indicate instances in time when the underlying dynamical system has a fundamentally different characteristic—which may be due to a change in the model parameters or due to intermittent phenomena arising from the same model. We propose two complementary approaches to achieve this, with the first devised using arguments from probabilistic unsupervised learning and the latter devised using supervised deep learning. To accelerate the deployment of transition detection algorithms in high-dimensional dynamical systems, we introduce dimensionality reduction techniques. Our experiments demonstrate that transitions can be detected efficiently, in real-time, for the two-dimensional forced Kolmogorov flow and the Rössler dynamical system, which are characterized by anomalous regimes in phase space where dynamics are perturbed off the attractor at potentially uneven intervals. Finally, we also demonstrate how variations in the frequency of detected changepoints may be utilized to detect a significant modification to the underlying model parameters by utilizing the Lorenz-63 dynamical system.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Application of machine learning approaches in the analysis of mass absorption cross-section of black carbon aerosols: Aerosol composition dependencies and sensitivity analyses

Physics-based models typically require an in-depth understanding of a phenomenon and assumptions of the underlying process(es), which are often hard to obtain in practice, whereas data-driven machine learning models learn the structure and patterns in the training data without any prior theoretical assumptions and then use inference to develop useful predictions. A novel machine learning-based algorithm has been previously developed for the prediction of black carbon mass absorption cross-section (MAC BC ) and applied to a variety of different atmospheric environments. In contrast to light scattering theories which require assumptions about the underlying physics, this algorithm uses time-series data of aerosol properties to estimate the temporally-varying MAC BC at 870 nm. Here, we analyze our algorithm and discuss the influence of aerosol optical properties (such as Ångström exponents and single scattering albedo) and chemical composition on the model outputs and the associated accuracy. Additionally, we conduct sensitivity analyses on our models to understand how the predictions change in response to different sets of input variables. Our support vector machine (SVM) for regression model is the least sensitive to variations in the input variables, although all models tend to exhibit a degradation to their accuracy when scattering Ångström exponents are less than one.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Global Characteristics of Observable Foreshocks for Large Earthquakes

Abstract Foreshocks are the only currently widely identified precursory seismic behavior, yet their utility and even identifiability are problematic, in part because of extreme variation in behavior. Here, we establish some global trends that help identify the expected frequency of foreshocks as well the type of earthquake most prone to foreshocks. We establish these tendencies using the global earthquake catalog of the U.S. Geological Survey National Earthquake Information Center with a completeness level of magnitude 5 and mainshocks with Mw≥7.0. Foreshocks are identified using three clustering algorithms to address the challenge of distinguishing foreshocks from background activity. The methods give a range of 15%–43% of large mainshocks having at least one foreshock but a narrower range of 13%–26% having at least one foreshock with magnitude within two units of the mainshock magnitude. These observed global foreshock rates are similar to regional values for a completeness level of magnitude 3 using the same detection conditions. The foreshock sequences have distinctive characteristics with the global composite population b-values being lower for foreshocks than for aftershocks, an attribute that is also manifested in synthetic catalogs computed by epidemic-type aftershock sequences, which intrinsically involves only cascading processes. Focal mechanism similarity of foreshocks relative to mainshocks is more pronounced than for aftershocks. Despite these distinguishing characteristics of foreshock sequences, the conditions that promote high foreshock productivity are similar to those that promote high aftershock productivity. For instance, a modestly higher percentage of interplate mainshocks have foreshocks than intraplate mainshocks, and reverse faulting events slightly more commonly have foreshocks than normal or strike-slip-faulting mainshocks. The western circum-Pacific is prone to having slightly more foreshock activity than the eastern circum-Pacific.

Geochemistry & Geophysics↗

NASA NTRS1988

Summary of along-track data from the earth radiation budget satellite for several representative ocean regions

For several days in January and August 1985, the Earth Radiation Budget Satellite, a component of the Earth Radiation Budget Experiment (ERBE), was operated in an along-track scanning mode. A survey of radiance measurements taken in this mode is given for five ocean regions: the north and south Atlantic, the Arabian Sea, the western Pacific north of the Equator, and part of the Intertropical Convergence Zone. Each overflight contains information about the clear scene and three cloud categories: partly cloudy, mostly cloudy, and overcast. The data presented include the variation of longwave and shortwave radiance in each scene classification as a function of viewing zenity angle during each overflight of one of the five target regions. Several features of interest in the development of anisotropic models are evident, including the azimuthal dependence of shortwave radiance that is an essential feature of shortwave bidirectional models. The data also demonstrate that the scene classification algorithm employed by the ERBE results in scene classifications that are a function of viewing geometry.

Brooks, David R.↗