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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

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 OSTI2015

The magnetic particle in a box: Analytic and micromagnetic analysis of probe-localized spin wave modes

The dipole field from a probe magnet can be used to localize a discrete spectrum of standing spin wave modes in a continuous ferromagnetic thin film without lithographic modification to the film. Obtaining the resonance field for a localized mode is not trivial due to the effect of the confined and inhomogeneous magnetization precession. We compare the results of micromagnetic and analytic methods to find the resonance field of localized modes in a ferromagnetic thin film, and investigate the accuracy of these methods by comparing with a numerical minimization technique that assumes Bessel function modes with pinned boundary conditions. We find that the micromagnetic technique, while computationally more intensive, reveals that the true magnetization profiles of localized modes are similar to Bessel functions with gradually decaying dynamic magnetization at the mode edges. We also find that an analytic solution, which is simple to implement and computationally much faster than other methods, accurately describes the resonance field of localized modes when exchange fields are negligible, and demonstrating the accessibility of localized mode analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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.↗

NASA NTRS1980

Fan noise-mode structure in a duct

Two computer programs help analyst meet low-noise limits on turbofan engines. Microphone Location Program computes optimum locations in turbofan duct for placement of microphones. After tests in first program are run, acoustic phase, amplitude, and pressure are used as inputs in Modal Calculation Program.

Love, R. A.↗

NASA NTRS1989

Recent developments in deployment analysis simulation using a multi-body computer code

Deployment is a candidate mode for construction of structural space systems components. By its very nature, deployment is a dynamic event, often involving large angle unfolding of flexible beam members. Validation of proposed designs and conceptual deployment mechanisms is enhanced through analysis. Analysis may be used to determine member loads thus helping to establish deployment rates and deployment control requirements for a given concept. Futhermore, member flexibility, joint free-play, manufacturing tolerances, and imperfections can affect the reliability of deployment. Analyses which include these effects can aid in reducing risks associated with a particular concept. Ground tests which can play a similar role to that of analyses are difficult and expensive to perform. Suspension systems just for vibration ground tests of large space structures in a 1 g environment present many challenges. Suspension of a structure which spatially expands is even more challenging. Analysis validation through experimental confirmation on relatively small simple models would permit analytical extrapolation to larger more complex space structures.

Housner, Jerrold M.↗

DOE OSTI2022

AI and extreme scale computing to learn and infer the physics of higher order gravitational wave modes of quasi-circular, spinning, non-precessing black hole mergers

We use artificial intelligence (AI) to learn and infer the physics of higher order gravitational wave modes of quasi-circular, spinning, non precessing binary black hole mergers. We trained AI models using 14 million waveforms, produced with the surrogate model NRHybSur3dq8, that include modes up to $\ell$ ≤ 4 and (5,5), except for (4,0) and (4,1), that describe binaries with mass-ratios $\textit{q}$ ≤ 8, individual spins $s^z_{\{1,2\}} \in$[–0.8,0.8], and inclination angle $θ \in$ [0,π]. Our probabilistic AI surrogates can accurately constrain the mass-ratio, individual spins, effective spin, and inclination angle of numerical relativity waveforms that describe such signal manifold. We compared the predictions of our AI models with Gaussian process regression, random forest, k-nearest neighbors, and linear regression, and with traditional Bayesian inference methods through the PyCBC Inference toolkit, finding that AI outperforms all these approaches in terms of accuracy, and are between three to four orders of magnitude faster than traditional Bayesian inference methods. Our AI surrogates were trained within 3.4 hours using distributed training on 1,536 NVIDIA V100 GPUs in the Summit supercomputer.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1996

The Integrated Mode Management Interface

Mode management is the processes of understanding the character and consequences of autoflight modes, planning and selecting the engagement, disengagement and transitions between modes, and anticipating automatic mode transitions made by the autoflight system itself. The state of the art is represented by the latest designs produced by each of the major airframe manufacturers, the Boeing 747-400, the Boeing 777, the McDonnell Douglas MD-11, and the Airbus A320/A340 family of airplanes. In these airplanes autoflight modes are selected by manipulating switches on the control panel. The state of the autoflight system is displayed on the flight mode annunciators. The integrated mode management interface (IMMI) is a graphical interface to autoflight mode management systems for aircraft equipped with flight management computer systems (FMCS). The interface consists of a vertical mode manager and a lateral mode manager. Autoflight modes are depicted by icons on a graphical display. Mode selection is accomplished by touching (or mousing) the appropriate icon. The IMMI provides flight crews with an integrated interface to autoflight systems for aircraft equipped with flight management computer systems (FMCS). The current version is modeled on the Boeing glass-cockpit airplanes (747-400, 757/767). It runs on the SGI Indigo workstation. A working prototype of this graphics-based crew interface to the autoflight mode management tasks of glass cockpit airplanes has been installed in the Advanced Concepts Flight Simulator of the CSSRF of NASA Ames Research Center. This IMMI replaces the devices in FMCS equipped airplanes currently known as mode control panel (Boeing), flight guidance control panel (McDonnell Douglas), and flight control unit (Airbus). It also augments the functions of the flight mode annunciators. All glass cockpit airplanes are sufficiently similar that the IMMI could be tailored to the mode management system of any modern cockpit. The IMMI does not replace the functions of the FMCS control and display unit. The purpose of the INMI is to provide flight crews with a shared medium in which they can assess the state of the autoflight system, take control actions on it, reason about its behavior, and communicate with each other about its behavior. The design is intended to increase mode awareness and provide a better interface to autoflight mode management. This report describes the IMMI, the methods that were used in designing and developing it, and the theory underlying the design and development processes.

Hutchins, Edwin↗

DOE OSTI2014

Integrated modeling of temperature profiles in L-mode tokamak discharges

Simulations of doublet III-D, the joint European tokamak, and the tokamak fusion test reactor L-mode tokamak plasmas are carried out using the PTRANSP predictive integrated modeling code. The simulation and experimental temperature profiles are compared. The time evolved temperature profiles are computed utilizing the Multi-Mode anomalous transport model version 7.1 (MMM7.1) which includes transport associated with drift-resistive-inertial ballooning modes (the DRIBM model [T. Rafiq et al., Phys. Plasmas 17, 082511 (2010)]). The tokamak discharges considered involved a broad range of conditions including scans over gyroradius, ITER like current ramp-up, with and without neon impurity injection, collisionality, and low and high plasma current. The comparison of simulation and experimental temperature profiles for the discharges considered is shown for the radial range from the magnetic axis to the last closed flux surface. The regions where various modes in the Multi-Mode model contribute to transport are illustrated. In the simulations carried out using the MMM7.1 model it is found that: The drift-resistive-inertial ballooning modes contribute to the anomalous transport primarily near the edge of the plasma; transport associated with the ion temperature gradient and trapped electron modes contribute in the core region but decrease in the region of the plasma boundary; and neoclassical ion thermal transport contributes mainly near the center of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

Expansion and orthogonalization of measured modes for structure identification

The purpose was to investigate a new simultaneous expansion/orthogonalization method in comparison with two previously published expansion methods and a widely used orthogonalization technique. Each expansion method uses data from an analytical model of the structure to complete the estimate of the mode shape vectors. Berman and Nagy used Guyan expansion in their work with improving analytical models. In this method, modes are expanded one at a time, producing a set not orthogonal with respect to the mass matrix. Baruch and Bar Itzhack's optimal orthogonalization procedure was used to subsequently adjust the expanded modes. A second expansion technique was presented by O'Callahan, Avitabile, and Reimer and separately by Kammer. Again, modes are expanded individually and orthogonalized after expansion with the same optimal technique as above. Finally, a simultaneous expansion/orthogonalization method was developed from the orthogonal Procrustes problem of computational mathematics. In this method modes are optimally expanded as a set and orthogonal with respect to the mass matrix as a result. Two demonstation problems were selected for the comparison of the methods described. The first problem is an 8 degree of freedom spring-mass problem first presented by Kabe. Several conditions were examined for expansion method including the presence of errors in the measured data and in the analysis models. As a second demonstration problem, data from tests of laboratory scale model truss structures was expanded for system identification. Tests with a complete structure produced a correlated analysis model and the stiffness and mass matrices. Tests of various damaged configurations produced measured data for 6 modes at 14 dof locations.

Smith, Suzanne Weaver↗

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 OSTI2025

Gate-Based Quantum Simulation of Gaussian Bosonic Circuits on Exponentially Many Modes

We introduce a framework for simulating, on an ( n + 1 )-qubit quantum computer, the action of a Gaussian bosonic (GB) circuit on a state over 2 n modes. Specifically, we encode the initial bosonic state’s expectation values over quadrature operators (and their covariance matrix) as an input qubit state. This is then evolved by a quantum circuit that effectively implements the symplectic propagators induced by the GB gates. We find families of GB circuits and initial states leading to efficient quantum simulations. For this purpose, we introduce a dictionary that maps between GB and qubit gates such that particle- (non-particle-) preserving GB gates lead to real- (imaginary-) time evolutions at the qubit level. For the special case of particle-preserving circuits, we present a bounded-error-quantum-polynomial time (BQP)-complete GB decision problem, indicating that GB evolutions of Gaussian states on exponentially many modes are as powerful as universal quantum computers. We also perform numerical simulations of an interferometer on ∼ 8 × 10 9 modes, illustrating the power of our framework. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1978

User's Manual for Thermal Analysis Program of Axially Grooved Heat Pipe (HTGAP)

A computer program that numerically predicts the steady state temperature distribution inside an axially grooved heat pipe wall for a given groove geometry and working fluid under various heat input and output modes is described. The program computes both evaporator and condenser film coefficients. The program is able to handle both axisymmetric and nonaxisymmetric heat transfer cases. Non-axisymmetric heat transfer results either from non-uniform input at the evaporator or non-uniform heat removal from the condenser, or from both. The presence of a liquid pool in the condenser region under one-g condition also causes non-axisymmetric heat transfer, and its effect on the pipe wall temperature distribution is included in the present program. The hydrodynamic aspect of an axially grooved heat pipe is studied in the Groove Analysis Program (GAP). The present thermal analysis program assumes that the GAP program (or other similar programs) is run first so that the heat transport limit and optimum fluid charge of the heat pipe are known a priori.

Kamotani, Y.↗

NASA NTRS1987

Tidal estimation in the Atlantic and Indian Oceans, 3 deg x 3 deg solution

An estimation technique was developed to extrapolate tidal amplitudes and phases over entire ocean basins using existing gauge data and the altimetric measurements provided by satellite oceanography. The technique was previously tested. Some results obtained by using a 3 deg by 3 deg grid are presented. The functions used in the interpolation are the eigenfunctions of the velocity (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 can provide the simulated data for the testing of the method and serve as a guide in choosing the most energetic functions for the interpolation.

Sanchez, Braulio V.↗