Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reduced order modeing”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Optimal Recursive Digital Filters for Active Bending Stabilization

In the design of flight control systems for large flexible boosters, it is common practice to utilize active feedback control of the first lateral structural bending mode so as to suppress transients and reduce gust loading. Typically, active stabilization or phase stabilization is achieved by carefully shaping the loop transfer function in the frequency domain via the use of compensating filters combined with the frequency response characteristics of the nozzle/actuator system. In this paper we present a new approach for parameterizing and determining optimal low-order recursive linear digital filters so as to satisfy phase shaping constraints for bending and sloshing dynamics while simultaneously maximizing attenuation in other frequency bands of interest, e.g. near higher frequency parasitic structural modes. By parameterizing the filter directly in the z-plane with certain restrictions, the search space of candidate filter designs that satisfy the constraints is restricted to stable, minimum phase recursive low-pass filters with well-conditioned coefficients. Combined with optimal output feedback blending from multiple rate gyros, the present approach enables rapid and robust parametrization of autopilot bending filters to attain flight control performance objectives. Numerical results are presented that illustrate the application of the present technique to the development of rate gyro filters for an exploration-class multi-engined space launch vehicle.

Orr, Jeb S.↗
DOE OSTI2026

Coherent anti-Stokes Raman scattering with squeezed light: CARS for quantum-enhanced spectroscopy and imaging

We theoretically investigate quantum-enhanced coherent anti-Stokes Raman scattering (CARS) using squeezed light to amplify vibrational transition rates at low photon flux. Quantum sensing approaches are needed for nondestructive nanometrology such as in bioimaging where reduced photodamage is desired while retaining resolution and sensitivity. We analyze both single-mode squeezing applied to the pump field and two-mode squeezing between the pump and Stokes fields. We also show that the ordering of displacement and squeezing operations—whether displacement precedes squeezing or squeezing precedes displacement—has an impact on the resulting CARS transition amplitudes due to a difference in the photon number and the quantum-enhancement coefficients, with the latter offering a stronger enhancement in the case of two modes squeezing of the pump and Stokes under experimentally accessible conditions. Furthermore, our calculations capture these quantum enhancements through the intrinsic photon-number correlations of squeezed light, eliminating the need for interferometric detection or higher pump powers that are otherwise required to reach comparable sensitivities in classical CARS. Finally, we outline a quantum plasmonic extension of our model in which local field enhancements caused by surface plasmon excitation in metallic nanoparticles can be incorporated via mode-selective field amplification factors, offering a pathway toward combining squeezed-light quantum optics with surface-enhanced nanoscale spectroscopy and imaging.

Atomic & molecular structure↗
NASA NTRS1973

Resonant oscillations of intermediate frequency in a stratified atmosphere.

A class of solutions to a model of forced oscillations in a rotating stratified atmospheric layer is derived and analyzed. The basic model is found to reduce to a boundary value problem with a second-order linear partial differential equation of the hyperbolic type for this range of forcing frequencies. The forced solutions are shown to exhibit resonances with the normal modes of oscillation of the layer. The characteristics of the resonant modes are analyzed in terms of mean tropospheric values of temperature, temperature lapse, wind speed, horizontal and vertical wind shears, latitude, and the frequency and horizontal wavelength of the forcing mechanism. These solutions are compared with solutions to the model for a different (subinertial) range of forcing frequencies. This comparison leads to an elliptic boundary value problem. The solutions in that case do not exhibit the same type of resonance and generally decay away from the region of forcing.

Thompson, O. E.↗
NASA NTRS2007

Optimization of Acoustic Pressure Measurements for Impedance Eduction

As noise constraints become increasingly stringent, there is continued emphasis on the development of improved acoustic liner concepts to reduce the amount of fan noise radiated to communities surrounding airports. As a result, multiple analytical prediction tools and experimental rigs have been developed by industry and academia to support liner evaluation. NASA Langley has also placed considerable effort in this area over the last three decades. More recently, a finite element code (Q3D) based on a quasi-3D implementation of the convected Helmholtz equation has been combined with measured data acquired in the Langley Grazing Incidence Tube (GIT) to reduce liner impedance in the presence of grazing flow. A new Curved Duct Test Rig (CDTR) has also been developed to allow evaluation of liners in the presence of grazing flow and controlled, higher-order modes, with straight and curved waveguides. Upgraded versions of each of these two test rigs are expected to begin operation by early 2008. The Grazing Flow Impedance Tube (GFIT) will replace the GIT, and additional capabilities will be incorporated into the CDTR. The current investigation uses the Q3D finite element code to evaluate some of the key capabilities of these two test rigs. First, the Q3D code is used to evaluate the microphone distribution designed for the GFIT. Liners ranging in length from 51 to 610 mm are investigated to determine whether acceptable impedance eduction can be achieved with microphones placed on the wall opposite the liner. This analysis indicates the best results are achieved for liner lengths of at least 203 mm. Next, the effects of moving this GFIT microphone array to the wall adjacent to the liner are evaluated, and acceptable results are achieved if the microphones are placed off the centerline. Finally, the code is used to investigate potential microphone placements in the CDTR rigid wall adjacent to the wall containing an acoustic liner, to determine if sufficient fidelity can be achieved with 32 microphones available for this purpose. Initial results indicate 32 microphones can provide acceptable measurements to support impedance eduction with this test rig.

Jones, M. G.↗
DOE OSTI2020

Modular high temperature gas reactor core modeling with RELAP5-3D/PHISICS – Optimization schemes for load following

The objective of this study is to develop a 3D/NK model to study the load following operation of a Modular High-Temperature Gas-Cooled Reactor (MHTGR). The selected design is the 350 MW prismatic, graphite moderated, helium cooled thermal reactor based on the MHTGR-350 transient benchmark led by the High-Temperature Gas-Cooled Reactor (HTGR) Methods Core Simulation Group at Idaho National Laboratory (INL) in the framework of the Organization for Economic Cooperation and Development Nuclear Energy Agency (OECD/NEA). This kind of reactor usually reacts very slowly to the perturbation of the core parameters due to the large amount of graphite in the core. This kind of behavior promotes the reactor stability but on the other hand limits the reactor load-following operability. Recently advances in gas reactor designs have made load-following in gas reactors a feasible and achievable goal. Modern nuclear reactors (such as the MHTGR-350) are designed to withstand the structural stresses associated with load-following. Operating the reactor in load-following mode will inevitably reduce the load factor. Although a higher load factor means more revenue and the best usage of the fuel, following the variable energy demand will increase the maximum achievable fraction of generated nuclear power, being no more limited to the base load power generation. In general, if low impact on material aging and the safe operability are demonstrated the plant economics will be minimally affected. In order to ensure that the system can be safely operated in a load-following mode an extensively study has been carried out. Some tests with linear coolant mass flow rate change demonstrated that the reactor behavior is suitable for operation in load following mode. This assumption has been confirmed by a 4 days load following transient test in which the reactor supplied the requested power with a negligible error.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1983

State vector formulation of substructure coupling for damped systems

A generalized substructure coupling procedure in state vector form is derived for a complex system with general viscous damping. This first-order differential equation formulation is used in order to permit complex substructure modes to be employed easily. The free-interface normal (complex) modes and rigid-body modes of the substructure are defined. Complex residual attachment modes which result from static approximation of neglected higher modes are derived. The motion of each substructure is represented by a selected set of component modes. Equations of interface compatibility are employed to obtain an independent set of system equations of motion. A new method which employs incomplete complex normal modes in conjunction with the complex residual attachment modes to account for the contribution of neglected higher order modes is derived. Examples are employed to indicate how a damped structure may be analyzed by including the effect of residual attachment modes. Numerical results indicate that the new component mode synthesis method provides sytem equations of motion with reduced number of degrees of feedom leading to more accurate approximations to the system frequencies and damping factors than are obtained by pure mode truncation.

Chung, Y.-T.↗
DOE OSTI2025

Laser wavelength dependence of laser imprint

In laser direct-drive inertial confinement fusion, laser imprint is one of the major causes of degradation in target performance through its seeding of hydrodynamic instabilities. Early experiments and simulations have shown that laser imprint could be mitigated with a longer laser wavelength because of its lower critical density and longer conduction zone. Building upon this work, we explore a scenario where the laser wavelength during the picket pulse differs from that of the main pulse in order to gain the benefit of reduced imprint, while avoiding losses in drive coupling efficiency or an increase in laser-plasma instabilities. A series of 2D radiation-hydrodynamic simulations, which test three different laser wavelengths for the picket pulse, has been performed, where the intensity of the picket pulse is adjusted in order to maintain the same implosion adiabat. A detailed analysis of the growth of the mass density modulations at the ablation front over a large range of mode numbers confirms that the laser imprint can be mitigated with a picket pulse operating at a longer wavelength than the main pulse because of the longer conduction zone and enhanced thermal smoothing. The amplitude of the ablation front modulations is found to be lower for all mode numbers, which reduces the seeding of the Rayleigh–Taylor instability without affecting the mode growth rates.

Fourier analysis↗
NASA NTRS1988

Axisymmetric perturbations of thin gaseous disks. I - Unstable convective modes and their consequences for the solar nebula

The axisymmetric perturbations of a thin, differentially rotating gas disk in which the vertical temperature stratification is superadiabatic are analyzed. The growth rate of adiabatic convective normal modes is calculated, departures from simple polytropic disk models are briefly discussed, and the detailed vertical structure of the eigenfunctions is analyzed. The stabilizing effect of radiative diffusion on convective modes is considered. It is found that rotation and compressibility tend to reduce the rate of growth of the disturbances, and that the growth rate increases without limit and is proportional to the square root of the radial wavenumber in polytropic equilibria. The maximum radial size of convective eddies scales like the square root of the degree of superadiabaticity times the size of the convective zone. There is significant convective penetration into radiatively stable layers only for the fundamental and low-order harmonic modes, whose vertical wavelength is comparable to the size of the convective layer.

Ruden, Steven P.↗
DOE OSTI2021

Role of microtearing mode in DIII-D and future high- β p core plasmas

This paper provides modeling evidence for a new variant of microturbulence instability that regulates the electron thermal fluxes in the transport barrier region of DIII-D "high β p " discharges. This scenario is being investigated as a possible scenario for future steady-state reactors, and thus it is essential to understand the physics processes which govern its performance in order to confidently extrapolate to future devices. Additional analysis examining the extrapolation of this scenario to ITER shows that while the new microtearing mode identified here is unlikely to be unstable in that regime (due to increased temperatures and reduced collisionality), the large shift in magnetic axis (Shafranov shift) that arises in these plasmas may still favor stabilization of common ballooning instabilities which peak on the low field side, but destabilize electrostatic slab modes on the high field side.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fracture length data for geothermal applications

Fracture lengths govern permeability and are unknowns in geothermal assessment. Along their lengths, fracture widths vary due to growth by linkage. Under the influence of diagenesis, narrow widths seal, breaking porosity continuity and reducing open length. The largest range of widths and thus susceptibility to fill occurs where fractures are linked by narrow segments. Outcrops of a geothermal target, Cambrian Potsdam quartz arenite, contain opening-mode fractures having lengths spanning five orders of magnitude from 0.082 mm to 17.9 m. Combined lengths measured at a range of scales can be described by power laws, but at a given image resolution, lengths are best fit by exponential functions. Owing to preferential sealing of small fractures, open fractures follow exponential functions, but values depend on rules for designating fractures as continuous. En échelon segments, offset 10 mm, are connected by narrow fractures or microfractures (hard linked) not evident on outcrop 1 m-elevation LiDAR or 30 m-height drone images. A rule that identifies where narrow and likely connected segments are located can yield lengths meaningful for flow simulation. Depending on diagenesis, continuity rules can halve or double average and maximum length values. Length values from outcrop for geothermal applications should be adjusted based on wellsite-specific diagenesis information.

15 GEOTHERMAL ENERGY↗
NASA NTRS1974

Transmission of sound through annular ducts of varying cross sections and sheared mean flow

An asymptotic solution is presented for the transmission and attenuation of acoustic waves in an annular duct of slowly varying cross section which carries a sheared mean flow. The analysis also takes into account the growth of the boundary layer as well as any slow variations in the acoustic liner properties. The problem is reduced to the solution of a first-order ordinary differential equation for the amplitude. The analysis is used to estimate the effects of the variations of the duct cross section and boundary layer growth on the different acoustic modes.

Nayfeh, A. H.↗
NASA NTRS1980

Identification of flexible aircraft from flight data

This paper describes a frequency-domain method for the estimation of the aeroelastic equations of motion of an aircraft using discrete sinusoidal inputs to the control surface actuators. The theory of estimation of the signal inputs and outputs is described as well as the method of obtaining accurate transfer functions from the resulting Bode plots. Methods of quantifying the accuracy of the results from more than one control input are discussed along with ways to obtain state-space and reduced-order models from the transfer function estimates. Presented is a summary of the results obtained in using the technique to develop a model of the USAF Total In-Flight Simulator (TIFS) aircraft which includes five symmetrical elastic modes of motion of the aircraft.

Eulrich, B. J.↗
NASA NTRS1989

Control/structure interaction methods for space station power systems

The Structural Dynamics Research Corporation and the NASA Lewis Research Center have been working together to develop tools and methods for the analysis of control/structure interaction problems related to the space station power systems. Flexible modes of the solar arrays below 0.1 Hz, suggest that even for relatively slow control systems, the potential for control/structure interaction exists. The emphasis of the effort has been to develop tools which couple NASTRAN's powerful capabilities in structural dynamics with EASY5's powerful capabilities in control systems analysis. One product is an interface software package called CO-ST-IN for COntrol-STructure-INteraction. CO-ST-IN acts to translate data between NASTRAN and EASY5, facilitating the analysis of complex coupled problems. Interfaces to SDRC I-DEAS and MATRIXx are also offered. Beside transferring standard modal information, CO-ST-IN implements a number of advanced methods. These include a modal ordering algorithm that helps eliminate uncontrollable or unobservable modes from the analysis, an implementation of the more accurate mode acceleration algorithm for recovery of element forces and stresses directly in EASY5 and an implementation of fixed interface modes in NASTRAN, which reduces the error in the closed-loop model due to the use of truncated mode sets.

Blelloch, Paul↗
DOE OSTI2023

Mshock HED Council Request FY24 [Slides]

ICF and HEDP experiments experience complicated loadings including successive shocks, but modeling and analytic work has mostly engaged the simpler shock or shock-reshock case. This is because the co propagating case is difficult to achieve with conventional (non-HED) drivers. Successive shocks will be challenging to model in BHR or modal model interpretations: 1) Ex: Consider the case where long wavelength modes re-invert coherently, but short wavelength nonlinear modes are spun up turbulently. 2) Necessitates more advanced diagnostics and analysis than simple “mixing layers.” High resolution spectral information (concomitant with Ω-EP campaigns) and higher-order moment analyses will be necessary to stress and validate the reduced BHR-type models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1985

UNAERO: A package of FORTRAN subroutines for approximating unsteady aerodynamics in the time domain

This report serves as an instruction and maintenance manual for a collection of CDC CYBER FORTRAN IV subroutines for approximating the unsteady aerodynamic forces in the time domain. The result is a set of constant-coefficient first-order differential equations that approximate the dynamics of the vehicle. Provisions are included for adjusting the number of modes used for calculating the approximations so that an accurate approximation is generated. The number of data points at different values of reduced frequency can also be varied to adjust the accuracy of the approximation over the reduced-frequency range. The denominator coefficients of the approximation may be calculated by means of a gradient method or a least-squares approximation technique. Both the approximation methods use weights on the residual error. A new set of system equations, at a different dynamic pressure, can be generated without the approximations being recalculated.

Dunn, H. J.↗
NASA NTRS1989

Algorithm For Optimal Control Of Large Structures

Cost of computation appears competitive with other methods. Problem to compute optimal control of forced response of structure with n degrees of freedom identified in terms of smaller number, r, of vibrational modes. Article begins with Hamilton-Jacobi formulation of mechanics and use of quadratic cost functional. Complexity reduced by alternative approach in which quadratic cost functional expressed in terms of control variables only. Leads to iterative solution of second-order time-integral matrix Volterra equation of second kind containing optimal control vector. Cost of algorithm, measured in terms of number of computations required, is of order of, or less than, cost of prior algoritms applied to similar problems.

Salama, Moktar A.↗
NASA NTRS1990

Cyclotron maser and plasma wave growth in magnetic loops

Cyclotron maser and plasma wave growth which results from electrons accelerated in magnetic loops are studied. The evolution of the accelerated electron distribution is determined by solving the kinetic equation including Coulomb collisions and magnetic convergence. It is found that for modest values of the column depth of the loop the growth rates of instabilities are significantly reduced and that the reduction is much larger for the cyclotron modes than for the plasma wave modes. The large decrease in the growth rate with column depth suggests that solar coronal densities must be much lower than commonly accepted in order for the cyclotron maser to operate. The density depletion has to be similar to that which occurs during auroral kilometric radiation events in the magnetosphere. The resulting distributions are much more complicated than the idealized distributions used in many theoretical studies, but the fastest growing mode can still simply be determined by the ratio of electron plasma to gyrofrequency, U=omega(sub p)/Omega(sub e). However, the dominant modes are different than for the idealized situations with growth of the z-mode largest for U approximately less than 0.5, and second harmonic x-mode (s=2) or fundamental o-mode (s=1) the dominant modes for 0.5 approximately less than U approximately less than 1. The electron distributions typically contain more than one inverted feature which could give rise to wave growth. It is shown that this can result in simultaneous amplification of more than one mode with each mode driven by a different feature and can be observed, for example, by differences in the rise times of the right and left circularly polarized components of the associated spike bursts.

Hamilton, Russell J.↗