Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Robust stability”

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 1,171 records · Page 65

Stability of a tokamak plasma with diffuse toroidal rotation

The present work considers the stability of a high- $\beta$ , large aspect ratio, circular plasma with diffuse profiles for the safety factor and the angular toroidal frequency (López & Guazzotto, Phys. Plasmas , vol. 24, 032501). An application of the Frieman–Rotenberg formalism results in a system of scalar eigenmode equations whose coupling is retained at the plasma–vacuum transition but is disregarded across the plasma column, which is a standard practice. The solution technique consists of a multidimensional shooting method for the poloidal harmonics; robust initial guesses are constructed by solving the dispersion relation in the static scenario with vanishing magnetic shear. Flow shear appears as a high- $\beta$ toroidal contribution, and we illustrate its destabilizing influence on $n=1$ external kink modes in the presence of ideal and resistive walls. Internal resonances are avoided by means of the selection of appropriate equilibrium parameters. The stabilizing influence of a finite positive average magnetic shear is also exemplified.

Physics↗

Deep Reinforcement Learning-Based Control of Energy Storage for Interarea Oscillation Damping

With the increasing electricity consumption and lack of transmission investment, today's power systems are operated much closer to their limits, raising concerns of inter-area oscillations that deteriorate the system stability. Here, this article presents a novel energy storage placement and control approach for enhanced damping of interarea oscillations. Combining the residual analysis and dominant mode analysis, we are able to identify the advantageous locations for placing energy storage that achieve improved damping performance. To overcome the challenges, such as fixed control parameters and insufficient damping, we propose to use a deep reinforcement learning-based approach for energy storage control. A state-of-the-art guided surrogate-gradient-based evolutionary strategy is used to train a learning agent in a robust, efficient, and reproducible manner. Parallel computing is also adopted to speed up the training process. The proposed strategy has been tested on both medium and large-scale systems. The proposed methods have demonstrated their effectiveness in mitigating various interarea oscillations within a timeframe of 20 s, thereby averting system collapse and enhancing power grid stability effectively.

25 ENERGY STORAGE↗

Long-range magnetic order and relaxor ferroelectricity in a hexagonal high-entropy ferrite

Multiferroics that combine ferroelectricity and magnetic order are attractive for electronic and spintronic technologies, yet chemical disorder that promotes relaxor ferroelectricity usually suppresses long-range magnetic order. Here, we report entropy-stabilized relaxor multiferroicity in epitaxial hexagonal (Tb0.2Dy0.2Ho0.2Lu0.2Yb0.2)FeO3 thin films. Structural, magnetic, dielectric, and synchrotron spectroscopic measurements show the coexistence of relaxor ferroelectricity and long-range ferromagnetic order. We find that improper ferroelectricity remains robust against A-site configurational disorder, while the Fe sublattice preserves magnetic exchange. This separation of the microscopic origins of the polar and magnetic responses enables chemically disordered multiferroicity. Our results establish entropy engineering in hexagonal ferrites as a route toward multifunctional oxide thin films and provide a general design strategy for high-entropy multiferroics.

Miertschin, Duncan [Baylor University]↗

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗

Entropy Stable Wall Boundary Conditions for the Three-Dimensional Compressible Navier-Stokes Equations

Non-linear entropy stability and a summation-by-parts framework are used to derive entropy stable wall boundary conditions for the three-dimensional compressible Navier-Stokes equations. A semi-discrete entropy estimate for the entire domain is achieved when the new boundary conditions are coupled with an entropy stable discrete interior operator. The data at the boundary are weakly imposed using a penalty flux approach and a simultaneous-approximation-term penalty technique. Although discontinuous spectral collocation operators on unstructured grids are used herein for the purpose of demonstrating their robustness and efficacy, the new boundary conditions are compatible with any diagonal norm summation-by-parts spatial operator, including finite element, finite difference, finite volume, discontinuous Galerkin, and flux reconstruction/correction procedure via reconstruction schemes. The proposed boundary treatment is tested for three-dimensional subsonic and supersonic flows. The numerical computations corroborate the non-linear stability (entropy stability) and accuracy of the boundary conditions.

Parsani, Matteo↗

Transient Stability-Enhancing Method for Grid-Forming Inverters Under Current Limiting

The primary controller of a grid-forming (GFM) inverter governs the internal reference voltage and angle, which enables the inverter to naturally synchronize and share power with the connected grid; however, during disturbances on the grid caused by, e.g., line faults, overloading, or frequency and phase shifts, a GFM inverter can be pushed into a current-limited operation to protect the device against thermal damage. This can lead to a windup of the internal reference frequency generated by the primary controller and cause a total loss of synchronism with the grid. Here, to make the primary controller of GFM inverters more robust against transient instability, this paper proposes a method that manipulates the reference signals governed by the primary controller by introducing the concept of fictitious power. The fictitious power counteracts the integrator windup in the primary controller, which enables the current-limited inverter to conceive a stable operating point. This method is benchmarked against five state-of-the-art GFM inverter stabilization methods to evaluate its stability-enhancing performance. Full-order electromagnetic-transient simulations and experimental results of a hardware inverter test setup validate the benefits of the method to improve the transient stability of GFM inverters in the face of grid disturbances.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗

Simultaneous removal of C 2 H 2 and C 2 H 6 for C 2 H 4 purification by robust MOFs featuring a high density of heteroatoms

Simultaneous removal of C 2 H 6 and C 2 H 2 from C 2 H 4 streams is of great importance in the petrochemical industry but remains a challenging task. To address this challenge, we have selected three isoreticular MOFs with high stability, low cost, and desirable scale-up ability, namely, MOF-303, MIL-160, and CAU-23 and assessed their potential in simultaneous removal of acetylene and ethane for ethylene purification. Each MOF exhibits desirable C 2 H 2 and C 2 H 6 uptake capacity (>5.5 mmol g –1 and >4 mmol g –1 , respectively), as well as good C 2 H 2 /C 2 H 4 selectivity (>2) and C 2 H 6 /C 2 H 4 selectivity (>1.5). Notably, MOF-303 takes up 4.96 mmol g –1 C 2 H 6 at 298 K and 1 bar, the highest value among the three MOFs, with C 2 H 6 /C 2 H 4 selectivity in the range of 1.55–2.47. MIL-160 possesses a very high C 2 H 2 uptake (9.1 mmol g –1 ) and C 2 H 2 /C 2 H 4 selectivity, 10.6 (1 : 1, v/v), at 298 K, much higher than those of all other MOFs tested to date for simultaneous removal of C 2 H 6 and C 2 H 2 from C 2 H 4 . The results from breakthrough experiments confirm that all three MOFs demonstrate excellent performance for C 2 H 4 purification in a ternary mixture of C 2 H 6 /C 2 H 4 /C 2 H 2 (1 : 1 : 1, v/v/v). Here, for MOF-303, MIL-160, and CAU-23, polymer-grade C 2 H 4 up to 0.164, 0.21, and 0.181 mmol g –1 can be obtained from the equimolar ternary mixture in a single separation step from the breakthrough experiment. Additionally, DFT calculations have been performed to further investigate the mechanism of adsorption/separation for C 2 H 6 , C 2 H 4 , and C 2 H 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1], [2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems with moderate coupling strength. They include a stiff linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1], [2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

97 MATHEMATICS AND COMPUTING↗

A flexible and scalable scheme for mixing computed formation energies from different levels of theory

Abstract Computational materials discovery efforts are enabled by large databases of properties derived from high-throughput density functional theory (DFT), which now contain millions of calculations at the generalized gradient approximation (GGA) level of theory. It is now feasible to carry out high-throughput calculations using more accurate methods, such as meta-GGA DFT; however recomputing an entire database with a higher-fidelity method would not effectively leverage the enormous investment of computational resources embodied in existing (GGA) calculations. Instead, we propose here a general procedure by which higher-fidelity, low-coverage calculations (e.g., meta-GGA calculations for selected chemical systems) can be combined with lower-fidelity, high-coverage calculations (e.g., an existing database of GGA calculations) in a robust and scalable manner. We then use legacy PBE(+ U ) GGA calculations and new r 2 SCAN meta-GGA calculations from the Materials Project database to demonstrate that our scheme improves solid and aqueous phase stability predictions, and discuss practical considerations for its implementation.

36 MATERIALS SCIENCE↗

High-throughput computation of electric polarization in solids via Berry flux diagonalization

Electric polarization in the absence of an externally applied electric field is a key property of polar materials, but the standard interpolation-based ab initio approach to compute polarization differences within the modern theory of polarization presents challenges for automated high-throughput calculations. Berry flux diagonalization [J. Bonini et al., Phys. Rev. B 102, 045141 (2020)] has been proposed as an efficient and reliable alternative, though it has yet to be widely deployed. Here, we assess Berry flux diagonalization using ab initio calculations of a large set of materials, introducing and validating heuristics that ensure branch alignment with a minimal number of intermediate interpolated structures. Our automated implementation of Berry flux diagonalization succeeds in cases where prior interpolation-based workflows fail due to band-gap closures or branch ambiguities. Benchmarking with ab initio calculations of 176 candidate ferroelectrics, we demonstrate the efficacy of the approach on a broad range of insulating materials and obtain accurate effective polarization values with fewer interpolated structures than prior automated interpolation-based workflows. Our real-space heuristics that can predict gauge stability a priori from ionic displacements enable a general automated framework for reliable polarization calculations and efficient high-throughput screening of chemically and structurally diverse polar insulators. These results establish Berry flux diagonalization as a robust and efficient method to compute the effective polarization of solids and to accelerate the data-driven discovery of functional polar materials.

Poteshman, Abigail N. [University of Chicago, IL (↗

Squeezing the periodicity of Néel-type magnetic modulations by enhanced Dzyaloshinskii-Moriya interaction of 4d electrons

In polar magnets, such as GaV 4 S 8 , GaV 4 Se 8 and VOSe 2 O 5 , modulated magnetic phases namely the cycloidal and the Néel-type skyrmion lattice states were identified over extended temperature ranges, even down to zero Kelvin. Our combined small-angle neutron scattering and magnetization study shows the robustness of the Néel-type magnetic modulations also against magnetic fields up to 2 T in the polar GaMo 4 S 8 . In addition to the large upper critical field, enhanced spin-orbit coupling stabilize cycloidal, Néel skyrmion lattice phases with sub-10 nm periodicity and a peculiar distribution of the magnetic modulation vectors. Moreover, we detected an additional single-q state not observed in any other polar magnets. Thus, our work demonstrates that non-centrosymmetric magnets with 4d and 5d electron systems may give rise to various highly compressed modulated states.

36 MATERIALS SCIENCE↗

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Single-Phase to Split-Phase Inverters with Advanced Grid Support Functions for Grid-Interactive Applications

This work presents a cost-effective single-phase to split-phase inverter with a reduced switch count, achieving grid interactive performance while maintaining operational efficiency. The proposed system integrates an Andronov-Hopf oscillator based secondary controller, which inherently embeds a nonlinear resistive droop architecture, ensuring rapid dynamic response. A Lyapunov energy function-based primary control enhances transient stability and regulation, while an internal model-based point of common coupling voltage estimation enables cost optimization without additional sensors. Equipped with advanced grid support functionalities, the inverter facilitates seamless distribution system operation with enhanced robustness. The effectiveness of the proposed architecture and control strategy is validated through MATLAB/Simulink and PLECS simulations, demonstrating its feasibility for high-performance grid-supportive applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Nonlinear Fuzzy Hybrid Control of Spacecraft

This paper describes a new approach for intelligent control of a spacecraft with large angle maneuvers. This new approach, based on fuzzy logic, determines the required torque to achieve a robust, high performance attitude response. This scheme extends the robustness, performance and portability of the existing linear or nonlinear attitude controllers. Formulations are presented for attitude-control but can be extended to other applications. A simulation study, which uses the new control strategy to stabilize the motion of the Microwave Anisotropy Probe spacecraft in the presence of disturbances and saturations, demonstrates the merits of the proposed scheme.

Mason, Paul A. C.↗

The NASA Constellation University Institutes Project: Thrust Chamber Assembly Virtual Institute

This paper documents key aspects of the Constellation University Institutes Project (CUIP) Thrust Chamber Assembly (TCA) Virtual Institute (VI). Specifically, the paper details the TCA VI organizational and functional aspects relative to providing support for Constellation Systems. The TCA VI vision is put forth and discussed in detail. The vision provides the objective and approach for improving thrust chamber assembly design methodologies by replacing the current empirical tools with verified and validated CFD codes. The vision also sets out ignition, performance, thermal environments and combustion stability as focus areas where application of these improved tools is required. Flow physics and a study of the Space Shuttle Main Engine development program are used to conclude that the injector is the key to robust TCA design. Requirements are set out in terms of fidelity, robustness and demonstrated accuracy of the design tool. Lack of demonstrated accuracy is noted as the most significant obstacle to realizing the potential of CFD to be widely used as an injector design tool. A hierarchical decomposition process is outlined to facilitate the validation process. A simulation readiness level tool used to gauge progress toward the goal is described. Finally, there is a description of the current efforts in each focus area. The background of each focus area is discussed. The state of the art in each focus area is noted along with the TCA VI research focus in the area. Brief highlights of work in the area are also included.

Tucker, P. Kevin↗

Spacecraft Modeling, Attitude Determination, and Control: Quaternion-Based Approach

Spacecraft attitude determination and control is an important part of a spacecraft to achieve its designed mission. As of today, many spacecrafts have been successfully launched, and most of them have performed well as they were designed. Many research papers have been published to address the attitude determination and control design problems. Several text books are available for students to learn the technology and for engineers to use as references. The most popular spacecraft models for attitude determination algorithms and control design methods are the Euler angle models and the quaternion models . The Euler angle models have been proved very efficient as the linearized models are controllable, and all standard linear control system design methods are directly applicable. The drawbacks related to the Euler angle methods are (a) the designs based on linearized models may not globally stabilize the original nonlinear spacecraft, i.e., the design may not work when the attitude of the spacecraft is far away from the point where the linearization is performed; (b) the models depend on the rotational sequences, this can be error prone if several teams work on the same project and they use different rotational sequences; (c) for any rotational sequence, there is a singular point where the model is not applicable; and (d) since most attitude determination methods use quaternion to represent the spacecraft attitude, there is a need to transform quaternion into Euler angles. On the other hand, for quaternion models, people have found controllers that can globally stabilize nonlinear spacecraft systems; the models do not depend on rotational sequences and they have no singular point; and the quaternion is provided by attitude determination system and ready to use. The main problem with the quaternion model based control system design is that the linearized quaternion model is not controllable. Therefore, most published design methods heavily rely on Lyapunov functions for the nonlinear spacecraft system. But there is no systematic way to obtain a desired Lyapunov functions. Moreover, the Lyapunov function based designs focus on the closed-loop system stability but pay little attention to the closed-loop system performance. In a series of papers, the author proposed some reduced quaternion models which lead to some controllable linearized spacecraft models. Therefore, all standard linear system theory can be directly applied to analyze and design the spacecraft control systems. We showed that, in some cases, the designed control system is not only optimal for the linearized system, but also globally stabilize the original nonlinear system . Clearly, the reduced quaternion models do not depend on rotational sequences. Due to the special structure of the linearized spacecraft model, some most important design methods, such as LQR design and robust pole assignment design are very simple, enjoy the analytical solutions for some problems, have direct connection to the performance measures, such as settling time, rising time , and percentage of overshoot . All these features are attractive for high quality control system designs. The idea mentioned above is then extended to more spacecraft control problems using specific actuators such as magnetic torque bars and control momentum gyroscopes. These types of actuators may not provide exactly desired torques. Most existing methods use different conversions to get approximate solutions, meaning that these actuators may generate a torque close to but not equal to the desired one. Using the reduced quaternion models that incorporate the actuators into the system model, the control inputs are not torques but the operational parameters. The main benefit of this idea is that the control actions are not approximate but accurate. As all actuators have their operational limit, design with input constraints are also considered in this book by using recently developed interior-point optimization techniques. This book grows up from my research on the spacecraft attitude determination and control design methods in more than a decade which is focused on using reduced quaternion models because of their merits stated above. The book provides all necessary background materials on orbital dynamics, rotations and quaternion, frequently used reference frames, transformations between reference frames, space environment and disturbance torques, ephemeris astronomical vector calculations and measurement instruments, spacecraft control actuators and their models, so that the readers will get a global picture and can apply all these information into the spacecraft system modeling, attitude determination, and spacecraft control system designs, which is the main purpose of this book. This book is different from existing books in that we focus on quaternion based spacecraft control system designs and we consider only attitude control system design related problems, from spacecraft modeling, to attitude determination and estimation, to control system design method selection, to control algorithm development, and to the simulation of the control system designs. Moreover, this book addresses different attitude control tasks in the spacecraft life cycle, including spacecraft maneuver, orbit raising, attitude control, and rendezvous. Finally, this book emphasizes the state space design methods rather than the classical frequency design methods.

Attitude determination and control↗

Enabling Structure/Interface Regulation for High Performance Ni-Rich Cathodes

Here, further commercialization of Ni-rich layered cathodes is hindered by severe structure/interface degradation and kinetic hindrance that occur during electrochemical operation, which leads to safety risks and reduced range in electric vehicles (EVs). Herein, by selecting elements with different solubility properties, a multifunctional strategy that synchronously fabricates perovskite-type SrZrO 3 coating and Sr/Zr co-doping is employed to strengthen the structure/interface stability and the Li + transport mobility of LiNi 0.85 Co 0.10 Mn 0.05 O 2 (NCM). Perovskite-type SrZrO 3 protective layers formed on the particle surface can substantially mitigate the unexpected interfacial side reactions and surface phase transitions. In addition, a robust crystal framework is constructed by optimizing local O coordination through the introduction of strong Zr-O bonds. Notably, Li + diffusion kinetics is effectively improved due to expanded cell parameters and O-Li-O slab spacing with the incorporation of large-diameter Sr pillar ions, as revealed by X-ray diffraction. As a result, the Sr/Zr-modified NCM achieves a remarkable capacity retention of 99.4% after 200 cycles at 1 C, and a high rate capacity of 168.9 mAh g -1 at 10 C. This work opens new avenues to develop high-performance NCM cathodes with high energy and high power for EVs with long calendar life.

25 ENERGY STORAGE↗