Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical analysis”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Numerical methods for nonlocal and fractional models

Partial differential equations (PDEs) are used with huge success to model phenomena across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDEs fail to adequately model observed phenomena, or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. Here, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis and of specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modelling and algorithmic extensions, which serve to show the wide applicability of nonlocal modelling.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

Characterization of the Heat Extraction Capability of a Compliant, Sliding, Thermal Interface for Use in a High Temperature, Vacuum Microgravity Furnace

Much of the material science gained in microgravity research requires processing a sample in a high performance furnace. One such furnace currently being designed is the Quench Module Insert (QMI). The Principle Investigators, for whom the furnace is proposed, require high temperature gradients in their cylindrical samples to achieve the science objectives. One of the components critical to achieving high sample axial temperature gradients in the Quench Module Insert is a high performance cold zone to extract the heat from the sample. This cold zone employs a compliant, sliding thermal interface based on a Vel-Therm felt. This felt provides a conductive path between the Sample Cartridge Assembly (SCA) exterior surface and the interior surface of the water cooled chill block while allowing movement of the sample relative to the chill block. The Vel-Therm felt is composed of long polymer-based fibers affixed to a thin flexible substrate layer. The fibers are oriented perpendicular to this substrate giving the felt the appearance of a velvet fabric. The Vel-Therm felt heat extraction capability was quantified in earlier tests performed in an inert gas environment. The current activity, described in this paper, is intended to characterize the extraction capability of Vel-Therm felt in a vacuum environment similar to the QMI environment. This testing is necessary to quantify the thermal performance of the Vel-Therm felt and the sensitivity of that performance to key variables. The data derived from these tests will be incorporated into the current thermal models to improve the quality of the models and reduce uncertainty of the analytical results. In addition, the data will be used to help select the appropriate Vel-Therm felt and set proper operating limits as well as assess the performance range of the furnace. The objective of this test is to measure the heat extraction rate of the Vel-Therm felt as specified by the effective heat transfer coefficient. Therefore, the test setup was designed to force the bulk of the heat transfer through the area where the Vel-Therm felt was applied. A heat source, consisting of a ceramic heating element encased in a Copper (Cu) housing is mounted on four isolated support rods. A 6-layer molybdenum radiation shield is used to insulate against heat loss from the heater and prevent heat exchange between the hot and cold sides of the test apparatus. The Vel-Therm felt is affixed to the surface of the cold sink, a water-cooled Cu chill block. An adjustable plate supports and isolates the cold sink from above and is used to control the amount of the deflection of the Vel- Therm when in contact with the Cu heating element housing. The primary means of establishing the power being conducted through the felt is to measure the energy being transferred to the water passing through the chill block. Analysis was performed to support the assumption that the source and sink surfaces were approximately isothermal under the specific test conditions. As a check on the amount of power passing through the felt, the power supplied to the heater was also measured. Thermocouples were strategically located throughout the test apparatus for measurement purposes. A bell jar was lowered over the assembly to impose vacuum conditions. Currently, variables tested have been fiber compression and fiber type and surface temperatures (both the hot and cold side temperatures are hypothesized to be important to the performance of the Vel-Therm.) Selected runs were repeated to ensure consistency and repeatability. Results obtained thus far reveal that Vel-Therm performance is significantly degraded by fibers being exposed to high compression. It also shows that performance is somewhat negatively impacted by previous compression, thereby, raising the question of repeatability. In addition, early results show a significant dependence on temperature. A computer aided mathematical analysis of the test setup is ongoing. The results will be correlated to actual results. The correlation will examine such details as parasitic loses, conduction down the power leads and many other concerns.

Bellomy-Ezell, Jenny↗
NASA NTRS1992

Digital tracking loops for a programmable digital modem

In this paper, an analysis and hardware emulation of the tracking loops for a very flexible programmable digital modem (PDM) will be presented. The modem is capable of being programmed for 2, 4, 8, 16-PSK, 16-QAM, MSK, and Offset-QPSK modulation schemes over a range of data rates from 2.34 to 300 Mbps with programmable spectral occupancy from 1.2 to 1.8 times the symbol rate; these operational parameters are executable in burst or continuous mode. All of the critical processing in both the modulator and demodulator is done at baseband with very high-speed digital hardware and memory. Quadrature analog front-ends are used for translation between baseband and the IF center frequency. The modulator is based on a table lookup approach, where precomputed samples are stored in memory and clocked out according to the incoming data pattern. The sample values are predistorted to counteract the effects of the other filtering functions in the link as well as any transmission impairments. The demodulator architecture was adapted from a joint estimator-detector (JED) mathematical analysis. Its structure is applicable to most signalling formats that can be represented in a two-dimensional space. The JED realization uses interdependent, mutually aiding tracking loops with post-detection data feedback. To expedite and provide for more reliable synchronization, initial estimates for these loops are computed in a parallel acquisition processor. The cornerstone of the demodulator realization is the pre-averager received data filter which allows operation over a broad range of data rates without any hardware changes and greatly simplifies the implementation complexity. The emulation results confirmed tracking loop operation over the entire range of operational parameters listed above, as well as the capability of achieving and maintaining synchronization at BER's in excess of 10(exp -1). The emulation results also showed very close agreement with the tracking loop analysis, and validated the resolution apportionment of the various hardware elements in the tracking loops.

Poklemba, John J.↗
DOE OSTI2021

Learning the solution operator of parametric partial differential equations with physics-informed DeepONets

Partial differential equations (PDEs) play a central role in the mathematical analysis and modeling of complex dynamic processes across all corners of science and engineering. Their solution often requires laborious analytical or computational tools, associated with a cost that is markedly amplified when different scenarios need to be investigated, for example, corresponding to different initial or boundary conditions, different inputs, etc. In this work, we introduce physics-informed DeepONets, a deep learning framework for learning the solution operator of arbitrary PDEs, even in the absence of any paired input-output training data. We illustrate the effectiveness of the proposed framework in rapidly predicting the solution of various types of parametric PDEs up to three orders of magnitude faster compared to conventional PDE solvers, setting a previously unexplored paradigm for modeling and simulation of nonlinear and nonequilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Defect and structural evolution under high-energy ion irradiation informs battery materials design for extreme environments

Understanding defect evolution and structural transformations constitutes a prominent research frontier for ultimately controlling the electrochemical properties of advanced battery materials. Herein, for the first time, we utilize in situ high-energy Kr ion irradiation with transmission electron microscopy to monitor how defects and microstructures evolve in Na- and Li-layered cathodes with 3d transition metals. Our experimental and theoretical analyses reveal that Li-layered cathodes are more resistant to radiation-induced structural transformations, such as amorphization than Na-layered cathodes. The underlying mechanism is the facile formation of Li-transition metal antisite defects in Li-layered cathodes. The quantitative mathematical analysis of the dynamic bright-field imaging shows that defect clusters preferentially align along the Na/Li ion diffusion channels ( a-b planes), which is likely governed by the formation of dislocation loops. Our study provides critical insights into designing battery materials for extreme irradiation environments and understanding fundamental defect dynamics in layered oxides.

25 ENERGY STORAGE↗
DOE OSTI2025

Simple analytic fusion hot spot models for fusion reaction history

The measured fusion reaction history is a combination of the temporal evolution of the fusion hot spot temperature, mass, and volume. Depending on the mechanism of evolution in inertial confinement fusion implosions—shocks, compression, convergence, mass ablation, ignition—the evolution of the reaction history varies. Here, we derive and catalog a set of simplified inertial confinement fusion hot spot models with analytic solutions to infer the evolution of the fusion reaction history for each mechanism. The models give valuable insight into the meaning and cause of fusion reaction history and nuclear burnwidth measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Arctic soil patterns analogous to fluid instabilities

Significance Slow-moving arctic soils form patterns resembling those found in common fluids, such as paint and cake icing drips. Inspired by fluid instabilities, we develop a conceptual model for soil patterns and use mathematical analysis to predict their wavelength. In particular, we propose that soil patterns arise due to competition between gravity and cohesion or the “stickiness” of soil grains. We compare our theoretical predictions with a dataset of soil features from Norway, finding that soil patterns are controlled by fluid-like properties as well as climate. Our work provides a physical explanation for a common pattern on both Earth and Mars, with implications for our understanding of landscapes and complex materials composed of both granular and fluid components.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Enhanced beam-beam modeling to include longitudinal variation during weak-strong simulation

Beam-beam interactions pose substantial challenges in the design and operation of circular colliders, significantly affecting their performance. In particular, the weak-strong simulation approach is pivotal for investigating single-particle dynamics during the collider design phase. This paper evaluates the limitations of existing models in weak-strong simulations, noting that while they accurately account for energy changes due to slingshot effects, they fail to incorporate longitudinal coordinate changes ( z variation). To address this gap, we introduce two novel transformations that enhance Hirata’s original framework by including both z variation and slingshot effect-induced energy changes. Through rigorous mathematical analysis and extensive weak-strong simulation studies, we validate the efficacy of these enhancements in achieving a more precise simulation of beam-beam interactions. Our results reveal that although z variation constitutes a higher-order effect and does not substantially affect the emittance growth rate within the specific design parameters of the Electron-Ion Collider, the refined model offers improved accuracy, particularly in scenarios involving the interaction between beam-beam effects and other random diffusion processes, as well as in simulations incorporating realistic lattice models. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this Letter, we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between nonorthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. Further, while eigenvector continuation is a nonperturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Reconfigurable Intelligent Surfaces Empowered THz Communication in LEO Satellite Networks

Massive swarms of low Earth orbit (LEO) satellites are poising to serve for high-speed and low-latency ubiquitous connectivity with almost global coverage. Broadband inter-satellite communication is one of the key elements of satellite communication systems that orchestrate massive satellite swarms in cooperation. Thanks to technological advancements in microelectronics and micro-systems, the terahertz (THz) band has emerged as a strong candidate for inter-satellite links (ISLs) due to its promise of wideband communication. Especially, multi-antenna systems can improve the system performance along with the wideband supported by the THz band. However, multi-antenna systems should be reconsidered due to their size, weight, and price/power (SWaP) constraints. On the other hand, as a state-of-art multi-antenna technology, reconfigurable intelligent surface (RIS) is able to relax SWaP constraints because of its passive component-based structures. However, as similar reflection characteristic throughout wideband is difficult to meet, it is possible to observe beam misalignment. In this work, we first provide an assessment of the use of the THz band for ISLs and quantify the impact of misalignment fading on error performance. Then, in order to compensate for the high path loss associated with high carrier frequencies, and to further improve the signal-to-noise ratio (SNR), we propose the use of RISs mounted on neighboring satellites to enable signal propagation. Based on a mathematical analysis of the problem, we present the error rate expressions for RIS-assisted ISLs with misalignment fading. Also, numerical results show that RIS can leverage the error rate performance and achievable capacity of THz ISLs.

reconfigurable intelligent surfaces (RISs)↗
DOE OSTI2023

Gradient-Descent-Based Learning Gain for Backstepping Controller and Disturbance Observer of Nonlinear Systems

This paper proposes a gradient-descent-based learning (GL) gain for backstepping controller and disturbance observer (DOB) of nonlinear system. The proposed method consists of the GL gain update law, controller, and DOB. The GL gain update law is proposed to adapt the control gain and DOB gain according to the direction that minimizes the cost function. The mathematical analysis reveals that the GL gain always has a positive sign and upper bound. The controller is designed via a backstepping procedure to track the desired output with GL control gain. The DOB is designed to estimate the unknown external disturbance with the GL DOB gain. Because the control and DOB gains are simultaneously tuned to achieve improved performance, the time consumption for tuning can be reduced. In addition, the peaking phenomenon can be avoided initially by a small initial value of GL gains. The stability of the closed-loop system is guaranteed using the input-to-state stability property. The performance of the proposed method was validated via simulations and experiments using a DC motor.

42 ENGINEERING↗
DOE OSTI2024

Structural Impact of Grid-Forming Inverters on Power System Coherency

This paper addresses the following fundamental research question: how does the integration of grid-forming inverters (GFMs) replacing conventional synchronous generators (SGs) impact the slow coherent eigen-structure and the low frequency oscillatory behavior of future power systems? Due to time-scale separated dynamics, generator states inside a coherent area synchronize over a fast time-scale due to stronger coupling, while the areas themselves synchronize over a slower time scale. Our mathematical analysis shows that due to the large-scale integration of GFMs, the weighted Laplacian structure of the frequency dynamics is preserved, however, the entries of the Laplacian may be significantly modified based on the location and penetration levels of the GFMs. This can impact and potentially significantly alter the coherency structure of the system. We have validated our findings with numerical results using the IEEE 68-bus test system.

Mukherjee, Sayak [BATTELLE (PACIFIC NW LAB)]↗
DOE OSTI2025

Flexibility Options: A Proposed Product for Managing Imbalance Risk

The presence of variable renewable energy resources with uncertain outputs in day-ahead electricity markets results in additional balancing needs in real-time. Addressing those needs cost-effectively and reliably within a competitive market with unbundled products is challenging as both the demand for and the availability of flexibility depends on day-ahead energy schedules. Existing approaches for reserve procurement usually rely either on oversimplified demand curves that do not consider how system conditions that particular day affect the value of flexibility, or on bilateral trading of hedging instruments that are not co-optimized with day-ahead schedules. This article proposes a new product, ‘Flexibility Options', to address these two limitations. The demand for this product is endogenously determined in the day-ahead market and it is met cost-effectively by considering real-time supply curves for product providers, which are co-optimized with the energy supply. As we illustrate with numerical examples and mathematical analysis, the product addresses the hedging needs of participants with imbalances cost-effectively, provides a less intermittent revenue stream for participants with flexible outputs, promotes value-driven pricing of flexibility, and ensures that the system operator is revenue-neutral. This article provides a comprehensive design that can be further tested and applied in large-scale systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Fermilab Booster Beam Emittances from Quadrupole Modes Measured BPMs

The measurement of beam emittances by extracting the quadrupole mode signal from a 4 plate beam position monitor (BPM) was published at least 40 years ago. Unfortunately, in practice, this method suffers from poor signal to noise ratio and requires a lot of tuning to extract out the emittances. In this paper, an improved method where multiple BPMs are used together with better mathematical analysis is described. The BPM derived emittances are then compared with those measured by the Ion Profile Monitor (IPM). Surprisingly, the BPM measured emittances behave very well and are more realistic than those measured by the IPM.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

An Optimal Power Control Strategy for Grid-Following Inverters in a Synchronous Frame

This work proposes a power control strategy based on the linear quadratic regulator with optimal reference tracking (LQR-ORT) for a three-phase inverter-based generator (IBG) using an LCL filter. The use of an LQR-ORT controller increases robustness margins and reduces the quadratic value of the power error and control inputs during transient response. A model in a synchronous reference frame that integrates power sharing and voltage–current (V–I) dynamics is also proposed. This model allows for analyzing closed-loop eigenvalue location and robustness margins. The proposed controller was compared against a classical droop approach using proportional-resonant controllers for the inner loops. Mathematical analysis and hardware-in-the-loop (HIL) experiments under variations in the LCL filter components demonstrate fulfillment of robustness and performance bounds of the LQR-ORT controller. Experimental results demonstrate accuracy of the proposed model and the effectiveness of the LQR-ORT controller in improving transient response, robustness, and power decoupling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Modeling, Load Profile Validation, and Assessment of Solar-Rooftop Energy Potential for Low-and-Moderate-Income Communities in the Caribbean

This document presents the modeling of load profile consumption for Low-and-Moderate-Income (LMI) communities in the Caribbean Islands, as well as an assessment of the solar-rooftop energy potential. In this work, real data, together with synthetic and electricity bill data, were collected to validate and improve the load profile models. The solar-rooftop energy potential was obtained through a National Renewable Energy Laboratory (NREL) software called the PVWatts calculator, and mathematical analysis. The analysis of rooftop solar energy potential was conducted to enable the minimum size of solar power systems to fit the energy demand in the community. The results obtained allow estimation of the capacity of the energy system for each house or an entire community.

14 SOLAR ENERGY↗
DOE OSTI2021

The impact of entrained air on ocean waves

We make a physical–mathematical analysis of the implications that the presence of a large number of tiny bubbles may have, when present, on the thin upper layer of the sea. In our oceanographic example, the bubbles are due to intense rain. It was found that the bubbles increase momentum dissipation in the near surface and affect the surface tension force. For short waves, the implications of increased vorticity are momentum exchanges between wave and mean flow and modifications to the wave dispersion relation. For the direct effect we have analyzed, the implications are estimated to be non-significant when compared to other processes of the ocean. However, we hint at the possibility that our analysis may be useful in other areas of research or practical application.

58 GEOSCIENCES↗