SEARCH · Engineering Papers

Results for “local motion stabilization”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Convergence Acceleration of the Navier-Stokes Equations Through Time-Derivative Preconditioning

Chorin's method of artificial compressibility is extended to both compressible and incompressible fluids by using physical arguments to define artificial fluid properties that make up a local preconditioning matrix. In particular, perturbation expansions are used to provide appropriate temporal derivatives for the equations of motion at both low speeds and low Reynolds numbers. These limiting forms are then combined into a single function that smoothly merges into the physical time derivatives at high speeds so that the equations are left unchanged at transonic, high Reynolds number conditions. The effectiveness of the resulting preconditioning procedures for the Navier-Stokes equations is demonstrated for a wide speed and Reynolds number ranges by means of stability results and computational solutions. Nevertheless, the preconditioned equations sometimes fail to provide a solution for applications for which the non-preconditioned equations converge. Often this is because the reduced dissipation in the preconditioned equations results in an unsteady solution while the more dissipative non-preconditioned equations result in a steady state. Problems of this type represent a computational challenge; it is important to distinguish between non-convergence of algorithms, and the non-existence of steady state solutions.

Merkle, Charles L.↗

NASA NTRS1975

Effect of damping on excitability of high-order normal modes

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

Merchant, D. H.↗

NASA NTRS1975

The effects of localized damping on structural response

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

Merchant, D. H.↗

NASA NTRS1986

Jeans criterion in a turbulent medium

According to the classical Jeans analysis, all the molecular clouds of mass larger than a few 100 M(solar), size larger than about 1pc and kinetic temperature Tk less than 30K are gravitationally unstable. We have shown that in clouds supported by internal supersonic motions, local gravitational instabilities may appear within molecular clouds which are globally stable. The argument is threefold: (1) when the turbulent kinetic energy is included into the internal energy term, the virial equilibrium condition shows that molecular clouds such as those observed, which are gravitationally unstable according to the Jeans criterion, are indeed globally stable if supported by a turbulent velocity field of power spectrum steeper than 3; (2) 2D compressible hydrodynamical simulations show that a supersonic turbulent velocity field generates a turbulent pressure within clouds, the gradients of which stabilize the unstable scales (i.e., the largest scales and the cloud itself) against gravitational collapse; (3) an analysis similar to the Jeans approach but including the turbulent pressure gradient term, gives basically the same results as those given in (1). Clouds of mean density lower than a critical value are found to be stable even though more massive than their Jeans mass. In clouds of mean density larger than that critical value, the gravitational instability appears only over a range of scales smaller than the cloud size, the largest scales being stable. In practice, the observed mean densities are lower than this critical value: the observation of a small number of cores and stars of a few solar masses embedded in clouds of several hundred solar masses can only be understood in terms of small scale density fluctuations of large amplitude generated by the supersonic turbulence which would occasionally overtake the limit of gravitational stability.

Bonazzola, S.↗

NASA NTRS1990

Buoyancy-driven instabilities and the nonlinear breakup of a sheared magnetic layer

Motivated by problems concerning the storage and subsequent escape of the solar magnetic field, a study has been made of how a magnetic layer embedded in a convectively stable atmosphere evolves due to axisymmetric instabilities driven by magnetic buoyancy. The initial equilibrium consists of a toroidal field sheared by a weaker poloidal component. The linear stability problem is investigated for both ideal and resistive MHD, and the nonlinear evolution is followed by numerical integration of the equations of motion. In all cases, the instability is greatly affected by the distribution and strength of the poloidal field. In particular, both the horizontal and vertical scales of the motions are controlled by the location of the surface on which the poloidal field vanishes: the resonant surface. In the nonlinear regime, a resonant surface close to the interface between the magnetized and field-free fluid leads to the localization of the instability, so that only a fraction of the magnetic region is disrupted by the motions. By contrast, a deeply seated resonant surface leads to the complete disruption of the layer and to the formation of large, helical magnetic fragments whose identity is preserved for the entire simulation.

Cattaneo, Fausto↗

NASA NTRS1993

Two and three-dimensional prediffuser combustor studies with air-water mixture

Two- and three-dimensional gas turbine prediffuser-combustor sectors were experimentally studied under a number of mixture and flow conditions in a tunnel operating with a two-phase, air-liquid film-droplet mixture. It is concluded that water vaporization in the combustor causes changes in both local gas temperature and state of vitiation and reduces reaction rates. Substantial accumulation of water and water vapor takes place in pocket over the combustor volume, even when the air-water mixture is steady in time. The accuracy of determining combustor performance changes increases with a better knowledge of the state of the air-water mixture in the primary zone. To establish flame-out conditions it is considered to be necessary to combine the prediction of detailed flowfield and chemical activity with that of flame stability and motion characteristics.

Laing, Peter↗

NASA NTRS1993

Boundary layer receptivity and control

Receptivity processes initiate natural instabilities in a boundary layer. The instabilities grow and eventually break down to turbulence. Consequently, receptivity questions are a critical element of the analysis of the transition process. Success in modeling the physics of receptivity processes thus has a direct bearing on technological issues of drag reduction. The means by which transitional flows can be controlled is also a major concern: questions of control are tied inevitably to those of receptivity. Adjoint systems provide a highly effective mathematical method for approaching many of the questions associated with both receptivity and control. The long term objective is to develop adjoint methods to handle increasingly complex receptivity questions, and to find systematic procedures for deducing effective control strategies. The most elementary receptivity problem is that in which a parallel boundary layer is forced by time-harmonic sources of various types. The characteristics of the response to such forcing form the building blocks for more complex receptivity mechanisms. The first objective of this year's research effort was to investigate how a parallel Blasius boundary layer responds to general direct forcing. Acoustic disturbances in the freestream can be scattered by flow non-uniformities to produce Tollmien-Schlichting waves. For example, scattering by surface roughness is known to provide an efficient receptivity path. The present effort is directed towards finding a solution by a simple adjoint analysis, because adjoint methods can be extended to more complex problems. In practice, flows are non-parallel and often three-dimensional. Compressibility may also be significant in some cases. Recent developments in the use of Parabolized Stability Equations (PSE) offer a promising possibility. By formulating and solving a set of adjoint parabolized equations, a method for mapping the efficiency with which external forcing excites the three-dimensional motions of a non-parallel boundary layer was developed. The method makes use of the same computationally efficient formulation that makes the PSE currently so appealing. In the area of flow control, adjoint systems offer a powerful insight into the effect of control forces. One of the simplest control strategies for boundary layers involves the application of localized mean wall suction.

Hill, D. C.↗

DOE OSTI2014

Molecular reordering processes on ice (0001) surfaces from long timescale simulations

We report results of long timescale adaptive kinetic Monte Carlo simulations aimed at identifying possible molecular reordering processes on both proton-disordered and ordered (Fletcher) basal plane (0001) surfaces of hexagonal ice. The simulations are based on a force field for flexible molecules and span a time interval of up to 50 μs at a temperature of 100 K, which represents a lower bound to the temperature range of earth's atmosphere. Additional calculations using both density functional theory and an ab initio based polarizable potential function are performed to test and refine the force field predictions. Several distinct processes are found to occur readily even at this low temperature, including concerted reorientation (flipping) of neighboring surface molecules, which changes the pattern of dangling H-atoms, and the formation of interstitial defects by the downwards motion of upper-bilayer molecules. On the proton-disordered surface, one major surface roughening process is observed that significantly disrupts the crystalline structure. Despite much longer simulation time, such roughening processes are not observed on the highly ordered Fletcher surface which is energetically more stable because of smaller repulsive interaction between neighboring dangling H-atoms. However, a more localized process takes place on the Fletcher surface involving a surface molecule transiently leaving its lattice site. The flipping process provides a facile pathway of increasing proton-order and stabilizing the surface, supporting a predominantly Fletcher-like ordering of low-temperature ice surfaces. Our simulations also show that eventual proton-disordered patches on the surface may induce significant local reconstructions. Further, a subset of the molecules on the Fletcher surface are susceptible to forming interstitial defects which might provide active sites for various chemical reactions in the atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗

DOE OSTI2016

Evaluation and Compensation of Detector Solenoid Effects in the JLEIC

The JLEIC detector solenoid has a strong 3 T field in the IR area, and its tails extend over a range of several meters. One of the main effects of the solenoid field is coupling of the horizontal and vertical betatron motions which must be corrected in order to preserve the dynamical stability and beam spot size match at the IP. Additional effects include influence on the orbit and dispersion caused by the angle between the solenoid axis and the beam orbit. Meanwhile it affects ion polarization breaking the figure-8 spin symmetry. Crab dynamics further complicates the picture. All of these effects have to be compensated or accounted for. The proposed correction system is equivalent to the Rotating Frame Method. However, it does not involve physical rotation of elements. It provides local compensation of the solenoid effects independently for each side of the IR. It includes skew quadrupoles, dipole correctors and anti-solenoids to cancel perturbations to the orbit and linear optics. The skew quadrupoles and FFQ together generate an effect equivalent to adjustable rotation angle to do the decoupling task. Details of all of the correction systems are presented.

Wei, Guohui↗

NASA NTRS1995

An LDEF follow-on spacecraft concept

The successful flight, retrieval, and analyses of the Long Duration Exposure Facility (LDEF) experiments demonstrated the value of long duration space exposure for a broad spectrum of science and engineering investigations. The original LDEF was an excellent gravity gradient spacecraft, but because of its 9 m length and 9,700 kg mass it was difficult to manifest on the Shuttle, for either launch or retrieval, in conjunction with other payloads. This paper discusses an LDEF follow-on spacecraft concept whose short stowed length (approximately 3 m) greatly improves Shuttle manifesting opportunities while still providing very large surface area exposure for experiments. Deployable 'wings' on each end of the short, 'cylindrical' main body of this new spacecraft provide the means for gravity gradient stabilization while greatly increasing the spacecraft surface area. The center section of the spacecraft is oriented with the end faces of the twelve sided, 4.2 m diameter 'cylinder' perpendicular to the velocity vector thus providing large areas for experiments in both the ram and anti-ram directions as well as additional exposure area around the periphery of the cylinder. When deployed and properly oriented with the Shuttle's Remote Manipulator System (RMS), both wings of the spacecraft are oriented edge on to the direction of motion and lie in the plane which contains the local gravity vector. The relatively thin wings readily accommodate dual side exposure of glass plate stacks for cosmic ray detection. Flat surfaces mounted normal to and on the periphery of the wings provide additional areas in both the ram and anti-ram directions for cosmic dust, micrometeoroid, and orbital debris collection free of contamination from 'splatter' off secondary surfaces. The baseline concept provides enhancements not available on the original LDEF such as solar array generated electrical power and data telemetry. Status of the efforts to promote support for and ultimately space flight of this concept will be presented. Suggestions for improvements in the spacecraft design and proposed utilization are solicited.

Keller, Vernon↗

DOE OSTI2011

ORBITAL MOTION OF HR 8799 b, c, d USING HUBBLE SPACE TELESCOPE DATA FROM 1998: CONSTRAINTS ON INCLINATION, ECCENTRICITY, AND STABILITY

HR 8799 is currently the only multiple-planet system that has been detected with direct imaging, with four giant planets of masses 7-10 M{sub Jup} orbiting at large separations (15-68 AU) from this young late A star. Orbital motion provides insight into the stability and possible formation mechanisms of this planetary system. Dynamical studies can also provide constraints on the planets' masses, which help calibrate evolutionary models, yet measuring the orbital motion is a very difficult task because the long-period orbits (50-500 yr) require long time baselines and high-precision astrometry. This paper studies the three planets HR 8799b, c, and d in the archival data set of HR 8799 obtained with the Hubble Space Telescope (HST) NICMOS coronagraph in 1998. The detection of all three planets is made possible by a careful optimization of the Locally Optimized Combination of Images algorithm, and we used a statistical analysis of a large number of reduced images. This work confirms previous astrometry for planet b and presents new detections and astrometry for planets c and d. These HST images provide a ten-year baseline with the discovery images from 2008, and therefore offer a unique opportunity to constrain their orbital motion now. Recent dynamical studies of this system show the existence of a few possible stable solutions involving mean motion resonances (MMRs), where the interaction between c and d plays a major role. We study the compatibility of a few of these stable scenarios (1d:1c, 1d:2c, or 1d:2c:4d) with the new astrometric data from HST. In the hypothesis of a 1d:2c:4b MMR our best orbit fit is close to the stable solution previously identified for a three-planet system and involves low eccentricity for planet d (e{sub d} = 0.10) and moderate inclination of the system (i = 28.0 deg), assuming a coplanar system, circular orbits for b and c, and exact resonance with integer period ratios. Under these assumptions, we can place strong constraints on the inclination of the system (27.3-31.4 deg) and on the eccentricity for d e{sub d} < 0.46. Our results are robust to small departures from exact integer period ratios and consistent with previously published results based on dynamical studies for a three-planet system prior to the discovery of the fourth planet.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Sharpening of the front of the current through a cylindrical foil liner

Results are presented from experiments on the electromagnetic implosion of aluminum foil liners at the MIG generator with a current rise time of ≈80 ns. Plasma with a density of 10{sup 17} cm{sup –3} was preliminarily injected into the liner region by using a set of radial plasma guns. The Lorentz force J × B causes plasma acceleration in the radial direction. Since the magnetic field pressure is inversely proportional to the radius squared, the plasma displacement is maximum near the liner surface. As a result, plasma motion becomes two-dimensional, a gap appears between the plasma and the liner, and the generator current is switched over to the liner. The plasma velocity at the liner surface is close to the local Alfvén velocity, while the time during which the current is switched over to the liner is nearly equal to the ratio of the liner length to the Alfvén velocity. The proposed scheme allows one to decrease the rise time of the current through the liner to several nano-seconds and, as a result, to reduce the initial liner radius and improve the stability of liner implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

Experimental evaluation of active-member control of precision structures

The results of closed loop experiments that use piezoelectric active-members to control the flexible motion of a precision truss structure are described. These experiments are directed toward the development of high-performance structural systems as part of the Control/Structure Interaction (CSI) program at JPL. The focus of CSI activity at JPL is to develop the technology necessary to accurately control both the shape and vibration levels in the precision structures from which proposed large space-based observatories will be built. Structural error budgets for these types of structures will likely be in the sub-micron regime; optical tolerances will be even tighter. In order to achieve system level stability and local positioning at this level, it is generally expected that some form of active control will be required.

Fanson, James↗

NASA NTRS1991

Aeroacoustic and aerodynamic applications of the theory of nonequilibrium thermodynamics

Recent developments in the field of nonequilibrium thermodynamics associated with viscous flows are examined and related to developments to the understanding of specific phenomena in aerodynamics and aeroacoustics. A key element of the nonequilibrium theory is the principle of minimum entropy production rate for steady dissipative processes near equilibrium, and variational calculus is used to apply this principle to several examples of viscous flow. A review of nonequilibrium thermodynamics and its role in fluid motion are presented. Several formulations are presented of the local entropy production rate and the local energy dissipation rate, two quantities that are of central importance to the theory. These expressions and the principle of minimum entropy production rate for steady viscous flows are used to identify parallel-wall channel flow and irrotational flow as having minimally dissipative velocity distributions. Features of irrotational, steady, viscous flow near an airfoil, such as the effect of trailing-edge radius on circulation, are also found to be compatible with the minimum principle. Finally, the minimum principle is used to interpret the stability of infinitesimal and finite amplitude disturbances in an initially laminar, parallel shear flow, with results that are consistent with experiment and linearized hydrodynamic stability theory. These results suggest that a thermodynamic approach may be useful in unifying the understanding of many diverse phenomena in aerodynamics and aeroacoustics.

Horne, W. Clifton↗

DOE OSTI2018

Coupled Multimodal Dynamics of Hydrogen-Containing Ion Networks in Water-Deficient, Sodium Hydroxide-Aluminate Solutions

The (meta)stability of low water activity sodium hydroxide/aluminate (Na + OH–/Al(OH) 4 –) electrolytes dictates kinetics in the Bayer process for aluminum refining and high-level nuclear waste processing. We utilized quasi-elastic neutron scattering (QENS) and proton nuclear magnetic resonance spectroscopy ( 1 H NMR) in extremely concentrated sodium aluminate solutions to investigate the picosecond (ps) to microsecond (ms) timescale motions of H-bearing species (Al(OH) 4 – monomers/clusters, OH– and H 2 O). In the QENS data, in contrast to typical liquids, no short-time translational diffusion was observed at 293 K, but two types of localized motions were found: (i) local backbone tumbling or a formation of large hydrated ion clusters on the order of 40–60 ps; and (ii) much slower, complex, and collective dynamics of the ensemble of H-bearing species on the order of 350–750 ps. Variable temperature, pulsed field gradient, diffusion-ordered 1H NMR was used to determine the ensemble translational motion along with relaxometry to calculate rotational correlation coefficients. The ensemble rotational correlation times were on the order of 184–300 ps from 1 H NMR, which is consistent with the timescale of the QENS components. Complementary molecular dynamics simulation of NaOH solutions exhibit extensive ion networks potentially responsible for the observed dynamical coupling of water with the motion of large hydrated ion clusters. Understanding these collective motions will aid in predicting the behavior of complex solutions during aluminum production and during nuclear waste processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Nanoscale Electrochemical Phenomena of Polarization Switching in Ferroelectrics

Polarization switching is a fundamental feature of ferroelectric materials, enabling a plethora of applications and captivating the attention of the scientific community for over half a century. Many previous studies considered ferroelectric switching as a purely physical process, whereas polarization is fully controlled by the superposition of electric fields. However, screening charge is required for thermodynamic stability of the single domain state that is of interest in many technological applications. The screening process has always been assumed to be fast; thus, the rate-limiting phenomena were believed to be domain nucleation and domain wall dynamics. In this manuscript, we demonstrate that polarization switching under an atomic force microscopy tip leads to reversible ionic motion in the top 3 nm of PbZr 0.2 Ti 0.8 O 3 surface layer. Furthermore, this evidence points to a strong chemical component to a process believed to be purely physical and has major implications for understanding ferroelectric materials, making ferroelectric devices, and interpreting local ferroelectric switching.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Configurational stability of a crack propagating in a material with mode-dependent fracture energy – Part II: Drift of fracture facets in mixed-mode I+II+III

In earlier papers (Leblond et al., 2011; 2019), we presented linear stability analyses of the coplanar propagation of a crack loaded in mixed-mode I+III, based on a “double” propagation criterion combining (Griffith, 1920)’s energetic condition and (Goldstein and Salganik, 1974)’s principle of local symmetry. The difference between the two papers was that in the more recent one, the local value of the critical energy-release-rate was no longer considered as a constant, but heuristically allowed to depend upon the ratio of the local mode III to mode I stress intensity factors. This led to a much improved, qualitatively acceptable agreement of theory and experiments, for the “threshold” value of the ratio of the unperturbed mode III to mode I stress intensity factors, above which coplanar propagation becomes unstable. In this paper, the analysis is extended to the case where a small additional mode II loading component is present in the initially planar configuration of the crack, generating a small, general kink of this crack from the moment it is applied. The main new effect resulting from presence of such a loading component is that the instability modes present above the threshold must drift along the crack front during its propagation. This prediction may be useful for future theoretical interpretations of a number of experiments where such a drifting motion was indeed observed.

36 MATERIALS SCIENCE↗