Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear reduction”

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 451 records · Page 25

Role of cross-field current instability in substorm onsets and intensifications

A cross field current instability is investigated as a potential mechanism for current reduction/disruption during substorm onsets and intensifications. Linear stability analysis shows that, for sufficiently strong current density, the instability can occur in the plasma sheet with a growth time comparable to the substorm onset time and excite waves with a significant electromagnetic component. Nonlinear analysis shows that the wave growth reaches the nonlinear stage in less than one ion gyroperiod and can reduce the cross tail current by approximately 15 to 28 of its initial value at saturation. The resulting anomalous resistivity is 11 to 12 orders of magnitude above the classical value. For a typical current reduction, the plasma in the disruption region is subjected to an earthward force. The substorm development scenario constructed based on this instability can readily account for a large number of substorm features.

Lui, Anthony T. Y.↗

A Scalable Space-Time Domain Decomposition Approach for Solving Large Scale Nonlinear Regularized Inverse Ill Posed Problems in 4D Variational Data Assimilation

We address the development of innovative algorithms designed to solve the strong-constraint Four Dimensional Variational Data Assimilation (4DVar DA) problems in large scale applications. We present a space-time decomposition approach which employs the whole domain decomposition, i.e. both along the spacial and temporal direction in the overlapping case, and the partitioning of both the solution and the operator. Starting from the global functional defined on the entire domain, we get to a sort of regularized local functionals on the set of sub domains providing the order reduction of both the predictive and the Data Assimilation models. The algorithm convergence is developed. Performance in terms of reduction of time complexity and algorithmic scalability is discussed on the Shallow Water Equations on the sphere. The number of state variables in the model, the number of observations in an assimilation cycle, as well as numerical parameters as the discretization step in time and in space domain are defined on the basis of discretization grid used by data available at repository Ocean Synthesis/Reanalysis Directory of Hamburg University.

97 MATHEMATICS AND COMPUTING↗

Analysis of Particle Image Velocimetry (PIV) Data for Acoustic Velocity Measurements

Acoustic velocity measurements were taken using Particle Image Velocimetry (PIV) in a Normal Incidence Tube configuration at various frequency, phase, and amplitude levels. This report presents the results of the PIV analysis and data reduction portions of the test and details the processing that was done. Estimates of lower measurement sensitivity levels were determined based on PIV image quality, correlation, and noise level parameters used in the test. Comparison of measurements with linear acoustic theory are presented. The onset of nonlinear, harmonic frequency acoustic levels were also studied for various decibel and frequency levels ranging from 90 to 132 dB and 500 to 3000 Hz, respectively.

Blackshire, James L.↗

GIFTS SM EDU Level 1B Algorithms

The Geosynchronous Imaging Fourier Transform Spectrometer (GIFTS) SensorModule (SM) Engineering Demonstration Unit (EDU) is a high resolution spectral imager designed to measure infrared (IR) radiances using a Fourier transform spectrometer (FTS). The GIFTS instrument employs three focal plane arrays (FPAs), which gather measurements across the long-wave IR (LWIR), short/mid-wave IR (SMWIR), and visible spectral bands. The raw interferogram measurements are radiometrically and spectrally calibrated to produce radiance spectra, which are further processed to obtain atmospheric profiles via retrieval algorithms. This paper describes the GIFTS SM EDU Level 1B algorithms involved in the calibration. The GIFTS Level 1B calibration procedures can be subdivided into four blocks. In the first block, the measured raw interferograms are first corrected for the detector nonlinearity distortion, followed by the complex filtering and decimation procedure. In the second block, a phase correction algorithm is applied to the filtered and decimated complex interferograms. The resulting imaginary part of the spectrum contains only the noise component of the uncorrected spectrum. Additional random noise reduction can be accomplished by applying a spectral smoothing routine to the phase-corrected spectrum. The phase correction and spectral smoothing operations are performed on a set of interferogram scans for both ambient and hot blackbody references. To continue with the calibration, we compute the spectral responsivity based on the previous results, from which, the calibrated ambient blackbody (ABB), hot blackbody (HBB), and scene spectra can be obtained. We now can estimate the noise equivalent spectral radiance (NESR) from the calibrated ABB and HBB spectra. The correction schemes that compensate for the fore-optics offsets and off-axis effects are also implemented. In the third block, we developed an efficient method of generating pixel performance assessments. In addition, a random pixel selection scheme is designed based on the pixel performance evaluation. Finally, in the fourth block, the single pixel algorithms are applied to the entire FPA.

Tian, Jialin↗

Development of an hp-version finite element method for computational optimal control

The purpose of this research effort was to begin the study of the application of hp-version finite elements to the numerical solution of optimal control problems. Under NAG-939, the hybrid MACSYMA/FORTRAN code GENCODE was developed which utilized h-version finite elements to successfully approximate solutions to a wide class of optimal control problems. In that code the means for improvement of the solution was the refinement of the time-discretization mesh. With the extension to hp-version finite elements, the degrees of freedom include both nodal values and extra interior values associated with the unknown states, co-states, and controls, the number of which depends on the order of the shape functions in each element. One possible drawback is the increased computational effort within each element required in implementing hp-version finite elements. We are trying to determine whether this computational effort is sufficiently offset by the reduction in the number of time elements used and improved Newton-Raphson convergence so as to be useful in solving optimal control problems in real time. Because certain of the element interior unknowns can be eliminated at the element level by solving a small set of nonlinear algebraic equations in which the nodal values are taken as given, the scheme may turn out to be especially powerful in a parallel computing environment. A different processor could be assigned to each element. The number of processors, strictly speaking, is not required to be any larger than the number of sub-regions which are free of discontinuities of any kind.

Hodges, Dewey H.↗

Combustion Stability Evaluation of Artificial Disturbance Dynamic Response

General guidelines are provided in CPIA Publication 6551 for evaluation of the response due to an artificial disturbance, however the guideline also allows for ambiguous interpretation. Stability rating devices that produce an artificial disturbance are traditionally explosive bombs and pulse guns. This paper reviews a recently developed objective process that can be applied consistently in the reduction of artificial disturbance dynamic data. It also examines three methods of response evaluation. The first method examines the response of a specific mode of interest and requires data filtering encompassing that mode. The second method examines the response of a specific mode of interest and its nonlinear components and requires a more complex filtering scheme. The third method examines the response of the entire dynamic system and consists of examining a wide bandwidth consisting of multiple modes of interest. The evaluation process is described and the advantages and disadvantages of the evaluation methods are discussed. Signal processing is used as a tool in quantifying the assessment, clearly as an improvement from the subjective heritage approach consisting primarily of engineering judgement. Data for several engines and components have been compiled and evaluated using these methods. A summary of these combustion devices is provided and observations are discussed.

Casiano, M. J.↗

Parallel interior-point solver for block-structured nonlinear programs on SIMD/GPU architectures

Here, we investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses SIMD/GPU accelerators locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations. We demonstrate our method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the reduction method is 50x faster than the sparse linear solver HSL MA57 running in serial on the CPU, and 6x faster than Pardiso running in parallel on CPU on the same number of processes.

97 MATHEMATICS AND COMPUTING↗

Multilane Automated Driving With Optimal Control and Mixed-Integer Programming

Road vehicle lane changes often initiate traffic disturbances and can therefore impact road networks’ energy and time efficiency. Furthermore, unexpected changes in traffic conditions may also render lane changes counterproductive for the lane-changing vehicle. Vehicle-to-vehicle connectivity combined with anticipative control could address these challenges via improved lane change decisions by automated vehicles. In a move toward this objective, receding horizon control cast as a mixed-integer quadratic program is used to plan lane changing and acceleration in a coupled optimization. A long-term pacing module, based on Pontryagin’s minimum principle from optimal control theory, sets terminal and input references for receding horizon control to target a user’s expected travel time. To remove nonlinear vehicle dynamics from the receding horizon controller, lane change commands are passed to a pure pursuit steering module whose response is approximated by a second-order linear model. Here, comparison against a rule-based reactive algorithm in arterial and highway scenarios shows an 8.9%–13.7% reduction in energy consumption and a 5.2%–10.3% reduction in the travel time, along with navigational improvements.

33 ADVANCED PROPULSION SYSTEMS↗

The effect of repeated hurricanes on the age of organic carbon in humid tropical forest soil

Abstract Increasing hurricane frequency and intensity with climate change is likely to affect soil organic carbon (C) stocks in tropical forests. We examined the cycling of C between soil pools and with depth at the Luquillo Experimental Forest in Puerto Rico in soils over a 30‐year period that spanned repeated hurricanes. We used a nonlinear matrix model of soil C pools and fluxes (“soilR”) and constrained the parameters with soil and litter survey data. Soil chemistry and stable and radiocarbon isotopes were measured from three soil depths across a topographic gradient in 1988 and 2018. Our results suggest that pulses and subsequent reduction of inputs caused by severe hurricanes in 1989, 1998, and two in 2017 led to faster mean transit times of soil C in 0–10 cm and 35–60 cm depths relative to a modeled control soil with constant inputs over the 30‐year period. Between 1988 and 2018, the occluded C stock increased and δ 13 C in all pools decreased, while changes in particulate and mineral‐associated C were undetectable. The differences between 1988 and 2018 suggest that hurricane disturbance results in a dilution of the occluded light C pool with an influx of young, debris‐deposited C, and possible microbial scavenging of old and young C in the particulate and mineral‐associated pools. These effects led to a younger total soil C pool with faster mean transit times. Our results suggest that the increasing frequency of intense hurricanes will speed up rates of C cycling in tropical forests, making soil C more sensitive to future tropical forest stressors.

14C↗

Data-Driven Modeling and Correction of Vehicle Dynamics

We develop a data-driven framework for learning and correcting nonautonomous vehicle dynamics. Physics-based vehicle models are often simplified for tractability and therefore exhibit inherent model-form uncertainty, motivating the need for data-driven correction. Moreover, nonautonomous dynamics are governed by time-dependent control inputs, which pose challenges in learning predictive models directly from temporal snapshot data. To address these, we reformulate the vehicle dynamics via a local parameterization of the time-dependent inputs, yielding a modified system composed ofa sequence of local parametric dynamical systems. Here, we approximate these parametric systems using two complementary approaches. First, we employ the dimension reduction and interpolation in parameter space (DRIPS) methodology to construct efficient linear surrogate models, equipped with lifted observable spaces and manifold-based operator interpolation. This enables data-efficient learning of vehicle models whose dynamics admit accurate linear representations in the lifted spaces. Second, for more strongly nonlinear systems, we employ flow map learning (FML), a deep neural network (DNN) approach that approximates the parametric evolution map without requiring special treatment of nonlinearities. We further extend FML with a transfer-learning-based model correction procedure, enabling the correction of misspecified prior models using only a sparse set of high-fidelity or experimental measurements, without assuming a prescribed form for the correction term. Through a suite of numerical experiments on unicycle, simplified bicycle, and slip-based bicycle models, we demonstrate that DRIPS offers robust and highly data-efficient learning of nonautonomous vehicle dynamics, while FML provides expressive nonlinear modeling and effective correction of model-form errors under severe data scarcity.

data-driven modeling↗

A NASTRAN/TREETOPS solution to a flexible, multi-body dynamics and controls problem on a UNIX workstation

Demands for nonlinear time history simulations of large, flexible multibody dynamic systems has created a need for efficient interfaces between finite-element modeling programs and time-history simulations. One such interface, TREEFLX, an interface between NASTRAN and TREETOPS, a nonlinear dynamics and controls time history simulation for multibody structures, is presented and demonstrated via example using the proposed Space Station Mobile Remote Manipulator System (MRMS). The ability to run all three programs (NASTRAN, TREEFLX and TREETOPS), in addition to other programs used for controller design and model reduction (such as DMATLAB and TREESEL, both described), under a UNIX Workstation environment demonstrates the flexibility engineers now have in designing, developing and testing control systems for dynamically complex systems.

Benavente, Javier E.↗

Processing and Testing of Thermoplastic Composite Cylindrical Shells Fabricated by Automated Fiber Placement

Two 61-cm-diameter eight-ply quasi-isotropic IM7/PEEK cylindrical shells were fabricated by automated fiber placement the NASA Langley Research Center using only infrared radiant heat to preheat the substrate and incoming composite uni-tape. The shells were characterized by ultrasonic c-scans for overall consolidation quality, and by optical microscopy and acid digestion for void content. Compression tests were also performed. Although the material used in the study was of generally poor quality due to numerous splits and dry fiber regions, the process was able to achieve a net reduction in void content in the as-placed component. Microscopy of the composite shells revealed well-consolidated, void-free interfaces. The two cylinders were then tested in uni-axial compression in a 1334 kN-capacity hydraulic test machine until buckling occurred. A geometrically nonlinear finite element analysis was conducted, and the differences between the predicted and measured values were 18.0 and 25.8%, respectively. Inclusion of measured imperfections of the cylinder into the analysis is expected to reduce these differences.

Hulcher, Anthony Bruce↗

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Architecture and Performance Predictions

A peak-seeking control approach for real-time trim configuration optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control approach is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an F/A-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are controlled for optimization of fuel flow. This paper presents the design and integration of this peak-seeking controller on a modified NASA F/A-18 airplane with research flight control computers. A research flight was performed to collect data to build a realistic model of the performance function and characterize measurement noise. This model was then implemented into a nonlinear six-degree-of-freedom F/A-18 simulation along with the peak-seeking control algorithm. With the goal of eventual flight tests, the algorithm was first evaluated in the improved simulation environment. Results from the simulation predict good convergence on minimum fuel flow with a 2.5-percent reduction in fuel flow relative to the baseline trim of the aircraft.

performance optimization↗

Vehicle-to-manufacturing (V2M) system: A novel approach to improve energy demand flexibility for demand response towards sustainable manufacturing

The U.S. manufacturing sector accounts for 77% of the industrial energy consumption, but its participation in demand response (DR) programs is largely lagged behind. The limited flexibility in production scheduling under high capacity utilization and the lack of DR schemes that incorporate energy demand flexibility measures are deemed as major barriers. In this study, a framework of the interactive vehicles-to-manufacturing (V2M) energy sharing system is proposed to improve the energy demand flexibility of the aggregated system and enhance the DR effectiveness for manufacturers. The V2M-based DR scheme aims to reduce the energy cost by load shifting through joint production and energy sharing control, which can eventually promote manufacturing DR implementation even under high capacity utilization requirements and dynamic real-time electricity prices. The V2M system is modeled based on a discrete-Markov chain considering the complex interconnections among various manufacturing resources and multi-directional energy flows among manufacturing facilities, electric vehicles, and the power grid. Based on the system model, a mixed-integer nonlinear programming (MINLP) problem is formulated to identify the optimal DR scheme. The effectiveness of the proposed approach is validated through comparisons with traditional manufacturing DR schemes. Finally, the results show that a 2.1 to 6.5 times energy demand flexibility improvement and an additional 4.7% to 6.9% energy cost reduction can be achieved by the proposed approach.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simulated impact of fill tube geometry on recent high-yield implosions at the National Ignition Facility

Inertial confinement fusion capsules fielded at the National Ignition Facility are filled with deuterium and tritium fuel by means of a fill tube. The fill tube introduces a low-density pathway into the fuel region of the capsule that allows high Z contaminant to invade the hot spot during the course of the implosion. A recent series of nominally identical high-yield implosions on the NIF has exhibited significant variability in performance. We evaluate the impact of the fill tube in these implosions computationally to determine whether variations in fill tube geometry could have contributed to this variability. The main contrast between the fill tube geometry in the six shots was the outer diameter of the capsule bore hole, a conical hole into which the fill tube is inserted. In our simulations, the geometry of the bore hole can play a significant role in the development of nonlinear flows seeded by the fill tube. Additionally, we find that the amount of space between the bore hole and the fill tube is the primary factor that determines the amount of contaminant jetted into the hot spot by the fill tube and, in turn, the level of yield reduction due to the fill tube in our simulations. As a consequence, some capsules with 5 μm fill tubes are predicted to outperform capsules with 2 μm fill tubes. We also find that micrometer-scale changes to bore hole size can impact fusion yields by up to four times near the ignition threshold. Nevertheless, simulation trends do not reproduce experimental yield trends, suggesting that the fill tube geometry was not the primary factor contributing to the observed variability in performance and that the fill tube could be masking sensitivity to other asymmetries such as other micrometer-scale capsule defects like voids that were not included in our simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Rapid data comparison technique development for N-ray/X-ray analysis

An operational N-ray/X-ray computer analysis system for nondestructive testing problems was realized by developing: (1) more efficient data transfer routines for use in existing image processing systems; (2) the capability to store and analyze multiple N-ray/X-ray images; (3) an efficient 2-D Fast Fourier Transform; and (4) transfer functions compensation capability to correct for distortions and nonlinearities of N-rays/X-rays. Applications of the technique are summarized for the following: thermal protection system tile with hydraulic fluid; thermal protection system tile with corrosion; thermal protection system tile debonding; stainless steel tubing corrosion; defective electrical devices; fast film restoration and deblur; corrosion tracking over time; reduction of source vignetting, and identification of materials (0 rings, Teflon parts, soft goods) via table look up.

Butterfield, R. L.↗

Peak-Seeking Optimization of Trim for Reduced Fuel Consumption: Architecture and Performance Predictions

A peak-seeking control approach for real-time trim configuration optimization for reduced fuel consumption has been developed by researchers at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center to address the goals of the NASA Environmentally Responsible Aviation project to reduce fuel burn and emissions. The peak-seeking control approach is based on a steepest-descent algorithm using a time-varying Kalman filter to estimate the gradient of a performance function of fuel flow versus control surface positions. In real-time operation, deflections of symmetric ailerons, trailing-edge flaps, and leading-edge flaps of an FA-18 airplane (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) are controlled for optimization of fuel flow. This presentation presents the design and integration of this peak-seeking controller on a modified NASA FA-18 airplane with research flight control computers. A research flight was performed to collect data to build a realistic model of the performance function and characterize measurement noise. This model was then implemented into a nonlinear six-degree-of-freedom FA-18 simulation along with the peak-seeking control algorithm. With the goal of eventual flight tests, the algorithm was first evaluated in the improved simulation environment. Results from the simulation predict good convergence on minimum fuel flow with a 2.5-percent reduction in fuel flow relative to the baseline trim of the aircraft.

peak-seeking↗

Impact of dynamic Jahn-Teller effect on magnetic excitations, lattice vibration, and thermal conductivity in U 𝑥 ⁢T⁢h 1−𝑥 ⁡O 2 system

Vibrational and magnetic properties of single-crystal uranium-thorium dioxide (U 𝑥 ⁢T⁢h 1−𝑥⁡ O 2 ) with a full range of 0 < 𝑥 < 1 are investigated. Thorium dioxide is a diamagnet whose thermal properties are governed by lattice vibration. The addition of paramagnetic uranium ions leads to the emergence of magnetic effects that alter the thermophysical properties noticeably even at room temperature. The interaction of phonons with magnetic moments of uranium 5⁢𝑓 electrons mediated by magnetoelastic coupling results in an anomalous low-temperature thermal conductivity profile. Analysis of the magnetic susceptibility measurements indicates a departure from the Curie-Weiss relationship characteristic of noninteracting paramagnetic ions, previously associated with the dynamic Jahn-Teller (DJT) effect characterized by coupling between spin and the oxygen sublattice. The T 2⁢g Raman peak position follows a nonlinear trend as a function of uranium concentration and hints that these Raman active optical modes play a role in either DJT or in mediating quadrupole-quadrupole interactions. A first-principle-based thermal transport model is implemented to explain the low-temperature transport measurements, where the anomalous reduction is attributed to phonon-spin resonant scattering. The interplay between spins and phonons is also captured using high-resolution inelastic x-ray scattering (IXS) measurements of phonon linewidths. Our results provide insights into the phonon interactions with the magnetic excitations governing DJT effect and impacting the low-temperature thermal transport processes in this material system. Furthermore, these findings have implications for understanding low-temperature thermal transport and magnetic properties in advanced materials for information processing and energy applications.

36 - MATERIALS SCIENCE↗