SEARCH · Engineering Papers
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.
The eigenvalue method to calculate models for periodically pulsating stars is used to study the stability of the resulting oscillations. This paper discusses results for a series of RR Lyrae models.
A study is made of the thermomechanical buckling and postbuckling responses of flat unstiffened composite panels. The panels are subjected to combined temperature change and applied edge displacement. The analysis is based on a first-order shear deformation, von Karman type nonlinear plate theory. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. An efficient multiple-parameter reduction method is used in conjunction with mixed finite element models, for determining the stability boundary and postbuckling response. The reduction method is also used for evaluating the sensitivity coefficients which measure the sensitivity of the buckling and postbuckling responses to variations in the different lamination and material parameters of the panel. Numerical results are presented showing the effects of variations in the laminate stacking sequence, fiber orientation, number of layers and aspect ratio of the panels on their thermomechanical buckling and postbuckling responses and their sensitivity coefficients.
Here, this paper proposed an energy function-embedded quasi-steady-state model for efficient simulation of cascading outages on a power grid while addressing transient stability concerns. Compared to quasi-steady-state models, the proposed model incorporates short-term dynamic simulation and an energy function method to efficiently evaluate the transient stability of a power grid together with outage propagation without transient stability simulation. Cascading outage simulation using the proposed model conducts three steps for each disturbance such as a line outage. First, it performs time-domain simulation for a short term to obtain a post-disturbance trajectory. Second, along the trajectory, the system state with the local maximum potential energy is found and used as the initial point to search for a relevant unstable equilibrium by Newton's method. Third, the transient energy margin is estimated based on this unstable equilibrium to predict an out-of-step condition with generators. The proposed energy function-embedded quasi-steady-state model is tested in terms of its accuracy and time performance on an NPCC 140-bus power system and compared to a quasi-steady-state model embedding transient stability simulation.
HabEx is a concept study for a next generation space telescope having a 4 to 6.5 m primary mirror diameter. Primarily designed for exoplanet science, but also carrying general astrophysics instruments, it will have stability requirements orders of magnitude more stringent than any previously flown space-based observatory. Wavefront control, including active sensing/control of wavefront jitter and drift, will be essential to achieving and maintaining coronagraph performance. There are currently two complementary techniques that can be used to sense the wavefront drift. An attractive method is to use starlight rejected from the coronagraph to sense changes in the wavefront incident on the coronagraph mask. This method suffers from photon noise especially when imaging dim stars. It is also unavailable during slews between target stars. Laser metrology is a complementary system in which a laser “truss” enables rigid body displacements of the optics to be measured from changes in the separation of pairs of optical fiducials which can be separated by many meters. In closed loop with actuated optics, laser metrology actively maintains rigid body alignment of the front-end optics, thereby eliminating the dominant source of wavefront drift. This method can maintain rigid body stability of the telescope optics and thereby wavefront control, even during slews between targets.
Abstract High-power lasers are vital for particle acceleration, imaging, fusion and materials processing, requiring precise control and high-energy delivery. Laser plasma accelerators (LPAs) demand laser positional stability at focus to ensure consistent electron beams in applications such as X-ray free-electron lasers and high-energy colliders. Achieving this stability is especially challenging for the low-repetition-rate lasers in current LPAs. We present a machine learning method that predicts and corrects laser pointing instabilities in real-time using a high-frequency pilot beam. By preemptively adjusting a correction mirror, this approach overcomes traditional feedback limits. Demonstrated on the BELLA petawatt laser operating at the terawatt level (30 mJ amplification), our method achieved root mean square pointing stabilization of 0.34 and 0.59 $\unicode{x3bc} \mathrm{rad}$ in the x and y directions, reducing jitter by 65% and 47%, respectively. This is the first successful application of predictive control for shot-to-shot stabilization in low-repetition-rate laser systems, paving the way for full-energy petawatt lasers and transformative advances across science, industry and security.
Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.
An adaptive method in which cruise and nonlinear orbit determination problems can be solved using a single program is presented. It involves singular value decomposition augmented with an extended partial step algorithm. The extended partial step algorithm constrains the size of the correction to the spacecraft state and other solve-for parameters. The correction is controlled by an a priori covariance and a user-supplied bounds parameter. The extended partial step method is an extension of the update portion of the singular value decomposition algorithm. It thus preserves the numerical stability of the singular value decomposition method, while extending the region over which it converges. In linear cases, this method reduces to the singular value decomposition algorithm with the full rank solution. Two examples are presented to illustrate the method's utility.
This presents a method to improve the thruster vector servo elasticity (TVSE) stability. TVSE is the coupling between the TVC system engine dynamics and the vehicle flexible modes. It can lead to instability of the flexible modes due to the energy injection from the thrust. For large boosters with heavy TVC engine, TVSE is more likely to cause instability. This instability is part of the plant dynamics and is not detected by the flight control. In this paper, the TVSE is treated as a disturbance rejection problem with the inertial force of the flexible motion viewed as an external disturbance to the TVC actuator control system. By effectively rejecting this disturbance, the engine can be decoupled from the flexible dynamics thus improving TVSE stability. The actuator disturbance transfer function is modified by inserting a TVSE filter into the control system. The filter aims to reduce the magnitude and introduce phase lead of the disturbance transfer function over the TVSE sensitive frequency range while maintaining adequate stability margins of the actuator control system. A prefilter is also designed for good command tracking. The filter and prefilter resulting from this method are usually of low order and can be included to the existing actuator controller if the controller is programmable. This method was demonstrated for a hydraulic Simplex actuator model. Its applications to SLS B1 and B2 TVC actuator are also described. B2 actuator controller is modified based on this method. The tradeoff between TVSE stability and the control system stability can be balanced by adjusting the amount of phase lead introduced by the TVSE filter. This method provides a useful way to solve the TVSE instability issue without costly and time-consuming a TVC system redesign. This method is based on the basic control principles of loop shaping and disturbance rejection and has been proven to be effective under proper circum-stances.
Numerical methods are presented that exploit the polynomial properties of discretizations of the Euler equations. It is noted that most finite difference or finite volume discretizations of the steady-state Euler equations produce a polynomial system of equations to be solved. These equations are solved using classical polynomial elimination theory, with some innovative modifications. This paper also presents some preliminary results of a new non-linear stability analysis technique. This technique is applicable to determining the stability of polynomial iterative schemes. Results are presented for applying the elimination technique to a one-dimensional test case. For this test case, the exact solution is computed in three iterations. The non-linear stability analysis is applied to determine the optimal time step for solving Burgers' equation using the MacCormack scheme. The estimated optimal time step is very close to the time step that arises from a linear stability analysis.
This investigation is divided into two parts, the treatment of the stability problem of fluid flows on the one hand, and that of the turbulent motion on the other. The first part summarizes all previous investigations under a unified point of view, that is, sets up as generally as possible the conditions under which a profile possesses unstable or stable characteristics, and indicates the methods for solution of the stability equation for any arbitrary velocity profile and for calculation of the critical Reynolds number for unstable profiles. In the second part, under certain greatly idealizing assumptions, differential equations for the turbulent motions are derived and from them qualitative information about several properties of the turbulent velocity distribution is obtained.
Hybrid Li-Be ceramics, combining both tritium (T) breeder (Li) and neutron multiplier (Be) for use as part of the fuel cycle of future nuclear fusion reactors are proposed. The development of such hybrid materials may reduce thermal gradients through better matching of thermal properties, mitigating the detrimental effects that may accompany traditional breeder systems while maintaining acceptable neutron multiplication and T breeding. Here, first-principles methods are used to investigate stability and thermal transport properties of a set of compounds containing both Li and Be. It is demonstrated that BeLi 2 O 4 Ge and BeLi 2 O 4 Si are mechanically stable and have formation energies comparable with leading candidates for solid state breeder materials, Li 2 TiO 3 , Li 2 ZrO 3 . It is also demonstrated that similar to the leading candidates, these compounds are insulators with thermal transport defined by phonons. The calculated thermal conductivity of BeLi 2 O 4 Ge is slightly higher compared to Li 2 TiO 3 or Li 2 ZrO 3 while in the BeLi 2 O 4 Si it is almost three times higher compared to the rest of compounds due to a higher phonon group velocities and increased phonon lifetimes. These results indicate that hybrid Li-Be ceramics offer a potential route towards better matching of thermal properties with minimal functional property degradation, thereby offering better overall fuel cycle performance.
The CANDECOMP/PARAFAC (CP) decomposition is widely used for analyzing multidimensional data, and the alternating least squares (CP-ALS) algorithm is a common method for its computation. CP rounding is the problem of computing a lower-rank CP decomposition of an input already in a higher-rank CP format. While the normal equations (NE) approach in CP-ALS is efficient for the CP rounding problem and frequently used, it becomes unstable in the presence of ill-conditioned subproblems. This paper presents a new QR-based CP-ALS method for CP rounding that preserves both numerical stability and computational efficiency. Here, our experiments show that the proposed method offers significant speedup over a previous QR-based approach and the Tensor Toolbox's NE-based implementation, particularly for higher-order tensors. Furthermore, our approach demonstrates a marked reduction in error for ill-conditioned problems, with error reductions several orders of magnitude smaller compared to the NE-based method, while achieving faster convergence and more accurate solutions. By using a more numerically stable approach, we can solve more problems in reduced working precision, which enables further reduction in time to solution.
In the present model of multilayered stability stratification, aimed at obtaining the analytic eigensolutions of the vertical structure equation, each layer is characterized by its own static stability value. By requiring continuity of pressure and vertical velocity across each interface level, and by imposing suitable upper and lower boundary conditions, matching eigensolutions are obtained in terms of the Bessel functions. Attention is given to an explicit example of a double-layered stratified atmosphere which demonstrates the mathematical manipulations involved; the resultant vertical structure functions are used to check the accuracy of the numerical solutions by the finite difference and finite element methods.
The multiple-grid scheme of Ni (1981) for the solution of the unsteady Euler equations for quasi-one-dimensional transonic flow problems is analyzed according to its ability to accelerate convergence to a steady state solution, its applicability to an implicit scheme, its accuracy, and its stability limits. This scheme is applied successfully to McCormack's (1981) implicit method of solving quasi-one-dimensional transonic flow problems for CFL numbers significantly greater than unity. It is determined that the second order accuracy of this implicit method is preserved by employing the multiple gridding. However, stability limitations that are flow dependent are introduced into the implicit method by the use of the multiple-grid scheme which will curtail its use under certain flow conditions. A scaling modification for the multiple grid scheme is developed which further improves its performance when used with explicit and artificially-damped implicit methods without any increase in computational work.
A new and easily used state variable method of aircraft ride quality analysis is developed and its use is illustrated by a numerical example using DC-8 airplane equations of motion. The method readily allows any type of stability augmentation system to be included for good handling qualities. As part of the ride quality analysis, a particularly powerful method of specifying and achieving all roots of a desired closed-loop characteristic equation by use of full state variable feedback is described and applied to the DC-8 example.
Abstract Jefferson Lab operates two central helium liquefiers (CHLs) which both utilize full cold compression from the saturation pressure at operating temperature (approximately 2.1 K) to just over atmospheric pressure. The original plant, CHL1, was recently outfitted with a replacement subatmospheric cold box (SC1R) containing state-of-the-art cold compressor technology, while the newer plant, CHL2, uses an older cold compressor system. In both cases, the heat of compression is absorbed at low temperature at the expense of electrical power consumed by the warm compressors. Due to the superior efficiency and turndown capabilities of SC1R, a new operating mode has been identified for CHL1 in which the required number of operating warm compressors is reduced by one. A cost-based method for optimizing cold compressor stability and efficiency has been developed and applied to CHL2, improving its turndown and lowering the warm compressor discharge pressure. As a result of these efforts, power consumption of the combined CHLs during normal operations has been reduced by nearly 10%, or a total of 650 kW. The observed performance of CHL1 with SC1R, as well as the cold compressor optimization method, leading to this improved energy consumption rate will be discussed in detail.
Perovskite solar cells (PSCs) are a promising alternative solar technology, but the technical challenges of (1) stability/durability, (2) efficiency at scale, and (3) manufacturability must be overcome to achieve widespread PSC commercialization. The challenge of balancing solution ink formulation and scalable manufacturing is often overlooked in the literature, where focus is on adopting inks from processes that will not scale. In this study, we apply a classical roll-to-roll manufacturing perspective, utilizing both compositional engineering and intense pulsed light (IPL) annealing, to develop a mixed-cation perovskite ink with a robust coating window that simultaneously solve issues of stability and manufacturability for PSCs. Our method resulted in blade-coated, flexible, mixed-cation PSCs on ITO-PET substrates with a champion power conversion efficiency (PCE) efficiency of 16.7% using IPL annealing of the absorber layer and, to our knowledge, is one of the fastest processing methods for the perovskite layer. This overall reduction in processing time with a stable ink represents an advance toward the scaled production of perovskite solar cells on flexible substrates.
The algebraic and geometric structures of certain classes of nonlinear stochastic systems were exploited in order to obtain useful stability and estimation results. The class of bilinear stochastic systems (or linear systems with multiplicative noise) was discussed. The stochastic stability of bilinear systems driven by colored noise was considered. Approximate methods for obtaining sufficient conditions for the stochastic stability of bilinear systems evolving on general Lie groups were discussed. Two classes of estimation problems involving bilinear systems were considered. It was proved that, for systems described by certain types of Volterra series expansions or by certain bilinear equations evolving on nilpotent or solvable Lie groups, the optimal conditional mean estimator consists of a finite dimensional nonlinear set of equations. The theory of harmonic analysis was used to derive suboptimal estimators for bilinear systems driven by white noise which evolve on compact Lie groups or homogeneous spaces.