Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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 415 records · Page 23

Mission and Vehicle-Level Updates for the Parallel Electric-Gas Architecture with Synergistic Utilization Scheme (PEGASUS) Concept Aircraft

NASA created the PEGASUS concept with the goal of lowering mission energy (a surrogate for operating cost) compared to other regional aircraft by leveraging electrified aircraft propulsion (EAP). Since its inception, researchers have explored multiple facets of PEGASUS in varying fidelity but have not completed a rigorous, integrated design. The goal of this memorandum is to provide an updated design using recent studies and improved methods. This memorandum explores the initial vehicle concept and concept of operations, while considering ways to improve both the mission concept of operations and the integrated vehicle-level performance. Additionally, the design and analysis methodologies for EAP-enabled aircraft concepts are improved in several areas. This research incorporates new propulsion-airframe integration and wing weight surrogates to model the impacts of wingtip propulsors on the configuration. Detailed weight and balance calculations enable calculating dynamic stability and flight qualities within the conceptual design environment. Ultimately, the vehicle is optimized to reduce "well-to-wake" equivalent CO2, CO2e, rather than minimizing either fuel (or total energy) consumption or maximum takeoff weight. Using fuel/energy or takeoff weight leads to conflicting optima for hybrid-electric aircraft. Two aircraft are developed to provide points of comparison for PEGASUS: an advanced conventional turboprop vehicle and a hybrid-electric variant. The results show that the PEGASUS concept can reduce CO2e relative to the advanced turboprop or a hybrid-electric propulsion architecture, albeit with an increase in maximum takeoff weight. PEGASUS's maximum takeoff weight is 55% heavier than the advanced conventional turboprop but releases 18% less CO2e for a 400 nmi mission. Over the same mission, PEGASUS's maximum takeoff weight is 47% heavier than the comparator hybrid-electric vehicle but releases 12% less CO2e. This study shows that the PEGASUS configuration reduces CO2e through its use of wingtip propulsors and that its benefit is not solely a result of switching to a hybrid-electric propulsion architecture. PEGASUS achieves this reduction in CO2e while maintaining satisfactory Level 1 or 2 flight qualities for all of its longitudinal- and lateral-directional modes.

Hybrid-electric↗

Rotary Balance Wind Tunnel Testing for the FASER Flight Research Aircraft

Flight dynamics research was conducted to collect and analyze rotary balance wind tunnel test data in order to improve the aerodynamic simulation and modeling of a low-cost small unmanned aircraft called FASER (Free-flying Aircraft for Sub-scale Experimental Research). The impetus for using FASER was to provide risk and cost reduction for flight testing of more expensive aircraft and assist in the improvement of wind tunnel and flight test techniques, and control laws. The FASER research aircraft has the benefit of allowing wind tunnel and flight tests to be conducted on the same model, improving correlation between wind tunnel, flight, and simulation data. Prior wind tunnel tests include a static force and moment test, including power effects, and a roll and yaw damping forced oscillation test. Rotary balance testing allows for the calculation of aircraft rotary derivatives and the prediction of steady-state spins. The rotary balance wind tunnel test was conducted in the NASA Langley Research Center (LaRC) 20-Foot Vertical Spin Tunnel (VST). Rotary balance testing includes runs for a set of given angular rotation rates at a range of angles of attack and sideslip angles in order to fully characterize the aircraft rotary dynamics. Tests were performed at angles of attack from 0 to 50 degrees, sideslip angles of -5 to 10 degrees, and non-dimensional spin rates from -0.5 to 0.5. The effects of pro-spin elevator and rudder deflection and pro- and anti-spin elevator, rudder, and aileron deflection were examined. The data are presented to illustrate the functional dependence of the forces and moments on angle of attack, sideslip angle, and angular rate for the rotary contributions to the forces and moments. Further investigation is necessary to fully characterize the control effectors. The data were also used with a steady state spin prediction tool that did not predict an equilibrium spin mode.

Denham, Casey↗

Macroscopic balance model for wave rotors

A mathematical model for multi-port wave rotors is described. The wave processes that effect energy exchange within the rotor passage are modeled using one-dimensional gas dynamics. Macroscopic mass and energy balances relate volume-averaged thermodynamic properties in the rotor passage control volume to the mass, momentum, and energy fluxes at the ports. Loss models account for entropy production in boundary layers and in separating flows caused by blade-blockage, incidence, and gradual opening and closing of rotor passages. The mathematical model provides a basis for predicting design-point wave rotor performance, port timing, and machine size. Model predictions are evaluated through comparisons with CFD calculations and three-port wave rotor experimental data. A four-port wave rotor design example is provided to demonstrate model applicability. The modeling approach is amenable to wave rotor optimization studies and rapid assessment of the trade-offs associated with integrating wave rotors into gas turbine engine systems.

Welch, Gerard E.↗

Aeroelastic Computations of a Compressor Stage Using the Harmonic Balance Method

The aeroelastic characteristics of a compressor stage were analyzed using a computational fluid dynamic (CFD) solver that uses the harmonic balance method to solve the governing equations. The three dimensional solver models the unsteady flow field due to blade vibration using the Reynolds-Averaged Navier-Stokes equations. The formulation enables the study of the effect of blade row interaction through the inclusion of coupling modes between blade rows. It also enables the study of nonlinear effects of high amplitude blade vibration by the inclusion of higher harmonics of the fundamental blade vibration frequency. In the present work, the solver is applied to study in detail the aeroelastic characteristics of a transonic compressor stage. Various parameters were included in the study: number of coupling modes, blade row axial spacing, and operating speeds. Only the first vibration mode is considered with amplitude of oscillation in the linear range. Both aeroelastic stability (flutter) of rotor blade and unsteady loading on the stator are calculated. The study showed that for the stage considered, the rotor aerodynamic damping is not influenced by the presence of the stator even when the axial spacing is reduced by nearly 25 percent. However, the study showed that blade row interaction effects become important for the unsteady loading on the stator when the axial spacing is reduced by the same amount.

Reddy, T. S. R.↗

Parametric Identification of Nonlinear Dynamical Systems

In this project, we looked at the application of harmonic balancing as a tool for identifying parameters (HBID) in a nonlinear dynamical systems with chaotic responses. The main idea is to balance the harmonics of periodic orbits extracted from measurements of each coordinate during a chaotic response. The periodic orbits are taken to be approximate solutions to the differential equations that model the system, the form of the differential equations being known, but with unknown parameters to be identified. Below we summarize the main points addressed in this work. The details of the work are attached as drafts of papers, and a thesis, in the appendix. Our study involved the following three parts: (1) Application of the harmonic balance to a simulation case in which the differential equation model has known form for its nonlinear terms, in contrast to a differential equation model which has either power series or interpolating functions to represent the nonlinear terms. We chose a pendulum, which has sinusoidal nonlinearities; (2) Application of the harmonic balance to an experimental system with known nonlinear forms. We chose a double pendulum, for which chaotic response were easily generated. Thus we confronted a two-degree-of-freedom system, which brought forth challenging issues; (3) A study of alternative reconstruction methods. The reconstruction of the phase space is necessary for the extraction of periodic orbits from the chaotic responses, which is needed in this work. Also, characterization of a nonlinear system is done in the reconstructed phase space. Such characterizations are needed to compare models with experiments. Finally, some nonlinear prediction methods can be applied in the reconstructed phase space. We developed two reconstruction methods that may be considered if the common method (method of delays) is not applicable.

Feeny, Brian↗

Exploration of Defect Dynamics and Color Center Qubit Synthesis with Pulsed Ion Beams

Short-pulse ion beams have been developed in recent years and now enable applications in materials science. A tunable flux of selected ions delivered in pulses of a few nanoseconds can affect the balance of defect formation and dynamic annealing in materials. We report results from color center formation in silicon with pulses of 900 keV protons. G-centers in silicon are near-infrared photon emitters with emerging applications as single-photon sources and for spin-photon qubit integration. G-centers consist of a pair of substitutional carbon atoms and one silicon interstitial atom and are often formed by carbon ion implantation and thermal annealing. Here, we report on G-center formation with proton pulses in silicon samples that already contained carbon, without carbon ion implantation or thermal annealing. The number of G-centers formed per proton increased when we increased the pulse intensity from 6.9 × 10 9 to 7.9 × 10 10 protons/cm 2 /pulse, demonstrating a flux effect on G-center formation efficiency. We observe a G-center ensemble linewidth of 0.1 nm (full width half maximum), narrower than previously reported. Pulsed ion beams can extend the parameter range available for fundamental studies of radiation-induced defects and the formation of color centers for spin-photon qubit applications.

43 PARTICLE ACCELERATORS↗

Reinforcement Learning for Load-balanced Parallel Particle Tracing

We explore an online reinforcement learning (RL) paradigm to dynamically optimize parallel particle tracing performance in distributed-memory systems. Our method combines three novel components: (1) a work donation algorithm, (2) a high-order workload estimation model, and (3) a communication cost model. First, we design an RL-based work donation algorithm. Our algorithm monitors workloads of processes and creates RL agents to donate data blocks and particles from high-workload processes to low-workload processes to minimize program execution time. The agents learn the donation strategy on the fly based on reward and cost functions designed to consider processes' workload changes and data transfer costs of donation actions. Second, we propose a workload estimation model, helping RL agents estimate the workload distribution of processes in future computations. Third, we design a communication cost model that considers both block and particle data exchange costs, helping RL agents make effective decisions with minimized communication costs. We demonstrate that our algorithm adapts to different flow behaviors in large-scale fluid dynamics, ocean, and weather simulation data. Our algorithm improves parallel particle tracing performance in terms of parallel efficiency, load balance, and costs of I/O and communication for evaluations with up to 16,384 processors.

Distributed and parallel particle tracing↗

Full-Span Tiltrotor Aeroacoustic Model (TRAM) Overview and 40- by 80-Foot Wind Tunnel Test

Most helicopter data trends cannot be extrapolated to tiltrotors because blade geometry and aerodynamic behavior, as well as rotor and fuselage interactions, are significantly different for tiltrotors. A tiltrotor model has been developed to investigate the aeromechanics of tiltrotors, to develop a comprehensive database for validating tiltrotor analyses, and to provide a research platform for supporting future tiltrotor designs. The Full-Span Tiltrotor Aeroacoustic Model (FS TRAM) is a dual-rotor, powered aircraft model with extensive instrumentation for measurement of structural and aerodynamic loads. This paper will present the Full-Span TRAM test capabilities and the first set of data obtained during a 40- by 80-Foot Wind Tunnel test conducted in late 2000 at NASA Ames Research Center. The Full-Span TRAM is a quarter-scale representation of the V-22 Osprey aircraft, and a heavily instrumented NASA and U.S. Army wind tunnel test stand. Rotor structural loads are monitored and recorded for safety-of-flight and for information on blade loads and dynamics. Left and right rotor balance and fuselage balance loads are monitored for safety-of-flight and for measurement of vehicle and rotor aerodynamic performance. Static pressure taps on the left wing are used to determine rotor/wing interactional effects and rotor blade dynamic pressures measure blade airloads. All of these measurement capabilities make the FS TRAM test stand a unique and valuable asset for validation of computational codes and to aid in future tiltrotor designs. The Full-Span TRAM was tested in the NASA Ames Research Center 40- by 80-Foot Wind Tunnel from October through December 2000. Rotor and vehicle performance measurements were acquired in addition to wing pressures, rotor acoustics, and Laser Light Sheet (LLS) flow visualization data. Hover, forward flight, and airframe (rotors off) aerodynamic runs were performed. Helicopter-mode data were acquired during angle of attack and thrust sweeps for a variety of tunnel speeds. Wake geometry images were acquired using LLS photographs and suggest dual tip vortex formation at low thrust conditions. The full paper will include comparisons to isolated-rotor TRAM data acquired at the Duits-Nederlandse Windtunnel (DNW) in 1998. The FS TRAM has been established as a valuable national asset for tiltrotor research. Data reduction and analysis of the 40- by 80-Foot Wind Tunnel test results are underway. Follow-on testing of the FS TRAM is currently being planned for the NASA Ames 80- by 120-Foot Wind Tunnel in late 2001.

McCluer, Megan S.↗

Increasing the Dynamic Pressure Capability of the NASA Langley/ODU 6-inch MSBS

The 6-inch NASA/ODU Magnetic Suspension and Balance System (MSBS) has been configured for dynamic stability testing of blunt-body atmospheric entry capsules. Tests have been successfully accomplished in the low-speed, open-circuit wind tunnel, at speeds up to around 40 m/s. The wind tunnel is designed to reach around 150 m/s, resulting in dynamic pressures comparable to those projected to arise in a future supersonic MSBS facility. Extensive system upgrades are being undertaken to permit testing at higher speeds/dynamic pressures, including control system enhancements, revised position and attitude sensing, and activation of additional electromagnets in the existing array. This paper will review recent progress in all these areas.

Mark Schoenenberger↗

Systematic Differences between Satellite-Based Presipitation Climatologies over the Tropical Oceans

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. It is clear however that there still remain significant uncertainties with satellite precipitation retrievals which limit their usefulness for many purposes. Systematic differences i'A tropical precipitation estimates have been brought to light in comparison activities such as the GPCP Algorithm Intercomparison Project and more recent Wetnet Precipitation Intercomparison Project 3. These uncertainties are assuming more importance because of the demands for validation associated with global climate modeling and data assimilation methodologies. The objective of the present study is to determine the physical basis for systematic differences in spatial structure of tropical precipitation as portrayed by several different satellite-based data sets. The study is limited to oceanic regions only and deals primarily with aspects of spatial variability. We are specifically interested in why MSU channel 1 and GPI precipitation differences are so striking over the Eastern Pacific ITCZ and why they both differ from other microwave emission-based precipitation estimates from SSM/I and a scattering-based deep convective ice index from MSU channel 2. Our results to date have shown that MSU channel I precipitation estimates are biased high over the Eastern Pacific ITCZ because of two factors: (1) the hypersensitivity of this frequency to cloud water in contrast to falling rain drops, and (2) unaccounted for scattering effects by precipitation-size ice which depresses the signal of the liquid water emission. Likewise, cold cloud top climatologies such as the GPI show an excess (a deficit) in estimated rainfall over the E. Pacific ITCZ (Warm Pool region). We show that these algorithms need to account for regionally varying heights (or temperatures) at which tropical convection detrains to form cirrus shields. A second objective we pursue is to identify variations in the macroscale cloud physical and thermodynamic properties of precipitation regimes" and relate these differences to tropical dynamical mechanisms of tropical heat and moisture balance. Finally, we interpret the algorithm differences and their associations with tropical dynamics in terms of WCRP GPCP goals for constructing precipitation climatologies.

Robertson, Frankin R.↗

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗

TPSAS-NF1676L-15194-DND

Technical Motivation - 1950s: Indication of large lateral (roll/yaw) acceleration (“β-dot”) dynamic derivatives for configurations with strong vortex flows in the stall region (Lichtenstein, Campbell, Riley)! - 1960s: Increased emphasis on use of long-range missiles tempered by aircraft losses in close-in air combat - 1970s: - Effect of lateral acceleration derivatives on aircraft stability (Nguyen) - Characteristics and air combat advantages of rapid controlled large-amplitude maneuvering at low speeds near maximum lift shown in piloted simulations (Herbst) - 1980s: - Dynamic (forced-oscillation, rotary balance, curved- and rolling-flow) wind tunnel tests (Grafton, Lutze, Hultberg) - Incorporation of test data into computer simulations (Kalviste) - Technologies developed for controlled maneuvers at stall, post-stall angles of attack - Concern that conventional math modeling inadequately represented the dynamic component of these motions - Need to evaluate dependency of control laws on accurate computer simulations and motion prediction

Marilyn E Ogburn↗

Recent NASA Research on Aerodynamic Modeling of Post-Stall and Spin Dynamics of Large Transport Airplanes

A simulation study was conducted to investigate aerodynamic modeling methods for prediction of post-stall flight dynamics of large transport airplanes. The research approach involved integrating dynamic wind tunnel data from rotary balance and forced oscillation testing with static wind tunnel data to predict aerodynamic forces and moments during highly dynamic departure and spin motions. Several state-of-the-art aerodynamic modeling methods were evaluated and predicted flight dynamics using these various approaches were compared. Results showed the different modeling methods had varying effects on the predicted flight dynamics and the differences were most significant during uncoordinated maneuvers. Preliminary wind tunnel validation data indicated the potential of the various methods for predicting steady spin motions.

Murch, Austin M.↗

H ∞ Control for Energy Dispatch in Autonomous Nanogrid With Communication Delays

This paper proposes an optimal controller and estimator for energy dispatch to balance the power supply and demand considering communication delays. The proposed algorithm involves modeling an autonomous nanogrid (ANG) consisting of distributed energy resources, energy storage systems, loads, an $H$ ∞ controller with a reference power modulation technique, and a state estimator. The ANG was developed to express the dynamic supply-demand energy balance of a nanogird system. Reference power modulation was designed to generate the desired ESS power based on the imbalanced energy. Random communication delays were modeled using a stochastic variable satisfying the Bernoulli random binary distribution. The optimal $H$ ∞ controller and estimator were developed using a linear matrix inequality approach to exponentially stabilize the closed-loop system. Simulations were performed using real daily demand forecasts obtained from the Korea Meteorological Administration to demonstrate the effectiveness of the proposed real-time optimization algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Cascading toward a Kilometer‐scale GCM: Impacts of a Scale‐aware Convection Parameterization in the Goddard Earth Observing System GCM.

The NASA Goddard Earth Observing System global model is evaluated through a cascade of simulations with increasing horizontal resolution. This model employs a non-hydrostatic dynamical core and includes a scale-aware, deep convection parameterization (DPCP). 40-daysimulations at six resolutions (100km to 3km) with unvarying model formulation were produced. At the highest resolution, extreme experiments were carried out; one with no DPCP, and one with its scale-awareness eliminated. Simulated precipitation, radiative balance, and atmospheric thermodynamic and dynamical variables are well reproduced with respect to both observational and reanalysis data. As model resolution increases, the convective precipitation smoothly transitions from being mostly produced by the convection parameterization to the cloud microphysics parameterization. However, contrary to current thought, these extreme cases argue for maintaining, to some extent, the scale-aware deep convection parameterization even at 3kmscale, as the run relying solely on explicit grid-scale production of rainfall performs more poorly at this resolution.

Saulo Freitas↗

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗

Greenland Glacier Albedo Variability

The program for Arctic Regional Climate Assessment (PARCA) is a NASA-funded project with the prime goal of addressing the mass balance of the Greenland ice sheet. Since the formal initiation of the program in 1995, there has been a significant improvement in the estimates of the mass balance of the ice sheet. Results from this program reveal that the high-elevation regions of the ice sheet are approximately in balance, but the margins are thinning. Laser surveys reveal significant thinning along 70 percent of the ice sheet periphery below 2000 m elevations, and in at least one outlet glacier, Kangerdlugssuaq in southeast Greenland, thinning has been as much as 10 m/yr. This study examines the albedo variability in four outlet glaciers to help separate out the relative contributions of surface melting versus ice dynamics to the recent mass balance changes. Analysis of AVHRR Polar Pathfinder albedo shows that at the Petermann and Jakobshavn glaciers, there has been a negative trend in albedo at the glacier terminus from 1981 to 2000, whereas the Stor+strommen and Kangerdlugssuaq glaciers show slightly positive trends in albedo. These findings are consistent with recent observations of melt extent from passive microwave data which show more melt on the western side of Greenland and slightly less on the eastern side. Significance of albedo trends will depend on where and when the albedo changes occur. Since the majority of surface melt occurs in the shallow sloping western margin of the ice sheet where the shortwave radiation dominates the energy balance in summer (e.g. Jakobshavn region) this region will be more sensitive to changes in albedo than in regions where this is not the case. Near the Jakobshavn glacier, even larger changes in albedo have been observed, with decreases as much as 20 percent per decade.

Source record↗

Lee waves: Benign and malignant

The flow of an incompressible fluid over an obstacle will produce an oscillation in which buoyancy is the restoring force, called a gravity wave. For disturbances of this scale, the atmosphere may be treated as dynamically incompressible, even though there exists a mean static upward density gradient. Even in the linear approximation - i.e., for small disturbances - this model explains a great many of the flow phenomena observed in the lee of mountains. However, nonlinearities do arise importantly, in three ways: (1) through amplification due to the decrease of mean density with height; (2) through the large (scaled) size of the obstacle, such as a mountain range; and (3) from dynamically singular levels in the fluid field. These effects produce a complicated array of phenomena - large departure of the streamlines from their equilibrium levels, high winds, generation of small scales, turbulence, etc. - that present hazards to aircraft and to lee surface areas. The nonlinear disturbances also interact with the larger-scale flow in such a manner as to impact global weather forecasts and the climatological momentum balance. If there is no dynamic barrier, these waves can penetrate vertically into the middle atmosphere (30-100 km), where recent observations show them to be of a length scale that must involve the coriolis force in any modeling. At these altitudes, the amplitude of the waves is very large, and the phenomena associated with these wave dynamics are being studied with a view to their potential impact on high performance aircraft, including the projected National Aerospace Plane (NASP). The presentation shows the results of analysis and of state-of-the-art numerical simulations, validated where possible by observational data, and illustrated with photographs from nature.

Wurtele, M. G.↗