Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation 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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRS1989

Galaxy formation by dust

It has been known since the early 1940's that radiation can cause an instability in the interstellar medium. Absorbing dust particles in an isotropic radiation field shadow each other by a solid angle which is inversely proportional to the square of the distance between the two particles, leading to an inverse-square attractive force - mock gravity. The effect is largest in an optically thin medium. Recently Hogan and White (HW, hereafter) proposed that if the pre-galactic universe contained suitable sources of radiation and dust, instability in the dust distribution caused by mock gravity may have led to the formation of galaxies and galaxy clusters. In their picture of a well-coupled dust-gas medium, HW show that mock gravity begins to dominate gravitational instability when the perturbation becomes optically thin, provided that the radiation field at the time is strong enough. The recent rocket observation of the microwave background at submillimeter wavelengths by Matsumoto et al. might be from pre-galactic stars, the consequence of the absorption of ultraviolet radiation by dust, and infrared reemission which is subsequently redshifted. HW's analysis omits radiative drag, incomplete collisional coupling of gas and dust, finite dust albedo, and finite matter pressure. These effects could be important. In a preliminary calculation including them, the authors have confirmed that mock gravitational instability is effective if there is a strong ultraviolet radiation at the time, but any galaxies that form would be substantially enriched in heavy elements because the contraction of the dust is more rapid than that of the gas. Moreover, since the dust moves with supersonic velocity through the gas soon after the perturbation becomes optically thin, the sputtering of dust particles by gas is significant, so the dust could disappear before the instability develops significantly. They conclude that the mock gravity by dust is not important in galaxy formations.

Wang, Boqi↗
NASA NTRS1990

Computational analysis of flow in 3D propulsive transition ducts

A numerical analysis of fully three dimensional, statistically steady flows in propulsive transition ducts being considered for use in future aircraft of higher maneuverability is investigated. The purpose of the transition duct is to convert axisymmetric flow from conventional propulsion systems to that of a rectangular geometry of high aspect ratio. In an optimal design, the transition duct would be of minimal length in order to reduce the weight penalty, while the geometrical change would be gradual enough to avoid detrimental flow perturbations. Recent experiments conducted at the Propulsion Aerodynamics Branch have indicated that thrust losses in ducts of superelliptic cross-section can be surprisingly low, even if flow separation occurs near the divergent walls. In order to address the objective of developing a rational design procedure for optimal transition ducts, it is necessary to have available a reliable computational tool for the analysis of flows achieved in a sequence of configurations. Current CFD efforts involving complicated geometries usually must contend with two separate but interactive aspects: namely, grid generation and flow solution. The first two avenues of the present investigation were comprised of suitable grid generation for a class of transition ducts of superelliptic cross-section, and the subsequent application of the flow solver PAB3D to this geometry. The code, PAB3D, was developed as a comprehensive tool for the solution of both internal and external high speed flows. The third avenue of investigation has involved analytical formulations to aid in the understanding of the nature of duct flows, and also to provide a basis of comparison for subsequent numerical solutions. Numerical results to date include the generation of two preliminary grid systems for duct flows, and the initial application of PAB3D to the corresponding geometries, which are of the class tested experimentally.

Sepri, Paavo↗
DOE OSTI2022

Pure quark and gluon observables in collinear drop

We construct a class of pure quark and gluon observables by using the collinear drop grooming technique. The construction is based on linear combinations of multiple cumulative distributions of the jet mass in collinear drop, whose specific weights are fully predicted perturbatively. This yields observables which obtain their values purely from quarks (or purely from gluons) in a wide region of phase space. We demonstrate this by showing that these observables are effective in two phase space regions, one dominated by perturbative resummation and one dominated by nonperturbative effects. The nonperturbative effects are included using shape functions which only appear as a common factor in the linear combinations constructed. We test this construction using a numerical analysis with next-to-leading logarithmic resummation and various shape function models, as well as analyzing these observables with PYTHIA and VINCIA. Choices for the collinear drop parameters are optimized for experimental use.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Electron dynamics in extended systems within real-time time-dependent density-functional theory

Due to a beneficial balance of computational cost and accuracy, real-time time-dependent density-functional theory has emerged as a promising first-principles framework to describe electron real-time dynamics. Here we discuss recent implementations around this approach, in particular in the context of complex, extended systems. Results include an analysis of the computational cost associated with numerical propagation and when using absorbing boundary conditions. Here we extensively explore the shortcomings for describing electron–electron scattering in real time and compare to many-body perturbation theory. Modern improvements of the description of exchange and correlation are reviewed. In this work, we specifically focus on the Qb@ll code, which we have mainly used for these types of simulations over the last years, and we conclude by pointing to further progress needed going forward.

2D materials↗
NASA NTRS1990

Numerical experiments in supersonic boundary-layer stability

The three-dimensional (3-D) time-dependent compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation method to study the stability of supersonic flows over a flat plate. Several direct simulations carried out in this study suggest the existence of a secondary instability that might provide a route to transition. The interaction of the modes involved in the secondary instability is possibly amenable to a Floquet-type analysis. Pertinent differences between this instability and the analogous incompressible K-type secondary instability are pointed out. Some preliminary results of a 2-D direct simulation of the nonlinear evolution of a second mode perturbation wave are also discussed.

Erlebacher, Gordon↗
NASA NTRS1990

Low order dynamical models for mixing layers

The goal of this project was to investigate the feasibility of developing useful low order dynamical systems to allow one to determine cheaply the effect of various initial conditions on the evolution of mixing layers. A mathematically rigorous approach for a time evolving base flow as found in a mixing layer seems to be the singular value decomposition of the linear operator transforming some initial perturbations into the solution at a given later time. Unfortunately, this is very costly, if at all feasible, and only deals with the linear evolution of a disturbance. A different approach based on weakly nonlinear analysis was attempted, but numerous difficulties were encountered. Several formulations were proposed which may be ill-posed. The reasons for this possible failure of the attempted approach is related to the fundamental limitations of weakly nonlinear analysis. The clarification of these issues is attempted in this report.

Riahi, D. N.↗
NASA NTRS1992

Pressure oscillation in the leakage annulus between a shrouded impeller and its housing due to impeller-discharge-pressure disturbances

The perturbed flow in the leakage path between a shrouded-pump impeller and its housing is analyzed using experiences with the Space Shuttle Main Engine (SSME), high pressure fuel turbopump (HPFTP) wearing-ring seals. Analysis is based on a bulk-flow model which consists of the path-momentum, circumferential momentum, and continuity equations. The pressure oscillations in the leakage annulus are driven by a circumferential variation of the impeller discharge pressure. It is shown that the occurrence and nature of the pressure oscillations depend on the tangential-velocity ratio of the fluid entering the seal, the order of the Fourier coefficient, the closeness of the precessional frequency of the rotating pressure field to the first natural frequency of the fluid annulus, and the clearance of the wearing-ring seal. The results obtained may explain the internal melting observed on SSME HPFTP seal parts.

Childs, D. W.↗
NASA NTRS1993

Modulational stability of periodic solutions of the Kuramoto-Sivaskinsky equation

We study the long-wave, modulational, stability of steady periodic solutions of the Kuramoto-Sivashinsky equation. The analysis is fully nonlinear at first, and can in principle be carried out to all orders in the small parameter, which is the ratio of the spatial period to a characteristic length of the envelope perturbations. In the linearized regime, we recover a high-order version of the results of Frisch, She, and Thual, which shows that the periodic waves are much more stable than previously expected.

Papageorgiou, Demetrios T.↗
NASA NTRS1997

The Effect of Three-Dimensional Freestream Disturbances on the Supersonic Flow Past a Wedge

The interaction between a shock wave (attached to a wedge) and small amplitude, three-dimensional disturbances of a uniform, supersonic, freestream flow are investigated. The paper extends the two-dimensional study of Duck et al, through the use of vector potentials, which render the problem tractable by the same techniques as in the two-dimensional case, in particular by expansion of the solution by means of a Fourier-Bessel series, in appropriately chosen coordinates. Results are presented for specific classes of freestream disturbances, and the study shows conclusively that the shock is stable to all classes of disturbances (i.e. time periodic perturbations to the shock do not grow downstream), provided the flow downstream of the shock is supersonic (loosely corresponding to the weak shock solution). This is shown from our numerical results and also by asymptotic analysis of the Fourier-Bessel series, valid far downstream of the shock.

Duck, Peter W.↗
NASA NTRS2004

Atmospheric Models for Aerocapture

There are eight destinations in the solar System with sufficient atmosphere for aerocapture to be a viable aeroassist option - Venus, Earth, Mars, Jupiter, Saturn and its moon Titan, Uranus, and Neptune. Engineering-level atmospheric models for four of these targets (Earth, Mars, Titan, and Neptune) have been developed for NASA to support systems analysis studies of potential future aerocapture missions. Development of a similar atmospheric model for Venus has recently commenced. An important capability of all of these models is their ability to simulate quasi-random density perturbations for Monte Carlo analyses in developing guidance, navigation and control algorithm, and for thermal systems design. Similarities and differences among these atmospheric models are presented, with emphasis on the recently developed Neptune model and on planned characteristics of the Venus model. Example applications for aerocapture are also presented and illustrated. Recent updates to the Titan atmospheric model are discussed, in anticipation of applications for trajectory and atmospheric reconstruct of Huygens Probe entry at Titan.

Justus, C. G.↗
NASA NTRS2005

Computational Analysis of the Flow and Acoustic Effects of Jet-Pylon Interaction

Computational simulation and prediction tools were used to understand the jet-pylon interaction effect in a set of bypass-ratio five core/fan nozzles. Results suggest that the pylon acts as a large scale mixing vane that perturbs the jet flow and jump starts the jet mixing process. The enhanced mixing and associated secondary flows from the pylon result in a net increase of noise in the first 10 diameters of the jet s development, but there is a sustained reduction in noise from that point downstream. This is likely the reason the pylon nozzle is quieter overall than the baseline round nozzle in this case. The present work suggests that focused pylon design could lead to advanced pylon shapes and nozzle configurations that take advantage of propulsion-airframe integration to provide additional noise reduction capabilities.

Hunter, Craig A.↗
DOE OSTI2024

Null Raychaudhuri: canonical structure and the dressing time

Abstract We initiate a study of gravity focusing on generic null hypersurfaces, non-perturbatively in the Newton coupling. We present an off-shell account of the extended phase space of the theory, which includes the expected spin-2 data as well as spin-0, spin-1 and arbitrary matter degrees of freedom. We construct the charges and the corresponding kinematic Poisson brackets, employing a Beltrami parameterization of the spin-2 modes. We explicitly show that the constraint algebra closes, the details of which depend on the non-perturbative mixing between spin-0 and spin-2 modes. Finally we show that the spin zero sector encodes a notion of a clock, called dressing time, which is dynamical and conjugate to the constraint. It is well-known that the null Raychaudhuri equation describes how the geometric data of a null hypersurface evolve in null time in response to gravitational radiation and external matter. Our analysis leads to three complementary viewpoints on this equation. First, it can be understood as a Carrollian stress tensor conservation equation. Second, we construct spin-0, spin-2 and matter stress tensors that act as generators of null time reparametrizations for each sector. This leads to the perspective that the null Raychaudhuri equation can be understood as imposing that the sum of CFT-like stress tensors vanishes. Third, we solve the Raychaudhuri constraint non-perturbatively. The solution relates the dressing time to the spin-2 and matter boost charge operators. Finally we establish that the corner charge corresponding to the boost operator in the dressing time frame is monotonic. These results show that the notion of an observer can be thought of as emerging from the gravitational degrees of freedom themselves. We briefly mention that the construction offers new insights into focusing conjectures.

Physics↗
NASA NTRS1979

Stability analysis of a liquid fuel annular combustion chamber

The problems of combustion instability in an annular combustion chamber are investigated. A modified Galerkin method was used to produce a set of modal amplitude equations from the general nonlinear partial differential acoustic wave equation. From these modal amplitude equations, the two variable perturbation method was used to develop a set of approximate equations of a given order of magnitude. These equations were modeled to show the effects of velocity sensitive combustion instabilities by evaluating the effects of certain parameters in the given set of equations. By evaluating these effects, parameters which cause instabilities to occur in the combustion chamber can be ascertained. It is assumed that in the annular combustion chamber, the liquid propellants are injected uniformly across the injector face, the combustion processes are distributed throughout the combustion chamber, and that no time delay occurs in the combustion processes.

Mcdonald, G. H.↗
NASA NTRS1992

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗
DOE OSTI2021

On the Rate Constant for NH 2 +HO 2 and Third-Body Collision Efficiencies for NH 2 +H(+M) and NH 2 +NH 2 (+M)

In low-temperature flash photolysis of NH 3 /O 2 /N 2 mixtures, the NH 2 consumption rate and the product distribution is controlled by the reactions NH 2 + HO 2 → products (R1), NH 2 + H (+M) → NH 3 (+M) (R2), and NH 2 + NH 2 (+M) → N 2 H 4 (+M) (R3). In the present work, published flash photolysis experiments by, among others, Cheskis and co-workers, are re-interpreted using recent direct measurements of NH 2 + H (+N 2 ) and NH 2 + NH 2 (+N 2 ) from Altinay and Macdonald. To facilitate analysis of the FP data, relative third-body collision efficiencies compared to N 2 for R2 and R3 were calculated for O 2 and NH 3 as well as for other selected molecules. We report results were in good agreement with the limited experimental data. Based on reported NH 2 decay rates in flash photolysis of NH 3 /O 2 /N 2 , a rate constant for NH 2 + HO 2 → NH 3 + O 2 (R1a) of $k_{1\text{a}}$ = 1.5(±0.5) × 10 14 cm 3 mol –1 s –1 at 295 K was derived. This value is higher than earlier determinations based on the FP results but in good agreement with recent theoretical work. Kinetic modeling of reported N 2 O yields indicates that NH 2 + HO 2 → H 2 NO + O (R1c) is competing with R1a, but perturbation experiments with addition of CH4 indicate that it is not a dominating channel. Measured HNO profiles indicate that this component is formed directly by NH 2 + HO 2 → HNO + H 2 O (R1b), but theoretical work indicates that R1b is only a minor channel. Based on this analysis, we estimate $k_{1\text{c}}$ = 2.5 × 10 13 cm 3 mol –1 s –1 and $k_{1\text{b}}$ = 2.5 × 10 12 cm 3 mol –1 s –1 at 295 K, with significant uncertainty margins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Uncovering Hazards Using a Multi-Objective Optimization to Explore the Faulty State-Space

Considering resilience when designing complex engineered systems is crucial to ensure the system is safe under unexpected hazardous scenarios. Traditional risk-based approaches, such as Failure Modes and Effects Analysis (FMEA) are useful for designing the system to mitigate hazardous scenarios that can be identified by the designer, but often require experience or prior knowledge of system failures to generate. More recently, researchers have developed simulation tools that enable the designer to model large sets of hazardous scenarios (driven by both internal faults and external factors) through simulation. While these tools enable a wider scope of fault modes to be evaluated (e.g., by injecting combined set of fault modes or injecting modes at different times), the resulting assessments (like FMEA) still require knowledge of the specific modes to be evaluated. However, failure to analyze a wide variety of fault scenarios can lead to an incomplete picture of the system resilience, especially to "surprise events'' which may be difficult for the designer to identify and predict beforehand. To overcome this challenge, previous work developed a fault sampling approach for resilience simulations which would procedurally-generate a wide variety of potential faults by systematically perturbing the health states of the system. While the resulting fault modes generated covered a much larger space hazards than would be otherwise considered (and identified many unique failure trajectories which would not have otherwise been identified), it also significantly increased the computational cost of the analysis and resulted in the simulation and analysis of a large set of essentially duplicate scenarios. Additionally, as the number of dimensions in the faulty state-space increases, the full elaboration of possible modes becomes computationally infeasible, justifying the use of a more targeted search. To resolve this limitation, this work proposes the use of a multiobjective optimization algorithm to search the health state space for potential fault modes that are both (1) hazardous and (2) unique. To solve this type of problem, this work proposes the use of a cooperative co-evolutionary algorithm. To demonstrate this approach, it will be applied to a model of an autonomous rover which uses line markings to navigate, focusing on potential hazards in the drive system which could cause the rover to crash. To determine the merit of the approach, it will further be compared with the previously-presented range elaboration approach and a random mode generation approach on the basis of computational efficiency and found modes.

Resilience↗
DOE OSTI2020

High–level β–carotene production from xylose by engineered Saccharomyces cerevisiae without overexpression of a truncated HMG1 (t HMG1 )

Beta-carotene is a natural pigment and health-promoting metabolite, and has been widely used in the nutraceutical, feed, and cosmetic industries. Here, we engineered a GRAS yeast Saccharomyces cerevisiae to produce beta-carotene from xylose, the sec- ond most abundant and inedible sugar component of lignocellulose biomass. Specifically, a beta-carotene biosynthetic pathway containing crtYB, crtI, and crtE from Xanthophyllomyces dendrorhous was introduced into a xylose-fermenting S. cerevisiae. Here, the resulting strain produced beta-carotene from xylose at a titer threefold higher than from glucose. Interestingly, overexpression of tHMG1, which has been reported as a critical genetic perturbation to enhance metabolic fluxes in the mevalonate pathway and beta-carotene production in yeast when glucose is used, did not further improve the production of beta-carotene from xylose. Through fermentation profiling, metabolites analysis, and transcriptional studies, we found the advantages of using xylose as a carbon source, instead of glucose, for beta-carotene production to be a more respiratory feature of xylose consumption, a larger cytosolic acetyl!CoA pool, and an upregulated expression level of rate-limiting genes in the beta-carotene!producing pathway, including ACS1 and HMG1. As a result, 772.8 mg/L of beta-carotene was ob- tained in a fed-batch bioreactor culture with xylose feeding. Considering the in- evitable large scale production of xylose when cellulosic biomass!based bioeconomy is implemented, our results suggest xylose utilization is a promising strategy for overproduction of carotenoids and other isoprenoids in engineered S. cerevisiae.

09 BIOMASS FUELS↗