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 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Root damage analysis of aircraft engine blade subject to ice impact

The blade root response due to ice impact on an engine blade is simulated using the NASA in-house code BLASIM. The ice piece is modeled as an equivalent spherical object impacting on the leading edge of the blade and has the velocity opposite to that of the aircraft with direction parallel to the engine axis. The effect of ice impact is considered to be an impulse load on the blade with its amplitude computed based on the momentum transfer principle. The blade response due to the impact is carried out by modal superposition using the first three modes. The maximum dynamic stresses at the blade root are computed at the quarter cycle of the first natural frequency. A combined stress failure function based on modified distortion energy is used to study the spanwise bending damage response at the blade root. That damage function reaches maximum value for very low ice speeds and increases steeply with increases in engine speed.

Reddy, E. S.↗

NASA NTRS1987

Solar Dynamics and G-modes

Various properties of r-modes that would occupy the solar convection zone were computed. Values were calculated for the energy loss from the coefficient of eddy viscosity using the solar convective zone model. Viscous damping is strong and removes one-third or more of the mode energy each oscillation period. A fast Fourier transform (FFT) of daily Zurich sunspot number was made for eight sets of data, taken near peaks of solar cycles 14 to 21. The FFTs were evaluated with the Wolff and Hickey model. The horizontal component of F sub h of the force exerted by a planetary torque was calculated to see is it works on the photospheric layers at a rate comparable to 0.000001 solar luminosity. The time span of torque in the external sector of the Sun was calculated. The duration of the effect of two planets in conjunction on sunspot numbers was examined. The north-south asymmetry of solar activity was investigated.

Blizard, Jane B.↗

DOE OSTI2020

Computational micromechanics model based failure criteria for chopped carbon fiber sheet molding compound composites

Chopped carbon fiber sheet molding compound has a great potential in lightweight automotive, marine, and aerospace applications. One of the most challenging tasks is to predict the failure strength of the material due to its anisotropy and heterogeneity, as well as complex stress states in real-world working conditions. In this work, a novel constitutive model of carbon fiber chip is proposed to capture the pre- and post-failure behaviors under different loading modes. On this basis, we propose a new computational micromechanics model, which is calibrated and validated by uniaxial tensile, compressive, and in-plane shear experiments. Furthermore, a set of microstructures representative volume element (RVE) models under complex loading conditions are reconstructed to understand the relationship between the microstructure characteristics and the failure envelopes. Finally, several modified versions of classical failure criteria are proposed for anisotropic materials with consideration of the fiber orientation tensor. The modified Tsai-Wu failure criterion, which shows the best accuracy among all failure criteria, is highlighted in the comparative study.

36 MATERIALS SCIENCE↗

NASA NTRS1989

The inviscid axisymmetric stability of the supersonic flow along a circular cylinder

The supersonic flow past a thin straight circular cylinder is investigated. The associated boundary layer flow (i.e., the velocity and temperature field) is computed; the asymptotic, far downstream solution is obtained, and compared with the full numerical results. The inviscid, linear, axisymmetric (temporal) stability of this boundary layer is also studied. A so called doubly generalized inflexion condition is derived, which is a condition for the existence of so called subsonic neutral modes. The eigenvalue problem (for the complex wavespeed) is computed for two freestream Mach numbers (2.8 and 3.8), and this reveals that curvature has a profound effect on the stability of the flow. The first unstable inviscid mode is seen to rapidly disappear as curvature is introduced, while the second (and generally the most important) mode suffers a substantially reduced amplification rate.

Duck, Peter W.↗

DOE OSTI2024

Dynamics of qudit gates and effects of spectator modes on optimal control pulses

Qudit gates for high-dimensional quantum computing can be synthesized with high precision using numerical quantum optimal control techniques. Large circuits are broken down into modules and the tailored pulses for each module can be used as primitives for a qudit compiler. Application of the pulses of each module in the presence of extra modes may decrease their effectiveness due to crosstalk. In this paper, we address this problem by simulating qudit dynamics for circuit quantum electrodynamics systems. As a test case, we take pulses for single-qudit swap gates optimized in isolation and then apply them in the presence of spectator modes each of which are in Fock states. We provide an experimentally relevant scaling formula that can be used as a bound on the fidelity decay. Our results show that frequency shift from spectator mode populations has to be ≲ 0.1 % of the qudit's nonlinearity in order for high-fidelity single-qudit gates to be useful in the presence of occupied spectator modes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2015

Quantum dynamics on a three-sheeted six-dimensional ab initio potential-energy surface of the phosphine cation: Simulation of the photoelectron spectrum and the ultrafast radiationless decay dynamics

A diabatic three-sheeted six-dimensional potential-energy surface has been constructed for the ground state and the lowest excited state of the PH{sub 3}{sup +} cation. Coupling terms of Jahn-Teller and pseudo-Jahn-Teller origin up to eighth order had to be included to describe the pronounced anharmonicity of the surface due to multiple conical intersections. The parameters of the diabatic Hamiltonian have been optimized by fitting the eigenvalues of the potential-energy matrix to ab initio data calculated at the CASSCF/MRCI level employing the correlation-consistent triple-ζ basis. The theoretical photoelectron spectrum of phosphine and the non-adiabatic nuclear dynamics of the phosphine cation have been computed by propagating nuclear wave packets with the multiconfiguration time-dependent Hartree method. The theoretical photoelectron bands obtained by Fourier transformation of the autocorrelation function agree well with the experimental results. It is shown that the ultrafast non-radiative decay dynamics of the first excited state of PH{sub 3}{sup +} is dominated by the exceptionally strong Jahn-Teller coupling of the asymmetric bending vibrational mode together with a hyperline of conical intersections with the electronic ground state induced by the umbrella mode. Time-dependent population probabilities have been computed for the three adiabatic electronic states. The non-adiabatic Jahn-Teller dynamics within the excited state takes place within ≈5 fs. Almost 80% of the excited-state population decay to the ground state within about 10 fs. The wave packets become highly complex and delocalized after 20 fs and no further significant transfer of electronic population seems to occur up to 100 fs propagation time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1995

Parallelized modelling and solution scheme for hierarchically scaled simulations

This two-part paper presents the results of a benchmarked analytical-numerical investigation into the operational characteristics of a unified parallel processing strategy for implicit fluid mechanics formulations. This hierarchical poly tree (HPT) strategy is based on multilevel substructural decomposition. The Tree morphology is chosen to minimize memory, communications and computational effort. The methodology is general enough to apply to existing finite difference (FD), finite element (FEM), finite volume (FV) or spectral element (SE) based computer programs without an extensive rewrite of code. In addition to finding large reductions in memory, communications, and computational effort associated with a parallel computing environment, substantial reductions are generated in the sequential mode of application. Such improvements grow with increasing problem size. Along with a theoretical development of general 2-D and 3-D HPT, several techniques for expanding the problem size that the current generation of computers are capable of solving, are presented and discussed. Among these techniques are several interpolative reduction methods. It was found that by combining several of these techniques that a relatively small interpolative reduction resulted in substantial performance gains. Several other unique features/benefits are discussed in this paper. Along with Part 1's theoretical development, Part 2 presents a numerical approach to the HPT along with four prototype CFD applications. These demonstrate the potential of the HPT strategy.

Padovan, Joe↗

NASA NTRS1986

Acceleration effects in solid propellant rocket motors

The performance variations due to acceleration loads imposed on spinning solid propellant rocket motors are investigated. The four potentially most significant modes of acceleration-induced phenomena are identified from a study of the literature and modeled. The four modes are a mechanical mode which deals with deformations of the propellant and case: a thermodynamic mode which covers acceleration-induced combustion phenomena; a stress mode which covers the stressed propellant's effect on burn rate; and a gas dynamic mode which deals with changes in gas flow in the chamber and through the nozzle. Simplified models of each mode are developed or taken from the literature and are added to an internal ballistics evaluation computer program. The resulting analysis is the first to include all of the modes. In order to do this an original analysis of the mechanical and stress modes was necessary. However, the analysis shows that the stress mode is not important for the circular perforated grains studied. The other effects are shown to have a significant influence on solid rocket motor performance. The magnitude of the different mode effects are such that one may not be ignored over the others as has been done in the past. The results of the analysis are compared to published rocket motor data. The comparisons indicate an erosive burning effect that is a function of spin rate. A qualitative explanation of the erosive effect is presented.

Langhenry, M. T.↗

DOE OSTI2014

Towards the clinical implementation of iterative low-dose cone-beam CT reconstruction in image-guided radiation therapy: Cone/ring artifact correction and multiple GPU implementation

Purpose: Compressed sensing (CS)-based iterative reconstruction (IR) techniques are able to reconstruct cone-beam CT (CBCT) images from undersampled noisy data, allowing for imaging dose reduction. However, there are a few practical concerns preventing the clinical implementation of these techniques. On the image quality side, data truncation along the superior–inferior direction under the cone-beam geometry produces severe cone artifacts in the reconstructed images. Ring artifacts are also seen in the half-fan scan mode. On the reconstruction efficiency side, the long computation time hinders clinical use in image-guided radiation therapy (IGRT). Methods: Image quality improvement methods are proposed to mitigate the cone and ring image artifacts in IR. The basic idea is to use weighting factors in the IR data fidelity term to improve projection data consistency with the reconstructed volume. In order to improve the computational efficiency, a multiple graphics processing units (GPUs)-based CS-IR system was developed. The parallelization scheme, detailed analyses of computation time at each step, their relationship with image resolution, and the acceleration factors were studied. The whole system was evaluated in various phantom and patient cases. Results: Ring artifacts can be mitigated by properly designing a weighting factor as a function of the spatial location on the detector. As for the cone artifact, without applying a correction method, it contaminated 13 out of 80 slices in a head-neck case (full-fan). Contamination was even more severe in a pelvis case under half-fan mode, where 36 out of 80 slices were affected, leading to poorer soft tissue delineation and reduced superior–inferior coverage. The proposed method effectively corrects those contaminated slices with mean intensity differences compared to FDK results decreasing from ∼497 and ∼293 HU to ∼39 and ∼27 HU for the full-fan and half-fan cases, respectively. In terms of efficiency boost, an overall 3.1 × speedup factor has been achieved with four GPU cards compared to a single GPU-based reconstruction. The total computation time is ∼30 s for typical clinical cases. Conclusions: The authors have developed a low-dose CBCT IR system for IGRT. By incorporating data consistency-based weighting factors in the IR model, cone/ring artifacts can be mitigated. A boost in computational efficiency is achieved by multi-GPU implementation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1990

The inviscid axisymmetric stability of the supersonic flow along a circular cylinder

The supersonic flow past a thin straight circular cylinder is investigated. The associated boundary-layer flow (i.e. the velocity and temperature field) is computed; the asymptotic, far downstream solution is obtained, and compared with the full numerical results. The inviscid, linear, axisymmetric (temporal) stability of this boundary layer is also studied. A so-called 'doubly generalized' inflexion condition is derived, which is a condition for the existence of so-called 'subsonic' neutral modes. The eigenvalue problem (for the complex wavespeed) is computed for two free-stream Mach numbers (2.8 and 3.8), and this reveals that curvature has a profound effect on the stability of the flow. The first unstable inviscid mode is seen to disappear rapidly as curvature is introduced, while the second (and generally the most important) mode suffers a substantially reduced amplification rate.

Duck, Peter W.↗

NASA NTRS1982

Development of a locally mass flux conservative computer code for calculating 3-D viscous flow in turbomachines

The VANS successive approximation numerical method was extended to the computation of three dimensional, viscous, transonic flows in turbomachines. A cross-sectional computer code, which conserves mass flux at each point of the cross-sectional surface of computation was developed. In the VANS numerical method, the cross-sectional computation follows a blade-to-blade calculation. Numerical calculations were made for an axial annular turbine cascade and a transonic, centrifugal impeller with splitter vanes. The subsonic turbine cascade computation was generated in blade-to-blade surface to evaluate the accuracy of the blade-to-blade mode of marching. Calculated blade pressures at the hub, mid, and tip radii of the cascade agreed with corresponding measurements. The transonic impeller computation was conducted to test the newly developed locally mass flux conservative cross-sectional computer code. Both blade-to-blade and cross sectional modes of calculation were implemented for this problem. A triplet point shock structure was computed in the inducer region of the impeller. In addition, time-averaged shroud static pressures generally agreed with measured shroud pressures. It is concluded that the blade-to-blade computation produces a useful engineering flow field in regions of subsonic relative flow; and cross-sectional computation, with a locally mass flux conservative continuity equation, is required to compute the shock waves in regions of supersonic relative flow.

Walitt, L.↗

DOE OSTI2015

Reactor Pressure Vessel Fracture Analysis Capabilities in Grizzly

Efforts have been underway to develop fracture mechanics capabilities in the Grizzly code to enable it to be used to perform deterministic fracture assessments of degraded reactor pressure vessels (RPVs). Development in prior years has resulted a capability to calculate -integrals. For this application, these are used to calculate stress intensity factors for cracks to be used in deterministic linear elastic fracture mechanics (LEFM) assessments of fracture in degraded RPVs. The -integral can only be used to evaluate stress intensity factors for axis-aligned flaws because it can only be used to obtain the stress intensity factor for pure Mode I loading. Off-axis flaws will be subjected to mixed-mode loading. For this reason, work has continued to expand the set of fracture mechanics capabilities to permit it to evaluate off-axis flaws. This report documents the following work to enhance Grizzly’s engineering fracture mechanics capabilities for RPVs: • Interaction Integral and -stress: To obtain mixed-mode stress intensity factors, a capability to evaluate interaction integrals for 2D or 3D flaws has been developed. A -stress evaluation capability has been developed to evaluate the constraint at crack tips in 2D or 3D. Initial verification testing of these capabilities is documented here. • Benchmarking for axis-aligned flaws: Grizzly’s capabilities to evaluate stress intensity factors for axis-aligned flaws have been benchmarked against calculations for the same conditions in FAVOR. • Off-axis flaw demonstration: The newly-developed interaction integral capabilities are demon- strated in an application to calculate the mixed-mode stress intensity factors for off-axis flaws. • Other code enhancements: Other enhancements to the thermomechanics capabilities that relate to the solution of the engineering RPV fracture problem are documented here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1986

Ground Software Maintenance Facility (GSMF) user's manual

Instructions for the Ground Software Maintenance Facility (GSMF) system user is provided to operate the GSMF in all modes. The GSMF provides the resources for the Automatic Test Equipment (ATE) computer program maintenance (GCOS and GOAL). Applicable reference documents are listed. An operational overview and descriptions of the modes in terms of operator interface, options, equipment, material utilization, and operational procedures are contained. Test restart procedures are described. The GSMF documentation tree is presented including the user manual.

Aquila, V.↗

NASA NTRS2011

Method of performing computational aeroelastic analyses

Computational aeroelastic analyses typically use a mathematical model for the structural modes of a flexible structure and a nonlinear aerodynamic model that can generate a plurality of unsteady aerodynamic responses based on the structural modes for conditions defining an aerodynamic condition of the flexible structure. In the present invention, a linear state-space model is generated using a single execution of the nonlinear aerodynamic model for all of the structural modes where a family of orthogonal functions is used as the inputs. Then, static and dynamic aeroelastic solutions are generated using computational interaction between the mathematical model and the linear state-space model for a plurality of periodic points in time.

Silva, Walter A.↗

NASA NTRS2025

PSE-Based Aerodynamic Flow Transition Prediction Using Automated Unstructured CFD Integration

Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.

Transition Prediction↗

NASA NTRS2020

Streak Instabilities on HIFiRE-5 Elliptic Cone

The Hypersonic International Flight Research and Experimentation (HIFiRE) 5 flight experiment by the Air Force Research Laboratories and Australian Defense Science and Technology Organization yielded inflight, hypersonic transition data for a canonical, 3D configuration in the form of an elliptic cone at various combinations of angle-of-attack and yaw. Previous analysis has shown that boundary-layer stability calculations for the as-flown, zero-yaw flight conditions are able to capture the second mode dominated transition along the attachment line as well as the crossflow induced transition in between the attachment line and the minor axis symmetry plane. However, there is still a need to understand the transition in the vicinity of the minor axis plane of symmetry, where the boundary layer thickness exhibits an azimuthally compact bulge due to the collision of the secondary flow approaching from the attachment lines on both sides. This paper reports a computational analysis of the stability of the boundary layer flow along the bulge regions on both the leeward and windward sides of the vehicle surface. At a representative flight condition corresponding to a time instant of t = 514.813 seconds, computations are performed to identify several different unstable modes within the bulge regions, with the eventual goal of developing a linear stability correlation for the transition onset near the minor axis symmetry plane. The most amplified flow instabilities in the pretransitional region correspond to a “stem” mode that is concentrated within a pair of nearly vertical shear layers below the cap of the mushroom shaped bulge in the mean velocity contours. An excellent agreement is found between the predictions of instability amplification based on the plane marching parabolized stability equations (PSE) and the direct numerical solutions (DNS). Given the large number of unstable modes within the bulge region, the observed agreement demonstrates, for first time, that the PSE can provide accurate predictions for the development of flow instabilities in streak-like flow features over three-dimensional bodies. The peak N-factor of approximately 15 at the transition location estimated from the thermocouple measurements along the leeward centerline is close to the transition N-factor for Mack mode transition during the HIFiRE-1 circular cone flight experiment with nearly zero degrees angle of attack.

Meelan Choudhari↗

DOE OSTI2012

A computational approach to continuum damping of Alfven waves in two and three-dimensional geometry

While the usual way of calculating continuum damping of global Alfven modes is the introduction of a small artificial resistivity, we present a computational approach to the problem based on a suitable path of integration in the complex plane. This approach is implemented by the Riccati shooting method and it is shown that it can be transferred to the Galerkin method used in three-dimensional ideal magneto-hydrodynamics (MHD) codes. The new approach turns out to be less expensive with respect to resolution and computation time than the usual one. We present an application to large aspect ratio tokamak and stellarator equilibria retaining a few Fourier harmonics only and calculate eigenfunctions and continuum damping rates. These may serve as an input for kinetic MHD hybrid models making it possible to bypass the problem of having singularities on the path of integration on one hand and considering continuum damping on the other.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Higher order topological superconductivity in magnet-superconductor hybrid systems

Abstract Quantum engineering of topological superconductors and of the ensuing Majorana zero modes might hold the key for realizing topological quantum computing and topology-based devices. Magnet-superconductor hybrid (MSH) systems have proven to be experimentally versatile platforms for the creation of topological superconductivity by custom-designing the complex structure of their magnetic layer. Here, we demonstrate that higher order topological superconductivity (HOTSC) can be realized in two-dimensional MSH systems by using stacked magnetic structures. We show that the sensitivity of the HOTSC to the particular magnetic stacking opens an intriguing ability to tune the system between trivial and topological phases using atomic manipulation techniques. We propose that the realization of HOTSC in MSH systems, and in particular the existence of the characteristic Majorana corner modes, allows for the implementation of a measurement-based protocols for topological quantum computing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗