Search NASASearch

SEARCH · Search NASA

Results for “Time step adaptivity”

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 199 records · Page 11

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted

The Role of Adaptation in Body Load-Regulating Mechanisms During Locomotion

Body loading is a fundamental parameter that modulates motor output during locomotion, and is especially important for controlling the generation of stepping patterns, dynamic balance, and termination of locomotion. Load receptors that regulate and control posture and stance in locomotion include the Golgi tendon organs and muscle spindles at the hip, knee, and ankle joints, and the Ruffini endings and the Pacinian corpuscles in the soles of the feet. Increased body weight support (BWS) during locomotion results in an immediate reorganization of locomotor control, such as a reduction in stance and double support duration and decreased hip, ankle, and knee angles during the gait cycle. Previous studies on the effect during exposure to increased BWS while walking showed a reduction in lower limb joint angles and gait cycle timing that represents a reorganization of locomotor control. Until now, no studies have investigated how locomotor control responds after a period of exposure to adaptive modification in the body load sensing system. The goal of this research was to determine the adaptive properties of body load-regulating mechanisms in locomotor control during locomotion. We hypothesized that body load-regulating mechanisms contribute to locomotor control, and adaptive changes in these load-regulating mechanisms require reorganization to maintain forward locomotion. Head-torso coordination, lower limb movement patterns, and gait cycle timing were evaluated before and after a 30-minute adaptation session during which subjects walked on a treadmill at 5.4 km/hr with 40% body weight support (BWS). Before and after the adaptation period, head-torso and lower limb 3D kinematic data were obtained while performing a goal directed task during locomotion with 0% BWS using a video-based motion analysis system, and gait cycle timing parameters were collected by foot switches positioned under the heel and toe of the subjects shoes. Subjects showed adaptive modification in the body load-regulating mechanisms that included increased head movement amplitude, increased knee and ankle flexion, and increased stance, stride, and double support time, with no change in the performance of the task with respect to that measured before exposure to BWS. These changes in locomotor control are opposite to that reported during 40% BWS exposure and indicative of an after-effect after removal of the adaptive stimulus. Therefore, it is evident that just 30 minutes of 40% BWS during locomotion was sufficient to induce adaptive modifications in the body load sensing systems that contribute to reorganization of sensory contributions to stable locomotor control.

Ruttley, Tara

On the Application of Hilbert Spectral Analysis for Climate Studies

The Hilbert spectral analysis (Huang et al, 1998, Proceedings of the Royal Society of London, A 454, pp 903-995) consisted of two steps: First, the data has to be reduced into a finite number of Intrinsic Mode Function by the Empirical Mode Decomposition method, then the resulting Intrinsic Mode Functions are converted to time-frequency-energy distribution through Hilbert transform. In this approach, the Empirical Mode Functions served as the basis functions with which the data is expanded. This basis function is adaptive, and the decomposition is nonlinear. Furthermore, as the Hilbert transform is a singular transform, it retains a high degree of local information. The instantaneous frequency is determined by differentiation of the phase function; therefore, there is no restriction of the 'uncertainty principle' for all the time-frequency analysis resulting from a priori basis approach. With the adaptive basis and the instantaneous frequency, the Hilbert Spectral analysis can represent data from nonlinear and nonstationary processes without resorting to the harmonics. Another advantage of using instantaneous frequency is the ability to find out frequency from limited length of data, which is a critical problem in climate studies. As the processes driving the climate changes could be both nonlinear and nonstationary, the Hilbert Spectral Analysis could be of great use in examining the underlying mechanisms. A preliminary study based on the length of day data will be presented as example for the application of the Hilbert Spectral Analysis for climate study.

Huang, Norden E.

Software Reduces Radio-Interference Effects in Radar Data

A computer program suppresses the effects of narrow-band radio-frequency interference (RFI) on the data collected by a wide-band radar system. The need for this program arises because some advanced wide-band synthetic-aperture radar systems utilize frequency bands that include frequencies used by other radio services. In this program, the RFI environment is represented by an auto-regressive process, the frequency band of which is narrow relative to that of the radar. Most of the RFI signals, both narrow- and wide-band, are estimated in one pass of a least-mean-square (LMS) adaptive filter. The program implements three popular LMS algorithms: the time-domain LMS, the frequency-domain LMS, and the filter-bank LMS adaptive-filter algorithms. The program can be run in a manual or automatic mode. In the manual mode, the user selects the filter parameters prior to execution. In the automatic mode, the program utilizes median-filter and spectral-estimation techniques plus the variable-step-size LMS algorithm for automatic determination of filter parameters, and the parameters are adaptively changed as functions of the inputs, resulting in better overall performance.

Le, Charles T-C

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL

Robotic Processing Of Rocket-Engine Nozzles

Automated manufacturing cell containing computer-controlled robotic processing system developed to implement some important related steps in fabrication of rocket-engine nozzles. Performs several tedious and repetitive fabrication, measurement, adjustment, and inspection processes and subprocesses now performed manually. Offers advantages of reduced processing time, greater consistency, excellent collection of data, objective inspections, greater productivity, and simplified fixturing. Also affords flexibility: by making suitable changes in hardware and software, possible to modify process and subprocesses. Flexibility makes work cell adaptable to fabrication of heat exchangers and other items structured similarly to rocket nozzles.

Gilbert, Jeffrey L.

Evaluation of influence coefficient matrices for an adaptive-wall wind tunnel and applications to one-step convergence scheme

Within the framework of ventilated adaptive-wall wind tunnel testing, experiments were conducted to evaluate the influence coefficients from the side-wall pressure distribution, the streamwise velocity distribution, and the vertical velocity distribution on two control surfaces. Consideration is given to the effects of Mach number, superposition, suction, blowing, and the presence of a nonlifting model. It is shown that a single-step convergence scheme reduces the testing time by eliminating the intermediate experimental iterations; however, its success depends on the accuracy of the influence coefficients.

Celik, Z. Z.

Grid generation and adaptation for the Direct Simulation Monte Carlo Method

A grid generation and adaptation procedure based on the method of transfinite interpolation is incorporated into the Direct Simulation Monte Carlo Method of Bird. In addition, time is advanced based on a local criterion. The resulting procedure is used to calculate steady flows past wedges and cones. Five chemical species are considered. In general, the modifications result in a reduced computational effort. Moreover, preliminary results suggest that the simulation method is time step dependent if requirements on cell sizes are not met.

Olynick, David P.

Identification of integrated airframe-propulsion effects on an F-15 aircraft for application to drag minimization

The application of an adaptive real-time measurement-based performance optimization technique is being explored for a future flight research program. The key technical challenge of the approach is parameter identification, which uses a perturbation-search technique to identify changes in performance caused by forced oscillations of the controls. The controls on the NASA F-15 highly integrated digital electronic control (HIDEC) aircraft were perturbed using inlet cowl rotation steps at various subsonic and supersonic flight conditions to determine the effect on aircraft performance. The feasibility of the perturbation-search technique for identifying integrated airframe-propulsion system performance effects was successfully shown through flight experiments and postflight data analysis. Aircraft response and control data were analyzed postflight to identify gradients and to determine the minimum drag point. Changes in longitudinal acceleration as small as 0.004 g were measured, and absolute resolution was estimated to be 0.002 g or approximately 50 lbf of drag. Two techniques for identifying performance gradients were compared: a least-squares estimation algorithm and a modified maximum likelihood estimator algorithm. A complementary filter algorithm was used with the least squares estimator.

Schkolnik, Gerald S.

Identification of integrated airframe: Propulsion effects on an F-15 aircraft for application to drag minimization

The application of an adaptive real-time measurement-based performance optimization technique is being explored for a future flight research program. The key technical challenge of the approach is parameter identification, which uses a perturbation-search technique to identify changes in performance caused by forced oscillations of the controls. The controls on the NASA F-15 highly integrated digital electronic control (HIDEC) aircraft were perturbed using inlet cowl rotation steps at various subsonic and supersonic flight conditions to determine the effect on aircraft performance. The feasibility of the perturbation-search technique for identifying integrated airframe-propulsion system performance effects was successfully shown through flight experiments and postflight data analysis. Aircraft response and control data were analyzed postflight to identify gradients and to determine the minimum drag point. Changes in longitudinal acceleration as small as 0.004 g were measured, and absolute resolution was estimated to be 0.002 g or approximately 50 lbf of drag. Two techniques for identifying performance gradients were compared: a least-squares estimation algorithm and a modified maximum likelihood estimator algorithm. A complementary filter algorithm was used with the least squares estimator.

Schkolnik, Gerard S.

A General Method for Solving Systems of Non-Linear Equations

The method of steepest descent is modified so that accelerated convergence is achieved near a root. It is assumed that the function of interest can be approximated near a root by a quadratic form. An eigenvector of the quadratic form is found by evaluating the function and its gradient at an arbitrary point and another suitably selected point. The terminal point of the eigenvector is chosen to lie on the line segment joining the two points. The terminal point found lies on an axis of the quadratic form. The selection of a suitable step size at this point leads directly to the root in the direction of steepest descent in a single step. Newton's root finding method not infrequently diverges if the starting point is far from the root. However, the current method in these regions merely reverts to the method of steepest descent with an adaptive step size. The current method's performance should match that of the Levenberg-Marquardt root finding method since they both share the ability to converge from a starting point far from the root and both exhibit quadratic convergence near a root. The Levenberg-Marquardt method requires storage for coefficients of linear equations. The current method which does not require the solution of linear equations requires more time for additional function and gradient evaluations. The classic trade off of time for space separates the two methods.

Nachtsheim, Philip R.

The Optoretinogram at 38

As early as 1980, Donald MacLeod coined the term "optoretinogram" to describe a change in the near-infrared reflectance of the retina following visual stimulation. The hope that such signals might exist was stimulated by an observation of a change in the infrared transmission of isolated toad retinae. This signal showed a time course of several seconds, much slower than the electroretinogram, and was thought to be related to some chemical step in the phototransduction cascade. During the early 1980's, we made several attempts to observe such signals in vivo in humans, without success. In the intervening decades, however, researchers have succeeded in observing what are now known as "intrinsic optical signals" or IOSs. Responses from single photoreceptors (obtained using adaptive optics) are surprisingly heterogenous, with some cells increasing their reflectance in response to visual stimulation, while others decrease. This talk will review our early work, survey recent findings, and consider various theories for the origin of the effect.

noninvasive monitoring

Adaptive-wall wind-tunnel development for transonic testing

Experimental techniques for rapid assessment and correction of wall interference in an adaptive-wall wind tunnel are described. The experimental arrangement allows laser velocimetry measurements on two control surfaces and incorporates a dedicated computer for data processing. The apparatus and its instrumentation are described, and typical results from an experiment on a nonlifting NACA 0012 airfoil at M = 0.78 are discussed. It is concluded that the time to acquire laser Doppler velocimeter data should be decreased, and the possibility of using one-step algorithms should be investigated.

Satyanarayana, B.

LightForce Photon-Pressure Collision Avoidance: Updated Efficiency Analysis Utilizing a Highly Parallel Simulation Approach

This paper provides an updated efficiency analysis of the LightForce space debris collision avoidance scheme. LightForce aims to prevent collisions on warning by utilizing photon pressure from ground based, commercial off the shelf lasers. Past research has shown that a few ground-based systems consisting of 10 kilowatt class lasers directed by 1.5 meter telescopes with adaptive optics could lower the expected number of collisions in Low Earth Orbit (LEO) by an order of magnitude. Our simulation approach utilizes the entire Two Line Element (TLE) catalogue in LEO for a given day as initial input. Least-squares fitting of a TLE time series is used for an improved orbit estimate. We then calculate the probability of collision for all LEO objects in the catalogue for a time step of the simulation. The conjunctions that exceed a threshold probability of collision are then engaged by a simulated network of laser ground stations. After those engagements, the perturbed orbits are used to re-assess the probability of collision and evaluate the efficiency of the system. This paper describes new simulations with three updated aspects: 1) By utilizing a highly parallel simulation approach employing hundreds of processors, we have extended our analysis to a much broader dataset. The simulation time is extended to one year. 2) We analyze not only the efficiency of LightForce on conjunctions that naturally occur, but also take into account conjunctions caused by orbit perturbations due to LightForce engagements. 3) We use a new simulation approach that is regularly updating the LightForce engagement strategy, as it would be during actual operations. In this paper we present our simulation approach to parallelize the efficiency analysis, its computational performance and the resulting expected efficiency of the LightForce collision avoidance system. Results indicate that utilizing a network of four LightForce stations with 20 kilowatt lasers, 85% of all conjunctions with a probability of collision Pc > 10 (sup -6) can be mitigated.

conjunction analysis

Estimating an executive summary of a time series: the tendency

In this paper, we revisit the problem of decomposing a signal into a tendency and a residual. The tendency describes an executive summary of a signal that encapsulates its notable characteristics while disregarding seemingly random, less interesting aspects. Building upon the Intrinsic Time Decomposition (ITD) and information-theoretical analysis, we introduce two alternative procedures for selecting the tendency from the ITD baselines. The first is based on the maximum extrema prominence, namely the maximum difference between extrema within each baseline. Specifically this method selects the tendency as the baseline from which an ITD step would produce the largest decline of the maximum prominence. The second method uses the rotations from the ITD and selects the tendency as the last baseline for which the associated rotation is statistically stationary. We delve into a comparative analysis of the information content and interpretability of the tendencies obtained by our proposed methods and those obtained through conventional low-pass filtering schemes, particularly the Hodrik–Prescott (HP) filter. Our findings underscore a fundamental distinction in the nature and interpretability of these tendencies, highlighting their context-dependent utility with emphasis in multi-scale signals. Through a series of real-world applications, we demonstrate the computational robustness and practical utility of our proposed tendencies, emphasizing their adaptability and relevance in diverse time series contexts.

Time series analysis

Adaptive temperature profile control of a multizone crystal growth furnace

An intelligent measurement system is described which is used to assess the shape of a crystal while it is growing inside a multizone transparent furnace. A color video imaging system observes the crystal in real time, and determines the position and the shape of the interface. This information is used to evaluate the crystal growth rate, and to analyze the effects of translational velocity and temperature profiles on the shape of the interface. Creation of this knowledge base is the first step to incorporate image processing into furnace control.

Batur, C.

Adaptive Graph Learning with Transformer for Multi-Reservoir Inflow Prediction

Reservoir inflow prediction is crucial for water resource management, yet existing approaches mainly focus on single-reservoir models that ignore spatial dependencies among interconnected reservoirs. We introduce AdaTrip as an adaptive, time-varying graph learning framework for multi-reservoir inflow forecasting. AdaTrip constructs dynamic graphs where reservoirs are nodes with directed edges reflecting hydrological connections, employing attention mechanisms to automatically identify crucial spatial and temporal dependencies. Evaluation on thirty reservoirs in the Upper Colorado River Basin demonstrates superiority over existing baselines, with improved performance for reservoirs with limited records through parameter sharing. Additionally, AdaTrip provides interpretable attention maps at edge and time-step levels, offering insights into hydrological controls to support operational decision-making. Our code is available at https://github.com/humphreyhuu/AdaTrip.

Hu, Pengfei [ORNL] (ORCID:0009000367130950)

Liquid Crystal on Silicon Wavefront Corrector

A low cost, high resolution, liquid crystal on silicon, spatial light modulator has been developed for the correction of huge aberrations in an optical system where the polarization dependence and the chromatic nature are tolerated. However, the overall system performance suggests that this device is also suitable for real time correction of aberration in human eyes. This device has a resolution of 1024 x 768, and is driven by an XGA display driver. The effective stroke length of the device is 700 nm and 2000 nm for the visible and IR regions of the device, respectively. The response speeds are 50 Hz and 5 Hz, respectively, which are fast enough for real time adaptive optics for aberrations in human eyes. By modulating a wavefront of 2 pi, this device can correct for arbitrary high order wavefront aberrations since the 2-D pixel array is independently controlled by the driver. The high resolution and high accuracy of the device allow for diffraction limited correction of the tip and tilt or defocus without an additional correction loop. We have shown that for every wave of aberration, an 8 step blazed grating is required to achieve high diffraction efficiency around 80%. In light of this, up to 125 waves peak to valley of tip and tilt can be corrected if we choose the simplest aberration. Corrections of 34 waves of aberration, including high order Zernicke terms in a high magnification telescope, to diffraction limited performance (residual wavefront aberration less than 1/30 lambda at 632.8 nm) have been observed at high efficiency.

Pouch, John