Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational modes”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Hypersonic Second-Mode Instability Response to Shaped Roughness

An experimental campaign was conducted on a 7-degree half-angle cone in the NASA Langley Research Center 20-Inch Mach 6 Wind Tunnel to examine the influence of arrays of regularly spaced roughness elements on instability growth and transition. The primary element shape was a pair of elliptical planform ramps that were inclined at equal and opposite angles with respect to the local streamwise direction. The element shapes were designed to induce transient growth disturbances that would lead to sustained azimuthal modulation of the boundary layer flow while limiting the nearfield disturbances to avoid an immediate, i.e., effective tripping of the boundary layer. The bulk of the run matrix consisted of testing different element height sat free stream unit Reynolds numbers ranging from 9.8 to 13.1 million per meter. Other element shapes previously designed for tripping hypersonic boundary layers were also implemented. The model was instrumented with surface mounted Kulite ® and PCB ® pressure transducers and thermocouples. Spectra from the PCBs ® indicated clear suppression of the second-mode instability; however, neither the PCB ® spectra nor the heat transfer data presented strong evidence for delayed turbulent flow. Complementary stability computations likewise demonstrated second-mode reduction, particularly just downstream of the roughness, but also revealed a rise in first mode (streak-instability) amplitudes from the baseline that was likely responsible for the earlier transition observed for taller roughness cases.

Hypersonics↗
NASA NTRS2001

Ab Initio and Improved Empirical Potentials for the Calculation of the Anharmonic Vibrational States and Intramolecular Mode Coupling of N-Methylacetamide

The second-order Moller-Plesset ab initio electronic structure method is used to compute points for the anharmonic mode-coupled potential energy surface of N-methylacetamide (NMA) in the trans(sub ct) configuration, including all degrees of freedom. The vibrational states and the spectroscopy are directly computed from this potential surface using the Correlation Corrected Vibrational Self-Consistent Field (CC-VSCF) method. The results are compared with CC-VSCF calculations using both the standard and improved empirical Amber-like force fields and available low temperature experimental matrix data. Analysis of our calculated spectroscopic results show that: (1) The excellent agreement between the ab initio CC-VSCF calculated frequencies and the experimental data suggest that the computed anharmonic potentials for N-methylacetamide are of a very high quality; (2) For most transitions, the vibrational frequencies obtained from the ab initio CC-VSCF method are superior to those obtained using the empirical CC-VSCF methods, when compared with experimental data. However, the improved empirical force field yields better agreement with the experimental frequencies as compared with a standard AMBER-type force field; (3) The empirical force field in particular overestimates anharmonic couplings for the amide-2 mode, the methyl asymmetric bending modes, the out-of-plane methyl bending modes, and the methyl distortions; (4) Disagreement between the ab initio and empirical anharmonic couplings is greater than the disagreement between the frequencies, and thus the anharmonic part of the empirical potential seems to be less accurate than the harmonic contribution;and (5) Both the empirical and ab initio CC-VSCF calculations predict a negligible anharmonic coupling between the amide-1 and other internal modes. The implication of this is that the intramolecular energy flow between the amide-1 and the other internal modes may be smaller than anticipated. These results may have important implications for the anharmonic force fields of peptides, for which N-methylacetamide is a model.

Gregurick, Susan K.↗
DOE OSTI2022

High-order multirate explicit time-stepping schemes for the baroclinic-barotropic split dynamics in primitive equations

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta approach, we propose two high-order multirate explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The proposed schemes allow for a large time step to be used for the three-dimensional baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast) mode, in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the baroclinic velocity is first computed by advancing the baroclinic mode and fluid thickness of the system with the large time-step and the assistance of some intermediate approximations of the barotropicmode obtained by substepping with the small time step; then the barotropic velocity is corrected by using the small time step to re-advance the barotropic mode un-der an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting are relieved via a reconciliation process with carefully calculated flux deficits. Here, two benchmark tests from the “MPAS-Ocean” platform are carried out to numerically demonstrate the performance and parallel scalability of the proposed SSPRK-SE schemes.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Effects of Nose Bluntness on Hypersonic Boundary-Layer Receptivity and Stability Over Cones

The receptivity to freestream acoustic disturbances and the stability properties of hypersonic boundary layers are numerically investigated for boundary-layer flows over a 5 straight cone at a freestream Mach number of 6.0. To compute the shock and the interaction of the shock with the instability waves, the Navier-Stokes equations in axisymmetric coordinates were solved. In the governing equations, inviscid and viscous flux vectors are discretized using a fifth-order accurate weighted-essentially-non-oscillatory scheme. A third-order accurate total-variation-diminishing Runge-Kutta scheme is employed for time integration. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. The appearance of instability waves near the nose region and the receptivity of the boundary layer with respect to slow mode acoustic waves are investigated. Computations confirm the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary-layer transition. The current solutions, compared with experimental observations and other computational results, exhibit good agreement.

Kara, Kursat↗
NASA NTRS1987

The configuration of slow-mode shocks

A simple method for computing blunt-body flows containing slow shocks is presented. The method is based on the solution of a free-boundary problem for the shock configuration that matches a postshock potential flow to a uniform preshock flow. The method is applied to shocks arising from flow about spherical and paraboloidal obstacles. The results of a parametric study of shock formation as a function of sonic and Alfven Mach numbers showed that when either Mach number decreased, the standoff distance increased. The results should give estimates of the location and the configuration of slow shocks, provided that two as yet unverified assumptions are correct: (1) that the standoff distance and the slow-shock configuration along the stagnation streamline are not affected significantly by portions of the shock still forming far upstream of the shock nose, and (2) that fast-mode waves propagating at small angles to the magnetic field do not significantly affect the formation of the slow shock.

Wolfson, Richard↗
NASA NTRS1990

An experimental modal testing/identification technique for personal computers

A PC-based system for mode shape identification is evaluated. A time-domain modal identification procedure is utilized to identify the mode shapes of a beam apparatus from discrete time-domain measurements. The apparatus includes a cantilevered aluminum beam, four accelerometers, four low-pass filters, and the computer. The method's algorithm is comprised of an identification algorithm: the Eigensystem Realization Algorithm (ERA) and an estimation algorithm called Minimum Model Error (MME). The identification ability of this algorithm is compared with ERA alone, a frequency-response-function technique, and an Euler-Bernoulli beam model. Detection of modal parameters and mode shapes by the PC-based time-domain system is shown to be accurate in an application with an aluminum beam, while mode shapes identified by the frequency-domain technique are not as accurate as predicted. The new method is shown to be significantly less sensitive to noise and poorly excited modes than other leading methods. The results support the use of time-domain identification systems for mode shape prediction.

Roemer, Michael J.↗
NASA NTRS2003

Inflow/Outflow Conditions for Unsteady Aerodynamics and Aeroacoustics in Nonuniform Flow

The effect of a nonuniform mean flow on the normal modes; the inflow/outflow nonreflecting boundary conditions; and the sound power are studied. The normal modes in an annular duct are computed using a spectral method in combination with a shooting method. The swirl causes force imbalance which couples the acoustic and vortical modes. The acoustic modes are distinguished from the vortical modes by their large pressure and small vorticity content. The mean swirl also produces a Doppler shift in frequency. This results in more counter-spinning modes cut-on at a given frequency than modes spinning with the swirl. Nonreflecting boundary conditions are formulated using the normal mode solutions. The inflow/outflow boundary conditions are implemented in a linearized Euler scheme and validated by computing the propagation of acoustic and vortical waves in a duct for a variety of swirling mean flows. Numerical results show that the evolution of the vortical disturbances is sensitive to the inflow conditions and the details of the wake excitations. All three components of the wake velocity must be considered to correctly compute the wake evolution and the blade upwash. For high frequencies, the acoustic-vortical mode coupling is weak and a conservation equation for the acoustic energy can be derived. Sound power calculations show significant mean flow swirl effects, but mode interference effects are small.

Atassi, Oliver V.↗
NASA NTRS1987

On computing eigensolution sensitivity data using free vibration solutions

A simplified method of computing eigensolution sensitivity derivatives in structural dynamics is developed. It is shown that if the elements of stiffness and mass matrices associated with a design variable are homogeneous functions of that design variable, then eigenvalue derivatives can be computed from element strain and kinetic energies. Furthermore, if cross-mode energies are known, eigensolution derivatives of modified systems can be computed approximately using assume mode reanalysis formulation. A ten bar truss example is used to illustrate the present formulations.

Wang, B. P.↗
NASA NTRS1976

The Application of Recent Techniques in Flight Flutter Testing

The relative merits of sinusoidal excitation versus random atmospheric turbulence was investigated. The randomdec and autocorrelation methods were used to analyze data from a Learjet flight test. A parameter identification digital program, using least squares approach, was developed to determine the aeroelastic characteristics of a two mode system. The flight test program, computer program, and data reduction procedure is presented. Final results of the two modes of excitation obtained by Randomdec method are discussed.

Abla, M. A.↗
NASA NTRS1995

PSC implementation and integration

Hardware and software design of the performance seeking control (PSC) for the NASA F-15 research aircraft are described. The hardware architecture, vehicle management system computer (VMSC), pilot interface, and PSC mode selection are discussed. The PSC software is distributed among the VMSC, central computer, digital electronic engine controls (DEEC's), and electronic air inlet controllers (EAIC's). The major PSC modules, VMSC logic, VMSC channel C memory requirements, VMSC channel C timing, and navigation control indicator (NCI) variables and where they are located are presented.

Nobbs, Steven↗
NASA NTRS1975

The effects of localized damping on structural response

The effect of localized structural damping on the excitability of higher order normal modes of the large space telescope was investigated. A preprocessor computer program was developed to incorporate Voigt structural joint damping models in a NASTRAN finite-element dynamic model. A postprocessor computer program was developed to select critical modes for low-frequency attitude control problems and for higher frequency fine-stabilization problems. The mode selection is accomplished by ranking the flexible modes based on coefficients for rate gyro, position gyro, and optical sensors, and on image-plane motions due to sinusoidal or random power spectral density force and torque inputs.

Merchant, D. H.↗
DOE OSTI2023

Pitfalls in the n -mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)—effectively a many-body expansion in the space of molecular vibrations—is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by “painting in” the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes—and restore SCF convergence—in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.

Chemistry↗
DOE OSTI2022

Insight into Pyroelectricity and Phase Transitions in Ferroelectrics from Nonequilibrium Approach: The Case of PbTiO 3

Computational methodology is proposed to simulate time dependent temperature variations via thermostat approach. Here, the methodology is applied to study pyroelectricity under nonequilibrium conditions in prototypical ferroelectric PbTiO 3 . It is found that at room temperature the nonequilibrium effects in pyroelectricity begin manifesting above 1.0 THz and originate from the overlap with the intrinsic soft mode frequency. Therefore, the proposed computational methodology can be used to predict equilibrium pyroelectric response by choosing a frequency of temperature variation below the intrinsic polar mode frequencies for the material. Probing the dynamics of ferroelectric phase transition with the proposed methodology revealed that it takes about a nanosecond for the computational supercell to undergo phase transition. It is concluded that computational techniques capable of simulating picoseconds can be used to compute pyroelectric response from nonequilibrium dynamical approach while those that can reach nanoseconds are suitable for simulating phase transitions in ferroelectrics.

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

Effects of symmetric and asymmetric modes on transonic aeroelastic characteristics of full-span wing-body configurations

To accurately study the transonic aeroelastic characteristics, it is important to model the full aircraft configuration, including asymmetry. Recently, an accurate method of computing unsteady transonic flows on full-span wing-body configurations was developed using the transonic potential flow theory. In this work, the method is further developed to account for the aeroelasticity of full-span wing-body configurations. This is accomplished by simultaneously integrating the unsteady aerodynamic forces and modal structural equations of the wing-body configurations. To validate the method, aeroelastic computations are made for a wing-body configuration with a rectangular wing. The aeroelastic responses of this configuration are correlated with the responses of a similar isolated wing. The comparisons are favorable. Aeroelastic computations associated with symmetric and asymmetric modes are also made to study the influence of modal asymmetry on responses. This new development is further illustrated by computing aeroelastic characteristics of a typical fighter aircraft. The results from this study will be useful in accurately computing the transonic flutter boundaries of aircraft, including those associated with asymmetric modes.

Guruswamy, Guru P.↗
DOE OSTI2022

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗
DOE OSTI2025

Light in the dark forest. Part I. An efficient optimal estimator for 3D Lyman-alpha forest power spectrum

The highly anisotropic nature of the Lyman-alpha (Lyα) forest data introduces a complex survey window function that complicates the measurement of the three-dimensional power spectrum ( P 3D ). In this paper, we present the first fully optimal estimator for P 3D , which exactly deconvolves the survey window function and marginalizes contaminated modes that distort the power spectrum. Our approach adapts optimal estimator techniques developed for the 2D cosmic microwave background data to the 3D case. To achieve computational feasibility, we employ the conjugate gradient method and implement the P 3 M formalism to handle large-scale and small-scale operations separately and efficiently. We validate our estimator using Monte Carlo mocks and Gaussian simulations, demonstrating its accuracy and computational efficiency. We confirm that mode marginalization eliminates distortions arising from quasar continuum errors and delivers robust power spectrum estimation, though it also inflates errors at large scales. This first implementation works in the flat-sky case; we discuss the remaining steps needed to generalize it to the curved-sky case. This formalism offers a foundation for the Lyα forest P 3D measurements and a new path toward cosmological constraints from the Lyα forest data.

Lyman alpha forest↗
NASA NTRS1986

Objective analysis of tidal fields in the Atlantic and Indian Oceans

An objective analysis technique has been developed to extrapolate tidal amplitudes and phases over entire ocean basins using existing gauge data and the altimetric measurements which are now beginning to be provided by satellite oceanography. The technique was previously tested in the Lake Superior basin. The method has now been developed and applied in the Atlantic-Indian ocean basins using a 6 deg x 6 deg grid to test its essential features. The functions used in the interpolation are the eigenfunctions of the velocity potential (Proudman functions) which are computed numerically from a knowledge of the basin's bottom topography, the horizontal plan form and the necessary boundary conditions. These functions are characteristic of the particular basin. The gravitational normal modes of the basin are computed as part of the investigation, they are used to obtain the theoretical forced solutions for the tidal constituents, the latter provide the simulated data for the testing of the method and serve as a guide in choosing the most energetic modes for the objective analysis. The results of the objective analysis of the M2 and K1 tidal constituents indicate the possibility of recovering the tidal signal with a degree of accuracy well within the error bounds of present day satellite techniques.

Sanchez, B. V.↗
DOE OSTI2021

Assessment of the Amplification Factor Transport Transition Model for High-Mach Number Flows

In its present state, the Amplification Factor Transport (AFT) model provides an efficient alternative for high-fidelity stability methods for complex geometries subject to low-speed conditions considered in the industrial-production CFD environment. Its performance has not yet been documented for the typical transition scenarios in hypersonic conditions. The goal of the present study is to apply the AFT model to a yawed- and flared-cone configuration, for which the laminar flows respectively support the crossflow and second-mode instability mechanisms, and to compare the amplification results against high-fidelity results computed with Linear Parabolized Stability Equations (LPSE). High-resolution viscous-flow solutions are provided with OVERFLOW and COFFE. For the yawed-cone case, the AFT model yields the largest amplification in a region of the flow where the expected crossflow instability is not amplified (in the leeward symmetry plane) according to LPSE. For the flared-cone case with a wall-temperature slightly lower than the adiabatic value, the AFT model yields extremely large amplification factors (in excess of 100), exceeding the computed LPSE amplification of second-mode instability results by a factor 14. Upon considering a cold-wall condition, the model shows much smaller amplification factors, which is the opposite with respect to the expected behavior of the second mode. The AFT model does not adequately capture the crossflow and second-mode instability mechanisms in its present state.

Groot, Koen J.↗