Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “process damping”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

External noise-induced transitions in a current-biased Josephson junction

We investigate noise-induced transitions in a current-biased and weakly damped Josephson junction in the presence of multiplicative noise. By using the stochastic averaging procedure, the averaged amplitude equation describing dynamic evolution near a constant phase difference is derived. Numerical results show that a stochastic Hopf bifurcation between an absorbing and an oscillatory state occurs. This means the external controllable noise triggers a transition into the non-zero junction voltage state. With the increase of noise intensity, the stationary probability distribution peak shifts and is characterised by increased width and reduced height. And the different transition rates are shown for large and small bias currents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Eigenvalue-Based Micromagnetic Analysis of Switching in Spin-Torque-Driven Structures

We present an eigenvalue-based approach for studying the magnetization dynamics in magnetic nanostructures driven by spintronic excitations, such as spin-transfer torque and spin-orbit torque. The approach represents the system dynamics in terms of normal oscillation modes (eigenstates) with corresponding complex eigenfrequencies. The dynamics is driven by a small number of active eigenstates and often considering just a single eigenstate is sufficient. We develop a perturbation theory that provides semianalytical dynamic solutions by using eigenstates for the case in the absence of damping and spintronic excitations as a basis. The approach provides useful insights into dynamics in such systems and allows solving several difficulties in their modeling, such as extracting the switching current in magnetic random-access memories and understanding switching mechanisms. We show that the presented approach directly predicts the critical switching current, i.e., switching current for an infinite time. The approach also provides solutions for the switching dynamics allowing the switching current to be obtained for a finite switching time, provided that the system symmetry is broken, e.g., by tilting the polarizer, so that switching by a finite pulse is possible.

42 ENGINEERING↗
DOE OSTI2024

Plasma-based particle sources

High-brightness particle beams generated by advanced accelerator concepts have the potential to become an essential part of future accelerator technology. In particular, high-gradient accelerators can generate and rapidly accelerate particle beams to relativistic energies. The rapid acceleration and strong confining fields can minimize irreversible detrimental effects to the beam brightness that occur at low beam energies, such as emittance growth or pulse elongation caused by space charge forces. Due to the high accelerating gradients, these novel accelerators are also significantly more compact than conventional technology. Advanced accelerators can be extremely variable and are capable of generating particle beams with vastly different properties using the same driver and setup with only modest changes to the interaction parameters. So far, efforts have mainly been focused on the generation of electron beams, but there are concepts to extend the sources to generate spin-polarized electron beams or positron beams. The beam parameters of these particle sources are largely determined by the injection and subsequent acceleration processes. Although, over the last decade there has been significant progress, the sources are still lacking a sufficiently high 6-dimensional (D) phase-space density that includes small transverse emittance, small energy spread and high charge, and operation at high repetition rate. This is required for future particle colliders with a sufficiently high luminosity or for more near-term applications, such as enabling the operation of free-electron lasers (FELs) in the X-ray regime. Major research and development efforts are required to address these limitations in order to realize these approaches for a front-end injector for a future collider or next-generation light sources. In particular, this includes methods to control and manipulate the phase-space and spin degrees-of-freedom of ultrashort plasma-based electron bunches with high accuracy, and methods that increase efficiency and repetition rate. These efforts also include the development of high-resolution diagnostics, such as full 6D phase-space measurements, beam polarimetry and high-fidelity simulation tools. A further increase in beam luminosity can be achieve through emittance damping. Emittance cooling via the emission of synchrotron radiation using current technology requires kilometer-scale damping rings. For future colliders, the damping rings might be replaced by a substantially more compact plasma-based approach. Here, plasma wigglers with significantly stronger magnetic fields are used instead of permanent-magnet based wigglers to achieve similar damping performance but over a two orders of magnitude reduced length.

43 PARTICLE ACCELERATORS↗
DOE OSTI2011

INCORPORATING KINETIC PHYSICS INTO A TWO-FLUID SOLAR-WIND MODEL WITH TEMPERATURE ANISOTROPY AND LOW-FREQUENCY ALFVEN-WAVE TURBULENCE

We develop a one-dimensional solar-wind model that includes separate energy equations for the electrons and protons, proton temperature anisotropy, collisional and collisionless heat flux, and an analytical treatment of low-frequency, reflection-driven, Alfven-wave (AW) turbulence. To partition the turbulent heating between electron heating, parallel proton heating, and perpendicular proton heating, we employ results from the theories of linear wave damping and nonlinear stochastic heating. We account for mirror and oblique firehose instabilities by increasing the proton pitch-angle scattering rate when the proton temperature anisotropy exceeds the threshold for either instability. We numerically integrate the equations of the model forward in time until a steady state is reached, focusing on two fast-solar-wind-like solutions. These solutions are consistent with a number of observations, supporting the idea that AW turbulence plays an important role in the origin of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

Gel stretch method: a new method to measure constitutive properties of cardiac muscle cells

Diastolic dysfunction is an important cause of congestive heart failure; however, the basic mechanisms causing diastolic congestive heart failure are not fully understood, especially the role of the cardiac muscle cell, or cardiocyte, in this process. Before the role of the cardiocyte in this pathophysiology can be defined, methods for measuring cardiocyte constitutive properties must be developed and validated. Thus this study was designed to evaluate a new method to characterize cardiocyte constitutive properties, the gel stretch method. Cardiocytes were isolated enzymatically from normal feline hearts and embedded in a 2% agarose gel containing HEPES-Krebs buffer and laminin. This gel was cast in a shape that allowed it to be placed in a stretching device. The ends of the gel were held between a movable roller and fixed plates that acted as mandibles. Distance between the right and left mandibles was increased using a stepper motor system. The force applied to the gel was measured by a force transducer. The resultant cardiocyte strain was determined by imaging the cells with a microscope, capturing the images with a CCD camera, and measuring cardiocyte and sarcomere length changes. Cardiocyte stress was characterized with a finite-element method. These measurements of cardiocyte stress and strain were used to determine cardiocyte stiffness. Two variables affecting cardiocyte stiffness were measured, the passive elastic spring and viscous damping. The passive spring was assessed by increasing the force on the gel at 1 g/min, modeling the resultant stress vs. strain relationship as an exponential [sigma = A/k(ekepsilon - 1)]. In normal cardiocytes, A = 23.0 kN/m2 and k = 16. Viscous damping was assessed by examining the loop area between the stress vs. strain relationship during 1 g/min increases and decreases in force. Normal cardiocytes had a finite loop area = 1.39 kN/m2, indicating the presence of viscous damping. Thus the gel stretch method provided accurate measurements of cardiocyte constitutive properties. These measurements have allowed the first quantitative assessment of passive elastic spring properties and viscous damping in normal mammalian cardiocytes.

Non-NASA Center↗
DOE OSTI2024

Superconducting microwave cavities and qubits for quantum information systems

Superconducting microwave cavities featuring ultrahigh Q-factors, which measure the efficiency of energy storage in relation to energy loss in a system, are revolutionizing quantum computing by providing long coherence times exceeding 1 ms, crucial for the development of scalable multi-qubit quantum systems with low error rates. Here, in this work, we provide an in-depth analysis of recent advances in ultrahigh Q-factor cavities, integration of Josephson junction-based qubits, and bosonic-encoded qubits in 3D cavities. We examine the sources of quantum state dephasing caused by damping and noise mechanisms in cavities and qubits, highlighting the critical challenges that need to be addressed to achieve even higher coherence times. We critically survey the latest progress made in implementing single 3D qubits using superconducting materials, normal metals, and multi-qubit and multi-state quantum systems. Our work sheds light on the promising future of this research area, including novel materials for cavities and qubits, modes with nontrivial topological properties, error correction techniques for bosonic qubits, and new light–matter interaction effects.

3D printing↗
NASA NTRS2011

Toward Supersonic Retropropulsion CFD Validation

This paper begins the process of verifying and validating computational fluid dynamics (CFD) codes for supersonic retropropulsive flows. Four CFD codes (DPLR, FUN3D, OVERFLOW, and US3D) are used to perform various numerical and physical modeling studies toward the goal of comparing predictions with a wind tunnel experiment specifically designed to support CFD validation. Numerical studies run the gamut in rigor from code-to-code comparisons to observed order-of-accuracy tests. Results indicate that this complex flowfield, involving time-dependent shocks and vortex shedding, design order of accuracy is not clearly evident. Also explored is the extent of physical modeling necessary to predict the salient flowfield features found in high-speed Schlieren images and surface pressure measurements taken during the validation experiment. Physical modeling studies include geometric items such as wind tunnel wall and sting mount interference, as well as turbulence modeling that ranges from a RANS (Reynolds-Averaged Navier-Stokes) 2-equation model to DES (Detached Eddy Simulation) models. These studies indicate that tunnel wall interference is minimal for the cases investigated; model mounting hardware effects are confined to the aft end of the model; and sparse grid resolution and turbulence modeling can damp or entirely dissipate the unsteadiness of this self-excited flow.

Kleb, Bil↗
DOE OSTI2016

DENSITY FLUCTUATIONS UPSTREAM AND DOWNSTREAM OF INTERPLANETARY SHOCKS

Interplanetary (IP) shocks as typical large-scale disturbances arising from processes such as stream–stream interactions or Interplanetary Coronal Mass Ejection (ICME) launching play a significant role in the energy redistribution, dissipation, particle heating, acceleration, etc. They can change the properties of the turbulent cascade on shorter scales. We focus on changes of the level and spectral properties of ion flux fluctuations upstream and downstream of fast forward oblique shocks. Although the fluctuation level increases by an order of magnitude across the shock, the spectral slope in the magnetohydrodynamic range is conserved. The frequency spectra upstream of IP shocks are the same as those in the solar wind (if not spoiled by foreshock waves). The spectral slopes downstream are roughly proportional to the corresponding slopes upstream, suggesting that the properties of the turbulent cascade are conserved across the shock; thus, the shock does not destroy the shape of the spectrum as turbulence passes through it. Frequency spectra downstream of IP shocks often exhibit “an exponential decay” in the ion kinetic range that was earlier reported at electron scales in the solar wind or at ion scales in the interstellar medium. We suggest that the exponential shape of ion flux spectra in this range is caused by stronger damping of the fluctuations in the downstream region.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Digital active material processing platform effort (DAMPER), SBIR phase 2

Applied Technology Associates, Inc., (ATA) has demonstrated that inertial actuation can be employed effectively in digital, active vibration isolation systems. Inertial actuation involves the use of momentum exchange to produce corrective forces which act directly on the payload being actively isolated. In a typical active vibration isolation system, accelerometers are used to measure the inertial motion of the payload. The signals from the accelerometers are then used to calculate the corrective forces required to counteract, or 'cancel out' the payload motion. Active vibration isolation is common technology, but the use of inertial actuation in such systems is novel, and is the focus of the DAMPER project. A May 1991 report was completed which documented the successful demonstration of inertial actuation, employed in the control of vibration in a single axis. In the 1 degree-of-freedom (1DOF) experiment a set of air bearing rails was used to suspend the payload, simulating a microgravity environment in a single horizontal axis. Digital Signal Processor (DSP) technology was used to calculate in real time, the control law between the accelerometer signals and the inertial actuators. The data obtained from this experiment verified that as much as 20 dB of rejection could be realized by this type of system. A discussion is included of recent tests performed in which vibrations were actively controlled in three axes simultaneously. In the three degree-of-freedom (3DOF) system, the air bearings were designed in such a way that the payload is free to rotate about the azimuth axis, as well as translate in the two horizontal directions. The actuator developed for the DAMPER project has applications beyond payload isolation, including structural damping and source vibration isolation. This report includes a brief discussion of these applications, as well as a commercialization plan for the actuator.

Blackburn, John↗
NASA NTRS2003

Fast-Running Aeroelastic Code Based on Unsteady Linearized Aerodynamic Solver Developed

The NASA Glenn Research Center has been developing aeroelastic analyses for turbomachines for use by NASA and industry. An aeroelastic analysis consists of a structural dynamic model, an unsteady aerodynamic model, and a procedure to couple the two models. The structural models are well developed. Hence, most of the development for the aeroelastic analysis of turbomachines has involved adapting and using unsteady aerodynamic models. Two methods are used in developing unsteady aerodynamic analysis procedures for the flutter and forced response of turbomachines: (1) the time domain method and (2) the frequency domain method. Codes based on time domain methods require considerable computational time and, hence, cannot be used during the design process. Frequency domain methods eliminate the time dependence by assuming harmonic motion and, hence, require less computational time. Early frequency domain analyses methods neglected the important physics of steady loading on the analyses for simplicity. A fast-running unsteady aerodynamic code, LINFLUX, which includes steady loading and is based on the frequency domain method, has been modified for flutter and response calculations. LINFLUX, solves unsteady linearized Euler equations for calculating the unsteady aerodynamic forces on the blades, starting from a steady nonlinear aerodynamic solution. First, we obtained a steady aerodynamic solution for a given flow condition using the nonlinear unsteady aerodynamic code TURBO. A blade vibration analysis was done to determine the frequencies and mode shapes of the vibrating blades, and an interface code was used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor was used to interpolate the mode shapes from the structural dynamic mesh onto the computational dynamics mesh. Then, we used LINFLUX to calculate the unsteady aerodynamic forces for a given mode, frequency, and phase angle. A postprocessor read these unsteady pressures and calculated the generalized aerodynamic forces, eigenvalues, and response amplitudes. The eigenvalues determine the flutter frequency and damping. As a test case, the flutter of a helical fan was calculated with LINFLUX and compared with calculations from TURBO-AE, a nonlinear time domain code, and from ASTROP2, a code based on linear unsteady aerodynamics.

Reddy, T. S. R.↗
NASA NTRS1996

Shape deformation of the organ of Corti associated with length changes of outer hair cell

Cochlear outer hair cells (OHC) are commonly assumed to function as mechanical effectors as well as sensory receptors in the organ of Corti (OC) of the inner ear. OHC in vitro and in organ explants exhibit mechanical responses to electrical, chemical or mechanical stimulation which may represent an aspect of their effector process that is expected in vivo. A detailed description, however, of an OHC effector operation in situ is still missing. Specifically, little is known as to how OHC movements influence the geometry of the OC in situ. Previous work has demonstrated that the motility of isolated OHCs in response to electrical stimulation and to K(+)-gluconate is probably under voltage control and causes depolarisation (shortening) and hyperpolarization (elongation). This work was undertaken to investigate if the movements that were observed in isolated OHC, and which are induced by ionic stimulation, could change the geometry of the OC. A synchronized depolarization of OHC was induced in guinea pig cochleae by exposing the entire OC to artificial endolymph (K+). Subsequent morphometry of mid-modiolar sections from these cochleae revealed that the distance between the basilar membrane (BM) and the reticular lamina (RL) had decreased considerably. Furthermore, in the three upper turns OHC had significantly shortened in all rows. The results suggest that OHC can change their length in the organ of Corti (OC) thus deforming the geometry of the OC. The experiments reveal a tonic force generation within the OC that may change the position of RL and/or BM, contribute to damping, modulate the BM-RL-distance and control the operating points of RL and sensory hair bundles. Thus, the results suggest active self-adjustments of cochlear mechanics by slow OHC length changes. Such mechanical adjustments have recently been postulated to correspond to timing elements of animal communication, speech or music.

Non-NASA Center↗
DOE OSTI2021

Electron acceleration at oblique angles via stimulated Raman scattering at laser irradiance >10 16 W cm –2 μm 2

The generation of hot, directional electrons via laser-driven stimulated Raman scattering (SRS) is a topic of great importance in inertial confinement fusion (ICF) schemes. Little recent research has been dedicated to this process at high laser intensity, in which back, side, and forward scatter simultaneously occur in high energy density plasmas, of relevance to, for example, shock ignition ICF. In this paper, we present an experimental and particle-in-cell (PIC) investigation of hot electron production from SRS in the forward and near-forward directions from a single speckle laser of wavelength λ 0 = 1.053μm, peak laser intensities in the range I 0 = 0.2–1.0×10 17 Wcm –2 and target electron densities between ne = 0.3–1.6% n c , where n c is the plasma critical density. As the intensity and density are increased, the hot electron spectrum changes from a sharp cutoff to an extended spectrum with a slope temperature T = 34 ± 1 keV and maximum measured energy of 350 keV experimentally. Multidimensional PIC simulations indicate that the high energy electrons are primarily generated from SRS-driven electron plasma wave phase fronts with k vectors angled ~50° with respect to the laser axis. These results are consistent with analytical arguments that the spatial gain is maximized at an angle which balances the tendency for the growth rate to be larger for larger scattered light wave angles until the kinetic damping of the plasma wave becomes important. The efficiency of generated high energy electrons drops significantly with a reduction in either laser intensity or target electron density, which is a result of the rapid drop in growth rate of Raman scattering at angles in the forward direction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Theoretical Plasma Physics

Lattice Boltzmann algorithms are a mesoscopic method to solve problems in nonlinear physics which are highly parallelized – unlike the direction solution of the original problem. These methods are applied to both fluid and magnetohydrodynamic turbulence. By introducing entropic constraints one can enforce the positive definiteness of the distribution functions and so be able to simulate fluids at high Reynolds numbers without numerical instabilities. By introducing a vector distribution function for the magnetic field one can enforce the divergence free condition on the magnetic field automatically, without the need of divergence cleaning as needed in most direct numerical solutions of the resistive magnetohydrodynamic equations. The principal reason for the high parallelization of lattice Boltzmann codes is that they consist of a kinetic collisional relaxation step (which is purely local) followed by a simple shift of the relaxed data to neighboring lattice sites. In large eddy simulations, the closure schemes are highly nonlocal – the most famous of these schemes is that due to Smagorinsky. Under a lattice Boltzmann representation the Smagorinsky closure is purely local – being simply a particular moment on the perturbed distribution fucntions. After nonlocal fluid moment models were discovered to represent Landau damping, it was found possible to model these fluid models using an appropriate lattice Boltzmann algorithm. The close to ideal parallelization of the lattice Boltzmann codes permitted us to be Gordon Bell finalists on using the Earth Simulation in Japan. We have also been involved in the radio frequency propagation of waves into a tokamak and into a spherical overdense tokamak plasma. Initially we investigated the use of a quasi-optical grill for the launching of lower hybrid waves into a tokamak. It was found that the conducting walls do not prevent the rods from being properly irradiated, the overloading of the quasi-optical grill is not severe with the electric field only being about three times higher than in the ideal case. Moreover, the quasi-optical grill was significantly fewer structural elements that the multijunction grill. Nevertheless there has not been much interest from experimental fusion groups to implementing these structures. Hence we have returned to optimizing the multijunction grill so that the large number of coupling matrix elements can be efficiently evaluated using symmetry arguments. In overdense plasmas, the standard electromagnetic waves cannot propagate into the plasma center, but are reflected at the plasma edge. By optimizing mode conversion processes (in particular, the O-X-B wave propagation of Ordinary Mode converting to an Extraordinary mode which then converts into an electrostatic Bernstein wave) one can excite within the plasma an electrostatic Bernstein wave that does not suffer density cutoffs and is absorbed on the electron cyclotron harmonics. Finally we have started looking at other mesoscopic lattice algorithms that involve unitary collision and streaming steps. Because these algorithms are unitary they can be run on quantum computers when they become available – unlike their computational cousin of lattice Boltzmann which is a purely classical code. These quantum lattice gas algorithms have been tested successfully on exact analytic soliton collision solution. These calculations are hoped to be able to study Bose Einstein condensed atomic gases and their ground states in an optical lattice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Symbolic generation of elastic rotor blade equations using a FORTRAN processor and numerical study on dynamic inflow effects on the stability of helicopter rotors

The process of performing an automated stability analysis for an elastic-bladed helicopter rotor is discussed. A symbolic manipulation program, written in FORTRAN, is used to aid in the derivation of the governing equations of motion for the rotor. The blades undergo coupled bending and torsional deformations. Two-dimensional quasi-steady aerodynamics below stall are used. Although reversed flow effects are neglected, unsteady effects, modeled as dynamic inflow are included. Using a Lagrangian approach, the governing equations are derived in generalized coordinates using the symbolic program. The program generates the steady and perturbed equations and writes into subroutines to be called by numerical routines. The symbolic program can operate on both expressions and matrices. For the case of hovering flight, the blade and dynamic inflow equations are converted to equations in a multiblade coordinate system by rearranging the coefficients of the equations. For the case of forward flight, the multiblade equations are obtained through the symbolic program. The final multiblade equations are capable of accommodating any number of elastic blade modes. The computer implementation of this procedure consists of three stages: (1) the symbolic derivation of equations; (2) the coding of the equations into subroutines; and (3) the numerical study after identifying mass, damping, and stiffness coefficients. Damping results are presented in hover and in forward flight with and without dynamic inflow effects for various rotor blade models, including rigid blade lag-flap, elastic flap-lag, flap-lag-torsion, and quasi-static torsion. Results from dynamic inflow effects which are obtained from a lift deficiency function for a quasi-static inflow model in hover are also presented.

Reddy, T. S. R.↗
DOE OSTI2013

MECHANISM FOR EXCITING PLANETARY INCLINATION AND ECCENTRICITY THROUGH A RESIDUAL GAS DISK

According to the theory of Kozai resonance, the initial mutual inclination between a small body and a massive planet in an outer circular orbit is as high as {approx}39. Degree-Sign 2 for pumping the eccentricity of the inner small body. Here we show that with the presence of a residual gas disk outside two planetary orbits, the inclination can be reduced to as low as a few degrees. The presence of the disk changes the nodal precession rates and directions of the planet orbits. At the place where the two planets achieve the same nodal processing rate, vertical secular resonance (VSR) occurs so that the mutual inclination of the two planets will be excited, which might further trigger the Kozai resonance between the two planets. However, in order to pump an inner Jupiter-like planet, the conditions required for the disk and the outer planet are relatively strict. We develop a set of evolution equations, which can fit the N-body simulation quite well but can be integrated within a much shorter time. By scanning the parameter spaces using the evolution equations, we find that a massive planet (10 M{sub J} ) at 30 AU with an inclination of 6 Degree-Sign to a massive disk (50 M{sub J} ) can finally enter the Kozai resonance with an inner Jupiter around the snowline. An inclination of 20 Degree-Sign of the outer planet to the disk is required for flipping the inner one to a retrograde orbit. In multiple planet systems, the mechanism can happen between two nonadjacent planets or can inspire a chain reaction among more than two planets. This mechanism could be the source of the observed giant planets in moderate eccentric and inclined orbits, or hot Jupiters in close-in, retrograde orbits after tidal damping.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

REVISITING ACCELERATION OF CHARGED GRAINS IN MAGNETOHYDRODYNAMIC TURBULENCE

We study the acceleration of charged grains by magnetohydrodynamic (MHD) turbulence in the interstellar medium (ISM). We begin with revisiting gyroresonance acceleration by taking into account the fluctuations of grain guiding center along a uniform magnetic field (i.e., nonlinear theory-NLT). We calculate grain velocities due to gyroresonance by fast MHD modes using the NLT for different phases of the ISM and compare them with results obtained using quasi-linear theory (QLT). We find for the parameters applicable to the typical ISM phases that the fluctuations of the grain guiding center reduce grain velocities by less than 15%, but they can be important for more special circumstances. We confirm that large grains can be accelerated to super-Alfvenic velocities through gyroresonance. For such super-Alfvenic grains, we investigate the effect of further acceleration via transit-time damping (TTD) by fast modes. We find that due to the broadening of the resonance condition in the NLT, the TTD acceleration is not only important for the cosines of grain pitch angle relative to the magnetic field {mu} > V{sub A}/v, but also for {mu} < V{sub A}/v where v is the grain velocity and V{sub A} is the Alfven speed. We show that the TTD acceleration is dominant over the gyroresonance for large grains, and can increase substantially grain velocities induced by gyroresonance acceleration. We quantify another stochastic acceleration mechanism arising from low-frequency Alfven waves. We discuss the range of applicability of the mechanisms and their implications.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

Evaluation of Aeroservoelastic Effects on Flutter

The HSCT Flight Controls Group is developing a longitudinal control law, known as Gamma-dot / V, for the NASA HSR program. Currently, this control law is based on a quasi-steady aeroelastic (QSAE) model of the vehicle. This control law was implemented into the p-k flutter analysis process for closed loop aeroservoelastic analysis. The available flexible models, developed for the TCA aeroelastic analysis, were used to assess the effect of control laws on flutter at several different Mach numbers and mass conditions. Significant structures and flight control system interaction was observed during the initial assessment. Figures 1 and 2 present a summary of the effect of total closed loop gain and phase on flutter mechanisms, based on ideal sensors and real sensors, for Mach 0.95 and mass M02 condition. Control laws based on ideal sensors gave rise to increased coupling between the rigid body short period mode and the first symmetric elastic mode. This reduced the stability margins for the first elastic mode and does not meet the required 6 dB gain margin requirement. The effect of "real" sensors significantly increased the structures and control system interactions. This caused the elastic,modes to be highly unstable throughout most of the flight envelope. State-space models were developed for several conditions and then MATLAB program was used for the aeroservoelastic stability analysis. These results provided an independent verification of the p-k flutter analysis findings. Good overall agreement was observed between the p-k flutter analysis and state-space model results for both damping and frequency comparisons. These results are also included in this document.

Nagaraja, K. S.↗
DOE OSTI2013

SUPERMASSIVE BLACK HOLE FORMATION AT HIGH REDSHIFTS VIA DIRECT COLLAPSE: PHYSICAL PROCESSES IN THE EARLY STAGE

We use numerical simulations to explore whether direct collapse can lead to the formation of supermassive black hole (SMBH) seeds at high redshifts. Using the adaptive mesh refinement code ENZO, we follow the evolution of gas within slowly tumbling dark matter (DM) halos of M{sub vir} {approx} 2 Multiplication-Sign 10{sup 8} M{sub Sun} and R{sub vir} {approx} 1 kpc. For our idealized simulations, we adopt cosmologically motivated DM and baryon density profiles and angular momentum distributions. Our principal goal is to understand how the collapsing flow overcomes the centrifugal barrier and whether it is subject to fragmentation which can potentially lead to star formation, decreasing the seed SMBH mass. We find that the collapse proceeds from inside out and leads either to a central runaway or to off-center fragmentation. A disk-like configuration is formed inside the centrifugal barrier, growing via accretion. For models with a more cuspy DM distribution, the gas collapses more and experiences a bar-like perturbation and a central runaway on scales of {approx}< 1-10 pc. We have followed this inflow down to {approx}10{sup -4} pc ({approx}10 AU), where it is estimated to become optically thick. The flow remains isothermal and the specific angular momentum, j, is efficiently transferred by gravitational torques in a cascade of nested bars. This cascade is triggered by finite perturbations from the large-scale mass distribution and by gas self-gravity, and supports a self-similar, disk-like collapse where the axial ratios remain constant. The mass accretion rate shows a global minimum on scales of {approx}1-10 pc at the time of the central runaway. In the collapsing phase, virial supersonic turbulence develops and fragmentation is damped. Models with progressively larger initial DM cores evolve similarly, but the timescales become longer. In models with more organized initial rotation-when the rotation of spherical shells is constrained to be coplanar-a torus forms on scales {approx}20-50 pc outside the disk, and appears to be supported by turbulent motions driven by accretion from the outside. The overall evolution of the models depends on the competition between two timescales, corresponding to the onset of the central runaway and of off-center fragmentation. In models with less organized rotation-when the rotation of spherical shells is randomized (but the total angular momentum remains unchanged)-the torus is greatly weakened, the central accretion timescale is shortened, and off-center fragmentation is suppressed-triggering the central runaway even in previously ''stable'' models. The resulting seed SMBH masses is found in the range M{sub .} {approx} 2 Multiplication-Sign 10{sup 4} M{sub Sun }-2 Multiplication-Sign 10{sup 6} M{sub Sun }, substantially higher than the mass range of Population III remnants. We argue that the above upper limit on M{sub .} appears to be more realistic, and lies close to the cutoff mass of detected SMBHs. Corollaries of this model include a possible correlation between SMBH and DM halo masses, and similarity between the SMBH and halo mass functions, at time of formation.

79 ASTRONOMY AND ASTROPHYSICS↗