Search NASA⌕ Search

SEARCH · Search NASA

Results for “Accelerated gradient momentum”

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 37 records · Page 2

All resistojet control of the NASA dual keel Space Station

This paper describes a control system that uses small (0.1-N) hydrogen-fueled, high-expansion-ratio resistojets to control a Space Station's position and attitude. The number and size of the resistojets is calculated based on requirements for stationkeeping, compensation for steady-state gravity-gradient and aerodynamic torques, momentum loading due to mobile remote manipulator motion, Shuttle Orbiter docking, and translational evasive maneuvering. The all-resistojet control system performs continuous orbit maintenance while keeping the linear accelerations on the spacecraft lower than they would be without any control system operating.

Paluszek, M. A.↗

Use of rotation to suppress thermosolutal convection in directionally solidified binary alloys

Effects of rotation on onset of convection during plane-front directional solidification of Pb-Sn and the pseudobinary system mercury cadmium telluride (Hg(1-x)Cd(x)Te), and on dendritic solidification of Pb-Sn have been studied by means of linear stability analysis. Incorporating Coriolis and centrifugal accelerations into the momentum equation of Coriell et al., we find that under realistic processing conditions, a large degree of stabilization can be achieved using modest rotation rates for both Pb-Sn and mercury cadmium telluride. At a growth velocity of 5 micron/sec and nominal liquid-side temperature gradient of 200 K/cm in Pb-Sn, rotation at 500 rpm results in a hundredfold increase in the critical Sn concentration. Large increases in the maximum allowable growth velocity at fixed melt composition are also attainable with modest rotation rates. The effect is amplified under conditions of reduced gravitational acceleration. For Hg(1-x)Cd(x)Te, we have also studied the nonrotating case. The key differences are due to the existence of a composition range for Hg(1-x)Cd(x)Te in which the melt density has a local maximum as a function of temperature. When the melt solidifies by cooling from below, the liquid density may initially increase with distance above the interface, before ultimately decreasing as the melt temperature increases above the value at which the local density maximum occurs. In contrast to the Pb-Sn case where density depends monotonically on temperature and composition, for Hg(1-x)Cd(x)Te there exists a critical value of the growth velocity above which plane-front solidification is unstable for all bulk CdTe mole fractions. Again, rotation leads to significant inhibition of onset. We identify the predicted stabilization with the Taylor-Proudman mechanism by which rotation inhibits thermal convection in a single-component fluid heated from below. In a binary liquid undergoing solidification, rotation inhibits the onset of buoyancy-driven convection, and has no effect on the short-wavelength morphological instability. At large growth velocities, the plane-front interface between liquid and solid becomes unstable with respect to a morphological instability and solidification occurs dendritically, with a mushy zone of dendrites and interdendritic fluid separating the solid from the melt. For the Pb-Sn system, rotation substantially suppresses the onset of convection in the mushy zone and in the overlying liquid, holding open the promise that rotation can suppress freckling and other macrosegregation defects.

Pearlstein, Arne J.↗

Direct Numerical Simulation of Smooth-Body Flow Separation Around a Ramp

Spanwise-periodic computation of a turbulent flow past a two-dimensional smooth ramp geometry is performed in the form of a direct numerical simulation. The Reynolds number based on the ramp height is about 147,000. A straight section that precedes the smooth ramp allows the incoming turbulent boundary layer to grow under a weak favorable pressure gradient. The boundary layer introduced at the domain inlet has a momentum-thickness based Reynolds number of 2000. The turbulent boundary layer nearing the ramp first interacts with a relatively stronger favorable pressure gradient, followed by a strong adverse pressure gradient. Consequently, the boundary layer experiences a modest acceleration before decelerating and separating. Analysis of the data over this region hints at the formation of an internal layer beneath the accelerated boundary layer. The analysis also reveals that this internal layer forms the origin of the free shear layer that emerges in the deceleration region and separates. The streamwise extent of the separated region is comparable to the ramp length, while the viscous layer thickness near reattachment is about the same as the ramp height; hence, the boundary layer undergoing separation and subsequent reattachment in the present configuration experiences its thickness being amplified by about tenfold. The reattached flow continues to develop further under a diminishing pressure gradient in the recovery region in a similar fashion to a zero pressure gradient turbulent boundary layer.

Turbulence↗

An experimental study of the properties of surface pressure fluctuations for separating turbulent boundary layers

Noise generated by helicopter and turbomachine rotors is a nuisance that designers would like to predict and to minimize within other design constraints. A key element for the noise calculation procedure is knowledge relating the flowfield structure to the surface pressure fluctuation structure. Surface pressure fluctuation data for zero-pressure-gradient and accelerating turbulent boundary layers were obtained. The zero-pressure-gradient turbulent boundary layers were examined with freestream velocities of 72 and 105 fps. Mean and fluctuation velocity profiles and streamwise velocity spectra and wavespeeds were obtained for momentum thickness Reynolds numbers up to 18000 for the zero-pressure-gradient case and up to 4000 for the favorable-pressure-gradient case. The wall shearing stress was estimated from a Clauser plot of the near wall data. It is clear that turbulent pressure fluctuations are produced by turbulent velocity fluctuations. Detailed simultaneous measurements of all of these fluctuations are needed to determine in more detail the structural relationships between velocity and pressure fields. Although some measurements were made for unseparated flows, none were made for separated flows.

Simpson, Roger L.↗

Space Station Control Requirements and Flywheel System Weights for Combined Momentum and Energy Storage

The specifications of the flywheel system for momentum storage and vehicle torquing are somewhat dependent upon the attitude control requirements of the space station in orbit. As a ground rule, the flywheel system will be sized large enough to provide all attitude maneuvers, if practical, to avoid or minimize turning on the reaction control system (RCS). The RCS, whenever used, expels expensive mass and tends to contaminate optical surfaces of the vehicle. The vehicle rate and acceleration specifications of 0.10 deg/sec and 0.01 deg/square sec are tentative, and may be reduced if lesser values are more practical for flywheel design. For local vertical attitude hold, the average attitude error should be zero, and not the classical 1 degree, since control moment gyro (CMG) gimbal angles provide an exact reference feedback for gravity gradient momentum. Docking presents a problem for docking transients and attitude alignment which will require use of the RCS.

Elam, F. M.↗

Simulating the Effect of Space Vehicle Environments on Directional Solidification of a Binary Alloy

Space microgravity missions are designed to provide a microgravity environment for scientific experiments, but these missions cannot provide a perfect environment, due to vibrations caused by crew activity, on-board experiments, support systems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps. Therefore, it is necessary to predict the impact of these vibrations on space experiments, prior to performing them. Simulations were conducted to study the effect of the vibrations on the directional solidification of a dendritic alloy. Finite element ca!cu!attie?ls were dme with a simd2titcr based on a continuum model of dendritic solidification, using the Fractional Step Method (FSM). The FSM splits the solution of the momentum equation into two steps: the viscous intermediate step, which does not enforce continuity; and the inviscid projection step, which calculates the pressure and enforces continuity. The FSM provides significant computational benefits for predicting flows in a directionally solidified alloy, compared to other methods presently employed, because of the efficiency gains in the uncoupled solution of velocity and pressure. finite differences, arises when the interdendritic liquid reaches the eutectic temperature and concentration. When a node reaches eutectic temperature, it is assumed that the solidification of the eutectic liquid continues at constant temperature until all the eutectic is solidified. With this approach, solidification is not achieved continuously across an element; rather, the element is not considered solidified until the eutectic isotherm overtakes the top nodes. For microgravity simulations, where the convection is driven by shrinkage, it introduces large variations in the fluid velocity. When the eutectic isotherm reaches a node, all the eutectic must be solidified in a short period, causing an abrupt increase in velocity. To overcome this difficulty, we employed a scheme to numerically predict a more accurate value for the rate of change of fraction of liquid as the liquid in an element solidifies. The new method enables us to contrast results of simulations in which the alloy is subjected to no gravity or a steady-state acceleration versus simulations when the alloy is subjected to vibration disturbances; therefore, the effect of vibration disturbances can be assessed more accurately. To assess the impact of these vibration-perturbations, transient accelerometer data from a space shuttle mission are used as inputs for the simulation model. These on-orbit acceleration data were obtained from the Microgravity Science Division at Glenn Research Center (GRC- MSD) and are applied to the buoyancy term of the momentum equation in a simulation of a Pb-5.8 wt. % Sb alloy that solidifies in a thermal gradient of 4000 K/m and a translation velocity of 3 p d s . Figure 2 shows the vertical velocity of a node that begins in the all-liquid region and subsequently solidifies; the vibrations are applied at 5000 seconds in this simulation. An important difficulty, common to all solidification models based on finite elements or 2 The magnitudes of the velocity oscillations that are vibration-induced are very small and acceptable. The biggest concern is whether the concentration of the liquid near the dendrite tips is distorted because of the vibration-induced perturbations. Results for this case show no concentration oscillations present in the all-liquid region.

Westra, D. G.↗

Accelerating Continuous Variable Coherent Ising Machines Via Momentum

The Coherent Ising Machine (CIM) is a non-conventional architecture that takes inspiration from physical annealing processes to solve Ising problems heuristically. Its dynamics are naturally continuous and described by a set of ordinary differential equations that have been proven to be useful for the optimization of continuous variables nonconvex quadratic optimization problems. The dynamics of such Continuous Variable CIMs (CV-CIM) encourage optimization via optical pulses whose amplitudes are determined by the negative gradient of the objective; however, standard gradient descent is known to be trapped by local minima and hampered by poor problem conditioning. In this work, we propose to modify the CV-CIM dynamics using more sophisticated pulse injections based on tried-and-true optimization techniques such as momentum and Adam. Through numerical experiments, we show that the momentum and Adam updates can significantly speed up the CV-CIM’s convergence and improve sample diversity over the original CV-CIM dynamics. We also find that the Adam-CV-CIM’s performance is more stable as a function of feedback strength, especially on poorly conditioned instances, resulting in an algorithm that is more robust, reliable, and easily tunable. More broadly, we identify the CIM dynamical framework as a fertile opportunity for exploring the intersection of classical optimization and modern analog computing.

Ising Model↗

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗

Doppler radar observations of the three-dimensional turbulent structure of a quasi-steady thunderstorm

The spatial patterns and temporal evolution of Doppler radar-inferred turbulence within a quasi-steady, left-moving thunderstorm are reported. An organized downdraft circulation was observed, which maintained a persistent, continuously-propagating, anticyclonic and turbulent updraft. The updrafts attained peak speeds of 25 m/s, and maximum estimated reflectivity factors exhibited peak magnitudes of 50-55 dBZ. Mechanical and buoyant production appeared to be the primary generators of turbulence, although buoyancy generation may have also significantly contributed, since turbulent and dynamic entrainment into the vertically accelerating updraft was appreciable below the strongest turbulent activity at 7 km. Patterns of the radial shears of radial velocity exhibited a significant decrease in strength and areal coverage as the storm weakened, which is attributed to decreased activity of large turbulent eddies and the diminution of large-scale mean wind gradients produced by vertical exchanges of horizontal momentum and strong updraft gradients.

Knupp, K. R.↗

Ocean-Atmosphere Interaction Over Agulhas Extension Meanders

Many years of high-resolution measurements by a number of space-based sensors and from Lagrangian drifters became available recently and are used to examine the persistent atmospheric imprints of the semi-permanent meanders of the Agulhas Extension Current (AEC), where strong surface current and temperature gradients are found. The sea surface temperature (SST) measured by the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and the chlorophyll concentration measured by the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) support the identification of the meanders and related ocean circulation by the drifters. The collocation of high and low magnitudes of equivalent neutral wind (ENW) measured by Quick Scatterometer (QuikSCAT), which is uniquely related to surface stress by definition, illustrates not only the stability dependence of turbulent mixing but also the unique stress measuring capability of the scatterometer. The observed rotation of ENW in opposition to the rotation of the surface current clearly demonstrates that the scatterometer measures stress rather than winds. The clear differences between the distributions of wind and stress and the possible inadequacy of turbulent parameterization affirm the need of surface stress vector measurements, which were not available before the scatterometers. The opposite sign of the stress vorticity to current vorticity implies that the atmosphere spins down the current rotation through momentum transport. Coincident high SST and ENW over the southern extension of the meander enhance evaporation and latent heat flux, which cools the ocean. The atmosphere is found to provide negative feedback to ocean current and temperature gradients. Distribution of ENW convergence implies ascending motion on the downwind side of local SST maxima and descending air on the upwind side and acceleration of surface wind stress over warm water (deceleration over cool water); the convection may escalate the contrast of ENW over warm and cool water set up by the dependence of turbulent mixing on stability; this relation exerts a positive feedback to the ENW-SST relation. The temperature sounding measured by the Atmospheric Infrared Sounder(AIRS) is consistent with the spatial coherence between the cloud-top temperature provided by the International Satellite Cloud Climatology Project (ISCCP) and SST. Thus ocean mesoscale SST anomalies associated with the persistent meanders may have a long-term effect well above the midlatitude atmospheric boundary layer, an observation not addressed in the past.

Atmospheric Infrared Sounder (AIRS)↗

Anti-solar acceleration of ionospheric plasma across the Venus terminator

It is demonstrated that the horizontal ionospheric particle pressure gradient across the Venus terminator is the principal body force accelerating the plasm to the observed anti-solar velocity. The large scale horizontal electromagnetic body force is typically an order of magnitude smaller than the particle pressure gradient. The viscous body forces above 300 km - drag or acceleration - are also an order of magnitude smaller than the pressure gradient body force. In the immediate vicinity of the ionopause where the ionospheric plasma density decreases below 1000 ions/cu cm and the magnetic field strength increases, the electromagnetic body force may become significant. The electromagnetic body force may also be significant in the nightside ionosphere. The bulk of the ionospheric flow momentum is not derived from the ionosheath momentum.

Knudsen, W. C.↗

Analysis of the low gravity tolerance of Bridgman-Stockbarger crystal growth. I - Steady and impulse accelerations

The effects of steady and impulse-type residual accelerations on dopant distributions during directional solidification in 2D and 3D 'generic' models of the Bridgman-Stockbarger technique are investigated using numerical methods. The calculations are based on the thermophysical properties of molten germanium doped with a low concentration of gallium. A Chebyshev collocation pseudospectral method is used for the solution of the governing momentum-, mass-, species-, and heat-transfer equations. Only convection caused by temperature gradients is considered. It is found that lateral nonuniformity in composition is very sensitive to the orientation of the steady component of the residual gravity vector and to the particular operating conditions under consideration. It is also found that laterally or radially averaged composition profiles are alone insufficient to describe the extent of residual convection in a spacecraft environment. The effects of impulse-type disturbances can be severe and can extend for times on the order of 1000 sec after the termination of the impulse.

Alexander, J. Iwan D.↗

Accuracy of the explicit energy-conserving particle-in-cell method for under-resolved simulations of capacitively coupled plasma discharges

The traditional explicit electrostatic momentum-conserving particle-in-cell algorithm requires strict resolution of the electron Debye length to deliver numerical stability and accuracy. The explicit electrostatic energy-conserving particle-in-cell algorithm alleviates this constraint with minimal modification to the traditional algorithm, retaining its simplicity, ease of parallelization, and acceleration on modern supercomputing architectures. In this article, we apply the algorithm to model a one-dimensional radio frequency capacitively coupled plasma discharge relevant to industrial applications. The energy-conserving approach closely matches the results from the momentum-conserving algorithm and retains accuracy even for cell sizes up to 8 times the electron Debye length. For even larger cells, the algorithm loses accuracy due to poor resolution of steep gradients within the radio frequency sheath. Accuracy can be recovered by adopting a non-uniform grid, which resolves the sheath and allows for cell sizes up to 32 times the electron Debye length in the quasi-neutral bulk of the discharge. The effect is an up to 8 times reduction in the number of required simulation cells, an improvement that can compound in higher-dimensional simulations. We therefore consider the explicit energy-conserving algorithm as a promising approach to significantly reduce the computational cost of full-scale device simulations and a pathway to delivering kinetic simulation capabilities of use to industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Scalar flux transport models for self-similar turbulent mixing

A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.

Reynolds-averaged Navier Stokes↗

Advances in Non-contact Measurement of Creep Properties

Our team has developed a novel approach to measuring creep at extremely high temperatures using electrostatic levitation (ESL). This method has been demonstrated on niobium up to 2300 C, while ESL has melted tungsten (3400 C). High-precision machined spheres of the sample are levitated in the NASA MSFC ESL, a national user facility, and heated with a laser. The laser is aligned off-center so that the absorbed photons transfer their momentum to the sample, causing it to rotate at up to 250,000+ RPM. The rapid rotation loads the sample through centripetal acceleration, causing it to deform. The deformation of the sample is captured on high-speed video, which is analyzed by machine-vision software from the University of Massachusetts. The deformations are compared to finite element models to determine the constitutive constants in the creep relation. Furthermore, the noncontact method exploits stress gradients within the sample to determine the stress exponent in a single test. This method was validated in collaboration with the University of Tennessee for niobium at 1985 C, with agreement within the uncertainty of the conventional measurements. A similar method is being employed on Ultra-High-Temperature ZrB2- SiC composites, which may see application in rocket nozzles and sharp leading edges for hypersonic vehicles.

Hyers, Robert W.↗

Human Outer Solar System Exploration via Q-Thruster Technology

Propulsion technology development efforts at the NASA Johnson Space Center continue to advance the understanding of the quantum vacuum plasma thruster (QThruster), a form of electric propulsion. Through the use of electric and magnetic fields, a Q-thruster pushes quantum particles (electrons/positrons) in one direction, while the Qthruster recoils to conserve momentum. This principle is similar to how a submarine uses its propeller to push water in one direction, while the submarine recoils to conserve momentum. Based on laboratory results, it appears that continuous specific thrust levels of 0.4 - 4.0 N/kWe are achievable with essentially no onboard propellant consumption. To evaluate the potential of this technology, a mission analysis tool was developed utilizing the Generalized Reduced Gradient non-linear parameter optimization engine contained in the Microsoft Excel® platform. This tool allowed very rapid assessments of "Q-Ship" minimum time transfers from earth to the outer planets and back utilizing parametric variations in thrust acceleration while enforcing constraints on planetary phase angles and minimum heliocentric distances. A conservative Q-Thruster specific thrust assumption (0.4 N/kWe) combined with "moderate" levels of space nuclear power (1 - 2 MWe) and vehicle specific mass (45 - 55 kg/kWe) results in continuous milli-g thrust acceleration, opening up realms of human spaceflight performance completely unattainable by any current systems or near-term proposed technologies. Minimum flight times to Mars are predicted to be as low as 75 days, but perhaps more importantly new "retro-phase" and "gravity-augmented" trajectory shaping techniques were revealed which overcome adverse planetary phasing and allow virtually unrestricted departure and return opportunities. Even more impressively, the Jovian and Saturnian systems would be opened up to human exploration with round-trip times of 21 and 32 months respectively including 6 to 12 months of exploration at the destinations. Finally, interstellar trip times are assessed at milli-g acceleration levels.

Joosten, B. Kent↗

Comprehensive Validation of an Intermittency Transport Model for Transitional Low-Pressure Turbine Flows

A transport equation for the intermittency factor is employed to predict transitional flows under the effects of pressure gradients, freestream turbulence intensities, Reynolds number variations, flow separation and reattachment. and unsteady wake-blade interactions representing diverse operating conditions encountered in low-pressure turbines. The intermittent behaviour of the transitional flows is taken into account and incorporated into computations by modifying the eddy viscosity, Mu(sub t), with the intermittency factor, gamma. Turbulent quantities are predicted by using Menter's two-equation turbulence model (SST). The onset location of transition is obtained from correlations based on boundary-layer momentum thickness, acceleration parameter, and turbulence intensity. The intermittency factor is obtained from a transport model which can produce both the experimentally observed streamwise variation of intermittency and a realistic profile in the cross stream direction. The intermittency transport model is tested and validated against several well documented low pressure turbine experiments ranging from flat plate cases to unsteady wake-blade interaction experiments. Overall, good agreement between the experimental data and computational results is obtained illustrating the predicting capabilities of the model and the current intermittency transport modelling approach for transitional flow simulations.

Suzen, Y. B.↗

The effects of laser polarization and wavelength on injection dynamics of a laser wakefield accelerator

Here, we investigate the effects of laser polarization and wavelength on electron injection dynamics in a laser wakefield accelerator. During the ionization process, electrons gain residual momentum and kinetic energy via above threshold ionization, which has a strong dependence on laser polarization. A circularly polarized laser pulse results in a much higher residual momentum and kinetic energy gain for the ionized electrons compared with the linearly polarized case. This residual momentum results in particle injection because of the sensitivity of particle trapping to the initial conditions and enhanced the total injected beam charge in both experiments and particle-in-cell simulations. Due to the strong correlation of above threshold ionization with laser wavelength, in this work we extended the investigation to long wavelength (up to 20 μm) drive pulses using particle-in-cell simulations. Owing to the gain in kinetic energy, it may be expected that the charge trapped would consistently increase for circular polarization with increasing laser wavelength, but this was not observed. Instead, there are oscillations with wavelength in the relative trapped charge between linear and circular polarization cases, which arise because of ionization and heating effects on the plasma. Our studies highlight the complex interplay between several different physical effects, including injection regimes—above threshold ionization assisted injection, wave-breaking injection by carrier-envelope-phase effects and ionization injection—ionization gradient induced laser pulse evolution, and thermal modifications to the wake structure that need considering when extrapolating laser wakefield acceleration to different wavelength regimes.

Ma, Y. (ORCID:0000000239194856)↗