Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
NASA NTRS1987

Perturbation analysis of internal balancing for lightly damped mechanical systems with gyroscopic and circulatory forces

Approximate expressions are developed for internally balanced singular values corresponding to the modes of mechanical systems with gyroscopic forces, light damping, and small circulatory forces. A brief overview is first given of the balanced realization model reduction method, including a discussion of recent work. The models considered are defined, and a perturbation analysis is used to show that the modal representation becomes asymptotically balanced as damping reduces to zero. The approximate balanced singular values are calculated, and a simple example of a flexible, dual-spin spacecraft is given as an illustration of the results.

Blelloch, P. A.↗
NASA NTRS1983

An algorithm for the empirical optimization of antenna arrays

A numerical technique is presented to optimize the performance of arbitrary antenna arrays under realistic conditions. An experimental-computational algorithm is formulated in which n-dimensional minimization methods are applied to measured data obtained from the antenna array. A numerical update formula is used to induce partial derivative information without requiring special perturbations of the array parameters. The algorithm provides a new design for the antenna array, and the method proceeds in an iterative fashion. Test case results are presented showing the effectiveness of the algorithm.

Blank, S.↗
NASA NTRS1975

Steady subsonic flow around finite-thickness wings

The general method for analyzing steady subsonic potential aerodynamic flow around a lifting body having arbitrary shape is presented. By using the Green function method, an integral representation for the potential is obtained. Under small perturbation assumption, the potential at any point, P, in the field depends only upon the values of the potential and its normal derivative on the surface of the body. Hence if the point P approaches the surface of the body, the representation reduces to an integral equation relating the potential and its normal derivative (which is known from the boundary conditions) on the surface. The question of uniqueness is examined and it is shown that, for thin wings, the operator becomes singular as the thickness approaches zero. This fact may yield numerical problems for very thin wings. However, numerical results obtained for a rectangular wing in subsonic flow show that these problems do not appear even for thickness ratio tau = .001. Comparison with existing results shows that the proposed method is at least as fast and accurate as the lifting surface theories.

Kuo, C. C.↗
DOE OSTI2025

Mechanism of the Non‐Kasha Fluorescence in Pyrene

The high-energy shoulder in the gas-phase fluorescence emission spectrum of pyrene is a well-known example of non-Kasha emission. Here, we comparatively assess two approaches, vibronic perturbation theory and nonadiabatic dynamics, in their ability to predict and explain the gas-phase fluorescence spectrum of pyrene. While both methods qualitatively capture the non-Kasha emission, they differ in their computational requirements, accuracy, and physical interpretation. Vibronic perturbation theory and nonadiabatic dynamics are complementary and can be combined in a two-step approach to non-Kasha fluorescence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

An analytical evaluation of mascon effects in a semi-analytic theory

A semianalytic method to predict the orbit of a spacecraft orbiting a planet and perturbed by oblateness and a gravity anomaly (mascon) is presented. The Hamiltonian is first expressed in closed form in terms of the nonsingular canonical set of Poincare elements where the mascon is represented by a point-mass. The planet global potential is limited to the second and third harmonics. The mascon potential in harmonics is limited to second degree and first order. The original Lie-Hori method is then applied to transform the canonical equations of motion from the osculating element space to the mean element space. Numerical integration of the new equations is used.

Salama, Ahmed H.↗
DOE OSTI2020

Sensor Distortion Effects in Photon Monte Carlo Simulations

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three-dimensional electrostatic simulations to parameterize the perturbed electric field profile for nonideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These nonideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the nonideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern; the distortion of the electron diffusion size and shape; and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread functions and the sublinear variance in flats. The method is implemented in the Photon Simulator and the code is publicly available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1993

Relativistic effects on the bonding and properties of the hydrides of platinum

The ground state of PtH2 and several low-lying states of PtH(+) and PtH have been studied at the all-electron self-consistent-field level of theory to examine the importance of relativistic effects. The results of calculations based on Dirac-Hartree-Fock theory, nonrelativistic theory, and the spin-free no-pair relativistic approximation of Hess are compared to separate the effects of the spin-free terms and the spin-orbit terms of the Hamiltonian on the relativistic corrections to the molecular properties. Comparison is also made between first-order perturbation theory including the one-electron spin-free terms and the method of Hess to determine the size of effects beyond first order. It is found that the spin-orbit interaction significantly affects the properties and energetics of these molecules because of the participation of the Pt 5d orbitals in the bonding, and that effects beyond first order in perturbation theory are large. Any treatment of Pt compounds will have to include both the spin-free and spin-orbit interactions for an accurate description.

Dyall, Kenneth G.↗
DOE OSTI2022

The WRF-Solar Ensemble Prediction System: Development, Test, and Validation

Providing reliable probabilistic solar radiation information is needed to improve management of the uncertainty and variability of solar generation. Thus, guidance on how to develop skillful and accurate ensemble forecasts is essential and it will ultimately contribute to integration of high amounts of solar energy on the grid. A team from the National Renewable Energy Laboratory and the National Center for Atmospheric Research had been collaborating to develop the WRF-Solar ensemble prediction system (WRF-Solar EPS) in the past three years to produce probabilistic solar irradiance forecasts and better predict solar energy by quantifying forecast uncertainty. The WRF-Solar EPS basically generates ensemble members for solar irradiance based on stochastic perturbations to provide intraday and day-ahead probabilistic forecasts. This study will present main research steps in developing the WRF-Solar EPS including: (a) tangent linear analysis for identifying key input variables of six WRF-Solar modules significantly related to predicting of cloud and solar irradiance, (b) combining stochastic perturbation technique with the WRF-Solar model, and (c) ensemble calibration method to decrease error and uncertainty of ensemble-based solar forecasts. The capability of WRF-Solar EPS is now updated to the most recent version of standard WRF model. This presentation will summarize comprehensive results from the evaluation of forecasts against the National Solar Radiation Data Base as well as ground-measured observations. Moreover, we will introduce the user's guide for WRF-Solar EPS (e.g., parameters to configure stochastic perturbations) and future extension of this research.

day-ahead forecast↗
NASA NTRS1986

A Petrov-Galerkin finite element method for the compressible Euler and Navier-Stokes equations

An overview of a new finite element method for the compressible Euler and Navier-Stokes equations is presented. The discretization is based on entropy variables. The method is developed within the framework of a Petrov-Galerkin formulation. Two perturbations are added to the weighting function; one is a generalization of the SUPG operator and the other is designed to enhance shock capturing capability. The treatment of boundary conditions and the consistent calculation of boundary fluxes are addressed. Results of numerical tests are presented which confirm the robustness and wide applicability of the method.

Franca, L. P.↗
NASA NTRS1992

IOPS advisor: Research in progress on knowledge-intensive methods for irregular operations airline scheduling

Our research focuses on the problem of recovering from perturbations in large-scale schedules, specifically on the ability of a human-machine partnership to dynamically modify an airline schedule in response to unanticipated disruptions. This task is characterized by massive interdependencies and a large space of possible actions. Our approach is to apply the following: qualitative, knowledge-intensive techniques relying on a memory of stereotypical failures and appropriate recoveries; and quantitative techniques drawn from the Operations Research community's work on scheduling. Our main scientific challenge is to represent schedules, failures, and repairs so as to make both sets of techniques applicable to the same data. This paper outlines ongoing research in which we are cooperating with United Airlines to develop our understanding of the scientific issues underlying the practicalities of dynamic, real-time schedule repair.

Borse, John E.↗
NASA NTRS2019

A Time-Distance Helioseismology Method for Quasi-Linear Geometries

Helioseismology is the study of the solar interior, through which we extract flow and wave-speed information from Doppler velocity observations at the surface. Local helioseismology involves the study of small regions on the solar disk and is used to create a detailed picture of the interior in that particular region. Perturbations in the flow and wave-speed results indicate, e.g. magnetic-flux or temperature variations. There are multiple methods used in local-helioseismic research, but all current local-helioseismic techniques assume a point-source perturbation. For this study, we develop a new time-distance (TD) helioseismic methodology that can exploit the quasi-linear geometry of an elongated feature, allowing us to i) improve the signal-to-noise ratio of the TD results, and ii) greatly decrease the number of calculations required and therefore the computing time of the TD analysis. Ultimately, the new method will allow us to investigate solar features with magnetic-field configurations previously unexplored. We validate our new technique using a simple F-mode wave simulation, comparing results of point-source and linear perturbations. Results indicate that local-helioseismic analysis is dependent on the geometry of the system and can be improved by taking the magnetic-field configuration into account.

Hess Webber, Shea A.↗
NASA NTRS1971

Application of the steepest ascent optimization method to a reentry trajectory problem

The direct optimization method is presented in detail. Nominal values of the control variables are input parameters. Perturbations are introduced into the control variables and the resulting first order predictions of changes in the payoff, and constraint functions are then determined. Through a sequence of prescribed cycles, a trajectory is eventually obtained which is reasonably close to the optimum. The method is successfully applied to an Apollo three-dimensional reentry problem. The study of this Apollo application problem has resulted in the development of a highly flexible computer program that can be modified to consider other trajectory optimization problems.

Junkin, B. G.↗
NASA NTRS1965

Lower Bounds for the Energy Levels of Anharmonic Oscillators

The lower‐bounds method presented by Löwdin in the preceding paper has been applied to oscillators perturbed by third‐ and fourth‐power terms in the potential‐energy expression. For favorable cases agreement between upper and lower bounds is easily carried to many more figures than are likely to be physically significant. In many cases the lower bounds agreed more closely to the true eigenvalue than did the corresponding upper bounds. For a given basis set, this method gives closer bounds than that of Bazley and Fox, except for energy levels too high to be satisfactorily treated in the given basis. The only disadvantage found was that for close bounds double precision proved necessary, indicating more than ordinary loss of computational accuracy.

Reid, Charles E.↗
NASA NTRS1965

Studies in Perturbation Theory. XI. Lower Bounds to Energy Eigenvalues, Ground State, and Excited States

The bracketing theorem in the partitioning technique for solving the Schrödinger equation may be used in principle to determine upper and lower bounds to energy eigenvalues. Practical lower bounds of any accuracy desired may be evaluated by utilizing the properties of ``inner projections'' on finite manifolds in the Hilbert space. The method is here applied to the ground state and excited states of a Hamiltonian H=H(sub 0)+V having a positive definite perturbation V. Even if inspiration is derived from the method of intermediate Hamiltonians, the final results are of bracketing type and independent of this approach. The method is numerically illustrated in some accompanying papers.

Loewdin, Per-Olov↗
DOE OSTI2022

Analysis of SCALE Criticality and Sensitivity Calculations for Reflected HEU Cylinders

The SCALE code package offers multiple nuclear data libraries and sensitivity and uncertainty (S/U) methods supporting and derived from Monte Carlo (MC) transport. The CSAS and TSUNAMI-3D sequences use KENO MC, utilizing either continuous-energy (CE) cross sections or multigroup (MG) cross section libraries. TSUNAMI-3D has two CE calculational methods: the iterated fission probability (IFP) method, and the Contribution-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method. Previous work has shown poor agreement between CLUTCH and confirmatory direct perturbation calculations in specific applications (e.g., fissionable and polyethylene reflectors). The HEU-MET-FAST-084 (HMF-084) International Criticality Safety Benchmark Evaluation Project evaluation consists of 27 cylindrical highly enriched uranium metal cores with 14 unique reflector materials of 0.5 and 1 in. thicknesses. Included in this list of reflector materials are natural uranium and polyethylene. This work utilized SCALE 6.2.4 models of the HMF-084 evaluation, with additional non-physical configurations to test both the MG bias and CLUTCH functionality across a variety of reflector material thicknesses. The evaluation’s use of concentric cylinders allowed for examination of several MG self-shielding methods: infinite homogenous, cylindrical, and spherical. The results indicate that the use of polyethylene reflectors with CLUTCH is not fundamentally impossible but sensitive to geometry. The poor performance of CLUTCH with fissionable reflectors was reaffirmed. The 2 in. and greater polyethylene-reflected calculations demonstrate the necessity of using the 302-group library for fast systems. The nickel MG bias was substantial, as discussed in a companion paper, as were cobalt and iron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Subspace Methods in Multi-Parameter Seismic Full Waveform Inversion

In full waveform inversion (FWI) high-resolution subsurface model parameters are sought. FWI is normally treated as a nonlinear least-squares inverse problem, in which the minimum of the corresponding misfit function is found by updating the model parameters. When multiple elastic or acoustic properties are solved for, simple gradient methods tend to confuse parameter classes. This is referred to as parameter cross-talk; it leads to incorrect model solutions, poor convergence and strong dependence on the scaling of the different parameter types. Determining step lengths in a subspace domain, rather than directly in terms of gradients of different parameters, is a potentially valuable approach to address this problem. The particular subspace used can be defined over a span of different sets of data or different parameter classes, provided it involves a small number of vectors compared to those contained in the whole model space. Additionally, in a subspace method, the basis vectors are defined first, and a local minimum is found in the space spanned by these. We examine the application of the subspace method within acoustic FWI in determining simultaneously updates for velocity and density. We first discuss the choice of basis vectors to construct the spanned space, from linear updates by distinguishing only the contributions of different parameter classes towards nonlinear updates by adding the contributions of higher-order perturbations of each parameter class. The numerical character of FWI solutions generated via subspace methods involving different basis vectors is then analyzed and compared with traditional FWI methods. The subspace methods can provide better reconstructions of the model, especially for the velocity, as well as improved convergence rates, while the computational costs are still comparable with the traditional FWI methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗