Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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 469 records · Page 26

A Mixed-Fidelity Approach for Design of Low-Boom Supersonic Aircraft

This paper documents a mixed-fidelity approach for the design of low-boom supersonic aircraft as a viable approach for designing a practical low-boom supersonic configuration. A low-boom configuration that is based on low-fidelity analysis is used as the baseline. Tail lift is included to help tailor the aft portion of the ground signature. A comparison of low- and high-fidelity analysis results demonstrates the necessity of using computational fluid dynamics (CFD) analysis in a low-boom supersonic configuration design process. The fuselage shape is modified iteratively to obtain a configuration with a CFD equivalent-area distribution that matches a predetermined low-boom target distribution. The mixed-fidelity approach can easily refine the low-fidelity low-boom baseline into a low-boom configuration with the use of CFD equivalent-area analysis. The ground signature of the final configuration is calculated by using a state-of-the-art CFD-based boom analysis method that generates accurate midfield pressure distributions for propagation to the ground with ray tracing. The ground signature that is propagated from a midfield pressure distribution has a shaped ramp front, which is similar to the ground signature that is propagated from the CFD equivalent-area distribution. This result confirms the validity of the low-boom supersonic configuration design by matching a low-boom equivalent-area target, which is easier to accomplish than matching a low-boom midfield pressure target.

Li, Wu↗

Anomaly Detection in Test Equipment via Sliding Mode Observers

Nonlinear observers were originally developed based on the ideas of variable structure control, and for the purpose of detecting disturbances in complex systems. In this anomaly detection application, these observers were designed for estimating the distributed state of fluid flow in a pipe described by a class of advection equations. The observer algorithm uses collected data in a piping system to estimate the distributed system state (pressure and velocity along a pipe containing liquid gas propellant flow) using only boundary measurements. These estimates are then used to further estimate and localize possible anomalies such as leaks or foreign objects, and instrumentation metering problems such as incorrect flow meter orifice plate size. The observer algorithm has the following parts: a mathematical model of the fluid flow, observer control algorithm, and an anomaly identification algorithm. The main functional operation of the algorithm is in creating the sliding mode in the observer system implemented as software. Once the sliding mode starts in the system, the equivalent value of the discontinuous function in sliding mode can be obtained by filtering out the high-frequency chattering component. In control theory, "observers" are dynamic algorithms for the online estimation of the current state of a dynamic system by measurements of an output of the system. Classical linear observers can provide optimal estimates of a system state in case of uncertainty modeled by white noise. For nonlinear cases, the theory of nonlinear observers has been developed and its success is mainly due to the sliding mode approach. Using the mathematical theory of variable structure systems with sliding modes, the observer algorithm is designed in such a way that it steers the output of the model to the output of the system obtained via a variety of sensors, in spite of possible mismatches between the assumed model and actual system. The unique properties of sliding mode control allow not only control of the model internal states to the states of the real-life system, but also identification of the disturbance or anomaly that may occur.

Solano, Wanda M.↗

Decentralized Data-Driven Estimation of Generator Rotor Speed and Inertia Constant Based on Adaptive Unscented Kalman Filter

This paper proposes an online decentralized and data-driven method for synchronous generator rotor speed and inertia estimation. First, an analytical relationship between the generator terminal voltage/current phasors and its internal rotor speed is derived using the Thevenin equivalent. The parameters of the latter are obtained via a recursive estimator, which allows us to estimate the generator's internal rotor speed from voltage/current phasors without the need for any generator information. Second, we reformulate the equivalent swing equation into the Kalman filter form, which enables us to use a priori-estimated internal rotor speed along with measured real and reactive power injections to obtain an estimate of the inertia for each online synchronous generator. This is achieved through an adaptive unscented Kalman filter. We demonstrate that the use of terminal bus frequency to approximate the rotor speed for inertia estimation by existing methods leads to large estimation biases/errors. Numerical simulations carried out on the IEEE 39-bus and Texas 2000-bus systems reveal that our proposed method is consistently more accurate than existing methods. Moreover, our proposed method is insensitive to both the type and location of system disturbances, and it is robust to different levels of measurement noise.

bus frequency↗

CFD-DEM Modeling of Autothermal Pyrolysis of Corn Stover with a Coupled Particle- and Reactor-Scale Framework

Autothermal operation of fast pyrolysis is an efficient process-intensification technique wherein exothermic oxidation reactions are used to overcome the heat-transfer bottleneck of conventional pyrolysis. The development of accurate, reliable modeling toolsets is imperative to generating a deeper understanding of biomass autothermal pyrolysis systems to support scale-up and industrial deployment. This modeling effort describes the development of single-particle and reactor models which incorporate detailed reaction schemes and simultaneous exothermic oxidation reactions. The particle-scale model was parameterized for corn stover feedstock with particle morphology, density, ash content, and biopolymer composition, all of which impact the emergent conversion characteristics during pyrolysis. Results were then used to parameterize a reactor-scale autothermal pyrolysis model, which was developed using a coarse-grained computational fluid dynamic-discrete element method. The simulation results compared well with experimental results, with the predicted bio-oil, light gas, and biochar yield within 3.0 wt% of the experimental yields. Further analyses were performed to test the influence of equivalence ratio, biomass injection position, and particle size distribution on autothermal pyrolysis. The analysis of the physio-chemical properties of the fluid and solid phase inside the reactor and at the reactor outlet help reveal important process interactions of autothermal pyrolysis.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Noise dynamics of quantum annealers: estimating the effective noise using idle qubits

Quantum annealing is a type of analog computation that aims to use quantum mechanical fluctuations in search of optimal solutions of QUBO (quadratic unconstrained binary optimization) or, equivalently, Ising problems. Since NP-hard problems can in general be mapped to Ising and QUBO formulations, the quantum annealing paradigm has the potential to help solve various NP-hard problems. Current quantum annealers, such as those manufactured by D-Wave Systems, Inc. have various practical limitations including the size (number of qubits) of the problem that can be solved, the qubit connectivity, and error due to the environment or system calibration, which can reduce the quality of the solutions. Typically, for an arbitrary problem instance, the corresponding QUBO (or Ising) structure will not natively embed onto the available qubit architecture on the quantum chip. Thus, in these cases, a minor embedding of the problem structure onto the device is necessary. However, minor embeddings on these devices do not always make use of the full sparse chip hardware graph, and a large portion of the available qubits stay unused during quantum annealing. In this work, we embed a disjoint random QUBO on the unused parts of the chip alongside the QUBO to be solved, which acts as an indicator of the solution quality of the device over time. Using experiments on three different D-Wave quantum annealers, we demonstrate that (i) long term trends in solution quality exist on the D-Wave device, and (ii) the unused qubits can be used to measure the current level of noise of the quantum system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring Transport Behavior in Hybrid Perovskites Solar Cells via Machine Learning Analysis of Environmental-Dependent Impedance Spectroscopy

Hybrid organic–inorganic perovskites are one of the promising candidates for the next-generation semiconductors due to their superlative optoelectronic properties. However, one of the limiting factors for potential applications is their chemical and structural instability in different environments. Herein, the stability of (FAPbI3)0.85(MAPbBr3)0.15 perovskite solar cell is explored in different atmospheres using impedance spectroscopy. An equivalent circuit model and distribution of relaxation times (DRTs) are used to effectively analyze impedance spectra. DRT is further analyzed via machine learning workflow based on the non-negative matrix factorization of reconstructed relaxation time spectra. This exploration provides the interplay of charge transport dynamics and recombination processes under environment stimuli and illumination. The results reveal that in the dark, oxygen atmosphere induces an increased hole concentration with less ionic character while ionic motion is dominant under ambient air. Under 1 Sun illumination, the environment-dependent impedance responses show a more striking effect compared with dark conditions. In this case, the increased transport resistance observed under oxygen atmosphere in equivalent circuit analysis arises due to interruption of photogenerated hole carriers. The results not only shed light on elucidating transport mechanisms of perovskite solar cells in different environments but also offer an effective interpretation of impedance responses.

36 MATERIALS SCIENCE↗

Application of system identification to analytic rotor modeling from simulated and wind tunnel dynamic test data, part 2

An introduction to aircraft state and parameter identification methods is presented. A simplified form of the maximum likelihood method is selected to extract analytical aeroelastic rotor models from simulated and dynamic wind tunnel test results for accelerated cyclic pitch stirring excitation. The dynamic inflow characteristics for forward flight conditions from the blade flapping responses without direct inflow measurements were examined. The rotor blades are essentially rigid for inplane bending and for torsion within the frequency range of study, but flexible in out-of-plane bending. Reverse flow effects are considered for high rotor advance ratios. Two inflow models are studied; the first is based on an equivalent blade Lock number, the second is based on a time delayed momentum inflow. In addition to the inflow parameters, basic rotor parameters like the blade natural frequency and the actual blade Lock number are identified together with measurement bias values. The effect of the theoretical dynamic inflow on the rotor eigenvalues is evaluated.

Hohenemser, K. H.↗

Interaction of feel system and flight control system dynamics on lateral flying qualities

An experimental investigation of the influence of lateral feel system characteristics on fighter aircraft roll flying qualities was conducted using the variable stability USAF NT-33. Forty-two evaluation flights were flown by three engineering test pilots. The investigation utilized the power approach, visual landing task and up-and-away tasks including formation, gun tracking, and computer-generated compensatory attitude tracking tasks displayed on the Head-Up Display. Experimental variations included the feel system frequency, force-deflection gradient, control system command type (force or position input command), aircraft roll mode time constant, control system prefilter frequency, and control system time delay. The primary data were task performance records and evaluation pilot comments and ratings using the Cooper-Harper scale. The data highlight the unique and powerful effect of the feel system of flying qualities. The data show that the feel system is not 'equivalent' in flying qualities influence to analogous control system elements. A lower limit of allowable feel system frequency appears warranted to ensure good lateral flying qualities. Flying qualities criteria should most properly treat the feel system dynamic influence separately from the control system, since the input and output of this dynamic element is apparent to the pilot and thus, does not produce a 'hidden' effect.

Bailey, R. E.↗

A note on the limitations of lattice least squares

This paper quantifies the known limitation of lattice least squares to ARX models in terms of the dynamic properties of the system being modeled. This allows determination of the applicability of lattice least squares in a given situation. The central result is that an equivalent ARX model exists for an ARMAX system if and only if the ARMAX system has no transmission zeros from the noise port to the output port. The technique used to prove this fact is a construction using the matrix fractional description of the system. The final section presents two computational examples.

Gillis, J. T.↗

Planetary Systems Around Neutron Stars

This project was initiated in 1993, about one year after the announcement of two planets around PSR B1257+12. Its goal was to investigate planetary systems around neutron stars using high precision timing of radio pulsars as a tool. A microsecond precision of the pulse timing analysis, which is equivalent to a millimeter-per-second radial velocity resolution, makes it possible to detect asteroid-mass bodies in orbit around pulsars and to study the dynamics of pulsar planetary systems. The project originally consisted of two longterm efforts: (i) routine observations and timing analysis of the millisecond pulsar PSR B1257+12 which was found to be orbited by at least two earth-mass bodies (Wolszczan and Frail, Nature, 355, 145) and (ii) a sensitive all-sky search for millisecond pulsars to detect further examples of neutron stars with planetary systems. In the third year of the project, it was expanded to include long-term timing observations of slow pulsars in search for planetary systems around these younger neutron stars. The instrumentation used to conduct these investigations included the 305-m Arecibo antenna with the Penn State Pulsar Machine (PSPM-1), the 100-m Effelsberg telescope with the local pulse timing hardware, and the 32-m paraboloid of the Torun Centre for Astronomy in Torun, Poland (TCFA) with the PSPM-2, the second pulsar machine built at Penn State. The PI's collaborators included pulsar groups led by D. Backer (Berkeley), R. Foster (NRL), S. Kulkarni (Caltech), J. Taylor (Princeton) and R. Wielebinski (Bonn). One postdoc (Stuart Anderson), one graduate student (Brian Cadwell) and several undergraduates have been engaged in various aspects of research related to this project.

Wolszczan, Alexander↗

A Mechanistic Study of Nucleate Boiling Heat Transfer Under Microgravity Conditions

Experimental studies of growth and detachment processes of a single bubble and multiple bubbles formed on a heated surface have been conducted in the parabola flights of KC-135 aircraft. Distilled water and PF5060 were used as the test liquids. A micro-fabricated test surface was designed and built. Artificial cavities of diameters 10 microns, 7 microns and 4 microns were made on a thin polished Silicon wafer that was electrically heated by a number of small heating elements on the back side in order to control the surface superheat. Bubble growth period, bubble size and shape from nucleation to departure were measured under subcooled and saturation conditions. Significantly larger bubble departure diameters and bubble growth periods than those at earth normal gravity were observed. Bubble departure diameters as large as 20 mm for water and 6 mm for PF5060 were observed as opposed to about 3 mm for water and less than 1 mm for PF5060 at earth normal gravity respectively. It is found that the bubble departure diameter can be approximately related to the gravity level through the relation D(sub d) proportional 1/g(exp 1/2). For water,the effect of wall superheat and liquid subcooling on bubble departure diameter is found to be small.The growth periods are found to be very sensitive to liquid subcooling at a given wall superheat. However,the preliminary results of single bubble dynamics using PF5060 showed that the departure diameter increases when wall superheat is elevated at the same gravity and subcooling. Growth period of single bubbles in water has been found to vary as t(sub g) proportional g(exp -.93). For water, when the magnitude of horizontal gravitational components was comparable to that of gravity normal to the surface, single bubbles slid along the heater surface and departed with smaller diameter at the same gravity level in the direction normal to the surface. For PF5060, even a very small horizontal gravitational component caused the sliding of bubble along the surface. The numerical simulation has been carried out by solving under the condition of axisymmetry, the mass, momentum, and energy equations for the vapor and the liquid phases. In the model the contribution of micro-layer has been included and instantaneous shape of the evolving vapor-liquid interface is determined from the analysis. Consistent with the experimental results, it is found that effect of reduced gravity is to stretch the growth period and bubble diameter It is found that effect of reduced gravity is to stretch the growth period and bubble diameter at departure. The numerical simulations are in good agreement with the experimental data for both the departure diameters and the growth periods. In the study on dynamics of multiple bubbles, horizontal merger of 2,3 4,and 5 bubbles was observed. It is found that after merger of 2 and 3 bubbles the equivalent diameter of the detached bubble is smaller than that of a single bubble departing at the same gravity level. During and after bubble merger, liquid still fills the space between the vapor stems so as to form mushroom type bubbles. The experimental and numerical studies conducted so far have brought us a step closer to prediction of nucleate boiling heat fluxes under low gravity conditions. Preparations for a space flight are continuing.

Dhir, V. K.↗

The Instrumented Walking and Turning Test to Evaluate Suited Gait Dynamics and Performance in Extravehicular Activity Training Environments

Background and aims: Walking will be required for many exploration tasks on the Moon during the Artemis program. Walking in a straight line on the confined floorspace of a testing area, and repetitive treadmill walking that requires no change in direction may not adequately reflect the balance and coordination required during ambulation. Also, performance of turning maneuvers may be affected differently in different extravehicular activity (EVA) training facilities that simulate partial gravity. For example, the Neutral Buoyancy Lab (NBL) simulates lunar gravity by adding weight to underwater subjects to alter buoyancy and achieve the equivalent ground reaction force of 1/6 of Earth’s gravity (1/6G), whereas the Active Response Gravity Offload System (ARGOS) uses a computer controlled overhead suspension system programmed to continuously offload a percentage of a subject’s weight to simulate 1/6G. The degree to which dynamic movements such as turning are comparable across these EVA training facilities has not yet been evaluated. The instrumented gait test helps NASA scientists and engineers evaluate gait dynamics and performance in suited conditions, and this test demonstrates the unique characteristics and limitations of EVA training facilities. We developed an instrumented walking and turning test using inertial measurement units (IMUs) and conducted the test at NASA’s EVA training facilities. Results were used to compare suited walking and turning characteristics in the ARGOS and the NBL. Methods: Subjects donned the Mark III space suit during offloading with the ARGOS spreader bar gimbal and donned the Z2.5 space suit while underwater in the NBL with weights and floatation added to achieve realistic suit center of gravity. The test team securely attached three Opal (APDM, OR, USA) wireless IMUs on the space suit for each test run: one on the middle of the hard upper torso, and one on the left and on the right ankle bearings. During the NBL tests, the IMUs were encased in a waterproof housing (GoPro) with foam added to create a tighter fit. At both testing facilities, 6.3 m x 1.0 m (LxW) walking lines were marked, and a cone for turning or walking around was located at the end of the walking path with another line on the other side of the cone to indicate the stopping point after walking around the cone. Under simulated 1/6G, subjects began by standing at the marked line with their arms folded across the chest, they then walked at a preferred speed along the straight walking path until they reached the end, turned 180 degrees around the cone, and finally stopped at the marked stopping point. All IMU data recorded during testing were automatically saved to the internal memory. Then, raw IMU signals were processed using custom MATLAB (Mathworks, MA, USA) code to compare gait parameters during both the walking and the turning components of the task. These parameters included time (s), speed (m/s for walking and rad/s for turning), step number (n) and walk:turn time ratio (% time spent straight walking versus turning). Results: Less time, faster gait, fewer steps, and higher walk:turn ratio during both walking and turning components were exhibited during tests performed at the ARGOS versus those performed at the NBL. During the NBL tests, the slower walking speed continued at the same rate throughout a U-shape turn. During the ARGOS tests, the subjects performed shorter and tighter turns at 4 times the speed of the NBL turns because they walked 30% faster and the vertical offloading system gave them more support. Conclusion: Our data show that the differences in walking and turning parameters during the NBL tests may be due to the high viscosity in the water environment where the motion of the lower limbs was slow and did not reach full flexion and extension. These tests improve the current knowledge of testing environments in preparation for EVAs on the lunar surface.

Kyoung Jae Kim↗

Characterization of dynamic stall phenomenon using two-dimensional unsteady Navier-Stokes equations

Among the new significant aspects of the present work are: (1) the treatment of the far-field boundary; (2) the use of C-grid topology, with the branch-cut singularity treated analytically; (3) evaluation of the effect of the envelope of prevailing initial states, and finally; (4) the ability to employ streakline/pathline 'visualization' to probe the unsteady features prevailing in vortex-dominated flows. The far-field boundary is placed at infinity, using appropriate grid stretching. This contributes to the accuracy of the solutions, but raised a number of important issues which needed to be resolved; this includes determining the equivalent time-dependent circulation for the pitching airfoil. A secondary counter-clockwise vortex erupts from within the boundary layer and immediately pinches off the energetic leading-edge shear layer which then, through hydrodynamic instability, rolls up into the dynamic stall vortex. The streakline/pathline visualization serves to provide information for insight into the physics of the unsteady separated flow.

Ghia, K. N.↗

A stochastic-dynamic model for global atmospheric mass field statistics

A model that yields the spatial correlation structure of atmospheric mass field forecast errors was developed. The model is governed by the potential vorticity equation forced by random noise. Expansion in spherical harmonics and correlation function was computed analytically using the expansion coefficients. The finite difference equivalent was solved using a fast Poisson solver and the correlation function was computed using stratified sampling of the individual realization of F(omega) and hence of phi(omega). A higher order equation for gamma was derived and solved directly in finite differences by two successive applications of the fast Poisson solver. The methods were compared for accuracy and efficiency and the third method was chosen as clearly superior. The results agree well with the latitude dependence of observed atmospheric correlation data. The value of the parameter c sub o which gives the best fit to the data is close to the value expected from dynamical considerations.

Ghil, M.↗

Development of a Novel Helmet Support Assembly for NASA Orion Crew Survival Suit

Orion Spacecraft water landing environments have necessitated the design and testing of a novel Helmet Support Assembly (HSA) as part of the NASA Orion Crew Survival Systems (OCSS) launch and re-entry spacesuit. A series of development sled tests using 5 th female, 50th male, and 95th male Hybrid III Anthropometric Test Devices (ATDs) simulated worst case water landings and identified the need for an integrated suit countermeasure to control the helmet during dynamic loading. Initial prototype countermeasures rigidly restrained the helmet keeping the helmet in place as desired, but led to large force and moment transmission to the upper neck. Further development led to a flexible HSA which used high-strength steel wire contoured to the size-matched ATD thorax and attached to both the front and rear of the neck ring. The selected wire diameter and contoured shape provided some ability to flex under the expected worst-case dynamic loads. The existing shoulder harness was used to restrain the helmet motion during “eyeballs out/down” loading while also attempting to prevent head-to-helmet contact in “eyeballs down” spinal-axis loading. A new small prototype helmet was also developed and tested for the 5th and 50th ATDs with both HSA devices. Dynamic impact tests were conducted on the HIA (Horizontal Impulse Accelerator) at Wright Patterson Air Force Base (WPAFB) under simulated off-nominal landing conditions using a representative Orion seat and 5-point harness. For each ATD size, peak Neck Injury Criteria (Nij) were compared to equivalent sled tests with an unsuited ATD configuration. Adding a helmet without attempting to control motion increased peak Nij values for all ATD sizes. The use of the rigid countermeasure showed decreased peak Nij results for the 5th female, but increased peak Nij values for both 50th configurations and the 95th male ATDs. By comparison, the flexible HSA showed reduced peak Nij values for all ATD sizes. Overall, this series of dynamic impact tests identified a risk of increased crew neck loading when under- or over-restraining the helmet during simulated Orion water landings and also demonstrated an effective strategy to mitigate those loads with a compliant HSA design.

Safety↗

A thermodynamic condition for ignition and burn-propagation in cryogenic layer inertially confined fusion implosions

A Lawson-like criterion for ignition (where self-heating dominates over all energy losses) in a dynamic implosion is developed, which accounts for asymmetry and for differences in an implosion x-ray confinement quality. It is shown that the thermodynamic ignition condition is equivalent to yield amplification levels of 16–32. Since negative pdV work of expansion after stagnation increases energy losses above that of x-ray and electron-conduction losses, the Lawson-like ignition criterion is necessary but not sufficient for igniting the hot spot to propagate into the DT fuel with sufficient vigor to generate high gain. Additionally, a higher dimensional generalization of the Lawson-like criterion, which includes the cooling of the implosion upon disassembly, does provide the needed criteria, and it shows that significantly higher temperature is needed for very high levels of yield amplification compared to what traditional ignition metrics imply.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Inlet Area Ratio on Operability of an Axial Air Inlet Rotating Detonation Combustor

This research presents the first results from a new Rotating Detonation Combustor (RDC) design with an axial air inlet and rear-facing angled fuel injection. This design features a variable air inlet that allows for the study of the effect of inlet area ratio on operation and to directly study the coupling between the air plenum and channel of the RDC while minimizing inlet pressure losses. This device maintains the same downstream channel configuration as those previously tested in our laboratory, which includes a 154 mm outer channel diameter and a constant area channel with a width of 7.6 mm. All tests were completed with a downstream hydrogen-air afterburner capable of flashing back in the channel to initiate detonation operation. Initial tests demonstrated successful operation at air mass flow rates from 160 to 770 g/s and equivalence ratios from 0.53 to 1.05 at multiple air inlet to channel area ratios, indicated that sensors installed within the air feed plenum are able to characterize the dynamic characteristics of waves within the detonation channel, and verified that the design goals for this configuration were achieved.

Shepard, Joshua↗

Automated design optimization of supersonic airplane wing structures under dynamic constraints.

The problems of the preliminary and first level detail design of supersonic aircraft wings are stated as mathematical programs and solved using automated optimum design techniques. The problem is approached in two phases: the first is a simplified equivalent plate model in which the envelope, plan form and structural parameters are varied to produce a design, the second is a finite element model with fixed configuration in which the material distribution is varied. Constraints include flutter, aeroelastically computed stresses and deflections, natural frequency and a variety of geometric limitations. The Phase I objective is a combination of weight and drag while Phase II is a weight minimization.

Fox, R. L.↗