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 343 records · Page 19

Stress analysis of circular semimonocoque cylinders with cutouts by a perturbation load technique

A method is presented for analyzing the stresses about a cutout in a circular cylinder of semimonocoque construction. The method involves the use of so-called perturbation solutions which are superposed on the stress distribution that would exist in the structure with no cutout in such a way as to give the effects of a cutout. The method can be used for any loading case for which the structure without the cutout can be analyzed and is sufficiently versatile to account for stringer and shear reinforcement about the cutout. (author)

Mccomb, Harvey G , Jr↗

GLOBAL REFERENCE ATMOSPHERIC MODELS FOR AEROASSIST APPLICATIONS

Aeroassist is a broad category of advanced transportation technology encompassing aerocapture, aerobraking, aeroentry, precision landing, hazard detection and avoidance, and aerogravity assist. The eight destinations in the Solar System with sufficient atmosphere to enable aeroassist technology are Venus, Earth, Mars, Jupiter, Saturn, Uranus, Neptune, and Saturn's moon Titan. Engineering-level atmospheric models for five of these targets - Earth, Mars, Titan, Neptune, and Venus - have been developed at NASA's Marshall Space Flight Center. These models are useful as tools in mission planning and systems analysis studies associated with aeroassist applications. The series of models is collectively named the Global Reference Atmospheric Model or GRAM series. An important capability of all the models in the GRAM series is their ability to simulate quasi-random perturbations for Monte Carlo analysis in developing guidance, navigation and control algorithms, for aerothermal design, and for other applications sensitive to atmospheric variability. Recent example applications are discussed.

Duvall, Aleta↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Ascent trajectory dispersion analysis for WTR heads-up space shuttle trajectory

The results of a Space Transportation System ascent trajectory dispersion analysis are discussed. The purpose is to provide critical trajectory parameter values for assessing the Space Shuttle in a heads-up configuration launched from the Western Test Range (STR). This analysis was conducted using a trajectory profile based on a launch from the WTR in December. The analysis consisted of the following steps: (1) nominal trajectories were simulated under the conditions as specified by baseline reference mission guidelines; (2) dispersion trajectories were simulated using predetermined parametric variations; (3) requirements for a system-related composite trajectory were determined by a root-sum-square (RSS) analysis of the positive deviations between values of the aerodynamic heating indicator (AHI) generated by the dispersion and nominal trajectories; (4) using the RSS assessment as a guideline, the system related composite trajectory was simulated by combinations of dispersion parameters which represented major contributors; (5) an assessment of environmental perturbations via a RSS analysis was made by the combination of plus or minus 2 sigma atmospheric density variation and 95% directional design wind dispersions; (6) maximum aerodynamic heating trajectories were simulated by variation of dispersion parameters which would emulate the summation of the system-related RSS and environmental RSS values of AHI. The maximum aerodynamic heating trajectories were simulated consistent with the directional winds used in the environmental analysis.

Source record↗

On the Search for Low W 0

The magnitude of the vacuum expectation value of the Gukov-Vafa-Witten superpotential |W o | plays a central role in the phenomenology of type IIB flux compactifications. Recent analytical constructions have shown that perturbatively flat vacua can be used to obtain very low values of |W o |. We present systematic algorithms to carry out exhaustive numerical searches for such vacua. Here we also analyse them in the statistical context, as part of the entire ensemble of type IIB flux vacua at low |W o |. Our preliminary analysis indicates that these perturbatively flat vacua are statistically sparse in the whole set of vacua at low |W o | as calculated by Denef and Douglas. Two-moduli examples are used to illustrate these more general findings in specific settings. We find that these simple (two moduli) cases are good examples for existence proofs but they do not feature a large statistical tuning freedom for phenomenological applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An ideal MHD δW stability analysis that bypasses the Newcomb equation

In the following work, we demonstrate the efficacy of a Riccati δW ideal MHD stability analysis that bypasses the numerically intractable integration of the Newcomb equation. By transforming the linear Newcomb equation into a quadratic Riccati equation, an accurate and equivalent δW analysis is performed that is shown to enjoy some numerical advantages. Additionally, we demonstrate that the Riccati approach is better conditioned than its Newcomb counterpart at the magnetic axis, and we apply dynamical systems insights to examine its behavior at singular points. We further discuss the constraints involved in solving for the extremal admissible perturbations of a δW analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Towards Realistic and High Fidelity Models for Nuclear Reactor Power Synthesis Simulation with Self-Powered Neutron Detectors

As presented in this report, a weighting function–based inferencing method is being applied to synthesize the power distribution in next-generation and university research reactors based on simulated self power neutron detector (SPND) responses. The overall goal is to assess the impacts of sensor uncertainty and true power distribution perturbations on the error in the synthesized power distribution. Regarding sensor uncertainty, the NuScale Small Modular Reactor (SMR) and the Westinghouse AP1000 serve as testbeds for analyzing the impact of varying the sensor uncertainty, as well as varying the number of sensors per sensor string in the reactor core. The reactor models are informed by Monte Carlo N-Particle (MCNP) neutron flux tallies. For the NuScale SMR and Westinghouse AP1000, the SPND response functions (i.e., the response of the SPNDs to individual segments of fuel) were determined homogeneously. Regarding an analysis of power distribution perturbation detection, the Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor was used as a demonstration case with one particular arrangement of SPNDs; the response functions for this reactor model were determined heterogeneously, making this a uniquely high-fidelity demonstration of perturbation detection. Finally, SPND models generated in the Geometry and Tracking 4 (Geant4) code have been generated and tested for comparison with traditionally implemented analytical SPND models, with the intent for Geant4 integration with the full methodological framework. SPND current outputs as a function of distance from some fuel assembly segment in the NuScale SMR are compared with the analytically determined currents. Results from the sensor uncertainty simulations for the NuScale SMR and AP1000 indicate that the average error in the inferred power distribution on the fuel assembly segment level is reasonably low, being slightly less than the random uncertainty applied to all respective SPNDs in both cores. For example, if all SPNDs in the core have a random uncertainty of 5%, then the corresponding fuel assembly segment level error (i.e. difference between the true and inferred local power) is ~2–3%. However, the maximum error in the inferred power distribution on the fuel assembly segment level can be considerably high (>15%) when SPND random uncertainties start to exceed ~3%. In general, the average and maximum errors in the inferred power distribution were slightly higher in the AP1000 as opposed to the NuScale SMR for the sensor string configurations considered herein. Another result determined from analysis of the sensor uncertainty simulations was that increasing the number of SPNDs per string does not clearly reduce inferred power distribution error and can in fact make the error large in some cases; however, this assessment may skewed due to imposed iteration limits. Results from the perturbation detection demonstration using the high-fidelity TAMU TRIGA model indicate that, given the arrangement of 17 SPND strings and 4 SPNDs per string considered herein, there is a clear, provable ability to infer a localized Gaussian-type peak perturbation in the 3D power distribution. Such a perturbation was detected with an average fuel assembly segment level error of 0.19%, and the general visualization of the detected perturbation clearly indicates that the magnitude and shape were appropriately resolved. Finally, the electrical current output generated by the Geant4 modeled SPND indicates significant magnitude differences than the analytically modeled SPND, demonstrating the need for accurate SPND models which account for finite sensor geometry effects to inform the power synthesis work described herein.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ALESEP: A computer program for the analysis of airfoil leading edge separation bubbles

The ALESEP program for the analysis of the inviscid/viscous interaction which occurs due to the presence of a closed laminar transitional separation bubble on an airflow is presented. The ALESEP code provides a iterative solution of the boundary layer equations expressed in an inverse formulation coupled to a Cauchy integral representation of the inviscid flow. This interaction analysis is treated as a local perturbation to a known solution obtained from a global airfoil analysis. Part of the required input to the ALESEP code are the reference displacement thickness and tangential velocity distributions. Special windward differencing may be used in the reversed flow regions of the separation bubble to accurately account for the flow direction in the discretization of the streamwise convection of momentum. The ALESEP code contains a forced transition model based on a streamwise intermittency function and a natural transition model based on a solution of the integral form of the turbulent kinetic energy equation. Instructions for the input/output, and program usage are presented.

Davis, R. L.↗

ALESEP. Part 2: A computer program for the analysis of leading edge separation bubbles on infinite swept wings

A program called ALESEP is presented for the analysis of the inviscid-viscous interaction which occurs due to the presence of a closed laminar-transitional separation bubble on an airfoil or infinite swept wing. The ALESEP code provides an iterative solution of the boundary layer equations expressed in an inverse formulation coupled to a Cauchy integral representation of the inviscid flow. This interaction analysis is treated as a local perturbation to a known solution obtained from a global airfoil analysis; hence, part of the required input to the ALESEP code are the reference displacement thickness and tangential velocity distributions. Special windward differencing may be used in the reversed flow regions of the separation bubble to accurately account for the flow direction in the discretization of the streamwise convection of momentum. The ALESEP code contains a forced transition model based on a streamwise intermittency function, a natural transition model based on a solution of the integral form of the turbulent kinetic energy equation, and an empirical natural transition model.

Davis, R. L.↗

Quasinormal Coupled-Mode Analysis of Dynamic Gain in Exceptional-Point Lasers

One of the key features of lasers operating near exceptional points (EPs) is that the gain medium can support an oscillating population inversion above a pump threshold, leading to self-modulated laser dynamics. This unusual behavior opens up new possibilities for frequency comb generation and temporal modulation. However, the dynamic population inversion couples signals with different frequencies and is difficult to capture using conventional temporal coupled-mode theory (TCMT) based on stationary saturable gain. In this paper, we develop a perturbative coupled-mode analysis framework to capture the spatial-temporal dynamics of near-EP lasers. By decomposing discrete frequency generation into multiple excitations of resonant modes, our analysis establishes a minimal physical model that translates the local distribution of dynamic population-inversion into a resonant modal interpretation of laser gain. Furthermore, this work enables the exploration of unique properties in this self-time-modulated systems, such as time-varying scattering and nonreciprocal transmission.

Absorption↗

Computational aspects of sensitivity calculations in linear transient structural analysis

A study was performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal of the study was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semi-analytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. In several cases this fixed mode approach resulted in very poor approximations of the stress sensitivities. Almost all of the original modes were required for an accurate sensitivity and for small numbers of modes, the accuracy was extremely poor. To overcome this poor accuracy, two semi-analytical techniques were developed. The first technique accounts for the change in eigenvectors through approximate eigenvector derivatives. The second technique applies the mode acceleration method of transient analysis to the sensitivity calculations. Both result in accurate values of the stress sensitivities with a small number of modes and much lower computational costs than if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, William H.↗

Effective field theory factorization for diffraction

We derive a factorization formula for coherent and incoherent ep diffraction using the soft collinear effective theory, utilizing multiple power expansion parameters to handle different kinematic regions. This goes beyond the known hard-collinear diffractive factorization to address the small-x Regge dynamics and Pomeron exchange from first principles. The effective field theory analysis also uncovers and factorizes an important irreducible incoherent background generated by color-nonsinglet exchange, dubbed “quasi-diffraction”, for which we calculate the associated Sudakov suppression. For unpolarized scattering we show that there are four diffractive structure functions at leading power, and point out the importance of studying ${F}_{3,4}^D$ through asymmetries, in addition to ${F}_{2,L}^D$. For the quasi-diffractive background, we make model independent predictions for ratios of the corresponding structure functions in a perturbative kinematic region. Our analysis also makes predictions for six leading-power spin-dependent structure functions. Finally, we provide connections to diffractive parton distributions, and assess the Ingelman-Schlein model. Our work lays a path for further QCD-based studies of diffraction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Composition Changes After the "Halloween" Solar Proton Event: The High-Energy Particle Precipitation in the Atmosphere (HEPPA) Model Versus MIPAS Data Intercomparison Study

We have compared composition changes of NO, NO2, H2O2,O3, N2O, HNO3 , N2O5, HNO4, ClO, HOCl, and ClONO2 as observed by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on Envisat in the aftermath of the "Halloween" solar proton event (SPE) in October/November 2003 at 25-0.01 hPa in the Northern hemisphere (40-90 N) and simulations performed by the following atmospheric models: the Bremen 2D model (B2dM) and Bremen 3D Chemical Transport Model (B3dCTM), the Central Aerological Observatory (CAO) model, FinROSE, the Hamburg Model of the Neutral and Ionized Atmosphere (HAMMONIA), the Karlsruhe Simulation Model of the Middle Atmosphere (KASIMA), the ECHAM5/MESSY Atmospheric Chemistry (EMAC) model, the modeling tool for SO1ar Climate Ozone Links studies (SOCOL and SOCOLi), and the Whole Atmosphere Community Climate Model (WACCM4). The large number of participating models allowed for an evaluation of the overall ability of atmospheric models to reproduce observed atmospheric perturbations generated by SPEs, particularly with respect to NOS, and ozone changes. We have further assessed the meteorological conditions and their implications on the chemical response to the SPE in both the models and observations by comparing temperature and tracer (CH4 and CO) fields. Simulated SPE-induced ozone losses agree on average within 5% with the observations. Simulated NO(y) enhancements around 1 hPa, however, are typically 30% higher than indicated by the observations which can be partly attributed to an overestimation of simulated electron-induced ionization. The analysis of the observed and modeled NO(y) partitioning in the aftermath of the SPE has demonstrated the need to implement additional ion chemistry (HNO3 formation via ion-ion recombination and water cluster ions) into the chemical schemes. An overestimation of observed H2O2 enhancements by all models hints at an underestimation of the OH/HO2 ratio in the upper polar stratosphere during the SPE. The analysis of chlorine species perturbations has shown that the encountered differences between models and observations, particularly the underestimation of observed ClONO2 enhancements, are related to a smaller availability of ClO in the polar night region already before the SPE. In general, the intercomparison has demonstrated that differences in the meteorology and/or initial state of the atmosphere in the simulations causes a relevant variability of the model results, even on a short timescale of only a few days.

Funke, B.↗

The condensational instability in the solar transition region and corona

The stability of plasmas at temperatures and densities typical of the solar transition region and corona is investigated using both a linear analysis and nonlinear time-dependent numerical simulations. Growth rates, decay rates, and oscillation frequencies of the perturbations determined from the linear analysis are in excellent agreement with the simulations. The nonlinear regime is characterized by a bifurcation of the plasma into a cool dense condensation surrounded by a hot tenuous corona. The condensation may then be accelerated by forces in the plasma such as those arising from gravity or differential heating. Finally, the results of the detailed simulation show that the transition region is a dynamically stable structure which is the result of the nonlinear evolution of the condensational instability.

Oran, E. S.↗

Aeroacoustic sensitivity analysis and optimal aeroacoustic design of turbomachinery blades

During the first year of the project, we have developed a theoretical analysis - and wrote a computer code based on this analysis - to compute the sensitivity of unsteady aerodynamic loads acting on airfoils in cascades due to small changes in airfoil geometry. The steady and unsteady flow though a cascade of airfoils is computed using the full potential equation. Once the nominal solutions have been computed, one computes the sensitivity. The analysis takes advantage of the fact that LU decomposition is used to compute the nominal steady and unsteady flow fields. If the LU factors are saved, then the computer time required to compute the sensitivity of both the steady and unsteady flows to changes in airfoil geometry is quite small. The results to date are quite encouraging, and may be summarized as follows: (1) The sensitivity procedure has been validated by comparing the results obtained by 'finite difference' techniques, that is, computing the flow using the nominal flow solver for two slightly different airfoils and differencing the results. The 'analytic' solution computed using the method developed under this grant and the finite difference results are found to be in almost perfect agreement. (2) The present sensitivity analysis is computationally much more efficient than finite difference techniques. We found that using a 129 by 33 node computational grid, the present sensitivity analysis can compute the steady flow sensitivity about ten times more efficiently that the finite difference approach. For the unsteady flow problem, the present sensitivity analysis is about two and one-half times as fast as the finite difference approach. We expect that the relative efficiencies will be even larger for the finer grids which will be used to compute high frequency aeroacoustic solutions. Computational results show that the sensitivity analysis is valid for small to moderate sized design perturbations. (3) We found that the sensitivity analysis provided important insight into how airfoils should be modified to improve their aeroelastic stability. Using this insight, we redesigned an aeroelastically unstable cascade to produce a stable cascade.

Hall, Kenneth C.↗

Optical Observations of Space Debris

This viewgraph presentation reviews the use of optical telescopes to observe space debris. .It will present a brief review of how the survey is conducted, and what some of the significant results encompass. The goal is to characterize the population of debris objects at GEO, with emphasis on the faint object population. Because the survey observations extend over a very short arc (5 minutes), a full six parameter orbit can not be determined. Recently we have begun to use a second telescope, the 0.9-m at CTIO, as a chase telescope to do follow-up observations of potential GEO debris candidates found by MODEST. With a long enough sequence of observations, a full six-parameter orbit including eccentricity can be determined. The project has used STK since inception for planning observing sessions based on the distribution of bright cataloged objects and the anti-solar point (to avoid eclipse). Recently, AGI's Orbit Determination Tool Kit (ODTK) has been used to determine orbits, including the effects of solar radiation pressure. Since an unknown fraction of the faint debris at GEO has a high area-to-mass ratio (A/M), the orbits are perturbed significantly by solar radiation. The ODTK analysis results indicate that temporal variations in the solar perturbations, possibly due to debris orientation dynamics, can be estimated in the OD process. Additionally, the best results appear to be achieved when solar forces orthogonal to the object-Sun line are considered. Determining the A/M of individual objects and the distribution of A/M values of a large sample of debris is important to understanding the total population of debris at GEO

Seitzer, Patrick↗

3D high-fidelity automated neutronics guided optimization of fusion blanket designs

The compact Fusion Pilot Plant (FPP) is defined in the recent National Academies of Sciences, Engineering, and Medicine report as the next step of fusion energy demonstration with a $50$ MWe peak net electricity production, $Q_e$ greater than $1$, and at least $3$ hours of continuous operation. This fusion pilot plant will be a test bed enabling materials, designs, and fuel management assessment, and it will represent an engineering challenge because of its high-fusion power and compact design targets. Previous reactor data is limited to experiments operating in different design space ranges. Therefore, design iterations and assessments should rely on high-fidelity first-principle theoretical and computational models. The high-fidelity integrated modeling of the plasma is a fundamental part of fusion energy research. However, the whole device modeling is often neglected, utilizing low-fidelity, system-level analysis. Recently, the need for high-fidelity multi-physics modeling was recognized, resulting in a selection of integrated tools. Further, autonomous design optimization requires a streamlined framework that perturbs the design point, reruns the analysis, and examines the outputs. However, high-fidelity analysis requires complex geometry specification that is difficult to perturb. This work presents the parametric CAD generation tool TRACER and a new neutronic workflow. TRACER allows the perturbation of the geometry representation, creating geometry files ready for further analysis. The streamlined neutronic workflow allows efficient and accurate calculations. The two new tools coupled together were used to perform a 3D high-fidelity multi-objective, multi-input optimization of an "ARC Class" compact tokamak design. The workflow was driven by an optimization driver for full automation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of Tibetan Plateau sensible heat on pre-monsoon dust burdens over South Asia: Observational and modeling evidence

Tibetan Plateau sensible heating (TPSH) effect is recognized as a key driver of the South Asian climate; however its role in regulating the regional dust burdens over South Asia during the pre-monsoon season remains unclear. Using long-term reanalysis and/or observational dust and TPSH data, we show that the dust optical depth at 550 nm (DOD 550 ) over South Asia is significantly positively correlated with TPSH at the interannual scale. DOD 550 over South Asia increases (decreases) by 20 % (12 %) in the strongest (weakest) TPSH years. The phenomenon of TPSH-dust relationships is reinforced by sensitivity experiments by altering TPSH from global climate models. Further mechanism analysis reveals that TPSH perturbation induces regional circulation anomalies over South Asia, which features a low-tropospheric cyclonic response around the Tibetan Plateau. This cyclonic anomaly strengthens the prevailing northwesterly winds, increasing regional dust emissions and transportation over South Asia. This novel mechanism helps understand the South Asian dust change.

54 ENVIRONMENTAL SCIENCES↗