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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation

In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2008

Improving the Optical Quality Factor of the WGM Resonator

Resonators usually are characterized with two partially dependent values: finesse (F) and quality factor (Q). The finesse of an empty Fabry-Perot (FP) resonator is defined solely by the quality of its mirrors and is calculated as F=piR(exp 1/2)/(1-R). The maximum up-to-date value of reflectivity R approximately equal to 1 - 1.6 x 10(exp -6) is achieved with dielectric mirrors. An FP resonator made with the mirrors has finesse F=1.9 x 10(exp 6). Further practical increase of the finesse of FP resonators is problematic because of the absorption and the scattering of light in the mirror material through fundamental limit on the reflection losses given by the internal material losses and by thermodynamic density fluctuations on the order of parts in 109. The quality factor of a resonator depends on both its finesse and its geometrical size. A one-dimensional FP resonator has Q=2 F L/lambda, where L is the distance between the mirrors and lambda is the wavelength. It is easy to see that the quality factor of the resonator is unlimited because L is unlimited. F and Q are equally important. In some cases, finesse is technically more valuable than the quality factor. For instance, buildup of the optical power inside the resonator, as well as the Purcell factor, is proportional to finesse. Sometimes, however, the quality factor is more valuable. For example, inverse threshold power of intracavity hyperparametric oscillation is proportional to Q(exp 2) and efficiency of parametric frequency mixing is proportional to Q(exp 3). Therefore, it is important to know both the maximally achievable finesse and quality factor values of a resonator. Whispering gallery mode (WGM) resonators are capable of achieving larger finesse compared to FP resonators. For instance, fused silica resonators with finesse 2.3 x 10(exp 6) and 2.8 x 10(exp 6) have been demonstrated. Crystalline WGM resonators reveal even larger finesse values, F=6.3 x 10(exp 6), because of low attenuation of light in the transparent optical crystals. The larger values of F and Q result in the enhancement of various nonlinear processes. Low-threshold Raman lasing, optomechanical oscillations, frequency doubling, and hyperparametric oscillations based on these resonators have been recently demonstrated. Theory predicts a possibility of nearly 10(exp 14) room-temperature optical Q-factors of optical crystalline WGM resonators, which correspond to finesse levels higher than 10(exp 9). Experiments have shown numbers a thousand times lower than that. The difference occurs due to media imperfections. To substantially reduce the optical losses caused by the imperfections, a specific, multi-step, asymptotic processing of the resonator is implemented. The technique has been initially developed to reduce microwave absorption in dielectric resonators. One step of the process consists of mechanical polishing performed after high temperature annealing. Several steps repeat one after another to lead to significant reduction in optical attenuation and, as a result, to the increase of Q-factor as well as finesse of the resonator which demonstrates a CaF2 WGM resonator with F greater than 10(exp 7) and Q greater than 10(exp 11).

Savchenkov, Anatoliy↗
DOE OSTI2025

Multi-Agent Hierarchical Deep Reinforcement Learning for HVAC Control With Flexible DERs

As electricity consumption in commercial and residential buildings continues to rise, reducing energy costs presents an increasing challenge. Heating, ventilating, and air-conditioning (HVAC) systems, which typically account for 40%-50% of a building's energy use, are prime targets for energy savings. Intelligent control of HVAC temperature through the exploitation of HVAC load flexibility brings significant potential to reduce energy consumption and electricity expenses. The nonlinear models of HVAC systems challenge traditional control methods, while the uncertainty introduced by HVAC load flexibility complicates distributed energy resource (DER) management using conventional optimal dispatch techniques. In response to these challenges, we propose a hierarchical multi-agent deep reinforcement learning (DRL) approach. The lower-level agents focus on balancing comfort and energy conservation, while the upper-level DRL agents optimize the use of DERs to reduce peak demand based on the control outcomes of the HVAC by the lower-level agents. Here, in the upper-level agents, we incorporate a multi-agent structure based on ensemble learning, which acts based on historical and current data without relying on precise load forecasting to address the delayed rewarding issue in DRL. This allows for the effective reduction of energy costs. The proposed method is tested using a real-world microgrid comprising 413 buildings in Southern California, and the results demonstrate that our approach can significantly reduce overall electricity bills while ensuring the comfort of consumers and residents.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
NASA NTRS2014

Elastomeric Structural Attachment Concepts for Aircraft Flap Noise Reduction - Challenges and Approaches to Hyperelastic Structural Modeling and Analysis

Airframe noise is a significant part of the overall noise of transport aircraft during the approach and landing phases of flight. Airframe noise reduction is currently emphasized under the Environmentally Responsible Aviation (ERA) and Fixed Wing (FW) Project goals of NASA. A promising concept for trailing-edge-flap noise reduction is a flexible structural element or link that connects the side edges of the deployable flap to the adjacent main-wing structure. The proposed solution is distinguished by minimization of the span-wise extent of the structural link, thereby minimizing the aerodynamic load on the link structure at the expense of increased deformation requirement. Development of such a flexible structural link necessitated application of hyperelastic materials, atypical structural configurations and novel interface hardware. The resulting highly-deformable structural concept was termed the FLEXible Side Edge Link (FLEXSEL) concept. Prediction of atypical elastomeric deformation responses from detailed structural analysis was essential for evaluating feasible concepts that met the design constraints. The focus of this paper is to describe the many challenges encountered with hyperelastic finite element modeling and the nonlinear structural analysis of evolving FLEXSEL concepts. Detailed herein is the nonlinear analysis of FLEXSEL concepts that emerged during the project which include solid-section, foamcore, hollow, extended-span and pre-stressed concepts. Coupon-level analysis performed on elastomeric interface joints, which form a part of the FLEXSEL topology development, are also presented.

Sreekantamurthy, Thammaiah↗
NASA NTRS1991

Nonlinear dynamics of global atmospheric and Earth-system processes

General Circulation Model (GCM) studies of the atmospheric response to change boundary conditions are discussed. Results are reported on an extensive series of numerical studies based on the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM) general circulation model. In these studies the authors determined the response to systematic changes in atmospheric CO2 ranging from 100 to 1000 ppm; to changes in the prescribed sea surface temperature (SST) in the Gulf of Mexico, such as occurred during the deglaciation phase of the last ice age; to changes in soil moisture over North America; and to changes in sea ice extent in the Southern Hemisphere. Study results show that the response of surface temperature and other variables is nearly logarithmic, with lower levels of CO2 implying greater sensitivity of the atmospheric state to changes in CO2. It was found that the surface temperature of the Gulf of Mexico exerts considerable control over the storm track and behavior of storm systems over the North Atlantic through its influence on evaporation and the source of latent heat. It was found that reductions in soil moisture can play a significant role in amplifying and maintaining North American drought, particularly when a negative soil moisture anomaly prevails late in the spring.

Saltzman, Barry↗
DOE OSTI2022

Experimental realization of one dimensional helium

As the spatial dimension is lowered, locally stabilizing interactions are reduced, leading to the emergence of strongly fluctuating phases of matter without classical analogues. Here we report on the experimental observation of a one dimensional quantum liquid of 4 He using nanoengineering by confining it within a porous material preplated with a noble gas to enhance dimensional reduction. The resulting excitations of the confined 4 He are qualitatively different than bulk superfluid helium, and can be analyzed in terms of a mobile impurity allowing for the characterization of the emergent quantum liquid beyond the Luttinger liquid paradigm. The low dimensional helium system offers the possibility of tuning via pressure—from weakly interacting, all the way to the super Tonks-Girardeau gas of strongly interacting hard-core particles.

36 MATERIALS SCIENCE↗
NASA NTRS2008

An Integrated Approach to Damage Accommodation in Flight Control

In this paper we present an integrated approach to in-flight damage accommodation in flight control. The approach is based on Multiple Models, Switching and Tuning (MMST), and consists of three steps: In the first step the main objective is to acquire a realistic aircraft damage model. Modeling of in-flight damage is a highly complex problem since there is a large number of issues that need to be addressed. One of the most important one is that there is strong coupling between structural dynamics, aerodynamics, and flight control. These effects cannot be studied separately due to this coupling. Once a realistic damage model is available, in the second step a large number of models corresponding to different damage cases are generated. One possibility is to generate many linear models and interpolate between them to cover a large portion of the flight envelope. Once these models have been generated, we will implement a recently developed-Model Set Reduction (MSR) technique. The technique is based on parameterizing damage in terms of uncertain parameters, and uses concepts from robust control theory to arrive at a small number of "centered" models such that the controllers corresponding to these models assure desired stability and robustness properties over a subset in the parametric space. By devising a suitable model placement strategy, the entire parametric set is covered with a relatively small number of models and controllers. The third step consists of designing a Multiple Models, Switching and Tuning (MMST) strategy for estimating the current operating regime (damage case) of the aircraft, and switching to the corresponding controller to achieve effective damage accommodation and the desired performance. In the paper present a comprehensive approach to damage accommodation using Model Set Design,MMST, and Variable Structure compensation for coupling nonlinearities. The approach was evaluated on a model of F/A-18 aircraft dynamics under control effector damage, augmented by nonlinear cross-coupling terms and a structural dynamics model. The proposed approach achieved excellent performance under severe damage effects.

Boskovic, Jovan D.↗
NASA NTRS1974

Low thrust space vehicle trajectory optimization using regularized variables

Optimizing the trajectory of a low thrust space vehicle usually means solving a nonlinear two point boundary value problem. In general, accuracy requirements necessitate extensive computation times. In celestial mechanics, regularizing transformations of the equations of motion are used to eliminate computational and analytical problems that occur during close approaches to gravitational force centers. It was shown in previous investigations that regularization in the formulation of the trajectory optimization problem may reduce the computation time. In this study, a set of regularized equations describing the optimal trajectory of a continuously thrusting space vehicle is derived. The computational characteristics of the set are investigated and compared to the classical Newtonian unregularized set of equations. The comparison is made for low thrust, minimum time, escape trajectories and numerical calculations of Keplerian orbits. The comparison indicates that in the cases investigated for bad initial guesses of the known boundary values a remarkable reduction in the computation time was achieved. Furthermore, the investigated set of regularized equations shows high numerical stability even for long duration flights and is less sensitive to errors in the guesses of the unknown boundary values.

Schwenzfeger, K. J.↗
NASA NTRS1991

Combined magnetic vector-scalar potential finite element computation of 3D magnetic field and performance of modified Lundell alternators in Space Station applications

A method of combined use of magnetic vector potential (MVP) based finite element (FE) formulations and magnetic scalar potential (MSP) based FE formulations for computation of three-dimensional (3D) magnetostatic fields is developed. This combined MVP-MSP 3D-FE method leads to considerable reduction by nearly a factor of 3 in the number of unknowns in comparison to the number of unknowns which must be computed in global MVP based FE solutions. This method allows one to incorporate portions of iron cores sandwiched in between coils (conductors) in current-carrying regions. Thus, it greatly simplifies the geometries of current carrying regions (in comparison with the exclusive MSP based methods) in electric machinery applications. A unique feature of this approach is that the global MSP solution is single valued in nature, that is, no branch cut is needed. This is again a superiority over the exclusive MSP based methods. A Newton-Raphson procedure with a concept of an adaptive relaxation factor was developed and successfully used in solving the 3D-FE problem with magnetic material anisotropy and nonlinearity. Accordingly, this combined MVP-MSP 3D-FE method is most suited for solution of large scale global type magnetic field computations in rotating electric machinery with very complex magnetic circuit geometries, as well as nonlinear and anisotropic material properties.

Wang, Ren H.↗
NASA NTRSDate not supplied

Computational Assessment of the Opening Mode Facesheet-to-Core Fracture Toughness of a Thin-Gage Facesheet Sandwich Structure

The single-cantilevered beam test (SCB) is currently being considered for standardization for measuring the opening mode fracture toughness of sandwich structures subjected to facesheet-to-core separation. At present, the Area method is considered the benchmark method for determining energy required to propagate the crack. This method is problematic since it requires periodic unloading and cannot be evaluated instantaneously. Therefore, a method is needed to accurately and instantaneously calculate opening mode fracture toughness. In previous work, SCB tests with thin-gauge facesheets were performed to compare the J-integral approach to data reduction methods including the Modified Beam Theory (MBT) and the Area method. Large deviations (over 10%) in the calculated fracture toughness were observed between the J-integral and MBT for large displacements. In this study, a simulated SCB experiment using a virtual crackclosure technique was performed to assess the previously observed discrepancies. The predicted load displacement response, the predicted angle at the load application point, and the predicted crack growth results are in good agreement (within 5%) with the corresponding experimental data. Using the predicted load and angles, the predicted J-integral is in good agreement (within 10%) with the experimentally obtained solution. Using the MBT approach, the previously observed deviations in the fracture toughness were determined to be associated with the empirical extrapolation procedure to estimate the crack-root rotation correction from nonlinear compliance curves. The nonlinear compliance is determined to be associated with large displacement effects observed within the predicted and experimental results.

Sandwich Composite↗
NASA NTRS1976

Optimization of multi-constrained structures based on optimality criteria

A weight-reduction algorithm is developed for the optimal design of structures subject to several multibehavioral inequality constraints. The structural weight is considered to depend linearly on the design variables. The algorithm incorporates a simple recursion formula derived from the Kuhn-Tucker necessary conditions for optimality, associated with a procedure to delete nonactive constraints based on the Gauss-Seidel iterative method for linear systems. A number of example problems is studied, including typical truss structures and simplified wings subject to static loads and with constraints imposed on stresses and displacements. For one of the latter structures, constraints on the fundamental natural frequency and flutter speed are also imposed. The results obtained show that the method is fast, efficient, and general when compared to other competing techniques. Extensions to the generality of the method to include equality constraints and nonlinear merit functions is discussed.

Rizzi, P.↗
NASA NTRS1981

A method for obtaining reduced-order control laws for high-order systems using optimization techniques

A method of synthesizing reduced-order optimal feedback control laws for a high-order system is developed. A nonlinear programming algorithm is employed to search for the control law design variables that minimize a performance index defined by a weighted sum of mean-square steady-state responses and control inputs. An analogy with the linear quadractic Gaussian solution is utilized to select a set of design variables and their initial values. To improve the stability margins of the system, an input-noise adjustment procedure is used in the design algorithm. The method is applied to the synthesis of an active flutter-suppression control law for a wind tunnel model of an aeroelastic wing. The reduced-order controller is compared with the corresponding full-order controller and found to provide nearly optimal performance. The performance of the present method appeared to be superior to that of two other control law order-reduction methods. It is concluded that by using the present algorithm, nearly optimal low-order control laws with good stability margins can be synthesized.

Mukhopadhyay, V.↗
DOE OSTI2021

Threshold Heat-Flux Reduction by Near-Resonant Energy Transfer

Near-resonant energy transfer to large-scale stable modes is shown to reduce transport above the linear critical gradient, contributing to the onset of transport at higher gradients. This is demonstrated for a threshold fluid theory of ion temperature gradient turbulence based on zonal-flow-catalyzed transfer. Furthermore, the heat flux is suppressed above the critical gradient by resonance in the triplet correlation time, a condition enforced by the wave numbers of the interaction of the unstable mode, zonal flow, and stable mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Global burned area increasingly explained by climate change

Fire behaviour is changing in many regions worldwide. However, nonlinear interactions between fire weather, fuel, land use, management and ignitions have impeded formal attribution of global burned area changes. Here, in this work, we demonstrate that climate change increasingly explains regional burned area patterns, using an ensemble of global fire models. The simulations show that climate change increased global burned area by 15.8% (95% confidence interval (CI) [13.1–18.7]) for 2003–2019 and increased the probability of experiencing months with above-average global burned area by 22% (95% CI [18–26]). In contrast, other human forcings contributed to lowering burned area by 19.1% (95% CI [21.9–15.8]) over the same period. Moreover, the contribution of climate change to burned area increased by 0.22% (95% CI [0.22–0.24]) per year globally, with the largest increase in central Australia. Our results highlight the importance of immediate, drastic and sustained GHG emission reductions along with landscape and fire management strategies to stabilize fire impacts on lives, livelihoods and ecosystems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1984

New Uses for Sensitivity Analysis: How Different Movement Tasks Effect Limb Model Parameter Sensitivity

Original results for a newly developed eight-order nonlinear limb antagonistic muscle model of elbow flexion and extension are presented. A wider variety of sensitivity analysis techniques are used and a systematic protocol is established that shows how the different methods can be used efficiently to complement one another for maximum insight into model sensitivity. It is explicitly shown how the sensitivity of output behaviors to model parameters is a function of the controller input sequence, i.e., of the movement task. When the task is changed (for instance, from an input sequence that results in the usual fast movement task to a slower movement that may also involve external loading, etc.) the set of parameters with high sensitivity will in general also change. Such task-specific use of sensitivity analysis techniques identifies the set of parameters most important for a given task, and even suggests task-specific model reduction possibilities.

Winters, J. M.↗
DOE OSTI2020

Cross-code comparison of the impact of the fill tube on high yield implosions on the National Ignition Facility

Fill tubes are used to inject deuterium and tritium fuel into inertial confinement fusion capsules fielded on the National Ignition Facility. These fill tubes have been shown to have a detrimental effect on capsule performance, primarily by introducing a low-density pathway into the central fuel region that enables the jetting of ablation material into the hot spot. Due to the complexity of the highly nonlinear flow associated with the fill tube and the challenge of diagnosing the evolution of the fill tube jet late in the implosion experiments, the uncertainty in how this perturbation source evolves is great. In this work, we report on the results of a detailed code comparison performed to understand uncertainties in computational modeling of the impact of fill tubes on implosion performance. The study employed two radiation-hydrodynamics codes, HYDRA and xRAGE, which employ very different meshing strategies and hydrodynamics solvers, as well as two radiation transport methodologies, discrete ordinates and multi-group diffusion. Our results demonstrate generally good agreement between codes through most of the implosion although they indicate sensitivity to opacity averaging methods. Late in the implosion, differences arise in the distribution and amount of contaminant although these differences have a remarkably small impact on the amount of yield reduction due to the fill tube. While these results demonstrate sensitivity in fill tube modeling to algorithmic choices, the observed differences between codes are small relative to known sensitivities due to expected variations in the fill tube geometry. Finally, we have developed a methodology for performing multi-group diffusion simulations that show good agreement with the more accurate discrete ordinates method

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Online Optimization of NSLS-II Dynamic Aperture and Injection Transient (Facility Improvement Project)

The Facility Improvement Project (FIP) “Methods of online optimization of NSLS-II storage ring concurrent with user operations” proposed and approved in 2018 is now complete. The first goal is the online optimization of nonlinear beam dynamics to increase the beam lifetime and injection efficiency. We have developed a model-independent optimization technique using advanced algorithms. Using this technique, we increased the NSLS-II dynamic aperture by 20% and reduced the amplitude-dependent tune shift by a factor of two. We applied sextupole optimization to the new high-chromaticity lattice, which has been developed to improve the beam stability and to increase the single-bunch beam intensity. We were able to double the injection efficiency and exceed 90%. The second goal of this project is to provide the top-off injection with minimized perturbations of the beamline user operations. To achieve the transparent injection, we optimize the matching of four storage ring injection kickers reducing the perturbation of the stored beam orbit. The first application of the online optimization resulted in a reduction of the amplitude of residual beam oscillation by a factor of six, from 1300 μm down to 200 μm. Further improvement was limited by the timing jitter of the trigger boards. Then, we replaced all trigger boards with a new design and reduced the timing jitter by a factor of five, from 5 ns down to 1 ns. This upgrade resulted in the injection transient reduction from 200 μm to 120 μm. To optimize the full set of kicker parameters, including the trigger timing, amplitude, and pulse width, we upgraded all kicker power supplies with the capability of tunable waveform width. As a result, we have reduced the injection transient to the limit of 60 μm.

43 PARTICLE ACCELERATORS↗