Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “successive linear approximation”

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 163 records · Page 9

Advanced data assimilation in strongly nonlinear dynamical systems

Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Miller, Robert N.↗

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage (CRADA Report)

NREL and Element 16 collaborated on sulfur thermal energy storage modeling using NREL’s high performance computing (HPC) resources to assist its application in industrial processes. Industrial process heat (IPH) accounts for ~70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately, 1500 TWht (~60% Terawatt hour thermal) of IPH demand is in the temperature range of 100-300°C. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. This HPC4EI project facilitated Element 16’s development of low-cost and high-impact molten sulfur TES for dispatchable IPH. The development of a high-fidelity model validated by experimental data and HPC simulations enabled the successful resolution of the complex interplay between fluid dynamics and heat transfer processes during transient operation of sulfur TES, overcoming the numerical challenges posed by the non-linear temperature-dependent physical properties of sulfur. The project helped accelerate Element 16’s molten sulfur TES product design and support its broad applications.

25 ENERGY STORAGE↗

A reliable algorithm for optimal control synthesis

In recent years, powerful design tools for linear time-invariant multivariable control systems have been developed based on direct parameter optimization. In this report, an algorithm for reliable optimal control synthesis using parameter optimization is presented. Specifically, a robust numerical algorithm is developed for the evaluation of the H(sup 2)-like cost functional and its gradients with respect to the controller design parameters. The method is specifically designed to handle defective degenerate systems and is based on the well-known Pade series approximation of the matrix exponential. Numerical test problems in control synthesis for simple mechanical systems and for a flexible structure with densely packed modes illustrate positively the reliability of this method when compared to a method based on diagonalization. Several types of cost functions have been considered: a cost function for robust control consisting of a linear combination of quadratic objectives for deterministic and random disturbances, and one representing an upper bound on the quadratic objective for worst case initial conditions. Finally, a framework for multivariable control synthesis has been developed combining the concept of closed-loop transfer recovery with numerical parameter optimization. The procedure enables designers to synthesize not only observer-based controllers but also controllers of arbitrary order and structure. Numerical design solutions rely heavily on the robust algorithm due to the high order of the synthesis model and the presence of near-overlapping modes. The design approach is successfully applied to the design of a high-bandwidth control system for a rotorcraft.

Vansteenwyk, Brett↗

Incorporation of star measurements for the determination of orbit and attitude parameters of a geosynchronous satellite: An iterative application of linear regression

Currently on NOAA/NESS's VIRGS system at the World Weather Building star images are being ingested on a daily basis. The image coordinates of the star locations are measured and stored. Subsequently, the information is used to determine the attitude, the misalignment angles between the spin axis and the principal axis of the satellite, and the precession rate and direction. This is done for both the 'East' and 'West' operational geosynchronous satellites. This orientation information is then combined with image measurements of earth based landmarks to determine the orbit of each satellite. The method for determining the orbit is simple. For each landmark measurement one determines a nominal position vector for the satellite by extending a ray from the landmark's position towards the satellite and intersecting the ray with a sphere with center coinciding with the Earth's center and with radius equal to the nominal height for a geosynchronous satellite. The apparent motion of the satellite around the Earth's center is then approximated with a Keplerian model. In turn the variations of the satellite's height, as a function of time found by using this model, are used to redetermine the successive satellite positions by again using the Earth based landmark measurements and intersecting rays from these landmarks with the newly determined spheres. This process is performed iteratively until convergence is achieved. Only three iterations are required.

Phillips, D.↗

Accretion shock geometries in the magnetic variables

The first self consistent shock models for the AM Herculis-type systems successfully identified the dominant physical processes and their signatures. These homogenous shock models predict unpolarized, Rayleigh-Jeans optical spectra with sharp cutoffs and rising polarizations as the shocks become optically thin in the ultraviolet. However, the observed energy distributions are generally flat with intermediate polarizations over a broad optical band. These and other observational evidence support a non-homogenous accretion profile which may extend over a considerable fraction of the stellar surface. Both the fundamental assumptions underlying the canonical 1-D shock model and the extension of this model to inhomogenous accretion shocks were identified, for both radial and linear structures. The observational evidence was also examined for tall shocks and little evidence was found for relative shock heights in excess of h/R(1) greater than or equal to 0.1. For several systems, upper limits to the shock height can be obtained from either x ray or optical data. These lie in the region h/R(1) is approximately 0.01 and are in general agreement with the current physical picture for these systems. The quasi-periodic optical variations observed in several magnetic variables may eventually prove to be a major aid in further understanding their accretion shock geometries.

Stockman, H. S.↗

Randomized Algorithms for Symmetric Nonnegative Matrix Factorization

Symmetric Nonnegative Matrix Factorization (SymNMF) is a technique in data analysis and machine learning that approximates a matrix with a product of a nonnegative, low-rank matrix and it transpose. To design faster and more scalable algorithms for SymNMF we develop two randomized algorithms for its computation. The first method uses randomized matrix sketching to compute an initial low-rank approximation to the input matrix and proceeds to uses this as a low-rank input to rapidly compute a SymNMF. The second methods uses randomized leverage score sampling to approximately solve constrained least squares problems. Many successful methods for SymNMF rely on (approximately) solving sequences of constrained least squares problems. Here, we prove theoretically that leverage score sampling can approximately solve constrained least squares problems to e-accuracy. Finally we demonstrate both methods work in practice by applying them to graph clustering tasks on large real world data sets. These experiments show that our methods approximately maintain solution quality and achieve significant speed ups for both large dense and large sparse problems.

97 MATHEMATICS AND COMPUTING↗

Comparison of dynamical approximation schemes for non-linear gravitational clustering

We have recently conducted a controlled comparison of a number of approximations for gravitational clustering against the same n-body simulations. These include ordinary linear perturbation theory (Eulerian), the adhesion approximation, the frozen-flow approximation, the Zel'dovich approximation (describable as first-order Lagrangian perturbation theory), and its second-order generalization. In the last two cases we also created new versions of approximation by truncation, i.e., smoothing the initial conditions by various smoothing window shapes and varying their sizes. The primary tool for comparing simulations to approximation schemes was crosscorrelation of the evolved mass density fields, testing the extent to which mass was moved to the right place. The Zel'dovich approximation, with initial convolution with a Gaussian e(exp -k(exp 2)/k(exp 2, sub G)) where k(sub G) is adjusted to be just into the nonlinear regime of the evolved model (details in text) worked extremely well. Its second-order generalization worked slightly better. All other schemes, including those proposed as generalizations of the Zel'dovich approximation created by adding forces, were in fact generally worse by this measure. By explicitly checking, we verified that the success of our best-choice was a result of the best treatment of the phases of nonlinear Fourier components. Of all schemes tested, the adhesion approximation produced the most accurate nonlinear power spectrum and density distribution, but its phase errors suggest mass condensations were moved to slightly the wrong location. Due to its better reproduction of the mass density distribution function and power spectrum, it might be preferred for some uses. We recommend either n-body simulations or our modified versions of the Zel'dovich approximation, depending upon the purpose. The theoretical implication is that pancaking is implicit in all cosmological gravitational clustering, at least from Gaussian initial conditions, even when subcondensations are present.

Melott, Adrian L.↗

Multi-Parameter Optical Fiber Sensing of Humidity, CH4, CO2, and Corrosion

In this work, the previously demonstrated capability of the optical fiber sensor (OFS) for successful monitoring of humidity has been extended to monitor the humidity, CH4, and CO2 with different gas composition based on the strain produced along the single mode fiber (SMF) sensor. This is enabled by absorption of H2O/gases on to the commercially available polyacrylate coated jacketed portion of the fiber resulting a change in strain. Under equilibrium, a differential microstrain was observed along the jacketed portion of the SMF with N2, CH4, and CO2 at different humidity conditions, while the unjacketed portion of the fiber was used only for sensing pressure/temperature induced strain. In case of N2 at 800 psig pressure, the observed microstrain was approximately 80, 65, 50 and 35 µε at 100, 75.0, 46.8, and 23.4 RH% respectively. Comparatively, a microstrain of approximately 95 µε was observed with 100% RH CH4 which demonstrates that SMF produces a measurable CH4 response alongside water. Similarly, the observed microstrain with CO2 was approximately 105, 90, 85, 80, and 70 µε at 100, 75.0, 46.8, 23.4, and 0 RH% respectively. Linear regression and principal component analysis of these dataset provided deconvolution of the impact of strain from H2O, CH4, and CO2 enabling good cross sensitivity. Additionally, modified OFS comprising of Fe coated fiber section was employed to monitor corrosion, using Fe as corrosion proxy, under harsh corrosive environment.

Mainali, Badri↗

Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water

Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Particle Filtering for Model-Based Anomaly Detection in Sensor Networks

A novel technique has been developed for anomaly detection of rocket engine test stand (RETS) data. The objective was to develop a system that postprocesses a csv file containing the sensor readings and activities (time-series) from a rocket engine test, and detects any anomalies that might have occurred during the test. The output consists of the names of the sensors that show anomalous behavior, and the start and end time of each anomaly. In order to reduce the involvement of domain experts significantly, several data-driven approaches have been proposed where models are automatically acquired from the data, thus bypassing the cost and effort of building system models. Many supervised learning methods can efficiently learn operational and fault models, given large amounts of both nominal and fault data. However, for domains such as RETS data, the amount of anomalous data that is actually available is relatively small, making most supervised learning methods rather ineffective, and in general met with limited success in anomaly detection. The fundamental problem with existing approaches is that they assume that the data are iid, i.e., independent and identically distributed, which is violated in typical RETS data. None of these techniques naturally exploit the temporal information inherent in time series data from the sensor networks. There are correlations among the sensor readings, not only at the same time, but also across time. However, these approaches have not explicitly identified and exploited such correlations. Given these limitations of model-free methods, there has been renewed interest in model-based methods, specifically graphical methods that explicitly reason temporally. The Gaussian Mixture Model (GMM) in a Linear Dynamic System approach assumes that the multi-dimensional test data is a mixture of multi-variate Gaussians, and fits a given number of Gaussian clusters with the help of the wellknown Expectation Maximization (EM) algorithm. The parameters thus learned are used for calculating the joint distribution of the observations. However, this GMM assumption is essentially an approximation and signals the potential viability of non-parametric density estimators. This is the key idea underlying the new approach.

Solano, Wanda↗

Transportation analysis and related design optimization of the Fermilab high-beta 650 MHz cryomodule

The Proton Improvement Plan-II (PIP-II) at Fermi National Accelerator Laboratory (FNAL) will create a new and vastly improved accelerator, which will be the source of high-energy particles for the experiments taking place at FNAL. The new linear accelerator (LINAC) contains several types of cryomodules, which are individual particle accelerators. The last cryomodule in the LINAC will be the High-Beta 650 (HB650), which will operate at 650 MHz. Each module is approximately 15 meters in length and 1.5 meters wide, weighs 13 tonnes, and shares many design features with three of the other cryomodules. The HB650 consists of two primary sections, the outer vacuum vessel and the components that reside within it, which will be cooled to cryogenic temperatures – the cold mass. Once assembly is complete at FNAL, it will be transported to one of three places: another location at FNAL, another national laboratory within the continental United States, or to a scientific partner in Europe. Any excitation the module experiences, such as going over rough roads when on a semitrailer, can create high stresses in components and cause failure if severe enough. Additionally, any delicate components with low resonant frequencies that match the excitation spectrum could achieve resonance, potentially causing a fatigue failure by repeated flexing. To ensure the successful transportation of the HB650, a transport analysis utilizing ANSYS simulation software has been performed on major subassemblies, as well as analysis of the subassemblies as a combined system. The optimization of components, design of support structures, and overall increasing of the lowest resonant frequencies have resulted in a cryomodule that will be better equipped to handle transportation and any large impact loads that come with it.

43 PARTICLE ACCELERATORS↗

Fuzzy self-learning control for magnetic servo system

It is known that an effective control system is the key condition for successful implementation of high-performance magnetic servo systems. Major issues to design such control systems are nonlinearity; unmodeled dynamics, such as secondary effects for copper resistance, stray fields, and saturation; and that disturbance rejection for the load effect reacts directly on the servo system without transmission elements. One typical approach to design control systems under these conditions is a special type of nonlinear feedback called gain scheduling. It accommodates linear regulators whose parameters are changed as a function of operating conditions in a preprogrammed way. In this paper, an on-line learning fuzzy control strategy is proposed. To inherit the wealth of linear control design, the relations between linear feedback and fuzzy logic controllers have been established. The exercise of engineering axioms of linear control design is thus transformed into tuning of appropriate fuzzy parameters. Furthermore, fuzzy logic control brings the domain of candidate control laws from linear into nonlinear, and brings new prospects into design of the local controllers. On the other hand, a self-learning scheme is utilized to automatically tune the fuzzy rule base. It is based on network learning infrastructure; statistical approximation to assign credit; animal learning method to update the reinforcement map with a fast learning rate; and temporal difference predictive scheme to optimize the control laws. Different from supervised and statistical unsupervised learning schemes, the proposed method learns on-line from past experience and information from the process and forms a rule base of an FLC system from randomly assigned initial control rules.

Tarn, J. H.↗

A Data Processing Workflow for Fixed-wing Drone Based Radiation Mapping in the Chornobyl Exclusion Zone (CEZ) - 20119

April 2020 marks the 34. anniversary of the high-profile radiological release from the Chornobyl Nuclear Power Plant (ChNPP). The release of radioactive material from reactor number four began on the 26 April 1986 and continued over a period of about 10 days, releasing approximately 1700 PBq of radioactive material (including 85 PBq of 137-Cs) into the environment. To this day, the accident remains the most significant release of radioactive material since civil nuclear power generation began. In the years since the accident, automated and remote radiation monitoring technologies have advanced significantly in their capabilities. One such example of this is the use of unmanned aerial vehicles (UAVs) in radiation mapping investigations. In April 2019, a team of scientists from the University of Bristol showcased a novel radiation mapping system within the CEZ, specifically aiming to map radiation over a large portion of the area immediately surrounding the ChNPP. Over six days of data collection, the system flew a total distance of 583.8 km, covering an area of 14.6 sq.km with an exceptional spatial resolution (sub 20 m/pixel). The work presented herein outlines and explains the data processing procedure to convert the raw data into 137-Cs activity (kBq/sq.m) and cesium-equivalent dose-rate (CED) at 1 m above ground level (μSv/hr). A demonstration of the validity of the method is demonstrated through the successful reduction of the raw data into a single linear relationship between the measured {sup 137}Cs net peak intensity and the {sup 137}Cs activity. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reactive Modified Epoxy Resin and Its Miscible Blends Based on Recycled Oligomers from Solvolysis

Chemical depolymerization of fully cured epoxy resin with 20% reactive modifier was successfully performed via a solvent-assisted solvolysis process into low molecular weight recyclable oligomers (RO) at 240 °C in a pressure vessel at 650 psi for 4 h. The thermoset epoxy resin was depolymerized into transparent brown viscous fluid with a higher viscosity than the uncured epoxy resin with approximately 93% yield. Different concentrations of the RO were homogeneously mixed with the pure epoxy resin, and their curing kinetics, viscosity, FTIR, mechanical properties, DMA, and cross-link density were investigated. The curing kinetics of the pure reactive modified epoxy resin (baseline) and its mixtures with RO of different concentrations were investigated under both isothermal and nonisothermal conditions using small amplitude oscillatory shear flow. The elastic and viscous moduli (G′ and G″), complex viscosity (η*), and tan δ values were evaluated at different curing times and temperatures. The G′, G″, and η* increased dramatically, while tan δ decreased strongly by several orders of magnitude at the gel point. The zero-shear viscosity (η 0 ) was determined from the angular frequency dependent on η* based on the Cross model for different blend compositions in the liquid state before curing. The composition dependence of η 0 showed a positive deviation from the linear mixing rule and was well described by the Lecyar model. Here, the apparent activation energy of curing (E a ) was also evaluated according to the Arrhenius equation and was found to be 46 ± 2 kJ/mol regardless of the different contents of RO. For all blends up to 40 wt % RO, only one tan δ peak systematically shifting to lower temperatures with increasing content of RO was observed in the DMA measurements, indicating that the epoxy resin and the RO are miscible with up to 40 wt % RO.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Consistent Pl Analysis of Aqueous Uranium-235 Critical Assemblies

The lethargy-dependent equations of the consistent Pl approximation to the Boltzmann transport equation for slowing down neutrons have been used as the basis of an IBM 704 computer program. Some of the effects included are (1) linearly anisotropic center of mass elastic scattering, (2) heavy element inelastic scattering based on the evaporation model of the nucleus, and (3) optional variation of the buckling with lethargy. The microscopic cross-section data developed for this program covered 473 lethargy points from lethargy u = 0 (10 Mev) to u = 19.8 (0.025 ev). The value of the fission neutron age in water calculated here is 26.5 square centimeters; this value is to be compared with the recent experimental value given as 27.86 square centimeters. The Fourier transform of the slowing-down kernel for water to indium resonance energy calculated here compared well with the Fourier transform of the kernel for water as measured by Hill, Roberts, and Fitch. This method of calculation has been applied to uranyl fluoride - water solution critical assemblies. Theoretical results established for both unreflected and fully reflected critical assemblies have been compared with available experimental data. The theoretical buckling curve derived as a function of the hydrogen to uranium-235 atom concentration for an energy-independent extrapolation distance was successful in predicting the critical heights of various unreflected cylindrical assemblies. The critical dimensions of fully water-reflected cylindrical assemblies were reasonably well predicted using the theoretical buckling curve and reflector savings for equivalent spherical assemblies.

Fieno, Daniel↗

Planetary circulations in the presence of transient and self-induced heating

The linearized global primitive equation model of Salby and Garcia (1987) was modified to study tropospheric circulations and their interaction with tropical convection. The vertical resolution was increased to approximately 1.5 km, with the effective vertical domain extending the approximately 40 km (higher altitudes in the calculation are used to ensure the radiation upper boundary condition is satisfied). A formal wave conditional instability of the second kind (CISK) parameterization was introduced. That parameterization allows low-frequency disturbances to interact with convection by organizing surface moisture. Boundary layer friction has two important effects in these calculations: it exaggerates surface convergence near the equator by driving the flow out of geostrophic balance; and it introduces a phase shift between the convergence pattern and the circulation, which drives the heating out of quadrature with the temperature field and allows a positive feedback to take place. The observational complement of this theoretical work involves Global Cloud Imagery (GCI) composited from six satellites in the International Satellite Cloud Climatology Project (ISCCP) which were simultaneously observing the Earth's cloud field. Having global coverage, horizontal resolution of half a degree, and produced synoptically every 3 hours, the GCI affords an unprecendented view of global convection. A full year (1983 to 1984) of GCI has now been successfully created. Among other phenomena, that imagery reveals a regular progression of easterly waves in the convective pattern of the Western Pacific. Those disturbances have odd symmetry about the equator and organize convection into a pattern with similar symmetry. In time-mean maps of cloud cover, these propagating disturbances produce a split Inter-Tropical Convergence Zone (ITCZ) in the western Pacific.

Salby, Murry L.↗

Hybrid Thermally Efficient Core (HyTEC) HyTEC Phase 1 – Advanced Aerodynamics Final Report

The objective of the Hybrid Thermally Efficient Core (HyTEC) – Advanced High Pressure Turbine (HPT) Aerodynamics project is to develop technology for a compact core that contributes to significant fuel burn reductions of 5-10% over current generation technologies. To accomplish this, the HPT is incorporating a range of aerodynamic features and technologies to improve component efficiency and provide favorable systems level trades. In particular, this project explored low solidity airfoils, advanced tip treatments, platform contouring, and advanced ceramic matrix composite (CMC) Stage 2 Nozzle (S2N) airfoils to eliminate the need for post throat cooling in a compact core environment. The maturation of these technology areas is expected to provide a significant improvement in component efficiencies, and consequently reductions in fuel burn, over the current state of the art (SoA). To mature these technologies to Technology Readiness Level (TRL) 4, a test campaign was performed that consisted of four tests at three facilities. TRL 3 testing was performed in the CW22 linear cascade at NASA Glenn for blade and nozzle technologies, TRL 4 nozzle testing was performed at GE Aerospace (GEA) Test Cell A8, and TRL 4 blade testing was performed at the Notre Dame Turbomachinery Laboratory (NDTL) using the Transonic Research Turbine (TRT) rig. Low solidity was successfully demonstrated to a TRL 4 level. Low solidity nozzles showed benefits in line with pre-project expectations, while low solidity blades were shown to have an aerodynamic penalty. Crucially, this program only considered the aerodynamic losses, and systems trades such as reductions in cooling flows are expected to continue to make low solidity blades a net positive. By successfully quantifying the aerodynamic performance in this project, these trades can be conducted to determine where in the engine architecture low solidity blades will contribute positively to system operation. Platform contouring was demonstrated to a TRL 4 level, with performance in line with the lower end of the pre-project expected range. Advanced tip treatments performance levels were indeterminate, showing the expected improvements to flow physics but with a performance level confounded by several rig issues including whirl mode induced variation in tip clearance. The elimination of post-throat cooling on the S2N was successfully demonstrated to provide a performance benefit, however that benefit was approximately half the level that was expected in pre-project predictions. Overall, the technology maturation plan for HyTEC Phase 1 was successful, bringing the suite of technologies to TRL 4.

High pressure turbine↗